Question

我的目标是从x，y和z坐标的数据框创建一个三维多面体。我希望值类似于SpatialPolygons对象，但欢迎其他建议。然后，我想将gBuffer函数扩展为3D多面体。 gBuffer坚持从原始几何体到新几何体的距离，是我最终目标的理想特征。

下面的代码说明了2D中的情况：

library(sp)
library(rgeos)
#Create a SpatialPolygons object from set of x-y coordinates
xy2SP <- function(xy, ID=NULL) {
if(is.null(ID)) ID <- sample(1e12, size=1)
SpatialPolygons(list(Polygons(list(Polygon(xy)), ID=ID)),
                proj4string=CRS("+proj=merc"))
}
xy <- data.frame(x=c(0,2,3,1,0), y=c(0,0,2,2,0))
plot(xy)

sp1 <- xy2SP(xy)
sp2 <- gBuffer(sp1, width=.2)
sp3 <- gBuffer(sp1, width=.4)
sp4 <- gBuffer(sp1, width=.6)

plot(sp1,xlim=c(-1,4),ylim=c(-1,3))
plot(sp2,add=T,border=2)
plot(sp3,add=T,border=2)
plot(sp4,add=T,border=2)

我的第一次尝试是将我的函数xy2SP用于3坐标而不是2坐标，但很快就会被挫败。在任何情况下，都有一个数据框的例子，它表示3D多面体的顶点。

library(scatterplot3d)
xyz <- data.frame(x=c(0,2,3,1,0,2,3,1),y=c(0,0,2,2,0,0,2,2),z=c(0,0,0,0,2,2,2,2))
sp<-scatterplot3d(xyz,pch=16)   
sp1 <- xy2SP(xyz)

从顶点坐标和偏移多面体创建多面体

0 个答案: