Создание ортофотоплана из облака точек *.las в R

Есть ли способ экспортировать облако точек *.las в R в ортофотоплан? Я загрузил свой las-файл, содержащий точки, с пакетом lidR. Я хочу экспортировать tif, который показывает облако точек сверху в RGB, аналогично тому, как будет выглядеть ортофото. Данные были собраны с помощью наземного лазерного сканера.

облако точек

это то, что я хочу


r lidar tiff lidar-data las
person Zoe    schedule 28.03.2021    source источник
comment
Эта глава книги может помочь вам jean-romain.github. io/lidRbook/. Не существует единого способа добиться этого, поэтому вы должны быть более конкретными в том, как вы хотите агрегировать свое облако точек.   -  person JRR    schedule 28.03.2021

Ответы (1)


arrow_upward
1
arrow_downward

Итак, я понял, как это сделать, хотя это не очень элегантно:

# 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)
person Zoe    schedule 29.03.2021
comment
Вы можете вычислить which.max(z) только один раз, и ваше решение будет элегантным ;-) - person JRR; 29.03.2021