技术领域
[0001] 本
发明涉及
图像处理技术领域，尤其涉及一种基于
线密度的不规则轮廓匹配算法。
背景技术
[0002] 碎纸拼接是图像处理领域中的一个典型分支，它被广泛应用于军事情报获取、司法取证和考古文物修复等领域。对于考古文物的碎片数量较大时，若采用人工拼接的方法不仅会使拼接工作较为繁琐复杂，消耗大量的人
力、物力，还可能对出土文物造成二次损坏。因此采用计算机辅助完成碎片拼接是非常必要的。
[0003] 碎纸拼接技术是通过运用图像处理技术和计算机辅助技术，将原本
破碎的目标经过图像获取、图像预处理、轮廓匹配、图像拼接几个步骤恢复出它的原貌。其中，轮廓匹配技术是匹配拼接技术中最关键的步骤之一，其核心技术是轮廓特征描述和相似性度量。特征描述的方法有很多，总的来说可分为基于空间域和基于变换域两种。目前，基于空间域的研究方法包括链码、多边形逼近、样条等；基于变换域的研究方法包括傅里叶变换描述符、
小波变换描述符等，其基本思想是对物体边缘点作傅里叶变换或者小波变换，将变换系数作为形状的描述符。相似性度量方法主要有欧氏距离、闵可夫斯基距离、汉明距离等。最常用的相似性度量方法是欧式距离，其优点是物理意义清晰，计算简单，适用范围比较广，缺点是与人的视觉
感知有较大差异。现有的轮廓匹配技术匹配准确度和速度较低，且计算量比较大。
发明内容
[0004] 本发明针对现有的轮廓匹配技术匹配准确度和速度较低，且计算量较大的问题，提出基于线密度的不规则轮廓匹配算法。
[0005] 为了实现上述目的，本发明采用以下技术方案：
[0006] 一种基于线密度的不规则轮廓匹配算法，包括：
[0007] 步骤1：获取碎纸的图像信息，从所述碎纸的图像信息中检测出目标碎纸图像及待匹配碎纸图像的边缘轮廓曲线；
[0008] 步骤2：对所述待匹配碎纸图像的边缘轮廓曲线进行匹配预处理；
[0009] 步骤3：利用线密度算法对所述目标碎纸图像的边缘轮廓曲线及匹配预处理后的待匹配碎纸图像的边缘轮廓曲线进行处理，获得目标碎纸图像及待匹配碎纸图像的边缘轮廓曲线的线密度；
[0010] 步骤4：通过目标碎纸图像与待匹配碎纸图像的边缘轮廓曲线的线密度之间的距离衡量所述待匹配碎纸图像的边缘轮廓曲线与所述目标碎纸图像的边缘轮廓曲线之间的匹配程度，距离越大，轮廓匹配程度越小，距离越小，轮廓匹配程度越大。
[0011] 进一步地，所述步骤1包括：
[0012] 步骤1.1：获取碎纸的图像信息，所述碎纸的图像信息包括目标碎纸图像信息及待匹配碎纸图像信息；
[0013] 步骤1.2：对所述目标碎纸图像信息及待匹配碎纸图像信息进行图像去噪处理；
[0014] 步骤1.3：对图像去噪处理后的所述目标碎纸图像信息及待匹配碎纸图像信息进行图像灰度化处理；
[0015] 步骤1.4：对图像灰度化处理后的所述目标碎纸图像信息及待匹配碎纸图像信息进行图像二值化处理，得到目标碎纸及待匹配碎纸的二值化图像；
[0016] 步骤1.5：对所述二值化图像进行二值形态学法轮廓的提取，提取出的目标碎纸图像或待匹配碎纸图像的边缘轮廓曲线为：
[0017] C＝I-IθS′ (1)
[0018] 其中，C为目标碎纸图像或待匹配碎纸图像的边缘轮廓曲线，I为目标碎纸或待匹配碎纸的二值化图像，S′为结构元素， IθS′为S′对I的
腐蚀。
[0019] 进一步地，所述步骤2包括：
[0020] 步骤2.1：保证每一张待匹配碎纸图像与目标碎纸图像
像素尺寸大小相同的情况下，获取待匹配碎纸图像；
[0021] 步骤2.2：通过公式(1)得出各待匹配碎纸图像的边缘轮廓曲线的像素点，通过MATLAB对待匹配碎纸图像进行细化并提取骨架，生成待匹配碎纸图像的边缘轮廓曲线；
[0022] 步骤2.3：比对待匹配碎纸图像的边缘轮廓曲线的像素点数，设置
阈值，以获取的所有碎片像素点尽可能相等为原则，筛选出在阈值范围内出现的待匹配碎纸图像的边缘轮廓曲线，得到符合阈值筛选的待匹配碎纸图像的边缘轮廓曲线数n。
[0023] 进一步地，所述步骤3包括：
[0024] 步骤3.1：将目标碎纸图像的边缘轮廓曲线及匹配预处理后的待匹配碎纸图像的边缘轮廓曲线划分为多个包含X个像素点数的曲线段；
[0025] 步骤3.2：将目标碎纸图像的边缘轮廓曲线及匹配预处理后的待匹配碎纸图像的边缘轮廓曲线进行数学建模，得出所述线密度的求取公式为：
[0026] Linedensityj＝Lj/U (2)
[0027] 其中，Linedensityj为第j个包含X个像素点数的曲线段的线密度，Lj为第j个包含X个像素点数的曲线段的长度，U为单位弧长；
[0028] 步骤3.3：通过公式(2)得出目标碎纸图像的边缘轮廓曲线的线密度W及第i个待匹配碎纸图像的边缘轮廓曲线的线密度Di，其中，W和Di为一维数组，数组的长度等于目标碎纸图像的边缘轮廓曲线或待匹配碎纸图像的边缘轮廓曲线包含X个像素点数的曲线段的个数。
[0029] 进一步地，所述步骤4包括：
[0030] 步骤4.1：通过所述W及所述Di得出目标碎纸图像与第i个待匹配碎纸图像的边缘轮廓曲线的线密度之间的距离Zi：
[0031] Zi＝|W-Di| (3)
[0032] 步骤4.2：通过所述Zi求取目标碎纸图像与待匹配碎纸图像的边缘轮廓曲线的线密度之间的距离的均值M和方差S：
[0033]
[0034]
[0035] 步骤4.3：为所述均值M和方差S分别设定
门限值threshold1、threshold2，若M小于threshold1且S小于threshold2，则第i个待匹配碎纸图像的边缘轮廓曲线与目标碎纸图像的边缘轮廓曲线相匹配，否则匹配失败。
[0036] 与
现有技术相比，本发明具有的有益效果：
[0037] 本发明利用二值形态学法对二维碎纸图像进行
边缘检测，提取出的轮廓曲线较为完整，连续性很好；通过匹配预处理，缩小匹配目标的待匹配范围；运用线密度算法，将碎纸轮廓的匹配与曲线的相似性联系起来，达到了较好的匹配效果；充分利用轮廓的几何特征，大大减少了运行时间，提高了运行结果的准确率。
附图说明
[0038] 图1为本发明
实施例一种基于线密度的不规则轮廓匹配算法的基本
流程图；
[0039] 图2为本发明实施例另一种基于线密度的不规则轮廓匹配算法的基本流程图；
[0040] 图3为本发明实施例一种基于线密度的不规则轮廓匹配算法的轮廓曲线检测流程图；
[0041] 图4为本发明实施例一种基于线密度的不规则轮廓匹配算法的匹配预处理流程图；
[0042] 图5为本发明实施例一种基于线密度的不规则轮廓匹配算法的待匹配碎纸图像的边缘轮廓曲线示意图；
[0043] 图6为本发明实施例一种基于线密度的不规则轮廓匹配算法的线密度得出流程图；
[0044] 图7为本发明实施例一种基于线密度的不规则轮廓匹配算法的构造
控制器初始
位置流程图；
[0045] 图8为本发明实施例一种基于线密度的不规则轮廓匹配算法的缘轮廓曲线匹配流程图。
具体实施方式
[0046] 下面结合附图和具体的实施例对本发明做进一步的解释说明：
[0047] 实施例一：
[0048] 如图1所示，一种基于线密度的不规则轮廓匹配算法，包括以下步骤：
[0049] 步骤S101：获取碎纸的图像信息，从所述碎纸的图像信息中检测出目标碎纸图像及待匹配碎纸图像的边缘轮廓曲线；
[0050] 步骤S102：对所述待匹配碎纸图像的边缘轮廓曲线进行匹配预处理；
[0051] 步骤S103：利用线密度算法对所述目标碎纸图像的边缘轮廓曲线及匹配预处理后的待匹配碎纸图像的边缘轮廓曲线进行处理，获得目标碎纸图像及待匹配碎纸图像的边缘轮廓曲线的线密度；
[0052] 步骤S104：通过目标碎纸图像与待匹配碎纸图像的边缘轮廓曲线的线密度之间的距离衡量所述待匹配碎纸图像的边缘轮廓曲线与所述目标碎纸图像的边缘轮廓曲线之间的匹配程度，距离越大，轮廓匹配程度越小，距离越小，轮廓匹配程度越大。
[0053] 本发明利用二值形态学法对二维碎纸图像进行边缘检测，提取出的轮廓曲线较为完整，连续性很好；通过匹配预处理，缩小匹配目标的待匹配范围；运用线密度算法，将碎纸轮廓的匹配与曲线的相似性联系起来，达到了较好的匹配效果；充分利用轮廓的几何特征，大大减少了运行时间，提高了运行结果的准确率。
[0054] 实施例二：
[0055] 如图2所示，另一种基于线密度的不规则轮廓匹配算法，包括以下步骤：
[0056] 步骤S201：获取碎纸的图像信息，采用二值形态学方法从所述碎纸的图像信息中检测出目标碎纸图像及待匹配碎纸图像的边缘轮廓曲线；
[0057] 具体地，如图3所示，所述步骤S201包括：
[0058] 步骤S2011：获取碎纸的图像信息，所述碎纸的图像信息包括目标碎纸图像信息及待匹配碎纸图像信息；
[0059] 步骤S2012：对所述目标碎纸图像信息及待匹配碎纸图像信息进行图像去噪处理；
[0060] 步骤S2013：对图像去噪处理后的所述目标碎纸图像信息及待匹配碎纸图像信息进行图像灰度化处理；
[0061] 步骤S2014：对图像灰度化处理后的所述目标碎纸图像信息及待匹配碎纸图像信息进行图像二值化处理，得到目标碎纸及待匹配碎纸的二值化图像；
[0062] 步骤S2015：对所述二值化图像进行二值形态学法轮廓的提取，提取出的目标碎纸图像或待匹配碎纸图像的边缘轮廓曲线为：
[0063] C＝I-IθS′ (1)
[0064] 其中，C为目标碎纸图像或待匹配碎纸图像的边缘轮廓曲线，I为目标碎纸或待匹配碎纸的二值化图像，S′为结构元素， IθS′为S′对I的腐蚀。
[0065] 二值形态学法轮廓提取的思想是用具有一定形态的结构元素腐蚀碎纸的二值化图像，并用碎纸的二值化图像减去腐蚀后的图像，得到的结果即为此碎纸图像的边缘轮廓曲线，此方法虽对噪声不是很敏感，但提取出的轮廓曲线较为完整，连续性很好。
[0066] 步骤S202：对所述待匹配碎纸图像的边缘轮廓曲线进行匹配预处理；在大批量进行匹配之前，以获取的所有碎片像素点尽可能相等为原则，通过匹配预处理的方法，在待匹配轮廓中筛选出一部分，进而缩小匹配目标的待匹配范围。
[0067] 具体地，如图4所示，所述步骤S202包括：
[0068] 步骤S2021：保证每一张待匹配碎纸图像与目标碎纸图像像素尺寸大小相同的情况下，获取待匹配碎纸图像；
[0069] 步骤S2022：通过公式(1)得出各待匹配碎纸图像的边缘轮廓曲线的像素点，通过MATLAB对待匹配碎纸图像进行细化并提取骨架，生成待匹配碎纸图像的边缘轮廓曲线；
[0070] 步骤S2023：比对待匹配碎纸图像的边缘轮廓曲线的像素点数，设置阈值，以获取的所有碎片像素点尽可能相等为原则，筛选出在阈值范围内出现的待匹配碎纸图像的边缘轮廓曲线，得到符合阈值筛选的待匹配碎纸图像的边缘轮廓曲线数n。作为一种可实施方式，本实施例中，匹配预处理之前共得出6条待匹配碎纸图像的边缘轮廓曲线，匹配预处理之后，得出4条待匹配碎纸图像的边缘轮廓曲线，如图5所示，匹配预处理之前共得出轮廓曲线1至轮廓曲线6共6条待匹配碎纸图像的边缘轮廓曲线，因为轮廓曲线3及轮廓曲线6与其他轮廓曲线的长度(像素点数)相差较大，故通过匹配预处理对其进行删除，即匹配预处理之后，得出轮廓曲线1、轮廓曲线2、轮廓曲线4及轮廓曲线5共4条待匹配碎纸图像的边缘轮廓曲线，即n为4。
[0071] 步骤S203：利用线密度算法对所述目标碎纸图像的边缘轮廓曲线及匹配预处理后的待匹配碎纸图像的边缘轮廓曲线进行处理，获得目标碎纸图像及待匹配碎纸图像的边缘轮廓曲线的线密度；线密度(Linedensity)是指在提取出的碎纸片边缘曲线上，求出每隔相同的X个像素点数的弧长L所包含的单位弧长(根据实验求得的最佳单位弧长)U的数量。通过计算分析轮廓曲线的线密度，从而获取碎纸边缘用于区分彼此的关键参数。该算法中最核心的步骤是求出轮廓曲线的线密度，在此可将轮廓曲线建模成数学中的函数曲线。
[0072] 具体地，如图6所示，所述步骤S203包括：
[0073] 步骤S2031：将目标碎纸图像的边缘轮廓曲线及匹配预处理后的待匹配碎纸图像的边缘轮廓曲线划分为多个包含X个像素点数的曲线段；
[0074] 步骤S2032：将目标碎纸图像的边缘轮廓曲线及匹配预处理后的待匹配碎纸图像的边缘轮廓曲线进行数学建模，得出所述线密度的求取公式为：
[0075] Linedensity＝L/U (2)
[0076] 其中，Linedensityj为第j个包含X个像素点数的曲线段的线密度，Lj为第j个包含X个像素点数的曲线段的长度，U为单位弧长；
[0077] 步骤S2033：通过公式(2)得出目标碎纸图像的边缘轮廓曲线的线密度W及第i个待匹配碎纸图像的边缘轮廓曲线的线密度Di，其中，W和Di为一维数组，数组的长度等于目标碎纸图像的边缘轮廓曲线或待匹配碎纸图像的边缘轮廓曲线包含X个像素点数的曲线段的个数。具体地，轮廓曲线1、轮廓曲线2、轮廓曲线4及轮廓曲线5共4条待匹配碎纸图像的边缘轮廓曲线的线密度图如图7所示。
[0078] 步骤S204：通过目标碎纸图像与待匹配碎纸图像的边缘轮廓曲线的线密度之间的距离衡量所述待匹配碎纸图像的边缘轮廓曲线与所述目标碎纸图像的边缘轮廓曲线之间的匹配程度，距离越大，轮廓匹配程度越小，距离越小，轮廓匹配程度越大。
[0079] 具体地，如图8所示，所述步骤S204包括：
[0080] 步骤S2041：通过所述W及所述Di得出目标碎纸图像与第i个待匹配碎纸图像的边缘轮廓曲线的线密度之间的距离Zi：
[0081] Zi＝|W-Di| (3)
[0082] 步骤S2042：通过所述Zi求取目标碎纸图像与待匹配碎纸图像的边缘轮廓曲线的线密度之间的距离的均值M和方差S：
[0083]
[0084]
[0085] 步骤S2043：将均值M和方差S称为轮廓匹配度的参数，并用M和S表示两个曲线的匹配度，为所述均值M和方差S分别设定门限值threshold1、threshold2，若M小于threshold1且S小于threshold2，则第i个待匹配碎纸图像的边缘轮廓曲线与目标碎纸图像的边缘轮廓曲线相匹配，否则匹配失败。具体地，threshold1和threshold2的值通过实验得出，从而实现精准匹配，针对不同批次的碎纸的图像，均值M和方差S的门限值thresho1d1和threshold2也将发生改变。作为一种可实施方式，本实施例中，轮廓曲线5为最终的目标碎纸图像的边缘轮廓曲线匹配结果。通过大量实验可以得出，本发明匹配的准确率可达到90％。
[0086] 本发明利用二值形态学法对二维碎纸图像进行边缘检测，提取出的轮廓曲线较为完整，连续性很好；通过匹配预处理，缩小匹配目标的待匹配范围；运用线密度算法，将碎纸轮廓的匹配与曲线的相似性联系起来，达到了较好的匹配效果；充分利用轮廓的几何特征，大大减少了运行时间，提高了运行结果的准确率。
[0087] 以上所示仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
