안녕하세요? 이번 글은 행정구역별로 토지피복 분류별 면적을 산정하는 방법을 정리해 보겠습니다.
대분류 토지피복지도를 사용할 건데요, 지도 다운로드, 편집, 스타일링 방법은 아래 글을 보시면 됩니다.
환경공간정보서비스: 토지피복지도 다운로드 / 병합 / 스타일링
안녕하세요? 이번 글은 환경공간정보서비스에서 우리나라 대분류 토지피복지도를 내려받은 후 래스터 병합, 스타일 적용 과정을 정리해 보겠습니다. 일단, 환경공간정보서비스에 접속해 보겠
blog.daum.net
아래와 같이 대분류 토지피복지도가 있습니다.
우리는 여기에 행정구역 경계를 얹은 후, 각 행정구역 단위로 토지피복 분류별 면적을 계산하고자 합니다.
시군구 경계는 국가공간정보포털 오픈마켓(data.nsdi.go.kr/dataset/15144)에서 강원 자료를 확보하였습니다.
내려받은 데이터를 올려보니 아래와 같이 좌표계 편차가 존재합니다.
시군구 경계의 추가 정보를 통해 해당 파일이 Bessel 중부원점(TM), EPSG:5174로 정의됨을 확인합니다.
레이어 좌표계 정의를 변경하니 좌표계 편차가 해소되었습니다.
이제 토지피복 분류별 면적을 계산해 보겠습니다.
해당 기능은 '공간 처리 툴박스 > 구획 히스토그램'을 이용하시면 됩니다.
실행 창에서 래스터 레이어는 토지피복지도를, 밴드 번호는 밴드 1을, 구역을 담고 있는 벡터 레이어는 시군구 경계를, 그리고 출력 열 접두어는 'L1_'를 입력합니다. 출력 구역의 경로, 이름을 설정하고 '실행' 버튼을 클릭합니다(면적 값만을 얻고자 한다면 확장자는 *.xlsx로 설정하시면 됩니다).
결과는 다음과 같습니다. L1_0, ..., L1_7은 대분류 토지피복지도 픽셀 값 0부터 7까지 각각의 개수를 나타냅니다. 대분류 토지피복지도는 30x30m 해상도로 정의되었기 때문에 픽셀 하나 당 면적은 900m로 산정하시면 됩니다. 참고로 L1_NODATA는 픽셀 값이 존재하지 않는 구역으로 L1_0과 대비되는 값입니다. 이상으로 행정구역으로 토지피복 분류별 면적을 계산하는 과정을 정리해 봤습니다.
