Loop 0 is not valid: Edge 4441 crosses edge 4445

[PanfengZhang]

library(sf)
cmap <- read_sf("https://geo.datav.aliyun.com/areas_v3/bound/100000.json")
cmap2 <- st_crop(cmap, xmin = 72, ymin = 17, xmax = 136, ymax = 55)

When I run the above code, the error is as follows：

Error in s2_geography_from_wkb(x, oriented = oriented, check = check) :
Evaluation error: Found 1 feature with invalid spherical geometry.
[1] Loop 0 is not valid: Edge 4441 crosses edge 4445.

How to resolve the problem?

[rsbivand]

The input geometry is mis-specified:

> library(sf)
Linking to GEOS 3.10.2, GDAL 3.4.1, PROJ 8.2.1; sf_use_s2() is TRUE
> cmap <- read_sf("https://geo.datav.aliyun.com/areas_v3/bound/100000.json")
> cmap2 <- st_crop(cmap, xmin = 72, ymin = 17, xmax = 136, ymax = 55)
Error in s2_geography_from_wkb(x, oriented = oriented, check = check) : 
  Evaluation error: Found 1 feature with invalid spherical geometry.
[1] Loop 0 is not valid: Edge 4441 crosses edge 4445.
> st_is_longlat(cmap)
[1] TRUE
> sf_use_s2()
[1] TRUE
> sf_use_s2(FALSE)
Spherical geometry (s2) switched off
> cmap2 <- st_crop(cmap, xmin = 72, ymin = 17, xmax = 136, ymax = 55)
although coordinates are longitude/latitude, st_intersection assumes that they are planar
Error in CPL_geos_op2(op, x, y) : 
  Evaluation error: TopologyException: Input geom 0 is invalid: Self-intersection at 114.320037 22.381126999999999.
> sf_use_s2(TRUE)
Spherical geometry (s2) switched on
> cmap1 <- st_make_valid(cmap)
> cmap2 <- st_crop(cmap1, xmin = 72, ymin = 17, xmax = 136, ymax = 55)
Error in s2_geography_from_wkb(x, oriented = oriented, check = check) : 
  Evaluation error: Found 1 feature with invalid spherical geometry.
[1] Loop 136 is not valid: Edge 744 crosses edge 748.
> sf_use_s2(FALSE)
Spherical geometry (s2) switched off
> cmap2 <- st_crop(cmap1, xmin = 72, ymin = 17, xmax = 136, ymax = 55)
although coordinates are longitude/latitude, st_intersection assumes that they are planar
Error in CPL_geos_op2(op, x, y) : 
  Evaluation error: TopologyException: Input geom 0 is invalid: Self-intersection at 119.60970386388662 26.150518787877377.
> sf_use_s2(TRUE)
Spherical geometry (s2) switched on
> sf_use_s2(FALSE)
Spherical geometry (s2) switched off
> cmap1 <- st_make_valid(cmap)
> cmap2 <- st_crop(cmap1, xmin = 72, ymin = 17, xmax = 136, ymax = 55)
although coordinates are longitude/latitude, st_intersection assumes that they are planar
Warning message:
attribute variables are assumed to be spatially constant throughout all geometries 
> sf_use_s2(TRUE)

The only remedy is to turn off spherical geometry, rebuild the geometry with st_make_valid() as planar geometry, then crop as planar geometry, even when the geometries are really spherical, and the bounding box should be curvilinear, not rectangular. See: https://r-spatial.org/r/2020/06/17/s2.html and #1771.

[PanfengZhang]

Thank you very much for helping me solve this problem