0

在下面的示例中,我希望在相交的 2D 轮廓下方的 lat-long 平面上绘制地理地图。

不需要使用包“plot3D”或“maps”。

  library(plot3D)

  alt <- rev(seq(0, 500, 5))
  lat <- seq(30, 65, 0.25)
  long <- seq(-135, -65, 0.5)

  grid <- mesh(long, lat, alt)
  names(grid) <- c("long", "lat", "alt")
  colvar <- with(grid, dnorm(long, -100, 20) * 
                       dnorm(lat, 50, 10) *
                       dnorm(alt, 250, 100))

  colvar <- colvar/max(colvar)

  slice3D (long, lat, alt, 
           colvar = colvar,
           ys  = 50, xs = -100, zs = NULL,
           xlab = "long",
           ylab = "lat",
           zlab = "alt",
           theta = 60, phi = 40,
           ticktype = "detailed")

示例地图:

library(maps)
plot(NULL, xlim = range(long), ylim = range(lat), 
     xlab = "long", ylab = "lat")
map("world", add = TRUE)

在此处输入图像描述

在此处输入图像描述

4

0 回答 0