从 R 中的 *.las 点云创建正射镶嵌

Question

有没有办法将 R 中的 *.las 点云导出为正射镶嵌？我用包 lidR 加载了包含点的 las 文件。我想导出一个 tif，它以 RGB 从上方显示点云，类似于正射影像的样子。数据是使用地面激光扫描仪收集的。

point cloud

this is what I want

Answer 1

好的，所以我想出了怎么做，虽然它不是很优雅：

# load data
points <- readLAS(input_path)

# returns the RGB values for the highest points
RGBZ <- function(r,g,b,z) {
  bands = list(
    R = r[which.max(z)],
    G = g[which.max(z)],
    B = b[which.max(z)]
  )
  return(bands)
}

# create & save ortho
ortho <- grid_metrics(points, ~RGBZ(R,G,B,Z), res = 0.1)
writeRaster(ortho, output_path)