基于双周外部ct检测装置在线检测大口径管壁的方法及系统的制作方法
【专利摘要】本发明涉及一种基于双周外部CT检测装置在线检测大口径管壁的方法及系统，该系统包括射线源、面阵探测器Ⅰ和面阵探测器Ⅱ、控制及图像处理系统；所述射线源、面阵探测器的信号线路与控制及图像处理系统相连，圆形轨道固定设置在管道上，使得圆形轨道的中心与管道的某一横截面的圆心重合；射线源、面阵探测器Ⅰ和面阵探测器Ⅱ分别固定在圆形轨道上；所述面阵探测器Ⅰ和面阵探测器Ⅱ在管道的两侧对称地偏置放置，使得射线源产生的锥形射线束能够覆盖管道的外壁区域。本发明可用较低能量的射线及较小尺寸的探测器扫描重建较大直径的管道外部环形区域的图像，射线束张角较小，射线一致性好，可用于重建管道(含固定管道)管壁的外部环形区域的三维图像。
【专利说明】
基于双周外部CT检测装置在线检测大口径管壁的方法及系统
技术领域
[000。 本发明属于管道管壁缺陷检测技术领域，设及一种基于双周外部CT(Computed Tomogra地y，计算机断层成像)检测装置在线检测大口径管壁的方法及系统。
【背景技术】
[0002] 管道运输已经广泛应用于石油、化工、冶金、制药、能源和环保等行业，已经成为现 代工业和国民经济的命脉。管道运输担负着高溫、高压、易燃、易爆和有毒等介质的输送任 务，一旦因缺陷的存在而引发泄漏或爆炸，有可能导致灾难性的事故。在管道的生产，加工 和使用的过程中，不可避免地会存在腐蚀、裂纹、气孔等缺陷，因此采用无损检测技术对管 道实施周期性的主动检测，对一些严重缺陷进行及时维修就可W大量避免事故发生，同时 也能大大延长管道寿命，经济效益十分可观。
[0003] 现有技术中，CT能够无损、准确地检测物体的内部结构、材料和缺陷等等，是最常 用的无损检测技术之一。由于射线束剂量和探测器尺寸的限制，传统的CT成像系统因无法 穿透物体或者无法完全覆盖物体而只能得到穿过物体外部环形区域的投影数据。由于投影 数据的缺失，传统的重建算法无法得到高质量的重建结果。外部CT是处理运种情况的一种 有效手段。传统的外部CT，通常将探测器在物体的两侧对称放置，只扫描并重建感兴趣的物 体外部环形区域。运种扫描方式可在射线源和探测器围绕被扫描物体旋转(或物体旋转）的 一个圆周内完成扫描，扫描速度和普通CT相当。但运种扫描方式要求射线束的张角要足够 大、或射线源到旋转中屯、的距离足够大、或者被检测物体直径相对较小，才能保证物体两侧 的探测器都能采集到射线投影数据。限于在役大口径管壁的检测条件，射线源到旋转中屯、 的距离不够大而且被检测物体直径相对较大，如果射线束的张角较大则边缘射线强度与中 屯、射线强度差别较大，运些都难W满足传统外部CT的应用条件；如只将探测器在物体的一 侧放置，射线束的张角较小，但采集到射线投影数据比传统外部CT少一半，将可能影响物体 外部环形区域的图像重建质量。
[0004] 因此，需要一种外部CT检测装置，能够在线检测大口径管壁处的缺陷，扫描过程易 于机械实现，并且射线束张角较小，射线一致性好，能得到高质量的Ξ维重建图像。
【发明内容】
[0005] 有鉴于此，本发明的目的在于提供一种基于双周外部CT检测装置在线检测大口径 管壁的方法及系统，该方法和系统易于机械实现，射线束张角较小，射线一致性好，并且能 得到高质量的Ξ维重建图像。
[0006] 为达到上述目的，本发明提供如下技术方案：
[0007] -种基于双周外部CT检测装置在线检测大口径管壁的系统，该系统包括射线源 (1)、面阵探测器I(2a)和面阵探测器Π (2b)、控制及图像处理系统(5);
[0008] 所述射线源（1)、面阵探测器I(2a)和面阵探测器Π (2b)的信号线路与控制及图像 处理系统(5)相连，圆形轨道(3)固定设置在管道(4)上，使得圆形轨道(3)的中屯、与管道(4) 的某一横截面的圆屯、重合;射线源（1)、面阵探测器I(2a)和面阵探测器π (2b)分别固定在 圆形轨道上；所述面阵探测器I(2a)和面阵探测器Π (2b)在管道(4)的两侧对称地偏置放 置，使得射线源(1)产生的锥形射线束能够覆盖管道(4)的外壁区域。
[0009] 本发明还提供了一种基于双周外部CT检测装置在线检测大口径管壁的方法，在该 方法中采用了射线源（1)、面阵探测器I(2a)和面阵探测器Π (2b)、控制及图像处理系统 (5)，具体包括W下步骤：
[0010] S1:在控制与图像处理系统(5)的控制下，首先将射线源（1)对准面阵探测器I(2a) 的中屯、，射线源（1)和面阵探测器I(2a)沿圆形轨道(3)旋转一周，得到一部分投影数据，并 传送到控制与图像处理系统(5)中存储，然后将射线源（1)对准面阵探测器Π (2b)的中屯、， 射线源（1)和面阵探测器Π (2b)沿圆形轨道(3)再旋转一周，得到另一部分投影数据，并传 送到控制与图像处理系统巧）中存储；因为射线源（1)发出的锥形射线束不能完全覆盖大口 径管壁，偏置后的面阵探测器I(2a)和面阵探测器Π (2b)只能得到管道外部环形区域的投 影数据，双周扫描得到的投影数据是射线束完全覆盖管道时得到的投影数据的一部分。
[0011] S2:在控制与图像处理系统(5)中将存储的两部分投影数据按照旋转分度进行组 合，得到外部的投影数据；
[0012] S3:根据组合得到的投影数据重建管道外部环形区域的Ξ维图像；
[0013] S4:显示Ξ维重建结果及二维切片图像。
[0014] 进一步，所述步骤S3具体包括W下步骤：
[0015] S31:对组合得到的投影数据进行加权；
[0016] S32:结合斜坡滤波器对加权后的投影数据进行一维滤波；
[0017] S33:对滤波结果进行Ξ维加权反投影。
[0018] 进一步，在步骤S31中，面阵探测器实际采集到的组合后的投影数据为p(i3，a，b)， 则投影数据加权后为：
[0019]
[0020] 其中，β表示中屯、射线与y轴的夹角，a，b分别代表虚拟探测器上的水平位置和竖直 位置坐标，R表示射线源到旋转中屯、的距离。
[0021] 进一步，在步骤S32中，对加权后的投影数据p/(i3，a，b)进行一维滤波：
[0022]
[002引其中V'代表卷积算子
'q为积分变量。
[0024] 进一步，在步骤S33中，将卷积后的投影数据另化。，6)用Ξ维坐标表示为 P( A "托3'，风斬.V，Z, /?))，则待重建点(X，y，Z )处的重建公式为：
[0025]
[0026] 其中（x，y，z)表示待重建点的Ξ维坐标，f(x，y，z)表示在(x，y，z)处的图像灰度，U (x,y ,β) =R-xsin0+ycos0
表示对应射线在虚拟探测器上的水平位 置坐标
表示对应射线在虚拟探测器上的竖直位置坐标。
[0027] 本发明的有益效果在于:本发明利用两块小尺寸的面阵探测器和现有的CT机的机 构实现大口径管壁外部环形区域内裂纹和缺陷的检测，扫描前将射线源和两块面阵探测器 设置在围绕管道的固定圆形轨道上，并使两块探测器对称地偏置放置，扫描开始后，射线源 对准一侧的面阵探测器中屯、，绕待检测的大口径管壁旋转一周，得到一部分投影数据，然后 射线源对准另一侧的面阵探测器中屯、，绕待检测的大口径管壁再旋转一周，得到另一部分 投影数据，然后将两部分投影数据按照旋转分度进行组合，得到重建算法所需的外部投影 数据，扫描过程易于机械实现。然后利用FDK算法对管道进行Ξ维图像重建。本发明可用较 低能量的射线及较小尺寸的探测器扫描重建较大直径的管道外部环形区域的图像，射线束 张角较小，射线一致性好，可用于重建管道(含固定管道)管壁的外部环形区域的Ξ维图像。
【附图说明】
[0028] 为了使本发明的目的、技术方案和有益效果更加清楚，本发明提供如下附图进行 说明：
[0029] 图1为本发明的管道检测结构示意图；
[0030] 图2为本发明的待检测管道扫描时的横截面示意图；
[0031] 图3为本发明的圆周轨迹锥束外部CT的抑K重建算法几何结构示意图；
[0032] 图4为本发明的双周外部CT检测装置在线检测大口径管壁的方法流程图。
【具体实施方式】
[0033] 下面将结合附图，对本发明的优选实施例进行详细的描述。
[0034] 图1为本发明的管道检测结构示意图，图2为本发明的待检测管道扫描时的横截面 示意图，图3为本发明的圆周轨迹锥束外部CT的FDK重建算法几何结构示意图，如图所示:将 射线源（1)，两个面阵探测器:面阵探测器I(2a)和面阵探测器Π (2b)设置在绕管道(4)的圆 形轨道(3)上，W起始射线源（υ到管道(4)中屯、轴的垂足为坐标原点0建立空间直角坐标系 〇-xyz，x轴为原点与射线源（1)的连线并且正方向为从原点指向射线源（l)，y轴为沿管道 (4)横向并垂直于X轴的坐标轴（图2)，z轴为与管道(4)中屯、轴线重合的坐标轴并且W管道 一侧作为正方向（图1)。W坐标原点0为旋转中屯、，（X，y，Z)表示被重建点坐标，在过坐标原 点0并且垂直于中屯、射线的位置处引入虚拟探测器，射线源（1)位于s，s〇/表示锥形射线的 中屯、射线，SK表示经过被重建点的一条射线，r表示K点在探测器中屯、层(z = 0)的投影，被 重建点在直线SK/上的投影为M，中屯、射线与y轴成β角，K为射线SK的锥角（图3)。
[0035] 本发明公开的一种双周外部CT检测装置在线检测大口径管壁的方法的主要步骤 包括:将射线源和两个探测器设置在围绕待检测管道中屯、的圆形轨道上，探测器相对于管 道对称地偏置放置，扫描开始后，射线源分别对准两个面阵探测器保持相对位置并沿围绕 管道的轨道作两周圆周运动（图1中实线部分的射线和探测器表示第一周扫描，虚线部分的 射线和探测器表示第二周扫描），获得两部分管道外部环形区域的投影数据，然后按照旋转 分度进行组合，得到算法所需的外部投影数据。此类扫描方式的机械运动实现简单易行，能 够利用小尺寸的面板探测器在线检测大口径的管道壁处的缺陷，射线束张角较小，射线一 致性好。
[0036] 图4为本发明的双周外部CT检测装置在线检测大口径管壁的方法流程图，如图所 示:双周外部CT检测装置在线检测大口径管壁的方法包括W下步骤：
[0037] S1.检测装置安装:包括射线源（1)、面阵探测器I(2a)和面阵探测器Π (2b) W及控 制及图像处理系统(5)，所述射线源（1)、面阵探测器I(2a)和面阵探测器Π (2b)的信号线路 与控制及图像处理系统（5)相连，圆形轨道(3)固定设置在管道(4)上，使得圆形轨道(3)的 中屯、与管道(4)某一横截面的圆屯、重合，射线源（1)、面阵探测器I(2a)和面阵探测器Π (2b) 分别固定在圆形轨道上，面阵探测器I(2a)和面阵探测器Π (2b)在管道(4)的两侧对称地偏 置放置，使得射线源(1)产生的锥形射线束能够覆盖管道(4)的外壁区域；
[0038] S2.扫描:在控制与图像处理系统(5)的控制下，首先将射线源（1)对准面阵探测器 I(2a)的中屯、，射线源（1)和面阵探测器I(2a)沿圆形轨道(3)旋转一周，得到一部分投影数 据，并传送到控制与图像处理系统（5)中存储，然后将射线源（1)对准面阵探测器Π (2b)的 中屯、，射线源（1)和面阵探测器Π (2b)沿圆形轨道(3)旋转一周，得到另一部分投影数据，并 传送到控制与图像处理系统(5)中存储；
[0039] S3.投影数据组合:在控制与图像处理系统(5)中将存储的两部分投影数据按照旋 转分度进行组合，得到外部的投影数据；
[0040] S4.管道外部环形区域Ξ维图像的重建:根据组合得到的投影数据重建管道外部 环形区域的Ξ维图像，主要包含Ξ个步骤:S41.对组合得到的投影数据进行加权;S42.结合 斜坡滤波器对加权后的投影数据进行一维滤波；S43.对步骤S42的滤波结果进行Ξ维加权 反投影。
[0041 ] 在本实施例中，假设管道感兴趣区域(Region of Interest,R0I)为：
[0042]
[0043] 其中rminJmax分别为管道外部环形区域的内径和外径，Zmin，Zmax分别表示所要检测 管道Z坐标的初始值与结束值。在感兴趣区域外，重建图像的值为零。
[0044] 重建算法如下：
[0045] 1)面阵探测器实际采集到的组合后的投影数据为p(i3，a，b)，则投影数据加权后 为：
[0046]
[0047] 其中，β表示中屯、射线与y轴的夹角，a，b分别代表虚拟探测器上的水平位置和竖直 位置坐标，R表示射线源到旋转中屯、的距离。
[004引2)对加权后的投影数据p/ W，a，b)进行一维滤波
[0049]
[0化0]其中V'代表卷积算子
^斜坡滤波函数，q为积分变量。
[0051 ] 3)将卷积后的投影数据パ(/?，。，W用Ξ维坐标表示为：パ'ク，。(W'，ク)，W.γ，,v，z.ク))，则待重 建点(x，y，z)处的重建公式为：
[0化2]
[0053]其中（x，y，z)表示待重建点的Ξ维坐标，f(x，y，z)表示在(x，y，z)处的图像灰度，U (x,y ,β) =R-xsin0+ycos0，
表示对应射线在虚拟探测器上的水平位 置坐标
表示对应射线在虚拟探测器上的竖直位置坐标。
[0化4] S5.显示Ξ维重建结果及二维切片图像。
[0055]最后说明的是，W上优选实施例仅用W说明本发明的技术方案而非限制，尽管通 过上述优选实施例已经对本发明进行了详细的描述，但本领域技术人员应当理解，可W在 形式上和细节上对其作出各种各样的改变，而不偏离本发明权利要求书所限定的范围。
【主权项】
1. 一种基于双周外部CT检测装置在线检测大口径管壁的系统，其特征在于:该系统包 括射线源（1)、面阵探测器I(2a)和面阵探测器Π (2b)、控制及图像处理系统(5); 所述射线源（1)、面阵探测器I(2a)和面阵探测器Π (2b)的信号线路与控制及图像处理 系统(5)相连，圆形轨道(3)固定设置在管道(4)上，使得圆形轨道(3)的中心与管道(4)的某 一横截面的圆心重合;射线源（1)、面阵探测器I(2a)和面阵探测器Π (2b)分别固定在圆形 轨道上;所述面阵探测器I(2a)和面阵探测器Π (2b)在管道(4)的两侧对称地偏置放置，使 得射线源(1)产生的锥形射线束能够覆盖管道(4)的外壁区域。2. -种利用如权利要求1所述系统进行在线检测大口径管壁的方法，其特征在于:在该 方法中采用了射线源（1)、面阵探测器I(2a)和面阵探测器Π (2b)、控制及图像处理系统 (5)，具体包括以下步骤： S1:在控制与图像处理系统(5)的控制下，首先将射线源（1)对准面阵探测器I(2a)的中 心，射线源（1)和面阵探测器I(2a)沿圆形轨道(3)旋转一周，得到一部分投影数据，并传送 到控制与图像处理系统(5)中存储，然后将射线源（1)对准面阵探测器Π (2b)的中心，射线 源（1)和面阵探测器Π (2b)沿圆形轨道(3)再旋转一周，得到另一部分投影数据，并传送到 控制与图像处理系统(5)中存储； S2:在控制与图像处理系统(5)中将存储的两部分投影数据按照旋转分度进行组合，得 到外部的投影数据； S3:根据组合得到的投影数据重建管道外部环形区域的三维图像； S4:显示三维重建结果及二维切片图像。3. 根据权利要求2所述的在线检测大口径管壁的方法，其特征在于:所述步骤S3具体包 括以下步骤： S31:对组合得到的投影数据进行加权； S32:结合斜坡滤波器对加权后的投影数据进行一维滤波； S33:对滤波结果进行三维加权反投影。4. 根据权利要求3所述的在线检测大口径管壁的方法，其特征在于:在步骤S31中，面阵 探测器实际采集到的组合后的投影数据为P(i3，a，b)，则投影数据加权后为：其中，β表示中心射线与y轴的夹角，a，b分别代表虚拟探测器上的水平位置和竖直位置 坐标，R表示射线源到旋转中心的距离。5. 根据权利要求4所述的在线检测大口径管壁的方法，其特征在于:在步骤S32中，对加 权后的投影数据P' (i3，a，b)进行一维滤波： p(P，a,b)=-^(p'{P，a,b)*g(a)) 其中代表卷积算子，咖)=?IW如，q为积分变量。 -QD6. 根据权利要求5所述的在线检测大口径管壁的方法，其特征在于:在步骤S33中，将卷 积后的投影数据；^(Αβ,?)用三维坐标表示为P(/?，七%]', /?)』(W，/?))，则待重建点（X，y，z)处的 重建公式为：其中（x，y，z)表示待重建点的三维坐标，f (x，y，z)表示在(x，y，z)处的图像灰度，U(x， y，β) = R-x表示对应射线在虚拟探测器上的水平位置 坐标：丨表示对应射线在虚拟探测器上的竖直位置坐标。
【文档编号】G01N23/04GK106066335SQ201610353562
【公开日】2016年11月2日
【申请日】2016年5月25日 公开号201610353562.3, CN 106066335 A, CN 106066335A, CN 201610353562, CN-A-106066335, CN106066335 A, CN106066335A, CN201610353562, CN201610353562.3
【发明人】曾理, 郭雨濛
【申请人】重庆大学