GGRS: Geoscience, GIS, & Remote Sensing

QGIS 2.18: AW3D30 데이터셋을 이용한 경사방향-경사 지도(Aspect-Slope Map) 제작

작성일 작성자 유병혁

이번 글은 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.txt

  • 0-5 : 10

  • 5-20 : 20

  • 20-40 : 30

  • 40-highest value : 40



이번에는 경사방향 래스터를 재분류해 볼까요?! 재분류 규칙은 아래와 같습니다.

aspect_recl.txt

  • 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)'를 지정합니다.

slope-aspect_color.qml


아래와 같이 경사방향-경사 래스터에 스타일까지 적용되었습니다.

ArcGIS의 *.lyr 파일처럼 QGIS도 '레이어 정의 파일로 저장'할 수 있는데요, 이 때 확장자는 'QGIS 레이어 정의 파일(*.qlr)'입니다.




맨위로
통합 검색어 입력폼