问题标签 [lidar]

我正在使用 liblas for python 来读取 .las 文件。当我输入：

它给了我：

没有名为 liblas 的模块。

我已经在系统中设置了 las 库路径，lasinfo工作正常。谁能告诉我如何在 Python 中导入 las 库？顺便说一句，我正在使用 Ubuntu。

I am a new MATLAB user with little programming experience (I have a mechanical engineering background) so I apologise in advance if this is a simple question!

I am trying to import a large point cloud file (.pts file extension) into MATLAB for processing. I'm lead to believe that the file contains a text header and 3 columns of integer data (x, y and z coordinates) - I managed to open the first part of the file as a text file and this is the case.

I cannot import the file directly into MATLAB as it is too large (875 million points) and can only import it 9000000 rows at a time, therefore I have written the script below to import the file (and consequently save) as 9000000x3 blocks, saved as MATLAB files (or another appropriate format).

Script:

The script runs fine and cycles through 5 iterations, importing 5 blocks of the data. Then, as it attempts to import the 6th chunk I get the following error:

I am unsure about what is causing the error, I believe it is relating to fread command size, as I have experimented with various values such as 3, which enables just one block to be imported before the dimension mismatch error occurs.

Once more I apologise if I am missing something very basic, my understanding of programming techniques is very limited only having been introduced to it a couple of months ago.

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

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

我正在尝试使用libLAS Python API将点数据集写入*.las文件。但是我遇到了一些浮点值被四舍五入的问题

如果我设置，pt.raw_x而不是pt.x我看不到四舍五入的问题，但没有las写入文件。

我不确定我错过了什么。我将不胜感激。

我试图在分类（植被）.las 数据集的点云中提取点的裸坐标。

使用 Arcgis 并查看属性表不是我想要的，所以我问自己，如果我可以将 .las“反转换”为 ASCII 左右，以获得这些坐标。我希望有人能理解我的问题。

编辑：我设法通过一个名为 Arcmap 10.2 的简单工具集获得了我想要的东西Featureclass Z to ASCII (3d Analyst)

我必须生成 3D 点云的网格。所以我用这个delaunay函数来执行点的三角剖分。激光雷达数据是人头扫描的结果。

当我使用delaunay两个输入时，它给了我输出但并不完美。所以我使用了三个输入 ( X,Y,Z)

但现在结果变得更糟了。我不知道是什么问题？

这是我编写的完整代码：

Head_cloud_point是具有 X,Y,Z 坐标的文件。我必须使用这些坐标生成网格。

我有一个使用 FUSION（激光雷达处理）命令执行以下工作流程的脚本：

1. 将 LiDAR LAS 文件剪辑到多边形边界，为每个多边形创建一个唯一的 LAS 文件
2. 运行cloudmetrics以基于每个剪辑的 LAS 文件记录点云信息

这是我为完成第 1 步和第 2 步而创建的脚本：

.bat如果我通过双击文件将此脚本作为文件运行，它会在polyclipdata命令中创建某种无限循环。但是，如果我在命令行上逐行执行脚本，工作流将执行所需的任务。这里哪里出错了？

我lasinfo使用 .las 文件运行并收到以下错误消息。

错误：主版本超出范围

我想知道是什么原因，我该如何解决？

我想用 3D Delaunay Triangulation 来表示一个表面。顶点必须是我的原始输入数据，来自市区的 LiDAR 点云。因此，表面必须调整/适应输入信息。

实际上，我需要做的是：

• 我有一个来自市区的 3D 点云 (x, y, z)；
• 我需要代表这个区域的表面；
• 我想做一个 3D delaunay 三角剖分（我用 CGAL 做过，我得到了四面体）并只识别代表表面的三角形（用 CGAL 我有 4 个顶点，我无法定义哪些代表表面三角形）；
• 由于这些表面三角形是已知的，因此我需要给出一个点并获取包含该给定点的三角形。

我想知道我需要哪个功能。我看到了“3D 表面网格生成”、“3D Delaunay 三角剖分”、“点集的表面重建”、“3D 网格生成”和“3D Alpha 形状”。但我不确定哪一个更适合我的工作/需要。什么函数可以给我表面三角形（由我的初始数据集中的点表示）作为输出。

有人有建议吗？

我用 CGAL 执行了 3D alpha 形状函数，得到了意想不到的结果。我的输入数据是一组 3D 点 (x, y, z)，它们代表平坦区域中的一个建筑物（盒子）（坐标中有一些噪音 - 小的噪音）。我想我只会得到代表建筑物（墙壁和屋顶）和地面的表面三角形。

但是，结果我得到了三角形，形成了表面的凸包。我试图改变“最佳阿尔法值”，但它是一样的。

是否有任何过滤过程或参数可以设置为仅获取表面三角形？