到 2026年,一个国家项目旨在绘制整个法国领土的三维地图。 美国国家地理和林业信息研究所 (IGN) 是这个名为“Lidar HD”的“项目”的幕后黑手,其目标是获得该地区非常精确的 3D 描述。
NSDT工具推荐: Three.js AI纹理开发包 - YOLO合成数据生成器 - GLTF/GLB在线编辑 - 3D模型格式在线转换 - 可编程3D场景编辑器 - REVIT导出3D模型插件 - 3D模型语义搜索引擎 - Three.js虚拟轴心开发包 - 3D模型在线减面 - STL模型在线切割
3D地图将满足多领域观察和分析的需求:
- 预防自然风险。
- 林业资源知识。
- 共同农业政策的监测。
- 能源转型和可持续城市发展。
- 区域规划。
- 内部安全。
- 揭示考古遗迹。
该项目于2021年启动,预计于2026年完成,需要飞行7000小时才能覆盖全境。 目前,法国南部地区可供下载:
在本文中,我们将从尼姆(竞技场周围)下载 LIDAR 数据,在 PDAL 中处理数据并在 Cesium 3D Tiles 中可视化。 你已经可以在这里观察结果。
1、数据下载
在这里可以从地图获取激光雷达数据:
对于尼姆,有文件 LIDARHD_1-0_LAZ_NP-0808_6305-2021.7z (405MB)。
下载并解压后有4个LAZ文件:
- Semis_2021_0808_6304_LA93_IGN69.laz (93MB)
- Semis_2021_0808_6305_LA93_IGN69.laz (107MB)
- Semis_2021_0809_6304_LA93_IGN69.laz (98MB)
- Semis_20 21_0809_6305_LA93_IGN69.laz (96MB)。
使用PDAL我们可以得到一些统计数据:
$ pdal info Semis_2021_0808_6304_LA93_IGN69.laz --all
显示了很多统计数据,例如面积(1121100)、计数(19148173)、创建年份(2021)。
2、数据处理
我们将使用 PDAL 合并 4 个文件 ( filters.merge
),选择竞技场周围的区域 ( filters.crop
),根据 z 值对激光雷达点进行分类并写入结果。
处理管道如下所示:
有关处理管线的代码,请参阅这里。
我们可以使用以下方式运行管道:
$ pdal pipeline pipeline.json
结果是文件 result.las (32MB)。 我们在 CloudCompare 中加载此文件:
下一步是创建 3D Tiles 以在网络上可视化激光雷达数据。 我们将使用工具Go Cesium Point Cloud Tiler :
$ gocesiumtiler -srid 2154 -input result.las -output 3dtiles -zoffset 60
作为“srid”,我们使用 2154 法国坐标系。 作为“zoffset”,我们使用 60 米,以便能够同时显示地形和激光雷达数据。
将创建一个包含大量文件的“3dtiles”文件夹。 LIDAR 数据位于 pnts 文件中,规格请参见这里。切片的空间放置在tileset.json 文件中定义,使用八叉树切片方案。
3、在 Cesium 中可视化
最后一步我们可以在 Cesium 中添加 3D Tileset,代码参见这里。这里我们根据 z 值定义分类点的各种颜色。
可以从这里查看最终的结果:
4、结束语
在本博客中,我们能够从法国下载原始激光测距数据并创建 3D 点云图块。
看到下载的 zip 文件包含多个 LAZ 文件有点令人惊讶,但我们可以使用 PDAL 合并这些文件。
另一个问题是未使用“分类”字段(至少在这些文件中),也许可以添加道路、水、建筑等分类。
作为下一步,我们可以更多地分析 LIDAR 数据,例如通过检测/提取对象......但那是另一个博客......
到目前为止,法国的第一个LiDAR下载看起来很有希望。
原文链接:LiDAR点云转3DTiles - BimAnt