3D визуализация облака точек отличается - PullRequest
Новая
       68

3D визуализация облака точек отличается

0 голосов
/ 15 февраля 2019

Я пытаюсь выполнить некоторую постобработку в облаке точек. Когда я пытаюсь построить облако точек с помощью Axis3D, я получаю другой вид облака точек, оно растягивается вдоль оси.Когда я пробую другой вид облака точек, я получаю совершенно другую сетку.

Вот мой код и картинка

enter image description here

def main():

    #process first point cloud
    f3data = np.loadtxt(r'c:\ahmed\Waltersdorf_F3.csv', delimiter=',',
                        dtype=[('id', np.str_, 20), ('x1', np.float32), ('x2', np.float32), ('x3', np.float32),
                               ])


    ptcloud_1 = np.vstack((f3data['x1'], f3data['x2'], f3data['x3'])).transpose()



    pc_1 = pcl.PointCloud.PointXYZ(ptcloud_1)
    X = np.vstack((pc_1.x, pc_1.y, pc_1.z)).transpose()

    import matplotlib.pyplot as plt

    import mpl_toolkits.mplot3d as m3d

    ax = m3d.Axes3D(plt.figure())
    ax.scatter3D(*X.T)
    plt.show()

    # pc_1 = pclpy.radius_outlier_removal(pc_1, search_radius=1.0, min_neighbors=4)
    #plotScatterRot3(f3rot)
    pc_type = utils.get_point_cloud_type(pc_1)
    seg = getattr(pcl.segmentation.SACSegmentation, pc_type)()
    seg.setOptimizeCoefficients(True)
    model = getattr(pcl.sample_consensus, "SACMODEL_" + 'PLANE'.upper())
    seg.setModelType(model)
    seg.setMethodType(pcl.sample_consensus.SAC_RANSAC)
    seg.setDistanceThreshold(1)
    seg.setInputCloud(pc_1)
    coefficients = pcl.ModelCoefficients()
    inliers = pcl.PointIndices()
    arr = np.ndarray(shape=(3,1), dtype= np.float32)
    arr[0] = 0
    arr[1] = 0
    arr[2] = 1
    seg.setAxis(arr)
    seg.segment(inliers, coefficients)

    inla = pc_1.xyz[inliers.indices]

    cloud_filtered = pcl.PointCloud.PointXYZ(inla)
    X = np.vstack((cloud_filtered.x, cloud_filtered.y, cloud_filtered.z)).transpose()


    pca = decomposition.PCA(n_components=3)
    pca.fit(X)
    X = pca.transform(X)
    final_cloud = pcl.PointCloud.PointXYZ(X)

    pcl.io.savePCDFile("original_first.txt", final_cloud)

    import mpl_toolkits.mplot3d as m3d
    ax = m3d.Axes3D(plt.figure())
    ax.scatter3D(*X.T)
    plt.show()

    pptk.viewer(X)

Вот набор данных

http://www.mediafire.com/file/f8tz1zbpxvyvko7/Waltersdorf_F3.csv/file

...