如何获取LIDAR块的唯一“点源ID”？

Question

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

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

这就是我到目前为止所尝试的：

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)

不幸的是，整个过程相当慢，对于更大的数据集，我得到了堆栈溢出错误，最终永远不会到达“写文件”部分。

Answer 1

这个过程比预期的要慢，因为您是在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)