CN1861001A - 一种适用于三维ct扫描系统投影坐标原点的标定方法 - Google Patents
一种适用于三维ct扫描系统投影坐标原点的标定方法 Download PDFInfo
- Publication number
- CN1861001A CN1861001A CN 200610066252 CN200610066252A CN1861001A CN 1861001 A CN1861001 A CN 1861001A CN 200610066252 CN200610066252 CN 200610066252 CN 200610066252 A CN200610066252 A CN 200610066252A CN 1861001 A CN1861001 A CN 1861001A
- Authority
- CN
- China
- Prior art keywords
- coordinate
- subpoint
- point
- detector
- circle
- 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
- Analysing Materials By The Use Of Radiation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种适用于三维CT扫描系统投影坐标原点的标定方法,其主要针对基于FDK算法的圆轨迹扫描成像系统投影坐标原点的标定。由于射线源焦点P和探测器成像平面的准确空间位置无法直接测量得到,从而使得投影坐标原点位置难以精确测量。本发明提出一种基于二次成像的双圆最小二乘拟合方法实现对投影坐标原点的高精度标定,对空间双圆目标体进行二次成像,利用图像、图形处理方法和最小二乘拟合技术求取不同成像位置下的双圆目标体投影圆心,利用圆心坐标解算出射线源焦点P在成像平面上的投影点O坐标,将本发明的标定结果应用于三维CT扫描设备的图像重构单元和CT控制单元,图像重建质量达到了设备出厂时的指标。
Description
技术领域
本发明涉及一种适用于三维CT扫描系统投影坐标原点的标定方法,是针对基于面阵成像器的X-射线3D-CT(Three Dimensional Computed Tomography)扫描成像系统投影坐标原点的标定,可用于医学和工业领域射线数字成像(DR-DigitalRadiography)、三维计算机断层扫描(3D-CT)成像过程中的相关测量。
背景技术
近些年来,随着计算机技术的飞速发展和面阵探测器件的出现,3D-CT日益成为NDT(Non-destructive Testing)领域内的研究热点,它采用锥束射线开放式扫描,利用探测器采集的二维数字投影序列,快速重构出扫描区域内的所有断层,真正实现对被检对象的三维全息。
在众多的三维重建算法中,考虑到运算量和工程实现难度,FDK(Feldkamp-Davis-Kress)(Feldkamp L A,Davis L C,Kress J W.Practicalcone-beam algorithm[J].Opt.Soc.Am.Al,1984:612~619)类型的算法最为实用,也一直是实际工程应用中的主流,其扫描原理如图1所示,射线源2沿圆轨迹对物体9进行扫描等效于射线源2、探测器4不动,物体9绕轴线旋转,利用探测器4获取的二维投影序列进行三维重建,一般工业低能X射线CT系统的射线源2与探测器4之间的距离为1~2.5m。
在图1中,探测器4的成像平面5坐标系为xdydzd。当探测器4在成像系统中的位置固定后,该坐标系即可确定,其原点Od即为探测器4采集的射线数字图像(即DR图像)坐标系原点。FDK扫描方式的三维CT重建在锥束射线3和探测器4构造的几何坐标系中进行,该坐标系被称为投影坐标系,即图1中的xyz坐标系,该坐标系各个坐标轴与坐标系xdydzd的对应坐标轴平行。然而,在实际的物理成像系统中,如何对xyz坐标系进行标定一直没有很好的解决方法,其关键就是如何准确定位坐标系的原点O在探测器4的成像平面5坐标系xdydzd中的位置。坐标原点O的误差会造成重建图像出现伪影和图像的几何失真,同时影响断层轴向(即z方向)定位的精度,尤其对于显微CT(Micro-Computed Tomgraphy),射线源焦点尺寸仅为几微米到几十微米,而且在扫描成像中充分地用到了几何放大(即物体9尽可能靠近射线源2),使得三维重建对投影坐标系原点O的定位精度要求更高。在实际的成像系统中,投影坐标系原点O即为射线源焦点P在探测器4的成像平面5(xsOdzd)上的投影点,由于射线源2的焦点、探测器4的成像平面5的准确位置无法直接测量得到,从而使得射线源2的焦点P在成像平面5(xdOdzd)上的投影点O的位置难以精确测量。
在目前已公开的文献中,提出了非线性最小二乘估计的测量方法。其思想就是计算空间一质点在不同转角下的投影坐标,质点理论投影坐标与实际求得的投影坐标之间应满足误差最小,通过求解满足此条件的非线性最小二乘解,即得到投影坐标原点参数值。在参数向量初始值接近于真值的情况下,最小二乘估计所得参数值有较好的重复性,但在实际的实现过程中,会存在下述困难:(A)不易得到空间质点的投影数据,即使能够对类似质点的对象进行锥束扫描,但追踪其投影质心坐标时会遇到较为复杂的图像、图形处理技术;(B)非线性最小二乘估计结果受参数向量的初始值影响较大。
发明内容
本发明的目的是提供一种适用于三维CT扫描系统投影坐标原点的标定方法,该标定方法是基于二次成像的双圆最小二乘拟合方法实现对投影坐标原点O的高精度标定,对空间双圆目标体进行二次成像,利用图像、图形处理方法和双圆最小二乘拟合技术求取不同成像位置下双圆目标体圆心的投影点(即圆形物体的投影圆心),利用圆心坐标计算得到射线源焦点在成像平面中的投影原点O坐标。
本发明是一种适用于三维CT扫描系统投影坐标原点的标定方法,其通过在位于射线源和探测器之间的多自由度载物台上放置一板状目标;目标上任意安装两个圆形物体;当射线源射出的锥束射线照射到目标上时,通过移动多自由度载物台,探测器采集到两个成像位置的双圆物体的射线数字图像(即DR图像),通过最小二乘拟合法解算出投影坐标原点O的坐标(λx,o,λz);具体步骤如下:
(A)调整射线源的中心射线PO与探测器的成像平面垂直,所述成像平面的坐标系记为xdydzd,射线源焦点P(λx,λy,λz)在所述成像平面中的投影点记为O(λx,o,λz),该点即为本发明所要标定的三维CT扫描系统的投影坐标原点;
(B)调整目标所在平面与探测器的成像平面平行,多自由度载物台在成像位置A点时,探测器采集目标上双圆物体的DR图像;其所述圆形物体在坐标系xdydzd下的圆心位于位置点E1和位置点F1,所述位置点E1和位置点F1在探测器的成像平面上的对应投影分别记为投影点B1和投影点C1;
(C)将多自由度载物台沿坐标系xdydzd的yd方向在射线源与探测器之间平移距离d后到达成像位置B点,探测器再次采集目标的DR图像;其所述圆形物体在坐标系xdydzd下的圆心位于位置点E2和位置点F2,所述位置点E2和位置点F2在探测器(4)的成像平面(5)上的对应投影分别记为投影点B2和投影点C2;
(D)对上述成像位置A点的双圆物体的DR图像进行双圆最小二乘拟合,回归出所述圆形物体圆心位置点E1、位置点F1的投影点坐标B1(xb1,0,zb1)、C1(xc1,0,zc1);
(E)对上述成像位置B点的双圆物体的DR图像进行双圆最小二乘拟合,回归出所述圆形物体圆心位置点E2、位置点F2的投影点坐标B2(xb2,0,zb2)、C2(xc2,0,zc2);
(F)根据步骤(D)和步骤(E)得到的坐标值xb1、xc1、zb1、zc1、xb2、xc2、zb2、zc2解算出投影坐标原点O的坐标(λx,o,λz)。
所述的三维CT扫描系统投影坐标原点的标定方法,其二次成像双圆最小二乘拟合求解所得投影坐标原点O的xd轴坐标满足
yd轴坐标满足λy=0,zd轴坐标满足
本发明标定方法的优点在于:(1)通过采集同一目标在两个不同位置的DR图像,实现对投影坐标原点的标定,该操作方法简单易行;(2)成像目标体为同一平面上的双圆,比加工一近似质点目标体更为容易,其制造成本相对较低;(3)利用最小二乘拟合法拟合双圆圆心的投影坐标,坐标值精度达到亚像素级,提高了标定的精度。
附图说明
图1是基于FDK重建算法的圆轨迹扫描原理图。
图2是本发明设计的二次成像系统原理图。
图2A是成像位置A点的空间几何关系示意图。
图2B是成像位置B点的空间几何关系示意图。
图2C二次成像的空间解析图。
图3是双圆最小二乘拟合方法求取投影坐标原点的流程图。
图中:1.目标 2.射线源 3.锥束射线 4.探测器 5.成像平面 6.多自由度载物台 7.成像位置A 8.成像位置B 9.物体
具体实施方式
下面将结合附图对本发明作进一步的详细说明。
本发明提出的一种适用于三维CT扫描系统投影坐标原点的标定方法,是出厂前对三维CT扫描设备的投影坐标原点O进行标定,或者是三维CT扫描设备经一段时间使用后,对扫描系统投影坐标原点O的位置进行校正时而采用的一种较为简便、易操作的投影坐标原点的标定方法。一套三维CT扫描设备一般由硬件部分和软件部分组成,其中,硬件部分包括:射线源、多自由度载物台、探测器、控制器、PC机;软件部分包括:CT控制单元、图像重构单元、图像处理与可视化单元。本发明对射线源2的焦点P投射在探测器4的成像平面5上的投影点O的位置,采用基于二次成像的双圆最小二乘拟合法进行测量,该测量的结果应用于图像重构单元和CT控制单元中。
本发明是一种适用于三维CT扫描系统投影坐标原点的标定方法,通过在位于射线源2和探测器4(面阵探测器)之间的多自由度载物台6上放置一板状目标1,目标1上任意安装两个圆形物体,其结构原理如图2所示。当射线源2射出的锥束射线3照射到目标1上时,通过移动多自由度载物台6(载物台6是平行移动,即在射线源2与探测器4之间平行移动,其移动距离在本发明中为任意),探测器4采集到两个成像位置的双圆物体的射线数字图像(即DR图像),通过最小二乘拟合解算出投影坐标原点O的坐标(λx,o,λz)。具体步骤如下:
(A)调整射线源2的中心射线PO与探测器4的成像平面5垂直,所述成像平面5的坐标系记为xdydzd,射线源焦点P(λx,λy,λz)在所述成像平面5上的投影点记为O(λx,o,λz),该点即为本发明所要标定的三维CT扫描系统的投影坐标原点;
(B)(请参见图2A所示)调整目标1所在平面与探测器4的成像平面5平行。当多自由度载物台6在成像位置A点7时,探测器4采集目标1上双圆物体的DR图像;所述双圆物体在坐标系xdydzd下的圆心位于位置点E1和位置点F1,所述位置点E1和位置点F1在DR图像上的对应投影分别记为投影点B1和投影点C1,投影点B1和投影点C1亦为位置点E1和位置点F1在探测器4的成像平面5上的投影点;
(C)(请参见图2B)将多自由度载物台6沿坐标系xdydzd的yd方向(正负均可,在本发明中,向射线源2方向移动称为负,向探测器4方向移动称为正)在射线源2与探测器4之间平移任意距离d后到达成像位置B点8时,探测器4再次采集目标1的DR图像;所述双圆物体在坐标系xdydzd下的圆心位于位置点E2和位置点F2,所述位置点E2和位置点F2在DR图像上的对应投影分别记为投影点B2和投影点C2,投影点B2和投影点C2亦为位置点E2和位置点F2在探测器4的成像平面5上的投影点;
(D)(请参见图3)对成像位置A点7的双圆物体的DR图像进行双圆最小二乘拟合,回归出所述圆形物体圆心位置点E1、位置点F1的投影点坐标B1(xb1,0,zb1)、C1(xc1,0,zc1);
(E)(请参见图3)对成像位置B点8的双圆物体的DR图像进行双圆最小二乘拟合,回归出所述圆形物体圆心位置点E2、位置点F2的投影点坐标B2(xb2,0,zb2)、C2(xc2,0,zc2);
(F)根据步骤(D)和步骤(E)利用双圆最小二乘拟合得到的坐标值xb1、xc1、zb1、zc1、xb2、xc2、zb2、zc2解算得到投影坐标原点O的坐标(λx,o,λz)。
请参见图2C所示,本发明的二次成像原理是以射线源2的中心射线PO垂直于探测器4的成像平面5为前提,射线源焦点P(λx,λy,λz)的投影点为O(λx,o,λz)。当多自由度载物台6的位置设定在成像位置A点7时,连接目标1上的位置点E1和位置点F1,则有直线E1F1//XdOdZd(XdOdZd表示成像平面5),位置点E1在探测器4的成像平面5坐标系中对应的投影坐标为B1(xb1,0,zb1),位置点F1在探测器4的成像平面5坐标系中对应的投影坐标为C1(xc1,0,zc1);多自由度载物台6在射线源与探测器之间沿yd方向平移(正负方向均可)一段距离d到达成像位置B点8时,位置点E1移至位置点E2,位置点F1移至位置点F2。连接目标1上的位置点E2和位置点F2,则有直线E2F2//XdOdZd(XdOdZd表示成像平面5),位置点E2在探测器4的成像平面5坐标系中对应的投影坐标为B2(xb2,0,zb2),位置点F2在探测器4的成像平面5坐标系中对应的投影坐标为C2(xc2、0,zc2);当在成像位置A点7时,设位置点E1、位置点F1的yd轴坐标为d1,当在成像位置B点8时,设位置点E2、位置点F2的yd轴坐标为d2;
因此,在成像位置A点7的射线PB1、射线PC1的直线方程分别为:
由式(1)、式(2)得到位置点E1和位置点F1的坐标:
当移动至成像位置B点8时,位置点E2和位置点F2的坐标:
由式(5)、式(6)得到射线PB2、射线PC2的直线方程分别为:
联立式(7)和式(8)解算得到:
上述各式中引入字母的物理意义为:xb1表示投影点B1的xd轴坐标,zb1表示投影点B1的zd轴坐标,xc1表示投影点C1的xd轴坐标,zc1表示投影点C1的zd轴坐标,xb2表示投影点B2的xd轴坐标,zb2表示投影点B2的zd轴坐标,xc2表示投影点C2的xd轴坐标,zc2表示投影点C2的zd轴坐标,d1表示位置点E1和位置点F1的yd轴坐标,d2表示位置点E2和位置点F2的yd轴坐标,λx表示射线源焦点P的xd轴坐标,λy表示射线源焦点P的yd轴坐标,λz表示射线源焦点P的zd轴坐标。
在实际应用中,将由步骤(D)和(E)回归得到的xb1、xc1、zb1、zc1、xb2、xc2、zb2、zc2代入式(9)、式(10)即可得到投影坐标原点O坐标值(λx,o,λz)。
图3为本发明采用的最小二乘拟合方法求取投影坐标原点的流程图,探测器4采集到双圆目标体的DR图像亦为双圆,对双圆DR图像进行边缘检测、阈值分割、轮廓细化与轮廓追踪,以得到DR图像中双圆的轮廓点坐标,然后利用最小二乘拟合法回归出DR图像中双圆的圆心坐标。回归出的圆心坐标即为计算投影坐标原点O的坐标(λx,o,λz)时所必须的xb1、xc1、zb1、zc1、xb2、xc2、zb2、zc2坐标值。该流程图中边缘检测、阈值分割、轮廓细化与轮廓追踪、最小二乘拟合求取圆心坐标模块为已公开的图像、图形处理方法。
利用计算机模拟对该标定方法的测量精度进行考核,模拟了λx=256,λz=256,λy=1000情况下的双圆DR图像,且图像叠加1%的高斯噪声。两次成像位置之差Δd=|d1-d2|分别为100、200、300、400时的计算结果见表1。Δλx、Δλz为绝对误差,err_λx、err_λz为相对误差,λx、λz的平均相对误差分别为0.0481%和0.0322%。
表1利用计算机模拟对参数(λx,λz)的计算结果(单位:像素)
Δd | λx | Δλx | λz | Δλz | err_λx | err_λz |
100 | 256.1432 | 0.1432 | 256.1747 | 0.1747 | 0.05559% | 0.0682% |
200 | 256.1248 | 0.1248 | 256.1328 | 0.1328 | 0.0488% | 0.0519% |
300 | 256.1092 | 0.1092 | 256.0041 | 0.0041 | 0.0427% | 0.0016% |
400 | 256.1153 | 0.1153 | 256.0177 | 0.0177 | 0.0450% | 0.0069% |
E(λx)=256.1231,E(λz)=256.0823,|E(λx)-λx|/λx=0.0481%,|E(λz)-λz|/λz=0.0322% |
注:E(λx)表示λx的平均值,E(λz)表示λz的平均值。
为了校正使用一段时间后的三维CT扫描设备投影坐标原点的偏差,本发明人根据图2所示的原理结构制作了双圆目标体,并通过二次成像获得该目标体的DR图像(参见图2A、图2B所示),利用该DR图像运用如图3所示的求取投影坐标原点的流程图,获得如表2所示的标定结果。
表2本发明标定方法对实际成像系统坐标原点的标定结果
λx | 969.5816 | 971.1252 | 970.3814 | 970.5180 | 971.5055 | 970.6422 | 970.2056 |
E(λx)=970.566,D(λx)=0.624 | |||||||
λz | 923.8883 | 924.1236 | 925.1794 | 924.2718 | 924.6505 | 924.3364 | 924.4902 |
E(λz)=924.4200,D(λz)=0.415 |
注:Δd=20mm,E(λx)表示λx的平均值,E(λz)表示λz的平均值,D(λx)为λx的标准方差,D(λz)为λz的标准方差。
将表2的标定结果应用于三维CT扫描设备的图像重构单元和CT控制单元,图像重建质量达到了设备出厂时的指标。
Claims (3)
1、一种适用于三维CT扫描系统投影坐标原点的标定方法,其特征在于:通过在位于射线源(2)和探测器(4)之间的多自由度载物台(6)上放置一板状目标(1);目标(1)上任意安装两个圆形物体;当射线源(2)射出的锥束射线(3)照射到目标(1)上时,通过移动多自由度载物台(6),探测器(4)采集到两个成像位置的双圆物体的DR图像,通过最小二乘拟合法联立解算出投影坐标原点O的坐标(λx,o,λz);具体步骤如下:
(A)调整射线源(2)的中心射线PO与探测器(4)的成像平面(5)垂直,所述成像平面(5)的坐标系记为xdydzd,射线源焦点P(λx,λy,λz)在所述成像平面(5)上的投影点记为O(λx,o,λz);
(B)调整目标(1)所在平面与探测器(4)的成像平面(5)平行,多自由度载物台(6)在成像位置A点(7)时,探测器(4)采集目标(1)上双圆物体的DR图像;其所述圆形物体在坐标系xdydzd下的圆心位于位置点E1和位置点F1,所述位置点E1和位置点F1在探测器(4)的成像平面(5)上的对应投影分别记为投影点B1和投影点C1;
(C)将多自由度载物台(6)沿坐标系xdydzd的yd方向在射线源(2)与探测器(4)之间平移距离d后到达成像位置B点(8)时,探测器(4)再次采集目标(1)上双圆物体的DR图像;其所述圆形物体在xdydzd坐标系下的圆心位于位置点E2和位置点F2,所述位置点E2和位置点F2在探测器(4)的成像平面(5)上的对应投影分别记为投影点B2和投影点C2;
(D)对上述成像位置A点(7)的双圆物体的DR图像进行双圆最小二乘拟合,回归出所述圆形物体的圆心位置点E1、位置点F1的投影点坐标B1(xb1,0,zb1)、C1(xc1,0,zc1);
(E)对上述成像位置B点(8)的双圆物体的DR图像进行双圆最小二乘拟合,回归出所述圆形物体的圆心位置点E2、位置点F2的投影点坐标B2(xb2,0,zb2)、C2(xc2,0,zc2);
(F)根据步骤(D)和步骤(E)得到的坐标值xb1、xc1、zb1、zc1、xb2、xc2、zb2、zc2解算出投影坐标原点O的坐标(λx,o,λz)。
2、根据权利要求1所述的三维CT扫描系统投影坐标原点的标定方法,其特征在于:所述投影坐标原点O的xd轴坐标满足
yd轴坐标满足λy=0,zd轴坐标满足
式中,xb1表示投影点B1的xd轴坐标,zb1表示投影点B1的zd轴坐标,xc1表示投影点C1的xd轴坐标,zc1表示投影点C1的zd轴坐标,xb2表示投影点B2的xd轴坐标,zb2表示投影点B2的zd轴坐标,xc2表示投影点C2的xd轴坐标,zc2表示投影点C2的zd轴坐标。
3、根据权利要求1所述的三维CT扫描系统投影坐标原点的标定方法,其特征在于:所述探测器(4)为平面阵列探测器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB200610066252XA CN100382763C (zh) | 2006-03-31 | 2006-03-31 | 一种适用于三维ct扫描系统投影坐标原点的标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB200610066252XA CN100382763C (zh) | 2006-03-31 | 2006-03-31 | 一种适用于三维ct扫描系统投影坐标原点的标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1861001A true CN1861001A (zh) | 2006-11-15 |
CN100382763C CN100382763C (zh) | 2008-04-23 |
Family
ID=37388598
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB200610066252XA Expired - Fee Related CN100382763C (zh) | 2006-03-31 | 2006-03-31 | 一种适用于三维ct扫描系统投影坐标原点的标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100382763C (zh) |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101832954A (zh) * | 2010-03-15 | 2010-09-15 | 中国工程物理研究院应用电子学研究所 | 锥束xct系统用移动组件以及用其进行图像重建坐标系原点标定的方法 |
CN101936720A (zh) * | 2010-07-30 | 2011-01-05 | 北京航空航天大学 | 一种适用于锥束xct系统的探测器扭转角的标定方法 |
CN101750021B (zh) * | 2009-12-04 | 2011-05-11 | 深圳先进技术研究院 | Ct系统中几何参数的标定方法、装置 |
CN101515370B (zh) * | 2009-03-06 | 2011-05-18 | 北京航空航天大学 | 三维显微ct扫描系统中射线源焦点的投影坐标的标定方法 |
CN103006251A (zh) * | 2012-12-06 | 2013-04-03 | 深圳先进技术研究院 | 用于ct系统中标定几何参数的标定体模、标定装置及标定方法 |
CN101842807B (zh) * | 2007-05-31 | 2013-07-31 | 伊利克塔股份有限公司 | Ct扫描中运动伪影的减少 |
CN104132950A (zh) * | 2014-07-18 | 2014-11-05 | 中国特种设备检测研究院 | 基于原始投影信息的cl扫描装置投影旋转中心标定方法 |
CN104257397A (zh) * | 2014-09-30 | 2015-01-07 | 清华大学 | 基于层析成像的x光机与探测器几何位置关系的标定方法 |
WO2015074254A1 (zh) * | 2013-11-19 | 2015-05-28 | 中国科学院过程工程研究所 | Ct探测器偏转角的确定方法和装置 |
CN105997126A (zh) * | 2016-05-25 | 2016-10-12 | 重庆大学 | 一种锥束ct系统几何参数校正模型及方法 |
CN107036563A (zh) * | 2016-02-03 | 2017-08-11 | 上海西门子医疗器械有限公司 | 确定投影角精度和投影角的方法和计算机断层扫描设备 |
WO2018218611A1 (zh) * | 2017-05-27 | 2018-12-06 | 广州华端科技有限公司 | 锥形束计算机断层成像系统的几何参数确定方法 |
CN108982556A (zh) * | 2018-08-22 | 2018-12-11 | 武汉科技大学 | 一种ct参数标定体膜、ct参数标定系统及ct参数标定方法 |
CN109146992A (zh) * | 2018-08-24 | 2019-01-04 | 吴玉新 | 一种断层扫描平面重现的方法及装置 |
CN109682843A (zh) * | 2019-02-13 | 2019-04-26 | 重庆交通大学 | 一种对ct系统的参数标定方法 |
CN110503713A (zh) * | 2019-07-03 | 2019-11-26 | 西安电子科技大学 | 一种基于轨迹平面法向量和圆心结合的旋转轴估计方法 |
CN110517330A (zh) * | 2019-08-07 | 2019-11-29 | 北京航空航天大学 | 一种偏置扫描模式下的工业锥束ct重建方法 |
CN111458739A (zh) * | 2020-03-16 | 2020-07-28 | 北京东软医疗设备有限公司 | 成像方法、装置及系统 |
CN113587810A (zh) * | 2021-07-20 | 2021-11-02 | 苏州工业园区智在天下科技有限公司 | 一种生成光源位置的方法和装置 |
CN113963056A (zh) * | 2021-09-07 | 2022-01-21 | 于留青 | Ct图像重建方法、装置、电子设备以及存储介质 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2700909B1 (fr) * | 1993-01-27 | 1995-03-17 | Gen Electric Cgr | Dispositif et procédé automatique de calibration géométrique d'un système d'imagerie par rayons X. |
JP3548306B2 (ja) * | 1995-12-22 | 2004-07-28 | 株式会社日立メディコ | X線断層撮影装置 |
FR2822273B1 (fr) * | 2001-03-13 | 2003-07-11 | Ge Med Sys Global Tech Co Llc | Procede d'etalonnage pour la reconstruction de modelisations tri-dimensionnelles a partir d'images obtenues par tomographie |
DE10140867B4 (de) * | 2001-08-21 | 2005-08-18 | Siemens Ag | Kalibrierphantom für projektive Röntgensysteme |
DE10215808B4 (de) * | 2002-04-10 | 2005-02-24 | Siemens Ag | Verfahren zur Registrierung für navigationsgeführte Eingriffe |
EP1654516B1 (en) * | 2003-08-08 | 2012-11-28 | University Health Network | Method and system for calibrating a source and detector instrument |
CN100381103C (zh) * | 2005-03-22 | 2008-04-16 | 东软飞利浦医疗设备系统有限责任公司 | 一种ct机几何参数校正方法 |
-
2006
- 2006-03-31 CN CNB200610066252XA patent/CN100382763C/zh not_active Expired - Fee Related
Cited By (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101842807B (zh) * | 2007-05-31 | 2013-07-31 | 伊利克塔股份有限公司 | Ct扫描中运动伪影的减少 |
CN101515370B (zh) * | 2009-03-06 | 2011-05-18 | 北京航空航天大学 | 三维显微ct扫描系统中射线源焦点的投影坐标的标定方法 |
CN101750021B (zh) * | 2009-12-04 | 2011-05-11 | 深圳先进技术研究院 | Ct系统中几何参数的标定方法、装置 |
CN101832954A (zh) * | 2010-03-15 | 2010-09-15 | 中国工程物理研究院应用电子学研究所 | 锥束xct系统用移动组件以及用其进行图像重建坐标系原点标定的方法 |
CN101936720B (zh) * | 2010-07-30 | 2012-03-21 | 北京航空航天大学 | 一种适用于锥束xct系统的探测器扭转角的标定方法 |
CN101936720A (zh) * | 2010-07-30 | 2011-01-05 | 北京航空航天大学 | 一种适用于锥束xct系统的探测器扭转角的标定方法 |
CN103006251A (zh) * | 2012-12-06 | 2013-04-03 | 深圳先进技术研究院 | 用于ct系统中标定几何参数的标定体模、标定装置及标定方法 |
WO2015074254A1 (zh) * | 2013-11-19 | 2015-05-28 | 中国科学院过程工程研究所 | Ct探测器偏转角的确定方法和装置 |
CN104132950B (zh) * | 2014-07-18 | 2016-07-06 | 中国特种设备检测研究院 | 基于原始投影信息的cl扫描装置投影旋转中心标定方法 |
CN104132950A (zh) * | 2014-07-18 | 2014-11-05 | 中国特种设备检测研究院 | 基于原始投影信息的cl扫描装置投影旋转中心标定方法 |
CN104257397A (zh) * | 2014-09-30 | 2015-01-07 | 清华大学 | 基于层析成像的x光机与探测器几何位置关系的标定方法 |
CN104257397B (zh) * | 2014-09-30 | 2016-08-24 | 清华大学 | 基于层析成像的x光机与探测器几何位置关系的标定方法 |
CN107036563A (zh) * | 2016-02-03 | 2017-08-11 | 上海西门子医疗器械有限公司 | 确定投影角精度和投影角的方法和计算机断层扫描设备 |
CN105997126A (zh) * | 2016-05-25 | 2016-10-12 | 重庆大学 | 一种锥束ct系统几何参数校正模型及方法 |
CN105997126B (zh) * | 2016-05-25 | 2019-04-02 | 重庆大学 | 一种锥束ct系统几何参数校正模型及方法 |
WO2018218611A1 (zh) * | 2017-05-27 | 2018-12-06 | 广州华端科技有限公司 | 锥形束计算机断层成像系统的几何参数确定方法 |
CN108982556B (zh) * | 2018-08-22 | 2020-11-13 | 武汉科技大学 | 一种ct参数标定体膜、ct参数标定系统及方法 |
CN108982556A (zh) * | 2018-08-22 | 2018-12-11 | 武汉科技大学 | 一种ct参数标定体膜、ct参数标定系统及ct参数标定方法 |
CN109146992A (zh) * | 2018-08-24 | 2019-01-04 | 吴玉新 | 一种断层扫描平面重现的方法及装置 |
CN109682843A (zh) * | 2019-02-13 | 2019-04-26 | 重庆交通大学 | 一种对ct系统的参数标定方法 |
CN109682843B (zh) * | 2019-02-13 | 2021-07-06 | 重庆交通大学 | 一种对ct系统的参数标定方法 |
CN110503713A (zh) * | 2019-07-03 | 2019-11-26 | 西安电子科技大学 | 一种基于轨迹平面法向量和圆心结合的旋转轴估计方法 |
CN110503713B (zh) * | 2019-07-03 | 2023-08-15 | 西安电子科技大学 | 一种基于轨迹平面法向量和圆心结合的旋转轴估计方法 |
CN110517330A (zh) * | 2019-08-07 | 2019-11-29 | 北京航空航天大学 | 一种偏置扫描模式下的工业锥束ct重建方法 |
CN111458739A (zh) * | 2020-03-16 | 2020-07-28 | 北京东软医疗设备有限公司 | 成像方法、装置及系统 |
CN111458739B (zh) * | 2020-03-16 | 2021-11-30 | 北京东软医疗设备有限公司 | 成像方法、装置及系统 |
CN113587810A (zh) * | 2021-07-20 | 2021-11-02 | 苏州工业园区智在天下科技有限公司 | 一种生成光源位置的方法和装置 |
CN113963056A (zh) * | 2021-09-07 | 2022-01-21 | 于留青 | Ct图像重建方法、装置、电子设备以及存储介质 |
CN113963056B (zh) * | 2021-09-07 | 2022-08-26 | 于留青 | Ct图像重建方法、装置、电子设备以及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN100382763C (zh) | 2008-04-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1861001A (zh) | 一种适用于三维ct扫描系统投影坐标原点的标定方法 | |
US20150260859A1 (en) | Method and device for correcting computed tomographiy measurements, comprising a coordinate measuring machine | |
CN107328798B (zh) | 一种新型icl系统及实现方法 | |
JP2016538552A (ja) | コンピュータ断層撮影の較正装置および方法 | |
CN1888814A (zh) | 三维主动视觉传感器的多视点姿态估计和自标定方法 | |
CN107796834B (zh) | 一种正交电子直线扫描cl成像系统及方法 | |
WO2018126335A1 (zh) | 基于小球模体的锥束ct系统几何参数评价及校正方法 | |
JP2011516870A (ja) | 機械的ワークを断層撮影法によって測定するための方法 | |
Bircher et al. | CT geometry determination using individual radiographs of calibrated multi-sphere standards | |
Dewulf et al. | Enhanced dimensional measurement by fast determination and compensation of geometrical misalignments of X-ray computed tomography instruments | |
CN103226113A (zh) | 锥束3d-ct扫描系统重建体素尺寸的自动标定方法 | |
Muralikrishnan et al. | X-ray computed tomography instrument performance evaluation, Part I: Sensitivity to detector geometry errors | |
CN101013065A (zh) | 一种基于像素频率的星敏感器高精度校准方法 | |
Ferrucci | Towards traceability of CT dimensional measurements | |
Ferrucci | Systematic approach to geometrical calibration of X-ray computed tomography instruments | |
CN114596222A (zh) | 适用于一般轨迹锥束ct系统几何校正的模体与标定方法 | |
Ferrucci et al. | Measurement of sample stage error motions in cone-beam x-ray computed tomography instruments by minimization of reprojection errors | |
Kang et al. | Geometric qualification for robot CT with flexible trajectories | |
CN101615293B (zh) | Vct系统参数标定装置及方法 | |
Zhang et al. | Precise on-line non-target pose measurement for cylindrical components based on laser scanning | |
CN112581524A (zh) | 基于点云的slam移动机器人机场道路检测数据获取方法 | |
JP5022820B2 (ja) | 放射線検査装置及びそれを用いた配管検査方法 | |
Duan et al. | Knowledge-based self-calibration method of calibration phantom by and for accurate robot-based CT imaging systems | |
Blumensath et al. | Calibration of robotic manipulator systems for cone-beam tomography imaging | |
US8897534B2 (en) | Method and evaluation device for determining the position of a structure located in an object to be examined by means of X-ray computer tomography |
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: 20080423 |