一种适用于锥束ＸＣＴ系统的探测器扭转角的标定方法技术方案

本发明专利技术公开了一种适用于锥束ＸＣＴ系统的探测器扭转角的标定方法，该标定采用了球状目标体在锥束射线场中绕ｚ轴旋转一周，其在成像坐标系ｘｄｏｄｚｄ中的投影质心轨迹为一近似椭圆，利用最小二乘拟合法将投影质心的轨迹拟合为一椭圆方程ｘｄ２＋ａｚｄ２＋ｂｘｄｚｄ＋ｃｘｄ＋ｅｚｄ＋ｆ＝０，椭圆的长轴与ｘｄ轴的夹角即为探测器扭转角ξ。本发明专利技术利用椭圆的长轴与ｘｄ轴的夹角来表征探测器的扭转角ξ，使得标定过程中计算量少，不易出现假解，重复精度高。

【技术实现步骤摘要】

本专利技术涉及一种锥束XCT(X-ray Computed Tomography)系统，更特别地说，是指一种锥束XCT系统的探测器扭转角标定方法。
技术介绍
近些年来，随着计算机技术的飞速发展和面阵探测器的出现，锥束XCT(Cone-beam X-ray Computed Tomography)日益成为NDT(Non-destructive Testing)领域内的研究热点。在众多的CT重建算法中，考虑到运算量和工程实现难度，FDK(Feldkamp-Davis-Kress)类型的算法最为实用，也一直是实际工程应用中的主流。锥束XCT系统的扫描原理如图1所示，即射线源1发出的锥束射线2对多自由度载物台3上的被测物体4进行透照，被测物体4在多自由度载物台3的带动下绕轴线旋转，面阵探测器5采集被测物体4在不同视角下的DR投影(DR-Digital Radiography，射线数字成像)，最后PC机中的图像重构单元利用这些二维DR投影序列进行三维重建。如图1A所示，FDK重建算法是在射线源和探测器构造的坐标系xyz(也称为重建坐标系)中进行，探测器5的成像坐标系为xdodzd，FDK重建算法理论上要求这两个坐标系的关系为：odzd轴平行于oz轴，odxd轴平行于ox轴。然而在实际的锥束XCT系统安装中，不可避免地存在着机械安装误差，导致odzd轴不平行于oz轴，odxd轴不平行于ox轴，相当于坐标系xdodzd绕y轴旋转了一定的角度，该角度即为探测器的扭转角ξ。扭转角ξ的存在影响了二维DR投影序列重建图像的精度，造成伪影的产生，从而影响重建图像的分辨力、以及重建图像细节的有效检出。
技术实现思路
本专利技术的目的是提出一种对锥束XCT扫描系统中探测器扭转角ξ的标定方法，该标定采用了球状目标体在锥束射线场中绕z轴旋转一周后，其在成像坐标系xdodzd中的DR投影质心轨迹为一近似椭圆，并利用最小二乘拟合法将DR投影质心的轨迹拟合为xd2+azd2+bxdzd+cxd+ezd+f＝0的椭圆方程；在DR投影质心椭圆轨迹中，其椭圆长轴与xd轴的夹角即为探测器扭转角ξ；在本专利技术中为了实现球状目标体绕z轴旋转，设计了偏心滑动支撑架，并将球状成像目标体(被测物体)安装在该偏心滑动支撑架上，而偏心滑动支撑架安装在多自由度载物台上；当多自由度载物台绕轴线旋转360度时，探测器将采集到球状目标体在不同视角下的DR投影序列。本专利技术的一种适用于锥束XCT系统的探测器扭转角的标定方法，具体步骤如下：第一步：调整球状目标体成像在探测器中的位置将多自由度载物台3置于射线源1与面阵探测器5之间的任意位置，偏心滑动支撑架的纵向套筒6插入多自由度载物台3的中心孔里，将一个球状目标体10固定在滑块9的立柱上，通过顶紧螺钉61调节升降杆7在轴线上的高度，并移动滑块9在横向导杆8上的位置；调整球状目标体成像在探测器中的位置是保证球状目标体10被锥束射线2照射，-->以及球状目标体成像能够被探测器5成像面所采集；第二步：调整旋转角度启动多自由度载物台3，并使多自由度载物台3绕轴线在360度范围内旋转，在每间隔5度～15度的旋转角度下，面阵探测器5采集得到球状目标体10的DR投影图像；第三步：获取二值图像在图像处理与可视化单元中对第二步得到的每幅DR投影图像进行阈值分割处理，从而得到二值DR图像；在所述二值DR图像中，球状目标体的投影值记为1，其余区域的投影值记为0；第四步：求DR投影质心坐标(xd-i，zd-i)在图像重构单元中对第三步中得到的每幅二值DR图像进行DR投影质心坐标求取，即DR投影质心坐标(xd-i，zd-i)中M表示DR投影的长度；N表示DR投影的高度；f(xd，zd)表示二值DR图像的二维函数，其中，xd表示成像坐标系xdodzd下xd轴上的坐标变量，zd表示成像坐标系xdodzd下zd轴上的坐标变量；xd-i表示第i幅二值DR图像中球状目标体投影质心在xd轴上的坐标；zd-i表示第i幅二值DR图像中球状目标体投影质心在zd轴上的坐标；第五步：拟合DR投影质心轨迹在图像重构单元中对第四步中得到的每一质心坐标利用最小二乘拟合法进行轨迹拟合，即将所有DR投影质心点坐标(xd-i，zd-i)拟合为xd2+azd2+bxdzd+cxd+ezd+f＝0的椭圆方程，依据该椭圆方程能够得到所有DR投影质心点坐标(xd-i，zd-i)形成的椭圆轨迹，所述椭圆轨迹的斜率为依据该椭圆轨迹斜率从而得到面阵探测器在安装时产生的扭转角ξ＝arctg(k)；xd表示成像坐标系xdodzd下xd轴上的坐标变量；zd表示成像坐标系xdodzd下zd轴上的坐标变量；a表示椭圆方程中变量zd的二次项系数；b表示椭圆方程中变量xdzd的二次项系数；c表示椭圆方程中变量xd的一次项系数；e表示椭圆方程中变量zd的一次项系数；f表示椭圆方程的常数项。本专利技术标定方法的优点：1)通过将球状目标体安装在偏心滑动支撑架上，多自由度载物台带动偏心滑动支撑架绕z轴旋转一周，球状目标体在成像坐标系xdodzd中形成近似椭圆的投影轨迹。通过球状目标体在滑动支撑架上的移动，可调节近似椭圆轨迹的形状，有利于得到最佳的椭圆轨迹。-->2)利用最小二乘法对近似椭圆轨迹进行拟合得到的椭圆方程xd2+azd2+bxdzd+cxd+ezd+f＝0，能够满足误差最小原则，从而保证了标定结果的精度。3)利用椭圆的长轴与xd轴的夹角来表征探测器的扭转角ξ，使得标定过程中计算量少，不易出现假解，重复精度高。4)本专利技术标定方法容易实现，原理简单，只需要将球状目标体在射线源与探测器之间旋转360度即可。附图说明图1是锥束XCT系统的扫描原理图。图1A是安装锥束XCT系统时出现的探测器扭转角机械安装误差示意图。图2是本专利技术的偏心滑动支撑架的结构图。图2A是本专利技术偏心滑动支撑架中横向导杆的结构图。图2B是本专利技术偏心滑动支撑架中滑块的结构图。图3是将本专利技术的偏心滑动支撑架置于锥束XCT系统中的位置摆放示意图。图4是将球状目标体的24幅DR图像合成为一幅DR图像的结果图。图4A是拟合得到的椭圆示意图。图4B是图4A的局部放大图。图5是未标定前的CT重建图像。图5A是采用本专利技术方法标定后的CT重建图像。具体实施方式下面将结合附图和实施例对本专利技术做进一步的详细说明。本专利技术提出的一种适用于锥束XCT系统的探测器扭转角的标定方法，是出厂前对锥束XCT系统的探测器扭转角ξ进行标定，或者是锥束XCT系统经一段时间使用后，对探测器扭转角ξ进行校正时而采用的一种较为简便、易操作的标定方法。一套锥束XCT系统一般由硬件部分和软件部分组成，其中，硬件部分包括：射线源、多自由度载物台、控制器、PC机、面阵探测器；软件部分包括：CT控制单元、图像重构单元、图像处理与可视化单元。参见图3所示，本专利技术对面阵探测器5的成像坐标系xdodzd绕y轴产生一定的旋转角度ξ进行标定，采用在轴线方向上能够升降、旋转方向上能够随多自由度载物台3转动的偏心滑动支撑架(如图2、图2A、图2B所示)来安装球状成像目标体10(被测物体)，使得球状成像目标体10在锥束射线2的照射下所成的投影被面阵探测器5采集。在本专利技术中，首先设计偏心滑动支撑架如图2、图2A、图2B所示，本文档来自技高网...

