我尝试使此代码打开LAS文件进行更改,然后将结果保存到新的LAS文件中。 las文件包含具有坐标(X和Y)和值(如高程的Z)的点。不幸的是,如果我没有将保存部分放入其中,代码可以工作,但是当我在下面这样做时,我收到以下错误:
Traceback (most recent call last):
File "C:\Users\Geri\Desktop\Sync\pythonlas\envisecond.py", line 33, in <module>
listx1=p.append(listx1, float[newx])
AttributeError: 'float' object has no attribute 'append'
enter code here
到目前为止,我知道我需要一个数组来将值保存到代码末尾的las文件中:
import laspy
import laspy.file
import numpy as p
header = laspy.header.Header()
inFile2 = laspy.file.File("C:\\Users\\Geri\\Desktop\\Sync\\pythonlas\\mapperclip\\2clip.las", mode = "r")
inFile3 = laspy.file.File("C:\\Users\\Geri\\Desktop\\Sync\\pythonlas\\mapperclip\\3clip.las", mode = "r")
point_records = inFile2.points
point_records = inFile3.points
t=0
listx1=p.array([], dtype=float)
listy1=p.array([], dtype=float)
listz1=p.array([], dtype=float)
while t < 415:
z=0
q=0
p=0.1
while z==0:
xmin=inFile3.x[t]-p
ymin=inFile3.y[t]-p
xmax=inFile3.x[t]+p
ymax=inFile3.y[t]+p
n=0
for points in inFile2.points:
ax=inFile2.x[n]
ay=inFile2.y[n]
if ax > xmin and ax < xmax and ay < ymax and ay > ymin:
newx = [inFile3.x[t]-((inFile3.x[t]-inFile2.x[n])/2)]
newy = [inFile3.y[t]-((inFile3.y[t]-inFile2.y[n])/2)]
newz = [inFile3.z[t]-((inFile3.z[t]-inFile2.z[n])/2)]
listx1=p.append(listx1, float[newx])
listy1=p.append(listy1, float[newy])
listz1=p.append(listz1, float[newz])
print listx1
print n
n+=1
q+=1
t+=1
else:
n+=1
if q>0:
z+=1
else:
p+=0.1
outfile = laspy.file.File("C:\\Users\\Geri\\Desktop\\Sync\\pythonlas\\mapperclip\\output2.las", mode="w", header=header)
outfile.X = [listx1]
outfile.Y = [listy1]
outfile.Z = [listz1]
outfile.close()
答案 0 :(得分:5)
您在外部p
循环开始时将float
设置为while
:
p=0.1
掩盖了顶部的numpy
导入:
import numpy as p
所以在while
循环中p
不再是模块,它是float
对象,调用p.append()
调用将失败。
为模块或float
值使用其他名称。