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 |
댓글