我使用的是多脉冲激光雷达数据,这些数据沿飞行路径内的多条线收集点。我正在尝试确定las文件中各个行的名称和数量。我在python中使用liblas模块。在
我找到了this documentation,它解释了las文件中存储的不同字段。它在页面底部提到了一个数据字段(get_data和set_data)。在
标头中的“点数据格式”和“点数据记录长度”为该“数据”字段留出了空间。我的头文件说我为数据字段预留了28个字节,数据字段中存储了28个值。第19个值(至少在来自两个不同传感器的两个数据集中)指的是行号。我在单脉冲数据中有一个值,在多脉冲数据中有4个值。在
我想知道是否有一个标准来存储这个字段中的内容,或者它是否是专有的。在
另外,为了得到每个扫描行的名称,我写了以下代码:
import liblas
from liblas import file as lasfile
# Get parameters
las_file = r"E:\Testing\00101.las"
f = lasfile.File(las_file, mode='r')
line_list = []
counter = 0
for p in f:
line_num = p.data[18]
if line_num not in line_list:
line_list.append(line_num)
counter += 1
print line_list
它会导致以下错误:
^{pr2}$有没有人知道更多关于las点/标题中存储的行号的信息?有人能解释一下这个错误吗?在我得到错误之前,它似乎分配了将近2gb的ram。我在WinXP上,所以我猜这是内存错误,但我不明白为什么访问这个“数据”字段占用内存。非常感谢任何帮助。在
我并不假装是这方面的专家,但我对地理信息系统的数据很感兴趣,所以这引起了我的兴趣。我在我的fedora19系统上安装了
liblas
及其依赖项,并使用了liblas
附带的示例数据文件。在在使用你的代码时,我也遇到了同样的问题,看着我所有的记忆都被吞噬了。我不知道为什么会发生这种情况——也许是不必要的引用挂在周围,阻止垃圾收集器像我们希望的那样工作。这可能是可以修复的,但我不会尝试。在
我确实注意到了
liblas
模块的一些有趣的特性,并决定尝试它们。我相信你能得到你想要的数据。在打开文件后,查看头中的XML描述。在
很难看到(随意玩儿xml.dom.minidom或者lxml.etree),但在我的示例文件中,它显示了点数据的字节布局(我的也有28个字节)。在我的例子中,偏移量18是分配给点源ID的一个短(2个字节)。您应该能够用}来检索这个值。不幸的是,
^{pr2}$p.data[18:19]
、p.get_data()[18:19]
、p.point_source_id
、或{data
引用耗尽了内存,p.point_source_id
有一个bug(bug修复请求提交给开发人员)。如果我们改变你的代码使用最后一个访问方法,一切似乎都很好。因此,请在for
循环中尝试以下操作:请注意
如果你想设置唯一的ID值。。。在
希望您的数据值也可用作
p.get_point_source_id()
。请在评论中告诉我它是如何为你工作的。干杯!在相关问题 更多 >
编程相关推荐