7.3. Lesson: 지형 분석
래스터가 표현하고 있는 지형에 대해 더 깊이 이해할 수 있게 해주는 특정한 래스터 형식이 있습니다. 수치 표고 모델(DEM)이 이런 점에서 특히 유용합니다. 이 강의에서는 지형 분석을 이용해서, 이전 모듈에서 다뤘던 주거 구역 개발에 대해 더 자세히 알아보도록 하겠습니다.
이 강의의 목표: 지형 분석을 이용해서 지형에 대해 더 많은 정보를 끌어내기.
7.3.1. Follow Along: 음영기복도 계산
We are going to use the same DEM layer as in the previous lesson.
If you are starting this chapter from scratch, use the
Browser panel and load the
The DEM layer shows you the elevation of the terrain, but it can sometimes seem a little abstract. It contains all the 3D information about the terrain that you need, but it doesn’t look like a 3D object. To get a better impression of the terrain, it is possible to calculate a hillshade, which is a raster that maps the terrain using light and shadow to create a 3D-looking image.
We are going to use algorithms in themenu.
Click on themenu
The algorithm allows you to specify the position of the light source: Azimuth has values from 0 (North) through 90 (East), 180 (South) and 270 (West), while the Vertical angle sets how high the light source is (0 to 90 degrees).
We will use the following values:
Z factor at
Azimuth (horizontal angle) at
Vertical angle at
Save the file in a new folder
exercise_data/raster_analysis/with the name
Finally click on Run
이제 다음과 같이 hillshade 라는 새 레이어가 생성되었습니다.
멋진 3D로 보입니다만, 좀 더 향상시킬 수 있을까요? 음영기복도 그 자체로는 마치 석고 모형처럼 보입니다. 다른, 좀 더 다채로운 래스터들과 어떻게 함께 사용할 수는 없을까요? 물론 할 수 있습니다. 음영기복도를 오버레이로서 사용한다면 말이죠.
7.3.2. Follow Along: 음영기복도를 오버레이로 사용
음영기복도는 낮 동안의 특정 시간의 햇빛에 대한 유용한 정보를 제공할 수 있습니다. 하지만 맵을 더 보기 좋게 하기 위한 심미적 목적을 위해서도 사용할 수 있습니다. 이렇게 하려면 음영기복도를 거의 투명하게 설정하면 됩니다.
Change the symbology of the original srtm_41_19 layer to use the Pseudocolor scheme as in the previous exercise
Hide all the layers except the srtm_41_19 and hillshade layers
Click and drag the srtm_41_19 to be beneath the hillshade layer in the Layers panel
Set the hillshade layer to be transparent by clicking on the Transparency tab in the layer properties
Set the Global opacity to
You’ll get a result like this:
Switch the hillshade layer off and back on in the Layers panel to see the difference it makes.
음영기복도를 이런 식으로 사용하면 경관의 지형을 향상시킬 수 있습니다. 만약 그 효과가 미미하다고 느껴진다면, hillshade 레이어의 투명도를 조절할 수 있습니다. 하지만 물론, 음영기복도가 밝아질수록 그 아래의 색상은 침침해질 것입니다. 여러분이 만족할 수 있는 균형을 찾아야 합니다.
Remember to save the project when you are done.
7.3.3. Follow Along: Finding the best areas
Think back to the estate agent problem, which we last addressed in the Vector Analysis lesson. Let us imagine that the buyers now wish to purchase a building and build a smaller cottage on the property. In the Southern Hemisphere, we know that an ideal plot for development needs to have areas on it that:
with a slope of less than 5 degrees
But if the slope is less than 2 degrees, then the aspect doesn’t matter.
Let’s find the best areas for them.
7.3.4. Follow Along: 경사도 계산
Slope informs about how steep the terrain is. If, for example, you want to build houses on the land there, then you need land that is relatively flat.
To calculate the slope, you need to use thealgorithm of the .
Open the algorithm
Choose srtm_41_19 as the Elevation layer
Keep the Z factor at
Save the output as a file with the name
slopein the same folder as the
Click on Run
Now you’ll see the slope of the terrain, each pixel holding the corresponding slope value. Black pixels show flat terrain and white pixels, steep terrain:
7.3.5. Try Yourself Calculating the aspect
Aspect is the compass direction that the slope of the terrain faces. An aspect of 0 means that the slope is North-facing, 90 East-facing, 180 South-facing, and 270 West-facing.
Since this study is taking place in the Southern Hemisphere, properties should ideally be built on a north-facing slope so that they can remain in the sunlight.
Use the Aspect algorithm of the
to get the
layer saved along with the
7.3.6. Follow Along: Finding the north-facing aspect
Now, you have rasters showing you the slope as well as the aspect, but you have no way of knowing where ideal conditions are satisfied at once. How could this analysis be done?
그 답은 Raster calculator 가 가지고 있습니다.
QGIS has different raster calculators available:
Each tool is leading to the same results, but the syntax may be slightly different and the availability of operators may vary.
We will usein the Processing Toolbox
Open the tool by double clicking on it.
The upper left part of the dialog lists all the loaded raster layers as
nameis the name of the layer and
Nis the band.
In the upper right part you will see a lot of different operators. Stop for a moment to think that a raster is an image. You should see it as a 2D matrix filled with numbers.
North is at 0 (zero) degrees, so for the terrain to face north, its aspect needs to be greater than 270 degrees or less than 90 degrees. Therefore the formula is:
aspect@1 <= 90 OR aspect@1 >= 270
Now you have to set up the raster details, like the cell size, extent and CRS. This can be done manually or it can be automatically set by choosing a
Reference layer. Choose this last option by clicking on the … button next to the Reference layer(s) parameter.
In the dialog, choose the aspect layer, because we want to obtain a layer with the same resolution.
Save the layer as
The dialog should look like:
Finally click on Run.
결과는 다음과 같을 것입니다.
The output values are
What does it mean?
For each pixel in the raster, the formula we wrote returns whether it matches
the conditions or not.
Therefore the final result will be False (0) and True (1).
7.3.7. Try Yourself More criteria
Now that you have done the aspect, create two new layers from the DEM.
The first shall identify areas where the slope is less than or equal to
The second is similar, but the slope should be less than or equal to
Save them under
7.3.8. Follow Along: 래스터 분석 결과 결합
Now you have generated three raster layers from the DEM:
aspect_north: terrain facing north
slope_lte2: slope equal to or below 2 degrees
slope_lte5: slope equal to or below 5 degrees
Where the condition is met, the pixel value is
Elsewhere, it is
Therefore, if you multiply these rasters, the pixels that have a value
1 for all of them will get a value of
1 (the rest will get
The conditions to be met are:
at or below 5 degrees of slope, the terrain must face north
at or below 2 degrees of slope, the direction that the terrain faces does not matter.
Therefore, you need to find areas where the slope is at or below five
AND the terrain is facing north,
OR the slope is at or
below 2 degrees. Such terrain would be suitable for development.
이 기준들을 만족하는 지역을 계산하려면,
Open the Raster calculator again
Use this expression in Expression:
( aspect_north@1 = 1 AND slope_lte5@1 = 1 ) OR slope_lte2@1 = 1
Set the Reference layer(s) parameter to
aspect_north(it does not matter if you choose another - they have all been calculated from
Save the output under
결과물은 다음과 같습니다.
The previous steps could have been simplified using the following command:
((aspect@1 <= 90 OR aspect@1 >= 270) AND slope@1 <= 5) OR slope@1 <= 2
7.3.9. Follow Along: 래스터 단순화
As you can see from the image above, the combined analysis has left us with many, very small areas where the conditions are met (in white). But these aren’t really useful for our analysis, since they are too small to build anything on. Let us get rid of all these tiny unusable areas.
Open the Sieve tool ( in the Processing Toolbox)
Set the Input file to
all_conditions, and the Sieved to
Set the Threshold to 8 (minimum eight contiguous pixels), and check Use 8-connectedness.
Once processing is done, the new layer will be loaded.
What is going on? The answer lies in the new raster file’s metadata.
View the metadata under the Information tab of the Layer Properties dialog. Look the
This raster, like the one it is derived from, should only feature the values
0, but it has also a very large negative number. Investigation of the data shows that this number acts as a null value. Since we are only after areas that weren’t filtered out, let us set these null values to zero.
Open the Raster Calculator, and build this expression:
(all_conditions_sieve@1 <= 0) = 0
This will maintain all non-negative values, and set the negative numbers to zero, leaving all the areas with value
Save the output under
산출물은 다음과 같습니다.
우리가 기대한 산출물입니다. 이전 결과물을 단순화시킨 버전이죠. 어떤 도구의 결과물이 기대했던 것과 다를 경우, 문제를 해결 하는 데 메타데이터(및, 적용 가능할 경우 벡터 속성)을 살펴보는 작업이 필수적일 수도 있습니다.
7.3.10. Follow Along: Reclassifying the Raster
We have used the Raster calculator to do calculations on raster layers. There is another powerful tool that we can use to extract information from existing layers.
Back to the
We know now that it has numerical values within a range from 0 through
What we want to do is to reclassify this layer to other discrete
values (from 1 to 4), depending on the aspect:
1 = North (from 0 to 45 and from 315 to 360);
2 = East (from 45 to 135)
3 = South (from 135 to 225)
4 = West (from 225 to 315)
This operation can be achieved with the raster calculator, but the formula would become very very large.
The alternative tool is the Reclassify by table tool in in the Processing Toolbox.
Open the tool
Choose aspect as the
Input raster layer
Click on the … of Reclassification table. A table-like dialog will pop up, where you can choose the minimum, maximum and new values for each class.
Click on the Add row button and add 5 rows. Fill in each row as the following picture and click OK:
The method used by the algorithm to treat the threshold values of each class is defined by the Range boundaries.
Save the layer as
Click on Run
If you compare the native aspect layer with the
reclassified one, there are not big differences.
But by looking at the legend, you can see that the values go from
Let us give this layer a better style.
Open the Layer Styling panel
Choose Paletted/Unique values, instead of Singleband gray
Click on the Classify button to automatically fetch the values and assign them random colors:
The output should look like this (you can have different colors given that they have been randomly generated):
With this reclassification and the paletted style applied to the layer, you can immediately differentiate the aspect areas.
7.3.11. Follow Along: Querying the raster
Unlike vector layers, raster layers don’t have an attribute table. Each pixel contains one or more numerical values (singleband or multiband rasters).
All the raster layers we used in this exercise consist of just one band. Depending on the layer, pixel values may represent elevation, aspect or slope values.
Select the tool from the Attributes toolbar.
Click on a random location of the srtm_41_19 layer. Identify Results will appear with the value of the band at the clicked location:
You can change the output of the Identify Results panel from the current
treemode to a
tableone by selecting Table in the View menu at the bottom of the panel:
Clicking each pixel to get the value of the raster could become annoying after a while. We can use the Value Tool plugin to solve this problem.
In the All tab, type
value tin the search box
Select the Value Tool plugin, press Install Plugin and then Close the dialog.
The new Value Tool panel will appear.
If you close the panel you can reopen it by enabling it in theor by clicking on the icon in the toolbar.
To use the plugin just check the Enable checkbox and be sure that the
srtm_41_19layer is active (checked) in the Layers panel.
Move the cursor over the map to see the value of the pixels.
But there is more. The Value Tool plugin allows you to query all the active raster layers in the Layers panel. Set the aspect and slope layers active again and hover the mouse on the map:
7.3.12. In Conclusion
You’ve seen how to derive all kinds of analysis products from a DEM. These include hillshade, slope and aspect calculations. You’ve also seen how to use the raster calculator to further analyze and combine these results. Finally you learned how to reclassify a layer and how to query the results.
7.3.13. What’s Next?
이제 여러분은 두 가지 분석을 할 수 있습니다. 벡터 분석은 잠재적으로 적합한 계획을 보여주며, 래스터 분석은 잠재적으로 적합한 지형을 보여줍니다. 이들을 어떻게 결합해야 문제에 대한 최종 결론을 얻을 수 있을까요? 이것이 다음 모듈에서 시작하는 강의의 주제입니다.