买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：浙江大学

摘要：本发明公开了一种基于垂直纵横线轮廓特征的复杂结构铸件潜在缺陷检测方法。采集复杂结构铸件的X射线图像，采用反锐化掩模算法对图像进行锐化，接着进行灰度级形态学顶帽运算；采用Canny边缘检测，标记连通域，再通过连通域像素面积初步筛选获得缺陷可能区域；提取感兴趣区域，求取八个方向上过感兴趣区域中心的纵横线灰度信息，筛选出两条相互垂直的纵横线；分析两条纵横线灰度图像轮廓变化特征并消除背景影响，根据轮廓变化特征判定是否为真实缺陷。本发明的突出了人眼难以分辨的模糊缺陷区域，改善了缺陷边缘检测效果，减小铸件角落区域结构的影响，有效提高最终缺陷检测率。

主权项：1.一种基于垂直纵横线轮廓特征的复杂结构铸件潜在缺陷检测方法，其特征在于包括以下步骤：1采集复杂结构铸件的X射线图像，采用反锐化掩模算法对图像进行锐化，接着进行灰度级形态学顶帽运算，实现图像预处理；2采用Canny边缘检测，标记连通域，再通过连通域像素面积初步筛选获得缺陷可能区域；3提取感兴趣区域，求取八个方向上过感兴趣区域中心的纵横线灰度信息，筛选出两条相互垂直的纵横线；所述的步骤3具体为：3.1针对每个缺陷可能区域在原X射线图像上截取建立一个矩形的感兴趣区域，感兴趣区域的中心为缺陷可能区域的几何中心，感兴趣区域的长为感兴趣区域沿x轴方向最大尺寸的两倍，感兴趣区域的宽设为感兴趣区域沿y轴方向最大尺寸的两倍，并将感兴趣区域的尺寸调整为32×32像素；3.2经过感兴趣区域的中心建立八个方向，八个方向分别为与图像x轴或者y轴的正方向呈22.5°×k角度的八个直线所在方向，k取0、1、2、3、4、5、6、7；3.3提取感兴趣区域在各方向直线上各个像素点的灰度信息，从八个方向中选取两端点灰度相差最小的方向及与其垂直的方向，分别作为第一纵横线和第二纵横线，将沿第一纵横线和第二纵横线上感兴趣区域各个像素点的灰度信息构建获得第一灰度向量P1和第二灰度向量P2；4分析两条纵横线灰度图像轮廓变化特征并消除背景影响，根据轮廓变化特征判定是否为真实缺陷；所述的步骤4具体为：4.1采用以下方式针对第一灰度向量P1处理获得第一灰度值：绘制第一纵横线的灰度-坐标图，获得第一纵横线所在方向灰度的变化斜率，表达式如下： 式中，P132和P11分别为第一纵横线的灰度-坐标图中坐标最大和最小点对应的灰度值，k0表示第一纵横线的灰度-坐标图中坐标最大和最小点连线的斜率；然后采用以下公式针对感兴趣区域中第一纵横线上的像素点进行计算获得各点的第一灰度值：Q1j＝P1j-P11-k0*j-1式中，Q1j为第一纵横线的灰度-坐标图中第j个点消除背景灰度影响后的灰度值，作为第一灰度值，P1j表示第一纵横线的灰度-坐标图中第j个点的灰度值，j表示灰度-坐标图中点的序号，j取1至32的整数；4.2采用以下方式针对第一灰度向量P1和第二灰度向量P2处理获得第二灰度值：取第一灰度向量P1中的第5至28个元素组成第一背景向量PP1，取第二灰度向量P2中的第5至28个元素组成第二背景向量PP2，绘制第一背景向量PP1和第二背景向量PP2的灰度-坐标图，对第一背景向量PP1和第二背景向量PP2采用以下公式计算处理，求得各自的斜率k1和k2，表达式如下： 式中，PP124、PP11和PP224、PP21分别为PP1和PP2灰度-坐标图中坐标最大和最小点对应的灰度值，斜率k1和k2分别表示第一背景向量PP1和第二背景向量PP2的两端点之间的斜率；然后采用以下公式消除第一背景向量PP1和第二背景向量PP2背景影响，获得各点的第二灰度值和第三灰度值，表达式如下：Q2i＝PP1i-PP11-k1*i-1Q3i＝PP2i-PP21-k2*i-1式中，PP1i表示第一背景向量PP1的灰度-坐标图中第i个坐标点的灰度值，Q2i为第一背景向量PP1的灰度-坐标图中第i个坐标点消除背景灰度影响后的灰度值，作为第二灰度值；PP2i表示第二背景向量PP2的灰度-坐标图中第i个坐标点的灰度值，Q3i为第二背景向量PP2的灰度-坐标图中第i个坐标点消除背景灰度影响后的灰度值，作为第三灰度值，i取1至24的整数；4.3提取轮廓特征计算获得两条纵横线的轮廓特征，轮廓特征包括将所有第一灰度值Q1j经快速傅里叶变换FFT得到的一次谐波幅值F0、将所有第二灰度值Q2i经快速傅里叶变换FFT得到的一次谐波幅值F1和将所有第三灰度值Q3i经快速傅里叶变换FFT得到的一次谐波幅值F2；然后利用一次谐波幅值F1和一次谐波幅值F2采用以下公式计算求得所有第二灰度值Q2i和所有第三灰度值Q3i的一次谐波幅值之间的相似度数S0，表达式如下： 式中，F1和F2分别为所有第二灰度值Q2i和所有第三灰度值Q3i的快速傅里叶变换一次谐波幅值；4.4然后根据幅值F0和相似度数S0各自的上下限预设阈值判断获得感兴趣区域是否为真实缺陷的结果：若幅值F0和相似度数S0均在各自的上限预设阈值和下限预设阈值之间，则感兴趣区域为真实缺陷；否则感兴趣区域不为真实缺陷。

