CN100522063C - 一种双排或多排螺旋ct中的图象重建方法 - Google Patents
一种双排或多排螺旋ct中的图象重建方法 Download PDFInfo
- Publication number
- CN100522063C CN100522063C CNB2004100265969A CN200410026596A CN100522063C CN 100522063 C CN100522063 C CN 100522063C CN B2004100265969 A CNB2004100265969 A CN B2004100265969A CN 200410026596 A CN200410026596 A CN 200410026596A CN 100522063 C CN100522063 C CN 100522063C
- Authority
- CN
- China
- Prior art keywords
- detector
- data
- row
- weights
- projection
- 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 - Lifetime
Links
- 238000000034 method Methods 0.000 title claims abstract description 85
- 238000005070 sampling Methods 0.000 claims abstract description 18
- 230000008569 process Effects 0.000 claims abstract description 11
- 230000008707 rearrangement Effects 0.000 claims abstract description 7
- 238000000205 computational method Methods 0.000 claims description 17
- 238000012545 processing Methods 0.000 claims description 8
- 210000003746 feather Anatomy 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000007704 transition Effects 0.000 claims description 3
- 230000004044 response Effects 0.000 abstract description 7
- 230000009977 dual effect Effects 0.000 abstract 1
- 238000013170 computed tomography imaging Methods 0.000 description 12
- 230000035945 sensitivity Effects 0.000 description 9
- 230000005540 biological transmission Effects 0.000 description 6
- 210000001519 tissue Anatomy 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000010521 absorption reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000007792 addition Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000002184 metal Substances 0.000 description 2
- NJPPVKZQTLUDBO-UHFFFAOYSA-N novaluron Chemical compound C1=C(Cl)C(OC(F)(F)C(OC(F)(F)F)F)=CC=C1NC(=O)NC(=O)C1=C(F)C=CC=C1F NJPPVKZQTLUDBO-UHFFFAOYSA-N 0.000 description 2
- 230000005693 optoelectronics Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000004308 accommodation Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000001915 proofreading effect Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000004846 x-ray emission Methods 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明的一种双排或多排螺旋CT中的图象重建方法，其包括以下步骤：对来自于探测器的原始数据进行校正；进行重建位置上的数据生成，对重建位置前后的采集数据进行纵向插值，近似出重建位置上的数据；进行重排过程，将扇束投影转化为平行投影；对重排后的数据进行滤波；对滤波后的数据进行反投影生成CT图象。本发明由于采用了在应用两点插值方法计算中间采样位置上的数据时，不光考虑了探测器距离采样平面的距离，而且还考虑了探测器在投影中所处的位置因素，这样加重距离中间采样位置近的投影数据的权重，而削减距离中间采样位置远的投影数据的权重，从而改善层灵敏度响应曲线的宽度，提高纵向分辨率。
Description
技术领域
本发明涉及一种双排甚至多排螺旋CT中的图像重建方法，尤其是涉及一种在纵向上应用两点非线性插值方法，生成重建位置上重建数据的方法，这种方法可应用于基于扩展了的高分辨率重建方法的双排或多排螺旋CT中，用于在重建位置上进行重建数据的生成。
背景技术
随着电子技术的发展，医疗器械中出现了无须解剖人体即可成像的CT设备。CT设备的原理是当高速电子流撞击金属时，将产生X-射线，它能穿透人体组织，从而产生透射效果。但X-射线在透射的过程中，不同的组织对射线的吸收程度不尽相同，例如骨骼相比于肌肉，对X-射线的吸收程度要高。为了描述不同的组织对X-射线的吸收程度，每一种组织对应着相应的吸收系数，从而根据位置r上的吸收系数μ(r)，就可知道此位置上的组织情况。为了了解某一层面上μ(r)的分布情况，CT成象系统是根据此层面中的许多方位上，μ(r)沿此方位路径的投影(即透射结果)间接得到的，其中μ(r)沿此方位路径的投影是指细束X-射线沿此路径穿过人体组织后，用所接收到的射线强度去除射线穿入人体组织之前的强度，并将这个比值取自然对数后得到的结果。
为了实现上述过程，在目前的单排CT成象系统上，产生X-射线的球管中的阴极与阳极之间加有很高的电压，阴极端产生的电子束经过高压场之后，以很快的速度打到阳极金属板上，从而产生扇束的X-射线。而与球管相对的位置上安装着一排探测器，用以接收从球管发出的X-射线。为了采样密度的需要，这排探测器的数目通常达到了上千个。这样，当从球管发出的扇束射线到达探测器上时，每个探测器检测出位于此探测器上经过衰减的射线强度，进行光电转换为模拟输出电流。而模数转换电路则将从每个探测器输出的模拟电流信号转换为数字信号，并将此信号送入滤波电路板进行滤波，从而得到当前方位上的采集信号。
不过，仅根据一个方位上的采集信号是不够了解一个层面上μ(r)的分布情况。为此，通过不断旋转装有球管及探测器的DISC装置，从而使得CT成象系统能在足够多的方位上采集到信号。这样，在获取这些方位角上的数据之后，CT成象系统对这些数据进行图象重建，从而得到此层面中μ(r)的分布图。
但根据CT成象系统的扫描方式不同，CT成象系统对采集数据进行图象重建的方法也不一样。在平扫方式下，数据采集过程中，病床是静止不动的，只有当前层面上的数据全部采集完之后，病床才移动到下一个确定的位置。因此，以这种方式采集到的数据是位于同一层面上，CT成象系统能够直接对这些数据进行校正、滤波、反投影等处理过程，从而生成当前层面上的图象，这种图象生成方法称为滤波反投影方法。然而，对螺旋扫描方式来说，在数据采集的同时，病床仍匀速运动。由于这种方式中的数据是连续采集，从而螺旋扫描能在较短的时间内获取全部的数据，大大提高工作效率，这促进了螺旋扫描的广泛应用。但是，以这种方式采集到的不同时刻的数据不位于同一层面上。为了生成一层面上的正确图象，CT成象系统必须根据位于此层面前后的数据，进行纵向插值，从而近似出此层面上的重建数据，之后再对此近似数据利用平扫方式下的图象重建方法生成图象。采用不同的纵向插值方法，所生成的图象效果是不同的，其中重要的一个方面是对纵向分辨率的影响。
在单螺旋扫描中，快速获取体数据与得到较高的纵向分辨率(Longitudinal Resolution)是很难同时兼顾的。为此，人们开发了双排及多排螺旋CT。多排螺旋CT中在病床移动的方向上安装有许多排紧密相连的探测器，每排探测器的宽度，根据厂商的需要，可以设计为相等，也可以依照一定的规律而改变。通过对这些探测器进行排之间的组合，CT成象系统就可以在数据采集中的同一时刻，采集到几个纵向位置上的数据。利用这种方法，加之采用滑环及并行处理技术，在快速获取体数据的同时，纵向分辨率可以得到很大程度的提高。
在多排螺旋CT中，为了得到重建位置(或层面)上的图象，通常的做法是根据多排螺旋CT所采集的数据，插值出此位置上的重建数据，之后，再利用平扫方式下的重建方法生成图象。因此，生成重建位置上的重建数据是非常重要的，它直接影响到重建图象的质量。一般来说，评估重建方法的指标有许多项，如层灵敏度曲线SSP(Slice Sensitivity Profile)、图象噪声、伪影等。其中，根据层灵敏度曲线中的半值全宽FWHM(Full Width atHalf Maximum)，可以反映出图象所利用的实际有效层厚(Effective SliceThickness)，进而了解纵向分辨率的情况。
对于基于扩展了的高分辨率重建方法的双排及多排螺旋CT成象系统来说，为生成目标重建位置上的重建数据，CT成象系统首先利用两点线性插值方法，生成位于目标重建位置所对应的滤波宽度(filter width)范围内一些中间采样位置上的数据，然后再对这些中间采样位置上的数据进行加权，从而得到目标位置上的最终重建数据，请见2001年9月份东芝医学杂志，Y.Ogawa的“纵向分辨率-从单排螺旋CT到多排螺旋CT”(Y.Ogawa，“Helical Reconstruction-From Single-slice CT to Multislice CT，”ToshibaMedical Review，Sept.2001)。在生成某个中间采样位置上的数据时，两点线性插值方法所利用的数据，是位于此中间采样位置两边，距离此采样位置最近的实际接收数据。然而，由于这种多排螺旋CT成象系统，所利用的两点线性插值方法对应的层灵敏度响应曲线较宽，因此降低了纵向分辨率。
发明内容
本发明的目的是提出一种双排或多排螺旋CT中的图象重建方法，基于扩展了的高分辨率重建方法，改进层灵敏度响应曲线，从而提高纵向分辨率。
为实现上述目的，本发明的技术方案如下：
一种双排或多排螺旋CT中的图象重建方法，其包括以下步骤：
a)对来自于探测器的原始数据进行校正；
b)进行重建位置上的数据生成，对重建位置前后的采集数据进行纵向插值，近似出重建位置上的数据，步骤b)中设置有为生成重建位置上的数据而对投影数据计算的第一权值和第二权值，以及综合该第一权值和第二权值的综合系数，处理过程如下：
b1)运用线性插值方法计算第一权值；
b2)根据探测器数据所处的射线源旋转角度β，及射线源到此探测器连线与射线源到扇束旋转中心连线间的夹角α，计算出此探测器数据的第二权值；
b3)根据探测器序号计算对第一及第二权值进行综合的系数；
b4)计算探测器数据的最终插值之权重系数；
c)进行重排过程，将扇束投影转化为平行投影；
d)对重排后的数据进行滤波；
e)对滤波后的数据进行反投影生成CT图象。
所述的图象重建方法，其中，所述步骤b1)中的第一权值计算方法为：
所述探测器为双排探测器；
所述双排探测器中第一及第二排探测器位于重建位置z0时，所对应的旋转角度分别为β1及β2，其中β1<β2，将重建平面上2π的方位角划分成两部分角度范围，设第一排探测器旋转角度为β1时，重建平面上所对应的方位角为0，由此该第二排探测器旋转角度为β2，重建平面上对应的方位角为β2-β1，为生成重建平面上0～β2-β1方位角范围内的重建数据，第一排与第二排探测器都利用β1到β2之间的数据，利用线性插值方法计算的第一排探测器的第一权值为whi2＝1.0-(β-β1)/(β2-β1)，及第二排探测器的第一权值为whi2＝(β-β1)/(β2-β1)；而为生成重建平面上β2-β1到2π方位角范围内的重建数据，第一排探测器利用β2-2π到β1之间的数据，利用线性插值方法计算的第一权值为whi2＝[β-(β2-2π)]/[β1-(β2-2π)]，而第二排探测器利用β2到β1+2π之间的数据，其第一权值为whi2＝1.0-(β-β2)/(β1+2π-β2)。
所述的图象重建方法，其中，所述步骤b2)中的第二权值计算方法为：
对位于采样位置一边的探测器数据，如纵向坐标小于重建位置的数据，在很大的角度变化范围内，其对应的第二个权重系数为零，而对于权重系数为非零时，第二权值随着β及α值的递增而递增，当第一排探测器旋转角度位于β2-2π到β1之间与第二排探测器旋转角度位于β1到β2之间时，其纵向坐标值都小于重建平面的Z值，该第一排探测器的第二权值的一种计算方法为whh＝[(β+2π)-β1-π+2α]/(π+2α)，其中当whh为负时，whh置为0，及第二排探测器的第二权值的计算方法为whh＝(β-β1-π+2α)/(π+2α)，其中当whh为负时，whh置为0；而对于位于采样位置另一边的探测器数据来说，如纵向坐标大于重建位置的数据，当其对应的第二个权重系数为非零时，随着β及α值的递增而递减，当第一排探测器旋转角度位于β1到β2之间与第二排探测器旋转角度位于β2到β1+2π之间时，其纵向坐标值都大于重建平面的Z值，所述第一排与第二排探测器的第二权值的计算方法同为whh＝(β-β1-π+2α)/(-π+2α)，其中当whh为负时，whh置为0。
所述的图象重建方法，其中，所述步骤b3)中的综合系数的计算方法为：
依探测器在投影中所处的位置，一排探测器被划分成四个子区域，如果系数f所对应的探测器位于A1～A2及A3～num区间内，那么系数f的值分别被置0与1.0，而对位于其它两个区域内的情形，f值的计算是根据此探测器距离当前子区域一端的长度与此子区域长度的比值x，依如下三次多项式f＝3*x*x-2*x*x*x计算而成的，对位于0～A1区域内的情况，比值x的计算为x＝(A1-d)/feather，其中feather为当前子区域的长度，d为探测器的序号；并且所述步骤b4)中计算探测器数据的最终插值权重为2.0*[f*whi2+(1.0-f)*whh]，其中，当探测器位于0～A1区域时，whi2重置为0，并且其中上述num值为每排探测器的个数，A3的值由基于扩展了的高分辨率重建方法中探测器偏移度n的值所决定，即等于旋转角度相差π的两组平行投影所重合的区域大小，数据A1的取值可根据实际情况进行调整，而A2的取值依赖于A1的值，两过渡子区域0～A1与A2～A3的大小相等。
所述的图象重建方法，其中，对于多排螺旋CT所述的第一及第二权值的计算方法为：将重建平面上2π的方位角被划分成更多的子方位角范围γ0～γ1，γ1～γ2，γ2～γ3，...，γm～2π，其中γ0＝0，为生成重建平面上γk～γk+1方位角内的重建数据，设此时距离重建平面最近的两排探测器序号分别为r，s，其中第r排探测器对应的投影数据的Z值小于重建平面的Z值，第s排探测器对应的投影数据的Z值大于重建平面的Z值，并且第r排投影数据的旋转角度变化范围为βr1～βr2，第s排投影数据的旋转角度变化范围为βs1～βs2，那么此时第r排投影数据的第一、第二权值分别为whi2＝(β-βr1)/(βr2-βr1)，及whh＝(γ-π+2α)/(π+2α)，其中γ为第r排投影数据的旋转角度为β时，重建平面上所对应的方位角，并且当whh为负数时，置whh为0；所述第s排探测器投影数据的第一、第二权值分别为whi2＝1.0-(β-βs1)/(βs2-βs1)，及whh＝(γ-π+2α)/(-π+2α)，其中γ为第s排投影数据的旋转角度为β时，重建平面上所对应的方位角，并且当whh为负数时，置whh为0。
本发明提供的一种双排或多排螺旋CT中的图象重建方法，与已有的技术相比，本发明由于采用了在应用两点插值方法计算中间采样位置上的数据时，不光考虑了探测器距离采样平面的距离(体现为第一权值)，而且还考虑了探测器在投影中所处的位置因素(体现为第二权值)，这样加重距离中间采样位置近的投影数据的权重，而削减距离中间采样位置远的投影数据的权重，从而改善层灵敏度响应曲线的宽度，提高纵向分辨率。
附图说明
下面结合附图对本发明的实施例作详细描述。
附图中，
图1示出的是本发明的双排或多排螺旋CT中的图象重建方法的双排螺旋CT系统示意图；
图2A和图2B分别示出的是本发明方法中经过重排后平行投影在重建平面上的投影的探测器分布图和旋转角度相差π的两组平行投影在重建平面上的投影相组合的示意图；
图3示出的是本发明方法的图象重建的流程示意图；
图4为本发明方法的生成重建位置上重建数据中插值数据权重的计算方法；
图5为本发明方法的双排螺旋CT的数据采集方式示意图；
图6为本发明方法与现有技术的效果对比图。
具体实施方式
下文，将详细描述本发明的一较佳实施例。
本发明的一种双排或多排螺旋CT中的图象重建方法，设计了一种两点间进行非线性插值的方法，以生成位于目标重建位置所对应的滤波宽度范围内，中间采样位置上的数据。本发明的计算探测器数据插值权重的技术实现方案为：首先利用线性插值方法，计算出探测器数据的第一个权重系数。位于同一组扇束投影内的所有探测器数据，此系数值相同。然后，再根据探测器数据所处的射线源旋转角度β，及射线源到此探测器连线与射线源到扇束旋转中心连线间的夹角α，计算出此探测器数据的第二个权重系数。设β值位于0到2π之间，对位于采样位置一边的探测器数据来说(比如纵向z坐标值小于采样位置值的数据)，只有当β+2α的值大于π时，其对应的第二个权重系数才为非零，并且随着β及α值的递增而递增。而对于位于采样位置另一边的探测器数据来说(比如纵向z坐标值大于采样位置值的数据)，第二个权重系数的变化情况正好相反，即只有当β+2
α的值小于π时，其对应的第二个权重系数才为非零，并且随着β及α值的递增而递减。这样，在获取这两个加权系数之后，将根据探测器序号而计算出的，对这两权值进行综合的系数f及1-f，分别作为第一及第二权值的修正系数，并将这两个修正后的权值相加，作为此探测器的最终权重。由于扩展了的高分辨率重建方法的特点，旋转角度相差π的两组平行投影(扇束投影可以通过重排算法，转化为平行投影的情况)所采集的区域有一部分是不重合的，而对于重合的区域，又可分为两端的过渡区域及中间区域。这样，一个投影区域可分成四个子区域，根据探测器所处的子区域，对第一及第二权值进行综和的系数f的计算方法也不一样。如果探测器位于重合区域的两端，那么将此探测器距离此子区域一端的长度，与此子区域整个长度的比值，作为一个三次多项式的自变量值，进行计算后的多项式结果便为综合系数f的值。如果探测器位于重合区域的中间段，那么f的值置为0。如果探测器位于不重合区域，那么f的值为1。此外，如果探测器位于重合区域内，且此区域位于整个投影区域的一端，那么第一权值重置为0。
本发明的一种双排或多排螺旋CT中的图象重建方法中，在本较佳实施例是基于扩展了的高分辨率重建方法的双排螺旋CT系统，对由数据采集子系统得到的投影数据进行图象重建时，生成重建位置上重建数据的方法。
本发明方法的所述双排螺旋CT系统的框图如图1所示，由一操作台110控制一主控系统120，该主控系统120与一CT主体系统130和一图象重建子系统140通讯连接，所述CT主体系统130包括一基座131以及在该基座131上的框架装置132，在该框架装置132内部为可旋转的DISC旋转装置133，在该DISC旋转装置133的内部容置空间中设置有一X射线球管134以及与该X射线球管134相对设置的数据采集系统135，所述X射线球管134产生X-射线，透过位于所述DISC旋转装置的容置腔内的被扫描者，由与该X射线球管134相对的安装有两排探测器的所述数据采集系统135进行信号接收。
整个CT系统的操作是通过所述操作台110控制所述主控系统120进行的。而所述主控系统120负责整个扫描系统的工作，其中包括启动与关闭所述X射线球管134中的X-射线发射、启动与停止装有所述X射线球管134与所述数据采集系统135的所述DISC旋转装置133的旋转、指挥所述图象重建子系统140接收来自所述数据采集系统133中的原始数据及进行图象重建、图象到所述操作台110的传输等。所述数据采集系统133负责检测经过衰减了的X-射线的强度，并进行光电转换为模拟电流信号，然后经过模数转换为计算机能处理的数字信号，并把该数字信号传输给所述图象重建子系统140。其中所述数据采集系统135的探测器并不是以射线源与旋转中心的连线呈对称分布的，而是偏移n+1/4个探测器夹角增量，其中n为整数，这样使得扇束投影在经过重排过程转换为平行投影后，一组平行投影在重建平面上的投影的探测器分布如图2A所示的，旋转角度相差π的两组投影在重建平面上的投影相互错交，从而在横向上不仅增加了一倍的采样密度，而且采样的范围也扩大了，如图2B所示的。其中num的值为每排探测器的个数，A3的值由上述n的值所决定，数据A1的取值可根据实际情况进行调整，而A2的取值依赖于A1的值，两过渡子区域0～A1与A2～A3的大小是相等的。
本发明方法的所述主控系统120启动所述图象重建子系统140进行图象重建过程如图3所示：
步骤1：对原始数据进行校正；
由于磁泻、空气、非线性、余辉等原因，所述图象重建子系统接收到的原始数据并不准确，必须进行校正。
步骤2：进行重建位置上的数据生成；
如前所述，双排螺旋CT采集的体数据不位于同一层面上，而为了生成重建位置上的图象，就必须通过对重建位置前后的采集数据进行纵向插值，近似出重建平面上的数据，否则图象将有伪影。
步骤3：进行重排过程；
通过重排过程，将扇束投影转化为平行投影，使得之后的处理能按照平行投影的情形进行，简化处理过程。
步骤4：滤波过程；
按照滤波反投影的重建方法，在数据进行反投影之前，必须进行滤波过程，以能得到正确的图象。
步骤5：反投影过程；
经过反投影过程之后，即生成了最终的CT图象。
由于所述步骤2的实现方法对层灵敏度响应曲线的影响非常大，因此如图4所示的给出实现此步骤中数据权重的计算方法：
步骤21：将校正后的体数据运用线性插值方法计算第一个权重；
双排螺旋CT的数据采集方式如图5所示，不失一般性，设第一排与第二排探测器位于重建位置z0时，射线源所对应的旋转角度分别为β1及β2(β1及β2的含义与图1中旋转角度β的相同)，其中图中的螺距p大于1.0及β1<β2。这时，重建平面上2π的方位角被划分成两部分角度范围，如图5所示设第一排探测器旋转角度为β1时，重建平面上所对应的方位角为0，那么当第二排探测器旋转角度为β2时，重建平面上对应的方位角为β2-β1。此外，设每排探测器的层厚为d，螺距比为p，则旋转一周后每排探测器在纵向上的位置将增加或减少p*d。为生成重建平面上0～β2-β1这部分方位角范围内的重建数据，第一排与第二排探测器都利用了β1到β2之间的数据，这时，对第一排探测器利用线性插值方法计算的第一权值为
whi2＝1.0-(β-β1)/(β2-β1)，
及对第二排探测器的第一权值为
whi2＝(β-β1)/(β2-β1).
而为生成重建平面上β2-β1到2π方位角范围内的重建数据，第一排探测器利用了β2-2π到β1之间的数据，利用线性插值方法计算的第一权值为
whi2＝[β-(β2-2π)]/[β1-(β2-2π)]，
第二排探测器利用了β2到β1+2π之间的数据，其第一权值为
whi2＝1.0-(β-β2)/(β1+2π-β2).
步骤22：根据探测器数据所处的射线源旋转角度β，及射线源到此探测器连线与射线源到扇束旋转中心连线间的夹角α，如图1所示，计算出此探测器数据的第二个权重系数。
该第二权值不仅与探测器的旋转角度有关，而且与探测器在扇束中的位置也有关。这个权值加重了距离重建平面近的投影数据的影响，而对距离重建平面远的投影数据，则正好相反。对位于采样位置一边的探测器数据来说，如纵向坐标小于重建位置的数据，在很大的角度变化范围内，其对应的第二个权重系数为零，而对于权重系数为非零时，随着β及α值的递增而递增。如图5所示，当第一排探测器旋转角度位于β2-2π到β1之间与第二排探测器旋转角度位于β1到β2之间时，其纵向坐标值都小于重建平面的Z值。这时，第一排探测器的第二权值的一种计算方法为
whh＝[(β+2π)-β1-π+2α]/(π+2α)，
其中当whh为负时，whh置为0(以下处理相同)，及第二排探测器的第二权值的计算方法为
whh＝(β-β1-π+2α)/(π+2α).
而对于位于采样位置另一边的探测器数据来说，如纵向坐标大于重建位置的数据，当其对应的第二个权重系数为非零时，随着β及α值的递增而递减。根据图5，当第一排探测器旋转角度位于β1到β2之间与第二排探测器旋转角度位于β2到β1+2π之间时，其纵向坐标值都大于重建平面的Z值。这时，第一排与第二排探测器的第二权值的计算方法同为
whh＝(β-β1-π+2α)/(-π+2α).
步骤23：根据探测器序号计算对第一及第二权值进行综合的系数；
将两种权值进行综合的系数f的计算方法依赖于探测器在投影中所处的位置，如附图2B所示，一组探测器被划分成四个子区域，如果系数f所对应的探测器位于A1～A2-1及A3～num区间内，那么系数f的值分别被置0与1.0，而对位于其它两个区域内的情形，f值的计算是根据此探测器距离当前子区域一端的长度与此子区域长度的比值x，依如下三次多项式
f＝3*x*x-2*x*x*x
计算而成的。对位于0～A1-1区域内的情况，比值x的计算为
x＝(A1-d)/feather
其中feather为当前子区域的长度，d为探测器的序号。
步骤24：计算探测器数据的最终插值权重；
根据以上步骤的计算结果，插值数据的最终的权重为
2.0*[f*whi2+(1.0-f)*whh]
其中，当探测器位于0～A1-1区域时，whi2重置为0。
对于多排螺旋CT，重建平面上2π的方位角则被划分成更多的子方位角范围(γ0＝0)～γ1，γ1～γ2，γ2～γ3，...，γm～2π，如同在双排螺旋CT中2π的方位角被划分成0～(γ1＝β2-β1)与(γ1＝β2-β1)～2π两个子方位角范围一样，如步骤21所述。为生成重建平面上γk～γk+1方位角内的重建数据，设此时距离重建平面最近的两排探测器序号分别为r，s，其中第r排探测器对应的投影数据的Z值小于重建平面的Z值，第s排探测器对应的投影数据的Z值大于重建平面的Z值，并且第r排投影数据的旋转角度变化范围为βr1～βr2，第s排投影数据的旋转角度变化范围为βs1～βs2，那么此时第r排投影数据的第一、第二权值分别为
whi2＝(β-βr1)/(βr2-βr1)，
whh＝(γ-π+2α)/(π+2α)，
其中γ为第r排投影数据的旋转角度为β时，重建平面上所对应的方位角，并且当whh为负数时，置whh为0(对第s排探测器中的whh处理相同)。而此时第s排投影数据的第一、第二权值分别为
whi2＝1.0-(β-βs1)/(βs2-βs1)，
whh＝(γ-π+2α)/(-π+2α)，
其中γ为第s排投影数据的旋转角度为β时，重建平面上所对应的方位角。之后，再根据上述步骤23、24计算根据探测器序号对第一及第二权值进行综合的系数及探测器数据的最终插值权重。
如图6所示显示了利用本发明方法两点非线性内插方法生成重建数据的层灵敏度曲线与现有技术相比较的效果示意图，这里的螺距比p为1.2，其中的响应值相对于中心位置上的响应值进行了归一化处理。由该图可见，本发明方法能够改善层灵敏度曲线，并且由图中反映有效层厚的半值全宽，可以看出本发明方法提高了纵向分辨率。根据实验结果，通常的线性内插方法相对于本发明方法来说，其对应的半值全宽，扩大了大约11％。
应当理解的是，对本发明的本领域普通技术人员来说，可以根据本发明的技术构思进行各种可能的改变或替换，而所有这些改变或替换都应属于本发明后附权利要求的保护范围。
Claims (5)
1、一种双排或多排螺旋CT中的图象重建方法，其包括以下步骤：
a)对来自于探测器的原始数据进行校正；
b)进行重建位置上的数据生成，对重建位置前后的采集数据进行纵向插值，近似出重建位置上的数据，步骤b)中设置有为生成重建位置上的数据而对投影数据计算的第一权值和第二权值，以及综合该第一权值和第二权值的综合系数，处理过程如下：
b1)运用线性插值方法计算第一权值；
b2)根据探测器数据所处的射线源旋转角度β，及射线源到此探测器连线与射线源到扇束旋转中心连线间的夹角α，计算出此探测器数据的第二权值；
b3)根据探测器序号计算对第一及第二权值进行综合的系数；
b4)计算探测器数据的最终插值之权重系数；
c)进行重排过程，将扇束投影转化为平行投影；
d)对重排后的数据进行滤波；
e)对滤波后的数据进行反投影生成CT图象。
2、根据权利要求1所述的图象重建方法，其特征在于，所述步骤
b1)中的第一权值计算方法为：所述探测器为双排探测器；
所述双排探测器中第一及第二排探测器位于重建位置z0时，所对应的旋转角度分别为β1及β2，其中β1<β2，将重建平面上2π的方位角划分成两部分角度范围，设第一排探测器旋转角度为β1时，重建平面上所对应的方位角为0，由此该第二排探测器旋转角度为β2，重建平面上对应的方位角为β2-β1，为生成重建平面上0～β2-β1方位角范围内的重建数据，第一排与第二排探测器都利用β1到β2之间的数据，利用线性插值方法计算的第一排探测器的第一权值为whi2＝1.0-(β-β1)/(β2-β1)，及第二排探测器的第一权值为whi2＝(β-β1)/(β2-β1)；而为生成重建平面上β2-β1到2π方位角范围内的重建数据，第一排探测器利用β2-2π到β1之间的数据，利用线性插值方法计算的第一权值为whi2＝[β-(β2-2π)]/[β1-(β2-2π)]，而第二排探测器利用β2到β1+2π之间的数据，其第一权值为whi2＝1.0-(β-β2)/(β1+2π-β2)。
3、根据权利要求2所述的图象重建方法，其特征在于，所述步骤b2)中的第二权值计算方法为：
对位于采样位置一边的探测器数据，如纵向坐标小于重建位置的数据，在很大的角度变化范围内，其对应的第二个权重系数为零，而对于权重系数为非零时，第二权值随着β及α值的递增而递增，当第一排探测器旋转角度位于β2-2π到β1之间与第二排探测器旋转角度位于β1到β2之间时，其纵向坐标值都小于重建平面的Z值，该第一排探测器的第二权值的一种计算方法为whh＝[(β+2π)-β1-π+2α]/(π+2α)，其中当whh为负时，whh置为0，及第二排探测器的第二权值的计算方法为whh＝(β-β1-π+2α)/(π+2α)，其中当whh为负时，whh置为0；而对于位于采样位置另一边的探测器数据来说，如纵向坐标大于重建位置的数据，当其对应的第二个权重系数为非零时，随着β及α值的递增而递减，当第一排探测器旋转角度位于β1到β2之间与第二排探测器旋转角度位于β2到β1+2π之间时，其纵向坐标值都大于重建平面的Z值，所述第一排与第二排探测器的第二权值的计算方法同为whh＝(β-β1-π+2α)/(-π+2α)，其中当whh为负时，whh置为0。
4、根据权利要求3所述的图象重建方法，其特征在于，所述步骤b3)中的综合系数的计算方法为：
依探测器在投影中所处的位置，一排探测器被划分成四个子区域，如果系数f所对应的探测器位于A1～A2及A3～num区间内，那么系数f的值分别被置0与1.0，而对位于其它两个区域内的情形，f值的计算是根据此探测器距离当前子区域一端的长度与此子区域长度的比值x，依如下三次多项式f＝3*x*x-2*x*x*x计算而成的，对位于0～A1区域内的情况，比值x的计算为x＝(A1-d)/feather，其中feather为当前子区域的长度，d为探测器的序号；并且所述步骤b4)中计算探测器数据的最终插值权重为2.0*[f*whi2+(1.0-f)*whh]，其中，当探测器位于0～A1区域时，whi2重置为0，并且其中上述num值为每排探测器的个数，A3的值由基于扩展了的高分辨率重建方法中探测器偏移度n的值所决定，即等于旋转角度相差π的两组平行投影所重合的区域大小，数据A1的取值可根据实际情况进行调整，而A2的取值依赖于A1的值，两过渡子区域0～A1与A2～A3的大小相等。
5、根据权利要求1所述的图象重建方法，其特征在于，对于多排螺旋CT所述的第一及第二权值的计算方法为：将重建平面上2π的方位角划分成更多的子方位角范围γ0～γ1，γ1～γ2，γ2～γ3，...，γm～2π，其中γ0＝0，为生成重建平面上γk～γk+1方位角内的重建数据，设此时距离重建平面最近的两排探测器序号分别为r，s，其中第r排探测器对应的投影数据的Z值小于重建平面的Z值，第s排探测器对应的投影数据的Z值大于重建平面的Z值，并且第r排投影数据的旋转角度变化范围为βr1～βr2，第s排投影数据的旋转角度变化范围为βs1～βs2，那么此时第r排投影数据的第一、第二权值分别为whi2＝(β-βr1)/(βr2-βr1)，及whh＝(γ-π+2α)/(π+2α)，其中γ为第r排投影数据的旋转角度为β时，重建平面上所对应的方位角，并且当whh为负数时，置whh为0；所述第s排探测器投影数据的第一、第二权值分别为whi2＝1.0-(β-βs1)/(βs2-βs1)，及whh＝(γ-π+2α)/(-π+2α)，其中γ为第s排投影数据的旋转角度为β时，重建平面上所对应的方位角，并且当whh为负数时，置whh为0。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2004100265969A CN100522063C (zh) | 2004-03-19 | 2004-03-19 | 一种双排或多排螺旋ct中的图象重建方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2004100265969A CN100522063C (zh) | 2004-03-19 | 2004-03-19 | 一种双排或多排螺旋ct中的图象重建方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1669528A CN1669528A (zh) | 2005-09-21 |
| CN100522063C true CN100522063C (zh) | 2009-08-05 |
Family
ID=35041115
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2004100265969A Expired - Lifetime CN100522063C (zh) | 2004-03-19 | 2004-03-19 | 一种双排或多排螺旋ct中的图象重建方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN100522063C (zh) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100495439C (zh) * | 2005-11-21 | 2009-06-03 | 清华大学 | 采用直线轨迹扫描的图像重建系统和方法 |
| US8781198B2 (en) * | 2008-10-10 | 2014-07-15 | Koninklijke Philips N.V. | High contrast imaging and fast imaging reconstruction |
| CN101694647B (zh) * | 2009-10-10 | 2013-01-23 | 中国科学院化学研究所 | 一种图、谱线实时转换分析方法及系统 |
| CN102727230B (zh) * | 2011-04-02 | 2014-06-04 | 沈阳东软医疗系统有限公司 | Ct扫描图像重建方法及装置 |
| CN103927768B (zh) * | 2013-01-16 | 2015-05-13 | 上海联影医疗科技有限公司 | Ct图像重建方法 |
| CN104167007B (zh) * | 2013-05-17 | 2018-02-27 | 上海联影医疗科技有限公司 | 基于部分扫描的ct图像重建方法、装置及ct设备 |
| WO2015053787A1 (en) * | 2013-10-11 | 2015-04-16 | Analogic Corporation | Tomosynthesis imaging |
| CN108460740B (zh) * | 2018-03-06 | 2021-12-14 | 赛诺威盛科技(北京)股份有限公司 | Ct螺旋重建图像伪影去除方法 |
| CN109685867B (zh) * | 2018-12-29 | 2023-05-30 | 上海联影医疗科技股份有限公司 | 一种ct系统、图像重建方法、装置、设备及介质 |
| CN112884855A (zh) * | 2021-01-13 | 2021-06-01 | 中广核贝谷科技有限公司 | 一种安检ct重建图像的处理方法和装置 |
| CN113298903A (zh) * | 2021-05-25 | 2021-08-24 | 苗闯 | 一种大螺距螺旋ct的重建方法、装置、设备和介质 |
-
2004
- 2004-03-19 CN CNB2004100265969A patent/CN100522063C/zh not_active Expired - Lifetime
Non-Patent Citations (4)
| Title |
|---|
| Multi-slice helical CT:Scan and reconstruction. Hui hu.Medical Physics,Vol.26 No.1. 1999 |
| Multi-slice helical CT:Scan and reconstruction. Hui hu.Medical Physics,Vol.26 No.1. 1999 * |
| 扇束工业CT滤波反投影重构算法的快速实现. 傅健，路宏年.计算机应用研究，第3期. 2003 |
| 扇束工业CT滤波反投影重构算法的快速实现. 傅健，路宏年.计算机应用研究，第3期. 2003 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1669528A (zh) | 2005-09-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6492005B2 (ja) | Ｘ線ｃｔ装置、再構成演算装置、及び再構成演算方法 | |
| US7272429B2 (en) | Methods and apparatus for facilitating a reduction in artifacts | |
| US5663995A (en) | Systems and methods for reconstructing an image in a CT system performing a cone beam helical scan | |
| JP5142664B2 (ja) | Ｘ線コンピュータ断層撮影装置 | |
| CN100522063C (zh) | 一种双排或多排螺旋ct中的图象重建方法 | |
| US8175218B2 (en) | Interpolation interlacing based data upsampling algorithm for cone-beam x-ray CT flying focal spot projection data | |
| JPH0614912A (ja) | コンピュータ断層撮影システムで画像を作成する方法と装置 | |
| JPH1128204A (ja) | Ｘ線ｃｔ装置 | |
| CN104597061B (zh) | 基于虚拟探测器的大视场锥束ct成像方法 | |
| JP2001506143A (ja) | 非対称検出システムを備えたらせんコンピュータ断層撮影法 | |
| WO2012049940A1 (ja) | 医用画像処理装置、ｘ線コンピュータ断層撮影装置および医用画像処理方法 | |
| Hagen et al. | Cycloidal computed tomography | |
| US6332013B1 (en) | Methods and apparatus for tilted helical reconstruction multislice CT | |
| CN110916708A (zh) | 一种ct扫描投影数据伪影校正方法、ct图像重建方法 | |
| US5469487A (en) | CT system with twin fan beam helical scan | |
| WO2000062674A1 (en) | Half field of view reduced-size ct detector | |
| JP2005524455A (ja) | 高分解能ｃｔスキャナ | |
| JP5329204B2 (ja) | Ｘ線ｃｔ装置 | |
| JP2002136510A (ja) | カバー範囲を拡大させたミリメートル未満のｃｔスライスを得るための方法及び装置 | |
| US6535572B2 (en) | Methods and apparatus for compensating computed tomographic channel ganging artifacts | |
| US8705688B2 (en) | Medical image processing apparatus, X-ray computed tomography apparatus, and medical image processing method | |
| JPH0998968A (ja) | 物体の断層写真像を作成する方法及び装置 | |
| CN109448071B (zh) | 一种能谱图像重建方法及系统 | |
| JP2013512725A (ja) | ４分の１検出器オフセットを使用する高解像度ニューテートスライス再構成の方法及びシステム | |
| JP2001128965A (ja) | マルチスライス式ｘ線コンピュータ断層撮影システムにおけるデータ補間のためのシステム及び方法 |
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 | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: Image reconstruction method in double-line or multi-line helical CT Effective date of registration: 20141226 Granted publication date: 20090805 Pledgee: Shenzhen SME credit financing guarantee Group Co.,Ltd. Pledgor: SHENZHEN ANKE HIGH-TECH Co.,Ltd. Registration number: 2014990001139 |
|
| PLDC | Enforcement, change and cancellation of contracts on pledge of patent right or utility model | ||
| PC01 | Cancellation of the registration of the contract for pledge of patent right |
Date of cancellation: 20160119 Granted publication date: 20090805 Pledgee: Shenzhen SME credit financing guarantee Group Co.,Ltd. Pledgor: SHENZHEN ANKE HIGH-TECH Co.,Ltd. Registration number: 2014990001139 |
|
| PLDC | Enforcement, change and cancellation of contracts on pledge of patent right or utility model | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: Image reconstruction method in double-line or multi-line helical CT Effective date of registration: 20160301 Granted publication date: 20090805 Pledgee: Shenzhen SME financing Company limited by guarantee Pledgor: SHENZHEN ANKE HIGH-TECH Co.,Ltd. Registration number: 2016990000154 |
|
| PLDC | Enforcement, change and cancellation of contracts on pledge of patent right or utility model | ||
| PC01 | Cancellation of the registration of the contract for pledge of patent right |
Date of cancellation: 20170330 Granted publication date: 20090805 Pledgee: Shenzhen SME financing Company limited by guarantee Pledgor: SHENZHEN ANKE HIGH-TECH Co.,Ltd. Registration number: 2016990000154 |
|
| PC01 | Cancellation of the registration of the contract for pledge of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: Image reconstruction method in double-line or multi-line helical CT Effective date of registration: 20170419 Granted publication date: 20090805 Pledgee: Shenzhen SME financing Company limited by guarantee Pledgor: SHENZHEN ANKE HIGH-TECH Co.,Ltd. Registration number: 2017990000306 |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| PC01 | Cancellation of the registration of the contract for pledge of patent right |
Date of cancellation: 20180522 Granted publication date: 20090805 Pledgee: Shenzhen SME financing Company limited by guarantee Pledgor: SHENZHEN ANKE HIGH-TECH Co.,Ltd. Registration number: 2017990000306 |
|
| PC01 | Cancellation of the registration of the contract for pledge of patent right | ||
| CP02 | Change in the address of a patent holder |
Address after: No.2 workshop, Lingya Industrial Park, No.1 Tangtou Road, Tangtou community, Shiyan street, Bao'an District, Shenzhen City, Guangdong Province Patentee after: SHENZHEN ANKE HIGH-TECH Co.,Ltd. Address before: 518067, No. 26, Mount Shun Road, Shekou, Guangdong, Shenzhen Patentee before: SHENZHEN ANKE HIGH-TECH Co.,Ltd. |
|
| CP02 | Change in the address of a patent holder | ||
| CX01 | Expiry of patent term |
Granted publication date: 20090805 |
|
| CX01 | Expiry of patent term |