jjzjj

Point Cloud Library(PCL开源库)学习一

doll ~CJ 2024-02-08 原文

一、PCL库简介(包含下载与配置方法)

        点云库(PCL)是一个开源算法库,用于点云处理任务和3D几何处理。该库包含用于点云滤波、特征点估计、表面重建、3D配准、模型拟合、对象识别、分割和可视化的算法。PCL库有自己存储点云的数据格式——PCD,但也允许以部分其它格式加载和保存数据集。PCL库是基于C++编写的,并在BSD许可下发布[1]。


        (一)C++版本PCL库下载与配置

        PCL库C++下载与配置方法详见参考资料[4]。

        (二)Python版本PCL库下载与配置

        PCL库anaconda+pycharm+windows下载与配置方法见参考资料[5][6]。其中,值得注意的是在Pycharm添加解释器时,Conda可执行文件为如下:

        在完成Python3.6版本的Anaconda虚拟环境配置与对应版本PCL库的安装后,在命令框内(Command)通过pip install laspy下载laspy外部库以配合PCL库读入和处理LAS格式点云数据,注意将对应库文件复制到虚拟环境的site-packages文件夹下(或直接在CMD.exe Prompt中激活虚拟环境直接pip install)。

        PCL库python下载与配置方法详见参考资料[5][6][7]。

        若python-pcl库配置仍存在问题,换种思路,在不影响项目工程要求的情况下换用pclpy库包。


          常用的点云数据处理工具补充:CloudCompare、Geomagic Studio、LAStools、matlab 、ENVI LiDAR等。

二、点云数据格式介绍(LAS格式文件、PCD格式文件为主)

         常见点云文件格式[12]:

LASLAZ(LiDAR数据的工业标准格式,属于二进制文件格式)

PCD(PCL库官方指定格式)

OBJ(是由Alias|Wavefront Technologies公司从几何学上定义的3D模型文件格式,属于文本文件)

PCAP(Velodyne公司出品的激光雷达默认采集数据的二进制文件格式)

PLY(斯坦福大学的Turk等人设计开发的多边形文件格式,斯坦福三角格式有文本和二进制两种)[5]

PTS(简便点云格式,属于文本格式,只包含点坐标信息,按XYZ顺序存储)

XYZTXT(文本格式)


(一)LAS格式文件

        LAS文件按每条扫描线排列方式存放数据,大致包括激光点的三维坐标、多次回波信息、强度信息、扫描角度、分类信息、飞行航带信息、飞行姿态信息、项目信息、GPS信息、数据点颜色信息等,LAS格式定义中用到的数据类型遵循1999年ANSI(美国国家标准化协会)C语言标准。LAS文件主要由Header、VLRS(变长记录,Variable Length Record)、Point Records三部分组成,LAZ是LAS格式的无损压缩版本。

        LAS文件格式详解见[12][13][19], 基于LAStools工具查看数据信息示例:

                                       

 图一  LAS点云原始数据(带地物标签)                                  图二 对应的光学影像

                 图三 LAS格式文件数据信息

(二)PCD格式文件

        PCD文件以ASCII字符编码,文件格式头(File Format Header)说明文件中存储的点云数据的格式,每个格式声明及点云数据之间用换行符(\n)隔开。PCD文件格式详解见[15]。

        基于记事本查看数据信息示例:

        通过代码格式转换的方式,我们可以控制我们想要的输出数据存储编码类型,而CloudCompare存储输出PCD格式的Data编码类型为Binary(二进制)。

                  

图四   代码转的PCD格式文件数据信息            图五 CloudCompare软件转的PCD格式文件数据信息

三、点云数据格式转换(LAS转PCD,基于Python)

import laspy
import numpy as np
from pclpy import pcl

filePath = r"G:\LiDAR\DATA\copy.las"
las = laspy.read(filePath)
point_format = las.point_format
# 查看点云数据类型信息
print("点云数据信息维名称")
print(list(point_format.dimension_names))
# intensity = las.intensity
# print(max(las.intensity))
# 通过max(las.intensity)==0,已知该las数据未存储有效的回波强度信息
x, y, z = np.float32(las.X), np.float32(las.Y), np.float32(las.Z)
# 构造二维矩阵(行代表各点,列代表点的X、Y、Z),若数据回波强度有效修改点云构造类即可,例如PointXYZI
Pmatrix = np.vstack((x,y,z)).transpose()
cloud = pcl.PointCloud.PointXYZ.from_array(Pmatrix)
pcl.io.savePCDFileASCII(r"G:\LiDAR\DATA\las2pcd.pcd",cloud)

       :利用CloudCompare软件也可以进行多种点云数据文件格式之间的转换,将数据导入后以所需格式保存即可。

四、点云数据的C++显示与python显示

