发明名称 工业CT检测中小细节特征尺寸测量方法
摘要 本发明涉及一种工业CT检测中小细节特征尺寸测量方法,在扫描获取的检测工件CT图像中,选取包含细节特征的细节特征区域以及无明显缺陷的参考区域,分别对细节特征区域和参考区域内的像素CT值进行计算,然后将细节特征区域内每个像素的CT值相对于参考区域内所有像素的平均CT值的差,与细节特征区域内每个像素的等效标准偏差值进行比例计算,当比例值大于3时,则认定该像素为细节特征点,进而累加细节特征点的数量,从而获取细节特征点的面积,经过换算后获取等效圆直径,减去有效射束宽度值即可获得细节特征区域的直径。该工业CT检测中小细节特征尺寸测量方法应用灵活,识别度和测量精度高。
申请公布号 CN105092616A 申请公布日期 2015.11.25
申请号 CN201510562522.5 申请日期 2015.09.07
申请人 中国兵器科学研究院宁波分院 发明人 张维国;倪培君;郭智敏;徐向群;齐子诚;乔日东;左欣
分类号 G01N23/04(2006.01)I 主分类号 G01N23/04(2006.01)I
代理机构 宁波诚源专利事务所有限公司 33102 代理人 袁忠卫;王莹
主权项 一种工业CT检测中小细节特征尺寸测量方法,其特征在于包括如下步骤:步骤一、开启工业CT系统并进行校准标定;步骤二、设定工业CT系统中的工艺扫描参数,包括重建矩阵H×H、重建范围D、射线源焦点尺寸a、单个探测器通道尺寸d、探测器通道间隔h、射线源到探测器距离L,射线源到转台中心距离q,工件材料等效钢系数α、射线等效能量数据,根据设定的工艺扫描参数对检测工件进行扫描;步骤三、工业CT系统生成检测工件的CT图像,观察CT图像中的各种细节特征,确定需要进行尺寸测量的细节特征;步骤四、在CT图像上选定能够包含一个细节特征且包括边界在内的细节特征区域W,统计计算细节特征区域内每个像素的CT值f(x<sub>i</sub>,y<sub>j</sub>),其中i,j∈W,确定细节特征区域内像素的最小CT值Min对应的像素位置(x<sub>min</sub>,y<sub>min</sub>);在CT图像中细节特征的附近,且无明显缺陷的工件图像上选择一参考区域C,统计计算参考区域内所有像素的CT值,确定参考区域内像素的最大CT值及最小CT值,计算参考区域内所有像素的平均CT值Aver,进而计算参考区域内像素的噪声标准偏差值S<sub>d</sub>;步骤五、识别CT图像中的工件边界,以细节特征区域中像素的最小CT值Min对应的像素位置(x<sub>min</sub>,y<sub>min</sub>)为中心,自标定的0°方向开始旋转一周至360°方向,利用等效钢系数α、射线等效能量数据通过试验的方法,计算确定各角度A方向上与工件的等效钢厚度相关的标准偏差系数T<sub>A</sub>,0°≤A<360°,进而形成等效钢厚度数组;步骤六、按照步骤五中角度方向的标定,计算细节特征区域W内每个像素位置与标定0°方向间的夹角的角度A(x<sub>i</sub>,y<sub>j</sub>),然后计算细节特征区域W内每个像素对应的噪声参数值K:<img file="FDA0000797118060000011.GIF" wi="484" he="162" />其中i,j∈W,<img file="FDA0000797118060000013.GIF" wi="152" he="78" />为像素位置为(x<sub>i</sub>,y<sub>j</sub>)对应的角度方向A(x<sub>i</sub>,y<sub>j</sub>)上与等效钢厚度相关的标准偏差系数;设定缺陷点数的初始值num<sub>0</sub>=0,当K>n时,缺陷点数num累计加一,否则丢弃,其中n为经验噪声参数值;步骤七、根据工艺扫描参数计算单个像素面积的大小S以及有效射束宽度BW;<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>S</mi><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mi>D</mi><mi>H</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>;</mo></mrow>]]></math><img file="FDA0000797118060000012.GIF" wi="249" he="154" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>B</mi><mi>W</mi><mo>&ap;</mo><mfrac><msqrt><mrow><msup><mi>d</mi><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>&lsqb;</mo><mi>a</mi><mrow><mo>(</mo><mi>M</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup></mrow></msqrt><mi>M</mi></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000797118060000021.GIF" wi="540" he="159" /></maths>其中M=L/q;进而计算获取细节特征区域的直径T<sub>d</sub>:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>T</mi><mi>d</mi></msub><mo>=</mo><mn>2</mn><mo>*</mo><msqrt><mfrac><mrow><mi>n</mi><mi>u</mi><mi>m</mi><mo>*</mo><mi>S</mi></mrow><mi>&pi;</mi></mfrac></msqrt><mo>-</mo><mi>B</mi><mi>W</mi><mo>.</mo></mrow>]]></math><img file="FDA0000797118060000022.GIF" wi="536" he="148" /></maths>
地址 315103 浙江省宁波市高新区凌云路199号