Rotational Projection Statistics for 3D Local Surface Description and Object Recognition

        本博客是针对论文《Rotational Projection Statistics for 3D Local Surface Description and Object Recognition》的一篇学习笔记。该论文介绍了一种在存在噪声和遮挡的情况下进行特征提取并创建描述符Rotational Projection Statistics (RoPS)的方法。该方法包括三个模块：定义局部参考帧（LRF），RoPS特征描述和3D对象识别。其中实现部分作者分别采用如下方法：

  1. 通过计算位于局部表面上的所有点的散射矩阵来定义LRF
  2. 通过将特征点的相邻点旋转投影到2D平面上并计算这些投影点的分布的一组统计（包括低阶中心矩和熵）来获得RoPS特征描述符。
  3. 使用所提出的LRF和RoPS描述符建立了分层三维物体识别算法。

一、物体识别

        目的：正确识别场景中存在的物体并恢复它们的姿势（即位置和方向）。

        存在问题：在有噪声、网格分辨率变换、遮挡和杂乱的情况下识别来自距离图像的对象仍然是一项具有挑战性的任务。

        目前基于3D对象识别的算法主要包括：基于全局特征和基于局部特征。 基于全局特征的算法构造了一组特征，这些特征对整个3D对象的几何属性进行编码。结合全局特征来进行3D对象识别主要是使用深度学习来完成，例如PointNet、PointNet++、PointCNN等方法。而基于局部特征的算法定义了一组特征，这些特征编码特征点的局部邻域的特征。基于局部特征的算法对于遮挡和杂乱是鲁棒的，因此可以很好地进行物体识别。

        目前提出的基于局部特征的算法论文主要有Point Signature、Using spin images for efficient object recognition in cluttered 3D scenes、Three-dimensional model-based object recognition and segmentation in cluttered scenes、3D geometric scale variability in range images: Features and descriptors。这些算法都包括三个阶段：特征匹配、假设生成和验证、姿态细化。对于局部特征描述符，根据是否使用本地参考帧（LRF）可以分为两类：（没有任何LRF的特征描述符情况下）使用直方图、局部几何信息（例如，法线，曲率）的统计来形成特征描述符。

        这篇论文中提出的Rotational Projection Statistics (RoPS)方法是一种用于局部特征描述和物体识别的方法，通过特征分解和符号重构来完成的LRF构造。

二、局部曲面特征描述

（一）Local Reference Frame

        作者在此通过对位于局部表面上的所有点的散射矩阵执行特征值分解来确定LRF的三个轴。 通过将方向与点散射的区域对齐来消除每个轴的符号。作者在这里用球体的局部表面网格和网格中的点来进行的例证：

        给定特征点p和支撑半径r，对于以p为中心、半径为r的球体，从距离图像中裁剪包含N个三角形和M个顶点的局部表面网格S。 对于具有顶点p_{i1}，p_{i2}和p_{i3}的第i个三角形，位于三角形内的任意一点可以表示为：p_{i}(s,t)=p_{i1}+s(p_{i2}-p_{i1})+t(p_{i3}-p_{i1})    ①。 

        位于第i个三角形内的所有点的散射矩阵C_{i}可以计算为：②。将①代入公式②中可以得到散射矩阵C_{i}表示为：③。对于整个局部表面网格S，它的总散射矩阵C被计算为所有三角形的散射矩阵的加权和，即：④，其中N是S中三角形的个数，w_{i1}是第i个三角形的面积和局部表面S面积之间的比率，即：⑤，其中“×”代表叉积的运算；w_{i2}是与特征点到第i个三角形的质心的距离相关的权重，即：⑥。在这里，第一权重w_{i1}将会提升LRF在面对变换的网格分辨率的表现，而第二权重w_{i2}则是用来改善LRF面对遮挡和杂波的滤波性能。

        为了构造LRF，在这里作者首先对散射矩阵进行特征值分解，即：CV=EV ⑦。其中，E是矩阵C的特征值\left \{ \lambda_{1},\lambda_{2},\lambda_{3}\right \}的对角矩阵，V包含着三个正交向量\left \{ v_{1},v_{2},v_{3}\right \}，且特征向量对应的特征值从大到小排列。

        为了更好的构造LRF（由于局部曲面的对称性，所以这里只用特征分解的方法并不能消除LRF的符号模糊性），作者在这里又引入了一种新的符号消歧技术，在下面会剖析一下作者所引出的符号消歧技术。

