发明名称 一种基于高斯金字塔分解的不规则圆检测方法
摘要 本发明涉及一种基于高斯金字塔分解的不规则圆检测方法。现有的圆检测算法在高分辨率、多圆的检测环境中会发生误检或漏检。本发明方法针对数字图像中的圆目标,通过高斯图像金字塔分解方法将图像分解成不同分辨率的图像,针对不同等级分辨率,采用不同圆检测方法,在低分辨率图像上进行粗略检测,之后在高分辨率图像进行精细圆检测,且在检测的各步骤中,相应的阈值给出了明确的设定方法。本发明提高了对不规则圆的检测精度,能够分析高噪声、多圆的复杂图像。
申请公布号 CN102129685A 申请公布日期 2011.07.20
申请号 CN201110072275.2 申请日期 2011.03.24
申请人 杭州电子科技大学 发明人 范影乐;赵磊;王海玲
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 杜军
主权项 1.一种基于高斯金字塔分解的不规则圆检测方法,其特征在于该方法包括以下步骤:步骤(1)用高斯图像金字塔将原始图像分解为低分辨率图像,所述的低分辨率图像分解级数为<img file="2011100722752100001DEST_PATH_IMAGE002.GIF" wi="14" he="21" />;在分析一组原始图像时,首先选择该组中圆半径差距最大的图像,该组图像的分解级数<img file="664870DEST_PATH_IMAGE002.GIF" wi="14" he="21" />具体确定方法分为以下两种情况:①当原始图像最大圆半径小于等于最小圆半径10倍时,分别以0、1、2、3、4、5为分解级数分析所选图像,然后根据各个图像分解级数下的分析速度,选择分析速度最快的分解级数为该组图像的分解级数<img file="881088DEST_PATH_IMAGE002.GIF" wi="14" he="21" />;②当原始图像最大圆半径大于最小圆半径10倍时,首先设原始图像中最大的圆半径为<img file="2011100722752100001DEST_PATH_IMAGE004.GIF" wi="32" he="25" />,将待检测目标圆的半径范围划分成<img file="2011100722752100001DEST_PATH_IMAGE006.GIF" wi="22" he="18" />个半径范围区域<img file="2011100722752100001DEST_PATH_IMAGE008.GIF" wi="21" he="25" />,<img file="2011100722752100001DEST_PATH_IMAGE010.GIF" wi="100" he="22" />,其中<img file="2011100722752100001DEST_PATH_IMAGE012.GIF" wi="98" he="26" />,<img file="950544DEST_PATH_IMAGE008.GIF" wi="21" he="25" />的范围为<img file="2011100722752100001DEST_PATH_IMAGE014.GIF" wi="165" he="26" />,此时可通过情况①中的方法确定每个半径范围区域<img file="739246DEST_PATH_IMAGE008.GIF" wi="21" he="25" />的最佳分解级数,然后对原始图像从最大的半径区域范围开始进行<img file="416215DEST_PATH_IMAGE006.GIF" wi="22" he="18" />次多分辨率的圆检测,每次检测只检测半径范围<img file="854150DEST_PATH_IMAGE008.GIF" wi="21" he="25" />的圆,检测到的圆在下次检测时应从分解后的图像上擦除;步骤(2)、运用改进的hough圆检测算法对低分辨率图像进行圆检测,检测出低分辨率图像中圆的个数以及各个圆的在低分辨率图像上的圆心和半径,具体是:1) 对图像的边缘点集<i>E</i>中的点跟据其连通性进行分类,分别得到拥有<img file="2011100722752100001DEST_PATH_IMAGE016.GIF" wi="14" he="16" />条连续边缘的边缘点集<i>E</i>的<img file="415712DEST_PATH_IMAGE016.GIF" wi="14" he="16" />个边缘点子集<img file="2011100722752100001DEST_PATH_IMAGE018.GIF" wi="20" he="25" />,<img file="2011100722752100001DEST_PATH_IMAGE020.GIF" wi="88" he="22" />,其中边缘点子集<img file="2011100722752100001DEST_PATH_IMAGE022.GIF" wi="20" he="25" />、<img file="2011100722752100001DEST_PATH_IMAGE024.GIF" wi="21" he="25" />……<img file="2011100722752100001DEST_PATH_IMAGE026.GIF" wi="21" he="25" />按子集中点的个数降序排列;2)对上述各子集按照<img file="928428DEST_PATH_IMAGE022.GIF" wi="20" he="25" />、<img file="41877DEST_PATH_IMAGE024.GIF" wi="21" he="25" />……<img file="904791DEST_PATH_IMAGE026.GIF" wi="21" he="25" />的排列顺序分别运用hough圆检测算法进行检测,具体是从边缘点集<img file="129099DEST_PATH_IMAGE018.GIF" wi="20" he="25" />中,等间隔取三个点,起始位置及间隔距离由边缘点集<img file="393858DEST_PATH_IMAGE018.GIF" wi="20" he="25" />中点的个数决定;设置阈值<img file="2011100722752100001DEST_PATH_IMAGE028.GIF" wi="22" he="25" />,当采样次数大于<img file="114427DEST_PATH_IMAGE028.GIF" wi="22" he="25" />时,则放弃该边缘所有点的采样权;如果检测到存在的圆,则记录圆的圆心和半径,同时取消该圆上所有边缘点的采样权;步骤(3)将低分辨率图像中检测到的每个圆根据其圆心和半径逐一映射到原始图像与之对应的圆存在区域;步骤(4)在上述每个圆存在区域利用最小二乘法圆参数拟合算法逐一检测原始图像中圆的圆心和半径;所述的最小二乘法圆参数拟合算法步骤为:Ⅰ)利用迭代剔除的最小二乘法圆参数拟合算法得到边缘点集<img file="2011100722752100001DEST_PATH_IMAGE030.GIF" wi="24" he="18" />,其步骤为:A)将原始图像中与经由改进的hough变换在低分辨率图像中检测到的圆A的边缘所对应的区域记为P,设<img file="2011100722752100001DEST_PATH_IMAGE032.GIF" wi="21" he="18" />为P中的边缘点集,<img file="199058DEST_PATH_IMAGE032.GIF" wi="21" he="18" />中点的个数为<img file="2011100722752100001DEST_PATH_IMAGE034.GIF" wi="18" he="18" />,将该点集中点的坐标参数拟合到最小二乘法的圆心、半径计算公式中得到圆B;B)设置阈值<img file="2011100722752100001DEST_PATH_IMAGE036.GIF" wi="18" he="25" />,遍历<img file="102423DEST_PATH_IMAGE032.GIF" wi="21" he="18" />中的点,记点到圆心的距离与半径的平方差为<img file="2011100722752100001DEST_PATH_IMAGE038.GIF" wi="21" he="25" />,<img file="2011100722752100001DEST_PATH_IMAGE040.GIF" wi="110" he="22" />当<img file="2011100722752100001DEST_PATH_IMAGE042.GIF" wi="56" he="28" />时,从<img file="657907DEST_PATH_IMAGE032.GIF" wi="21" he="18" />中剔除该点;C)重复执行步骤(A)和(B);设经过剔除后,<img file="50842DEST_PATH_IMAGE032.GIF" wi="21" he="18" />中剩余点数为<img file="2011100722752100001DEST_PATH_IMAGE044.GIF" wi="20" he="20" />,设置阈值<img file="2011100722752100001DEST_PATH_IMAGE046.GIF" wi="22" he="25" />、<img file="2011100722752100001DEST_PATH_IMAGE048.GIF" wi="22" he="25" />,记<img file="324566DEST_PATH_IMAGE032.GIF" wi="21" he="18" />中的剩余点到圆心的距离与半径的平方差为<img file="2011100722752100001DEST_PATH_IMAGE050.GIF" wi="17" he="25" />,<img file="2011100722752100001DEST_PATH_IMAGE052.GIF" wi="104" he="22" />记<img file="2011100722752100001DEST_PATH_IMAGE054.GIF" wi="76" he="46" />,当<img file="2011100722752100001DEST_PATH_IMAGE056.GIF" wi="20" he="18" />小于<img file="703726DEST_PATH_IMAGE046.GIF" wi="22" he="25" />,或剩余边缘点个数小于<img file="739815DEST_PATH_IMAGE048.GIF" wi="22" he="25" />时,停止迭代,输出剩余的边缘点集<img file="365969DEST_PATH_IMAGE030.GIF" wi="24" he="18" />;Ⅱ)将边缘点集<img file="929586DEST_PATH_IMAGE030.GIF" wi="24" he="18" />中的点拟合到最小二乘法圆参数拟合算法中求得原始图像中圆的圆心和坐标。
地址 310018 浙江省杭州市下沙高教园区2号大街