본문 바로가기
Python Tips

Proj Transformer 좌표 변환

by suminhan 2020. 4. 8.
from pyproj import Proj, transform

inProj = Proj('+proj=tmerc +lat_0=38 +lon_0=127.0028902777778 +k=1 +x_0=200000 +y_0=500000 +ellps=bessel +units=m +no_defs +towgs84=-115.80,474.99,674.11,1.16,-2.31,-1.63,6.43')
outProj = Proj(init='epsg:4326')
x1,y1 = yldf.iloc[0].x, yldf.iloc[0].y
x2,y2 = transform(inProj,outProj,x1,y1)
print (y2,',',x2)

간단하게는 이렇게 바꾸거나

 

from pyproj import Transformer
transformer = Transformer.from_proj(inProj, outProj)
points = []
for index, item in yldf[['x', 'y']].iterrows():
    #x1,y1 = item.x, item.y
    points.append((item.x, item.y))
    
latlist, lnglist = [], []
for pt in transformer.itransform(points):
    latlist.append(pt[1])
    lnglist.append(pt[0])
    
yldf['lat'] = latlist
yldf['lng'] = lnglist

이런식으로 위도 경도를 한꺼번에 변환 가능하다.

'Python Tips' 카테고리의 다른 글

CountVectorizer + pyLDAvis (N-gram LDA)  (0) 2020.06.10
Python MapBox 지도 visualize  (0) 2020.05.11
주로 쓰는 좌표계 정리  (0) 2020.04.08
GeoPandas 에서 GeoWithin Query하기  (0) 2020.03.17
dictionary 값으로 정렬하기  (0) 2020.03.17

댓글