专利名称：改进的Zernike矩边缘检测方法
技术领域：
本发明涉及一种图像处理方法，尤其是涉及一种针对交点亚像素定位的改进的Zernike矩边缘检测方法。
背景技术：
工业CT(Industry Computerized Tomogragh)广泛应用于航空、航天、汽车等等工业中精密工件内部结构的测量及其缺陷检测领域。工业CT图像尺寸测量方法是一种以工业CT图像为信息载体并从中提取定量数据的测量方法，它以工业CT断层图像为研究对象，计算目标的长度、宽度、高度、壁厚等各种几何参数。与其它测量方法相比，其既能测得实体复杂的内部几何形状(尤其是封闭内腔)的参数，又不破坏实体。边缘提取是工业CT图像 测量的基本任务，它对测量结果的精度有很大的影响。传统的边缘检测算子如Sobel，Canny等只能提取像素级边缘不能满足定位的高精度需求。Hueckel首先提出亚像素边缘检测算法，到目前为止，亚像素级边缘检测方法主要分为四大类插值法、拟合法、基于小波变换的方法、基于矩的方法、基于活动轮廓模型的方法。其中基于矩的方法具有快速、高精度且有一定抗噪能力的优势。Ghosal和Mehrotal首次提出了利用Zernike正交矩来检测亚像素边缘(参见Ghosal S，Mehrotra R. Orthogonal moment operators for sub-pixel edge detection. PatternRecognition, 1993, 26(2) : 295-306),在他们的算法中建立了理想的阶跃灰度模型，通过图像的三个不同阶次Zernike矩计算模型的4个参数，以这4个参数为判断边缘的依据并确定图像中物体边缘的位置。分析发现，Ghosal提出的Zernike矩方法是基于直线边缘的，当检测目标为交点附近像素点时，算法将每个模板内的折线边缘作为直线处理，亚像素检测效果较差。因此,需要一种改进的Zernike矩边缘检测算法,解决Ghosal算法在利用Zernike矩进行亚像素边缘检测时交点附近的边缘检测误差大的问题。
发明内容
有鉴于此，本发明所要解决的技术问题是提供一种改进的Zernike矩边缘检测方法，对Ghosal算法在利用Zernike矩进行亚像素边缘检测时交点附近的亚像素边缘误差大的不足进行判断和修正。本发明的目的是这样实现的本发明提供的改进的Zernike矩边缘检测方法，包括以下步骤I)获取图像，提取图像的边缘；2)对边缘进行细化处理；3)获取边缘方向参数分布和边缘方向参数差分值a (J)；4)根据预设阈值T来判断边缘像素点是否为交点附近像素点，如果是，则该像素点为交点附近像素点；
5)如果否，则该像素点为非交点附近像素点，判断为待修正像素点，并利用非交点附近像素点进行最小二乘拟合，得到直线方程P ；6)利用直线方程P来计算待修正像素点的坐标。进一步，所述步骤2)中边缘点细化是对图像中的所有边缘点做边缘梯度方向上边缘距离参数I的非极小值抑制来细化边缘。进一步，所述步骤3)中对所有边缘点采用8链码跟踪获取边缘方向参数Φ分布。进一步，所述步骤3)中采用以下公式来计算检测到边缘点的边缘方向参数Φ的向前差分值Δ φ
Δ φ = φ (η+1) - Φ (η), η = I, 2, 3. . . (NUM-1),其中，Δ φ表示边缘方向参数向前差分值，Φ (η+1)表示后一个边缘点的边缘方向参数，Φ (η)表示当前边缘点的边缘方向参数，NUM为检测到的像素点总数。进一步，所述步骤4)中判断像素点是否为交点附近像素点包括以下步骤如果边缘方向参数差分值满足关系式I Λ φ I彡T时，则边缘像素点为交点附近待修正像素点，如果边缘方向参数差分值满足关系式I Δ φ I < Τ,则边缘像素点为非交点附近像素点。进一步，所述步骤5)中最小二乘拟合直线方程根据以下公式进行
二 I其中，F (X，P) =Ax+By+C = 0表示拟合的直线方程，P表示直线参数P= [A, B, C],X = [X，y, 1]τ表示坐标向量,对于直线上任意边缘点(U,.)，则F(X，P)表示该点到F(X，P)=O的代数距离,根据目标函数求解直线参数P。进一步,所述步骤6)中交点附近的像素点的精确坐标是通过直线参数P来确定的。本发明的优点在于本发明与传统的Zernike矩边缘检测方法相比，针对交点附近像素点亚像素定位误差较大进行改进，采用Zernike矩方法对工业CT图像进行亚像素边缘提取，并细化边缘，然后采用8链码跟踪方法获取边缘方向参数Φ分布，并计算边缘方向参数差分值Λ φ，再选取阈值T，判断边缘点是否为交点附近的像素点，若为非交点附近的边缘点，进行最小二乘拟合从而得到直线方程，再通过直线方程计算交点附近待修正边缘点的亚像素坐标。改进Zernike矩边缘检测方法很好地弥补了 Ghosal算法在交点附近像素点定位误差较大的不足，对高精度面积测量和高精度三维测量有重要意义。本发明的其它优点、目标和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。本发明的目标和其它优点可以通过下面的说明书，权利要求书，以及附图中所特别指出的结构来实现和获得。
为了使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作进一步的详细描述，其中
图I为本发明实施例提供的改进的Zernike矩边缘检测方法流程图；图2为本发明实施例的矩形数字模拟二值图；图3为本发明实施例的边缘方向参数$分布图；图4为本发明实施例的边缘方向参数差分值A (]5分布图；图5为本发明实施例的亚像素边缘定位结果图。
具体实施例方式以下将结合附图，对本发明的优选实施例进行详细的描述；应当理解，优选实施例仅为了说明本发明，而不是为了限制本发明的保护范围。图I为本发明实施例提供的改进的Zernike矩边缘检测方法流程图，实施例I本发明提供的改进的Zernike矩边缘检测方法，包括以下步骤I)获取图像，提取图像的边缘；本实施例中采用Zernike矩方法对图像进行边缘提取，本实施例中的图像为工业CT图像。2)对边缘进行细化处理；本实施例的细化边缘是通过对图像中的每一个边缘点在梯度最大方向上做边缘距离参数I的非极小值抑制来进行的。3)获取边缘方向参数0分布和边缘方向参数差分值A (]5 ;本实施例采取8链码跟踪获取边缘方向参数0分布；本实施例采用以下公式来计算检测到边缘点的边缘方向参数0的向前差分值A (J)A (J) = (J) (n+1) - (n), n = I, 2, 3. . . (NUM-I),其中，A (J)表示边缘方向参数向前差分值，(n+1)表示后一个边缘点的边缘方向参数，Mn)表示当前边缘点的边缘方向参数，NUM为检测到的像素点总数。4)根据预设阈值T来判断边缘像素点是否为交点附近像素点；判断像素点是否为交点附近像素点如下步骤如果边缘方向参数差分值满足关系式I A 4) I彡T时,则边缘像素点为交点像素点，如果边缘方向参数差分值满足关系式I A (J) I < T,则边缘像素点为待修正像素点,如果是，则该像素点为交点附近像素点；5 )如果否,则该像素点为非交点附近像素点,判断为待修正像素点，并利用非交点附近像素点进行最小二乘拟合，得到直线方程P ;最小二乘拟合直线方程根据以下公式进行
权利要求
1.改进的Zernike矩边缘检测方法，其特点在于，包括以下步骤 1)获取图像并提取图像的边缘； 2)对边缘进行细化处理； 3)获取边缘方向参数分布和边缘方向参数差分值A(J)； 4)根据预设阈值T来判断边缘像素点是否为交点附近像素点，如果是，则该像素点为交点附近像素点； 5)如果否，则该像素点为非交点附近像素点，判断为待修正像素点，并利用非交点附近像素点进行最小二乘拟合，得到直线方程P ; 6)利用直线方程P来计算待修正像素点的坐标。
2.根据权利要求I所述改进的Zernike矩边缘检测方法，其特点在于，所述步骤2)中边缘点细化是对图像中的所有边缘点做边缘方向上边缘距离参数I的非极小值抑制来细化边缘。
3.根据权利要求I所述改进的Zernike矩边缘检测方法，其特点在于，所述步骤3)中对所有边缘点采用8链码跟踪获取边缘方向参数0分布。
4.根据权利要求I所述改进的Zernike矩边缘检测方法，其特点在于，所述步骤3)中采用以下公式来计算检测到边缘点的边缘方向参数0的向前差分值△小 A (J) = (J) (n+1) ~<i> (n), n = I, 2, 3. . . (NUM-1), 其中，△ 4)表示边缘方向参数向前差分值，(n+1)表示后一个边缘点的边缘方向参数，(n)表示当前边缘点的边缘方向参数，NUM为检测到的像素点总数。
5.根据权利要求I所述改进的Zernike矩边缘检测方法，其特点在于，所述步骤4)中判断像素点是否为交点附近像素点包括以下步骤 如果边缘方向参数差分值满足关系式I A Cj5 I彡T时，则边缘像素点为交点附近待修正像素点，如果边缘方向参数差分值满足关系式I A Cj5 I < T，则边缘像素点为非交点附近像素点。
6.根据权利要求I所述改进的Zernike矩边缘检测方法，其特点在于，所述步骤5)中最小二乘拟合直线方程根据以下公式进行
7.根据权利要求I所述改进的Zernike矩边缘检测方法，其特点在于，所述步骤6)中交点附近的像素点的精确坐标是通过直线参数P来确定的。
全文摘要
改进的Zernike矩边缘检测方法，涉及一种图像处理方法。提供一种针对交点附近像素点亚像素定位的改进的Zernike矩边缘检测方法。采用Zernike矩方法对工业CT图像进行亚像素边缘提取，并细化边缘，再采用8链码跟踪方法，获取边缘方向参数φ分布，计算边缘方向参数差分值Δφ，选取阈值T，判断边缘点是否为交点附近的像素点，若不是交点附近边缘点，则进行最小二乘拟合，得到直线方程，再通过直线方程计算交点附近待修正边缘点的亚像素坐标。改进Zernike矩边缘检测方法很好地弥补了Ghosal算法在交点附近像素点定位误差较大的不足，对高精度面积测量和高精度三维测量有重要意义。
文档编号G06T7/00GK102637300SQ20121013624
公开日2012年8月15日 申请日期2012年4月26日 优先权日2012年4月26日
发明者王慧倩, 王珏, 罗珊, 蔡玉芳 申请人:重庆大学