本发明属于数字图像处理技术领域，特别涉及一种基于轮廓匹配优化的矩形检测提取方法。
背景技术：
大量人工物体形状呈现为矩形特征，例如自然场景图像中的路牌、车牌，高分遥感影像中的汽车、建筑物、盐田。这种矩形特征与自然场景目标中具有分型特征的形状具有显著的区别，为自动检测、提取这些目标提供了极佳的特征。尤其是伴随着我国资源三号、高分二号等一批高空间分辨率遥感卫星发射，获取了大量急需实时处理的高分遥感影像。然而，当前以分割为基础的面向对象方法在提取这些目标时，常常会遇到边缘模糊、断裂、弯曲、内部差异大等状况，导致分割效果较差，进而影响最终的提取结果。因此，有针对性的矩形检测提取方法研究十分必要。
现有的矩形检测方法多是采用自底向上的先分解后组合的策略。这种策略首先从影像中提取边缘或者直线，然后由矩形对边平行、邻边垂直的约束条件，对提取的边缘或直线进行组合、分类，进而实现最终的矩形检测提取。例如目前引用最多的基于hough直线变换的矩形检测方法(jung,claudiorosito,androdrigoschramm.2004.rectangledetectionbasedonawindowedhoughtransform.paperpresentedatthecomputergraphicsandimageprocessing,xviibraziliansymposium.)，在边缘提取后，根据矩形的平行、垂直约束条件，在hough变换空间中寻找满足这些条件的极值点对，实现矩形的检测提取。但由于所采用的方法中包含逐像素滑窗策略，使得这一方法速度较慢。此外，边缘弯曲等状况也会导致直线检测的失败，进而导致矩形的漏检。其它一些依赖于直线检测结果并组合的策略(河谷大和，小坂清人，吴藤弘康.图像处理装置和矩形检测方法:中国，201210439013.x[p].2012-11-06.)，通常也会遇到类似边缘弯曲、断裂的漏检状况。
技术实现要素：
针对上述问题，本发明提出了一种基于轮廓匹配优化的矩形检测提取方法，采用矩形整体检测策略，通过对矩形轮廓与检测边缘匹配度的优化，实现矩形的检测提取。
本发明的具体技术方案是一种基于轮廓匹配优化的矩形检测提取方法，其特征在于，包括以下步骤：
步骤1：检测和提取原始影像边缘轮廓，得到边缘影像；
步骤2：计算边缘影像中每一个非边缘像素距离最近边缘像素的距离d(p)，然后进行归一化距离变换，得到归一化后的距离s(p)，p为非边缘像素；
步骤3：采用规则网格覆盖边缘影像，并对规则网格中的矩形赋初始值；
步骤4：计算每一个规则网格中的矩形及该矩形领域集中的矩形的综合轮廓匹配度，选取轮廓匹配度最大的矩形作为提取的初始矩形；
步骤5：去除重复矩形和非矩形，得到最终的边缘影像的提取矩形。
更进一步地，所述的步骤1中检测和提取原始影像边缘轮廓的方法为采用sobel算子、prewitt算子、canny算子的方法或用图像分割算法计算闭合轮廓的方法。
更进一步地，采用倒角距离变换的方法，计算所述的步骤2中每一个非边缘像素距离最近边缘像素的距离d(p)。
更时一步地，按下式(i)对d(p)进行归一化距离变换，
其中，σd为尺度因子参数，用于降低噪声的影响。
更进一步地，所述的步骤3中赋有初始值的规则网格rij表达为下式(ii)，
rij＝(i·(1-ξ)h0,j·(1-ξ)w0,w0,h0,0)……(ii)
其中，i∈{i|i≥0,i·(1-ξ)h0≤nr}，
j∈{j|j≥0,j·(1-ξ)w0≤nc}，w0,h0为规则网格中赋有初始值的矩形的宽和高，ξ是规则格网覆盖的重叠率，取值范围为[0,1]，nr和nc分别是原始影像的行数和列数。
更进一步地，所述的步骤4中计算每一个规则网格中的矩形及该矩形领域集中的矩形的轮廓匹配度，选取轮廓匹配度最大的矩形作为提取的初始矩形的方法包括以下步骤：
1)给需计算的规则格网中一个矩形设置一个存储空间r＝(x,y,w,h,θ)，设置一个迭代优化算法的尺度参数的存储空间t，给t设置初始尺度1，其中，x,y初始为该矩形的中心点坐标，w,h初始为该矩形的宽和高，θ初始为该矩形的夹角，
2)按照下式(iii)计算矩形r＝(x,y,w,h,θ)的邻域集n，
n＝{nr|nr＝r+k1a1+k2a2+k3a3+k4a4+k5a5,k1,k2,k3,k4,k5∈{-1,0,1}}……(iii)
其中，a1＝0.5tw0(cosθ,-sinθ,0,0,0)，a2＝0.5tw0(sinθ,cosθ,0,0,0)，a3＝tw0(0,0,1,0,0)，a4＝tw0(0,0,0,1,0)，a5＝45t(0,0,0,0,1)，
3)将矩形r按下式(iv)栅格化，得到栅格化的矩形表达式er，
其中，
c＝(x,y)，m为整数集，
4)按下式(v)计算轮廓匹配度外力fe，
5)按下式(vi)计算轮廓匹配度内力fi，
其中，
6)按下式(vii)计算综合轮廓匹配度f，
7)如果t的数值小于预设的精度，则跳转至步骤8)，否则，提取t的数值除以2后，再次给t赋值，按照步骤2)更新邻域集，然后任意取出领域集中的一个矩形赋值给存储空间r＝(x,y,w,h,θ)，按步骤3)－6)计算综合轮廓匹配度，如果当前计算的矩形综合轮廓匹配度大于前一个矩形的综合轮廓匹配度，则认为当前矩形为最匹配矩形，否则，前一个矩形为最匹配矩形，重复步骤7)，
8)终止迭代，以最匹配矩形为提取的初始矩形。
更进一步地，所述步骤5中去除重复矩形和非矩形的方法是，在所有规则格网中的矩形对应的初始矩形中，选取任意一对初始矩形，如果
更进一步地，所述的第一阈值为3，第二阈值为0.7。
本发明的有益效果是：本发明的方法通过归一化距离变化，使得轮廓匹配度外力的计算能够尽可能多的反应靠近轮廓的边缘像素，并且为后续的重复矩形移除和非矩形移除提供了关键指标；通过轮廓匹配度的外力引入，使得矩形的检测提取能够提取符合预定尺度的矩形，防止矩形的过度收缩或者膨胀；多尺度的邻域集的迭代优化，以较快的速度达到预设的精度；规则格网和近似矩形的初始化策略，配合多尺度邻域迭代优化方法，实现了多个矩形目标的检测提取。
本发明对矩形目标的边缘模糊、断裂、弯曲、局部缺失等状况，具有较好的鲁棒性，并具有多个矩形目标的检测提取能力。
附图说明
图1为本发明的基于轮廓匹配优化的矩形检测提取方法流程图；
图2为实施例采用的原始影像；
图3为影像边缘提取结果图；
图4为归一化距离变换结果图；
图5是规则格网生成结果图；
图6是匹配度优化迭代结果图；
图7是重叠矩形和非矩形去除结果图。
具体实施方式
下面结合说明书附图对本发明的具体技术方案作进一步地描述。
如附图1所示，本实施例中的矩形目标提取方法主要分成以下五个主要部分：
如附图2所示，在本实施例中，所采用的原始影像为高空间分辨率遥感影像，分辨率达1米，图像尺寸为1772行×2342列，待检测提取的目标是盐田。
步骤1：检测和提取原始影像边缘轮廓，得到边缘影像。常规的sobel算子、prewitt算子、canny算子，以及图像分割算法的闭合轮廓，均可应用于边缘提取。由于现有软件或者计算机视觉开源库中均提供了大量的边缘检测提取算法，例如matlab以及opencv，此处不再赘述。本实施例中，使用canny算子进行提取，结果如附图3所示；
步骤2：计算边缘影像中每一个非边缘像素距离最近边缘像素的距离d(p)，然后进行归一化距离变换，得到归一化后的距离s(p)，p为非边缘像素。通过对边缘影像进行倒角距离变换，计算边缘影像中每一个非边缘像素距离其最近边缘像素的距离，如附图4所示。鉴于倒角距离变换技术比较成熟，现有软件或者计算机视觉开源库中均提供了此算法，例如matlab以及opencv，此处不再赘述。
按下式(i)对d(p)进行归一化距离变换，
其中，σd为尺度因子参数，用于降低噪声的影响，σd根据处理的原始影像的尺寸、待提取的目标大小而设定，本例中设置为10；
步骤3：采用规则网格覆盖边缘影像，并对规则网格中的矩形赋初始值。
赋有初始值的规则网格rij表达为下式(ii)，
rij＝(i·(1-ξ)h0,j·(1-ξ)w0,w0,h0,0)……(ii)
其中，i∈{i|i≥0,i·(1-ξ)h0≤nr}，
j∈{j|j≥0,j·(1-ξ)w0≤nc}，w0,h0为规则网格中赋有初始值的矩形的宽和高，ξ是规则格网覆盖的重叠率，取值范围为[0,1]，本例设置为0.55，nr和nc分别是原始影像的行数和列数，即1772和2342。生成的网格如附图5所示；
步骤4：计算每一个规则网格中的矩形及该矩形领域集中的矩形的综合轮廓匹配度，选取轮廓匹配度最大的矩形作为提取的初始矩形。具体方法包括以下步骤：
1)给需计算的规则格网中一个矩形设置一个存储空间r＝(x,y,w,h,θ)，设置一个迭代优化算法的尺度参数的存储空间t，采用的迭代优化算法为多尺度爬山算法，给t设置初始尺度1，其中，x,y初始为该矩形的中心点坐标，w,h初始为该矩形的宽和高，θ初始为该矩形的夹角，
2)按照下式(iii)计算矩形r＝(x,y,w,h,θ)的邻域集n，
n＝{nr|nr＝r+k1a1+k2a2+k3a3+k4a4+k5a5,k1,k2,k3,k4,k5∈{-1,0,1}}……(iii)
其中，a1＝0.5tw0(cosθ,-sinθ,0,0,0)，a2＝0.5tw0(sinθ,cosθ,0,0,0)，a3＝tw0(0,0,1,0,0)，a4＝tw0(0,0,0,1,0)，a5＝45t(0,0,0,0,1)，
3)将矩形r按下式(iv)栅格化，得到栅格化的矩形表达式er，
其中，
c＝(x,y)，m为整数集，
4)按下式(v)计算轮廓匹配度外力fe，
5)按下式(vi)计算轮廓匹配度内力fi，
其中，
6)按下式(vii)计算综合轮廓匹配度f，
7)如果t的数值小于预设的精度，本实施例的预设的精度为0.005，则跳转至步骤8)，否则，提取t的数值除以2后，再次给t赋值，按照步骤2)更新邻域集，然后任意取出领域集中的一个矩形赋值给存储空间r＝(x,y,w,h,θ)，按步骤3)－6)计算综合轮廓匹配度，如果当前计算的矩形综合轮廓匹配度大于前一个矩形的综合轮廓匹配度，则认为当前矩形为最匹配矩形，否则，前一个矩形为最匹配矩形，重复步骤7)，
8)终止迭代，以最匹配矩形为提取的初始矩形。
通过以上方法，可以求出规则网格中每一个矩形对应的初始矩形，如附图6所示。
步骤5：去除重复矩形和非矩形，得到最终的边缘影像的提取矩形。具体方法是，在所有规则格网中的矩形对应的初始矩形中，选取任意一对初始矩形，如果
则视为这两个矩形为重复矩形，保留其中fe较大的矩形，移除较小的矩形，将所有初始矩形两两比较完毕后，再将fe小于设定的第二阈值0.7的矩形视为非矩形予以移除，最终剩余的初始矩形为最终的边缘影像的提取矩形，如附图7所示。这个阈值是根据前面的参数σd的设置情况而设定的，通常取
本发明的方法利用轮廓匹配度的优化，实现了具有矩形形状的盐田高精度提取，所得的结果是严格的矩形形状，能够满足gis数据库直接入库的要求。
显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。