我有一个数据帧,其中包含包含x坐标和y坐标的列。我想将它们转换为各自的纬度和经度值。我一直在阅读,但找不到一组可以使用的公式。
我正在使用python。我正在使用圣安东尼奥市的 311 呼叫开放数据集。
有人可以指导我如何将 x、y 坐标转换为纬度和经度吗?
任何帮助将不胜感激。
def xy_to_lonlat(x, y):
proj_latlon = pyproj.Proj(proj='latlong',datum='WGS84')
proj_xy = pyproj.Proj(proj="utm", zone=33, datum='WGS84')
lonlat = pyproj.transform(proj_xy, proj_latlon, x, y)
return lonlat[0], lonlat[1]
此函数应完成工作。 如果要执行反向转换,请使用以下命令:
def lonlat_to_xy(lon, lat):
proj_latlon = pyproj.Proj(proj='latlong',datum='WGS84')
proj_xy = pyproj.Proj(proj="utm", zone=33, datum='WGS84')
xy = pyproj.transform(proj_latlon, proj_xy, lon, lat)
return xy[0], xy[1]