r - raster:: extract 生成一个空列表

Question

我正在尝试提取数字化树冠内 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 和多边形位于同一位置。我能做些什么来解决这个问题？

标签： rrasterlidar

Answer 1

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

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