工业数字射线成像检测散射处理方法技术

本发明专利技术公开了一种工业数字射线成像检测散射处理方法。它解决了现有散射处理方法复杂、运行速度慢，依赖具体成像条件，不适于工程应用的问题，具有运行速度快、算法有效性强、工程应用价值大，适于工业Ｘ－射线数字成像检测等优点。其方法为：（１）获取原始图像Ａ（ｍ，ｎ），其中ｍ、ｎ为像素点的横纵坐标值，并读取图像长度ｈ和图像宽度ｗ；（２）在空间域建立一个模板Ｍ，该模板与图像点对点平滑运算，得到散射图像，其中模板尺寸为：ｈ↓［０］×ｗ↓［０］＝ｈ／１６×ｗ／１６（ｈ↓［０］，ｗ↓［０］通常取二者中最小的奇数，如为偶数，则近似取），ｈ为图像长度，ｗ为图像宽度，模板系数采用等平均权重分布；（３）定义被检测工件厚度不均匀补偿系数ｒ，与散射图像Ｓ’的灰度值相乘，再与原始图像相减即可。

【技术实现步骤摘要】

专利技术涉及一种成像散射处理方法，尤其涉及一种。
技术介绍
目前在工业X-射线的数字成像检测中，散射是降低图像质量的一个重要因素。采用图像处理手段减低图像散射是数字化检测技术的优势。在现有的处理方法中，通常是通过获取系统的点扩展函数PSF(Point Spread Function)，建立系统散射模型，这种处理方法复杂、运行速度慢、依赖于具体的成像条件，不适于工程应用。工业中非常需要一种高速、有效的散射处理方法。
技术实现思路
本专利技术的目的就是为了解决现有的散射处理方法复杂、运行速度慢，依赖具体成像条件，不适于工程应用的问题，提供一种，它是一种基于空间域的工业X-射线数字成像检测散射处理方法，具有运行速度快、算法有效性强、工程应用价值大，适于工业X-射线数字成像检测等优点。为实现上述目的，本专利技术采用如下技术方案一种工业X-射线数字成像检测散射处理方法，它的方法为，(1)获取的原始图像A(m，n)，其中m、n为像素点的横纵坐标值，并读取图像长度h和图像宽度w；(2)在空间域建立一个模板M，该模板与图像点对点平滑运算，得到散射图像，其中模板尺寸为h0×w0=h16×w16]]>(h0，w0通常取二者中最小的奇数，如为偶数，则近似取)，h为图像长度，w为图像宽度，模板系数采用等平均权重分布；(3)定义被检测工件厚度不均匀补偿系数r，与散射图像S’的灰度值相乘，再与原始图像相减，便可达到散射处理的目的。所述步骤(2)中，模板系数为在模板尺寸的窗口内，采用等平均权重分布 所述步骤(2)中，模板M与原始图像的每一点A(m，n)卷积运算，运算公式为S′(m,n)=1h0×w0Σm=1h0Σn=1w0A(m,n).]]> 所述步骤(3)中，被检测工件厚度补偿系数r是一个在0～1之间的数，厚度均匀的工件，r在0.3～0.5中取值；厚度明显不均匀，则r在0.8～0.9之间取值。所述步骤(3)中根据被检测工件厚度，确定散射补偿系数r，与散射图像S’中的每一点相乘，运算公式为S＝r×S′(m，n)。所述步骤(3)，散射图像S与原始图像中的每一点进行C(m，n)＝A(m，n)-S(m，n)运算，完成散射处理。本专利技术处理方法包括一个空间域的模板M，被检测工件厚度补偿系数r、卷积运算、乘法运算和减法运算。空间域的模板有两个参数模板尺寸和模板系数。模板尺寸基于试验数据，设置为h0×w0=h16×w16]]>(h为图像长度，w为图像宽度，h0，w0通常取二者中最小的奇数，如为偶数，则近似取)；模板系数在模板尺寸的窗口内，采用等平均权重分布，即为 被检测工件厚度补偿系数r是一个在0～1之间的数，由于工件厚度的均匀分布程度，会引起图像中散射信息的变化。对于厚度分布均匀的被检测工件，如板类(焊缝等)，r在0.3～0.5中取值；厚度明显不均匀，则r在0.8～0.9之间取值。空间域的模板和原图像A实行卷积运算，得到散射图像S’，即S′＝AM被检测工件厚度分布补偿系数r与散射图像S’进行乘法运算，即S＝r×S′散射图像S与原始图像A进行减法运算，即可实现了散射处理，即C＝A-S本专利技术的有益效果由于采用了上述的处理方法，散射得到很好的处理；被检测工件厚度散射补偿系数的设置，使算法的边缘失真降低；空间域的卷积和四则运算结构，更是提高了运行速度，满足了工业中高效、高质的要求。附图说明图1为原始图像A；图2为散射图像S’；图3为散射处理后图像C。具体实施例方式下面结合附图与实施例对本专利技术作进一步说明。首先获取的原始图像A，读出图像的尺寸h、w，根据h、w计算模板的尺寸h0×w0=h16×w16,]]>确定模板中的等权重系数 然后确定了尺寸和系数的模板M与原始图像的每一点A(m，n)卷积运算，运算公式为S′(m,n)=1h0×w0Σm=1h0Σn=1w0A(m,n);]]>再根据被检测工件厚度特点，确定散射补偿系数r，与散射图像S’中的每一点相乘，运算公式为S＝r×S′(m，n)；最后，散射图像S与原始图像中的每一点进行C(m，n)＝A(m，n)-S(m，n)运算，便完成了的散射处理算法。如图像A，见图1为在基于Paxscan4030探测器的成像系统上获得的航空发动机涡轮叶片图像(透照参数为管电压110Kv，管电流1.8mA，焦距1160mm，滤波参数1mm黄铜等)。图像尺寸h×w＝393像素×350像素，散射处理的模板尺寸为h0×w0=h16×w16=24×21]]>(单位为像素)，按照步骤(2)，h0×w0＝21×21，模板系数为 按照S′(m,n)=1h0×w0Σm=1h0Σn=1w0A(m,n)]]>运算，得到的图像为S’，见图2，由于航空发动机涡轮叶片厚度分布不均匀，厚度最小处为2mm，最大处40mm，散射补偿系数r取0.9，按照S＝r×S′(m，n)运算图像，得到散射图像S，按照C(m，n)＝A(m，n)-S(m，n)，便完成了散射处理，得到了抑制散射后的图像，见图3。权利要求1.一种，其特征是(1)获取的原始图像A(m，n)，其中m、n为像素点的横纵坐标值，并读取图像长度h和图像宽度w；(2)在空间域建立一个模板M，该模板与图像点对点平滑运算，得到散射图像，其中模板尺寸为h0×w0=h16×w16,]]>其中h0，w0取二者中最小的奇数，如为偶数，则近似取，h为图像长度，w为图像宽度，模板系数采用等平均权重分布；(3)定义被检测工件厚度不均匀补偿系数r，与散射图像S’的灰度值相乘，再与原始图像相减，便可达到散射处理的目的。2.根据权利要求1所述的工业X-射线数字成像检测散射处理方法，其特征是所述步骤(2)中，模板系数为在模板尺寸的窗口内，采用等平均权重分布3.根据权利要求1所述的工业X-射线数字成像检测散射处理方法，其特征是所述步骤(2)中，模板M与原始图像的每一点A(m，n)卷积运算，运算公式为S′(m,n)=1h0×w0Σm=1h0Σn=1w0A(m,n).]]>4.根据权利要求1所述的工业X-射线数字成像检测散射处理方法，其特征是所述步骤(3)中，被检测工件厚度补偿系数r是一个在0～1之间的数，厚度均匀的工件，r在0.3～0.5中取值；厚度明显不均匀，则r在0.8～0.9之间取值。5.根据权利要求1所述的工业X-射线数字成像检测散射处理方法，其特征是所述步骤(3)中根据被检测工件厚度，确定散射补偿系数r，与散射图像S’中的每一点相乘，运算公式为S＝r×S′(m，n)。6.根据权利要求1所述的工业X-射线数字成像检测散射处理方法，其特征是所述步本文档来自技高网...

【技术保护点】
一种工业数字射线成像检测散射处理方法，其特征是：（１）获取的原始图像Ａ（ｍ，ｎ），其中ｍ、ｎ为像素点的横纵坐标值，并读取图像长度ｈ和图像宽度ｗ；（２）在空间域建立一个模板Ｍ，该模板与图像点对点平滑运算，得到散射图像，其中模板 尺寸为：ｈ↓［０］×ｗ↓［０］＝ｈ／１６×ｗ／１６，其中ｈ↓［０］，ｗ↓［０］取二者中最小的奇数，如为偶数，则近似取，ｈ为图像长度，ｗ为图像宽度，模板系数采用等平均权重分布；（３）定义被检测工件厚度不均匀补偿系数ｒ，与散射图像Ｓ’的 灰度值相乘，再与原始图像相减，便可达到散射处理的目的。

【技术特征摘要】

【专利技术属性】
技术研发人员：孔凡琴，
申请(专利权)人：孔凡琴，陈仁富，
类型：发明
国别省市：88[中国|济南]