set z_range and m_range in st_as_sf.data.frame()

fixes geoarrow/geoarrow-r#75

Author: edzer

NEWS.md

@@ -1,3 +1,7 @@
+# version 1.0-23
+
+* `st_as_sf.data.frame()` sets `z_range` and `m_range`; https://github.com/geoarrow/geoarrow-r/issues/75
+
 # version 1.0-22
 
 * `st_combine()` on `POINT` geometries ignores empty points; #2551

R/sf.R

@@ -60,20 +60,24 @@ st_as_sf.data.frame = function(x, ..., agr = NA_agr_, coords, wkt,
 				do.call(c, apply(cc, 1, fn))
 			} else { # points:
 				structure(points_rcpp(cc, dim),
-				n_empty = 0L, precision = 0, crs = NA_crs_,
-				bbox = structure(
-					c(xmin = min(cc[,1], na.rm = TRUE),
-					ymin =   min(cc[,2], na.rm = TRUE),
-					xmax =   max(cc[,1], na.rm = TRUE),
-					ymax =   max(cc[,2], na.rm = TRUE)), class = "bbox"),
-				class =  c("sfc_POINT", "sfc" ), names = NULL)
+					n_empty = 0L, precision = 0, crs = NA_crs_,
+					bbox = structure(
+						c(xmin = min(cc[,1], na.rm = TRUE),
+						ymin =   min(cc[,2], na.rm = TRUE),
+						xmax =   max(cc[,1], na.rm = TRUE),
+						ymax =   max(cc[,2], na.rm = TRUE)), class = "bbox"),
+					class =  c("sfc_POINT", "sfc" ), names = NULL)
 			}
 
 		if (remove) {
 			if (is.character(coords))
 				coords = match(coords, names(x))
 			x = x[-coords]
 		}
+		if (grepl("Z", dim))
+			attr(x[[sf_column_name]], "z_range") = compute_z_range(x[[sf_column_name]])
+		if (grepl("M", dim))
+			attr(x[[sf_column_name]], "m_range") = compute_m_range(x[[sf_column_name]])
 
 	}
 	st_sf(x, ..., agr = agr, sf_column_name = sf_column_name)