QGIS에서 다중 래스터를 클리핑하고 밴드 셋(band set)을 생성하는 과정을 학습해 보도록 하겠습니다.
여기서는 소백산국립공원 경계에 5km 버퍼를 적용한 후, 해당 지역을 촬영한 센티널 2호(Sentinel-2A)
위성영상일 일괄적으로 잘라낸 후에, 결과 래스터들을 하나의 밴드 셋(band set)으로 생성해 보겠습니다.
먼저, 아래 그림과 같이 소백산국립공원 경계 파일을 추가합니다.
'벡터 > 공간 연산 도구 > 버퍼(들)'을 클릭하고,
버퍼 거리를 5000m로 설정하고 출력 Shape 파일을 지정합니다.
아래와 같이 버퍼 결과가 생성되었습니다.
이번에는 '벡터 > 조사 도구 > 레이어 범위로 폴리곤 만들기'를 클릭하고,
버퍼 결과를 입력 레이어로 선택한 후, 출력 폴리곤 Shape 파일도 아래와 같이 지정해 줍니다.
이렇게 센티널 2호 위성영상을 잘라낼 연구지역 폴리곤이 생성되었습니다. 이제 위성영상을 추가해볼까요?!
아래 표는 센티널 2호 위성영상의 밴드 목록표인데요, 1번, 9번, 10번 60m급
해상도 3개 밴드를 제외하고, 나머지 10개 밴드를 연구지역으로 잘라보겠습니다.
자, 일단 센티널 2호 위성영상 10개 밴드를 추가합니다.
여기서 문제가 하나 발견되었는데요, 센티널 2호 위성영상 영역이 연구지역 전체를 커버하지는 않습니다.
즉, 밴드들을 잘라낼 때 최대 y 좌표는 현재 연구지역 ymax과 센티널 영상 ymax값 차이만큼 보정해야겠습니다.
이제 영상을 자르기 위해 'SCP > Pre processing > Clip multiple rasters'를 클릭합니다.
Refresh list를 클릭하여 현재 Layers Panel에 등록된 밴드들을 추가하고 Select all을 클릭합니다.
다음으로 Clip coordinates에서 UL(upper left) 값과 LR(lower right) 값을 조정해 주면 되겠는데요,
Clip raster by extent에서 제공하는 클리핑 영역 값을 참고하여 UL과 LR 값을 채우면 되겠습니다.
xmin, xmas, ymin 값들은 연구지역 폴리곤의 것을 사용하고 ymax는 위성영상 값을 사용하면 되겠죠?!
아래와 같이 설정이 모두 끝났으면 Clip selected rasters 버튼을 클릭하여 파일들을 잘라내면 되겠습니다.
아래와 같이 clip가 접두사가 되어 클리핑 결과 래스터들이 생성되었습니다.
이제 'SCP > Band set'을 클릭하고,
아래와 같이 Band set definition에서 밴드 순서와 중앙 파장(Central Wavelength) 값까지 입력한 후,
Create raster of band set(stack bands)를 클릭하여 밴드 셋 생성을 마무리하면 되겠습니다. 간단하죠?!
아래와 같이 밴드 셋을 RGB 합성(composites)으로 처리한 결과입니다.
영상 일부를 확대해 본 결과인데요, SCP Toolbar를 이용하면 RGB 합성을 보다 쉽게 적용하실 수 있습니다.
8번(NIR), 4번(Red), 3번(Green) 밴드를 RGB 합성한 결과입니다. 이제 ROI 생성에 관해 살펴보도록 하겠습니다.
