CN202049120U - 一种消除ct图像中的几何伪影的系统 - Google Patents
一种消除ct图像中的几何伪影的系统 Download PDFInfo
- Publication number
- CN202049120U CN202049120U CN2011200556334U CN201120055633U CN202049120U CN 202049120 U CN202049120 U CN 202049120U CN 2011200556334 U CN2011200556334 U CN 2011200556334U CN 201120055633 U CN201120055633 U CN 201120055633U CN 202049120 U CN202049120 U CN 202049120U
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- math
- mover
- 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.)
- Expired - Fee Related
Links
- 238000002591 computed tomography Methods 0.000 title claims abstract description 158
- 238000012545 processing Methods 0.000 claims abstract description 38
- 238000000034 method Methods 0.000 claims description 33
- 230000014509 gene expression Effects 0.000 claims description 29
- 238000003672 processing method Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 description 8
- 238000012937 correction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 239000008188 pellet Substances 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
Images
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本实用新型涉及一种消除CT图像中的几何伪影的系统。CT包括发射X射线的X射线源;带动物体旋转的转台;探测到达自身的X射线强度的X射线探测器;处理X射线探测器所探测的X射线强度从而得到待成像物体的扇束投影数据的处理装置;该系统进一步包括模体;处理装置根据得到的模体质心在X射线探测器上的投影坐标确定X射线源的焦点、转台的旋转中心、X射线探测器间相对位置的几何参数,将其代入含有几何参数的重建公式,利用重建公式来处理扇束投影数据从而得到待成像物体不包含几何伪影的CT图像。利用本实用新型的技术方案,能精确得到没有几何伪影的CT图像。
Description
技术领域
本实用新型涉及CT应用领域,特别是涉及一种消除CT图像中的几何伪影的系统。
背景技术
CT(X射线计算机断层成像,computed tomography)是一种重要的无损检测技术,它的工作原理为:利用不同物体对X射线的吸收和透过率的不同,用X射线源发射X射线穿透物体,然后用灵敏度极高的探测器对穿透物体后的X射线强度进行检测,由处理装置对获取的强度数据进行处理,就可以得到物体的断面或立体图像,通过观察该图像,就可以了解物体内部的结构。在医学领域,利用CT对人体进行检测,就可以得到人体受检查部位的断面或立体图像,进而发现体内任何部位的细小病变。
但是,作为一种利用X射线进行检测的仪器,其各构件之间的位置关系必须精确满足一定的条件,例如,为了防止出现由CT构件的几何参数误差所引起的CT图像的几何伪影,必须使CT构件满足理想的几何成像关系,即支撑待成像物体的转轴必须与探测器所在直线平行,该转轴在探测器上的投影点必须位于探测器的中心,X射线源和转轴所在的平面必须与探测器所在平面垂直,另外,在检测中还需要精确已知X射线源焦点到待成像物体的质心的距离以及X射线源焦点到探测器的距离。然而实际使用的CT很难保证满足上述条件,由此产生的几何参数误差导致用CT得到的图像有几何伪影,这大大降低了CT图像的分辨率,也影响了对CT图像的正确判读,进而影响了对病情的判断。
现有技术是用手动方式调节CT中的某些构件,依靠对CT图像质量的主观判断来消除CT图像中的几何伪影的,这种校正方式比较粗糙,主要依靠调节人员的经验来判断校正效果,因而校正质量没有保证,还可能会发生因为调节失误而误差加大且无法复原的情况。
实用新型内容
本实用新型所要解决的技术问题是提供一种消除CT图像中的几何伪影的系统,能精确得到没有几何伪影的CT图像。
本实用新型解决上述技术问题的技术方案如下:一种消除CT图像中的几何伪影的系统,所述CT包括用于发射X射线的X射线源、带动待成像物体旋转的转台、探测到达自身的X射线强度的X射线探测器、对所述X射线探测器所探测到的X射线的强度进行处理从而得到所述待成像物体的扇束投影数据的处理装置,该系统进一步包括模体;
所述X射线源用于,发射对模体进行CT扫描的X射线;
所述转台用于,带动模体旋转;
所述处理装置用于,对所述X射线探测器探测到的X射线的强度进行处理,得到所述模体的质心在所述X射线探测器上的投影坐标;根据所述模体的质心在X射线探测器上的投影坐标,确定几何参数,并将所述几何参数代入重建公式,其中,所述几何参数为用于确定X射线源焦点、转台旋转中心、X射线探测器之间相对位置的几何参数;利用所述重建公式处理所述待成像物体的扇束投影数据,得到所述待成像物体的CT图像数据,则用所述待成像物体的CT图像数据所得到的CT图像中没有几何伪影。
本实用新型的有益效果是:本实用新型中,由于处理装置用CT扫描的方式得到模体的质心在X射线探测器上的投影坐标,处理装置根据该投影坐标即可确定出几何参数,这些几何参数可以用于确定CT中的X射线源焦点、转台旋转中心、X射线探测器之间的相对位置,这样,将这些几何参数代入重建公式,然后用重建公式对待成像物体的扇束投影数据进行处理,用所得到的待成像物体的CT图像数据作为CT图像的灰度,即可精确消除CT图像中的几何伪影,得到没有几何伪影的CT图像;另外,该方法不需人工参与,因而也消除了人工调节误差的存在以及人工误差无法复原的可能性。
在上述技术方案的基础上,本实用新型还可以做如下改进:
进一步,所述处理装置用于,对所述X射线探测器探测到的X射线的强度进行处理,得到所述模体的投影数据;用图像处理方法从所述模体的投影数据中得到所述模体的质心在X射线探测器上的投影坐标。
进一步,所述模体的质心在X射线探测器上的投影坐标为:所述模体的质心在X射线探测器上n个不同角度的投影坐标,其中,n为正整数;
所述模体的质心在X射线探测器上第k个角度处的投影坐标为uk,k为不小于零的整数,所述第k个角度为γk=β+45*k度,β为0度至45度之间的任一角度。
进一步,所述几何参数包括:转台旋转中心在X射线探测器上的投影点与X射线探测器的中心点之间的距离h、中心射线与X射线探测器所在直线的夹角α、转台旋转中心在X射线探测器上的投影点与X射线源焦点之间的距离SDD,其中,所述中心射线为穿过所述转台旋转中心的X射线。
进一步,所述重建公式为:
其中,f(r0,θ0)为所述重建图像中极坐标为(r0,θ0)的点的重建图像数据;p(r,θ)为所述待成像物体上极坐标为(r,θ)的点的扇束投影数据;
进一步,h为:利用以下四个表达式中的任一个计算得到的h,或,利用以下四个表达式中的任意两个以上所计算得到的h的平均值,所述四个表达式为:
表达式一:
表达式二:
表达式三:
表达式四:
进一步, 其中,
或,
进一步,所述几何参数进一步包括转台旋转中心与X射线源焦点之间的距离SOD。
进一步,所述模体包括第一个模体和第二个模体;则
进一步,第一个模体的质心和第二个模体的质心在X射线探测器上第k个角度处的投影坐标分别为u0k和u1k,所述第k个角度为γk=β+45*k度,其中,k为0与n之间的所有整数,n为正整数,β为0度至45度之间的任一角度;则和可利用下列任一组关系式得到:
第一组关系式:
第二组关系式:
第三组关系式:
第四组关系式:
进一步,所述模体为柱形物体或球形物体。
附图说明
图1为利用本实用新型消除CT图像中的几何伪影的方法流程图;
图2为本实用新型提供的消除CT图像中的几何伪影的系统结构图;
图3为模体为细丝时利用各点投影数据得到的细丝的CT图像,其中,图3a为其整体图像,图3b为图3a中方框部分的图像;
图4为应用本实用新型的CT各构件的位置关系图;
图5为利用现有技术对CT进行校正后扫描两个易拉罐模体得到的CT图像,其中,图5a为整体图像,图5b和图5c分别为图5a中曲线501和502上的方框部分的放大图;
图6为假设其他构件的位置正确,手工调节模体质心在X射线探测器上的投影点与X射线探测器的中心点之间距离得到的两个易拉罐模体的CT图像,其中,图6a为整体图像,图6b和图6c分别为图6a中曲线601和602方框部分的放大图;
图7为利用本实用新型对CT进行校正后扫描两个易拉罐模体得到的CT图像,其中,图7a为整体图像,图7b和图7c分别为图7a中曲线701和702上的方框部分的放大图。
具体实施方式
以下结合附图对本实用新型的原理和特征进行描述,所举实例只用于解释本实用新型,并非用于限定本实用新型的范围。
图1为利用本实用新型消除CT图像中的几何伪影的方法流程图。如图1所示,该方法包括:
步骤101:对模体进行CT扫描,获得模体的质心在X射线探测器上的投影坐标。
这里,CT包括X射线源、转台和X射线探测器,其中,X射线源可发射穿过模体的X射线,转台带动待成像物体旋转,X射线探测器可精确检测穿过模体到达自身的X射线的强度。本步骤中用X射线源发射的X射线对模体进行CT扫描,以及步骤103中对待成像物体进行CT扫描,都是通过X射线探测器与模体或待成像物体之间的相对运动,利用X射线来实现对模体或待成像物体的扫描。
本步骤中对模体进行CT扫描后,X射线探测器得到的穿过模体的X射线强度数据经过CT中的处理装置的处理,就成为模体的投影数据。从这些投影数据中可以分辨出模体的质心在X射线探测器上的投影坐标。当然,也可以从这些投影数据中分辨出模体上质心之外的其他点在X射线探测器上的投影坐标,例如,在模体为细丝的情况下,本实用新型可以获得细丝的质心在X射线探测器上的投影坐标,也可以获得细丝的某一个端点在X射线探测器上的投影坐标。
这里所说的模体为结构简单、易于判断质心的物体,它可以为柱形物体,如细丝等,也可以为球形物体,如小球等,当然,模体也可以为其他物体。
模体位于CT中的转台上,转台具有旋转中心。
步骤102:根据模体的质心在X射线探测器上的投影坐标,确定几何参数,并将几何参数代入重建公式,其中,几何参数为用于确定CT中的X射线源焦点、转台旋转中心、X射线探测器之间相对位置的几何参数。
这里,几何伪影的产生,主要是由于CT中的X射线源焦点、转台旋转中心、X射线探测器之间的相对位置不能精确满足预定的位置条件,具有位置参数偏差。而本步骤中的几何参数是指可以确定CT中的X射线源焦点、转台旋转中心、X射线探测器之间相对位置的几何参数,即本实用新型用几何参数来表征X射线源焦点、转台旋转中心、X射线探测器之间相对位置的偏差,将几何参数代入重建公式,然后用重建公式对待成像物体的扇束投影数据进行校正,即可得到消除了几何伪影的CT图像中。
几何参数可以为一个参数,也可以为多个参数,例如,几何参数可以包括:转台旋转中心在X射线探测器上的投影点与X射线探测器的中心点之间的距离h、中心射线与X射线探测器所在直线的夹角α、转台旋转中心在X射线探测器上的投影点与X射线源焦点之间的距离SDD,其中,中心射线为穿过转台旋转中心的X射线,还可以进一步包括转台旋转中心与X射线源焦点之间的距离SOD,这几个几何参数对于利用重建公式消除CT图像中的几何伪影所发挥的作用是相互独立的。当然,几何参数还可以为其他形式,例如,可以为由上述h、α、SDD和SOD推导出的其他几何参数,这些几何参数对于消除CT图像中的几何伪影所发挥的作用可以是相互独立的,当然,各几何参数之间也可以是不相互独立的。
几何参数是根据步骤101得到的模体的质心在X射线探测器上的投影坐标来确定的,将几何参数代入重建公式之后,就可以利用重建公式来消除CT图像中的几何伪影了。
步骤103:对待成像物体进行CT扫描,得到待成像物体的扇束投影数据;利用重建公式处理待成像物体的扇束投影数据,得到待成像物体的CT图像数据,则利用待成像物体的CT图像数据所得到的CT图像中没有几何伪影。
这里,待成像物体的扇束投影数据是由具有几何位置偏差的CT对待成像物体进行扫描后得到的,用步骤102得到的代入了几何参数的重建公式对扇束投影数据进行处理,即可得到待成像物体的CT图像数据,然后利用待成像物体的CT图像数据可以得到没有几何伪影的CT图像,这是用数学的方式来精确消除CT图像中的几何伪影,由此可见,本实用新型是用间接的方式来消除几何伪影的,并没有采用现有技术直接消除几何伪影的方式。
利用重建公式处理待成像物体的扇束投影数据,该处理的方式可以为将待成像物体的扇束投影数据代入重建公式,从而得到待成像物体的CT图像数据,然后利用这些CT图像数据来得到CT图像,这里,利用CT图像数据来得到CT图像的方法随着所要得到的CT图像的不同而有所不同,例如,要得到灰度图,就可以对CT图像数据进行量化等处理,要得到其他种类的CT图像,例如高度图,则要进行其他形式的处理,不管得到何种CT图像,只要利用的是重建公式处理后所得到的CT图像数据,这样形成的CT图像中都不包含几何伪影。
该方法利用步骤101和102得到代入了几何参数的重建公式之后,可以保持该重建公式不变,稳定工作一段时间,即多次利用该重建公式进行步骤103所述的重建运算,而无需在每次执行步骤103之前都重新进行步骤101和102。这样,该方法可以在计算一次几何参数的基础上进行多次CT图像重建,大大节约了CT扫描的人力、物力和时间。
由此可见,本实用新型中,由于用CT扫描的方式得到模体的质心在X射线探测器上的投影坐标,根据该投影坐标即可确定出几何参数,这些几何参数可以用于确定CT中的X射线源焦点、转台旋转中心、X射线探测器之间的相对位置,这样,将这些参数代入重建公式,然后用重建公式对待成像物体的扇束投影数据进行处理,用所得到的待成像物体的CT图像数据作为CT图像的灰度,即可精确消除CT图像中的几何伪影,得到没有几何伪影的CT图像;另外,该方法不需人工参与,因而也消除了人工调节误差的存在以及人工误差无法复原的可能性。
步骤101中,对模体进行CT扫描,获得模体的质心在X射线探测器上的投影坐标的方法为:对模体进行CT扫描,得到模体的投影数据;用图像处理方法从模体的投影数据中获得模体的质心在X射线探测器上的投影坐标。
这里,从模体的投影数据中获得模体的质心在X射线探测器上的投影坐标的图像处理方法多种多样,可以为质心法,也可以为二值法、阈值法等。
模体的质心在X射线探测器上的投影坐标可以包括:模体的质心在X射线探测器上n个不同角度的投影坐标,其中,n为正整数,模体的质心在X射线探测器上第k个角度处的投影坐标为uk,k为不小于零的整数,第k个角度为γk=β+45*k度,β为0度至45度之间的任一角度。
下面均以n等于8为例进行说明,n为大于8的整数的情况可以认为是在n等于8的情况下进行的延伸,增加了(n-8)个角度处的投影坐标,从而使本实用新型对几何伪影的校正效果更好。n为小于8的正整数的情况也可以类推。
在n等于8的情况下,模体的质心在X射线探测器上第k个角度处的投影坐标为uk,第k个角度为γk=β+45*k度,其中,k为0与7之间的所有整数,β为0度至45度之间的任一角度。
由此可见,这8个不同的角度是具有一定的对应关系的,由于相邻角度之间相差45度,因此,这8个角度可分为四组,每组包括两个相差180度的角度,而且由于β为0度至45度之间的任一角度,因此,第k个角度为γk=β+45*k度意味着这8个角度为0度至360度中满足相邻角度差为45度且组内角度差为180度条件的任意四组角度。
uk可以通过下式得到:
式中的Maxnum为X射线探测器的通道总数;pi,k为X射线探测器的第i个通道所测得的第k个角度处的投影坐标。
步骤102所确定的几何参数有多种表示形式,例如,这些几何参数可以包括:转台旋转中心在X射线探测器上的投影点与X射线探测器的中心点之间的距离h、中心射线与X射线探测器所在直线的夹角α、转台旋转中心在X射线探测器上的投影点与X射线源焦点之间的距离SDD,其中,中心射线为穿过转台旋转中心的X射线。
这样,步骤102所述的重建公式为:
其中,f(r0,θ0)为重建图像中极坐标为(r0,θ0)的点的重建图像数据;p(r,θ)为待成像物体上极坐标为(r,θ)的点的扇束投影数据;
由此可见,该重建公式是利用h、α和SDD这三个几何参数,对CT得到的待成像物体上极坐标为(r,θ)的点的扇束投影数据p(r,θ)进行处理,得到CT图像数据,从而修正由几何位置误差所引起的几何伪影,这样得到的以CT图像数据为灰度的CT图像就没有几何伪影了。
这些几何参数中,h为利用以下四个表达式中的任一个计算得到的h,或为利用以下四个表达式中的任意两个以上所计算得到的h的平均值,这四个表达式为:
表达式一:
表达式二:
表达式三:
表达式四:
α和SDD可以分别通过以下方法得到:
其中,n1和n2可以通过以下两种方法之一来确定:
第一种方法为:
第二种方法为:
这里的几何参数还可以为其他形式,例如,几何参数包括:转台旋转中心在X射线探测器上的投影点与X射线探测器的中心点之间的距离h以及参数n1、n2,其中,
h为利用以下四个解析式中的任一个计算得到的h,或为利用以下四个解析式中的任意两个以上所计算得到的h的平均值,这四个解析式为:
解析式一:
解析式二:
解析式三:
解析式四:
n1和n2为利用以下两组解析式中的任一组计算得到的n1和n2,这两组解析式为:
第一组解析式:
第二组解析式:
这种情况下,重建公式可以为:
其中,f(x0,y0)为重建图像中坐标为(x0,y0)的点的重建图像数据;p(x,y)为待成像物体上坐标为(x,y)的点的扇束投影数据;
可见,该重建公式是利用h、n1和n2这三个几何参数对待成像物体上坐标为(x,y)的点的扇束投影数据p(x,y)进行处理,得到CT图像数据,从而修正由CT构件的几何位置误差所引起的几何伪影,这样得到的以CT图像数据为灰度的CT图像就没有几何伪影了。
由以上表达式可以看出,h、α和SDD这三个几何参数与h、n1和n2这三个几何参数之间是可以相互推导得到的,即这两组几何参数在本质上是同一组参数。另外,利用这两组几何参数中的任一组还可以推导出其他形式的几何参数,均在本实用新型的保护范围之内。
另外,步骤102确定的几何参数进一步包括转台旋转中心与X射线源焦点之间的距离SOD。SOD不会对重建图像的质量造成影响,即SOD不会对消除CT图像中的几何伪影提供帮助,而只与标定CT图像的实际尺寸有关,即SOD不同,则最终在步骤103中得到的CT图像的尺寸是不同的。
SOD可以利用两个模体用以下方法来确定:
这里,模体包括第一个模体和第二个模体;则SOD通过以下方法得到:
设第一个模体的质心和第二个模体的质心在X射线探测器上第k个角度处的投影坐标分别为u0k和u1k,第k个角度为γk=β+45*k度,其中,k为0与n之间的所有整数,n为正整数,β为0度至45度之间的任一角度;则和可利用下列任一组关系式得到:
第一组关系式:
第二组关系式:
第三组关系式:
第四组关系式:
这里的u0k和u1k可以分别用上面所述的uk的计算公式得到。
本实用新型采用间接方法来消除CT图像中的几何伪影,只需对模体进行一次或多次扫描,即可得到代入了几何参数的重建公式,然后利用该重建公式对CT扫描得到的待成像物体的扇束投影数据进行运算处理,以得到的CT图像数据为灰度,即可得到没有几何伪影的CT图像,由此可见,本实用新型无需任何手工操作,在CT中的构件具有几何位置偏差的情况下亦可正常使用,得到没有几何伪影的CT图像,相对于现有技术需要手工操作,费时、费力且精确度不高的情况,本实用新型节省时间、操作简单且结果准确。
图2为本实用新型提供的消除CT图像中的几何伪影的系统结构图。如图2所示,CT包括用于发射X射线的X射线源201、带动待成像物体旋转的转台206、探测到达自身的X射线的强度的X射线探测器203、对X射线探测器203所探测到的X射线的强度进行处理从而得到待成像物体205的扇束投影数据的处理装置204,其中,X射线源201具有焦点,其发射的X射线可用于对待成像物体205进行CT扫描;转台206是围绕旋转轴来旋转的,其旋转轴与X射线源201所发射的X射线具有一个交点,该交点为转台206的旋转中心;处理装置204利用扇束投影数据即可进一步得到待成像物体205的CT图像,本实用新型所解决的问题,即为该处处理装置204利用扇束投影数据所得到的CT图像含有几何伪影的问题;该系统进一步包括模体202,这里的模体202具有质心;
X射线源201用于,发射对模体202进行CT扫描的X射线;
转台206用于,带动模体202旋转;
处理装置204用于,对X射线探测器203探测到的X射线的强度进行处理,得到模体202的质心在X射线探测器203上的投影坐标;根据模体202的质心在X射线探测器203上的投影坐标,确定几何参数,并将几何参数代入重建公式,其中,几何参数为用于确定X射线源焦点、转台旋转中心、X射线探测器203之间相对位置的几何参数;利用重建公式处理待成像物体205的扇束投影数据,得到待成像物体205的CT图像数据,则用待成像物体205的CT图像数据所得到的CT图像中没有几何伪影。
由此可见,本实用新型中,由于X射线源201发射的X射线扫描模体202,处理装置204得到模体202的质心在X射线探测器203上的投影坐标,根据该投影坐标即可确定出几何参数,这些几何参数可以用于确定CT中的X射线源201的焦点、转台旋转中心、X射线探测器203之间的相对位置,这样,处理装置204将这些几何参数代入重建公式,然后用重建公式对待成像物体205的扇束投影数据进行处理,得到待成像物体的CT图像数据,然后以CT图像数据未灰度,即可得到精确消除了几何伪影的CT图像。
由于处理装置用CT扫描的方式得到模体的质心在X射线探测器上的投影坐标,处理装置根据该投影坐标即可确定出几何参数,这些几何参数可以用于确定CT中的X射线源焦点、转台旋转中心、X射线探测器之间的相对位置,这样,将这些几何参数代入重建公式,然后用重建公式对待成像物体的扇束投影数据进行处理,用所得到的待成像物体的CT图像数据作为CT图像的灰度,即可精确消除CT图像中的几何伪影,得到没有几何伪影的CT图像;另外,该方法不需人工参与,因而也消除了人工调节误差的存在以及人工误差无法复原的可能性。
该系统中的模体可以为多种形式,通常采用结构简单、易于判断质心的物体,它可以为柱形物体,如相对于X射线束以及探测器而言比较细的细丝等,也可以为球形物体,如相对于X射线束以及探测器而言直径比较小的小球等,当然,模体也可以为专门制作的模体或其他物体。
该系统中,处理装置204用于,对X射线探测器203探测到的X射线的强度进行处理,得到模体202的投影数据;用图像处理方法从模体202的投影数据中得到模体202的质心在X射线探测器203上的投影坐标。
这里,图像处理方法很多,例如二值法、阈值法、质心法等。
模体202的质心在X射线探测器203上的投影坐标为:模体202的质心在X射线探测器203上n个不同角度的投影坐标;这里,n为正整数,模体的质心在X射线探测器上第k个角度处的投影坐标为uk,k为不小于零的整数,第k个角度为γk=β+45*k度,β为0度至45度之间的任一角度。
以n等于8为例,n大于8的情况可以看做是在n等于8的情况的基础上进一步的延伸,n小于8的情况可以类推。
在n等于8的情况下,模体202的质心在X射线探测器203上第k个角度处的投影坐标为第k个角度为γk=β+45*k度;其中,k为0与7之间的所有整数,β为0度至45度之间的任一角度,Maxnum为X射线探测器203的通道总数,pi,k为由X射线探测器203的第i个通道所测得的第k个角度处的投影坐标。
处理装置204所确定的几何参数也有多种形式,只要将其代入相应的重建公式中,对待成像物体205上各点的扇束投影数据进行处理,得到没有几何伪影的CT图像,都在本实用新型的保护范围之内。
例如,几何参数可以包括:转台旋转中心在X射线探测器203上的投影点与X射线探测器203的中心点之间的距离h、中心射线与X射线探测器203所在直线的夹角α、转台旋转中心在X射线探测器203上的投影点与X射线源焦点之间的距离SDD,其中,中心射线为穿过转台旋转中心的X射线。
这样,重建公式为:
其中,f(r0,θ0)为重建图像中极坐标为(r0,θ0)的点的重建图像数据;p(r,θ)为待成像物体205上极坐标为(r,θ)的点的扇束投影数据;
这一组几何参数和重建公式中,h可以为利用以下四个表达式中的任一个计算得到的h,或为利用以下四个表达式中的任意两个以上所计算得到的h的平均值,这四个表达式为:
表达式一:
表达式二:
表达式三:
表达式四:
该组几何参数中, 其中,n1和n2可由以下两组计算公式中的任一组计算得到:
第一组计算公式为:
第二组计算公式为:
几何参数还可以包括:转台旋转中心在X射线探测器203上的投影点与X射线探测器203的中心点之间的距离h以及参数n1、n2,其中,
h为利用以下四个解析式中的任一个计算得到的h,或为利用以下四个解析式中的任意两个以上所计算得到的h的平均值,这四个解析式为:
解析式一:
解析式二:
解析式三:
解析式四:
n1和n2为利用以下两组解析式中的任一组计算得到的n1和n2,两组解析式为:
第一组解析式:
第二组解析式:
利用这组几何参数的重建公式为:
其中,f(x0,y0)为重建图像中坐标为(x0,y0)的点的重建图像数据;p(x,y)为待成像物体205上坐标为(x,y)的点的扇束投影数据;
由上面所述的各式可以看出,几何参数可以为h、α和SDD,也可以为h、n1和n2,这两组几何参数之间是可以相互推导出来的,即这两组几何参数在本质上是同一组几何参数,另外,本实用新型所述的几何参数还可以为其他几何参数,以及这些几何参数能够推导出来的几何参数,只要在将其代入相应的重建公式之后,得到的重建公式中没有几何伪影,都在本实用新型的保护范围之内。
本实用新型所述的系统中,几何参数进一步包括转台旋转中心与X射线源焦点之间的距离SOD。
设第一个模体202和第二个模体202的质心在X射线探测器203上第k个角度处的投影坐标分别为u0k和u1k,u0k和u1k的计算公式如上述uk的计算公式所示,这里,第k个角度为γk=β+45*k度,其中,k为0与n之间的所有整数,n为正整数,β为0度至45度之间的任一角度;则和可利用下列任一组关系式得到:
第一组关系式:
第二组关系式:
第三组关系式:
第四组关系式:
图3为模体为细丝时利用各点投影数据得到的细丝的CT图像,其中,图3a为其整体图像,图3b为图3a中方框部分的图像。如图3a所示,该整体图像是利用CT对细丝进行0-360度的精确扫描得到的1800个投影数据所形成的图像,其表现为一条类似抛物线形状的曲线301。而图3b是对图3a中方框部分的局部放大图,如图3b所示,对图3a中方框部分进行放大后,各点的投影数据所形成的图像是成带状分布的,而不是一条纯粹的曲线,这意味着图3a中的曲线301也是一条带状分布的曲线带,图3b中的曲线302是用质心法从细丝上各点的投影数据中获得的细丝质心在X射线探测器上的投影坐标所组成的曲线,而其他的黑色部分则是质心之外的部分的投影坐标所形成的图像。质心法意味着用于确定几何参数的投影坐标为模体,此处为细丝,的质心在X射线探测器上的投影坐标。
图4为应用本实用新型的CT各构件的位置关系图。如图4所示,S表示X射线源焦点,O表示转台旋转中心;O’表示X射线探测器的中心点;虚线表示X射线探测器所在的直线的理想位置,由于构件的几何位置偏差,该X射线探测器所在直线的实际位置为穿过O’点的u(n1,n2)轴,从而与理想位置有一个角度偏差α,即为本实用新型所述的中心射线与X射线探测器所在直线的夹角,同时,该几何位置偏差也造成转台旋转中心在X射线探测器上的投影点OD与X射线探测器的中心点O’之间具有一定的距离h;横轴,即图4中的ξ轴,为穿过O且平行于虚线的轴;纵轴,即η轴为穿过O和O’的轴;这样,就可以确定出模体上某一点的坐标为P0(ξ0,η0);另外,本实用新型中所述的几何参数还包括转台旋转中心O在X射线探测器上的投影点与X射线源焦点S之间的距离SDD,以及转台旋转中心O与X射线源焦点S之间的距离SOD。
图5为利用现有技术对CT进行校正后扫描两个易拉罐模体得到的CT图像,其中,图5a为整体图像,图5b和图5c分别为图5a中曲线501和502上的方框部分的放大图。图5为利用现有技术,采用手工方法调节CT中的构件,用肉眼来判断CT图像的优劣,从而检验手工调节的质量,这样来消除CT图像中的几何伪影,如图5a所示,曲线501和502分别为靠近X射线探测器中心和远离X射线探测器中心的两个易拉罐的圆形切面的CT图像,对这两条曲线中的方框部分分别进行放大,如图5b和5c所示,两条放大曲线均由两条几何伪影组成,靠近X射线探测器中心的几何伪影之间的距离要比远离X射线探测器中心的几何伪影之间的距离小一些。
图6为假设其他构件的位置正确,手工调节转台旋转中心在X射线探测器上的投影点与X射线探测器的中心点之间的距离而得到的两个易拉罐模体的CT图像,其中,图6a为整体图像,图6b和图6c分别为图6a中曲线601和602方框部分的放大图。图6消除几何伪影的方法为:假设其他构件的位置正确,只有转台旋转中心在X射线探测器上的投影点与X射线探测器的中心点之间的距离不为零,即有距离偏差,这样,采用手工方法调节CT中的构件,用肉眼来判断CT图像的优劣,从而检验手工调节的质量,这样来消除CT图像中的几何伪影。如图6a所示,曲线601和602分别为靠近X射线探测器中心和远离X射线探测器中心的两个易拉罐的圆形切面的CT图像,对这两条曲线中的方框部分分别进行放大,如图6b和6c所示,其中,靠近X射线探测器中心位置的易拉罐的CT图像已消除了几何伪影,而远离X射线探测器中心位置的易拉罐的CT图像则仍由两条几何伪影构成。
综合图5和图6的描述可以看出,利用现有技术的手工调节方式,很难消除CT图像中的几何伪影。
图7为利用本实用新型对CT进行校正后扫描两个易拉罐模体得到的CT图像,其中,图7a为整体图像,图7b和图7c分别为图7a中曲线701和702上的方框部分的放大图。如图7所示,利用本实用新型确定出几何参数,将几何参数代入重建公式,然后利用重建公式对两个易拉罐的扇束投影数据进行处理后,得到图7a的整体图像,靠近X射线探测器中心和远离X射线探测器中心位置的两个易拉罐的重建图像分别为曲线701和702,对这两条曲线上的方框部分进行放大后,得到图7b和图7c,可以看出,图7b和图7c中均没有几何伪影,即利用本实用新型对CT进行校正后得到的CT图像没有几何伪影。另外,通过对图像数据进行分析也可知,图7中各曲线的精确度在一个像素之内,这从根本上保证了曲线上任何一点均不会产生几何伪影。
由此可见,本实用新型具有以下优点:
(1)本实用新型中,由于处理装置用CT扫描的方式得到模体的质心在X射线探测器上的投影坐标,处理装置根据该投影坐标即可确定出几何参数,这些几何参数可以用于确定CT中的X射线源焦点、转台旋转中心、X射线探测器之间的相对位置,这样,将这些几何参数代入重建公式,然后用重建公式对待成像物体的扇束投影数据进行处理,用所得到的待成像物体的CT图像数据作为CT图像的灰度,即可精确消除CT图像中的几何伪影,得到没有几何伪影的CT图像。
(2)本实用新型不需人工参与,因而也消除了人工调节误差的存在以及人工误差无法复原的可能性。
(3)本实用新型采用间接方式来消除CT图像中的几何伪影,只需对模体进行一次或多次扫描,即可得到代入了几何参数的重建公式,然后利用该重建公式对CT扫描得到的待成像物体的扇束投影数据进行运算处理,即可得到没有几何伪影的CT图像,由此可见,本实用新型无需任何手工操作,在CT中的构件具有几何位置偏差的情况下亦可正常使用,得到没有几何伪影的CT图像,相对于现有技术需要手工操作,费时、费力且精确度不高的情况,本实用新型节省时间、操作简单且结果准确。
以上所述仅为本实用新型的较佳实施例,并不用以限制本实用新型,凡在本实用新型的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本实用新型的保护范围之内。
Claims (11)
1.一种消除CT图像中的几何伪影的系统,所述CT包括用于发射X射线的X射线源、带动待成像物体旋转的转台、探测到达自身的X射线强度的X射线探测器、对所述X射线探测器所探测到的X射线的强度进行处理从而得到所述待成像物体的扇束投影数据的处理装置,其特征在于,该系统进一步包括模体;
所述X射线源用于,发射对模体进行CT扫描的X射线;
所述转台用于,带动模体旋转;
所述处理装置用于,对所述X射线探测器探测到的X射线的强度进行处理,得到所述模体的质心在所述X射线探测器上的投影坐标;根据所述模体的质心在X射线探测器上的投影坐标,确定几何参数,并将所述几何参数代入重建公式,其中,所述几何参数为用于确定X射线源焦点、转台旋转中心、X射线探测器之间相对位置的几何参数;利用所述重建公式处理所述待成像物体的扇束投影数据,得到所述待成像物体的CT图像数据,则用所述待成像物体的CT图像数据所得到的CT图像中没有几何伪影。
2.根据权利要求1所述的系统,其特征在于,所述处理装置用于,对所述X射线探测器探测到的X射线的强度进行处理,得到所述模体的投影数据;用图像处理方法从所述模体的投影数据中得到所述模体的质心在X射线探测器上的投影坐标。
3.根据权利要求1或2所述的系统,其特征在于,所述模体的质心在X射线探测器上的投影坐标为:所述模体的质心在X射线探测器上n个不同角度的投影坐标,其中,n为正整数;
所述模体的质心在X射线探测器上第k个角度处的投影坐标为uk,k为不小于零的整数,所述第k个角度为γk=β+45*k度,β为0度至45度之间的任一角度。
4.根据权利要求1或2所述的系统,其特征在于,所述几何参数包括:转台旋转中心在X射线探测器上的投影点与X射线探测器的中心点之间的距离h、中心射线与X射线探测器所在直线的夹角α、转台旋转中心在X射线探测器上的投影点与X射线源焦点之间的距离SDD,其中,所述中心射线为穿过所述转台旋转中心的X射线。
5.根据权利要求4所述的系统,其特征在于,所述重建公式为:
其中,f(r0,θ0)为所述重建图像中极坐标为(r0,θ0)的点的重建图像数据;p(r,θ)为所述待成像物体上极坐标为(r,θ)的点的扇束投影数据;
6.根据权利要求4所述的系统,其特征在于,h为:利用以下四个表达式中的任一个计算得到的h,或,利用以下四个表达式中的任意两个以上所计算得到的h的平均值,所述四个表达式为:
表达式一:
表达式二:
表达式三:
表达式四:
7.根据权利要求6所述的系统,其特征在于,
或,
8.根据权利要求7所述的系统,其特征在于,所述几何参数进一步包括转台旋转中心与X射线源焦点之间的距离SOD。
11.根据权利要求1或2所述的系统,其特征在于,所述模体为柱形物体或球形物体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011200556334U CN202049120U (zh) | 2011-03-04 | 2011-03-04 | 一种消除ct图像中的几何伪影的系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011200556334U CN202049120U (zh) | 2011-03-04 | 2011-03-04 | 一种消除ct图像中的几何伪影的系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN202049120U true CN202049120U (zh) | 2011-11-23 |
Family
ID=44989425
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011200556334U Expired - Fee Related CN202049120U (zh) | 2011-03-04 | 2011-03-04 | 一种消除ct图像中的几何伪影的系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN202049120U (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102652674A (zh) * | 2011-03-04 | 2012-09-05 | 首都师范大学 | 一种消除ct图像中的几何伪影的方法和系统 |
CN106706675A (zh) * | 2015-07-16 | 2017-05-24 | 中国科学院高能物理研究所 | 一种基于计算机分层扫描成像cl系统的校正方法 |
CN107684435A (zh) * | 2017-08-16 | 2018-02-13 | 深圳先进技术研究院 | 锥束ct系统几何校准方法及其校准装置 |
CN108876730A (zh) * | 2018-05-24 | 2018-11-23 | 沈阳东软医疗系统有限公司 | 校正运动伪影的方法、装置及设备和存储介质 |
CN109272472A (zh) * | 2018-10-15 | 2019-01-25 | 天津大学 | 面向医用能谱ct图像的噪声及伪影消除方法 |
CN110461236A (zh) * | 2017-06-27 | 2019-11-15 | 深圳市奥沃医学新技术发展有限公司 | Ct系统参数确定方法及装置 |
CN112748135A (zh) * | 2020-12-28 | 2021-05-04 | 北京航空航天大学 | 一种外部ct扫描模式下转台旋转中心标定方法及装置 |
-
2011
- 2011-03-04 CN CN2011200556334U patent/CN202049120U/zh not_active Expired - Fee Related
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102652674A (zh) * | 2011-03-04 | 2012-09-05 | 首都师范大学 | 一种消除ct图像中的几何伪影的方法和系统 |
CN102652674B (zh) * | 2011-03-04 | 2014-02-19 | 首都师范大学 | 一种消除ct图像中的几何伪影的方法和系统 |
CN106706675A (zh) * | 2015-07-16 | 2017-05-24 | 中国科学院高能物理研究所 | 一种基于计算机分层扫描成像cl系统的校正方法 |
CN110461236A (zh) * | 2017-06-27 | 2019-11-15 | 深圳市奥沃医学新技术发展有限公司 | Ct系统参数确定方法及装置 |
CN110461236B (zh) * | 2017-06-27 | 2022-05-17 | 深圳市奥沃医学新技术发展有限公司 | Ct系统参数确定方法及装置 |
CN107684435A (zh) * | 2017-08-16 | 2018-02-13 | 深圳先进技术研究院 | 锥束ct系统几何校准方法及其校准装置 |
CN108876730A (zh) * | 2018-05-24 | 2018-11-23 | 沈阳东软医疗系统有限公司 | 校正运动伪影的方法、装置及设备和存储介质 |
CN108876730B (zh) * | 2018-05-24 | 2022-03-04 | 东软医疗系统股份有限公司 | 校正运动伪影的方法、装置及设备和存储介质 |
CN109272472A (zh) * | 2018-10-15 | 2019-01-25 | 天津大学 | 面向医用能谱ct图像的噪声及伪影消除方法 |
CN109272472B (zh) * | 2018-10-15 | 2022-07-15 | 天津大学 | 面向医用能谱ct图像的噪声及伪影消除方法 |
CN112748135A (zh) * | 2020-12-28 | 2021-05-04 | 北京航空航天大学 | 一种外部ct扫描模式下转台旋转中心标定方法及装置 |
CN112748135B (zh) * | 2020-12-28 | 2022-03-15 | 北京航空航天大学 | 一种外部ct扫描模式下转台旋转中心标定方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102652674B (zh) | 一种消除ct图像中的几何伪影的方法和系统 | |
CN202049120U (zh) | 一种消除ct图像中的几何伪影的系统 | |
US7978895B2 (en) | X-ray CT system | |
US7672424B2 (en) | Image reconstruction with voxel dependent interpolation | |
CN106667512B (zh) | X射线成像设备的几何校正方法、乳腺断层成像设备 | |
CN103961125B (zh) | 一种用于锥束ct的ct值校正方法 | |
US20090202127A1 (en) | Method And System For Error Compensation | |
WO2018218611A1 (zh) | 锥形束计算机断层成像系统的几何参数确定方法 | |
US20070116183A1 (en) | Tomograph | |
CN103479379B (zh) | 一种倾斜螺旋扫描的图像重建方法及装置 | |
EP2469472A1 (en) | Computed tomography (ct) image reconstruction method and apparatus | |
CN105608720B (zh) | 计算机断层成像伪影校正方法及装置 | |
CN111553849B (zh) | 基于局部特征匹配的锥束ct几何伪影去除方法及装置 | |
CN106651977A (zh) | 基于重建图像梯度的l0范数最小化的锥束ct旋转中心标定方法 | |
CN105844594A (zh) | 一种工业ct扫描系统旋转中心校正方法 | |
CN102488528B (zh) | 一种层析成像几何参数的校准方法 | |
CN105869130B (zh) | Ct系统几何校正效果的检验方法及装置 | |
Sun et al. | Reduction of artifacts induced by misaligned geometry in cone-beam CT | |
JP2004113785A (ja) | コンピュータ断層撮影法における画像形成方法およびこの方法を実施するためのct装置 | |
Duan et al. | Knowledge-based self-calibration method of calibration phantom by and for accurate robot-based CT imaging systems | |
CN107016655A (zh) | 锥束cl几何全参数迭代校正方法 | |
US7394887B2 (en) | Method and apparatus for reconstruction of tilted cone beam data | |
CN109919868A (zh) | 一种锥束ct射束硬化曲线侦测及投影加权校正方法 | |
Aloisi et al. | Investigation on the effects of X-ray CT system geometrical misalignments on dimensional measurement errors | |
CN107730569B (zh) | 一种医学图像伪影校正方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20111123 Termination date: 20150304 |
|
EXPY | Termination of patent right or utility model |