raster:: extract 生成一个空列表答案

【问题标题】：raster:: extract produces an empty listraster:: extract 生成一个空列表
【发布时间】：2020-12-15 17:18:28
【问题描述】：

我正在尝试提取数字化树冠内 DSM(CHM) 中的像素值。

首先我将工作目录设置为在 shapefile 和光栅中读取。

TreeCrowns <-shapefile("plot1sag_shape/plot1sag.shp")
CHM <- raster('272280split4.tif')

然后我尝试提取像素值

pixel <- raster::extract(CHM, TreeCrowns, method= 'simple', weights=FALSE, fun=NULL)

但我得到一个空列表，其中每个多边形的所有 NULL 值。我已经确认 CHM 和多边形位于同一位置。我该怎么做才能解决这个问题？

【问题讨论】：

【解决方案1】：

由于您的 shapefile 由多边形组成，extract() 函数需要知道如何通过 fun= 参数汇总整个多边形的像素值。由于您提供了fun=NULL，因此该函数解释为返回NULL 值以汇总像素值。

试试fun=mean 或fun=sum（它们的含义不同，所以看看哪一个适合你）。

【讨论】：

【解决方案2】：

这可能是因为多边形和栅格不重叠。你能show(CHM) 和TreeCrowns 吗？看过了吗

plot(CHM)
lines(TreeCrowns)

或者您的多边形相对于栅格单元是否非常小？在这种情况下尝试参数small=TRUE

【讨论】：