提取包含在空间面中的栅格像元数 [英] Extract number of raster cells included in a SpatialPolygons
本文介绍了提取包含在空间面中的栅格像元数的处理方法,对大家解决问题具有一定的参考价值,需要的朋友们下面随着小编来一起学习吧!
问题描述
我想创建一个函数来计算SpatialPolygonsDataframe对象的多边形内的栅格单元格的数量,并在不使用循环的情况下将值添加为新列。 我找不到该怎么做...
以下是我的代码:
library(sp)
library(raster)
# Create a SpatialPolygonsDataframe and a raster objets to overlay
# Polygons
p1 <- rbind(c(-180,-20), c(-140,55), c(-50, 0), c(-140,-60), c(-180,-20))
p2 <- rbind(c(-10,0), c(140,60), c(160,0), c(140,-55), c(-10,0))
p3 <- rbind(c(-125,0), c(0,60), c(40,5), c(15,-45), c(-125,0))
polys <- as(spPolygons(p1, p2, p3), "SpatialPolygonsDataFrame")
# Raster
p4 <- rbind(c(-180,10), c(0,90), c(40,90), c(145,-10),
c(-25, -15), c(-180,0), c(-180,10))
rpol <- spPolygons(p4)
r <- raster(ncol=90, nrow=45)
grd <- rasterize(rpol, r, fun=sum)
# Function to count the share of occupied cells per polygon
myFunction <- function(polys,grd){
over <- crop(mask(grd, polys), polys)
share <- length(over[!is.na(over)]) / ncell(over) * 100
return(share)
}
myFunction(polys,grd)
非常感谢您的帮助!
推荐答案
可以使用函数raster::extract
。默认情况下,这将提取SpatialPolygonDataFrame的每个多边形的单元信息。输出是一个列表,这样您就可以计算列表中的对象数。
library(sp)
library(raster)
# Create a SpatialPolygonsDataframe and a raster objets to overlay
# Polygons
p1 <- rbind(c(-180,-20), c(-140,55), c(-50, 0), c(-140,-60), c(-180,-20))
p2 <- rbind(c(-10,0), c(140,60), c(160,0), c(140,-55), c(-10,0))
p3 <- rbind(c(-125,0), c(0,60), c(40,5), c(15,-45), c(-125,0))
polys <- as(spPolygons(p1, p2, p3), "SpatialPolygonsDataFrame")
# Create raster
r <- raster(ncol=90, nrow=45)
values(r) <- 1:ncell(r)
# Extract values
r.ext <- extract(r, polys)
# Count number of cell for each polygon
lengths(r.ext)
这篇关于提取包含在空间面中的栅格像元数的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持IT屋!
查看全文