#include <iostream>
#include <pcl/point_types.h>
#include <pcl/io/pcd_io.h>
#include <pcl/features/normal_3d.h>
int main()
{
// 加载点云数据
pcl::PointCloud<pcl::PointXYZ>::Ptr cloud_in(new pcl::PointCloud<pcl::PointXYZ>);
if (pcl::io::loadPCDFile<pcl::PointXYZ>("input.pcd", *cloud_in) == -1)
{
PCL_ERROR("Couldn't read file input.pcd \n");
return (-1);
}
// 创建法线存储对象
pcl::PointCloud<pcl::Normal>::Ptr cloud_normals(new pcl::PointCloud<pcl::Normal>);
// 创建KdTree对象,用于在点云中搜索邻近点
pcl::search::KdTree<pcl::PointXYZ>::Ptr tree(new pcl::search::KdTree<pcl::PointXYZ>());
tree->setInputCloud(cloud_in);
// 创建法线估计对象
pcl::NormalEstimation<pcl::PointXYZ, pcl::Normal> ne;
ne.setInputCloud(cloud_in);
ne.setSearchMethod(tree);
// 设置搜索半径
ne.setRadiusSearch(0.03);
// 计算法线
ne.compute(*cloud_normals);
// 输出法线
std::cout << "Point Cloud normals with " << cloud_normals->points.size() << " data points." << std::endl;
for (size_t i = 0; i < cloud_normals->points.size(); ++i)
{
std::cout << " point " << cloud_in->points[i].x << " "
<< cloud_in->points[i].y << " "
<< cloud_in->points[i].z << " has normal "
<< cloud_normals->points[i].normal_x << " "
<< cloud_normals->points[i].normal_y << " "
<< cloud_normals->points[i].normal_z << std::endl;
}
// 如果你想把带有法线的点云保存为PCD文件
pcl::PointCloud<pcl::PointNormal>::Ptr cloud_with_normals(new pcl::PointCloud<pcl::PointNormal>);
pcl::concatenateFields(*cloud_in, *cloud_normals, *cloud_with_normals);
pcl::io::savePCDFileASCII("output_with_normals.pcd", *cloud_with_normals);
return 0;
}
得到点的法向量存储输出
最新推荐文章于 2024-08-03 18:36:56 发布