래스터 붙이기(Mosaicing)와 자르기(Clipping)

이 예제에서는 QGIS에서 기본적인 래스터 다루기로 래스터 보기, 붙이기 그리고 부속 세팅(subsetting)을 알아봅니다.

작업 개요

이 예제에서는 브라질의 래스터데이터를 공공도메인으로 다운로드하고 그것을 QGIS에서 볼 것입니다. 다음으로 이러한 데이터를 합쳐 붙여서 하나로 만들고 경계 데이터를 이용하여 잘라내어서 끊김없는 단일 브라질 데이터셋을 얻을 것입니다.

기타 스킬

  • 공공 도메인의 실시간 인공위성 영상을 찾고 다운로드합니다.

  • 벡터 레이어로부터 한가지 피쳐를 선택하고 새로운 shapefile 로 저장합니다.

데이터 획득

래스터 데이터를 잘라낼 브라질 경계 데이터가 필요합니다. Natural Earth에서 Admin 0 - Countries shapefile을 얻을 수 있습니다.

`NASA/GSFC, Rapid Response site <https://earthdata.nasa.gov/data/near-real-time-data/rapid-response/modis-subsets>`_는 실시간 인공위성영상을 잘 모아둔 곳입니다. `USDA Foreign Agricultural Service (FAS) subsets. <http://lance-modis.eosdis.nasa.gov/imagery/subsets/?project=fas>`_는 우수한 지역적 결과물입니다. 이 예제에서는 브라질의 2km 해상도의 FAS subsets 을 사용합니다.

해당 데이터를 어떻게 찾고 다운로드하는 가는 다음과 같습니다.

  1. Open the South America region subsets. Find the Brazil subsets shown in the FAS Subsets section. Click on any one of them.
../_images/1109.png
  1. 상세 페이지에서 선택한 결과물 아래의 2km 링크를 클릭합니다. 여기서 :guilabel:`NDVI`결과물을 다운로드 할 것입니다. NDVI에 대해 더 알아봅시다. `Learn more about NDVI <http://en.wikipedia.org/wiki/Normalized_Difference_Vegetation_Index>`_

../_images/268.png
  1. 래스터 이미지를 다운로드 하기 위해 :guilabel:`Download GeoTIFF file`를 클릭합니다.

../_images/337.png

브라질의 7개 FAS subsets에 대해 과정을 반복합니다.

For convenience, you can directly download sample data used in this tutorial from links below.

ne_10m_admin_0_countries.zip

FAS_Brazil1.2013363.aqua.ndvi.2km.tif

FAS_Brazil2.2013363.terra.ndvi.2km.tif

FAS_Brazil3.2013363.aqua.ndvi.2km.tif

FAS_Brazil4.2013363.aqua.ndvi.2km.tif

FAS_Brazil5.2013363.aqua.ndvi.2km.tif

FAS_Brazil6.2013363.terra.ndvi.2km.tif

FAS_Brazil7.2013363.aqua.ndvi.2km.tif

데이터 출처: [LANCE] [NATURALEARTH]

과정

  1. QGIS를 열고 메뉴 레이어 –> 레이어 추가 –> 래스터 레이어 추가 :menuselection:`Layer –> Add Raster Layer..`로 갑니다.

../_images/427.png
  1. 개별 영상들의 목록을 찾습니다. 컨트롤 :kbd:`Ctrl`키를 누른 채로 이미지파일을 클릭하여 다중선택을 합니다. :guilabel:`Open`을 클릭합니다.

../_images/528.png
  1. 왼쪽 패널의 Table of Content에 이미지들이 올라온 것을 보게될 것입니다. 이제 개별 이미지들을 합쳐서 하나의 `Mosaic`이미지를 만듭시다. 메뉴에서 래스터 –> 기타 사항 –> 병합 :menuselection:`Raster –> Miscellaneous –> Merge`을 클릭합니다.

주석

QGIS의 메뉴에서 래스터 :guilabel:`Raster`는 :guilabel:`GdalTools`라는 핵심 플러그인에서 비롯된 것입니다. 만약 래스터 :guilabel:`Raster`가 보이지 않는다면 :guilabel:`GdalTools`플러그인을 플러그인 –> 플러그인 관리 및 설치 –> 설치 :menuselection:`Plugins –> Manage and install plugins –> Installed`로 가서 활성화 하십시오. 보다 자세한 사항은 :doc:`using_plugins`를 참조하기 바랍니다.

