QGIS GRASS commands를 이용하여 SPOT-VGT DN 값을 월 최대 식생지수(MNDVI)로 합성하는 방법을 학습해 보도록 하겠습니다.
이번 학습은 아래 FOSS4G KOREA 2014 구두발표에서 소개한 내용을 글로 정리한 것입니다.
QGIS를 이용한 보호지역의 원격탐사 활용사례 소개(FOSS4G KOREA 2014)
http://blog.daum.net/geoscience/679
SPOT-VGT DN값을 Real NDVI로 변환하는 방법은 아래 글에서 별도로 정리하였습니다.
QGIS 래스터 계산기를 이용하여 SPOT-VGT DN 값을 Real NDVI로 변환하기
http://blog.daum.net/geoscience/681
이번 글에서는 10일 단위로 합성된 SPOT-VGT S10 DN값을 각각 Real NDVI로 변환하고,
이것을 중첩되는 셀의 최대값으로 합성하여 월 최대 식생지수(Monthly Maximum NDVI)로 합성해 보도록 하겠습니다.
여기서는 2012년 1월에 합성된 3장의 영상을 사용해 보도록 하겠습니다.
먼저 '레스터 레이어 추가' 버튼을 클릭하고,
아래와 같이 3장의 영상을 추가합니다.
이제 이 3장의 영상을 합성해 보도록 하겠습니다.
여기서는 GRASS commands의 r.mapcalculator를 이용해 계산해보도록 하겠습니다.
r.mapcalculator 실행 창은 아래와 같습니다.
Raster layer A, B, C를 각각 영상으로 지정하고 Formula를 입력하면 되겠습니다.
r.mapcalculator로 이용 가능한 함수는 아래와 같습니다.
여기서 type F는 결과값이 부동 소수점(floating point)으로, I는 정수(integer)로,
*는 정의된 수식 중 인수 하나가 부동 소수점일 경우에는 부동 소수점으로, 인수가 모두 정수일 경우에는 정수로 계산됩니다.
여기서는 아래 수식을 적용하기 위해 max 함수를 사용하였습니다.
수식을 아래와 같이 지정하고 결과값을 확인해 보면,
아래와 같이 월 최대 식생지수(MNDVI)가 계산되었습니다.
Greens 색상표를 이용하여 스타일을 아래와 같이 변경해 본 결과입니다. 어렵지 않죠?!
