根据R中重叠的栅格单元的数量来划分多边形值

Question

目前我正在使用raster包。我有不同的多边形具有某些值（假设为100），我设法rasterize。问题是，当我rasterize每个栅格单元格结果的值为100，但我想要多边形值（100）时，每个单元格会覆盖多边形。例如，如果多边形覆盖100个栅格单元格，我希望每个栅格单元格的值为1，而不是100.任何人都可以帮助我吗？ 这里的光栅是我创建的shp：https://drive.google.com/drive/folders/0B6-UFgI67v99c3ZhUFp0eWpzOGM 我试着这样做：

ncell<-freq(union,digits=6)
ncell[,"value"]/ncell[,"count"]
new<-rep(c(union[,"value"],ncell[,"count"]))
union$new<-c(new)

但我无法加入我在栅格中将栅格单元格与新值相关联的列。

Answer 1

我能想到两种方式：

1. 在使用栅格化之前计算多边形所需的数字

2. 像你一样使用freq，但是然后使用subs

例如：

r <- subs(union, data.frame(ncell)) 
x <- union / r