1

我正在尝试使用 laspy 获取有关出现在 LIDAR 数据块中的独特航线的信息。

我已经尝试为整个块运行 lasInfo 模块，但我得到的只是一个最小和最大 point_source_ID 值，而不是我需要的单个飞行路线列表。

这是我迄今为止尝试过的：

import laspy
import glob

las_files_list = glob.glob(r'PATH\*.las')
print(las_files_list)

las_source_id_set = set()

for f in las_files_list:
    las_file = laspy.file.File(f, mode='r')
    las_source_id_list = las_file.pt_src_id


    for i in las_source_id_list:
        las_source_id_set.add(i)

    las_file.close()
    print(las_source_id_set,' ', f)

print(las_source_id_set)

with open('point_source_id.txt', 'w') as f:
    f.write(las_source_id_set)

不幸的是，整个过程相当缓慢，并且使用更大的数据集，我得到一个堆栈溢出错误，最终永远无法进入“写入文件”部分。

4

1 回答 1

0

该过程比它可能的要慢，因为您正在 Python 中执行循环。有一个 numpy 函数可用于加快处理速度：numpy.unique

您的脚本将变为：

import laspy
import glob
import numpy as np

las_files_list = glob.glob(r'PATH\*.las')
print(las_files_list)

las_source_id_set = set()

for f in las_files_list:
    with laspy.file.File(p) as las:
        las_source_id_set.update(np.unique(las.pt_src_id))


    print(las_source_id_set,' ', f)

print(las_source_id_set)

with open('point_source_id.txt', 'w') as f:
    f.write(las_source_id_set)
于 2019-04-26T10:53:38.483 回答