CN101832954A - 锥束xct系统用移动组件以及用其进行图像重建坐标系原点标定的方法 - Google Patents
锥束xct系统用移动组件以及用其进行图像重建坐标系原点标定的方法 Download PDFInfo
- Publication number
- CN101832954A CN101832954A CN 201010123892 CN201010123892A CN101832954A CN 101832954 A CN101832954 A CN 101832954A CN 201010123892 CN201010123892 CN 201010123892 CN 201010123892 A CN201010123892 A CN 201010123892A CN 101832954 A CN101832954 A CN 101832954A
- Authority
- CN
- China
- Prior art keywords
- guide rail
- axis
- array detector
- cone
- imaging plane
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种适用于锥束XCT系统的移动组件以及用其进行图像重建坐标系原点标定方法,其主要针对基于FDK算法的圆轨迹扫描成像系统投影坐标原点的标定。由于射线源焦点和探测器成像平面的准确空间位置无法直接测量得到,从而使得投影坐标原点难以精确测量。本发明提出利用单圆球目标体在锥束场不同位置多次成像获得的DR图像序列,基于图像、图形处理方法和最小二乘拟合技术,重构超定方程组,然后求解出射线源焦点在成像平面上的投影坐标点,利用此方法得到的参数值进行三维重建,得到了精确的重建结果。
Description
技术领域
本发明涉及一种锥束XCT(X-ray Computed Tomography)系统,更特别地说,是指一种锥束XCT系统的图像重建坐标系原点标定方法。
背景技术
近些年来,随着计算机技术的飞速发展和面阵探测器的出现,锥束XCT(Cone-beam X-ray Computed Tomography)日益成为NDT(Non-destructiveTesting)领域内的研究热点。在众多的CT重建算法中,考虑到运算量和工程实现难度,FDK(Feldkamp-Davis-Kress)类型的算法最为实用,也一直是实际工程应用中的主流。
根据FDK算法要求,需要精确知道图像重建坐标系原点的位置,即射线源的焦点在面阵探测器的成像平面上的投影点坐标。一般射线源与面阵探测器的成像平面之间的距离为1.5m~2m。但在实际的成像系统中,由于射线源的焦点和面阵探测器的成像平面的准确位置无法直接测量得到,从而使得焦点在成像平面上的投影点的位置难以精确测量。因此,投影点的误差会影响重建图像的精度,造成伪影的产生,影响图像的分辨力和细节的有效检出。
在“利用非线性最小二乘估计锥束扫描三维CT重建几何参数,北京航空航天大学学报,2005,31(10):1135-1139”中提出了一种非线性最小二乘估计的测量方法。其思想就是计算空间一质点在不同转角下的投影坐标,质点理论投影坐标与实际求得的投影坐标之间应满足误差最小,通过求解满足此条件的非线性最小二乘解,即得到投影坐标原点参数值。非线性最小二乘估计结果受参数向量的初始值影响较大,容易产生病态的解。
作者张全红在“X射线工业CT成像优化研究[D].北京:北京航空航天大学机械工程及自动化学院,2006”中提出了利用夹芯面包板测量图像重建坐标系原点的方法,夹芯面包板由长准直比的准直栅构成,当对面包板进行成像时,惟有焦点发出的主射线通过准直栅成像,通过对准直栅DR(Digital Radiography)图像的灰度分布统计和拟合,计算出射线源焦点投影的位置。该方法测量精度较低,操作起来难度较大。
发明内容
本发明的目的是提出一种适用于锥束XCT系统的图像重建坐标系原点标定方法、以及在标定过程中获取X轴向或Y轴向上位置的移动组件,该方法利用单球状目标体在锥束场不同位置下多次成像获得的DR投影图像序列,并将在同一X轴向位置下的、两个Y轴方向具有位移差的移动位置下的DR投影图像归为一组,然后采用图像处理方法提取出每组两个DR投影图像的中心坐标,连接两中心点构成一直线方程,以此得到的多个直线方程重构出超定方程组,然后求解超定方程组得到图像重建坐标系原点坐标。
本发明的一种适用于锥束XCT系统中提供位置点的移动组件,由Y轴向导轨(6)、X轴向导轨(7)、滑块(71)和支撑架(8)组成;支撑架(8)由升缩杆(81)和套筒(83)组成,升缩杆(81)的一端能够在套筒(83)内上下滑动,升缩杆(81)的另一端上用于放置球状目标体(3),支撑架(8)的底部设有导柱(82),导柱(82)能够在X轴向导轨(7)的滑道(72)中滑动;X轴向导轨(7)的上方设有滑道(72),滑块(71)安装在X轴向导轨(7)的下方;X轴向导轨(7)与Y轴向导轨(6)通过滑块(71)滑动连接;移动组件摆放在射线源(1)与面阵探测器(4)之间的任意位置,且移动组件中的Y轴向导轨(6)垂直于面阵探测器(4)的成像平面(5)。
本发明的一种适用于锥束XCT系统的图像重建坐标系原点标定方法,进行图像重建坐标系原点标定时执行包括下列步骤:
第一步:将移动组件摆放在射线源(1)与面阵探测器(4)之间的任意位置,且移动组件中的Y轴向导轨(6)垂直于面阵探测器(4)的成像平面(5);
然后将一个球状目标体(3)安装在升缩杆(81)的顶部;
第二步:开启射线源设备出射锥束射线(2),调整升缩杆(81)的高度,使锥束射线(2)能够照射球状目标体(3)在面阵探测器(4)的成像平面(5)上并有DR投影图像(31);
第三步:(A)调整导柱(82)在X轴向导轨(7)上的第一位置PX1,在该第一位置PX1下面阵探测器(4)采集到球状目标体(3)在成像平面(5)上的第一位置投影B11,第一位置投影B11的中心坐标记为Q11(x11,z11);
(B)在第一位置PX1下,将X轴向导轨(7)沿Y轴向导轨(6)从左向右移动一段距离Δy到达第一移动位置PY(X1)后,面阵探测器(4)采集到球状目标体(3)在成像平面(5)上的第一移动位置投影B12,第一移动位置投影B12的中心坐标记为Q12(x12,z12);
第四步:(A)再调整导柱(82)在X轴向导轨(7)上的第二位置PX2,在该第二位置PX2下面阵探测器(4)采集到球状目标体(3)在成像平面(5)上的第二位置投影B21,第二位置投影B21的中心坐标记为Q21(x21,z21);
(B)在第二位置PX2下,将X轴向导轨(7)沿Y轴向导轨(6)从左向右移动一段距离Δy到达第二移动位置PX(X2)后,面阵探测器(4)采集到球状目标体(3)在成像平面(5)上的第二移动位置投影B22,第二移动位置投影B22的中心坐标记为Q22(x22,z22);
第五步:依次调整导柱(82)在X轴向导轨(7)上的位置,分别记为位置PX3,……,PXN,然而面阵探测器(4)分别采集到球状目标体(3)在成像平面(5)上的DR投影图像分别记为B31,……,BN1;则B31的中心坐标为Q31(x31,z31),BN1的中心坐标为QN1(xN1,zN1);在位置PX3,……,PXN下,将X轴向导轨(7)沿Y轴向导轨(6)移动一段距离Δy后到达位置PY(X3),……,PY(XN),则面阵探测器(4)分别采集到球状目标体(3)在成像平面(5)上的DR投影图像分别记为B32,……,BN2;则B32的中心坐标为Q32(x32,z32),BN2的中心坐标为QN2(xN2,zN2);
(B)连接中心点Q21和Q22,得到第二直线方程
(C)连接中心点Q31和Q32得到第三直线方程
(D)连接中心点QN1和QN2得到第N直线方程
(E)联立在一次测量过程中所有的直线方程,得到超定方程组 且
附图说明
图1是本发明球状目标体在锥束射线场下的成像示意图。
图1A是本发明设计的支撑架结构。
图1B是本发明设计的二自由度移动组件结构图。
图2A和图2B是球状目标体在X轴向移动前后的成像示意图。
图3A和图3B是球状目标体在Y轴向移动前后的成像示意图。
图4是面阵探测器4采集到的DR投影合成图像。
图5A是160kV锥束XCT系统使用2年后的CT重建图像。
图5B是采用本发明方法重新进行160kV锥束XCT系统图像重建坐标系原点标定后的CT重建图像。
具体实施方式
下面将结合附图和实施例对本发明做进一步的详细说明。
本发明提出的一种适用于锥束XCT系统的图像重建坐标系原点标定方法,是出厂前对锥束XCT系统的图像重建坐标系原点Od(xo,zo)进行标定,或者是锥束XCT系统经一段时间使用后,对图像重建坐标系原点Od(xo,zo)的位置进行校正时而采用的一种较为简便、易操作的标定方法。一套锥束XCT系统一般由硬件部分和软件部分组成,其中,硬件部分包括有射线源、多自由度载物台、探测器、控制器和PC机;软件部分(存储在PC机的硬盘中)包括有CT控制单元、图像重构单元、图像处理与可视化单元。
在本发明中(参见图1、图1A、图1B所示),不采用锥束XCT系统中原有的多自由度载物台,而采用由支撑架8、X轴向导轨7和Y轴向导轨6构成的二自由度移动组件。支撑架8由升缩杆81和套筒83组成,升缩杆81的一端能够在套筒83内上下滑动,升缩杆81的另一端上用于放置球状目标体3,支撑架8的底部设有导柱82,导柱82能够在X轴向导轨7的滑道72中滑动,从而确定球状目标体3在X轴向上的位置;X轴向导轨7的上板面上设有滑道72,该滑道72用于安装支撑架8的导柱82,从而实现支撑架8沿X轴向导轨7的轴向滑动,X轴向导轨7的下方安装有滑块71,X轴向导轨7通过滑块71与Y轴向导轨6安装,且使X轴向导轨7沿Y轴向导轨6的轴向方向滑动;Y轴向导轨6的长度为1000~1200mm,X轴向导轨7的长度为500~800mm。支撑架8的一端设有升缩杆81,另一端设有导柱82,升缩杆81能够在套筒83内滑动,支撑架8的升缩杆81上放置球状目标体3。移动组件位于射线源1与面阵探测器4之间任意位置,并在移动组件的升缩杆81的顶部放置一球状目标体3。当射线源1射出的锥束射线2照射到球状目标体3上时,球状目标体3会在阵探测器4的成像平面5上形成投影图像31(即DR投影图像)。
本发明是对射线源1的焦点P投射在面阵探测器4的成像平面5上的投影点Od(xo,zo)的位置进行标定,采用单球状目标体3在锥束场下,在同一X轴向位置下、通过移动单球状目标体3在射线源1与面阵探测器4之间的位置(即沿着Y轴向导轨6移动),使得成像平面5上形成两个不同Y轴向位置下的DR投影图像;然后采用图像处理方法提取出两个DR投影图像的中心坐标,连接两个中心点构成直线方程重复进行不同X轴向位置下的两个不同Y轴向位置下的DR投影图像获取,以此得到的多个直线方程重构出超定方程组最后在图像处理与可视化单元中依据该超定方程组求解得到源焦点P在成像平面5上的投影坐标点Od(xo,zo),即得到图像重建坐标系原点坐标。该测量的结果应用于图像重构单元和CT控制单元中。
参见图2A、图2B所示的球状目标体3沿X轴向移动成像。将支撑架8的导柱82在X轴向导轨7上由后向前每次移动一段距离Δx(Δx=15mm~20mm)后到达的位置记为PXN,X表示沿X轴向移动,N表示沿X轴向的移动次数,N=5~7。球状目标体3沿X轴向移动的起始位置记为PX1,其他位置依次记为PX2,PX3,……,PXN。
参见图3A、图3B所示的球状目标体3沿Y轴向移动成像。当支撑架8在X轴向导轨7上的位置确定后,X轴向导轨7在Y轴向导轨6上由左向右每次移动一段距离Δy(Δy=100mm~150mm)后到达的位置记为PY(XN),Y表示沿Y轴向移动,XN表示球状目标体3在X轴向上的位置点,N表示沿X轴向的移动次数,N=5~7。
依据图2A、图2B、图3A、图3B所示,采用本发明的投影点位置标定方法,首先是支撑架8在X轴向导轨7上移动到某一位置PXN时,然后沿Y轴向导轨6由左向右移动一段距离Δy后到达一个位置。当球状目标体3在X轴向上的位置为PX1时,再沿Y轴移动一段距离Δy后到达的位置记为PY(X1),依次类推,则PX2对应PY(X2),PX3对应PY(X3),……,PXN对应PY(XN)。
本发明是一种适用于锥束XCT系统的图像重建坐标系原点标定方法,该图像重建坐标系原点标定包括有下列步骤:
第一步:将移动组件摆放在射线源1与面阵探测器4之间的任意位置,且移动组件中的Y轴向导轨6垂直于面阵探测器4的成像平面5,支撑架8的底部安装在X轴向导轨7与Y轴向导轨6的垂直位置,将一个球状目标体3安装在升缩杆81的顶部;如图1所示,图中移动组件是为了直观表达安装关系,采用了筒示的移动组件,对于移动组件的具体结构请参见图1A和图1B所示。
第二步:开启射线源设备出射锥束射线2,调整升缩杆81的高度,使锥束射线2能够照射球状目标体3在面阵探测器4的成像平面5上出现有DR投影图像31;
第三步:(A)调整导柱82在X轴向导轨7上的位置,记为位置PX1(第一位置PX1),在该第一位置PX1下面阵探测器4采集到球状目标体3在成像平面5上的DR投影图像记为B11(第一位置投影B11),第一位置投影B11的中心坐标记为Q11(x11,z11);
(B)在第一位置PX1下,将X轴向导轨7沿Y轴向导轨6从左向右移动一段距离Δy到达位置PY(X1)(第一移动位置PY(X1))后,面阵探测器4采集到球状目标体3在成像平面5上的DR投影图像记为B12(第一移动位置投影B12),第一移动位置投影B12的中心坐标记为Q12(x12,z12);
在本发明中,对DR投影图像进行边缘检测、阈值分割、轮廓细化与轮廓追踪,以得到DR投影图像的轮廓点坐标,然后利用最小二乘拟合法回归出DR投影图像的中心坐标。其中边缘检测、阈值分割、轮廓细化与轮廓追踪、最小二乘拟合求取中心坐标为已公开的图像、图形处理方法。
第四步:(A)再调整导柱82在X轴向导轨7上的位置,记为位置PX2(第二位置PX2),在该第二位置PX2下面阵探测器4采集到球状目标体3在成像平面5上的DR投影图像记为B21(第二位置投影B21),第二位置投影B21的中心坐标记为Q21(x21,z21);
(B)在第二位置PX2下,将X轴向导轨7沿Y轴向导轨6从左向右移动一段距离Δy到达位置PY(X2)(第二移动位置PY(X2))后,面阵探测器4采集到球状目标体3在成像平面5上的DR投影图像记为B22(第二移动位置投影B22),第二移动位置投影B22的中心坐标记为Q22(x22,z22);
第五步:依次调整导柱82在X轴向导轨7上的位置,分别记为位置PX3,……,PXN(N=5~7次),面阵探测器4分别采集到球状目标体3在成像平面5上的DR投影图像分别记为B31,……,BN1;则B31的中心坐标Q31(x31,z31),BN1的中心坐标QN1(xN1,zN1);在位置PX3,……,PXN下,将球状目标体3沿Y轴向导轨6移动一段距离Δy后到达位置PY(X3),……,PY(XN),则面阵探测器4分别采集到球状目标体3在成像平面5上的DR投影图像分别记为B32,……,BN2;则B32的中心坐标Q32(x32,z32),BN2的中心坐标QN2(xN2,zN2)。为了叙述方便,PX3称为第三位置,PY(X3)称为第三移动位置,PXN称为第N位置,PY(XN)称为第N移动位置,B31称为第三位置投影,B32称为第三移动位置投影,BN1称为第N位置投影,BN2称为第N移动位置投影。
(B)连接中心点Q21和Q22,得到第二直线方程
(C)连接中心点Q31和Q32得到第三直线方程
……;
(D)连接中心点QN1和QN2得到第N直线方程
(E)联立在一次测量过程(N=5~7次)中所有的直线方程,得到超定方程组 且
在本发明中,所重构的方程组的个数大于未知数的个数,因而是一超定方程组。该方程组的系数矩阵的行数大于列数,此时,该方程组没有确定意义下的解,只存在最小二乘意义下的最逼近解。本发明采用方程组残差的2-范数取极小值时的最小二乘解为该超定方程组的最逼近解。最小二乘法求解超定方程组为已公开数学解算方法。目前常用的数学计算软件Matlab即可实现。
实施例子
为了验证本发明测量方法的有效性,发明人设计了移动组件,将直径为20mm的钢珠安装在移动组件的升缩杆81上,并将移动组件和钢珠一起放置在160kV锥束XCT系统的射线源与面阵探测器之间,分别进行在X轴上的七个位置点(N=7次,在选取出第一个X轴上的位置后,以后的6个位置点的移动间距为20mm)的七组位置测量。每组采集两个沿Y轴方向不同位置(Y轴上第一位置点与第二位置点相距120mm)的钢珠的DR投影图像,图4为将7组DR投影图像合成为一幅DR投影图像的结果。利用160kV锥束XCT系统中的图像处理技术,提取出每组钢珠DR投影图像的中心坐标,两个中心坐标相连成一直线,那么7组数据就重构出7条直线方程,7条直线方程就重构成一个超定方程组为然后采用最小二乘法求解该方程组得出图像重建坐标系原点坐标为(1025.255,769.700)。
选用的160kV锥束XCT系统的主要硬件配置如下:
(1)射线源:德国YXLON公司产MG165/2.25X-射线源,焦点尺寸0.2mm;
(2)移动组件:由支撑架8、X轴向导轨7和Y轴向导轨6构成;
(3)面阵探测器:美国Varian公司产PaxScan4030CB型非晶硅面阵探测器,成像面积400×300mm2,探元尺寸0.194mm。
该设备图像重建坐标系原点坐标的出厂标定值为(1022,777),在使用2年后,图像重建坐标系原点的实际位置会偏离出厂值。利用该值进行CT重建,结果如图5A所示,图像边缘产生重影。用本发明的测量值进行重建,结果如图5B所示,图像边缘清晰。
实验结果表明,本发明的标定方法具有以下优点:
1)能够适用于普通焦点和微焦点的锥束XCT系统中图像重建坐标系原点的标定。
2)采用一个球状目标体多次成像重构超定方程组,求解超定方程组得到最小二乘解,求解原点坐标值精度达到亚像素级,从而降低了测量中随机误差造成的影响,保证了测量的重复精度。
3)对使用一段时间后的锥束XCT系统,采用本发明的方法进行图像重建坐标系原点标定,球状目标体3在面阵探测器4的成像平面5上的投影点坐标能够恢复至出厂时的位置。
4)测量方法容易实现,操作简单,只需要将球体在射线源与探测器之间沿导轨移动不同距离成像即可。
Claims (10)
1.一种锥束XCT系统用移动组件,其特征在于由Y轴向导轨(6)、X轴向导轨(7)、滑块(71)和支撑架(8)组成;支撑架(8)由升缩杆(81)和套筒(83)组成,升缩杆(81)的一端能够在套筒(83)内上下滑动,升缩杆(81)的另一端上用于放置球状目标体(3),支撑架(8)的底部设有导柱(82),导柱(82)能够在X轴向导轨(7)的滑道(72)中滑动;X轴向导轨(7)的上方设有滑道(72),滑块(71)安装在X轴向导轨(7)的下方;X轴向导轨(7)与Y轴向导轨(6)通过滑块(71)滑动连接;移动组件摆放在射线源(1)与面阵探测器(4)之间的任意位置,且移动组件中的Y轴向导轨(6)垂直于面阵探测器(4)的成像平面(5)。
2.根据权利要求1所述的移动组件,其特征在于:Y轴向导轨(6)的长度为1000mm~1200mm,X轴向导轨(7)的长度为500mm~800mm。
3.应用权利要求1所述的移动组件在锥束XCT系统中进行图像重建坐标系原点标定的方法,其特征在于图像重建坐标系原点标定包括有下列步骤:
第一步:将移动组件摆放在射线源(1)与面阵探测器(4)之间的任意位置,且移动组件中的Y轴向导轨(6)垂直于面阵探测器(4)的成像平面(5);
然后将一个球状目标体(3)安装在升缩杆(81)的顶部;
第二步:开启射线源设备出射锥束射线(2),调整升缩杆(81)的高度,使锥束射线(2)能够照射球状目标体(3)在面阵探测器(4)的成像平面(5)上并有DR投影图像(31);
第三步:(A)调整导柱(82)在X轴向导轨(7)上的第一位置PX1,在该第一位置PX1下面阵探测器(4)采集到球状目标体(3)在成像平面(5)上的第一位置投影B11,第一位置投影B11的中心坐标记为Q11(x11,z11);
(B)在第一位置PX1下,将X轴向导轨(7)沿Y轴向导轨(6)从左向右移动一段距离Δy到达第一移动位置PY(X1)后,面阵探测器(4)采集到球状目标体(3)在成像平面(5)上的第一移动位置投影B12,第一移动位置投影B12的中心坐标记为Q12(x12,z12);
第四步:(A)再调整导柱(82)在X轴向导轨(7)上的第二位置PX2,在该第二位置PX2下面阵探测器(4)采集到球状目标体(3)在成像平面(5)上的第二位置投影B21,第二位置投影B21的中心坐标记为Q21(x21,z21);
(B)在第二位置PX2下,将X轴向导轨(7)沿Y轴向导轨(6)从左向右移动一段距离Δy到达第二移动位置PY(X2)后,面阵探测器(4)采集到球状目标体(3)在成像平面(5)上的第二移动位置投影B22,第二移动位置投影B22的中心坐标记为Q22(x22,z22);
第五步:依次调整导柱(82)在X轴向导轨(7)上的位置,分别记为位置PX3,……,PXN,然而面阵探测器(4)分别采集到球状目标体(3)在成像平面(5)上的DR投影图像分别记为B31,……,BN1;则B31的中心坐标为Q31(x31,z31),BN1的中心坐标为QN1(xN1,zN1);在位置PX3,……,PXN下,将X轴向导轨(7)沿Y轴向导轨(6)移动一段距离Δy后到达位置PY(X3),……,PY(XN),则面阵探测器(4)分别采集到球状目标体(3)在成像平面(5)上的DR投影图像分别记为B32,……,BN2;则B32的中心坐标为Q32(x32,z32),BN2的中心坐标为QN2(xN2,zN2);
(D)连接中心点QN1和QN2得到第N直线方程
(E)联立在一次测量过程中所有的直线方程,得到超定方程组 且
第七步:采用最小二乘法解求解超定方程组得到图像重建坐标系原点坐标点Qd(xO,zO)。
4.根据权利要求3所述的适用于锥束XCT系统的图像重建坐标系原点标定方法,
其特征在于:支撑架(8)的导柱(82)在X轴向导轨(7)上由后向前每次移动距离Δx=15mm~20mm。
5.根据权利要求3所述的适用于锥束XCT系统的图像重建坐标系原点标定方法,其特征在于:当支撑架(8)在X轴向导轨(7)上的位置确定后,X轴向导轨(7)在Y轴向导轨(6)上由左向右每次移动距离Δy=100mm~150mm。
6.根据权利要求3所述的适用于锥束XCT系统的图像重建坐标系原点标定方法,其特征在于:沿X轴向的移动次数N=5~7次。
7.根据权利要求3所述的适用于锥束XCT系统的图像重建坐标系原点标定方法,其特征在于:面阵探测器(4)为平面阵列探测器。
8.根据权利要求3所述的适用于锥束XCT系统的图像重建坐标系原点标定方法,其特征在于:球状目标体(3)的直径为10mm~30mm。
9.根据权利要求3所述的适用于锥束XCT系统的图像重建坐标系原点标定方法,其特征在于:能够适用于普通焦点和微焦点的锥束XCT系统中图像重建坐标系原点的标定,且求解原点坐标值精度达到亚像素级。
10.根据权利要求3所述的适用于锥束XCT系统的图像重建坐标系原点标定方法,其特征在于:对使用一段时间后的锥束XCT系统,经所述的图像重建坐标系原点标定,能够使球状目标体(3)在面阵探测器(4)的成像平面(5)上的投影点坐标能够恢复至出厂时的位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010123892 CN101832954B (zh) | 2010-03-15 | 2010-03-15 | 锥束xct系统中进行图像重建坐标系原点标定的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010123892 CN101832954B (zh) | 2010-03-15 | 2010-03-15 | 锥束xct系统中进行图像重建坐标系原点标定的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101832954A true CN101832954A (zh) | 2010-09-15 |
CN101832954B CN101832954B (zh) | 2012-06-20 |
Family
ID=42717097
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010123892 Expired - Fee Related CN101832954B (zh) | 2010-03-15 | 2010-03-15 | 锥束xct系统中进行图像重建坐标系原点标定的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101832954B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102253061A (zh) * | 2011-04-19 | 2011-11-23 | 东南大学 | 立式锥束计算机断层成像校准系统及应用该系统的方法 |
CN103226113A (zh) * | 2013-03-29 | 2013-07-31 | 中国计量科学研究院 | 锥束3d-ct扫描系统重建体素尺寸的自动标定方法 |
CN103759679A (zh) * | 2013-11-04 | 2014-04-30 | 杭州长庚测量技术有限公司 | 一种锥束ct系统角度偏差测量方法 |
CN104665862A (zh) * | 2015-02-16 | 2015-06-03 | 清华大学 | 在cbct中消除几何伪影的方法以及使用该方法的cbct系统 |
CN107831180A (zh) * | 2016-09-14 | 2018-03-23 | 奚岩 | X射线原位成像方法及系统 |
CN107945132A (zh) * | 2017-11-29 | 2018-04-20 | 深圳安科高技术股份有限公司 | 一种基于神经网络的ct图像的伪影校正方法及装置 |
CN109187591A (zh) * | 2018-06-04 | 2019-01-11 | 东南大学 | 一种x射线超分辨成像方法及其应用 |
CN111175322A (zh) * | 2020-01-09 | 2020-05-19 | 江苏科技大学 | 基于x射线数字成像技术的浆锚搭接连接质量检测方法 |
CN111272782A (zh) * | 2020-03-13 | 2020-06-12 | 彭晟 | 自适应视场的x射线装置及自适应调整视场的方法 |
CN111458739A (zh) * | 2020-03-16 | 2020-07-28 | 北京东软医疗设备有限公司 | 成像方法、装置及系统 |
CN113029053A (zh) * | 2021-04-06 | 2021-06-25 | 中国科学技术大学 | 通用的ct对轴方法 |
CN113137934A (zh) * | 2021-04-23 | 2021-07-20 | 中国工程物理研究院流体物理研究所 | 一种单探头射线方程的标定系统及标定方法 |
CN115060747A (zh) * | 2022-04-13 | 2022-09-16 | 中国兵器科学研究院宁波分院 | 一种工业ct系统用x射线源的焦点尺寸定量方法 |
CN116892893A (zh) * | 2023-09-11 | 2023-10-17 | 上海福柯斯智能科技有限公司 | 工业ct锥束中心投影点量测方法及存储介质 |
CN115060747B (zh) * | 2022-04-13 | 2024-09-24 | 中国兵器科学研究院宁波分院 | 一种工业ct系统用x射线源的焦点尺寸定量方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1392402A (zh) * | 2001-06-20 | 2003-01-22 | 株式会社日立工程服务 | X射线ct装置 |
CN2711727Y (zh) * | 2004-07-13 | 2005-07-20 | 丁莉 | 微型γ射线工业计算机断层扫描成像装置 |
JP2006078251A (ja) * | 2004-09-08 | 2006-03-23 | National Research Institute For Cultural Properties | 木材又は木造文化財の年輪幅又は密度測定装置並びに測定方法 |
CN1861001A (zh) * | 2006-03-31 | 2006-11-15 | 北京航空航天大学 | 一种适用于三维ct扫描系统投影坐标原点的标定方法 |
CN101515370A (zh) * | 2009-03-06 | 2009-08-26 | 北京航空航天大学 | 三维显微ct扫描系统中射线源焦点的投影坐标的标定方法 |
-
2010
- 2010-03-15 CN CN 201010123892 patent/CN101832954B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1392402A (zh) * | 2001-06-20 | 2003-01-22 | 株式会社日立工程服务 | X射线ct装置 |
CN2711727Y (zh) * | 2004-07-13 | 2005-07-20 | 丁莉 | 微型γ射线工业计算机断层扫描成像装置 |
JP2006078251A (ja) * | 2004-09-08 | 2006-03-23 | National Research Institute For Cultural Properties | 木材又は木造文化財の年輪幅又は密度測定装置並びに測定方法 |
CN1861001A (zh) * | 2006-03-31 | 2006-11-15 | 北京航空航天大学 | 一种适用于三维ct扫描系统投影坐标原点的标定方法 |
CN101515370A (zh) * | 2009-03-06 | 2009-08-26 | 北京航空航天大学 | 三维显微ct扫描系统中射线源焦点的投影坐标的标定方法 |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102253061A (zh) * | 2011-04-19 | 2011-11-23 | 东南大学 | 立式锥束计算机断层成像校准系统及应用该系统的方法 |
CN103226113A (zh) * | 2013-03-29 | 2013-07-31 | 中国计量科学研究院 | 锥束3d-ct扫描系统重建体素尺寸的自动标定方法 |
CN103226113B (zh) * | 2013-03-29 | 2015-10-21 | 中国计量科学研究院 | 锥束3d-ct扫描系统重建体素尺寸的自动标定方法 |
CN103759679A (zh) * | 2013-11-04 | 2014-04-30 | 杭州长庚测量技术有限公司 | 一种锥束ct系统角度偏差测量方法 |
CN103759679B (zh) * | 2013-11-04 | 2016-08-17 | 杭州长庚测量技术有限公司 | 一种锥束ct系统角度偏差测量方法 |
CN104665862A (zh) * | 2015-02-16 | 2015-06-03 | 清华大学 | 在cbct中消除几何伪影的方法以及使用该方法的cbct系统 |
CN104665862B (zh) * | 2015-02-16 | 2017-05-03 | 清华大学 | 在cbct中消除几何伪影的方法以及使用该方法的cbct系统 |
CN107831180A (zh) * | 2016-09-14 | 2018-03-23 | 奚岩 | X射线原位成像方法及系统 |
CN107945132A (zh) * | 2017-11-29 | 2018-04-20 | 深圳安科高技术股份有限公司 | 一种基于神经网络的ct图像的伪影校正方法及装置 |
CN107945132B (zh) * | 2017-11-29 | 2022-10-04 | 深圳安科高技术股份有限公司 | 一种基于神经网络的ct图像的伪影校正方法及装置 |
CN109187591B (zh) * | 2018-06-04 | 2020-10-02 | 东南大学 | 一种x射线超分辨成像方法及其应用 |
CN109187591A (zh) * | 2018-06-04 | 2019-01-11 | 东南大学 | 一种x射线超分辨成像方法及其应用 |
CN111175322B (zh) * | 2020-01-09 | 2022-08-16 | 昆山市建设工程质量检测中心 | 基于x射线数字成像技术的浆锚搭接连接质量检测方法 |
CN111175322A (zh) * | 2020-01-09 | 2020-05-19 | 江苏科技大学 | 基于x射线数字成像技术的浆锚搭接连接质量检测方法 |
CN111272782A (zh) * | 2020-03-13 | 2020-06-12 | 彭晟 | 自适应视场的x射线装置及自适应调整视场的方法 |
CN111458739A (zh) * | 2020-03-16 | 2020-07-28 | 北京东软医疗设备有限公司 | 成像方法、装置及系统 |
CN111458739B (zh) * | 2020-03-16 | 2021-11-30 | 北京东软医疗设备有限公司 | 成像方法、装置及系统 |
CN113029053A (zh) * | 2021-04-06 | 2021-06-25 | 中国科学技术大学 | 通用的ct对轴方法 |
CN113029053B (zh) * | 2021-04-06 | 2022-05-13 | 中国科学技术大学 | 通用的ct对轴方法 |
CN113137934A (zh) * | 2021-04-23 | 2021-07-20 | 中国工程物理研究院流体物理研究所 | 一种单探头射线方程的标定系统及标定方法 |
CN113137934B (zh) * | 2021-04-23 | 2022-10-28 | 中国工程物理研究院流体物理研究所 | 一种单探头射线方程的标定系统及标定方法 |
CN115060747A (zh) * | 2022-04-13 | 2022-09-16 | 中国兵器科学研究院宁波分院 | 一种工业ct系统用x射线源的焦点尺寸定量方法 |
CN115060747B (zh) * | 2022-04-13 | 2024-09-24 | 中国兵器科学研究院宁波分院 | 一种工业ct系统用x射线源的焦点尺寸定量方法 |
CN116892893A (zh) * | 2023-09-11 | 2023-10-17 | 上海福柯斯智能科技有限公司 | 工业ct锥束中心投影点量测方法及存储介质 |
CN116892893B (zh) * | 2023-09-11 | 2023-12-12 | 上海福柯斯智能科技有限公司 | 工业ct锥束中心投影点量测方法及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN101832954B (zh) | 2012-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101832954B (zh) | 锥束xct系统中进行图像重建坐标系原点标定的方法 | |
CN101726503B (zh) | 用于x射线相衬层析成像的系统和方法 | |
CN101515370B (zh) | 三维显微ct扫描系统中射线源焦点的投影坐标的标定方法 | |
CN104997529B (zh) | 基于对称重复的模板校正锥束ct系统几何失真的方法 | |
CN201199234Y (zh) | 分层摄影检测系统 | |
CN108444449B (zh) | 一种对具有平行线特征的目标空间姿态测量方法 | |
CN102155940B (zh) | 用于双目视觉定位跟踪系统的立体靶标 | |
CN100382763C (zh) | 一种适用于三维ct扫描系统投影坐标原点的标定方法 | |
CN102346922B (zh) | 航天遥感载荷成像几何畸变三维可视化仿真方法 | |
CN103226113B (zh) | 锥束3d-ct扫描系统重建体素尺寸的自动标定方法 | |
CN107367229A (zh) | 自由双目立体视觉转轴参数标定方法 | |
CN104390633B (zh) | 一种非接触式机构空间运动测量装置及实现方法 | |
CN103759679B (zh) | 一种锥束ct系统角度偏差测量方法 | |
CN106289283A (zh) | 利用激光扫描仪写入占用网格地图的系统和方法 | |
TW201514446A (zh) | 三維量測模擬取點系統及方法 | |
Ferrucci et al. | Measurement of the X-ray computed tomography instrument geometry by minimization of reprojection errors—Implementation on experimental data | |
CN103697811B (zh) | 一种相机与结构光源结合获取物体轮廓三维坐标的方法 | |
CN104154955B (zh) | 贮箱液体推进剂液面形貌与剂量动态测量方法和系统 | |
CN107991676A (zh) | 星载单航过InSAR系统对流层误差校正方法 | |
CN104048600A (zh) | 基于光耦探测器x射线三维显微镜重建体素尺寸标定方法 | |
Zhang et al. | Global Measurement Method for Large‐Scale Components Based on a Multiple Field of View Combination | |
CN103792567A (zh) | 一种实用的板状构件层析扫描装置射束倾角标定方法 | |
CN105469391A (zh) | 一种云阴影检测方法及系统 | |
CN201765199U (zh) | 一种适用于锥束xct系统中的二自由度转台 | |
CN105387818A (zh) | 一种基于一维图像序列的大尺度三维形貌测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120620 Termination date: 20130315 |