今天要介绍的文章是边佳旺在CVPR2017以及IJCV2020发表的GMS，论文名是”GMS: Grid-based Motion Statistics for Fast, Ultra-robust Feature Correspondence”。同之前介绍的外点滤除算法AdaLAM的目的相同，该算法能够实现对初始匹配的筛选，以减少错误匹配。
(a)图是ORB+Lowe Ratio的结果，有很多错误的匹配；(b)图是ORB+Lowe Ratio+GMS的结果，匹配效果明显变好。 首先祭出代码 Repo：
原有算法复杂，为了达到好的效果，设计的算法比较复杂，导致运行速度慢，性价比低，无法部署在实时性要求较高的系统中。
首先给出第一个假设：
运动平滑性：真实匹配的小邻域内匹配通常是在三维空间中的同一块区域。同样地，一个错误匹配的邻域内的匹配通常是几何上不同的三维位置。
这个假设告诉我们：正确匹配的邻域内有多个支持它的匹配，而错误匹配的邻域内支持它的匹配是很少的。这里其实隐含着一个逻辑：作者通过观察发现，正确匹配的邻域内匹配的数量多的概率就会很大，反之则少；那么根据贝叶斯法则，那些匹配对邻域内有很多匹配的匹配关系有很大概率是正确匹配。 一句话：正确匹配周围会有较多的匹配去支持它，而错误的匹配周围支持它的匹配很少。
输入图像为，它们分别有个特征匹配；
表示图到图的最近邻匹配；
其中表示像素点与的匹配；定义的邻域为：
以及它的相似邻域可以表示为：
其中表示两点之间的欧式距离，和为距离阈值，对于：
它表示在中元素的个数，即匹配的运动支持。 如下图所示：黄圈区域即表示邻域。支持它的匹配对为，错误匹配没有支持它的匹配对，即。那么我们就可以通过统计的方式知道哪些匹配是正确的，哪些是错误的。
由上文我们知道：统计支持域内的匹配数量就能够判别正确/错误匹配。于是我们可以对建模，如下：
其中表示二项分布：
表示邻域内匹配对的数量，和分别表示正确/错误匹配被其某个邻域窗口匹配支持的概率。
进一步解释：是受特征质量的控制，即接近匹配的正确率。通常很小，因为错误匹配几乎随机分布在图像范围内。上式的期望为：
方差为：
于是真假匹配的可区分度被表示为如下形式（均值差异比方差之和）：
上式告诉我们：
如果时，。也就是说，当图像中的匹配点越多，此时正确与错误匹配的可区分度就越强！同时，即使是在比稍微大一丁点的情况下，也可以通过增多特征点的数量来弥补正确匹配与错误匹配区分程度不足的缺点；当然，也可以通过提高特征质量来增加正确/错误匹配之间的区分度。
上图展示了作者在Oxford Affine Dataset上验证模型合理性的示意图。利用SIFT特征进行匹配，根据真值标记出正确以及错误匹配。统计每个匹配所在小邻域内的匹配数量。可以发现，正确匹配的支持域得分明显高于错误匹配，即使在非常困难的匹配序列上，该现象仍然存在。
因此，可通过对设置阈值来判定是正确或者错误匹配：
上式中，以及分别表示正确/错误匹配集合，的阈值被设置为：
其中为超参数，经验值为4~6。
看到这里，大家肯定有一个疑问：如何高效地实现上述算法呢？难道对每一个匹配画个圈圈，然后统计圈圈内的匹配数？当然不是，本文设计了一种基于划分网格的算法对上述算法进行加速。将两幅图像划分为多个不重合的网格：以及，表示为落到网格和中的一个匹配对，如上图中的一个红线段。于是邻域（表示在网格网格中的匹配）被重新定义为：
相似邻域被重新定义为：
上面的两个式子中，表示某个网格，表示落在中的匹配对，表示同时落在和中的匹配对。换句话说，本文将落在同一个网格中的匹配当作邻域，将同时落在两个网格中的匹配称为相似邻域，即cell-pair。
注意：一个cell-pair中的匹配满足运动平滑性假设，所以在判断匹配关系是正确还是错误时，仅需要判定一个cell-pair中的所有匹配是否正确即可，无需逐个匹配判断。此外，不需要确定所有可能的cell-pair，而是只检查与第一个图像网格匹配数量最多的一个cell-pair即可。
如果网格很小，则很少邻域信息将被考虑， 这会降低算法性能。 但是，如果网格很大，则将包括更多不正确的对应关系。为解决这个矛盾，我们将网格大小设置为较小以提高准确性，并提出运动内核以考虑更多邻域信息。
考虑了原有网格邻域的8个网格，构成一个更大的网格。
于是我们从新定义了邻域：
其中； 重新定义相似邻域：
其中：。
为了应对匹配过程中的尺度与旋转问题，本文提出了多尺度以及多旋转策略。
固定一个图像网格，假设网格大小为，改变另一个图像的网格大小为，其中有5个候选的取值，分别为
在不同尺度的网格上运行GMS算法，然后统计最好的匹配结果。如果场景的尺度发生较大改变，此时可以设置更多的候选值或者增大。
利用旋转运动核模拟不同方向的旋转，如下图所示，固定，对按照顺时针旋转，这样可以得到8个运动核。然后利用GMS算法在所有的运动核上，然后选择最好的结果（选择匹配数量最多那个）。
综上，算法伪代码如下：
GMS能够在PC端速度2ms，multi-scale(GMS-S)以及multi-rotation(GMS-R)会增加一定的耗时。 
求解位姿速度快，且位姿精确。 
集成到SLAM也可获得较好的结果（该实验不够充分，仅测试了将GMS集成到SLAM初始化的阶段，为了展示GMS能够更快/更好的进行初始化，同时能够产生较多的地图点）。
matchGMS()，可直接调用；此内容由惯性聚合(RSS阅读器)自动聚合整理，仅供阅读参考。 原文来自 — 版权归原作者所有。