폴리곤내부의 점 분석

참고

In QGIS3, you can perform this analysis using the Join attributes by location (summary) tool in the Processing Toolbox. Use the polygon layer as the Input Layer, the point layer as the Join layer and count as the Summaries to calculate. See 공간 결합 수행 (QGIS3) tutorial for step-by-step instructions.

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

Data Sources: [NCEI] [NATURALEARTH]

과정

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

../_images/1202.png

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

참고

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

../_images/2166.png

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

../_images/3105.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/468.png

  1. 메뉴 벡터 –> 분석 도구 –> 폴리곤의 점 :menuselection:`Vector –> Analysis Tools –> Point in Polygon`을 클릭합니다.

../_images/562.png

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

참고

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

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

../_images/758.png

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

../_images/856.png

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

../_images/955.png

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

../_images/1063.png

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

../_images/11102.png

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

If you want to give feedback or share your experience with this tutorial, please comment below. (requires GitHub account)