函数在全局环境中作为对象工作，但不是来自包

Question

我有一个函数可以创建栅格图层的样本变异函数。该函数采用RasterStack 或RasterBrick 并返回一个列表，其中每层有一个变异函数。该函数在 .GlobalEnv 中加载时工作正常，但是当我从包中使用它时，我收到一条错误消息：

Error in (function (classes, fdef, mtable) : unable to find an inherited method for function ‘nlayers’ for signature ‘"numeric"’

回溯说出现问题我使用rasterToPoints 函数的行，以及我在dlist 上使用lapply 的行，但我不知道为什么要使用nlayers()，或者为什么在尝试时它只是一个问题使用包中的函数。

这里是 roxygen2 cmets 的函数：

#' Create sample variograms of CHIRPS data
#'
#' @name sampvar-chrp
#' @title Create sample variograms of CHIRPS data
#'
#' @param data RasterStack of CHIRPS images
#' @return List of sample variograms
#' @examples
#' data(chrprf)
#' sampvar_chrp(chrprf)
#' @export
sampvar_chrp <- function(data) {
  dlist <- as.list(data)
  spdf <- lapply(dlist, function(x){
    names(x) <- "rain"
    pts <- rasterToPoints(x, spatial = TRUE)
  })
  rsamp <- lapply(spdf, function(x){
    s <- x[sample(1:length(x), size = 150), ]
  })
  chrp_vario <- lapply(rsamp, function(x){
    v <- variogram(object = rain ~ 1, data = x)
  })
  return(chrp_vario)
}

感谢任何建议！

编辑：我认为（至少部分）问题是data 在运行包中的函数时没有成功转换为列表。我编辑了函数，重新记录，重新安装，得到了"error"打印语句。所以看起来as.list(data) 在全球环境中工作，但在我的包中不起作用。这是编辑后的函数：

#' Create sample variograms of CHIRPS data
#'
#' @name sampvar_chrp
#' @title Create sample variograms of CHIRPS data
#'
#' @param data RasterStack of CHIRPS images
#' @return List of sample variograms
#' @importFrom raster rasterToPoints
#' @examples
#' data(chrprf)
#' sampvar_chrp(chrprf)
#' @export
sampvar_chrp <- function(data) {
  dlist <- as.list(data)
  if (is.list(dlist) == TRUE) {
    spdf <- lapply(dlist, function(x){
      names(x) <- "rain"
      pts <- rasterToPoints(x, spatial = TRUE)
    })
    rsamp <- lapply(spdf, function(x){
      s <- x[sample(1:length(x), size = 150), ]
    })
    chrp_vario <- lapply(rsamp, function(x){
      v <- variogram(object = rain ~ 1, data = x)
    })
    return(chrp_vario)
  } else { print("error")}

}

Answer 1

如果您查看raster::rasterToPoints，您可以看到

head(rasterToPoints)

#1 function (x, fun = NULL, spatial = FALSE, ...)                                      
#2 {                                                                                   
#3     nl <- nlayers(x)           

它首先调用的是nlayers。错误消息表明它得到一个numeric 向量，而它需要一个Raster* 对象。

问题似乎是您的变量data 不是Raster* 对象。所以你需要（至少）退后一步看看发生了什么。