全文数据：基于垂直纵横线轮廓特征的复杂结构铸件潜在缺陷检测方法技术领域[0001]本发明涉及数字图像缺陷检测领域，具体涉及一种基于垂直纵横线轮廓特征的复杂结构铸件潜在缺陷检测方法。背景技术[0002]X射线检测是实现工业铸件无损检测的一种方式，通过X射线获得工业铸件的内部结构图像，再由人工检测确定不合格铸件。这种人工筛选的方式能够适应复杂结构的缺陷检测，但效率低并且对一些模糊的缺陷容易误判。实现工业零件缺陷的自动化检测，对于提高工业零件的生产效率、降低制造成本有着重要意义。[0003]图像处理的方法用于检测铸件缺陷可以大大提高生产效率、降低生产成本，同时满足自动化生产的要求。目前利用图像处理检测铸件内部缺陷的方法可以分为两类:一类为借助无缺陷的模板图像与有缺陷的实际图像差异检测缺陷的方法，这类方法在能够获得优秀模板图像的前提下有精确的检测结果，但要获得优秀模板十分困难;另一类为无模板图像的检测方法，这类方法大多采用目标分割，特征提取与匹配的模式，如通过阈值分割及形状、位置、宽度、对比度等特征进行缺陷识别，但往往对于铸件结构复杂以及亮度不均匀的铸件图像检测效果较差;缺陷追踪匹配法，缺陷可能区域需要在系列每张图中一一匹配，能有效解决亮度不均匀等随机干扰的影响，但计算复杂，且不能很好地解决铸件结构复杂的影响。发明内容[0004]本发明的目的是克服现有技术的不足，提出一种基于垂直纵横线轮廓特征的复杂结构铸件潜在缺陷检测方法，通过灰度级形态学顶帽和Canny边缘检测，提取缺陷可能区域;分析两条相互垂直的纵横线灰度曲线轮廓性质，提取特征。该方法受图像亮度不均匀和铸件复杂结构影响小，可靠性高，可对特征不明显的潜在缺陷进行精确检测。[0005]为实现上述目的，以下为该发明的技术方案：[0006]本发明包括以下步骤：[0007]1采集复杂结构铸件的X射线图像，采用反锐化掩模算法对图像进行锐化，接着进行灰度级形态学顶帽运算突出局部高亮度区域，实现图像预处理；[0008]2采用Canny边缘检测，标记连通域，再通过连通域像素面积初步筛选获得缺陷可能区域；[0009]3提取感兴趣区域，求取八个方向上过感兴趣区域中心的纵横线灰度信息，筛选出两条相互垂直的纵横线；[0010]4分析两条纵横线灰度图像轮廓变化特征并消除背景影响，根据轮廓变化特征判定是否为真实缺陷，具体是将轮廓变化相似且幅度在一定范围内的区域判定为真实缺陷。[0011]所述的复杂结构铸件是指汽车轮毂。[0012]所述的步骤1具体为：[0013]I.1提取图像高频成分，利用高斯算子对原图像高斯滤波，用原图像减去高斯滤波的结果得到图像高频成分，再在原图像上加上图像高频成分的倍数，实现图像锐化，得到锐化图像；[00M]1.2对锐化图像进行形态学顶帽处理:根据预设的潜在缺陷大小，设置与潜在缺陷大小相同的圆形结构体，用圆形结构体对锐化灰度图像进行开运算（先进行腐蚀，再膨胀），再用锐化图像减去开运算之后的图像，获得顶帽运算结果图像。[0015]所述的潜在缺陷的大小应大于最大缺陷，小于X射线图像的长宽尺寸。[0016]所述的步骤2具体为：[0017]2.1边缘检测采用Canny边缘检测方法，先用带有标准偏差的高斯滤波器平滑图像处理，然后提取梯度方向上梯度强度局部最大的点为边缘点，在边缘点形成的脊中，提取脊的顶部，对脊像素做阈值处理，获得图像边缘，结果为二值图，获得背景部分和边缘部分；[0018]2.2将背景部分分割连通域，以四连通区域为一个连通域，按顺序进行编号；[0019]2.3遍历每个连通域，求取每个连通域的面积，取面积在0.2A〜1.2A范围内的连通域作为候选的缺陷可能区域，其中A为预设潜在缺陷大小。[0020]所述的步骤3具体为：[0021]3.1针对每个缺陷可能区域在原X射线图像上截取建立一个矩形的感兴趣区域，感兴趣区域的中心为缺陷可能区域的几何中心，感兴趣区域的长为感兴趣区域沿X轴方向最大尺寸的两倍，感兴趣区域的宽设为感兴趣区域沿y轴方向最大尺寸的两倍，并将感兴趣区域的尺寸调整为32X32像素；[0022]3.2经过感兴趣区域的中心建立八个方向，八个方向分别为与图像X轴或者y轴的正方向呈22.5°父1^角度的八个直线所在方向沽取0，1，2,3,4,5,6,7;[0023]3.3提取感兴趣区域在各方向直线上各个像素点的灰度信息，从八个方向中选取两端点灰度相差最小的方向及与其垂直的方向，分别作为第一纵横线和第二纵横线，将沿第一纵横线和第二纵横线上感兴趣区域各个像素点的灰度信息构建获得第一灰度向量Pl和第二灰度向量P2。[0024]本发明筛选出两条相互垂直的纵横线，筛选条件为：其中一条纵横线在八个方向纵横线中两端点灰度相差最小，筛选的另一条纵横线为该纵横线的垂直方向纵横线。[0025]八个方向分别为与图像X轴或者y轴的正方向呈22.5°Xk角度的八个直线所在方向。[0026]所述的步骤4具体为：[0027]4.1采用以下方式针对第一灰度向量Pl处理获得第一灰度值以消除背景影响，消除背景影响包括消除背景变化趋势和消除基础背景：[0028]绘制第一纵横线的灰度-坐标图，即对第一灰度向量Pl处理，获得第一纵横线所在方向灰度的变化斜率，表达式如下：ΓΊ,P,32-Ri\032-1[0030]式中，P132⑴分别为第一纵横线的灰度-坐标图中坐标最大和最小点对应的灰度值，ko表不第一纵横线的灰度-坐标图中坐标最大和最小点连线的斜率；[0031]然后采用以下公式针对感兴趣区域中第一纵横线上的像素点进行计算获得各点的第一灰度值：[0033]式中，Q1i为第一纵横线的灰度-坐标图中第i个点消除背景灰度影响后的灰度值，作为第一灰度值，Pii表不第一纵横线的灰度-坐标图中第i个点的灰度值，i表不灰度-坐标图中点的序号，i取1至32的整数；[0034]4.2采用以下方式针对第一灰度向量Pl和第二灰度向量P2处理获得第二灰度值以消除背景影响，消除背景影响包括消除背景变化趋势和消除基础背景：[0035]取第一灰度向量Pl中的第5至28个元组成第一背景向量PPl，取第二灰度向量P2中的第5至28个元素组成第二背景向量PP2,绘制第一背景向量PPl和第二背景向量PP2的灰度-坐标图，对第一背景向量PPl和第二背景向量PP2采用以下公式计算处理，求得各自的斜率ki和k2,表达式如下：[0038]式中，PP124、PPi⑴和PP224、PP2⑴分别为PPl和PP2灰度-坐标图中坐标最大和最小点对应的灰度值，斜率1^和1«分别表示第一背景向量PPl和第二背景向量PP2的两端点之间的斜率；[0039]然后采用以下公式消除第一背景向量PPl和第二背景向量PP2背景影响，获得各点的第二灰度值和第三灰度值，表达式如下：[0042]式中，PP1⑴表示第一背景向量PPl的灰度-坐标图中第i个坐标点的灰度值，Q2⑴为第一背景向量PPi的灰度-坐标图中第i个坐标点消除背景灰度影响后的灰度值，作为第二灰度值;PP2i表示第二背景向量PP2的灰度-坐标图中第i个坐标点的灰度值，Q3i为第二背景向量PP2的灰度-坐标图中第i个坐标点消除背景灰度影响后的灰度值，作为第三灰度值，i取1至24的整数；[0043]4.3提取轮廓特征[0044]计算获得两条纵横线的轮廓特征，轮廓特征包括将所有第一灰度值Qoi经快速傅里叶变换（FFT得到的一次谐波幅值F0、将所有第二灰度值Q2i经快速傅里叶变换FFT—次谐波幅值Fl和将所有第三灰度值Q3i经快速傅里叶变换（FFT—次谐波幅值F2;[0045]然后利用一次谐波幅值Fl和一次谐波幅值F2采用以下公式计算求得所有第二灰度值Q2⑴和所有第三灰度值Q3⑴的一次谐波幅值之间的相似度数SO，表达式如下：[0047]式中，Fl和F2分别为所有第二灰度值出（i和所有第三灰度值Q3i的快速傅里叶变换一次谐波幅值；[0048]4.3然后根据幅值FO和相似度数SO各自的上下限预设阈值判断获得感兴趣区域是否为真实缺陷的结果：[0049]若幅值FO和相似度数SO均在各自的上限预设阈值和下限预设阈值之间，则感兴趣区域为真实缺陷；否则感兴趣区域不为真实缺陷。[0050]本发明的快速傅里叶变换一次谐波幅值FO主要表征了一定频率范围的信号幅值，相似度数SO主要表征以一定频率范围内的两个信号FFT—次谐波幅值评判的相似程度，能反映出图像在两个方向上的灰度变化相似性。[0051]本发明的有益效果主要表现在：[0052]1采用灰度级形态学顶帽运算，突出人眼难以分辨的模糊缺陷区域，同时去除大部分复杂铸件结构，可改善缺陷边缘检测效果[0053]2采用垂直纵横线灰度变化轮廓的相似度作为特征值，可筛去顶帽无法筛除的复杂结构的角落区域，减小铸件角落区域结构的影响，有效提高最终缺陷检测率。附图说明[0054]图1是本发明方法的流程图。[0055]图2是实施例采集到典型潜在缺陷的复杂结构铸件图像。[0056]图3是实施例反锐化掩模算法得到的锐化图像。[0057]图4是实施例灰度级开运算结果图。[0058]图5是实施例灰度级顶帽运算结果图。[0059]图6是实施例Canny边缘检测图像。[0060]图7是实施例缺陷可能区域图像。[0061]图8是实施例的8个方向所在直线示意图。[0062]图9是实施例8个方向所在直线灰度信息图。[0063]图10是实施例复杂结构铸件潜在缺陷检测结果图像。具体实施方式[0064]下面利用附图和实施例对本发明基于垂直纵横线轮廓特征的复杂结构铸件潜在缺陷检测方法作进一步说明。[0065]如图1所示，本发明的实施例如下：[0066]步骤1:采用反锐化掩模算法对原始图像进行锐化，并采用灰度级形态学顶帽运算突出局部高亮度区域，实现图像预处理。过程如下：[0067]1.1先提取图像高频成分，创建大小为·酌高斯算子，其中A为预设潜在缺陷大小，设置标准偏差为4,对原图像高斯滤波，用原图减去高斯滤波的结果得到图像高频成分;再在原图像上加上图像高频成分的3〜6倍，得到锐化图像，如图3所示；[0068]1.2对锐化图像进行形态学顶帽处理:统计模糊缺陷的大小，设置一定大小的圆形结构体，圆形结构体的大小设定为预设潜在缺陷大小，对锐化灰度图像先进行腐蚀，再膨胀，即进行开运算，效果如图4所示;接着，用锐化图像减去膨胀之后的图像，获得顶帽运算结果图像。实现的效果是:令尺寸大的大部分铸件结构被消除，仅留下尺寸小的缺陷和噪声等部分，减小铸件结构对边缘检测的干扰，顶帽运算结果图像如图5所示；[0069]步骤2:采用Canny边缘检测，对边缘二值图像取反，标记连通域，并通过连通域像素数目初步筛选缺陷可能区域。具体过程如下：[0070]2.1边缘检测采用Canny边缘检测方法，先用标准偏差为2.5的高斯滤波器平滑图像，减小噪声干扰的影响，然后提取梯度方向上梯度强度局部最大的点为边缘点，在边缘点形成的脊中，提取脊的顶部，对脊像素做阈值处理，获得图像边缘，结果为二值图，背景部分灰度为0，边缘灰度为1，如图6所示；[0071]2.2对边缘二值图像取反，使背景部分灰度为1，边缘灰度为0;标记连通域，连通域标记为灰度为1的区域编号，设定四连通区域为一个连通域，编号相同；[0072]2.3遍历每个连通域，求取每个连通域的面积，即编号相同的像素数量，设定阈值，仅取连通域面积在一定范围内的连通域作为候选的缺陷可能区域，并重新为缺陷可能区域编号，缺陷可能区域如图7所示；[0073]步骤3:提取感兴趣区域，求取8个方向上过感兴趣区域中心的纵横线灰度信息，筛选出两条相互垂直的纵横线，分析两条纵横线灰度图像轮廓变化特征，表征轮廓变化相似程度的特征值在〇-〇.57范围内且表征轮廓变化幅度的特征值在15-300范围内的区域判定为真实缺陷。详细过程如下：[0074]3.1遍历缺陷可能区域，提取感兴趣区域，原始图像中截取矩形区域中的图像形成一新图像，矩形区域的中心设为连通域的中心坐标，矩形区域长设为连通域沿X方向尺寸的两倍，矩形区域宽设为连通域沿y轴方向尺寸的两倍，截取出新矩形区域图像后，将感兴趣区域尺寸调整为32X32;[0075]3.2提取图像感兴趣区域各方向中心线上的灰度信息，经过32X32矩形区域的中心，绘制与X轴正方向呈22.5°父1^角度的8个方向的直线，其中1^取0，1，2,3,4,5,6,7,8条直线如图8所示；[0076]提取8条直线经过像素的灰度信息，由这些灰度信息绘制灰度-坐标图，如图9所不，从k个方向中选取两端点灰度相差最小的方向和其垂直的方向，分别为第一纵横线和第二纵横线，本实施例中分别对应图9中的k=7方向和k=0方向，其对应的保存灰度信息的向量分别为第一灰度向量Pl和第二灰度向量P2;[0077]3.3首先，针对第一灰度向量Pl处理获得第一灰度值:从第一纵横线的灰度-坐标图，获得第一灰度向量Pl灰度的变化斜率，然后消除第一灰度向量Pl的背景影响，获得各点的第一灰度值。[0078]例如，求图9中第一灰度向量Pl灰度的变化斜率：[0080]例如，求图9中第一灰度向量Pl的第i=2点的第一灰度值：[0081]Q1⑵=P1⑵-P1I-Ico*2-1=100.32-97.41-0.497*1=2.413[0082]接着，针对第一灰度向量Pl和第二灰度向量P2处理获得第二灰度值:取第一灰度向量Pl中的第5至28个元组成第一背景向量PPl，取第二灰度向量P2中的第5至28个元素组成第二背景向量PP2,绘制第一背景向量PPl和第二背景向量PP2的灰度-坐标图，对第一背景向量PPl和第二背景向量PP2计算处理，求得各自的斜率kdPk2;然后消除第一背景向量PPl和第二背景向量PP2背景影响，获得各点的第二灰度值和第三灰度值。[0083]例如，求图9中第一背景向量PPl和第二背景向量PP2的变化斜率：[0086]例如，分别求图9中第i=2点的第二灰度值和第三灰度值：[0089]最后，计算获得两条纵横线的轮廓特征，轮廓特征包括所有第一灰度值Q1i的快速傅里叶变换（FFT得到的一次谐波幅值FO、所有第二灰度值如（i的快速傅里叶变换FFT得到的一次谐波幅值Fl和所有第三灰度值Q3i的快速傅里叶变换FFT得到的一次谐波幅值F2;然后利用Fl和F2采用以下公式计算求得所有第二灰度值Q2i和所有第三灰度值Q3⑴的FFT—次谐波幅值之间的相似度数S0。[0090]图9所示实施例最终得到的相似度数SO是：[0092]图9所示最终实施例最终得到的FO为129.17[0093]5设置上下限预设阈值15〈F0〈300且S0〈0.57时，判定缺陷可能区域为真实缺陷，比较情况是:图8、图9中表示的缺陷可能区域被判定为真实缺陷。[0094]本发明共统计70个典型样本的两个特征值情况，设置较合适的阈值，两个特征值各自在一定阈值以内的可能缺陷判断为真实缺陷并标记，对每个缺陷可能区域进行计算和判定，获得最终结果，如图10所示。[0095]通过上述的实例分析总结出，本发明方法提供了一种基于垂直纵横线轮廓特征的复杂结构铸件潜在缺陷检测方法，能有效实现真实缺陷的检测。上述对典型实例的说明并不用以限制本发明，在不背离本发明的情况下可以进行修改和改进，其范围在权利要求书以及等同物中进行了限定。

