我有一个数据帧,结果如下:
我想将这些坐标列从wgs84转换为lon&lat,最后在我的数据帧中添加这些新列:
对于转换,我使用了以下代码,但我认为应该有更好的方法,而不必将坐标列转换为列表并创建新的一个数据帧。
import pyproj as pp
from mpl_toolkits.basemap import Basemap
import pandas as pd
cx =dfb.COORDENADA_X.tolist()
cy =dfb.COORDENADA_Y.tolist()
utm15_wgs84 = pp.Proj(init='epsg:32615')
for ix, iy in zip(cx, cy):
lon, lat = utm15_wgs84(ix, iy, inverse=True)
print(lon, lat)
有什么建议吗?
最佳答案
使用apply
pandas
中的DataFrame
功能例如
dfb[['wgs_x', 'wgs_y']] = dfb.apply(lambda row:utm15_wgs84(row['COORDENADA_X'], row['COORDENADA_Y'], inverse=True), axis=1).apply(pd.Series)