R 使用显式二元谓词子集化sfc对象时出错

这里出现问题是因为`[.sfc`在内部使用NextMethod().其工作方式是将显式传递给`[.sfc`的参数传递给列表的默认子设置方法.这是在C代码函数do_subset_dflt中实现的原语"["函数.

对于R函数来说，这是不寻常的，它忽略所有参数的名称(除了drop以外)，并且只按顺序执行:

my_list <- list(a = 1, b = 2)

my_list[the_name_of_this_argument_is_ignored = 2]
#> $b
#> [1] 2

my_list[the_name_of_this_argument_is_ignored = 2, drop = FALSE]
#> $b
#> [1] 2

这甚至导致了令人困惑的结果

my_mat <- matrix(1:4, 2, 2)

my_mat[i = 1, j = 2]
#> [1] 3

my_mat[j = 1, i = 2]
#> [1] 3

如果(非删除)参数的数量超过了对象的维度数量，我们就会得到您的错误:

my_list[2, op = sf::st_intersects]
#> Error in my_list[2, op = sf::st_intersects] : 
#>   incorrect number of dimensions

默认参数可以毫无问题地工作，这似乎令人困惑，但进行完全相同的函数调用但显式命名默认参数会引发错误.这是NextMethod工作方式的结果-它只传递您专门传递给初始调用的符号，而忽略任何默认参数.

通过创建一个具有自己的子集方法的新类，可以非常容易地显示整个情况:

test <- function(x) structure(matrix(x, 2, 2), class = "test")

`[.test` <- function(x, i, j, ..., nag = FALSE) {
  if(nag) cat("This won't work now\n")
  NextMethod()
}

x <- test(5)

x[1, 1]
#> [1] 5

x[1, 1, nag = TRUE]
#> This won't work now
#> Error in `[.default`(x, 1, 1, nag = TRUE): incorrect number of dimensions

x[1, 1, nag = FALSE]
#> Error in `[.default`(x, 1, 1, nag = FALSE): incorrect number of dimensions

没有任何直接的解决办法，因为[是一个原始函数.NextMethod的使用意味着命名自变量op实际上不可供终端用户使用.由于该函数没有导出并且没有文档记录，因此这并不是一个真正的错误，但将op作为参数传递而不是在函数体中硬编码`st_intersects似乎是一个奇怪的设计决定.

与此同时，最好避免使用op参数，而直接使用谓词函数:

library(sf)
#> Linking to GEOS 3.9.3, GDAL 3.5.2, PROJ 8.2.1; sf_use_s2() is TRUE

p1 <- st_sfc(st_point(c(0, 0)))
p2 <- st_sfc(st_point(c(1, 1)))

p1[lengths(st_intersects(p1, p2)) != 0]
#> Geometry set for 0 features 
#> Bounding box:  xmin: NA ymin: NA xmax: NA ymax: NA
#> CRS:           NA

^{创建于2024-01-21，包含reprex v2.0.2}