폴리곤내부의 점 분석

GIS의 힘은 여러 데이터 소스를 한꺼번에 같이 분석하는데 있습니다. 종종 찾고있는 답은 수많은 다른 레이어에 있고 이런 정보를 추출하고 컴파일하기 위한 분석을 할 필요가 있습니다. 그러한 유형의 분석이 **Points-in-Polygon**입니다. 폴리곤 레이어와 포인트 레이어를 가지고 있을때 그리고 얼만나 많이 혹은 어떤 저들이 각 폴리곤의 경계내에 있는지를 알고자 할때 이 분석방법을 사용할 수 있습니다.

과업 개요

모두에게 알려져 있는 주요 지진발생위치 데이터가 주어지고, 어느 나라가 가장 높은 빈도의 지진이 일어났는지 알아내려고 시도합니다.

데이터 획득

모든 주요한 지진을 나타내는 레이어로써 NOAA’s National Geophysical Data Center의 `Significant Earthquake Database <http://www.ngdc.noaa.gov/nndc/struts/form?t=101650&s=1&d=1>`_를 사용합니다. `tab-delimited earthquake data <http://www.ngdc.noaa.gov/nndc/struts/results?type_0=Exact&query_0=$ID&t=101650&s=13&d=189&dfn=signif.txt>`_를 다운로드 합니다.

Natural Earth는 Admin 0 - Countries 데이터셋을 가지고 있습니다. `countries <http://www.naturalearthdata.com/http//www.naturalearthdata.com/download/10m/cultural/ne_10m_admin_0_countries.zip>`_를 다운로드 합니다.

For convenience, you may directly download a copy of the dataset from the link below:

signif.txt

ne_10m_admin_0_countries.zip

데이터 출저: [NGDC] [NATURALEARTH]

과정

  1. 메뉴 레이어 –> 구분자로 분리된 텍스트 레이어 추가 :menuselection:`Layer –> Add Delimited Text Layer`를 열고 다운로드한 ``signif.txt``파일을 찾습니다.

../_images/1106.png
  1. 이것은 탭으로 구분된 파일 *tab-delimited file*이므로 파일 포맷 :guilabel:`File format`에 탭 :guilabel:`Tab`을 선택합니다. X필드 :guilabel:`X field`와 Y필드 :guilabel:`Y field`는 자동적으로 나타납니다. :guilabel:`OK`를 클릭합니다.

주석

파일은 QGIS로 불러오기를 시도하면서 에러메세지를 를 보게됩니다. 이것들은 유효한 에러이고 파일의 몇 행은 불러들여지지 않을 것입니다. 이 예제에서는 에러를 무시할 수 있습니다.

../_images/258.png
  1. 지진 데이터셋은 위도/경도 좌표를 가지므로 좌표계선택 Coordinate Reference System Selector 다이알로그에서 좌표계로 :guilabel:`WGS 84 EPSG:436`를 선택합니다.

../_images/334.png
  1. The earthquake point layer would now be loaded and displayed in QGIS. Let’s also open the Countries layer. Go to Layer ‣ Add Vector Layer. Browse to the downloaded ne_10m_admin_0_countries.zip file and click Open. Select the ne_10m_admin_0_countries.shp as the layer in the Select layers to add... dialog.
../_images/424.png
  1. 메뉴 벡터 –> 분석 도구 –> 폴리곤의 점 :menuselection:`Vector –> Analysis Tools –> Point in Polygon`을 클릭합니다.

../_images/525.png
  1. 팝업창에서 폴리곤 레이어와 포인트 레이어를 각각 선택합니다. 출력 레이어에 ``earthquake_per_coutry.shp``를 입력하고 :guilabel:`OK`를 클릭합니다.

주석

OK를 누른 후 기다립니다. QGIS는 결과를 계산하기까지 최대 10분이 소요될 수 있습니다.

  1. 레이어를 TOC에 추가할 것을 물어보면 :guilabel:`Yes`를 클릭하십시오.

../_images/722.png
  1. TOC에 새로운 렝어가 추가되는 것을 보게될 것입니다. 레이어를 우측클릭하여 속성 테이블을 열고 속성테이블 열기 :guilabel:`Open Attribute Table`를 선택합니다.

../_images/821.png
  1. 속성 테이블에서 ``PNTCNT``라는 이름의 새로운 필드명을 주목하게 됩니다. 이것은 각 폴리곤내에 포함된 지진레이어에서 점의 갯수를 계산한 것입니다.

../_images/919.png
  1. 답을 얻기 위해 ``PNTCNT``필드를 정렬하면 가장 수가 높은 나라가 답이 됩니다. 내림차순으로 정렬하기 위해 ``PNTCNT``열을 2번 클릭합니다. 그것을 선택하기 위하여 첫번째 행을 클릭하고 속성 테이블을 닫습니다.

../_images/1020.png
  1. QGIS 창으로 되돌아가면 한 객체가 노란색으로 밝게 빛나는 것을 보게됩니다. 이것이 속성테이블에서 가장 점이 많아서 선택한 행과 연결된 객체입니다. 객체 확인 :guilabel:`Identify`툴을 선택하고 폴리곤을 클릭합니다. 유효한 지진빈도가 가장 높은 국가가 중국 **China**임을 볼 수 있습니다.

../_images/1123.png

중국이 가장 많은 주요 지진 횟수를 가지고 있다는 것이 2가지 데이터셋의 간단한 분석에서 명확해 졌습니다. 국가의 인구와 면적을 고려하여 분석을 더욱 개선하고 반대로 주요 지진의 영향을 가장 덜 받는 국가를 찾아 낼 수도 있습니다.

comments powered by Disqus

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