（二）Sign Disambiguation

        为了消除LRF的符号模糊性，每个特征向量应指向散射矢量的主方向（从特征点开始，并指向位于局部表面上的点的方向）。 因此，每个特征向量的符号由特征向量的内积和散射向量的符号确定。因此，在这里作者采用将特征向量v_{1}模糊定义为：的方法。其中，函数sign(h)代表提取实数的符号，而h的计算为：

        同样，将特征向量v_{3}定义为：。在这里将v_{2}定义为\widetilde{v}_{3}\times \widetilde{v}_{1}。这里，p是原点，v1，v2和v3分别是x，y和z轴。此时生成LRF，使用此LRF可以生成唯一、姿势不变且高辨别性的局部特征描述符。

（三）对LRF的优化（实验部分）

        作者为了评估提出的LRF的可重复性和鲁棒性，计算了场景和模型中点对之间的LRF误差。

        实验中使用的六个场景是通过对原始的网格数据降采样了\frac{1}{2}，然后将高斯噪声和0.1的网格分辨率标准差加到原始数据中。在每个模型中随机选择1,000个点，并将这些点称为特征点。然后通过搜索与模型中的特征点的距离最小的点来获得场景中的对应点。对于每个点对\left ( p_{Si} ,p_{Mi}\right )，我们计算两个点之间的LRF，并将值定义为L_{Si}和L_{Mi}。对于第\large i个点对两个LRF之间的误差，作者提出的计算公式为：，其中\LARGE \epsilon _{i}代表两个LRF之间的旋转误差量，在没有错误的情况下为零。

        作者通过五次随机选择1,000个不同的点对来测试每个方法的总LRF，每种技术的总LRF误差显示为直方图。

       图中部分LRF误差大于80°，主要是由于场景中存在局部对称表面（例如，平坦或球形表面）。 对于局部对称曲面，由于其分布，其LRF存在固有的符号模糊性，因为点的分布在所有方向上几乎相同。 为了处理这种情况，作者采用了一种特征点选择技术，该技术使用特征值的比率来避免局部对称的表面。确定LRF后，下一步是定义局部表面描述符。

