使用 shapefile 多边形裁剪 .las（激光雷达）文件

Posted 2023-02-25

【中文标题】使用 shapefile 多边形裁剪 .las（激光雷达）文件

【英文标题】：Clipping .las (lidar) files with shapefile polygons

【发布时间】：2014-02-09 17:33:42

【问题描述】：

我有每个包含数百万点的激光雷达文件。我正在尝试使用包含多边形的 shapefile 过滤这些激光雷达文件，基本上是在每个点与多边形之间进行相交或包含。

我是激光雷达的新手，所以对于那些空间裁剪操作的开源库的每一个建议都足够好。

【问题讨论】：

我猜最常见的库是点云库PCL。我不确定这是否适用于 DEM 等地理空间数据，这取决于您的应用程序！

我会阅读有关 PCL 的信息，谢谢您的回复！

Here 我提供了一个详细的例子来解释如何做到这一点。希望你能理解。

【参考方案1】：

PDAL 包含一个过滤器，可以根据几何范围裁剪点。裁剪多边形预计位于WKT。 http://www.pdal.io/apps.html#translate-command 包含一个使用内联 WKT 多边形裁剪文件的示例。一个简单的例子如下所示：

pdal translate --polygon="((0 0, 10 0, 10 10, 0 10, 0 0))" infile.las outfile.las

您还可以在 WKT 中使用文件系统上的文件来代替字符串：

pdal translate --polygon=mywktfile.wkt infile.las outfile.las

PDAL是免费开源的，可以安装在各种系统上；不幸的是，安装文档有点，呃，sparse。我列出了我个人对在各种系统上启动和运行 PDAL 的最简单方法的看法：

Windows：OSGeo4W 在其 64 位安装程序中包含 PDAL 的最新版本 (1.0.0b1-1)（我没有亲自尝试过，因此无法说明它的有效性） OSX：使用homebrew、brew install pdal 其他系统必须build from source（请注意，这些安装说明可能已过时）

【讨论】：

手册页已移至 pdal.io/apps/translate.html，如果需要从 Geopackage 或 SHAPE 文件生成 WKT 字符串，OGR CSV 驱动程序使用 -lco GEOMETRY=AS_WKT 来完成（参见 gdal.org/drivers/vector/csv.html#vector-csv )。