问题标签 [pdal]

我正在使用带有 LAS 和 LAZ 格式的点云文件的PDAL命令行界面。我想知道 LAZ/LAS 文件中是否有索引机制。有没有办法使用 PDAL 来获得它。

我需要阅读 .las 文件并通过 google colab 运行代码。

我已经尝试了几个库，如 lilbas 和 pdal。但是，在尝试导入它们时，我遇到了一些难以解决的错误。

为了安装库，我运行了以下代码： !pip install liblas 或 !pip install pydal

但是我收到以下错误：对于 liblas：OSError：liblas_c.so.3：无法打开共享对象文件：没有这样的文件或目录

对于 pdal：使用缓存的https://files.pythonhosted.org/packages/e4/67/338c4164a02b1bb0ecb34adca254b29e88c6c3ed874402d2cd67f5c611a5/PDAL-2.1.8.tar.gz 错误：命令“python setup.py egg_info”在/tmp 中失败，错误代码为 1 /pip-install-tiflm77w/pdal/

你可以帮帮我吗？

我在 Conda Python3.5.6 中安装了 PDAL。在手动链接许多“.so”文件后，当我想将 LAS 转换为 LAZ 文件时仍然出错。

有人可以告诉我如何解决吗？

所以这是我第一次使用 PDAL。我使用 python 3.6 和 PDAL 1.9。

它显示错误，

我检查了网站上的示例代码，它看起来一样。只是不知道为什么它不起作用？

我正在使用几个 python 库（whitebox、PCL、PDAL）对 3D 点云进行分类。我的目标是对土壤进行分类。该数据集已被一家公司分类，因此我将其分类为基本事实。

目前我能够对土壤进行分类，为此我对数据集进行了解密并使用 PDAL 重新进行分类。现在我正处于比较两个数据集以查看我的分类质量的阶段。

我制作了一个脚本，它获取 2 个集合的 XYZ 坐标并将其放入一个列表中，然后我将它们逐个进行比较，但是数据集包含大约 500 万个点，并且在开始时需要 1 分钟乘以 5 个点。几分钟后，一切都崩溃了。谁能给我提示？这里是我的云的照片 Lets 的集合是基本事实，右边是我分类的那个

我想将 .las 点云文件中的所有数据写入 PostgreSQL/postgis 数据库。

我的管道文件如下：

这很好用，我已添加"output_dims":"X,Y,Z,RED,GREEN,BLUE"到作者中。 ST_AsEWKT(PC_EXPLODE(pa)::geometry)返回SRID=32635;POINT(580985.1431 4275315.1888 230.0292)，没有 RGB 数据。

返回：

颜色值看起来像是 longint，我想将所有这些值写入数据库。

我在我的 Ubuntu 18 中以这种方式构建 PDAL：

运行PGUSER=postgres PGPASSWORD=*** PGHOST=localhost PGPORT=5432 ctest可以确认一切正常。

但是当我尝试检查 LAZ 文件时，我收到了这个错误：

PDAL: readers.las: Can't read compressed file without LASzip or LAZperf decompression library.

这是我的管道文件：

我认为 lazperf 还可以，因为 pgpointcloud 不会抱怨PGUSER=postgres PGPASSWORD=**** PGHOST=localhost make installcheck并告诉我：

PDAL 测试没有告诉我任何关于压缩的信息。

我如何构建或告诉 PDAL 我的 LAZPerf 安装？

编辑 pdal info install/PDAL/test/data/las/autzen_trim.las一切正常。

我正在开发一个基于激光雷达数据从头开始编写 LAS 文件的应用程序。以前，该应用程序是使用 LibLAS 构建的，并且运行良好。但是，LibLAS 已停产，他们建议将其替换为 PDAL，后者是用于处理点云数据的更通用的库。

到目前为止，除了 LAS 标头之外，几乎所有内容都可以轻松地从 LibLAS 转换为 PDAL。以前，在 LibLAS 中，我们将定义一个 liblas::header，我们将在其中添加诸如 LAS 版本等数据。然后我们将在写入文件时使用该标头，它将部分定义文件的写入方式。像这样：

现在，使用 PDAL，我似乎找不到任何东西可以将标题添加到我正在编写的文件中。经过无数次谷歌搜索后，我似乎找不到任何可以帮助我解决这个问题的东西。

我找到了这个教程，但它没有提到标题：https ://pdal.io/development/writing.html

我找到了本教程，但它使用了正在读取的文件的标题，这不是我想要的：https ://pdal.io/api/transition/index.html

我正在使用一个从激光雷达数据生成数字高程模型的程序。代码是使用 PDAL 2.3.0 库在 python 3.8 中编写的。我需要将重投影从 EPSG 3395 更新到 EPSG 4326。当我使用 EPSG 4326 运行管道时，会生成一个全黑的 geotiff，但没有数据。不知道我在这里做错了什么。我已经附加了传递给 pdal 的 JSON 数据。

我有许多分类的激光雷达点云文件，我想将它们转换为 geotiff 光栅文件。为此，我编写了一个函数，该函数创建一个使用 PDAL 进行转换所需的 json-Pipeline 文件，然后执行该管道。

哪里classification是包含与地面/建筑物/植被相对应的数字的字典，所以我不必记住这些数字。

通过遍历tiles. 但是，由于我有很多文件，我想为此使用多个内核。如何拆分任务以至少使用我机器中的所有四个内核？我试图用以下方法做到这一点：

但这不起作用，因为我得到了AttributeError: 'list' object has no attribute 'split'. 但我没有通过列表，或者我是吗？这通常是并行化的方法吗？