[发明专利]一种基于轮廓随机采样的局部形状匹配方法

摘要

本发明公开一种基于轮廓随机采样的局部形状匹配方法，包括：采用几何圆半径方法提取待匹配目标轮廓上的所有子轮廓序列和模板轮廓上随机子轮廓；提取目标子轮廓和模板随机子轮廓的最小旋转外接矩形长宽比及轮廓面积特征，进行第一次匹配；对第一次匹配得到的轮廓点，提取其子轮廓基于角度分类的均值距离链码和均值角度链码特征；利用模板随机子轮廓得到的链码特征与第一次匹配后的子轮廓进行二次匹配，以得到少量最优匹配轮廓；计算模板随机子轮廓与匹配子轮廓的坐标变换矩阵，通过该变换矩阵将整个模板轮廓投影到目标轮廓上；更新子轮廓的坐标变换矩阵；识别最优匹配轮廓。本发明能够快速可靠地给出匹配结果，图像配准计算效率高，适应性强。

主权项

1.一种基于轮廓随机采样的局部形状匹配方法，其特征在于，该方法包括如下步骤：A、提取待匹配目标轮廓上的所有子轮廓序列和模板轮廓上随机子轮廓；其中，所述步骤A之前还包括：一、对相机采集得到的灰度图像进行预处理：(1)对灰度图像进行3×3模板的均值滤波，并采用最大类间方差法计算待匹配目标图像f(x,y)的全局最优阈值T，得到二值图像g(x,y)；(2)从待匹配目标图像序列中截取需匹配的目标区域作为模板图像，对该目标区域进行包括滤波、阈值分割在内的预处理，取得该模板图像的二值图像；二、提取二值图像连通区域的轮廓：(1)通过采用链码方法搜索二值图像g(x,y)的轮廓，并将链码编码中的所有点转换为坐标点，即得到一系列逆时针排序的坐标点对；使用极大极小值来限定区域的周长，遍历所有轮廓移除过长或过短的轮廓，保留待匹配的轮廓；(2)采用链码提取模板图像的二值图像的轮廓，即得到模板图像轮廓上的一系列坐标点；对位于模板图像四周边界的轮廓进行剔除，保留模板图像中的目标轮廓；所述步骤A具体包括：A1、采用几何圆半径范围的抽取方法提取待匹配目标轮廓每个坐标点上的子轮廓：首先计算所有轮廓点之间的欧式距离，并保存在相应的vector数组中，然后统计每个轮廓点处的前后小于圆半径长度radius的轮廓点个数，最后遍历轮廓上所有坐标点得到所有的子轮廓，并保存于vector数组；A2、采用随机采样的方法提取模板轮廓上的随机轮廓点几何圆半径构成的子轮廓：通过计算随机轮廓点与前后轮廓点之间的欧式距离，并统计前后小于圆半径长度radius的轮廓坐标点的个数，最终提取出该随机轮廓点上的圆半径子轮廓；所述步骤A2还包括：将模板轮廓的开始和末尾的在圆半径范围的点剔除；B、提取待匹配目标子轮廓和模板随机子轮廓的特征，其中，该特征包括但不限于待匹配目标子轮廓和模板随机子轮廓的最小旋转外接矩形的长宽比和轮廓面积两个特征；所述步骤B中待匹配目标子轮廓和模板随机子轮廓的最小旋转外接矩形的长宽比和轮廓面积两个特征的提取过程如下：B1、计算最小旋转外接矩形的长宽比特征Fratio：对给定的子轮廓坐标点集，按照水平轴逆时针旋转，每旋转一次，查找包围轮廓点集最小面积的矩形minAreaRect，即得到该矩形的长和宽，从而得到长宽比特征：B2、计算轮廓面积特征Farea：将子轮廓的两端连接，并与轮廓部分构成封闭区域，计算该区域的面积；C、进行第一次匹配：设定匹配幅度阈值Tmagnitude，当待匹配目标图像中子轮廓的所述两个特征同时满足在模板中随机子轮廓对应特征的±Tmagnitude幅度内，则认为所述待匹配目标图像中子轮廓为匹配点，从而剔除非匹配轮廓点，其中，所述待匹配目标图像中子轮廓简称为轮廓点；D、对第一次匹配得到的轮廓点，提取基于角度分类的均值距离链码和均值角度链码特征；所述步骤D具体包括：D1、计算最小旋转外接矩形的中心rect_center和轮廓的重心mass_center；通过计算最小旋转外接矩形得到该矩形的中心rect_center和四个角点的坐标rect_points[4]，通过零阶矩和一阶矩得到轮廓的重心，计算公式如下：D2、依据轮廓的重心将矩形中心偏移到矩形外侧，并将偏移中心位于矩形长边一侧：首先通过计算矩形四个角坐标点之间的欧式距离，锁定长边；然后计算矩形的两条长边和过中心点的平行线共三条直线的方程式ajx+bjy+c＝0其中，j＝1,2,3；最后，锁定重心位于过中心点的平行线的一侧方向，作该中心点的对称点offset，即完成对中心点的偏移；D3、以偏移点为中心，偏移点与轮廓起始点为边，计算轮廓上所有点角度，并搜索从轮廓起始点顺时针或逆时针出发的最大角度为参考角度angle，然后搜索在参考角度范围内的所有点角度和轮廓坐标点；D4、将参考角度划分为n份，以偏移点为中心，偏移点与轮廓起始点为边，按照平均角度步长进行顺时针或逆时针旋转将各个点角度和轮廓坐标点分别进行分类；D5、计算各个角度内的轮廓坐标点的欧式距离，并对各个角度内的距离求平均值以构成距离链码，即为Davg＝{d0,d1,…,dn‑1}，n为分类数；同理，对各个角度内的点角度求平均值以构成角度链码，即为Aavg＝{a0,a1,…,an‑1}；D6、利用如下公式，对均值距离链码和均值角度链码特征进行归一化处理：E、利用模板随机子轮廓得到的链码特征与第一次匹配后的子轮廓进行第二次匹配，对第一次匹配获得的轮廓点进行再次剔除；F、计算模板随机子轮廓与第二次匹配后的子轮廓的坐标变换矩阵，通过子轮廓坐标变换矩阵将整个模板轮廓投影到目标轮廓上；G、更新子轮廓的坐标变换矩阵：采用均匀采样的方法得到投影模板轮廓上的轮廓采样点，并遍历目标轮廓上的坐标点，求得投影模板轮廓采样点在目标轮廓上的最近邻坐标点，然后利用奇异值分解求得整个模板轮廓的坐标变换矩阵；H、识别最优匹配轮廓：采用均匀采样的方法将模板轮廓点投影到目标轮廓上，并计算这些采样点到目标轮廓与之相对应的最近邻坐标点的欧式距离之和，以其为评判匹配误差的标准，对第二次匹配后的子轮廓进行匹配度评价，选取距离和最小的轮廓即为最优匹配轮廓。