../_images/626.png
  1. 결합 Merge`다이알로그에서, 입력파일 :guilabel:`Input files 옆에 있는 선택 Select... 을 클릭하고 모든 geotiff파일이 포함되어 있는 디렉토리를 찾습니다. 컨트롤 Ctrl`키를 누른 채로 해당되는 모든 파일을 클릭합니다. 이제 출력파일  :guilabel:`Output file`에 있는 선택  :guilabel:`Select... 을 클릭하고 출력파일명을 ``Brazil_mosaic.tif``로 입력합니다. 밑에서 종료시 캔버스로 불러옴 Load into canvas when finished 옆에 있는 상자에 체크합니다. :guilabel:`OK`를 클릭합니다.

../_images/725.png
  1. You will get a pop-up message saying Processing complete, once the mosaic is created and loaded to the QGIS Canvas. You will see that the individual images and now combined and mosaiced into a single layer. You can now turn off individual layers by un-checking the box next to them.
../_images/824.png
  1. 여기서 할 수 있는 다른 래스터 작업으로 영상을서브셋을 하거나 ‘잘라내기’를 할 수 있습니다. 정확한 형태로 래스터를 잘라내기 위해서 벡터 레이어의 폴리곤을 사용할 수 있습니다. Natural Earth에서 다운로드한 국가 경계 shapefile을 불러옵니다. 메뉴 레이어 –> 벡터 레이어 추가 :menuselection:`Layer –> Add Vector Layer`로 갑니다.

../_images/922.png
  1. Select the ne_10m_admin_0_countries.zip file and click Open. When prompted to select the layer within the zip file, select ne_10m_admin_0_countries.shp.
../_images/1023.png
  1. 벡터레이어가 불러들여지면 브라질을 선택하여 폴리곤을 추출합니다. 툴바에서 단일클릭이나 영역으로 객체선택 :guilabel:`Select Single Feature`중 객체선택을 선택합니다.

../_images/1126.png
  1. 브라질 폴리곤 아무곳이나 클릭을 하면 선택이 됩니다.

../_images/1224.png
  1. ``ne_10m_admin_0_countries``레이어를 오른쪽 클릭하고 다른 이름으로 저장 :guilabel:`Save Selection As...`을 선택합니다.

../_images/1321.png
  1. 출력파일명으로 ``brazil_boundary.shp``를 입력하고 저장된 파일을 지도에 추가 :guilabel:`Add saved file to map`상자가 체크되어 있는지 확인합니다. :guilabel:`OK`를 클릭합니다.

../_images/1420.png
  1. QGIS에 브라질 경계 폴리곤이 불러들여집니다.

../_images/1519.png
  1. 이제 메뉴 래스터 –> 추출 –> 잘라내기 :menuselection:`Raster –> Extraction –> Clipper`로 갑니다.

../_images/1618.png
  1. 입력파일 input file (raster)``Brazil_mosaic``를 선택합니다. 출력 파일 Output file`에 ``Brazil_mosaic_clipped``라고 명칭을 입력합니다. :guilabel:`Clipping mode 섹션에서 :guilabel:`Mask layer`를 선택합니다. 마스크 레이어로써 새로 만들어진 ``brazil_boundary``를 선택합니다. 종료시 캔버스로 불러옴 :guilabel:`Load into canvas when finished`옆에 있는 상자를 체크합니다. :guilabel:`OK`를 클릭합니다.

../_images/1717.png
  1. 새롭게 자르기한 레이어가 QGIS로 불러들여집니다. 모자이크 주변의 검은색 픽셀을 보일 겁니다. ``Brazil_mosaic_clipped``레이어를 오른쪽 클릭하고 속성 :guilabel:`Properties`을 선택합니다.

../_images/1817.png
  1. 투명도 Transparency`탭으로 가서 :guilabel:`Additional no data value`에 `0 을 추가합니다.

../_images/1915.png
  1. 이제 국가경계로 잘려진 훌륭한 모자이크를 얻습니다. 이 모자이크는 배경이나 추가 분석시 프로젝트에서 사용할 수 있습니다.

../_images/2012.png
comments powered by Disqus

This work is licensed under a Creative Commons Attribution 4.0 International License