CN108201447A - 一种静态ct系统几何参数校正方法 - Google Patents
一种静态ct系统几何参数校正方法 Download PDFInfo
- Publication number
- CN108201447A CN108201447A CN201711165876.1A CN201711165876A CN108201447A CN 108201447 A CN108201447 A CN 108201447A CN 201711165876 A CN201711165876 A CN 201711165876A CN 108201447 A CN108201447 A CN 108201447A
- Authority
- CN
- China
- Prior art keywords
- detector
- geometric parameter
- radiographic source
- static state
- correction
- 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.)
- Pending
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/02—Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computerised tomographs
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating apparatus or devices for radiation diagnosis
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2560/00—Constructional details of operational features of apparatus; Accessories for medical measuring apparatus
- A61B2560/02—Operational features
- A61B2560/0223—Operational features of calibration, e.g. protocols for calibrating sensors
Abstract
本发明提出了一种静态CT系统几何参数校正方法,所述静态CT系统由多个射线源和一个以上探测器构成,所述几何参数为每个射线源与对应的探测器之间的几何参数,所述方法包括如下步骤:步骤a:将校正体模放置在需要校正几何参数的静态CT系统中,所述校正体模具有标记小球;步骤b:每个射线源触发一次射线,获得多幅投影;以及步骤c:根据每幅投影中标记小球的投影中心点的坐标,结合已知的标记小球在待校正系统中的坐标,计算获得整个静态CT系统的校正后的几何参数。本发明可以同时适用于离体CT系统和活体CT系统,并可以同时校正CT系统的全部10个几何参数。另外,本发明仅需要各射线源对校正体模进行1次投影,即可得到CT重建所需的全部几何参数。
Description
技术领域
本发明涉及医学和工业领域的计算机断层扫描(CT)系统,尤其涉及一种静态CT系统几何参数校正方法。
背景技术
计算机断层成像(CT)是通过无损方式获取物体内部结构信息的一种重要成像手段,它拥有高分辨率、高灵敏度以及多层次等众多优点,被广泛应用于各个医疗临床检查领域。相对于传统的CT系统,静态CT系统具有多个射线源,在数据的采集过程中不需要X射线源和X射线探测器等部件旋转,可以降低成本,缩短CT成像扫描时间的优势,并且由于静态CT在成像的过程中没有机械运动,就避免了运动伪影的产生。
几何参数校正是静态CT系统调试的重要组成部分,也是获得精确CT图像的前提条件。中国专利申请号201210148432.8提出了“一种X射线锥束计算机层析成像系统的几何参数校正方法”。该方法通过相机校正技术来复现锥束CT系统中的X射线源、平板探测器和转轴之间的几何位置关系,从而对锥束CT的几何参数及其误差进行直接求解,是一种系统化的测量求解方法,可以将锥束CT抽象为基本的相机系统模型,从而同时求解系统中相互关联的多个几何参数。
但是,上述技术存在如下缺点:1、不能同时适用于离体CT系统(样本旋转,射线源和探测器不动)和活体CT系统(样本不动,射线源和探测器旋转);2、不能同时校正CT系统的全部10个几何参数;3、参数校正过程较复杂(例如,需要对校正体模进行多次投影),难于多次重复实现。
因此,需要一种简单、有效、实用的静态CT系统快速几何参数校正方法。
发明内容
针对现有技术的上述问题,本发明提出了一种静态CT系统几何参数校正方法,可以同时提取成像系统的10个重要的系统几何参数,并设计了适用于CT系统校正的体模装置。
本发明的静态CT系统几何参数校正方法,所述静态CT系统由多个射线源和一个以上探测器构成,所述几何参数为每个射线源与对应的探测器之间的几何参数,所述方法包括如下步骤:
步骤a:将校正体模放置在需要校正几何参数的静态CT系统中,所述校正体模具有标记小球;
步骤b:每个射线源触发一次射线,获得多幅投影;以及
步骤c:根据每幅投影中标记小球的投影中心点的坐标,结合已知的标记小球在待校正系统中的坐标,计算获得整个静态CT系统的校正后的几何参数。
优选地,所述几何参数为所述静态CT系统的10个几何参数,包括:射线源S的坐标位置(Ox,Oy,Oz)、射线源S在探测器平面的投影中心点坐标(u0,v0)、探测器的扭转角θ、探测器的倾斜角φ、探测器的旋转角η、射线源到扫描物体中心的距离(SOD),射线源到探测器的距离(SDD)。
优选地,所述校正体模由支撑体和安装于所述支撑体上的标记小球构成,所述支撑体为低吸收材料,所述标记小球为高吸收材料。
优选地,所述支撑体为碳纤维材料或者聚甲基丙烯酸甲酯,所述标记小球为钢珠。
优选地,所述标记小球的数量、大小以及位置以校正过程中在探测器上的投影不造成重叠为宜。
优选地,所述步骤c进一步包括如下步骤:
步骤c1:对于每幅投影图,在小球投影坐标(u,v)和实际校正体模中的标记小球的三维坐标(x,y,z)之间利用一个3×4投影矩阵A建立映射关系:
[Ax,Ay,Az,A]T=[au,av,a]T (1)
其中,α为一比例参数,投影矩阵A的构成如下:
A=K[S|t] (2)
上述(2)式中,K为一个3×3上三角矩阵,S为一个3×3旋转矩阵,t为3×1位移向量,
上述(3)式中,u0和v0为射线源在探测器平面的投影点坐标,λ为探测器像素的大小,
公式(2)中旋转矩阵S可以进一步用三角函数表达为:
θ为探测器的扭转角、φ为探测器的倾斜角、η为探测器的旋转角,
公式(2)中位移向量t为:
步骤c2:由以上公式联立求解,得到静态CT系统的校正后的几何参数,包括:
探测器平面的投影中心点坐标(u0,v0):
u0=K13,v0=K23 (6)
射线源到探测器的距离(SDD):
SDD=K11λ (7)
探测器的扭转角θ、探测器的倾斜角φ、探测器的旋转角η分别为:
θ=Arctan 2(S32,S33)
η=Arctan2(S12,S11) (8)
射线源的位置坐标为:
o=[ox,oy,oz]T=-STt(9)
得出的位移向量t=(t1,t2,t3):
t1=A34
t2=(A24-K23A34)/K22
t3=(A14-t13A34-K2A12)/K11 (10)
射线源到扫描物体中心的距离(SOD):
优选地,所述静态CT成像系统的多个射线源各单独发射一次射线,利用发射的射线穿过所述校正体模在探测器上的各幅投影成像,可以得到整个静态CT系统的校正后的几何参数。
优选地,所述方法适合离体CT系统和活体CT系统。
优选地,所述多个射线源形成一个以上射线源阵列,所述一个以上射线源阵列和所述一个以上探测器布置为四边形或者六边形或者N边形方式。
优选地,在一个实施例中,所述校正体模的形状为立方体,所述标记小球分别安装于所述立方体的两个相对表面的角部。
有益效果
本发明的一种简单、有效、实用的静态CT系统快速几何参数校正方法,可以同时适用于离体CT系统和活体CT系统,并可以同时校正CT系统的全部10个几何参数。另外,参数校正过程较简单,仅需要各射线源对校正体模进行1次投影,即可得到CT重建所需的全部系统几何参数。
附图说明
图1是本发明一实施例的静态CT系统几何结构原理示意图。
图2是本发明的静态CT系统中的几何参数示意图。
图3是本发明一实施例的校正体膜的示意图,其中(a)是校正体模的结构示意图,(b)是该校正体模在探测器上的投影图。
图4是经过本发明的几何参数校正方法校正前后的重建图像质量比较图。
图5是将图4中所有图像和理想重建图像残差图像进行比较的示意图。
图6是本发明一实施例的静态CT系统几何参数校正方法的流程图。
具体实施方式
在下列说明中,为了提供对本发明的彻底了解而提出许多具体细节。本发明可在不具有部分或所有这些具体细节的情况下实施。在其他情况下,为了不对本发明造成不必要的混淆,不详述众所周知的过程操作。虽然本发明将结合具体实施例来进行说明,但应当理解的是,这并非旨在将本发明限制于这些实施例。
下面结合图1~6来说明本发明的静态CT系统几何参数校正方法。
本发明提供了一种静态CT系统几何参数校正方法,所述静态CT系统由多个射线源和一个以上探测器构成,所述几何参数为每个射线源与对应的探测器之间的几何参数,如图6所示,该方法包括如下步骤:
步骤a:将校正体模放置在需要校正几何参数的静态CT系统中,所述校正体模具有标记小球;
步骤b:每个射线源触发一次射线,获得多幅投影;以及
步骤c:根据每幅投影中标记小球的投影中心点的坐标,结合已知的标记小球在待校正系统中的坐标,计算获得整个静态CT系统的校正后的几何参数。
图1是本发明一实施例的的静态CT系统几何结构原理示意图。如图1所示,该静态CT系统由2个射线源阵列(各射线源阵列由11个射线源排列形成)和2个探测器构成,包括第一方向射线源阵列、第二方向射线源阵列、第一方向探测器和第二方向探测器。第一方向射线源阵列于第一方向探测器对应,第二方向射线源阵列与第二方向探测器对应。2个射线源阵列和2个探测器布置为四边形的形式。在其他实施例中,射线源阵列和探测器也可以布置为六边形或者N边形的形式。
本发明中所述几何参数为所述静态CT系统的10个几何参数,包括:射线源S的坐标位置(Ox,Oy,Oz)、射线源S在探测器平面的投影中心点坐标(u0,v0)、探测器的扭转角θ、探测器的倾斜角φ、探测器的旋转角η、射线源到扫描物体中心的距离(SOD),射线源到探测器的距离(SDD)。图2是本发明的静态CT系统中的几何参数示意图。如图2所示,S为射线源阵列(射线源阵列的坐标),定义坐标系XYZ,Y为校正体模的中心轴,(u0,v0)表示射线源S在探测器平面的投影中心点坐标。SOD为射线源到中心轴的距离,SDD为射线源到探测器的距离,U和V为探测器上像素的排列方向,定义探测器的左下角为坐标原点。
本发明中,所述校正体膜由支撑体和安装于所述支撑体上的标记小球构成,所述支撑体为低吸收材料,例如可以是碳纤维材料或者聚甲基丙烯酸甲酯,所述标记小球为高吸收材料。一般地,标记小球可以采用钢珠。本发明中,所述标记小球的数量、大小和位置不限制,以校正过程中在探测器上的投影不造成重叠为宜。
图3是本发明一实施例的校正体膜的示意图,其中(a)是校正体模的结构示意图,(b)是该校正体模在探测器上的投影图,白色圆点为标记小球的投影,单位为mm。本实施例中,所述校正体模的支撑体31的形状为立方体,所述标记小球32分别安装于所述立方体31的两个相对表面的角部。在其他实施例中,标记小球的位置不一定是在角部,位置没有限制,在其他任意位置都可以,只需要确保在实验的过程中其在待校正系统中的坐标是已知的。
本发明的几何参数校正方法中,步骤c进一步包括如下步骤c1和c2。
步骤c1:对于每幅投影图,在小球投影坐标(u,v)和实际校正体模中的标记小球的三维坐标(x,y,z)之间利用一个3×4投影矩阵A建立映射关系(该标记小球的三维坐标是已知的):
[Ax,Ay,Az,A]T=[au,av,a]T (1)
其中,α为一比例参数,投影矩阵A的构成如下:
A=K[S|t] (2)
上述(2)式中,K为一个3×3上三角矩阵,S为一个3×3旋转矩阵,t为3×1位移向量,
上述(3)式中,u0和v0为射线源在探测器平面的投影点坐标,λ为探测器像素的大小,
公式(2)中旋转矩阵S可以进一步用三角函数表达为:
θ为探测器的扭转角、φ为探测器的倾斜角、η为探测器的旋转角,
公式(2)中位移向量t为:
步骤c2:由以上公式(1)~(5)联立求解,可以得到静态CT系统的校正后的几何参数,包括:
探测器平面的投影中心点坐标(u0,v0):
u0=K13,v0=K23 (6)
射线源到探测器的距离(SDD):
SDD=K11λ (7)
探测器的扭转角θ、探测器的倾斜角φ、探测器的旋转角η分别为:
θ=Arc tan 2(S32,S33)
η=Arc tan2(S12,S11) (8)
射线源的位置坐标为:
o=[ox,oy,oz]T=-STt (9)
得出的位移向量t=(t1,t2,t3):
t1=A34
t2=(A24-K23A34)/K22
t3=(A14-t13A34-K2A12)/K11 (10)
射线源到扫描物体中心的距离(SOD):
在静态CT成像系统具有多个射线源的情况下,静态CT成像系统的多个射线源各单独发射一次射线,利用发射的射线穿过校正体模在探测器上的各幅成像,就可以得到整个静态CT系统的校正后的几何参数。获得的校正后的几何参数可用于系统机械安装误差的检测及修正,也可用于系统CT高质量成像重建。
本发明的静态CT系统几何参数校正方法适合离体CT系统和活体CT系统。离体CT系统是样本旋转,射线源和探测器不动的系统,活体CT系统是样本不动,射线源和探测器旋转的系统。另外,本发明的方法除应用于锥束CT系统,经适当变形后,也可应用于PET、SPECT系统。
图4是经过本发明的几何参数校正方法校正前后的重建图像质量比较图。(a)为理想无误差图像。(b)~(d)为不同探测器偏移情况下未经过参数校正直接重建的图像质量:其中,(b)偏移offset=1pixel(像素);(c)偏移offset=2pixels;(d)偏移offset=5pixels。(e)~(h)为不同探测器偏移情况下经过几何参数校正后再重建的图像质量:其中,(e)偏移offset=0pixel;(f)偏移offset=1pixel;(g)偏移offset=2pixels;(h)偏移offset=5pixels。从图中可以看到,经过几何参数校正后进行重建的图像成像质量较好,图像中的伪影远远小于没有经过参数校正的图像。
图5是将图4中所有图像和理想重建图像之间的残差图像进行比较的示意图。所谓残差图像即将图4中所有图像与图4(a)分别相减后的残余图像。将图4中所有图像和理想重建图像之间的残差图像进行比较,从图中可以看到,经过几何参数校正后进行重建的图像的残差远远小于没有经过几何参数校正的图像。这里的图像窗宽为[-0.13 0.07]。这也从另一方面证明,静态CT系统经过几何参数校正后进行重建的图像成像质量较好。
显然,本技术领域中的普通技术人员应当认识到,以上的实施例仅是用来说明本发明,而并非用作为对本发明的限定,只要在本发明的实质精神范围内,对以上所述实施例的变化、变型都将落在本发明的权利要求书范围内。
Claims (10)
1.一种静态CT系统几何参数校正方法,所述静态CT系统由多个射线源和一个以上探测器构成,所述几何参数为每个射线源与对应的探测器之间的几何参数,其特征在于,所述方法包括如下步骤:
步骤a:将校正体模放置在需要校正几何参数的静态CT系统中,所述校正体模具有标记小球;
步骤b:每个射线源触发一次射线,获得多幅投影;以及
步骤c:根据每幅投影中标记小球的投影中心点的坐标,结合已知的标记小球在待校正系统中的坐标,计算获得整个静态CT系统的校正后的几何参数。
2.根据权利要求1所述的静态CT系统几何参数校正方法,其特征在于,所述几何参数为所述静态CT系统的10个几何参数,包括:射线源S的坐标位置(Ox,Oy,Oz)、射线源S在探测器平面的投影中心点坐标(u0,v0)、探测器的扭转角θ、探测器的倾斜角φ、探测器的旋转角η、射线源到扫描物体中心的距离(SOD),射线源到探测器的距离(SDD)。
3.根据权利要求1所述的静态CT系统几何参数校正方法,其特征在于,所述校正体模由支撑体和安装于所述支撑体上的标记小球构成,所述支撑体为低吸收材料,所述标记小球为高吸收材料。
4.根据权利要求3所述的静态CT系统几何参数校正方法,其特征在于,所述支撑体为碳纤维材料或者聚甲基丙烯酸甲酯,所述标记小球为钢珠。
5.根据权利要求4所述的静态CT系统几何参数校正方法,其特征在于,所述标记小球的数量、大小以及位置以校正过程中在探测器上的投影不造成重叠为宜。
6.根据权利要求5所述的静态CT系统几何参数校正方法,其特征在于,所述步骤c进一步包括如下步骤:
步骤c1:对于每幅投影图,在小球投影坐标(u,v)和实际校正体模中的标记小球的三维坐标(x,y,z)之间利用一个3×4投影矩阵A建立映射关系:
[Ax,Ay,Az,A]T=[au,av,a]T (1)
其中,α为一比例参数,投影矩阵A的构成如下:
A=K[S|t] (2)
上述(2)式中,K为一个3×3上三角矩阵,S为一个3×3旋转矩阵,t为3×1位移向量,
上述(3)式中,u0和v0为射线源在探测器平面的投影点坐标,λ为探测器像素的大小,
公式(2)中旋转矩阵S可以进一步用三角函数表达为:
θ为探测器的扭转角、φ为探测器的倾斜角、η为探测器的旋转角,公式(2)中位移向量t为:
步骤c2:由以上公式联立求解,得到静态CT系统的校正后的几何参数,包括:
探测器平面的投影中心点坐标(u0,v0):
u0=K13,v0=K23 (6)
射线源到探测器的距离(SDD):
SDD=K11λ (7)
探测器的扭转角θ、探测器的倾斜角φ、探测器的旋转角η分别为:
θ=Arc tan 2(S32,S33)
η=Arc tan 2(S12,S11) (8)
射线源的位置坐标为:
o=[ox,oy,oz]T=-STt (9)
得出的位移向量t=(t1,t2,t3):
t1=A34
t2=(A24-K23A34)/K22
t3=(A14-t13A34-K2A12)/K11 (10)
射线源到扫描物体中心的距离(SOD):
7.根据权利要求1所述的静态CT系统几何参数校正方法,其特征在于,所述静态CT成像系统的多个射线源各单独发射一次射线,利用发射的射线穿过所述校正体模在探测器上的各幅投影成像,可以得到整个静态CT系统的校正后的几何参数。
8.根据权利要求1所述的静态CT系统几何参数校正方法,其特征在于,所述方法适合离体CT系统和活体CT系统。
9.根据权利要求1所述的静态CT系统几何参数校正方法,其特征在于,所述多个射线源形成一个以上射线源阵列,所述一个以上射线源阵列和所述一个以上探测器布置为四边形或者六边形或者N边形方式。
10.根据权利要求3所述的静态CT系统几何参数校正方法,其特征在于,在一个实施例中,所述校正体模的形状为立方体,所述标记小球分别安装于所述立方体的两个相对表面的角部。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711165876.1A CN108201447A (zh) | 2017-11-21 | 2017-11-21 | 一种静态ct系统几何参数校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711165876.1A CN108201447A (zh) | 2017-11-21 | 2017-11-21 | 一种静态ct系统几何参数校正方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108201447A true CN108201447A (zh) | 2018-06-26 |
Family
ID=62603810
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711165876.1A Pending CN108201447A (zh) | 2017-11-21 | 2017-11-21 | 一种静态ct系统几何参数校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108201447A (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108903964A (zh) * | 2018-07-09 | 2018-11-30 | 广州华端科技有限公司 | 计算机断层图像的散射校正方法以及装置 |
CN109452950A (zh) * | 2018-12-29 | 2019-03-12 | 上海联影医疗科技有限公司 | 一种模体及成像设备的几何标记方法 |
CN109658465A (zh) * | 2018-12-07 | 2019-04-19 | 广州华端科技有限公司 | 图像重建过程中的数据处理、图像重建方法和装置 |
CN109671128A (zh) * | 2018-12-07 | 2019-04-23 | 广州华端科技有限公司 | 图像重建过程中的数据处理、图像重建方法和装置 |
CN109875591A (zh) * | 2019-03-14 | 2019-06-14 | 浙江大学 | 一种锥束ct系统的几何标定方法 |
CN110353713A (zh) * | 2019-07-16 | 2019-10-22 | 上海联影医疗科技有限公司 | 几何校正模体、重建图像几何参数的确定方法和装置 |
CN110766629A (zh) * | 2019-10-17 | 2020-02-07 | 广州华端科技有限公司 | Cbct系统几何校正方法、装置、计算机设备和存储介质 |
CN111178457A (zh) * | 2020-04-09 | 2020-05-19 | 南京佗道医疗科技有限公司 | 一种精度校验器及三维影像重建精度评价方法 |
CN111265231A (zh) * | 2019-04-15 | 2020-06-12 | 清华大学 | 分布式光源ct图像重建方法与系统 |
CN111419257A (zh) * | 2020-04-07 | 2020-07-17 | 南京安科医疗科技有限公司 | 校准模体自动定位方法、扫描系统和可读存储介质 |
CN112294349A (zh) * | 2020-10-29 | 2021-02-02 | 上海联影医疗科技股份有限公司 | X射线设备校正方法、成像方法和成像系统 |
CN113643428A (zh) * | 2021-08-17 | 2021-11-12 | 北京唯迈医疗设备有限公司 | 一种适用于多自由度锥形束ct的全参数几何校准方法 |
CN113892960A (zh) * | 2021-10-09 | 2022-01-07 | 清华大学 | X射线自成像几何标定方法及装置 |
CN113939229A (zh) * | 2019-04-10 | 2022-01-14 | 皇家飞利浦有限公司 | 静态增益校准 |
CN115077436A (zh) * | 2022-07-21 | 2022-09-20 | 浙江双元科技股份有限公司 | 一种叠片电池的极片错位检测系统及方法 |
WO2022213274A1 (zh) * | 2021-04-06 | 2022-10-13 | 中国科学技术大学 | 通用的ct对轴方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5442674A (en) * | 1993-01-27 | 1995-08-15 | Ge Medical Systems | Device and automatic method for the geometrical calibration of an X-ray imaging system |
CN103006251A (zh) * | 2012-12-06 | 2013-04-03 | 深圳先进技术研究院 | 用于ct系统中标定几何参数的标定体模、标定装置及标定方法 |
CN103462630A (zh) * | 2013-09-13 | 2013-12-25 | 深圳先进技术研究院 | Ct系统及ct扫描方法 |
CN103549971A (zh) * | 2013-11-07 | 2014-02-05 | 北京航空航天大学 | 一种确定c型臂断层成像系统中几何标定参数的方法 |
CN104783824A (zh) * | 2014-01-20 | 2015-07-22 | 上海联影医疗科技有限公司 | X射线成像系统的校正方法 |
CN105931202A (zh) * | 2016-04-20 | 2016-09-07 | 广州华端科技有限公司 | 几何校正模体的校正方法和系统 |
CN105997126A (zh) * | 2016-05-25 | 2016-10-12 | 重庆大学 | 一种锥束ct系统几何参数校正模型及方法 |
CN106153648A (zh) * | 2016-06-22 | 2016-11-23 | 深圳先进技术研究院 | 一种静态高分辨显微ct成像系统及应用其的成像方法 |
CN106388848A (zh) * | 2016-10-18 | 2017-02-15 | 深圳先进技术研究院 | 对ct图像进行预处理的方法、系统及静态ct成像装置 |
CN106491152A (zh) * | 2016-10-18 | 2017-03-15 | 深圳先进技术研究院 | 一种基于光子探测器的静态ct成像系统 |
-
2017
- 2017-11-21 CN CN201711165876.1A patent/CN108201447A/zh active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5442674A (en) * | 1993-01-27 | 1995-08-15 | Ge Medical Systems | Device and automatic method for the geometrical calibration of an X-ray imaging system |
CN103006251A (zh) * | 2012-12-06 | 2013-04-03 | 深圳先进技术研究院 | 用于ct系统中标定几何参数的标定体模、标定装置及标定方法 |
CN103462630A (zh) * | 2013-09-13 | 2013-12-25 | 深圳先进技术研究院 | Ct系统及ct扫描方法 |
CN103549971A (zh) * | 2013-11-07 | 2014-02-05 | 北京航空航天大学 | 一种确定c型臂断层成像系统中几何标定参数的方法 |
CN104783824A (zh) * | 2014-01-20 | 2015-07-22 | 上海联影医疗科技有限公司 | X射线成像系统的校正方法 |
CN105931202A (zh) * | 2016-04-20 | 2016-09-07 | 广州华端科技有限公司 | 几何校正模体的校正方法和系统 |
CN105997126A (zh) * | 2016-05-25 | 2016-10-12 | 重庆大学 | 一种锥束ct系统几何参数校正模型及方法 |
CN106153648A (zh) * | 2016-06-22 | 2016-11-23 | 深圳先进技术研究院 | 一种静态高分辨显微ct成像系统及应用其的成像方法 |
CN106388848A (zh) * | 2016-10-18 | 2017-02-15 | 深圳先进技术研究院 | 对ct图像进行预处理的方法、系统及静态ct成像装置 |
CN106491152A (zh) * | 2016-10-18 | 2017-03-15 | 深圳先进技术研究院 | 一种基于光子探测器的静态ct成像系统 |
Non-Patent Citations (2)
Title |
---|
XINHUA LI,ET AL: "A generic geometric calibration method for tomographic imaging systems with flat-panel detectors—A detailed implementation guide", 《MEDICAL PHYSICS》 * |
YOUNGBIN CHO,ET AL: "Accurate technique for complete geometric calibration of cone-beam computed tomography systems", 《MEDICAL PHYSICS》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108903964A (zh) * | 2018-07-09 | 2018-11-30 | 广州华端科技有限公司 | 计算机断层图像的散射校正方法以及装置 |
CN109658465A (zh) * | 2018-12-07 | 2019-04-19 | 广州华端科技有限公司 | 图像重建过程中的数据处理、图像重建方法和装置 |
CN109671128A (zh) * | 2018-12-07 | 2019-04-23 | 广州华端科技有限公司 | 图像重建过程中的数据处理、图像重建方法和装置 |
CN109452950A (zh) * | 2018-12-29 | 2019-03-12 | 上海联影医疗科技有限公司 | 一种模体及成像设备的几何标记方法 |
CN109875591B (zh) * | 2019-03-14 | 2020-05-29 | 浙江大学 | 一种锥束ct系统的几何标定方法 |
CN109875591A (zh) * | 2019-03-14 | 2019-06-14 | 浙江大学 | 一种锥束ct系统的几何标定方法 |
CN113939229A (zh) * | 2019-04-10 | 2022-01-14 | 皇家飞利浦有限公司 | 静态增益校准 |
CN111265231A (zh) * | 2019-04-15 | 2020-06-12 | 清华大学 | 分布式光源ct图像重建方法与系统 |
CN111265231B (zh) * | 2019-04-15 | 2021-08-31 | 清华大学 | 分布式光源ct图像重建方法与系统 |
CN110353713A (zh) * | 2019-07-16 | 2019-10-22 | 上海联影医疗科技有限公司 | 几何校正模体、重建图像几何参数的确定方法和装置 |
CN110766629A (zh) * | 2019-10-17 | 2020-02-07 | 广州华端科技有限公司 | Cbct系统几何校正方法、装置、计算机设备和存储介质 |
CN110766629B (zh) * | 2019-10-17 | 2022-03-01 | 广州华端科技有限公司 | Cbct系统几何校正方法、装置、计算机设备和存储介质 |
CN111419257B (zh) * | 2020-04-07 | 2022-11-22 | 南京安科医疗科技有限公司 | 校准模体自动定位方法、扫描系统和可读存储介质 |
CN111419257A (zh) * | 2020-04-07 | 2020-07-17 | 南京安科医疗科技有限公司 | 校准模体自动定位方法、扫描系统和可读存储介质 |
CN111178457B (zh) * | 2020-04-09 | 2020-07-03 | 南京佗道医疗科技有限公司 | 一种精度校验器及三维影像重建精度评价方法 |
CN111178457A (zh) * | 2020-04-09 | 2020-05-19 | 南京佗道医疗科技有限公司 | 一种精度校验器及三维影像重建精度评价方法 |
CN112294349A (zh) * | 2020-10-29 | 2021-02-02 | 上海联影医疗科技股份有限公司 | X射线设备校正方法、成像方法和成像系统 |
WO2022213274A1 (zh) * | 2021-04-06 | 2022-10-13 | 中国科学技术大学 | 通用的ct对轴方法 |
CN113643428A (zh) * | 2021-08-17 | 2021-11-12 | 北京唯迈医疗设备有限公司 | 一种适用于多自由度锥形束ct的全参数几何校准方法 |
CN113892960A (zh) * | 2021-10-09 | 2022-01-07 | 清华大学 | X射线自成像几何标定方法及装置 |
CN115077436A (zh) * | 2022-07-21 | 2022-09-20 | 浙江双元科技股份有限公司 | 一种叠片电池的极片错位检测系统及方法 |
CN115077436B (zh) * | 2022-07-21 | 2022-11-22 | 浙江双元科技股份有限公司 | 一种叠片电池的极片错位检测系统及方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108201447A (zh) | 一种静态ct系统几何参数校正方法 | |
US7950849B2 (en) | Method and device for geometry analysis and calibration of volumetric imaging systems | |
US7307252B2 (en) | Detector head position correction for hybrid SPECT/CT imaging apparatus | |
US7424089B2 (en) | System and method for reconstructing image by using straight-line trajectory scan | |
Stevens et al. | Alignment of a volumetric tomography system | |
CN104107065B (zh) | 3d图像集在不同空间之间的最佳变换 | |
US20050265523A1 (en) | C-arm device with adjustable detector offset for cone beam imaging involving partial circle scan trajectories | |
US7103233B2 (en) | Methods and apparatus for determining component alignment | |
JP2003019128A (ja) | イメージング・システムを較正するシステム及び方法、並びに記憶媒体 | |
US9336569B2 (en) | Image compensation method and apparatus | |
US20100208969A1 (en) | Rotation center identifying method, ring artifact correction method, rotation center identifying apparatus, x-ray diagnostic apparatus, recording medium on which program for executing rotation center identification is recorded and recording medium on which program for executing ring artifact correction is recorded | |
US20150071515A1 (en) | Image reconstruction method and device for tilted helical scan | |
CN107157505B (zh) | 一种提高锥形束spect成像分辨率的方法 | |
CN105319225B (zh) | 一种实现板状样品高分辨率大视野cl成像的扫描方法 | |
CN109658465B (zh) | 图像重建过程中的数据处理、图像重建方法和装置 | |
US8023713B2 (en) | Method and system for reducing artifact due to time delay in data acquisition system in computer tomography | |
Ferrucci et al. | Measurement of sample stage error motions in cone-beam x-ray computed tomography instruments by minimization of reprojection errors | |
CN101447085B (zh) | Micro CT 系统中对投影图像旋转中心偏移的校正方法 | |
CN109146800A (zh) | 锥束计算机断层成像图像校正方法和系统 | |
CN111839568A (zh) | 一种新型大视场直线扫描ct系统及图像重建方法 | |
CN103606144A (zh) | Ct投影旋转中心的确定方法和装置 | |
US11175242B2 (en) | Geometric alignment, sample motion correction, and intensity normalization of computed tomography projections using pi-line optimization | |
Duan et al. | Knowledge-based self-calibration method of calibration phantom by and for accurate robot-based CT imaging systems | |
CN112603346B (zh) | 一种基于标记物成像的探测器偏转校正方法 | |
CN108364325A (zh) | 规则样品x射线ct投影图位置平移偏差检测及校正方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180626 |