GIS 위경도 변환

태평양돌고래·2021년 8월 5일
0

GIS

목록 보기
1/1

GIS

Geometry Information System의 약자로, "지리정보시스템" 이라고 굳이 번역할 수 있겠지만,
지리정보를 가지고 데이터를 분석하는 것을 GIS분야라고 통칭하는 듯 하다.

활용 예시

매출 분포


( 백화점을 어떻게 운영하고 의사결정 할지 그 다음 스텝이 보이는 GIS의 예 )

출처: http://www.biz-gis.com/index.php?mid=GIS_Essay&document_srl=5054

GIS에는 좌표계 라는 개념이 항상 따라오는데, 그것에 대해 궁금한 점이 생기거나, 잊어서 다시 찾는 블로그 링크는 여기


무언가 분석하려면 EPSG code가 4326일 경우가 가장 편리했다. 아직까지는..
하지만 우리나라 지리정보를 받을 경우에는 대부분 EPGS code는 5179인 경우가 많다.

따라서 EPSG:4326으로 변환하는 내용을 다시 구글링하지 않기 위해 여기에 정리해 둠.

# 파일 읽어오는 라이브러리 import
import pandas as pd
import numpy as np
import geopandas as gpd

내용에서 사용한 행정구역 경계선 좌표정보 다운로드는 여기 링크

df_region1 = gpd.read_file( '../TL_SCCO_CTPRVN.shp' , encoding = 'cp949')
# encoding = 'cp949'를 해주면 한글이 깨지지 않는다. 그래도 깨진다면 UTF-8로 시도. 그래도 깨진다면.. 엄청난 구글링을 해야 함 ㅠㅠ
df_region1.head()

결과

결과의 POLYGON은 사전적 의미로 다각형 이라는 뜻인데, GIS에서는 뚫려있지 않는 2차원 공간 지면이라고 생각하면 편하다.
여기의 POLYGON을 친숙한 위경도 좌표계로 바꾸고 싶다.

# 변환에 필요한 라이브러리 import
import pyproj
from fiona.crs import from_epsg

위의 라이브러리를 이용하여

df_region1_4 = df_region1.to_crs(epsg=4326)

해주면 끝이다. 예전에 구글링할때는 엄청 복잡하게 했던 것 같은데, 간편한 방법을 다시 찾아서 이렇게 정리해 둠...

df_region1_4.head()

결과

GIS를 다루다 보면 여기 저기서 데이터를 긁어모으는 경우가 많은데,
좌표계를 통일시키고 작업해야 한다.
그래서 이와 같이 EPSG code = 4326으로 해 두고 작업하는 것을 추천.
파이썬으로 GIS 정보 시각화 할 때 사용하는 folium 라이브러리에서는 EPSG code = 4326일 때만 인식하는 것 같다.

0개의 댓글