Error in MtrxSet ... polygons not (all) closed

Hi Micha,
Thank you for sharing the file.
I?ve noticed that this polygon has only 5 vertices, while the one in 
the previous message has 40, so this is not exactly the same object. 
Both seem OK, and I was able to download and crop satellite images 
with CDSE without problems.

Difficult to say more without a reproducible example raising the error.
Best,
Zivan


On 20 Aug 2024, at 09:25, Micha Silver <tsvibar at gmail.com> wrote:

Hi Zivan, Roger:
Thanks for the quick response.Here's the RDS of the polygon:
https://filebin.net/pdrazh0e4ej7el9f

and sessionInfo, etc:

sessionInfo()

R version 4.2.2 Patched (2022-11-10 r83330)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Debian GNU/Linux 12 (bookworm)

Matrix products: default
BLAS:?? /usr/lib/x86_64-linux-gnu/openblas-pthread/libblas.so 
<http://libblas.so/>.3
LAPACK: 
/usr/lib/x86_64-linux-gnu/openblas-pthread/libopenblasp-r0.3.21.so 
<http://libopenblasp-r0.3.21.so/>

locale:
?[1] LC_CTYPE=en_US.UTF-8?????? LC_NUMERIC=C LC_TIME=en_IL.UTF-8
?[4] LC_COLLATE=en_US.UTF-8???? LC_MONETARY=en_IL.UTF-8 
LC_MESSAGES=en_US.UTF-8
?[7] LC_PAPER=en_US.UTF-8?????? LC_NAME=C LC_ADDRESS=C
[10] LC_TELEPHONE=C???????????? LC_MEASUREMENT=en_US.UTF-8 
LC_IDENTIFICATION=C

attached base packages:
[1] stats???? graphics? grDevices utils???? datasets methods?? base

other attached packages:
?[1] tidyr_1.3.1?????? sf_1.0-16???????? dplyr_1.1.4 lubridate_1.9.3
?[5] RSQLite_2.3.7???? CDSE_0.2.0??????? rOPTRAM_0.3.0.000 ggplot2_3.5.1
?[9] knitr_1.47??????? terra_1.7-78????? devtools_2.4.5 usethis_2.2.3

loaded via a namespace (and not attached):
?[1] pkgload_1.3.4????? bit64_4.0.5 jsonlite_1.8.8???? here_1.0.1
?[5] shiny_1.8.1.1????? askpass_1.2.0 blob_1.2.4???????? remotes_2.5.0
?[9] sessioninfo_1.2.2? pillar_1.9.0 glue_1.7.0???????? digest_0.6.35
[13] promises_1.3.0???? colorspace_2.1-0 htmltools_0.5.8.1? httpuv_1.6.15
[17] pkgconfig_2.0.3??? s2_1.1.6 purrr_1.0.2??????? xtable_1.8-4
[21] scales_1.3.0?????? processx_3.8.4 later_1.3.2??????? openssl_2.2.0
[25] timechange_0.3.0?? tibble_3.2.1 proxy_0.4-27?????? generics_0.1.3
[29] ellipsis_0.3.2???? cachem_1.1.0 withr_3.0.0??????? cli_3.6.3
[33] lutz_0.3.2???????? magrittr_2.0.3 mime_0.12????????? memoise_2.0.1
[37] ps_1.7.6?????????? fs_1.6.4 fansi_1.0.6??????? class_7.3-22
[41] geojsonsf_2.0.3??? pkgbuild_1.4.4 httr2_1.0.1??????? profvis_0.3.8
[45] tools_4.2.2??????? lifecycle_1.0.4 stringr_1.5.1????? munsell_0.5.1
[49] callr_3.7.6??????? compiler_4.2.2 e1071_1.7-14?????? rlang_1.1.4
[53] classInt_0.4-10??? units_0.8-5 grid_4.2.2???????? rstudioapi_0.16.0
[57] rappdirs_0.3.3???? htmlwidgets_1.6.4 miniUI_0.1.1.1???? wk_0.9.2
[61] gtable_0.3.5?????? codetools_0.2-20 DBI_1.2.3????????? curl_5.2.1
[65] R6_2.5.1?????????? fastmap_1.2.0 bit_4.0.5????????? utf8_1.2.4
[69] rprojroot_2.0.4??? KernSmooth_2.23-24 desc_1.4.3???????? 
stringi_1.8.4
[73] Rcpp_1.0.12??????? vctrs_0.6.5 tidyselect_1.2.1?? xfun_0.45
[77] urlchecker_1.0.1

sf_extSoftVersion()

????????? GEOS?????????? GDAL???????? proj.4 GDAL_with_GEOS???? USE_PROJ_H
????? "3.11.1"??????? "3.6.2"??????? "9.1.1" "true"???????? "true"
????????? PROJ
?????? "9.1.1"

This error occurs, during a longer analysis, when I begin downloading 
imagery , and only with the above aoi polygon. So I'm trying to strip 
out all the steps in the code to identify what the cause is. No luck 
yet... Note: I don't think that the CDSE package is the cause.

Best regards, Micha


On 19/08/2024 21:29, Roger Bivand wrote:

Adding the output of sessionInfo() and sf_extSoftVersion() might also help. The difficulty Zivan points out is that the structure as read is valid, but before conversion to text by dput() may have actually had non-identical first and last coordinates, because of the number of digits used. Also try putting the file created by saveRDS(aoi, ...) on a website with a link here. It would be useful to see the code in the workflow at the point of failure (you point to Error in MtrxSet, maybe also run traceback() after the failure to give more precision in where this is occurring, maybe in stars?

Hop this helps,

Roger

--
Roger Bivand
Emeritus Professor
Norwegian School of Economics
Postboks 3490 Ytre Sandviken, 5045 Bergen, Norway
Roger.Bivand at nhh.no

-- 
Micha Silver
Remote Sensing Lab, Sde Boker
Ben Gurion University
+972-523-665918