이번 글은 QGIS에서 경사방향-경사 지도(Slope-Aspect Map) 제작과정을 정리해 보겠습니다.
Aspect-Slope Maps in QGIS | https://kingsgeocomputation.org/2016/03/16/aspect-slope-maps-in-qgis/
경사방향-경사 지도는 지형(또는 연속적인 면)의 경사방향(direction)과 경사(steepness)를 동시에 보여줍니다.
경사방향 범주는 색상(hue)으로, 경사도 등급은 채도(saturation)로 표현됩니다.
앞서 Esri 블로그에서는 ArcGIS를 이용하여 경사방향-경사 지도 제작과정을 아래와 같이 정리한 바가 있는데요,
제가 참조한 글은 이것을 QGIS에서 구현한 것입니다. 저는 여기에 QGIS 스타일 기능의 장점을 보완하였습니다.
Esri 블로그: Aspect-slope map (2008)
https://blogs.esri.com/esri/arcgis/2008/05/23/aspect-slope-map/
Esri 블로그: New Aspect-Slope Raster Function Now Available (2017)
https://blogs.esri.com/esri/arcgis/2017/03/28/new-aspect-slope-raster-function-now-available/
경사방향-경사 지도 제작을 위한 실험 데이터는 서울시 AW3D30 DSM 데이터셋입니다.
[1] 전세계 30m급 수치표면모델(DSM) 'AW3D30' 다운로드 소개 | http://blog.daum.net/geoscience/1232
[2] QGIS 2.18: AW3D30 데이터셋을 이용한 구역 통계(zonal statistics) | http://blog.daum.net/geoscience/1233
현재 위 데이터셋의 CRS는 EPSG:4326(지리 좌표계)으로 정의되어 있는데요, 이것을 EPSG:5186(투영 좌표계)로 변경하겠습니다.
일단 DSM 데이터를 경사방향(Aspect)과 경사(Slope) 데이터로 각각 변환해 보겠습니다.
'래스터 > 분석 > DEM (지형 모델)'을 클릭합니다.
아래와 같이 경사 래스터를 생성합니다. 이 때 경사는 도단위 대신 백분율로 표시합니다.
이번에는 경사방향 래스터를 생성했습니다.
자, 이번에는 앞서 생성한 경사, 경사방향 래스터를 규칙에 따라 재분류하도록 하겠습니다.
이를 위해, QGIS 공간처리 툴박스에서 GRASS GIS 7 commands > Raster (r.*) > r.reclass를 이용합니다.
r.reclass | https://grass.osgeo.org/grass70/manuals/r.reclass.html
먼저, 경사 래스터는 아래와 같은 재분류 규칙(reclass rules)을 갖습니다.
재분류 규칙을 새 txt 파일로 생성한 후, r.reclass에서 사용하는 문법으로 작성합니다.
r.reclass실행 시, 아래 재분류 규칙 txt 파일을 반영하여 경사 래스터를 재분류합니다.
slope_recl.txt0-5 : 10
5-20 : 20
20-40 : 30
40-highest value : 40
이번에는 경사방향 래스터를 재분류해 볼까요?! 재분류 규칙은 아래와 같습니다.
0-22.5 : 1
22.5-67.5 : 2
67.5-112.5 : 3
112.5-157.5 : 4
157.5-202.5 : 5
202.5-247.5 : 6
247.5-292.5 : 7
292.5-337.5 : 8
337.5-359.5 : 1
마찬가지로 r.reclass에 재분류 규칙 txt 파일을 반영하여 경사방향 래스터를 재분류합니다.
자, 이제 '래스터 > 래스터 계산기'에서 경사 재분류와 경사방향 재분류를 더해주면 '경사방향-경사 지도'가 제작됩니다.
아직까지 '그레이 레벨'이라 크게 와닿지 않는데요, 이 지도에 경사방향-경사 래스터의 색상표를 적용해 보겠습니다.
속성의 스타일 탭에서,
하단 '스타일'을 클릭하고 '스타일 불러오기'를 선택합니다.
아래 'QGIS 레이어 스타일 파일(*.qml)'를 지정합니다.
아래와 같이 경사방향-경사 래스터에 스타일까지 적용되었습니다.
ArcGIS의 *.lyr 파일처럼 QGIS도 '레이어 정의 파일로 저장'할 수 있는데요, 이 때 확장자는 'QGIS 레이어 정의 파일(*.qlr)'입니다.
