fixes #2320

Author: edzer

NEWS.md

@@ -1,5 +1,7 @@
 # version 1.0-16
 
+* `[.sfc` works when setting argument `op`; #2320
+
 * `st_sample` for polygons is sensitive to setting `oriented = TRUE` to prevent wrongly correcting ring directions; #2308
 
 * add support for the GDAL `footprint` utility (requiring GDAL >= 3.8.0) to `gdal_utils`; #2305, by @goergen95

R/sfc.R

@@ -143,8 +143,11 @@ st_sfc = function(..., crs = NA_crs_, precision = 0.0, check_ring_dir = FALSE, d
 "[.sfc" = function(x, i, j, ..., op = st_intersects) {
 	if (!missing(i) && (inherits(i, "sf") || inherits(i, "sfc") || inherits(i, "sfg")))
 		i = lengths(op(x, i, ...)) != 0
-	st_sfc(NextMethod(), crs = st_crs(x), precision = st_precision(x),
-		dim = if(length(x)) class(x[[1]])[1] else "XY")
+	precision = st_precision(x)
+	crs = st_crs(x)
+	dim = if (length(x)) class(x[[1]])[1] else "XY"
+	x = unclass(x)[i] # now a list
+	st_sfc(x, crs = crs, precision = precision, dim = dim)
 }