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

【问题标题】：Clipping .las (lidar) files with shapefile polygons使用 shapefile 多边形裁剪 .las（激光雷达）文件
【发布时间】：2014-02-09 17:33:42
【问题描述】：

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

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

【问题讨论】：

【解决方案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 的最简单方法的看法：

【讨论】：

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