【技术保护点】
一种适用于锥束ＸＣＴ系统的探测器扭转角的标定方法，其特征在于包括有下列标定步骤：第一步：调整球状目标体成像在探测器中的位置将多自由度载物台（３）置于射线源（１）与面阵探测器（５）之间的任意位置，偏心滑动支撑架的纵向套筒（６）插入多自由度载物台（３）的中心孔里，将一个球状目标体（１０）固定在滑块（９）的立柱上，通过顶紧螺钉（６１）调节升降杆（７）在轴线上的高度，并移动滑块（９）在横向导杆（８）上的位置；调整球状目标体成像在探测器中的位置是保证球状目标体（１０）被锥束射线（２）照ｄ］ｚ↓［ｄ］的二次项系数；ｃ表示椭圆方程中变量ｘ↓［ｄ］的一次项系数；ｅ表示椭圆方程中变量ｚ↓［ｄ］的一次项系数；ｆ表示椭圆方程的常数项。射，以及球状目标体成像能够被探测器（５）成像面所采集；第二步：调整旋转角度启动多自由度载物台（３），并使多自由度载物台（３）绕轴线在３６０度范围内旋转，在每间隔５度～１５度的旋转角度下，面阵探测器（５）采集得到球状目标体（１０）的ＤＲ投影图像；第三步：获取二值图像在图像处理与可视化单元中对第二步得到的每幅ＤＲ投影图像进行阈值分割处理，从而得到二值ＤＲ图像；在所述二值ＤＲ图像中，球状目标体的投影值记为１，其余区域的投影值记为０；第四步：求ＤＲ投影质心坐标（ｘ↓［ｄ－ｉ］，ｚ↓［ｄ－ｉ］）在图像重构单元中对第三步中得到的每幅二值ＤＲ图像进行ＤＲ投影质心坐标求取，即ＤＲ投影质心坐标（ｘ↓［ｄ－ｉ］，ｚ↓［ｄ－ｉ］）中ｘ↓［ｄ－ｉ］＝１／ＭＮ＊＊ｘ↓［ｄ］×ｆ（ｘ↓［ｄ］，ｚ↓［ｄ］），ｚ↓［ｄ－ｉ］＝１／ＭＮ＊＊ｚ↓［ｄ］×ｆ（ｘ↓［ｄ］，ｚ↓［ｄ］）；Ｍ表示ＤＲ投影的长度；Ｎ表示ＤＲ投影的高度；ｆ（ｘ↓［ｄ］，ｚ↓［ｄ］）表示二值ＤＲ图像的二维函数，其中，ｘ↓［ｄ］表示成像坐标系ｘ↓［ｄ］ｏ↓［ｄ］ｚ↓［ｄ］下ｘ↓［ｄ］轴上的坐标变量，ｚ↓［ｄ］表示成像坐标系ｘ↓［ｄ］ｏ↓［ｄ］ｚ↓［ｄ］下ｚ↓［ｄ］轴上的坐标变量；ｘ↓［ｄ－ｉ］表示第ｉ幅二值ＤＲ图像中球状目标体投影质心在ｘ↓［ｄ］轴上的坐标；ｚ↓［ｄ－ｉ］表示第ｉ幅二值ＤＲ图像中球状目标体投影质心在ｚ↓［ｄ］轴上的坐标；第五步：拟合ＤＲ投影质心轨迹在图像重构单元中对第四步中得到的每一质心坐标利用最小二乘拟合法进行轨迹拟合，即将所有ＤＲ投影质心点坐标（ｘ↓［ｄ－ｉ］，ｚ↓［ｄ－ｉ］）拟合为ｘ↓［ｄ］↑［２］＋ａｚ↓［ｄ］↑［２］＋ｂｘ↓［ｄ］ｚ↓［ｄ］＋ｃ...