（四）RoPS描述符

        作者对于RoPS描述符通过执行局部表面旋转、相邻点投影和统计计算进行求解。

        作者通过上图给出了整个RoPS方法的说明。首先，图a和图b分别给出了模型和局部表面重构的网格图像。由于已经为p定义了LRF，并且将局部表面S中三角形的顶点构成了一个点云\large Q=\left \{ q_{1},q_{2},......,q_{M} \right \}，然后相对于Q进行LRF变换以便实现旋转不变性，新的点云为\large Q'=\left \{ q'_{1},q'_{2},......,q'_{M} \right \}。将点云围绕x轴旋转角度\large \theta _{k}，从而产生旋转的点云\large Q\left (\theta _{k} \right )。然后将该点云\large Q\left (\theta _{k} \right )投影到三个坐标平面（即xy，xz和yz平面）上以获得三个投影点云\large Q_{i}\left (\theta _{k} \right )，i = 1,2,3，如图d。对投影点云\large Q_{i}\left (\theta _{k} \right )，获取其2D边界矩阵，然后计数落入每个箱中的点的数量以产生L×L矩阵D，如图e。矩阵D称为“归属矩阵”，它表示相邻点的2D分布。分布矩阵D被进一步归一化，使得所有区间的总和等于1，以便实现网格分辨率的变化的不变性。

        为了进一步压缩分布矩阵D中的信息以便实现计算和存储效率，作者从分布矩阵D中提取了一组统计数据，包括中心矩和香农熵，中心矩被用于它们的数学简单性和丰富的描述性，而香农熵被选择用于测量概率分布中包含的信息的功能。矩阵D的阶m + n的中心矩\large \mu _{mn}可以被定义为:，其中，。香农熵e计算为：。这些选定的中心矩与香农熵一起用于形成统计向量，即图f。然后将来自xy，xz和yz平面的三个统计矢量连接起来以形成子特征\large f_{x}\left ( \theta _{k} \right )，表示围绕x轴的第k次旋转的总统计量，即图g。然后通过将所有旋转的子特征连接成向量来生成整体特征描述符，即：

        RoPs描述符有4个参数：

  1. 统计数的组合（文章中通过列举不同的组合找到了最佳的组合方式）
  2. 分区L的数目（本文中推荐设置为5）
  3. 每个坐标系中的旋转系数T（本文中推荐设置为3）
  4. 支撑半径r（取15mr）

        对于这些参数的取值作者在论文中有所讲解，在这里就不一一陈述了。

        描述符通常通过准确率和召回率来进行评估，召回率的定义为：，精准度的定义为：。通过改变阈值，可以生成RP曲线。理想情况下，RP曲线将落在图的左上角，这意味着该特征可以获得高召回率和精度。

（五）物体识别流程图

 

id="article_content" class="article_content clearfix"> id="content_views" class="markdown_views prism-tomorrow-night">

在这里插入图片描述

class="table-box">
🚀 算法题 🚀

🌲 算法刷题专栏 | 面试必备算法 | 面试高频算法 🍀
🌲 越难的东西,越要努力坚持，因为它具有很高的价值，算法就是这样✨
🌲 作者简介：硕风和炜，CSDN-Java领域优质创作者🏆，保研|国家奖学金|高中学习JAVA|大学完善JAVA开发技术栈|面试刷题|面经八股文|经验分享|好用的网站工具分享💎💎💎
🌲 恭喜你发现一枚宝藏博主,赶快收入囊中吧🌻
🌲 人生如棋，我愿为卒，行动虽慢，可谁曾见我后退一步？🎯🎯

class="table-box">
🚀 算法题 🚀

在这里插入图片描述
在这里插入图片描述

class="toc">

🍔 目录

🚩 题目链接

⛲ 题目描述

请你来实现一个 myAtoi(string s) 函数，使其能将字符串转换成一个 32 位有符号整数。

函数 myAtoi(string s) 的算法如下：

  • 空格：读入字符串并丢弃无用的前导空格（" "）
  • 符号：检查下一个字符（假设还未到字符末尾）为 ‘-’ 还是 ‘+’。如果两者都不存在，则假定结果为正。
  • 转换：通过跳过前置零来读取该整数，直到遇到非数字字符或到达字符串的结尾。如果没有读取数字，则结果为0。
  • 舍入：如果整数数超过 32 位有符号整数范围 [−231, 231 − 1] ，需要截断这个整数，使其保持在这个范围内。具体来说，小于 −231 的整数应该被舍入为 −231 ，大于 231 − 1 的整数应该被舍入为 231 − 1 。
    返回整数作为最终结果。

示例 1：

输入：s = “42”

输出：42

解释：加粗的字符串为已经读入的字符，插入符号是当前读取的字符。

带下划线线的字符是所读的内容，插入符号是当前读入位置。
第 1 步：“42”（当前没有读入字符，因为没有前导空格）
^
第 2 步：“42”（当前没有读入字符，因为这里不存在 ‘-’ 或者 ‘+’）
^
第 3 步：“42”（读入 “42”）
^
示例 2：

输入：s = " -042"

输出：-42

解释：

第 1 步：" -042"（读入前导空格，但忽视掉）
^
第 2 步：" -042"（读入 ‘-’ 字符，所以结果应该是负数）
^
第 3 步：" -042"（读入 “042”，在结果中忽略前导零）
^
示例 3：

