-1

我有这个用于在 Las 文件中操作激光雷达点的小代码。目标是从文件 1 中获取 1 个点并在文件 2 中查找最近的点。然后使用这些点的坐标创建一个新点，然后保存。到目前为止我遇到了2个问题，解释如下：

import laspy
import laspy.file
import liblas
from liblas import header
h = header.Header()

inFile2 = laspy.file.File("C:\\Users\\Geri\\Desktop\\Sync\\Sync\\pythonlas\\mapperclip\\2clip.las", mode = "r")
inFile3 = laspy.file.File("C:\\Users\\Geri\\Desktop\\Sync\\Sync\\pythonlas\\mapperclip\\3clip.las", mode = "r")
point_records = inFile2.points
point_records = inFile3.points

t=0
for points in inFile3.points:
    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:          
                                f = liblas.file.File("C:\\Users\\Geri\\Desktop\\Sync\\Sync\\pythonlas\\mapperclip\\proba.las",mode='w', header= h)
                                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)
                                pt = liblas.point.Point(newx, newy, newz)

问题出在这条线之后。我想将 newx,y,z 保存为新 las 文件中新点的 XY 坐标和 z 值。

收到以下错误：

ctypes.Argumenterror: 参数 1: : 错误类型

                                f.write(pt)
                                f.close()
                                print t
                                print newx
                                print newy
                                print newz
                                print inFile2.y[n]
                                print inFile2.z[n]                                
                                print inFile3.z[t]
                                n+=1
                                q+=1
                                t+=1
            else:
                n+=1
        if q>0:            
            z+=1
        else:
            p+=0.1

当我尝试从代码中写入文件时，我得到以下信息：（脚本成功运行 415 次）

Traceback (most recent call last):
   File "........
       newx = inFile3.x[t]-((inFile3.x[t]-inFile2.x[n])/2)

IndexError：索引 416 超出轴 0 的范围，大小为 416

4

1 回答 1

1

所以这是我到目前为止得到的：

对于第一个问题，我只是省略了 liblas 模块并使用 laspy 来保存创建的新文件。

对于第二个问题，我更换了

for points in inFile3.points:

和

 while t < 415:
于 2015-09-30T01:08:21.353 回答