我的目标是在 R 中创建 3D 可视化。我有柏林市区 (Ortsteile) 的 shapefile,并希望将值(居民/平方公里)突出显示为 z 值。我已将 shapefile 实现到 R 中,并为密度值(“Einwohnerd”)着色,如下所示:
library(rgdal)
library(sp)
berlin=readOGR(dsn="C...etc.", layer="Ortsteile")
berlin@data
col <- rainbow(length(levels(berlin@data$Name)))
spplot(berlin, "Einwohnerd", col.regions=col, main="Ortsteil Berlins", sub="Datensatz der Stadt Berlin", lwd=.8, col="black")
如何将某个多边形(市区)引用到 z 值(居民/平方公里)以及如何突出显示该 z 值?
希望有人能解答!
最好的问候
SB
感谢您的回答,但我仍在寻找使用密度作为 z 值的最佳方法,以便我可以创建 3D 模型。我发现不可能使用形状的多边形,但可以光栅化多边形并使用矩阵进行不同的透视和旋转。
这是代码,但最终的 3D 可视化看起来不够清晰和良好。也许以另一种方式计算 z 值会更好,这样第一个值不会开始那么高,或者使用多边形的中心而不是在 z 方向上绘制一列:
library(rgdal)
library(sp)
setwd("C:\\...")
berlin=readOGR(dsn="C:\\...\\Ortsteile", layer="Ortsteile")
col <- rainbow(length(levels(berlin@data$Name)))
spplot(berlin, "Einwohnerd", col.regions=col, main="Ortsteil Berlins",
sub="Datensatz der Stadt Berlin", lwd=.8, col="black")
library(raster)
raster <- raster(nrows=100, ncols=200, extent(berlin))
test <- rasterize(berlin, raster, field="Einwohnerd")
persp(test, theta = 40, phi = 40, col = "gold", border = NA, shade = 0.5)
for(i in seq(0,90,10)){
persp(test, theta = 40, phi = i, col = "gold", border = NA, shade = 0.5)
}
library(rgl)
library(colorRamps)
mat <- matrix(test[], nrow=test@nrows, byrow=TRUE)
image(mat)
persp3d(z = mat, clab = "m")
persp3d(z = mat, col = rainbow(10),border = "black")
persp3d(z = mat, facets = FALSE, curtain = TRUE)