（多）多边形的反向坐标

Question

我正在处理多个国家/地区的分区的 shapefile，对于一个国家/地区（冰岛），X 和 Y 坐标似乎在 shapefile 中交换了位置。

数据可以在这里下载：shapefile data； IS_50V:mork_kjordaemi 是相关数据集，select 下载下拉菜单中的“shape-zip”选项。我一直在使用 R 中的“sf”包处理所有 shapefile 工作，它与我拥有的所有其他 shapefile 数据完美配合。

library(sf)

ic_2003 <- downloaded_data

st_crs(ic_2003) 给我

Coordinate Reference System:
  User input: ISN2016 
  wkt:
GEOGCRS["ISN2016",
    DATUM["Islands Net 2016",
        ELLIPSOID["GRS 1980",6378137,298.257222101,
            LENGTHUNIT["metre",1]]],
    PRIMEM["Greenwich",0,
        ANGLEUNIT["degree",0.0174532925199433]],
    CS[ellipsoidal,2],
        AXIS["geodetic latitude (Lat)",north,
            ORDER[1],
            ANGLEUNIT["degree",0.0174532925199433]],
        AXIS["geodetic longitude (Lon)",east,
            ORDER[2],
            ANGLEUNIT["degree",0.0174532925199433]],
    USAGE[
        SCOPE["unknown"],
        AREA["Iceland"],
        BBOX[59.96,-30.87,69.59,-5.55]],
    ID["EPSG",8086]]

head(ic_2003) 给我

Simple feature collection with 6 features and 15 fields
Geometry type: MULTIPOLYGON
Dimension:     XY
Bounding box:  xmin: 63.29577 ymin: -24.53268 xmax: 66.56644 ymax: -13.49462
Geodetic CRS:  ISN2016

我试过 ic_2003 <- st_transform(ic_2003, 4326) 但这并不能解决问题。

我也试过 ic_2003 <- st_transform(ic_2003, pipeline = "+proj=pipeline +step +proj=axisswap +order=2,1")，就像 here 那样，但这也没有解决问题。

如果我绘制数据

ggplot(ic_2003) +
 geom_sf() +
 coord_sf()

我得到了正确的形状，但旋转了 90 度并且在世界地图上的位置错误。

如果您能给我任何帮助，我们将不胜感激。

Answer 1

肯定有一种 sf 方法可以轻松做到这一点，但您也可以使用 purrr::modify（其工作原理类似于 map）来交换所有几何图形 lat/lon 列（列表中的列表中的列表中的矩阵）而不更改 sf 属性...

library(sf)
library(tidyverse)

ic_2003 <- st_read("mork_kjordaemiPolygon.shp") #from link above

ic_2003 <- ic_2003 %>% 
      mutate(geometry = modify(geometry, modify, ~list(.[[1]][,c(2,1)]))) 

ggplot(ic_2003) +
  geom_sf() +
  coord_sf()

（多）多边形的反向坐标

Reverse coordinates for (multi-) polygon

r

axes

shapefile

sf