[R-sig-Geo] Error in MtrxSet ... polygons not (all) closed

[Micha Silver]

Searching for the error in the subject line returns some discussions 
from a few years ago. But I'm not able to overcome this error.

Here's the polygon in question:

 > dput(aoi)

structure(list(x =

structure(list(structure(list(structure(c(11.0127043, 11.0127061, 
11.0134526, 11.0134554, 11.0134626, 11.0134692, 11.0134776, 11.0134813, 
11.0157343, 11.0157373, 11.0157437, 11.0157495, 11.0157546, 11.0157585, 
11.015761, 11.0157617, 11.0157616, 11.0157613, 11.0153992, 11.0153983, 
11.0153956, 11.0153922, 11.0153869, 11.0153804, 11.0153729, 11.0153694, 
11.0127478, 11.0127447, 11.0127379, 11.012732, 11.0127269, 11.0127232, 
11.0127212, 11.01272, 11.0127197, 11.0126956, 11.0126958, 11.0126965, 
11.0126978, 11.0127004, 11.0127043, 46.8484069, 46.8484048, 46.8476584, 
46.847656, 46.8476511, 46.8476481, 46.8476461, 46.8476456, 46.8475213, 
46.8475214, 46.847522, 46.8475235, 46.8475266, 46.8475311, 46.8475365, 
46.8475425, 46.8475489, 46.8475519, 46.8495623, 46.8495656, 46.8495729, 
46.8495788, 46.8495833, 46.8495856, 46.849587, 46.8495874, 46.8496081, 
46.8496079, 46.8496069, 46.849605, 46.8496014, 46.8495965, 46.8495906, 
46.8495839, 46.8495808, 46.8484309, 46.8484282, 46.8484218, 46.8484166, 
46.848412, 46.8484069), dim = c(41L, 2L), dimnames = list( NULL, c("X", 
"Y")))), class = c("XY", "POLYGON", "sfg"))), class = c("sfc_POLYGON", 
"sfc"), precision = 0, bbox = structure(c(xmin = 11.0126956, ymin = 
46.8475213, xmax = 11.0157617, ymax = 46.8496081), class = "bbox"), crs 
= structure(list( input = "EPSG:4326", wkt = "GEOGCRS[\"WGS 84\",\n 
ENSEMBLE[\"World Geodetic System 1984 ensemble\",\n MEMBER[\"World 
Geodetic System 1984 (Transit)\"],\n MEMBER[\"World Geodetic System 1984 
(G730)\"],\n MEMBER[\"World Geodetic System 1984 (G873)\"],\n 
MEMBER[\"World Geodetic System 1984 (G1150)\"],\n MEMBER[\"World 
Geodetic System 1984 (G1674)\"],\n MEMBER[\"World Geodetic System 1984 
(G1762)\"],\n MEMBER[\"World Geodetic System 1984 (G2139)\"],\n 
ELLIPSOID[\"WGS 84\",6378137,298.257223563,\n 
LENGTHUNIT[\"metre\",1]],\n ENSEMBLEACCURACY[2.0]],\n 
PRIMEM[\"Greenwich\",0,\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n 
CS[ellipsoidal,2],\n AXIS[\"geodetic latitude (Lat)\",north,\n 
ORDER[1],\n ANGLEUNIT[\"degree\",0.0174532925199433]],\n AXIS[\"geodetic 
longitude (Lon)\",east,\n ORDER[2],\n 
ANGLEUNIT[\"degree\",0.0174532925199433]],\n USAGE[\n SCOPE[\"Horizontal 
component of 3D system.\"],\n AREA[\"World.\"],\n 
BBOX[-90,-180,90,180]],\n ID[\"EPSG\",4326]]"), class = "crs"), n_empty 
= 0L)), row.names = 1L, class = c("sf", "data.frame"), sf_column = "x", 
agr = structure(integer(0), class = "factor", levels = c("constant", 
"aggregate", "identity"), names = character(0)))

This polygon is used as area of interest to crop Copernicus imagery that 
I am downloading using the {CDSE} package. (I have run this workflow 
successfully many times with several other aoi polygons)

I have tried:

1- sf::st_make_valid()

2- transforming to a UTM CRS

3- extracting coordinates and recreating the polygon (i.e.:
aoi_p <- st_polygon(list(st_coordinates(aoi)[,1:2]))
     aoi3 <- aoi_p |>
       st_sfc() |>
       st_as_sf()
     st_crs(aoi3) <- "EPSG:4326"

4- buffering by a small amount


The above error recurs in all cases (only with this problem polygon). 
Any suggestions?


Thanks


Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918
https://orcid.org/0000-0002-1128-1325

-- 
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918