权利要求：1.一种基于垂直纵横线轮廓特征的复杂结构铸件潜在缺陷检测方法，其特征在于包括以下步骤：1采集复杂结构铸件的X射线图像，采用反锐化掩模算法对图像进行锐化，接着进行灰度级形态学顶帽运算，实现图像预处理；2采用Canny边缘检测，标记连通域，再通过连通域像素面积初步筛选获得缺陷可能区域；3提取感兴趣区域，求取八个方向上过感兴趣区域中心的纵横线灰度信息，筛选出两条相互垂直的纵横线；4分析两条纵横线灰度图像轮廓变化特征并消除背景影响，根据轮廓变化特征判定是否为真实缺陷。2.根据权利要求1所述的基于垂直纵横线轮廓特征的复杂结构铸件潜在缺陷检测方法，其特征在于:所述的步骤1具体为：1.1提取图像高频成分，利用高斯算子对原图像高斯滤波，用原图像减去高斯滤波的结果得到图像高频成分，再在原图像上加上图像高频成分的倍数，实现图像锐化，得到锐化图像；1.2对锐化图像进行形态学顶帽处理:根据预设的潜在缺陷大小，设置与潜在缺陷大小相同的圆形结构体，用圆形结构体对锐化灰度图像进行开运算，再用锐化图像减去开运算之后的图像，获得顶帽运算结果图像。3.根据权利要求1所述的基于垂直纵横线轮廓特征的复杂结构铸件潜在缺陷检测方法，其特征在于:所述的步骤2具体为：2.1边缘检测采用Canny边缘检测方法，先用带有标准偏差的高斯滤波器平滑图像处理，然后提取梯度方向上梯度强度局部最大的点为边缘点，在边缘点形成的脊中，提取脊的顶部，对脊像素做阈值处理，获得图像边缘，结果为二值图，获得背景部分和边缘部分；2.2将背景部分分割连通域，以四连通区域为一个连通域，按顺序进行编号；2.3遍历每个连通域，求取每个连通域的面积，取面积在0.2A〜1.2A范围内的连通域作为候选的缺陷可能区域，其中A为预设潜在缺陷大小。4.根据权利要求1所述的基于垂直纵横线轮廓特征的复杂结构铸件潜在缺陷检测方法，其特征在于:所述的步骤3具体为：3.1针对每个缺陷可能区域在原X射线图像上截取建立一个矩形的感兴趣区域，感兴趣区域的中心为缺陷可能区域的几何中心，感兴趣区域的长为感兴趣区域沿X轴方向最大尺寸的两倍，感兴趣区域的宽设为感兴趣区域沿y轴方向最大尺寸的两倍，并将感兴趣区域的尺寸调整为32X32像素；3.2经过感兴趣区域的中心建立八个方向，八个方向分别为与图像X轴或者y轴的正方向呈22.5°父1^角度的八个直线所在方向沽取0，1，2,3,4,5,6,7;3.3提取感兴趣区域在各方向直线上各个像素点的灰度信息，从八个方向中选取两端点灰度相差最小的方向及与其垂直的方向，分别作为第一纵横线和第二纵横线，将沿第一纵横线和第二纵横线上感兴趣区域各个像素点的灰度信息构建获得第一灰度向量Pl和第二灰度向量P2。5.根据权利要求1所述的基于垂直纵横线轮廓特征的复杂结构铸件潜在缺陷检测方法，其特征在于:所述的步骤4具体为：4.1采用以下方式针对第一灰度向量Pl处理获得第一灰度值：绘制第一纵横线的灰度-坐标图，获得第一纵横线所在方向灰度的变化斜率，表达式如下：式中，Pi32和Pi1分别为第一纵横线的灰度-坐标图中坐标最大和最小点对应的灰度值，ko表示第一纵横线的灰度-坐标图中坐标最大和最小点连线的斜率；然后采用以下公式针对感兴趣区域中第一纵横线上的像素点进行计算获得各点的第一灰度值：Qii=Pi⑴-Pil-k〇*i-l式中，Q1i为第一纵横线的灰度-坐标图中第i个点消除背景灰度影响后的灰度值，作为第一灰度值,P1i表示第一纵横线的灰度-坐标图中第i个点的灰度值，i表示灰度-坐标图中点的序号，i取1至32的整数；4.2采用以下方式针对第一灰度向量Pl和第二灰度向量P2处理获得第二灰度值：取第一灰度向量Pl中的第5至28个元组成第一背景向量PPl，取第二灰度向量P2中的第5至28个元素组成第二背景向量PP2,绘制第一背景向量PPl和第二背景向量PP2的灰度-坐标图，对第一背景向量PPl和第二背景向量PP2采用以下公式计算处理，求得各自的斜率Iu和k2,表达式如下：式中，PPi24、PPi1和PP224、PP21分别为PPl和PP2灰度-坐标图中坐标最大和最小点对应的灰度值，斜率1^和1«分别表示第一背景向量PPl和第二背景向量PP2的两端点之间的斜率；然后采用以下公式消除第一背景向量PPl和第二背景向量PP2背景影响，获得各点的第二灰度值和第三灰度值，表达式如下：式中，PP1i表示第一背景向量PPl的灰度-坐标图中第i个坐标点的灰度值，Q2i为第一背景向量PPl的灰度-坐标图中第i个坐标点消除背景灰度影响后的灰度值，作为第二灰度值;PP2⑴表示第二背景向量PP2的灰度-坐标图中第i个坐标点的灰度值，Q3i为第二背景向量PP2的灰度-坐标图中第i个坐标点消除背景灰度影响后的灰度值，作为第三灰度值，i取1至24的整数；4.3提取轮廓特征计算获得两条纵横线的轮廓特征，轮廓特征包括将所有第一灰度值Qoi经快速傅里叶变换FFT得到的一次谐波幅值FO、将所有第二灰度值Q2i经快速傅里叶变换FFT—次谐波幅值Fl和将所有第三灰度值Q3⑴经快速傅里叶变换FFT—次谐波幅值F2;然后利用一次谐波幅值Fl和一次谐波幅值F2采用以下公式计算求得所有第二灰度值Q2⑴和所有第三灰度值Q3⑴的一次谐波幅值之间的相似度数S0,表达式如下：式中，Fl和F2分别为所有第二灰度值〇2i和所有第三灰度值Q3i的快速傅里叶变换一次谐波幅值；4.3然后根据幅值FO和相似度数SO各自的上下限预设阈值判断获得感兴趣区域是否为真实缺陷的结果：若幅值FO和相似度数SO均在各自的上限预设阈值和下限预设阈值之间，则感兴趣区域为真实缺陷；否则感兴趣区域不为真实缺陷。

百度查询： 浙江大学 基于垂直纵横线轮廓特征的复杂结构铸件潜在缺陷检测方法