REMOTE SENSING

SNAP에서 센티널 3호 위성영상 재투영하고 GeoTIFF로 내보내기

유병혁 2019. 8. 15. 13:11

안녕하세요? 이번 글은 SNAP에서 센티널 3호 위성영상을 재투영하고 GeoTIFF로 내보내는 과정을 정리해 보겠습니다.


일단 SNAP을 다운로드 받고 설치합니다.

SNAP Download | http://step.esa.int/main/download/snap-download/


설치화면은 아래와 같은데요, 현재 버전은 2019년 7월 22일 출시된 7.0입니다.


설치 구성요소는 SNAP과 Sentinel-1, Sentinel-2, Sentinel-3, Radarsat-2 도구상자(Toolbox)가 있습니다.


참고로 SNAP은 파이썬 개발자를 위한 SNAP API도 제공합니다: http://step.esa.int/main/community/developers/


SNAP 설치를 완료했습니다. 이제 SNAP에서 센티널 3호 위성영상을 열어볼까요?! 다운로드 방법은 아래 글을 참조하시기 바랍니다.

코페르니쿠스 오픈 엑세스 허브(Copernicus Open Access Hub)의 위성영상 다운로드 방법 소개 | http://blog.daum.net/geoscience/1399


이 글에서 사용하는 데이터는 센티널 3호 위성영상 중 SLSTR 센서 데이터입니다.

SLSTR 데이터는 총 11개 밴드로 구성되어 있으며, 이 중 가시 단파 적외선 채널 S1 - S6는

500m 공간 해상도로, 열 적외선 채널 S7 - S9과 산불 채널 F1, F2는 1km 공간 해상도로 저장됩니다.


SNAP 실행화면은 다음과 같습니다.


센티널 3호 SLSTR 데이터를 열어볼까요?! 아래 'Open Product' 버튼을 클릭하고,


다운로드 받고 압축 해제한 데이터 폴더 내에 위치한 xfdumanifest.xml 파일을 열어줍니다.


아래와 같이 Product Explorer 도구 창에 해당 데이터셋이 표시됩니다.


Bands 폴더 목록 중 LST(Land Surface Temperature) 데이터를 열어보겠습니다.


LST 화면은 아래와 같습니다. 재투영을 적용해 보겠습니다.


SNAP에서 재투영 기능은 상단 메뉴 'Raster > Geometric Operations > Reprojection'을 실행하시면 됩니다.


아래와 같이 Reprojection 창이 실행되었는데요,


'Reprojection Parameters' 탭을 보시면 기본 재투영이 'Geographic Lat/Lon (WGS 84)'임을 확인하실 수 있습니다.


아래와 같이 재투영이 적용되었습니다.


Pixel Info 도구 창을 통해 해당 밴드에서 커서가 위치한 곳의 화소 정보를 확인하실 수 있습니다.

아래와 같이 이 밴드는 2019년 8월 14일 오전 1시 39분(UTC 시간)에 취득되었으며 LST 밴드의

단위는 K(Kelvin)인 것을 확인하였습니다. 


이번에는 재투영된 데이터셋을 선택한 상태에서 GeoTIFF로 내보내는 방법을 학습해 보겠습니다.


상단 메뉴에서 'File > Export > GeoTIFF'를 실행하면,


아래와 같이 'Export Product' 창이 실행됩니다.


여기서 'Subset' 버튼을 클릭하고,


'Specify Product Subset' 창을 통해 내가 원하는 만큼의 공간, 밴드, 메타데이터를 서브셋(Subset)으로 정의합니다.


여기서는 'Band Subset' 탭을 선택하고, 


전체 밴드들 중 아래 3개 밴드만 선택해주도록 하겠습니다.


이제 서브셋 설정까지 마쳤으면 'Export Product' 버튼을 클릭합니다.


내보내기한 GeoTIFF를 QGIS 3.4에서 열어본 화면입니다. 참고로 배경지도는 Vworld Satellite 레이어를 사용했습니다.