1、道BICS 19. 100 J 04 中华人民共和国国家标准GB/T 29067-2012 工业计算机层析成像(CT)图像测量方法无损检测Non-destructive testing-Test method for measuring industrial computed tomography (CT) image 龟飞胁g 2013-10-01实施2012-12-31发布发布中华人民共和国国家质量监督检验检夜总局中国国家标准化管理委员会队孔子宫/、自量阳、曲曲层ff，哈制 、4 . GB/T 29067-2012 目。本标准按照GB/T1. 1-2009给出的规则起草。本标准由全国元损检
2、测标准化技术委员会(SAC/TC56)提出并归口。本标准起草单位z重庆真测科技股份有限公司、重庆大学ICT研究中心、上海泰司检测科技有限公司、中国石油集团济柴动力总厂、齐齐哈尔轨道交通装备有限责任公司、上海航天动力技术研究所。本标准主要起草人z曾理、王在、李树生、蔡玉芳、沈宽、邹永宁、王福全、段黎明、杨蓉、高金生、俞翔、沈伟民。I GB/T 29067-2012 无损检测工业计算机层析成像(CT)图像测量方法1 范围本标准规定了用于元损检测的工业计算机层析成像(CT)系统图像密度测量校准和尺寸测量校准的方法，以及从CT图像测量被测物体物质密度和尺寸的方法。本标准适用于测量被测物体内部各部分的密
3、度(包括由不同材料组成的被测物体内部各部分的密度)及被测物体内外结构(特别是封闭内腔)的几何尺寸。2 规范性引用文件下列文件对于本文件的应用是必不可少的。凡是注日期的引用文件，仅注日期的版本适用于本文件。凡是不注日期的引用文件，其最新版本(包括所有的修改单适用于本文件。GBjT 12604. 2元损检测术语射线照相检测3 术语和定义GBjT 12604. 2界定的术语和定义适用于本文件。4 应用基础4. 1 检测设备要求要求具备相关的检测设备(包括工业CT设备及相关图像处理软件)，并要求客户和供应商在相关方面达成一致。4.2 避免图像伪影对密度测量的影响密度测量方法要求:CT密度的变化反映被测
4、物体密度的真实变化而不是图像伪影。由于探测器响应不一致、射线散射、射线硬化等原因，CT图像可能呈现出一定程度的伪影，且伪影往往被误认为是密度变化。为确保该密度测量方法的有效性，工业CT设备制造商应对射线投影数据进行校正，以尽量减少CT图像伪影。在实际检测中，用户应该尽量避开受伪影影响的区域进行密度测量。4.3 避免图像伪影对尺寸测量的影响尺寸测量方法要求E被测物体的被测区域对应的CT图像应具有比较陡峭的边缘。由于探测器响应不一致、射线散射、射线硬化等原因，CT图像可能呈现出一定程度的伪影，加之被测物体的直线型边缘可能会造成条状伪影。为确保该尺寸测量方法的有效性，工业CT设备制造商应对射线投影数
5、据进行校正，以尽量减少CT图像伪影。在实际检测中，用户应该尽量避开受伪影影响的区域进行尺寸测量。1 GB/T 29067-2012 5 标准试件的制作5. 1 密度标准试件的制作5. 1. 1 若客户与供应商之间没有其他的协议，密度标准试件组应由多个不同密度的圆柱形密度标准试件组成，在基体上沿圆周均匀分布，如图1所示。截面./-A 圄1密度标准试件组5. 1.2 密度标准试件组的材料密度应涵盖被测物体材料密度范围，且每个标准试件的密度已知。每个标准试件应材质均匀，各个标准试件的尺寸规格应相同。5. 1.3 密度标准试件高度庇为CT切片厚度的3倍以上。5.1.4 密度标准试件圆柱直径应该确保每个
6、标准试件CT切片图像的像素个数不少于100。5.2 尺寸标准试件的制作5.2. 1 尺寸标准试件宜选择与被翻物体材料一致的试块，也可选用标准量块。5.2.2 尺寸标准试件的偏差应小于被测物体的测暨偏差要求至少一个数量级，例如被测物体的测量偏差要求小于0.1mm，则标准试件尺寸的偏差应小于0.01mm。5.2.3 尺寸标准试件的形状应为长方体。设系统的空间分辨力为K(单位为LP/mm)，尺寸标准试件的最长边的长度应大于50X才Kmm，且对应CT图像的最长边宜具有50个或以上的像素，以减少测量的系统误差。6 校准与测量6. 1 标准试件和被测物体的检测工艺要求6. 1. 1 将标准试件或被测物体安
7、置在CT系统转台上，使得试件中心轴垂直于扫描基准平面。6. 1.2 标准试件应尽量安置在转台上与被测物体的待测区域大致相同的位置，或与被测物体的待测区域到转轴的距离近似相同。4 GB/T 29067-2012 6. 1.3 标准试件扫描切片位置应大致选择在标准试件的中间位置，被测物体扫描切片位置根据需要而定。6. 1.4 标准试件与被测物体应采用相同的工艺参数(包括射线源参数、扫描参数、重建参数等及图像数据处理方法。6.2 密度测量校准6.2. 1 密度测量校准时宜将扫描参数、重建参数调整至系统密度分辨力最优的参数设置下。6.2.2 在CT图像上的每个密度标准试件中心部位分别选一圆形区域(每个
8、圆形区域像素数宜多于50个，但不大于整个标准试件区域的2/白，测得该区域内每一点的CT密度f(x，计算该区域的平均CT密度f。6.2.3 选取密度标准试件组中一个密度标准试件作为密度测量授准的基准标准试件。根据基准标准试件指定区域的平均CT密度70和CT系统设定的比例系数k(根据CT系统总体参数确定，一般为1 0005 000) ，对每个密度标准试件指定区域的平均CT密度f按式(1)进行归一化处理，获得该区域内的平均CT数ji:式中z(-f =J_JOXk fo f一一密度标准试件指定区域的平均CT密度FF一一密度标准试件指定区域的平均CT数zfO一-基准标准试件指定区域的平均CT密度zh一一
9、-CT系统设定的比倒系数。. ( 1 ) 6.2.4 利用各个密度标准试件指定区域的平均CT数豆和l该标准试件己知的密度p，采用多项式拟合校准平均CT数F与物质密度p之间的关系式:p=f(，推荐采用三次多项式拙合。6.3 被测物体的密度测量6.3. 1 选取与密度校准时采用的同一基准括准试件，将基准标准试件放在被割物体的待测区域附近与被测物体同时进行CT扫描z或将基准标准试件单独放在与被测物体的待测区域丁同一扫描位置，采用与被测物体扫描时相同的工艺参数，对密度标准试件单独进行扫描。6.3.2 从基准标准试件的CT图像中，获取指定区域(指定区域内像素数宜多于50个且不大于整个标准试件区域的2/3
10、)的平均CT密度foo6.3.3 从被测物体的CT图像中，获取待测区域区域内像素数宜多于50个的平均CT密度f。按式(2)进行归一化得到被测物体的待测区域的平均CT数L(-f, =J _JO Xk 式中zf一一被测物体的待测区域的平均CT密度55一一被测物体的待测区域的平均CT数EfO fO-基准标准试件的指定区域的平均CT密度zh一一CT系统设定的比例系数(与密度测量校准时相同)。.( 2 ) 6.3.4 根据被测物体的待测区域的平均CT数F和已校准的函数关系p=f(，求得该区域平均物质密度p.6.3.5 重复测量同一均匀材质不同区域至少5次，以5次测得密度的平均值作为此材质的密度值，并计算
11、密度测量的标准偏差。3 G/T 29067-2012 6.4 尺寸测量投准6.4. 1 尺寸标准试件校准时应将扫描参数、重建参数调整至系统空间分辨力最优的参数设置下。6.4.2 测试时，宜使尺寸标准试件对射线的衰减系数与实际被测物体接近，必要时可采用射线衰减补偿。6.4.3 如图2所示，设定尺寸标准试件1为校准测量的标准试件，尺寸标准试件2为测试测量误差的标准试件。在CT图像上尺寸标准试件1短边中心处，沿长边方向画一条直线。将直线所经过点的CT数绘制在坐标图上则形成该标准试件的CT数曲线。由于尺寸标准试件和空气的密度相差很大，因此该曲线会形成一个边缘陡峭的波峰，其理想曲线如图3所示。6.4.4
12、 由于噪声的干扰，标准试件的CT数曲线往往不能形成理想曲线。因此，建议沿长边方向画5条平行直线，将5条直线所经过点的CT数绘制在坐标图上，并进行最小二乘拟合，以获得近似的理想CT数曲线。6.4.5 在该波峰高度约一半的位置定位一测量直线P，该直线与CT数曲线相交于A、B两点，在坐标系中找出A、B点对应的横坐标值，由两者差的绝对值可得到尺寸标准试件1的图像长度D(单位z像素。6.4.6 根据标准试件1的图像长度D(单位z像素、图像矩阵的宽度W(单位z像素)、CT系统扫描视场成像直径V(单位z毫米)，由式(3)计算出标准试件1的CT几何长度测量值L(单位z毫米): L=DX号. . . . . .
13、 .( 3 ) 6.4.7 比较标准试件1的几何长度测量值L与实际长度值之间的差异，如果测量值大于实际值则将测量直线P向上移动直至两者相等，反之则将直线P向下移动。记录此时测量直线P所处位置的CT数件。尺寸标件2圈2尺寸标准试件组CT示意图4 单G/T 29067-2012 CT数RB - r -_ -_ 一一二可二-一一一一一一一一A P B D 4二l民一像索个数图3圄2中标准试件1直线位置的理想CT数曲线图6.4.8 计算曲线上低密度物质对应CT数的均值L即图3中RA和RC所指曲线段的CT数平均值;计算曲线上高密度位置对应CT数的均值阳，即图3中RB所指曲线段的CT数平均值;计算测量直线
14、P在CT数曲线陡坡上的相对位置比例Cp，Cp由式(4)计算得到z口二F C ( 4 ) H一严L6.4.9 在图像上尺寸标准试件2短边中心处，沿长边方向画一条直线。将该直线所经过各个点的CT值绘制在坐标图上则形成类似图3的CT数曲线。同样在该曲线上定位一直线Q，调整Q的位置使得Q所处位置的CT数等于阳，向通过式。)得到z Q-L十CpX (H一L)(5 ) 则直线Q与CT数曲线相交两点之间横坐标差的绝对值，即为本次测量的尺寸标准试件2的CT图像长度(单位z像素)。6.4. 10 重复测量3次，依次计算出标准试件2的图像长度Dl鸟，D3o由式(6)计算取三者的平均值作为尺寸标准试件2的CT图像像
15、素长度测量值D(单位z像素)，并计算标准偏差。4 D ( 6 ) 6.4. 11 依据式(3)计算出标准试件2的几何长度测量值L。由标准试件2的CT几何长度测量值与其实际长度之差的绝对值，得到标件2的CT长度L测量误差的绝对值:a=IL测量值一L实际值|。6.4. 12 由于该方法是基于图像中物质CT数曲线，因此在测量线上需消除或避开图像伪影。6.4. 13 为了克服本测量方法中由于人工划线带来的测量结果的不稳定性，图像处理系统可以采用软件自动测量只寸标准试件两条平行短边的距离，代替手动测量平均值。6.5 被测物体的人工尺寸测量6.5. 1 测量被测物体时，建议将被测物体与其密度相近的尺寸标准
16、试件一起做CT扫描或采用相同的工艺参数扫描。6.5.2 采用和前述尺寸标准试件技准时测量标准试件2尺寸类似的步骤和方法，测量被测物体尺寸。先测量标准试件的长度，确定标准试件测量线的相对位置比例Cpo然后测量被测物体待测区域的长5 GB/T 29067-2012 度，根据Cp的值调整被测物体测量线的位置。计算被测物体测量线与CT数曲线相交两点的横坐标差的绝对值，然后根据CT视场直径求出两点的实际长度。6.5.3 由于被测物体待测区域周围的材料密度可能不一致，导致上述CT数曲线左右高度可能不一致。此时可根据测量线的相对位置比例Cp，分别在CT数曲线的左右陡边中确定相交点，并计算两点的横坐标之差。6
17、.5.4 由于CT系统的固有特性和空间分辨力的限制，CT图像边缘存在扩散现象。该测量方法可有效地克服CT图像边缘扩散对测量精度的影响，但要求被测物体待测区域的长度应大于边缘扩散带宽的5倍以上。6.6 被测物体的自动尺寸测量6.6. 1 对管道(包括被测物体内腔及内部管道)壁厚的自动测量，可用图像处理软件自动提取内外壁边缘点并跟踪，自动计算内壁上各点到外壁的最短距离，从而得到壁厚变化曲线，最大、最小和平均壁厚等几何参数。6.6.2 对于已知CAD模型的被测物体，可采用(二维或三维)CAD模型和CT图像的尺寸比对自动测量技术，自动计算并直观显示被测物体的尺寸误差，以便于改进制造工艺。包括自动配准(
18、根据CT视场直径缩放CT图像或CAD图形，CT图像边缘提取，CT图像边缘与CAD图形线的平移旋转配准)和自动尺寸误差测量(计算CT图像边缘点到最近的CAD图形线的距离，并用不同的颜色或灰度显示被测物体边缘点的误差大小)等。7 精度和偏差7.1 密度测量的精度由密度测量最大误差除以系统密度测量量程，可得到系统密度测量(相对于量程的)精度(引用误差。密度测量精度与系统密度分辨力等因素有关。有实例显示z若系统密度分辨力达到千分之三，且CT图像元明显伪影，CT系统密度测量精度可达到1%左右。7.2 人工尺寸测量的偏差具体偏差和系统分辨力等因素有关。设系统的空间分辨力为K(单位为LP/mm)，则CT系统
19、尺寸测量偏差可小于l/(2XK)mmo系统的空间分辨力越高，CT尺寸测量偏差越小。人工划线结合CT数曲线自动定位的尺寸测量方法，可有效地消除数字图像按像素采样和CT图像边缘扩散带来的误差，使得测量偏差可小于一个像素代表的实际长度。7.3 自动尺寸测量的偏差自动尺寸测量(包括壁厚自动测量、CAD图形和CT图像的尺寸比对自动测量效率高、测量结果稳定。自动尺寸测量的偏差取决于边缘点的定位偏差等因素，而CT图像存在的像素采样误差和边缘扩散影响了边缘点的定位偏差。自动尺寸测量的偏差一般大于或等于一个像素代表的实际长度，采用亚像素的边缘点定位方法可减少边缘定位偏差和自动尺寸测量偏差。8 报告报告应包括重建
20、参数、显示参数、相关图像和测量数据。其中密度测量部分至少要包括指定材料的密度测量值，尺寸测量部分至少要包括指定位置(标示于CT图像上)的尺寸测量值。客户和供应商可协商其他需要包括的参数和数据。、-、. GB/T 29067-2012 9 测量实例9. 1 密度控准和测量实例根据图1制作密度标准试件组，密度标准试件物质分别为铜、钢、铁合金、铝合金、塑料和有机玻璃。根据各标准试件指定区域平均CT数户=6872，6385，3549，2302，1395，1011J及密度Pi= 8. 370, 7.850, 4.430 , 2.800 , 1.590, 1. 190J拟合多项式并进行测量，结果如表1所示
21、，系统测量精度(引用误差)在1%之内。表1实际材料密度测量对照表编号材料材料实际密度测量密度绝对误差精度(引用误差g/cm3 g/cm3 g/cm3 % 1 铜8.370 8.400 0.030 0.358 2 钢7.850 7.845 一0.005一0.0603 铁合金4.430 4.420 一0.010一0.1194 铝合金2.800 2. 770 一0.030一0.3585 塑料1. 590 1. 610 0.020 0.239 6 有机玻璃1. 190 1. 200 0.010 0.119 注z系统测量精度(引用误差计算是以标准试件组中最大密度材料的密度(铜:8.370g/cm3)为系
22、统密度测量量程。9.2 尺寸技准和测量实例以9MeV工业CT系统为例，该系统的空间分辨力为2LP/mm，分别选取长度为85.000mm和65.000 mm的两个钢质3级标准量块为标准试件1和标准试件2，选取视场直径为275.000mm，图像大小为1024X1 024，按照本标准方法进行CT图像尺寸校准和测量。首先确定标准试件1测量线位置P并计算出测量线的相对位置比例Cp，然后根据Cp计算出标准试件2测量线所在位置，从而计算出标准试件2的图像长度。3次测量标准试件2的图像长度分别为242.240，242.118，242.126(单位:像素)，3次测量均值为242.161(单位z像素)，转换后标准
23、试件2的CT几何长度测量值为(单位z毫米): 242.161 X告茫=65肌则该次测量的测量值与实际值的偏差为65H川阳刚mm)。NFON|homNH阁。,. . 华人民共和国家标准无损检测工业计算机层析成像(CT)图像测量方法GB/T 29067-2012 国中 中国标准出版社出版发行北京市朝阳区和平里西街甲2号(100013)北京市西城区三里河北街16号(100045)网址总编室:(010)64275323发行中心:(010)51780235读者服务部:(010)68523946中国标准出版社秦皇岛印刷厂印刷各地新华书店经销* 印张o.75 字数15千字2013年5月第一次印刷开本880X1230 1/16 2013年5月第一版 书号:155066. 1-47019定价16.00元如有印装差错由本社发行中心调换版权专有侵权必究举报电话:(010)68510107打印H期:2013年5月17日F002A
