안녕하세요? 이번 글은 QGIS SCP 플러그인으로 센티널 2호 위성영상으로부터 분광지수를 계산해 보겠습니다.
첫번째 계산해볼 분광지수는 가장 친숙한 '정규차식생지수(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.560㎛, 8번(NIR, 중심파장 0.842㎛), 11번(SWIR, 중심파장 1.610㎛)입니다.
이번에는 'SCP > Batch(일괄)' 기능을 실행합니다.
'Batch' 창은 아래와 같은데요,
우측 'Functions > !working_dir!'를 더블클릭하고,
작업 디렉터리를 센티널 2호 위성영상 밴드가 위치한 폴더로 지정해 줍니다.
다음으로는 'Functions > band_calc'를 더블 클릭합니다.
NDWI 계산을 위한 해당 구문을 아래와 같이 작성해 봅니다.
위 구문을 복사해서 NDMI를 계산하는 구문도 아래와 같이 작성해 봅니다.
이제 'RUN'을 실행하면, Band calc을 두 번 실행한 것처럼 NDWI와 NDMI가 일괄 생성됩니다. 본인이 선호하는 방법을 쓰시면 되겠습니다.