Создание отражений и изображений дальности из лидарных данных - PullRequest
1 голос
/ 14 мая 2019

Мне нужно преобразовать лидарные данные (3D Point Cloud) в карты отражения и дальности (изображения), чтобы использовать их в качестве каналов для CNN.

Любой, кто уже работал с лидарными данными и может мне помочь?

Заранее спасибо.

1 Ответ

0 голосов
/ 14 мая 2019

Продукты данных, на которые вы ссылаетесь, обычно называют «метриками сетки». Метрики сетки являются описательной статистикой для наборов данных высоты и интенсивности облака точек LiDAR, рассчитанных в пределах определенной пользователем сетки (например, растровой ячейки). Я использую FUSION (ссылка стр.72) , набор инструментов командной строки, и lidr (ссылка стр.29) , пакет обработки LiDAR в R, для получения метрик сетки LiDAR.

Вот пример из документации, использующей пакет lidr:

LASfile <- system.file("extdata", "Megaplot.laz", package="lidR")
las = readLAS(LASfile)
colors = height.colors(50)
# Canopy surface model with 4 m^2 cells
metrics = grid_metrics(las, max(Z), 2)
plot(metrics, col = colors)
# Mean height with 400-m^2 cells
metrics = grid_metrics(las, mean(Z), 20)
plot(metrics, col = colors)
# Define your own new metrics
myMetrics = function(z, i)
{
metrics = list(
zwimean = sum(z*i)/sum(i), # Mean elevation weighted by intensities
zimean = mean(z*i), # Mean products of z by intensity
zsqmean = sqrt(mean(z^2)) # Quadratic mean
)
return(metrics)
}
metrics = grid_metrics(las, myMetrics(Z, Intensity))
plot(metrics, col = colors)
plot(metrics, "zwimean", col = colors)
plot(metrics, "zimean", col = colors)
plot(metrics, "zsqmean", col = colors)

Вот синтаксис из FUSION документации:

GridMetrics [switches] groundfile heightbreak cellsize outputfile datafile1 [datafile2 ...
datafileN]
...