(一)基于pcl.visualization模块点云显示

from pclpy import pcl

filePath = r"G:\LiDAR\DATA\las2pcd.pcd"
# 加载点云数据
cloud = pcl.PointCloud.PointXYZ()
reader = pcl.io.PCDReader()
reader.read(filePath,cloud)
# 可视化点云数据
viewer = pcl.visualization.CloudViewer("viewer")
viewer.showCloud(cloud,"sample cloud")
while not viewer.wasStopped(10):
    pass

(二)基于C++的PCL库点云显示

       ① 在调试运行项目时,dll文件缺失解决方法可查看参考资料[22]。

        ②对于生成解决方案时产生Error:4996,可以选择通过如下方法解决:

        解决方案资源管理器 ->右击项目文件名 ->属性->C/C++ ->高级->禁用特定警告->确定。


         在具体调试测试时,我发现利用pclpy来显示上述代码转换格式数据CloudCompare转换数据均可以成功显示运行;但是在VS 2019上基于PCL库来显示代码转换格式数据就无法成功,会跳到PCL.exe面板卡住(此问题并未解决,还需要进一步学习)。


# include<iostream>
# include<pcl/visualization/cloud_viewer.h>
# include<pcl/io/io.h>
# include<pcl/io/pcd_io.h>
# include<pcl/io/ply_io.h>
# include<pcl/point_types.h>
# include<cstdlib>
# include <vtkAutoInit.h>


using namespace std;
using namespace pcl;
using namespace pcl::io;

int user_data;

void viewerOneOff(pcl::visualization::PCLVisualizer& viewer)
{
	//设置背景颜色
	viewer.setBackgroundColor(1.0, 0.5, 1.0);
	//添加圆球几何对象
	pcl::PointXYZ o;
	o.x = 1.0;
	o.y = 0;
	o.z = 0;
	viewer.addSphere(o, 0.25, "sphere", 0);
}

void viewerPsycho(pcl::visualization::PCLVisualizer& viewer)
{
	static unsigned count = 0;
	std::stringstream ss;
	ss << "Once per viewer loop:" << count++;
	viewer.removeShape("text", 0);
	viewer.addText(ss.str(), 200, 300, "text", 0);
	user_data++;
}

int main()
{
	// 创建点云渲染对象,导入待渲染文件
	pcl::PointCloud<pcl::PointXYZ>::Ptr cloud(new pcl::PointCloud<pcl::PointXYZ>);
	pcl::io::loadPCDFile("G:\\LiDAR\\DATA\\copy - Cloud.pcd",*cloud);
	// 创建点云渲染句柄
	pcl::visualization::CloudViewer viewer("Cloud Viewer");
	viewer.showCloud(cloud);
	//以下代码均为高级复杂操作代码,此处仅依据参考尝试调试查看(若需具体尝试可向其中添加相应代码)
	//该注册函数在可视化的时候只执行一次
	viewer.runOnVisualizationThreadOnce(viewerOneOff);
	//该注册函数在渲染输出时每次都调用
	viewer.runOnVisualizationThread(viewerPsycho);
	while (!viewer.wasStopped())
	{
		user_data++;
	}
	return 0;
}

五、Anaconda虚拟环境管理配置小Tips

        在面对实际项目时,我们通常会遇到需要在不同项目使用不同的python编程环境,如果同时安装多个python环境就可能造成许多混乱和错误。这时利用Anaconda的虚拟环境机制就可以方便安装不同的python版本,在不同的项目中使用不同的运行环境(因为某些项目涉及不同包版本的匹配等问题,如python-PCL库包需要对应支持的版本等),从而方便来回切换运行环境。

        若在具体实践项目时,我们的库包及模块出现报错“未能被找到”,可以利用Denpendency Walker工具检查模块及库包的依赖关系,查找缺少的动态链接库(DLL,Dynamic Link Library),从而具体分析具体解决。

        若下载库包时由于网络问题造成下载速度较慢或下载超时,我们可以采用基于镜像网页的下载方式去下载库包,格式:pip install 安装包名 -i 镜像网页。

        具体操作见参考资料[17][19]。

参考资料:

[1] https://en.wikipedia.org/wiki/Point_Cloud_Library(wiki)

[2]Introduction — Point Cloud Library 0.0 documentation(PCL库使用指南)

[3]Point Cloud Library | The Point Cloud Library (PCL) is a standalone, large scale, open project for 2D/3D image and point cloud processing.(PCL官网)

[4]win10系统下 VS2019点云库PCL1.12.0的安装与配置_苦瓜王子的博客-CSDN博客_pcl库下载

[5]【python3.6】python安装PCL(适用命令行或pycham中)_米码收割机的博客-CSDN博客_python pcl

