1

最近我明白了使用列表理解的巨大优势。我正在处理以特殊格式 *.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 中，我可以将列表理解用作：

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'

提前感谢您的帮助。詹尼

4

3 回答 3

3

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

于 2012-10-09T18:01:57.567 回答
1

尝试

import numpy as np
...
points = np.array([f.X, f.Y]).T
于 2012-10-09T18:26:42.600 回答
1

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

for p in f: p.make_nice()

现在您的列表组合应该可以工作（使用大写的 X 和 Y——请参阅下面的评论）。

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

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

相关来源：

def make_nice(self):
    '''Turn a point instance with the bare essentials (an unpacked list of data)
    into a fully populated point. Add all the named attributes it possesses, 
    including binary fields.
    '''
    i = 0
    for dim in self.reader.point_format.specs: 
            self.__dict__[dim.name] = self.unpacked[i]
            i += 1

    # rest of method snipped
于 2012-10-09T18:27:38.613 回答