两种ICP的改进算法：PLICP与NICP

作者：小L

编辑：3D视觉工坊

前言

在之前的文章中（ICP方法详细推导)，我们介绍了ICP的基本思想与详细的推导。本文将介绍ICP方法的两种改进，分别是：PLICP[1]与NICP[2]。本文将分别介绍两种改进的基本思想，具体算法以及一些补充说明。若有理解不到位和错误之处，请以论文原文为准。

第一部分 PLICP

一、基本思想

PLICP中的“PL”表示”Point to Line”，顾名思义，在匹配时是一个点与一个直线进行匹配，而不是传统方法的点与点进行匹配。之所以有这种思想，是因为我们认为每次扫描的数据是对真实物理世界的一个平面的采样，所以我们在匹配时应该尽可能与这个直线去匹配而不是具体的采样点。

（左图棕色曲线表示真实的物理面，蓝色的为带有噪声的采样点；中间表示传统ICP的点点距离，右图表示PLICP方法，匹配时是计算到平面的距离）

二、 算法描述

2.1 利用上一次迭代的变化参数（或初值），对当前采样（curr）的每个点进行变化；

2.2 变换后，寻找每个点在参考点云（ref）中的最近邻的两个点

2.3 使用论文[3]中提到的方法，去除离群点；

2.4 构建目标函数：

三、补充说明

3.1 论文[3]介绍了一种截断剔除离群点的方法，具体而言，在完成两组点云的匹配后，计算每组匹配点的欧氏距离，只保留距离最小的一定百分比的匹配，从而对噪声鲁棒；

3.2 对算法中的目标函数进一步解释：可以看出点最近直线的距离，采用了投影的思想，法向量点乘即为在法向量上的投影；

3.3 算法中对目标函数求解最佳变换参数的方法有很多，论文给出了一种闭式解的方法，具体请参考论文[1]的附录，这里不再展开。

3.4 PLICP相比于ICP而言，收敛速度更快（论文证明，ICP是一阶收敛，而PLICP是二阶收敛）。但更容易陷入局部极值，故一般使用时，多采用全局ICP方法（例如论文采用了GPM[4]）进行粗匹配，然后在使用PLICP进行精确计算。

3.5 作者给出了PLCIP方法的C语言实现：https://censi.science/software/csm/

第二部分NICP

一、基本思想

NICP的基本思想是，curr和ref的两个点在匹配时，不仅要距离接近，而且所在处的法向量方向也要相同。在匹配时，根据距离、曲率以及法向量进行筛选，并在优化变换参数时优化增加了法向量的参数。

（左侧图片表示采样点的曲率，越大的区域为红色；右图表示匹配，绿色和蓝色为两次扫描，紫红色线表示匹配的点，可以看出，右上角部分虽然蓝色和绿色的点在距离上重合，但由于法向量不同，并不会建立匹配关系）

二、算法描述

2.1 法向量计算方法

2.2 匹配原则

在进行匹配时，不同于传统ICP，距离大于一定阈值时剔除，NICP采用3个准则剔除错误匹配，分别是：1. 点距离超过阈值；2. 曲率接近；3. 法向量方向接近。

2.3 目标函数

2.4 优化求解

优化求解可以采用任何优化求解方法。论文采用了LM算法。

三、补充说明

1）NICP采用法向量进行扩充，包含了一定的语义成分；

2）对于协方差矩阵意义的个人理解

协方差矩阵求逆获得了信息矩阵，在对匹配点误差进行加权时，对不同方向上的误差进行了不同权重的约束。例如，如果某个点在一个平面上，那么对应的协方差矩阵的特征值最小值代表了法向量上的“厚度”，最小值越小，表示越接近于平面，那么在信息矩阵中对应的位置权重越大，放大了ICP时在法向量方向上的误差，避免在法向量方向上产生严重的“错位”。

备注：在公众号「3D视觉工坊」后台，回复「ICP」，即可获得以下四篇的参考文献合集。

参考文献

[1]. A. Censi, "An ICP variant using a point-to-line metric," 2008 IEEE International Conference on Robotics and Automation, Pasadena, CA, 2008, pp. 19-25, doi: 10.1109/ROBOT.2008.4543181.

[2] J. Serafin and G. Grisetti, "NICP: Dense normal based point cloud registration," 2015 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Hamburg, 2015, pp. 742-749, doi: 10.1109/IROS.2015.7353455.

[3]. D. Chetverikov, D. Svirko, D. Stepanov and P. Krsek, "The Trimmed Iterative Closest Point algorithm," Object recognition supported by user interaction for service robots, Quebec City, Quebec, Canada, 2002, pp. 545-548 vol.3, doi: 10.1109/ICPR.2002.1047997.

[4]. A. Censi, "Scan matching in a probabilistic framework," Proceedings 2006 IEEE International Conference on Robotics and Automation, 2006. ICRA 2006., Orlando, FL, 2006, pp. 2291-2296, doi: 10.1109/ROBOT.2006.1642044.

本文仅做学术分享，如有侵权，请联系删文。

下载1

在「3D视觉工坊」公众号后台回复：3D视觉，即可下载 3D视觉相关资料干货，涉及相机标定、三维重建、立体视觉、SLAM、深度学习、点云后处理、多视图几何等方向。

下载2

在「3D视觉工坊」公众号后台回复：3D视觉github资源汇总，即可下载包括结构光、标定源码、缺陷检测源码、深度估计与深度补全源码、点云处理相关源码、立体匹配源码、单目、双目3D检测、基于点云的3D检测、6D姿态估计汇总等。

下载3

在「3D视觉工坊」公众号后台回复：相机标定，即可下载独家相机标定学习课件与视频网址；后台回复：立体匹配，即可下载独家立体匹配学习课件与视频网址。

重磅！3DCVer-学术论文写作投稿 交流群已成立

扫码添加小助手微信，可申请加入3D视觉工坊-学术论文写作与投稿 微信交流群，旨在交流顶会、顶刊、SCI、EI等写作与投稿事宜。

同时也可申请加入我们的细分方向交流群，目前主要有3D视觉、CV&深度学习、SLAM、三维重建、点云后处理、自动驾驶、CV入门、三维测量、VR/AR、3D人脸识别、医疗影像、缺陷检测、行人重识别、目标跟踪、视觉产品落地、视觉竞赛、车牌识别、硬件选型、学术交流、求职交流等微信群。

一定要备注：研究方向+学校/公司+昵称，例如：”3D视觉 + 上海交大 + 静静“。请按照格式备注，可快速被通过且邀请进群。原创投稿也请联系。

▲长按加微信群或投稿

▲长按关注公众号

3D视觉从入门到精通知识星球：针对3D视觉领域的知识点汇总、入门进阶学习路线、最新paper分享、疑问解答四个方面进行深耕，更有各类大厂的算法工程人员进行技术指导。与此同时，星球将联合知名企业发布3D视觉相关算法开发岗位以及项目对接信息，打造成集技术与就业为一体的铁杆粉丝聚集区，近2000星球成员为创造更好的AI世界共同进步，知识星球入口：

学习3D视觉核心技术，扫描查看介绍，3天内无条件退款

圈里有高质量教程资料、可答疑解惑、助你高效解决问题