输入：s = “1337c0d3”

输出：1337

解释：

第 1 步：“1337c0d3”（当前没有读入字符，因为没有前导空格）
^
第 2 步：“1337c0d3”（当前没有读入字符，因为这里不存在 ‘-’ 或者 ‘+’）
^
第 3 步：“1337c0d3”（读入 “1337”；由于下一个字符不是一个数字，所以读入停止）
^
示例 4：

输入：s = “0-1”

输出：0

解释：

第 1 步：“0-1” (当前没有读入字符，因为没有前导空格)
^
第 2 步：“0-1” (当前没有读入字符，因为这里不存在 ‘-’ 或者 ‘+’)
^
第 3 步：“0-1” (读入 “0”；由于下一个字符不是一个数字，所以读入停止)
^
示例 5：

输入：s = “words and 987”

输出：0

解释：

读取在第一个非数字字符“w”处停止。

提示：

0 <= s.length <= 200
s 由英文字母（大写和小写）、数字（0-9）、’ ‘、’+‘、’-’ 和 ‘.’ 组成

🌟 求解思路&实现代码&运行结果

⚡ 模拟

🥦 求解思路
  1. 处理空白字符：跳过字符串开头的所有空白字符（如空格、制表符等）。
  2. 处理正负号：读取可选的正负号，并记录符号。
  3. 读取数字：从第一个非空白字符开始，读取连续的数字字符，并将其转换为整数。
  4. 处理溢出：在转换过程中，检查是否超出 32 位有符号整数范围。如果超出，则返回边界值。
  5. 返回结果：根据符号返回最终的整数。
  6. 有了基本的思路，接下来我们就来通过代码来实现一下。
🥦 实现代码
class Solution {
    public int myAtoi(String s) {
        if (s == null || s.length() == 0) {
            return 0; // 空字符串直接返回 0
        }

        int index = 0; // 当前字符的索引
        int sign = 1; // 符号，默认为正
        int result = 0; // 转换结果

        // 1. 跳过空白字符
        while (index < s.length() && s.charAt(index) == ' ') {
            index++;
        }

        // 2. 处理正负号
        if (index < s.length() && (s.charAt(index) == '+' || s.charAt(index) == '-')) {
            sign = s.charAt(index) == '-' ? -1 : 1;
            index++;
        }

        // 3. 读取数字字符
        while (index < s.length() && Character.isDigit(s.charAt(index))) {
            int digit = s.charAt(index) - '0'; // 当前字符对应的数字

            // 4. 处理溢出
            if (result > Integer.MAX_VALUE / 10 || 
                (result == Integer.MAX_VALUE / 10 && digit > Integer.MAX_VALUE % 10)) {
                return sign == 1 ? Integer.MAX_VALUE : Integer.MIN_VALUE;
            }

            result = result * 10 + digit; // 更新结果
            index++; // 移动到下一个字符
        }

        // 5. 返回结果
        return sign * result;
    }
}
 class="hljs-button signin active" data-title="登录复制" data-report-click="{"spm":"1001.2101.3001.4334"}">
class="hide-preCode-box">
🥦 运行结果

在这里插入图片描述

💬 共勉

class="table-box">
最后，我想和大家分享一句一直激励我的座右铭，希望可以与大家共勉！

在这里插入图片描述

在这里插入图片描述

data-report-view="{"mod":"1585297308_001","spm":"1001.2101.3001.6548","dest":"https://blog.csdn.net/Coder_ljw/article/details/145582156","extend1":"pc","ab":"new"}">>
注：本文转载自blog.csdn.net的我很懒但我很软乎的文章"https://blog.csdn.net/u014509577/article/details/88909050"。版权归原作者所有，此博客不拥有其著作权，亦不承担相应法律责任。如有侵权，请联系我们删除。
复制链接
复制链接
相关推荐
发表评论

评论记录：

未查询到任何数据！