R空间：具有相同CRS的2个shapefile之间的纬度/经度不匹配

如何解决R空间：具有相同CRS的2个shapefile之间的纬度/经度不匹配

我的目标是计算非洲次国家区域内极端天气事件的频率。为此，我设置了一个包含非洲省份的 shapefile，主要使用 GADM 数据和新的地理编码 EMDAT GDIS dataset 用于天气事件的点数据。

这是区域形状文件的样子：

    library(sf)

    st_geometry(africa_map)
    Geometry set for 796 features 
    Geometry type: GEOMETRY
    Dimension:     XY
    Bounding box:  xmin: -25.36042 ymin: -46.96575 xmax: 63.49391 ymax: 37.3452
    Geodetic CRS:  WGS 84
    First 5 geometries:
    MULTIPOLYGON (((-4.821226 24.99475,-4.821355 2...
    MULTIPOLYGON (((1.853562 35.8605,1.8424 35.865...
    MULTIPOLYGON (((-1.361976 35.3199,-1.358957 35...
    MULTIPOLYGON (((2.984874 36.81497,3.014171 36....
    MULTIPOLYGON (((7.262677 37.076,7.266449 37.07..

以及将经纬度转换为WGS 84后的GDIS数据集：

    gdis_africa_sf <- st_as_sf(x = gdis_africa,coords = c("longitude","latitude"),crs = "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0")

    st_geometry(gdis_africa_sf)
    Geometry set for 5171 features 
    Geometry type: POINT
    Dimension:     XY
    Bounding box:  xmin: -34.04233 ymin: -25.19619 xmax: 37.08849 ymax: 63.4228
    CRS:           +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0
    First 5 geometries:
    POINT (-17.09348 15.66576)
    POINT (-16.53153 15.77399)
    POINT (-16.20006 15.84419)
    POINT (-17.09348 15.66576)
    POINT (-16.53153 15.77399)

到现在为止，您已经可以看出有些不对劲，因为边界框根本不对应，即使投影看起来很合适。

    st_crs(africa_map)==st_crs(gdis_africa_sf)
[1] TRUE

将两者并排绘制时，无论我是使用新的 shapefile 还是仅应用数据框的经度和纬度，问题都会变得更加清晰。

    ggplot() +
      geom_sf(data = africa_map) +
      geom_sf(data = gdis_africa_sf)

Plot 1

    ggplot(data = africa_map) +
      geom_sf() +
      geom_point(data = gdis_africa,aes(x = longitude,y = latitude),color = "red",alpha = 0.3,size = 2,shape = 1) 

Plot 2

天气事件坐标似乎向西北移动了数千公里 - 但来源是什么？我该如何解决这个问题并使我的两个地理数据兼容？任何提示将不胜感激。

解决方法

您似乎误诊了问题，这不是翻译，而是对角线上的反映 - 您在站点点中交换了纬度和经度。试试：

 gdis_africa_sf <- st_as_sf(x = gdis_africa,coords = c("latitude","longitude"),crs = "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0")