Question

我具有提取点序列时间序列栅格值的功能。到现在为止，这一切都很好。但是，今天尝试使用它时，我经常收到错误消息。我收到的错误提示：

UseMethod（“ extract_”）中的错误：没有适用的方法 “ extract_”应用于类“ c（'RasterLayer'，'Raster'， 'BasicRaster'）“

请您帮忙。我已经尝试了很多，但无法解决此问题。

library(raster)
#--------------------------------------start myfunction-----------------------------

# extract raster values for each single point in each path and put it a dataframe
lst <- function (rstack.lst, points, df.name, IDstr){
  df.name <- list()
  ii <- 1
  #number of layers in a list of raster stacks
  n = Reduce(`+`, lapply(rstack.lst, nlayers))

  for (j in 1:length(rstack.lst)){
    df.name[[j]] <- as.data.frame(matrix(0, ncol = nlayers(rstack.lst[[j]])+3, nrow = nrow(points)))
    names(df.name[[j]]) <- append(c("coords.x","coords.y","ID"), substr(names(rstack.lst[[j]]),5,14), after = 3)
    #calculating x and y coordinates
    df.name[[j]][1:2] <- coordinates(points) 
    # setting up unique IDs 
    df.name[[j]]$ID <- paste0(IDstr, ".",seq(1:(nrow(df.name[[j]])))) 

    for (i in 1:nlayers(rstack.lst[[j]])){
      df.name[[j]][i+3] <- extract(rstack.lst[[j]][[i]], points)
      colnames(df.name[[j]][i+3]) <- substr(names(rstack.lst[[j]][[i]]),5,14)
      cat(paste0(round((ii/(n))*100), '% completed'))
      ii <- ii+1
      Sys.sleep(.05)
      if (ii == n) cat(': Done')
      else cat('\014')
    }

  }
  return(df.name)
}
#--------------------------------------end start myfunction-----------------------------

lsdata[[1]] <- lst(rstack.lst=r.lst, points=Pnt.shp, IDstr="P1")

Answer 1

我可能是错的，但是我认为extract函数没有故障，但是您可能正在加载其他一些具有完全相同功能的软件包（我的意思是相同的名称）。许多软件包可能具有相同名称的功能。例如，extract也是tidyr中的函数。为避免此类问题，建议您在指定函数时也添加程序包名称。您可以在这里执行：raster::extract。

要确定，只需执行extract并仔细检查它包含的内容。无论包含什么内容，可能都不会如下所示：

function (x, y, ...) 
standardGeneric("extract")
<bytecode: 0x000000001873c4d0>
<environment: 0x000000001855db68>
Methods may be defined for arguments: x, y
Use  showMethods("extract")  for currently available ones.

执行使用空间点提取时间序列值的函数时出现错误消息

1 个答案: