REMOTE SENSING

SM(Status Map, 상태지도)를 이용한 SPOT VEGETATION 영상 전처리

유병혁 2012. 10. 15. 04:15

SPOT VEGETATION S10 NDVI 값 영상을 다운로드 받으면, NDVI와 SM 두 개 영상이 포함되어 있습니다.

여기서 NDVI는 DN값으로 저장된 상태이고, SM은 Status Map(상태지도)로써 SPOT VGT 영상 처리에 유용한 정보들을 제공합니다.

 

SM 레이어에 대한 공식적인 설명은 아래와 같습니다.

 

 

 

그럼, SM을 이용하여 SPOT VEGETATION 영상을 전처리하는 방법을 살펴 보겠습니다.

처리 방법을 알아보기 위하여 2011년 7월 1일 촬영된 S-10 NDVI 영상을 다운 받았습니다(아래).

 

SM 영상은 아래와 같습니다.

 

참고로, SPOT VEGETATION 영상은 4개의 밴드를 가지고 있으며 파장(Wavelength)은 아래와 같습니다.

 

먼저, SM 밴드에서 RED 밴드의 방사 품질(radiometric quality)을 조회해 보도록 하겠습니다.

ENVI 프로그램에서 Basic Tools > Band Math를 실행합니다.

 

Bit NR 6에서 RED 밴드(B2)의 방사 품질이 양호(good)하면 값이 1, 불량(bad)하면 값이 0인 특성에 따라 수식은 아래와 같이 입력합니다.

(b1 and 2^6) ge 1

 

Variables to Bands Pairings에서 SM 밴드를 지정합니다.

 

ENVI Output to Memory Warning에서 출력 결과를 File 또는 Memery 중 선택합니다.

 

처리 결과는 아래와 같습니다.

 

 

같은 방식으로 NIR 밴드 방사 품질은 아래와 같이 계산합니다.

(b1 and 2^5) ge 1

 

Bit NR 3에서 육상(land)은 1, 해양(water)은 0으로 구분되며, 육상을 분리하기 위한 산식은 아래와 같습니다.

(b1 and 2^3) ge 1

 

날씨가 많은(clear) 경우는 Bit NR 1과 Nit NR 0에서 모두 0이므로, 아래 산식으로 계산됩니다.

((b1 and 2^1) eq 0) and ((b1 and 2^0) eq 0)

 

따라서 적색과 근적외 밴드의 방사 품질이 양호하고, 육상지역이며, 날씨가 많은 경우를 모두 종합하면 산식은 아래와 같습니다.

((b1 and 2^6) ge 1) and ((b1 and 2^5) ge 1) and ((b1 and 2^3) ge 1) and (((b1 and 2^1) eq 0) and ((b1 and 2^0) eq 0))

 

이렇게 생성된 매스크(Mask)자료를 가지고 NDVI 영상에서 분석 목적에 불필요한 화소를 제거해주면  되겠습니다.