pcl之open3d数据与pcl数据相互转换

30 篇文章 16 订阅
#include<pcl/point_types.h>
#include<pcl/io/pcd_io.h>
#include <open3d/Open3D.h>
#include <iostream>


using namespace open3d;
using namespace std;

std::shared_ptr<geometry::PointCloud> LoadPLY(const string &path) {
    auto cloud_ptr = std::make_shared<geometry::PointCloud>();
    open3d::io::ReadPointCloud(path, *cloud_ptr);
    return cloud_ptr;
}

pcl::PointCloud<pcl::PointXYZ>::Ptr CvtO3DToPCL(std::shared_ptr<geometry::PointCloud> pc) {
    pcl::PointCloud<pcl::PointXYZ>::Ptr cloud_ptr(new pcl::PointCloud<pcl::PointXYZ>);
    cloud_ptr->resize(pc->points_.size());
    for (int i = 0; i < pc->points_.size(); i++) {
        cloud_ptr->points[i].x = pc->points_[i].x();
        cloud_ptr->points[i].y = pc->points_[i].y();
        cloud_ptr->points[i].z = pc->points_[i].z();
    }
    return cloud_ptr;
}

std::shared_ptr<geometry::PointCloud> CvtPCLToO3D(pcl::PointCloud<pcl::PointXYZ>::Ptr pc) {

    auto open3d_cloud_ptr = std::make_shared<geometry::PointCloud>();
    open3d_cloud_ptr->points_.resize(pc->points.size());
    for (int i = 0; i < pc->points.size(); i++) {
        open3d_cloud_ptr->points_[i][0] = pc->points[i].x;
        open3d_cloud_ptr->points_[i][1] = pc->points[i].y;
        open3d_cloud_ptr->points_[i][2] = pc->points[i].z;
    }
    return open3d_cloud_ptr;
}


int main() {
    auto o3d_ply = LoadPLY("/home/rvbust/Documents/test1.ply");

    auto pcl_cloud_ptr = CvtO3DToPCL(o3d_ply);
    auto o3d_cloud_ptr = CvtPCLToO3D(pcl_cloud_ptr);

    visualization::DrawGeometries({o3d_cloud_ptr}, "PointCloud");

    return 0;
}

以下是使用pcl拟合曲线并使用open3d可视化的示例代码: ```python import open3d as o3d import numpy as np import pclpy from pclpy import pcl # 生成点云数据 points = np.array([[0, 0, 0], [0.1, 0.1, 0.1], [0.2, 0.2, 0.2], [0.3, 0.3, 0.3], [0.4,0.4, 0.4], [0.5, 0.5, 0.5], [0.6, 0.6, 0.6], [0.7, 0.7, 0.7], [0.8, 0.8, 0.8], [0.9, 0.9, 0.9]]) # 将点云数据转换pcl格式 cloud = pcl.PointCloud.PointXYZ(points) # 使用pcl进行曲线拟合 b_spline = pcl.surface.BSplineSurface() b_spline.setInputCloud(cloud) b_spline.setOrder(3) b_spline.setNumCtrlPoints(10, 1) b_spline.setNumKnots(10, 1) b_spline.setCtrlPoint(0, 0, pcl.PointXYZ(0, 0, 0)) b_spline.setCtrlPoint(1, 0, pcl.PointXYZ(0.1, 0.1, 0.1)) b_spline.setCtrlPoint(2, 0, pcl.PointXYZ(0.2, 0.2, 0.2)) b_spline.setCtrlPoint(3, 0, pcl.PointXYZ(0.3, 0.3, 0.3)) b_spline.setCtrlPoint(4, 0, pcl.PointXYZ(0.4, 0.4, 0.4)) b_spline.setCtrlPoint(5, 0, pcl.PointXYZ(0.5, 0.5, 0.5)) b_spline.setCtrlPoint(6, 0, pcl.PointXYZ(0.6, 0.6, 0.6)) b_spline.setCtrlPoint(7, 0, pcl.PointXYZ(0.7, 0.7, 0.7)) b_spline.setCtrlPoint(8, 0, pcl.PointXYZ(0.8, 0.8, 0.8)) b_spline.setCtrlPoint(9, 0, pcl.PointXYZ(0.9, 0.9, 0.9)) b_spline.initCompute() # 将拟合曲线转换为numpy数组 curve_points = np.zeros((100, 3)) for i in range(100): u = i / 99.0 curve_points[i] = b_spline.getPoint(u, 0) # 将拟合曲线可视化 curve = o3d.geometry.LineSet() curve.points = o3d.utility.Vector3dVector(curve_points) curve.lines = o3d.utility.Vector2iVector(np.array([[i, i+1] for i in range(99)])) o3d.visualization.draw_geometries([curve]) ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值