在Plotly(R API)

时间:2016-08-25 09:53:50

标签: r geolocation raster plotly

我想在R中使用Plotly来创建考古发掘的三维模式。我非常成功地绘制了点和曲面数据(例如:Vignette我正在研究的R package),但我还想添加沟槽地理参考剖面照片的栅格信息

我没有找到任何方法在Plotlys 3D环境中绘制栅格数据。到目前为止我提出的唯一解决方案(感谢this帖子)是使用SFM使用Photoscan创建配置文件的3D模型,将彩色网格导出为.ply文件,修复此文件的标题并将其导入R以使用以下示例代码进行绘图:

library(geomorph)
library(plotly)

#load data
mesh <- read.ply("plotly/expply8_corr.ply", ShowSpecimen = FALSE)

# extract vertex coordinates
x <- mesh$vb["xpts",]
y <- mesh$vb["ypts",]
z <- mesh$vb["zpts",]

# plot
plot_ly(
  x = x, y = y, z = z,
  i = mesh$it[1,]-1, j = mesh$it[2,]-1, k = mesh$it[3,]-1,
  facecolor = c(mesh$material$color[1, ]),
  type = "mesh3d"
)

您将找到示例数据here

不幸的是,这种情况非常糟糕。如果增加网格分辨率,一切都会变慢。我真的想添加一个简单的地理参考栅格来保持高性能,并避免创建配置文件的3D模型的必要性。是否有使用Plotly或其他绘图库实现此目的的工作流程?

1 个答案:

答案 0 :(得分:1)

我找到了包rgl的一个很好的解决方案。例如:

library(rgl)
library(jpeg)

# download and load picture
download.file(
  url = 'https://upload.wikimedia.org/wikipedia/en/6/6d/Chewbacca-2-.jpg',
  destfile = "chewbacca.jpg",
  mode = 'wb'
)

chewie <- readJPEG("chewbacca.jpg", native = TRUE)

# create some sample data
x <- sort(rnorm(1000))
y <- rnorm(1000)
z <- rnorm(1000) + atan2(x, y)

# plot sample data
plot3d(x, y, z, col = rainbow(1000), size = 5)

# add picture
show2d(
  # plot raster
  {
    par(mar = rep(0, 4))
    plot(
      0:1, 0:1, type="n",
      ann = FALSE, axes = FALSE,
      xaxs = "i", yaxs = "i"
    )
    rasterImage(chewie, 0, 0, 1, 1)
  },
  # image position and extent
  # coordinate order: lower left, lower right, upper right and upper left
  x = c(-2,  1,  1, -2),
  y = c(-1, -1,  1,  1),
  z = c(-3, -3,  2,  2)
)

图片必须与其他软件进行地理配准(GIS / CAD中的摄影测量)。如果你有地理参考栅格,你只需要角点的坐标来绘制它。