0

不幸的是，我从爱尔兰横轴墨卡托 (ITM) 到 WGS84 经纬度的投影似乎出错了，因为绘制的坐标与来自 CSO的都柏林地图不一致（见下文）。

我的转换坐标绘制在都柏林地图上

转换后的数据来自爱尔兰估价办公室，ITM X 和 Y 坐标被输入一个函数，该函数改编自先前的 stackoverflow 讨论，该讨论使用 geopandas 的内置points_from_xy方法来转换坐标参考系统之间的坐标：

def create_geodf_from_GPS (df, latitude, longitude, crs):

    locations = gpd.points_from_xy(longitude, latitude)
    geo_df = gpd.GeoDataFrame(df, geometry=locations)
    geo_df.crs = crs

    return geo_df


VO_geo = create_geodf_from_GPS(VO, VO[" X ITM"], VO[" Y ITM"], crs = 'epsg:2157')
VO_geo = VO_geo.to_crs('epsg:4326')

有谁知道这里可能出了什么问题？

4

1 回答 1

0

感谢@joris，非常简单的修复

使用 x & y 作为 gpd.points_from_xy 的参数而不是先前混合的经度和纬度的更改函数：

def create_geodf_from_GPS (df, x, y, crs):

   locations = gpd.points_from_xy(x, y)
   geo_df = gpd.GeoDataFrame(df, geometry=locations)
   geo_df.crs = crs

return geo_df

现在在 WGS84 纬度-经度中绘制数据按预期工作：

VO_geo = create_geodf_from_GPS(VO, x=VO[" X ITM"], y=VO[" Y ITM"], crs = 'epsg:2157')
VO_geo.to_crs('epsg:4326').plot()

注意：必须通过使用 geopandas 的 (gpd) 空间连接功能过滤掉非都柏林数据来清理数据以去除明显的异常值

VO_geo_clean = gpd.sjoin(VO_geo.to_crs('epsg:4326'), map_of_Dublin)

结果： VO 数据绘制在都柏林地图上

于 2020-03-21T11:11:01.847 回答