Question

最近我理解使用列表理解的巨大优势。我正在处理以特殊格式* .las文件存储的数百万个点（x，y，z）。在python中，使用这种格式有两种方法：

Liblas module  [http://www.liblas.org/tutorial/python.html][1] (in a C++/Python)
laspy module [http://laspy.readthedocs.org/en/latest/tut_part_1.html][2] (pure Python)

我对liblas有几个问题，我希望测试laspy。

在liblas中我可以使用list comprehension：

from liblas import file as lasfile
f = lasfile.File(inFile,None,'r') # open LAS
points = [(p.x,p.y) for p in f] # read in list comprehension

在laspy中我无法说明如何做同样的事情：

from laspy.file import File
f = file.File(inFile, mode='r')
f
<laspy.file.File object at 0x0000000013939080>
(f[0].X,f[0].Y)
(30839973, 696447860)

我尝试了几种组合：

points = [(p.X,p.Y) for p in f]

但是我收到了这条消息

Traceback (most recent call last):
 File "<interactive input>", line 1, in <module>
AttributeError: Point instance has no attribute 'x'

我尝试使用大写和非大写，因为Python区分大小写：

>>> [(p.x,p.y) for p in f]
Traceback (most recent call last):
  File "<interactive input>", line 1, in <module>
AttributeError: Point instance has no attribute 'x'
>>> [(p.X,p.Y) for p in f]
Traceback (most recent call last):
  File "<interactive input>", line 1, in <module>
AttributeError: Point instance has no attribute 'X'

这是交互式提示：

C:\Python27>python.exe
Python 2.7.3 (default, Apr 10 2012, 23:24:47) [MSC v.1500 64 bit (AMD64)] on win
32
Type "help", "copyright", "credits" or "license" for more information.
>>> from laspy.file import File
>>> inFile="C:\\04-las_clip_inside_area\\Ku_018_class.las"
>>> f = File(inFile, None, 'r')
>>> f
<laspy.file.File object at 0x00000000024D5E10>
>>> points = [(p.X,p.Y) for p in f]
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
AttributeError: Point instance has no attribute 'X'
>>>

列表后面的打印p是：

print dir(p) 
['__doc__', '__init__', '__module__', 'make_nice', 'pack', 'packer', 'reader', 'unpacked']

以循环格式我总是有同样的错误

>>> for p in f:
...     print dir(p)
...     print p.X,p.Y
...     
['__doc__', '__init__', '__module__', 'make_nice', 'pack', 'packer', 'reader', 'unpacked']
Traceback (most recent call last):
  File "<interactive input>", line 3, in <module>
AttributeError: Point instance has no attribute 'X'

使用nneonneo建议的代码

import numpy as np
for p in f:
...   points = np.array([f.X, f.Y]).T

我可以存储在数组中

points
array([[ 30839973, 696447860],
       [ 30839937, 696447890],
       [ 30839842, 696447832],
       ..., 
       [ 30943795, 695999984],
       [ 30943695, 695999922],
       [ 30943960, 695999995]])

但错过了创建列表理解的方法

points = [np.array(p.X,p.Y).T for p in f]
Traceback (most recent call last):
  File "<interactive input>", line 1, in <module>
AttributeError: Point instance has no attribute 'X'

提前感谢您的帮助。詹尼

Answer 1

Python区分大小写。对我来说，看起来你要求属性x，但它应该是一个大写的X。

Answer 2

尝试

import numpy as np
...
points = np.array([f.X, f.Y]).T

Answer 3

看起来Point有一个make_nice()方法可以显示更多属性。

for p in f: p.make_nice()

现在你的列表组合应该有效（大写的X和Y - 见下面的评论）。

[(p.X,p.Y) for p in f]

注意：此答案未经过测试。它基于阅读laspy.util.Point的来源。

Python：使用模块laspy进行列表理解的问题

3 个答案: