QGIS에서 SCP 플러그인을 이용하여 위성영상을 다운로드 받는 방법을 학습해보도록 하겠습니다.
SCP는 원격탐사를 위해 특화된 QGIS 플러그인으로, Semi-Automatic Classification Plugin의 약어입니다.
사이트: https://fromgistors.blogspot.com/p/semi-automatic-classification-plugin.html | 개발자: 루카 콘지도(Luca Congedo)
먼저, SCP를 설치해보도록 하겠습니다.
상단 메뉴에서 '플러그인 > 플러그인 관리 및 설치'를 클릭하고,
Semi-Automatic Classification Plugin을 검색하고 설치합니다.
같은 방식으로 구글지도, 빙맵, 오픈스트리트맵 등 레이어를 이용할 수 있는 'OpenLayers Plugin'도 설치하도록 하겠습니다.
저는 이전글에 이어 스리랑카 콜롬보지역의 위성영상을 SCP를 이용하여 다운로드 받아보도록 하겠습니다.
먼저, 아래와 같이 래스터 레이어 추가 버튼을 클릭하고,
앞서 랜드 뷰어(Land Viewer)에서 다운로드 받은 2016년 1월 27일 취득된 랜드셋 8호 위성영상 중,
(http://blog.daum.net/geoscience/992) 해상도가 가장 높은 15m급 8번 전정색 밴드를 열어보겠습니다.
아래와 같이 콜롬보 지역을 포함한 위성영상이 추가되었습니다.
전체 영상에서 콜롬보 지역이 어디쯤인지 확인하기 위해 OpenLayers를 이용해 보겠습니다.
상단 메뉴에서 '웹 > OpenLayers plugin > OpenStreepMap > OpenStreetMap'을 클릭해 보면,
아래와 같이 해당 지역의 오픈스트리트맵 레이어가 추가됩니다.
참고로, 콜롬보 지역은 스리랑카의 서부지역에 위치하고 있습니다.
콜롬보 지역을 확대했는데요, 이제 이 영역에 대한 위성영상을 조회하고 다운로드 받아보도록 하겠습니다.
상단 메뉴에서 'SCP > Download images > Landsat download'를 클릭합니다.
아래와 같이 SCP는 미국 NASA와 USGS가 협력 개발한 Landsat 위성영상 뿐만 아니라, 유럽 ESA의 Sentinel-2 위성영상,
NASA와 일본 경제산업성(METI, Ministry of Economy, Trade asnd Industry)이 협력한 ASTER 위성영상 다운로드를 제공합니다.
실행 창은 아래와 같이 Semi-Automatic Classfication Plugin 창이 실행되고, Landsat download 탭이 선택됩니다.
아래와 같이 Sentinel-2 download와 ASTER download 탭도 비슷한 기능 구성을 가지고 있습니다.
자, 그럼 Landsat download 탭을 자세히 살펴보도록 하겠습니다.
먼저, 해당 위성 아카이브 공식 사이트의 계정(https://ers.cr.usgs.gov)이 필요합니다.
Search area는 QGIS 화면상에서 상단 좌측(UL: Upper Left), 하단 우측(LR: Lower Right) 좌표를 설정하는 기능입니다.
화면 우측 '+' 버튼을 클릭한 후, UL은 마우스 좌클릭, LR은 마우스 우클릭을 통해 사각영역을 아래와 같이 지정할 수 있습니다.
아래와 같이 해당 영역의 좌표가 자동으로 입력되었는데요, Search area와 중복되는 위성영상들을 조회하게 됩니다.
Search는 검색할 위성영상의 범위를 지정하는 기능입니다.
Satellites를 보시면 Landsat 위성 시리즈를 선택하실 수 있는데요, 저는 L4-5 TM을 선택하도록 하겠습니다.
Date(취득날짜)는 2006년 1월 1일부터 12월 31일까지, Max cloud cover(최대 토지피복률)는 10%로 설정한 후 'Find'를 클릭해보겠습니다.
아래와 같이 검색결과가 조회되었는데요, CloudCover 순으로 조회한 후, 2006년 2월 7일 위성영상을 다운로드 받아보도록 하겠습니다.
화면 우측에 3개의 버튼이 있는데요, 위로부터 각각 미리보기(display preview), 선택된 행 삭제(delete row), 재설정(reset) 기능들입니다.
'미리보기'를 클릭하면 QGIS 상에 저해상도 레이어가 아래와 같이 추가됩니다. 위성영상 영역을 파악하는데 참고하시면 되겠습니다.
Download options에서는 밴드별 다운로드 선택을 하실 수 있습니다.
현재는 아마존 웹 서비스(AWS)를 통해 Landsat 8 bands들에 한해 적용되는 기능입니다.
이제 Image list에서 원하는 영상만 정리하고 다운로드를 진행해 보도록 하겠습니다.
아래 'Preprocess images'를 체크하시면 영상 DN값이 반사율(reflectance)로 자동 변환됩니다.
실행(Run) 버튼을 클릭하고,
위성영상을 다운로드 받을 폴더를 지정해 줍니다.
QGIS 화면 상에서 아래와 같이 Executing 바가 추가되어 진행상태(%)를 알려줍니다.
앞서 'Load bands in QGIS'를 체크해뒀기 때문에 아래와 같이 자동 계산된 반사율 영상이 추가됩니다.
다운로드 위치에 보시면 폴더가 2개 존재하는데요, 하나는 기본 DN값, _converted 폴더는 반사율 값이 저장되어 있습니다.
좌측은 DN값 폴더, 우측은 반사율값 폴더입니다.