【技术特征摘要】
1.一种适用于锥束XCT系统的探测器扭转角的标定方法，其特征在于包括有下列标定步骤：第一步：调整球状目标体成像在探测器中的位置将多自由度载物台(3)置于射线源(1)与面阵探测器(5)之间的任意位置，偏心滑动支撑架的纵向套筒(6)插入多自由度载物台(3)的中心孔里，将一个球状目标体(10)固定在滑块(9)的立柱上，通过顶紧螺钉(61)调节升降杆(7)在轴线上的高度，并移动滑块(9)在横向导杆(8)上的位置；调整球状目标体成像在探测器中的位置是保证球状目标体(10)被锥束射线(2)照射，以及球状目标体成像能够被探测器(5)成像面所采集；第二步：调整旋转角度启动多自由度载物台(3)，并使多自由度载物台(3)绕轴线在360度范围内旋转，在每间隔5度～15度的旋转角度下，面阵探测器(5)采集得到球状目标体(10)的DR投影图像；第三步：获取二值图像在图像处理与可视化单元中对第二步得到的每幅DR投影图像进行阈值分割处理，从而得到二值DR图像；在所述二值DR图像中，球状目标体的投影值记为1，其余区域的投影值记为0；第四步：求DR投影质心坐标(xd-i，zd-i)在图像重构单元中对第三步中得到的每幅二值DR图像进行DR投影质心坐标求取，即DR投影质心坐标(xd-i，zd-i)中M表示DR投影的长度；N表示DR投影的高度；f(xd，zd)表示二值DR图像的二维函数，其中，xd表示成像坐标系xdodzd下xd轴上的坐标变量，zd表示成像坐标系xdod...

【专利技术属性】
技术研发人员：杨民，刘永瞻，高海东，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：11[中国|北京]