안녕하세요? 이번 글은 국토정보플랫폼에서 정사영상을 다운로드 받아 QGIS에서 사용하는 과정을 정리해 보겠습니다. 이전 글(https://blog.daum.net/geoscience/1691)에서는 수치지도를 받았었는데요, 이번에는 같은 지역의 정사영상을 다운로드 받겠습니다.
'정사영상'을 클릭하고,
목록 중 '제주 권역(2019, 영암/34603099)'을 선택해 보겠습니다. '미리보기'를 클릭하면 해당 정사영상을 확인하실 수 있습니다.
'영상비교'를 클릭하면,
아래와 같이 동일 지역을 촬영한 다중 시기 정사영상(여기서는 2019년, 2017년, 2015년, 2010년)을 비교할 수 있습니다.
'메타정보'를 클릭하시면, 해당 정사영상의 메타데이터, 연혁, 식별, 좌표계, 배포, 데이터품질을 확인할 수 있습니다.
해당 정사영상을 다운로드한 후 레이어 추가해 보겠습니다.
현재 좌표 및 좌표계가 정의되어 있지 않은데요, 실습을 통해 지리참조된 정사영상으로 변환해 보겠습니다. 그래야 수치지도와 정사영상을 중첩해볼 수 있겠죠?!
먼저 수치지도 레이어의 공간범위를 확인해 보겠습니다. 앞서 작업한 '34603099_LN' 레이어 속성 창을 보시면, 공간범위를 쉽게 확인하실 수 있습니다.
공간범위 값은 다음과 같습니다:
<좌측 x>,<하단 y> : <우측 x>,<상단 y>
172530.4700000000011642,239394.7500000000000000 : 174827.1900000000023283,242174.7200000000011642
다음으로 '공간 처리 툴박스 > GDAL > 래스터 변환 > 변환 (포맷 변경)'을 클릭합니다.
아래와 같이 '변환 (포맷 변경)' 창이 실행되었는데요,
'입력 레이어'는 정사영상을, '산출 파일 용 투영체 무시'는 수치지도와 동일하게 EPSG:5186을 지정해 줍니다.
앞서 확인한 공간범위는 '추가 명령줄 파라미터'에서 적용하시면 됩니다.
아래 추가 명령줄을 통해 공간범위를 설정할 수 있습니다. 명령줄은 '-a_ullr <좌측 x> <상단 y> <우측 x> <하단 y>' 순으로 앞서 확인한 공간범위와는 순서 차이(<좌측 x>,<하단 y> : <우측 x>,<상단 y>)가 있음을 유의해야 합니다.
또한 정사영상은 수치지도 도엽단위로 분할하여 저장하며, 도곽보다 50m 크게 제작합니다. 따라서 '-a_ullr <좌측 x-50> <상단 y+50> <우측 x+50> <하단 y-50>'으로 값을 조정해야 합니다. 실습에서 조정 값은 아래와 같습니다:
-a_ullr 172480.4700000000011642 242224.7200000000011642 174877.1900000000023283 239344.7500000000000000
자, 이제 '추가 명령줄 파라미터'를 입력하고,
'변환한 산출물'을 지정한 후 결과를 확인해 봅니다.
결과는 다음과 같습니다. 한번 확대해 볼까요?!
자, 이제 수치지도와 지리참조된 정사영상을 중첩해서 사용할 수 있습니다!
