单个封闭轮廓内的R区域

r geospatial spatial contour polygons

95 观看

1回复

37 作者的声誉

使用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空间包。

作者: Peter 的来源 发布者: 2017 年 9 月 15 日

回应 1


0

37 作者的声誉

我能够通过转换为栅格并使用栅格包中的簇来解决此问题:

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)

这个例子直接来自栅格数据包-我已经发布了它,因为它可能与其他人有关!

作者: Peter 发布者: 2017 年 9 月 17 日
32x32