R Under development (unstable) (2024-11-03 r87286 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2024 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> suppressPackageStartupMessages(library(sf))
> round_trip = function(x, EWKB = FALSE, pureR = FALSE) {
+ 	if (inherits(x, "sfg"))
+ 		x = st_sfc(x)
+ 	wkb = st_as_binary(x, EWKB = EWKB, pureR = pureR)
+ 	class(wkb) = "WKB"
+ 	# print(wkb)
+ 	y = st_as_sfc(wkb, EWKB = EWKB, pureR = pureR)
+ 	a = all.equal(x, y)
+ 	if (length(a) == 1 && is.logical(a) && a)
+ 		TRUE
+ 	else {
+ 		print(x)
+ 		print(wkb)
+ 		print(y)
+ 		FALSE
+ 	}
+ }
> 
> p3 = st_point(c(0,0,0))
> p3m = st_point(c(0,0,0), "XYM")
> p4 = st_point(c(0,0,0,0))
> p2 = st_point(c(0,0))
> ls = st_linestring(matrix(1:6,3))
> mp = st_multipoint(matrix(1:6,3))
> 
> outer = matrix(c(0,0,10,0,10,10,0,10,0,0),ncol=2, byrow=TRUE)
> hole1 = matrix(c(1,1,1,2,2,2,2,1,1,1),ncol=2, byrow=TRUE)
> hole2 = matrix(c(5,5,5,6,6,6,6,5,5,5),ncol=2, byrow=TRUE)
> pts = list(outer, hole1, hole2)
> pl1 = st_polygon(pts)
> 
> pol1 = list(outer, hole1, hole2)
> pol2 = list(outer + 12, hole1 + 12)
> pol3 = list(outer + 24)
> mp1 = st_multipolygon(list(pol1,pol2,pol3))
> 
> ml1 = st_multilinestring(list(outer, hole1, hole2))
> gc = st_geometrycollection(list(p2, ls, pl1, mp1))
> 
> sapply(list(p3, p3m, p4, p2, ls, mp, pl1, mp1, ml1, gc), round_trip, EWKB = FALSE)
 [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
> sapply(list(p3, p3m, p4, p2, ls, mp, pl1, mp1, ml1, gc), round_trip, EWKB = FALSE, pureR = TRUE)
 [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
> sapply(list(p3, p3m, p4, p2, ls, mp, pl1, mp1, ml1, gc), round_trip, EWKB = TRUE)
 [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
> sapply(list(p3, p3m, p4, p2, ls, mp, pl1, mp1, ml1, gc), round_trip, EWKB = TRUE, pureR = TRUE)
 [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
> 
> rawToHex(st_as_binary(st_multipoint(matrix(1:6,3))))
[1] "0104000000030000000101000000000000000000f03f0000000000001040010100000000000000000000400000000000001440010100000000000000000008400000000000001840"
> rawToHex(st_as_binary(st_sfc(st_point(c(0,1)), st_multipoint(matrix(1:6,3)))))
[1] "01010000000000000000000000000000000000f03f"                                                                                                      
[2] "0104000000030000000101000000000000000000f03f0000000000001040010100000000000000000000400000000000001440010100000000000000000008400000000000001840"
> try(rawToHex("error"))
Error in rawToHex("error") : 
  not implemented for objects of class character
> 
> # debug roundtrips sf -> GDAL -> sf; 
> # the first WKT is what GDAL reports, and will lack M
> st_as_text(st_sfc(sf:::CPL_roundtrip(st_sfc(st_linestring(matrix(1:18,6,3),dim="XYZ")))))
LINESTRING (1 7 13,2 8 14,3 9 15,4 10 16,5 11 17,6 12 18)
[1] "LINESTRING Z (1 7 13, 2 8 14, 3 9 15, 4 10 16, 5 11 17, 6 12 18)"
> st_as_text(st_sfc(sf:::CPL_roundtrip(st_sfc(st_multipoint(matrix(1:18,6,3),dim="XYZ")))))
MULTIPOINT (1 7 13,2 8 14,3 9 15,4 10 16,5 11 17,6 12 18)
[1] "MULTIPOINT Z ((1 7 13), (2 8 14), (3 9 15), (4 10 16), (5 11 17), (6 12 18))"
> st_as_text(st_sfc(sf:::CPL_roundtrip(st_sfc(st_point(c(0,0,0), dim="XYZ")))))
POINT (0 0 0)
[1] "POINT Z (0 0 0)"
> 
> if (sf:::CPL_gdal_version() >= "2.1.0") { # address GDAL/Fedora (gdal 2.0.2) error:
+   st_as_text(st_sfc(sf:::CPL_roundtrip(st_sfc(st_linestring(matrix(1:18,6,3),dim="XYM")))))
+   st_as_text(st_sfc(sf:::CPL_roundtrip(st_sfc(st_multipoint(matrix(1:18,6,3),dim="XYM")))))
+   st_as_text(st_sfc(sf:::CPL_roundtrip(st_sfc(st_point(c(0,0,0), dim="XYM")))))
+ } else {
+ 	"(output expected when gdal <= 2.1.0, e.g. CRAN/fedora)"
+ }
LINESTRING (1 7,2 8,3 9,4 10,5 11,6 12)
MULTIPOINT (1 7,2 8,3 9,4 10,5 11,6 12)
POINT (0 0)
[1] "POINT M (0 0 0)"
> 
> proc.time()
   user  system elapsed 
   0.57    0.10    0.67