[6](详细安装python_pcl) python_pcl+windows+anaconda_没李不邢的博客-CSDN博客_python3.7 pcl

[7](详细PCL安装)PCL+python+windows+anaconda环境_没李不邢的博客-CSDN博客_anaconda配置pcl

[8]Large Geometric Models Archive(几何模型文件库)

[9]Point Cloud Library (PCL): Modules(PCL模块函数库)

[10]PCL(Point Cloud Library)学习记录(2022) · 语雀(PCL大佬的自建网站)

[11]PCL(Point Cloud Library)学习入门指南&代码实践(最新版)_哔哩哔哩_bilibili

[12]点云Las文件是什么?_程序媛一枚~的博客-CSDN博客_las文件

[13]点云数据格式说明_scott198510的博客-CSDN博客_点云数据格式

[14]Python,C++中点云 .las转.pcd_程序媛一枚~的博客-CSDN博客

[15]PCD文件格式(记录) - 哔哩哔哩

PCL——(3)PCD(点云数据)文件格式简介 - 一抹烟霞 - 博客园

[16]数据来源:美国休斯顿大学高光谱实验室(采集时间:2017年)

[17][转]使用Dependency Walker工具查看DLL文件_mahui85的博客-CSDN博客_查看dll(查找库包之间的依赖关系,用于针对module报错)

[18]ImportError:DLL load failed:找不到指定的模块 - CodeAntenna

[19]pyCharm下配置Anaconda虚拟环境_IT大佬啊的博客-CSDN博客_pycharm anaconda虚拟环境

[20]API Documentation — laspy 2.3.0 documentation

[21]https://github.com/maguangyan/pclpy_tutorial(pclpy大佬教程)

[22]vs2017+vtk由于找不到**.dll,无法执行代码。重新安装程序可能会解决此问题。_CHENxiaoni_的博客-CSDN博客

国内的优质镜像网页

[1](清华镜像) https://pypi.tuna.tsinghua.edu.cn/simple/

[2](豆瓣镜像)http://pypi.doubanio.com/simple/    (http)使用时需加入网页信任:--trusted-host pypi.doubanio.com

[3](阿里镜像)http://mirrors.aliyun.com/pypi/simple/

[4](中国科学技术大学镜像)https://pypi.mirrors.ustc.edu.cn/simple/

