单个闭合轮廓内的 R 区域答案

【问题标题】：R area within individual closed contours单个闭合轮廓内的 R 区域
【发布时间】：2018-02-24 05:47:55
【问题描述】：

使用 R - 我希望计算封闭在特定闭合轮廓级别内的区域（以形成多边形）。这已在其他地方得到解决，但如果这些闭合轮廓彼此隔离，我也希望单独计算面积。

一个例子：

library(splancs)
x <- c(1,0,0,1,1,1.5,1.5,3,3,1)
y <- c(0,0,1,1,0,0,2,2,0,0)
m <- cbind(x, y)

plot(m, type="b")
areapl(m)

这给出了 area=4 .... 但是我希望将输出作为与每个孤立多边形的单独区域相对应的区域 c(1,3) 的列表。

这只是一个玩具示例。在我的真实数据中，我首先不知道多边形的坐标，也不知道会有多少孤立的多边形。

任何想法都会很棒 - 我愿意使用其他 R 空间包。

【问题讨论】：

如果数据的格式可以很容易地从sp 包导入到 Spatial* 中，并且如果每个孤立的多边形都有自己的 id，rgeos::gArea(m,byID=T) 将满足您的需求。

【解决方案1】：

我能够通过转换为光栅并使用光栅包中的丛来解决这个问题：

r <- raster(ncols=12, nrows=12)
set.seed(0)
r[] <- round(runif(ncell(r))*0.7 )
rc <- clump(r) 
freq(rc)
plot(r)
plot(rc)

此示例直接来自光栅包 - 我已发布它，因为它可能与其他人相关！

【讨论】：