时间戳

首页 > 解决方案 > 在 Python 中将点云拆分为特定大小的单元格

python - 在 Python 中将点云拆分为特定大小的单元格

问题描述

我需要将点云拆分为特定的单元格大小，以便在该单元格中找到 Z min 和 Z max。我将点云转换为数组并创建网格网格，但我不知道如何将点云/数组拆分为由网格网格定义的扩展。建议？还是有更合适的方法？

import laspy
import numpy as np
inf = laspy.file.File('SAVA_000012.las', mode='r')
points = inf.points
points_xyz = np.array((points['point']['X'],
                   points['point']['Y'],
                   points['point']['Z'])).transpose()
x = np.arange(min(inf.X), max(inf.X), 200)
y = np.arange(min(inf.Y), max(inf.Y), 200)
print(x)
xx, yy = np.meshgrid(x, y, sparse=True, indexing= 'xy')

提前致谢

标签： pythonsplitpoint-clouds

解决方案

使用https://github.com/daavoo/pyntcloud：

from pyntcloud import PyntCloud

cloud = PyntCloud.from_file("SAVA_000012.las")

# Using 0.15 just for example. Omit `z` for 2D grid
voxelgrid_id = cloud.add_structure("voxelgrid", size_x=0.15, size_y=0.15)

voxelgrid = cloud.structures[voxelgrid_id]

z_max = voxelgrid.get_feature_vector(mode="z_max")

您可以在此处了解有关 VoxelGrid 的更多信息：

https://github.com/daavoo/pyntcloud/blob/master/examples/%5Bstructures%5D%20VoxelGrid.ipynb

z_min目前不支持，但它应该很容易实现（贡献很好^^）以下作为参考：

https://github.com/daavoo/pyntcloud/blob/master/pyntcloud/utils/numba.py#L19

用于：

https://github.com/daavoo/pyntcloud/blob/master/pyntcloud/structures/voxelgrid.py#L161

推荐阅读