从R中的*.las点云创建正交马赛克

Question

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

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)