有关Point Cloud Library(PCL开源库)学习一的更多相关文章

  1. LC滤波器设计学习笔记(一)滤波电路入门 - 2

    目录前言滤波电路科普主要分类实际情况单位的概念常用评价参数函数型滤波器简单分析滤波电路构成低通滤波器RC低通滤波器RL低通滤波器高通滤波器RC高通滤波器RL高通滤波器部分摘自《LC滤波器设计与制作》,侵权删。前言最近需要学习放大电路和滤波电路,但是由于只在之前做音乐频谱分析仪的时候简单了解过一点点运放,所以也是相当从零开始学习了。滤波电路科普主要分类滤波器:主要是从不同频率的成分中提取出特定频率的信号。有源滤波器:由RC元件与运算放大器组成的滤波器。可滤除某一次或多次谐波,最普通易于采用的无源滤波器结构是将电感与电容串联,可对主要次谐波(3、5、7)构成低阻抗旁路。无源滤波器:无源滤波器,又称

  2. CAN协议的学习与理解 - 2

    最近在学习CAN,记录一下,也供大家参考交流。推荐几个我觉得很好的CAN学习,本文也是在看了他们的好文之后做的笔记首先是瑞萨的CAN入门,真的通透;秀!靠这篇我竟然2天理解了CAN协议!实战STM32F4CAN!原文链接:https://blog.csdn.net/XiaoXiaoPengBo/article/details/116206252CAN详解(小白教程)原文链接:https://blog.csdn.net/xwwwj/article/details/105372234一篇易懂的CAN通讯协议指南1一篇易懂的CAN通讯协议指南1-知乎(zhihu.com)视频推荐CAN总线个人知识总

  3. 深度学习部署:Windows安装pycocotools报错解决方法 - 2

    深度学习部署:Windows安装pycocotools报错解决方法1.pycocotools库的简介2.pycocotools安装的坑3.解决办法更多Ai资讯:公主号AiCharm本系列是作者在跑一些深度学习实例时,遇到的各种各样的问题及解决办法,希望能够帮助到大家。ERROR:Commanderroredoutwithexitstatus1:'D:\Anaconda3\python.exe'-u-c'importsys,setuptools,tokenize;sys.argv[0]='"'"'C:\\Users\\46653\\AppData\\Local\\Temp\\pip-instal

  4. ruby - 我正在学习编程并选择了 Ruby。我应该升级到 Ruby 1.9 吗? - 2

    我完全不是程序员,正在学习使用Ruby和Rails框架进行编程。我目前正在使用Ruby1.8.7和Rails3.0.3,但我想知道我是否应该升级到Ruby1.9,因为我真的没有任何升级的“遗留”成本。缺点是什么?我是否会遇到与普通gem的兼容性问题,或者甚至其他我不太了解甚至无法预料的问题? 最佳答案 你应该升级。不要坚持从1.8.7开始。如果您发现不支持1.9.2的gem,请避免使用它们(因为它们很可能不被维护)。如果您对gem是否兼容1.9.2有任何疑问,您可以在以下位置查看:http://www.railsplugins.or

  5. ruby - 我如何学习 ruby​​ 的正则表达式? - 2

    如何学习ruby​​的正则表达式?(对于假人) 最佳答案 http://www.rubular.com/在Ruby中使用正则表达式时是一个很棒的工具,因为它可以立即将结果可视化。 关于ruby-我如何学习ruby​​的正则表达式?,我们在StackOverflow上找到一个类似的问题: https://stackoverflow.com/questions/1881231/

  6. 深度学习12. CNN经典网络 VGG16 - 2

    深度学习12.CNN经典网络VGG16一、简介1.VGG来源2.VGG分类3.不同模型的参数数量4.3x3卷积核的好处5.关于学习率调度6.批归一化二、VGG16层分析1.层划分2.参数展开过程图解3.参数传递示例4.VGG16各层参数数量三、代码分析1.VGG16模型定义2.训练3.测试一、简介1.VGG来源VGG(VisualGeometryGroup)是一个视觉几何组在2014年提出的深度卷积神经网络架构。VGG在2014年ImageNet图像分类竞赛亚军,定位竞赛冠军;VGG网络采用连续的小卷积核(3x3)和池化层构建深度神经网络,网络深度可以达到16层或19层,其中VGG16和VGG

  7. 机器学习——时间序列ARIMA模型(四):自相关函数ACF和偏自相关函数PACF用于判断ARIMA模型中p、q参数取值 - 2

    文章目录1、自相关函数ACF2、偏自相关函数PACF3、ARIMA(p,d,q)的阶数判断4、代码实现1、引入所需依赖2、数据读取与处理3、一阶差分与绘图4、ACF5、PACF1、自相关函数ACF自相关函数反映了同一序列在不同时序的取值之间的相关性。公式:ACF(k)=ρk=Cov(yt,yt−k)Var(yt)ACF(k)=\rho_{k}=\frac{Cov(y_{t},y_{t-k})}{Var(y_{t})}ACF(k)=ρk​=Var(yt​)Cov(yt​,yt−k​)​其中分子用于求协方差矩阵,分母用于计算样本方差。求出的ACF值为[-1,1]。但对于一个平稳的AR模型,求出其滞

  8. Unity Shader 学习笔记(5)Shader变体、Shader属性定义技巧、自定义材质面板 - 2

    写在之前Shader变体、Shader属性定义技巧、自定义材质面板,这三个知识点任何一个单拿出来都是一套知识体系,不能一概而论,本文章目的在于将学习和实际工作中遇见的问题进行总结,类似于网络笔记之用,方便后续回顾查看,如有以偏概全、不祥不尽之处,还望海涵。1、Shader变体先看一段代码......Properties{ [KeywordEnum(on,off)]USL_USE_COL("IsUseColorMixTex?",int)=0 [Toggle(IS_RED_ON)]_IsRed("IsRed?",int)=0}......//中间省略,后续会有完整代码 #pragmamulti_c

  9. ruby-on-rails - 这个 C 和 PHP 程序员如何学习 Ruby 和 Rails? - 2

    按照目前的情况,这个问题不适合我们的问答形式。我们希望答案得到事实、引用或专业知识的支持,但这个问题可能会引发辩论、争论、投票或扩展讨论。如果您觉得这个问题可以改进并可能重新打开,visitthehelpcenter指导。关闭9年前。我来自C、php和bash背景,很容易学习,因为它们都有相同的C结构,我可以将其与我已经知道的联系起来。然后2年前我学了Python并且学得很好,Python对我来说比Ruby更容易学。然后从去年开始,我一直在尝试学习Ruby,然后是Rails,我承认,直到现在我还是学不会,讽刺的是那些打着简单易学的烙印,但是对于我这样一个老练的程序员来说,我只是无法将它

  10. jquery - 使用 Rails 3 学习 Ajax 的资源 - 2

    按照目前的情况,这个问题不适合我们的问答形式。我们希望答案得到事实、引用或专业知识的支持,但这个问题可能会引发辩论、争论、投票或扩展讨论。如果您觉得这个问题可以改进并可能重新打开,visitthehelpcenter指导。关闭10年前。有没有学习Ajax(jQuery)和Rails3的好资源?

随机推荐