Could you please provide an example of code that actually raises this error?
Best regards,
Zivan
On 19 Aug 2024, at 16:23, Micha Silver <tsvi...@gmail.com> wrote:
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
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
