我需要将激光雷达数据(3D点云(转换为反射率和范围地图(图像(,以便将它们用作CNN的通道。
任何已经使用LiDAR数据并可以帮助我的人?
预先感谢您。
您所指的数据产品通常称为"网格指标"。网格指标是在用户定义的网格(例如栅格单元格(中计算出的LIDAR点云高度和强度数据集的描述性统计数据。我同时使用FUSION
(参考P.72(,一组命令行工具和R中的LIDAR处理软件包lidr
(参考P.29(来产生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
文档中的Sytax:
GridMetrics [switches] groundfile heightbreak cellsize outputfile datafile1 [datafile2 ...
datafileN]