在 R 中的 shapefile 中查找百分比土地覆盖类

【问题标题】：Finding percentage land cover classes in a shapefile in R在 R 中的 shapefile 中查找百分比土地覆盖类
【发布时间】：2020-12-03 02:34:49
【问题描述】：

我有 R 中的 shapefile 和土地覆盖类的 tiff 文件。如何在 shapefile 中找到土地覆盖类别的百分比？例如，该 shapefile 中的研究区域有多少土地？

#example data for cover
rat <- levels(raster2 )[[1]]
rat[["landcover"]] <- c("land","ocean/lake", "rivers")
levels(raster2 ) <- rat

#and a raster
raster1 <-  raster(matrix(c(1,1,1,2,3,4,5,6,7),ncol =3))

【问题讨论】：

标签： r raster tiff

【解决方案1】：

这是一个最小的、独立的、可重现的示例：

栅格数据

library(raster)
r <- raster(ncol=36, nrow=18)
set.seed(111)
values(r) <- sample(3, 18*36, replace=TRUE)

多边形

cds1 <- rbind(c(-180,-20), c(-160,5), c(-60, 0), c(-160,-60), c(-180,-20))
polys <- spPolygons(cds1)

解决方案

x <- extract(r, polys)
table(x)
#x
# 1  2  3 
#16  9 13

【讨论】：

我的价值观是像这些 124886 443141 2570 46059 395207 等，它们是区域的形式吗？因为我想知道每个土地覆盖类别的百分比。