光栅：在R中也出现计数零点答案

【问题标题】：raster: Count zeros points occurrence too in R光栅：在R中也出现计数零点
【发布时间】：2020-10-24 13:33:19
【问题描述】：

当我没有点时，我无法按像素计算光栅中的点数。在我的例子中：

#Create a raster and some points
xmn = 17;xmx = 23;ymn=42;ymx=49
my_ras<-raster(matrix(1:12,3,4), xmx=xmx, xmn=xmn,ymx=ymx,ymn=ymn)
my_points<-data.frame(x=runif(10,xmn,xmx), y=runif(10,ymn,ymx))

#Count the points
pointcount<- rasterize(my_points, my_ras, fun='count')

#Vizualize
plot(my_ras)
points(my_points, pch=16)
d <- data.frame(coordinates(pointcount), count=pointcount[])
summary(d)
#       x               y             count      
# Min.   :17.75   Min.   :43.17   Min.   :1.000  
# 1st Qu.:18.88   1st Qu.:43.17   1st Qu.:1.000  
# Median :20.00   Median :45.50   Median :1.500  
# Mean   :20.00   Mean   :45.50   Mean   :1.667  
# 3rd Qu.:21.12   3rd Qu.:47.83   3rd Qu.:2.000  
# Max.   :22.25   Max.   :47.83   Max.   :3.000  
#                                 NA's   :6

在这里，问题是我没有没有点的像素的 0 值。我无法制作零栅格，因为在我的栅格中的现实世界中，我在目标区域之外有 NA 值，这会导致实际零（无点）和区域边界外的零之间出现错误。

如果有（显然在这种情况下没有），我需要零和 NA 值，有什么想法吗？也许sf 包有这样的东西。

【问题讨论】：

标签： r raster sf

【解决方案1】：

示例数据

library(raster)
xmn = 17;xmx = 23;ymn=42;ymx=49
ras <-raster(matrix(1:12,3,4), xmx=xmx, xmn=xmn,ymx=ymx,ymn=ymn)
points <-data.frame(x=runif(10,xmn,xmx), y=runif(10,ymn,ymx))

您可以使用background=0 代替 NA

pointcount<- rasterize(points, ras, fun='count', background=0)

如果您现在想排除某些区域，可以使用raster::mask。或者，如果您有一个 NA 和零值的栅格，您可以使用 rasterize 和 update=TRUE

【讨论】：