QGIS 3.4 SCP 플러그인 - (3) 센티널 2호 위성영상 분광지수 계산

댓글수0 다음블로그 이동

REMOTE SENSING

QGIS 3.4 SCP 플러그인 - (3) 센티널 2호 위성영상 분광지수 계산

유병혁
댓글수0

안녕하세요? 이번 글은 QGIS 3.4 SCP 플러그인으로 센티널 2호 위성영상으로부터 분광지수를 계산해 보겠습니다.


분광지수 관련해서는 GIS 및 지구관측을 위한 온라인 학습플랫폼 'GEO University'의 글이 참고가 되실 것 같습니다.

Spectral Indices with multispectral satellite data | https://www.geo.university/pages/spectral-indices-with-multispectral-satellite-data


첫번째 계산해볼 분광지수는 가장 친숙한 '정규차식생지수(NDVI: Normalized Difference Vegetation Index)'입니다.

NDVI는 센티널 2호 영상의 4번 밴드(Red, 중심파장 0.665㎛)와 8번 밴드(NIR, 중심파장 0.842㎛)로 계산됩니다.


일단 QGIS 3.4에서 센티널 2호 위성영상의 4번과 8번 밴드 레이어를 추가합니다.


'SCP > Band calc'를 실행합니다.


'Band calc' 창의 Band list에서 각각의 레이어(Band name)가 변수(Variable)로 정의된 것을 확인하실 수 있습니다.


'expression'에서 NDVI 계산식을 아래와 같이 적용하고 'RUN' 버튼을 클릭해 봅니다.


결과는 다음과 같습니다. 북한산국립공원 영역이 잘 드러나죠?!


이번에는 SCP 플러그인의 'Batch()' 기능을 통해 한번에 여러 분광지수를 계산해 보겠습니다. 정규차수분지수

(NDWI: Normalized Difference Water Index)와 정규차식생수분지수(NDMI: Normalized Difference Moisture Index)

를 일괄 처리해 봅니다.


분광지수 계산에 필요한 밴드는 3번(Green, 중심파장 0.5608번(NIR, 중심파장 0.842㎛), 11번(SWIR, 중심파장 1.610㎛)입니다.


이번에는 'SCP > Batch(일괄)' 기능을 실행합니다.


'Batch' 창은 아래와 같은데요,


우측 'Functions > !working_dir!'를 더블클릭하고,


작업 디렉터리를 센티널 2호 위성영상 밴드가 위치한 폴더로 지정해 줍니다.


다음으로는 'Functions > band_calc'를 더블 클릭합니다.


NDWI 계산을 위한 해당 구문을 아래와 같이 작성해 봅니다.


위 구문을 복사해서 NDMI를 계산하는 구문도 아래와 같이 작성해 봅니다.


이제 'RUN'을 실행하면, Band calc을 두 번 실행한 것처럼 NDWI와 NDMI가 일괄 생성됩니다. 본인이 선호하는 방법을 쓰시면 되겠습니다.

맨위로

http://blog.daum.net/geoscience/1446

신고하기