CN103083031A - Ct设备的螺旋扫描图像重建方法、装置和计算机程序产品 - Google Patents

Abstract

本发明公开了CT设备的螺旋扫描图像重建方法、装置和计算机程序产品。在所述方法中，所述CT设备包括多行探测器，所述多行探测器包括具有较小行高的窄探测器和具有较大行高的宽探测器，所述方法包括：对从所述多行探测器得到的投影数据进行预加权处理，以使得从窄探测器得到的投影数据具有与从宽探测器得到的投影数据相同的切片方向灵敏度曲线；对预加权处理后的投影数据直接进行插值处理而忽略窄探测器和宽探测器之间的不同行高；以及对插值处理后的投影数据进行反投影处理以得到重建的CT图像。通过本发明，可以大大减少或消除螺旋扫描CT图像中的伪影。

Description

CT设备的螺旋扫描图像重建方法、装置和计算机程序产品

技术领域

本发明总体上涉及CT(Computed Tomography，计算机X射线断层扫描技术)领域，更特别地涉及一种用于CT设备的螺旋扫描图像的重建方法、装置和计算机程序产品。

背景技术

通常，CT设备用X射线束对人体某部位一定厚度的层面进行扫描，由探测器接收透过该层面的X射线，对所接收的X射线进行光电转换、模数转换，再经过图像重建后得到CT图像。

螺旋扫描是CT设备的常见扫描方式之一，其使得扫描速度大大提高。但是，与螺旋扫描相关联的是CT图像中出现伪影，例如风车状伪影，阶梯状伪影。

伪影的根本原因是z轴方向的插值误差。通常在CT图像重建过程中使用内插技术。在内插时一般选择和使用最接近图像z轴位置的数据来重建图像，而不管这些数据来自于探测器阵列的哪些排(即行)探测器。这样，重建一幅图像所需的数据来自于多排探测器，由于不同排探测器采集数据的角度不同，导致图像在重建时出现周期性误差。重建一幅图像所需要的探测器的排数由螺旋扫描的螺距决定，因此伪影的严重程度与螺距成正比，螺距越大，伪影越严重。

CT图像中各种伪影的存在严重影响了CT图像的质量。

发明内容

因此，需要一种能够减少或消除螺旋扫描CT图像中的伪影的CT图像重建方法。

根据本发明的一方面，本发明提供一种用于CT设备的螺旋扫描图像的重建方法，其中所述CT设备包括多行探测器，所述多行探测器包括具有较小行高的窄探测器和具有较大行高的宽探测器(即行高不是均一的，部分行的行高较宽，部分行的行高较窄)，所述方法包括：

对从所述多行探测器得到的投影数据进行预加权处理，以使得从窄探测器得到的投影数据具有与从宽探测器得到的投影数据相同的切片方向灵敏度曲线；

对预加权处理后的投影数据直接进行插值处理而忽略窄探测器和宽探测器之间的不同行高；以及

对插值处理后的投影数据进行反投影处理以得到重建的CT图像。

根据本发明的一个实施例，其中，所述预加权处理包括对每个投影数据执行以下操作：

如果待处理投影数据是从窄探测器得到的，则对该待处理投影数据以及从与该窄探测器在同一列但在不同行中的至少一个探测器得到的投影数据进行加权平均，并将加权平均值作为相应的经预加权处理后的投影数据；

如果待处理投影数据是从宽探测器得到的，则将该待处理投影数据直接作为相应的经预加权处理后的投影数据。

根据本发明的一个实施例，其中，所述对该待处理投影数据以及从所述至少一个探测器得到的投影数据进行加权平均的步骤进一步包括：

确定分别用于该待处理投影数据以及从所述至少一个探测器得到的投影数据的相应的加权系数；以及

将该待处理投影数据以及从所述至少一个探测器得到的投影数据乘以各自的加权系数，然后相加。

根据本发明的一个实施例，其中，确定加权系数的步骤包括确定所述相应的加权系数以使得从窄探测器得到的投影数据与从宽探测器得到的投影数据具有相同的切片曲线。

根据本发明的一个实施例，其中，确定加权系数的步骤包括迭代地计算所述相应的加权系数以使得窄探测器的经预加权处理后的投影数据的切片方向灵敏度曲线的半高全宽与宽探测器的相等。

根据本发明的一个实施例，其中，从所述至少一个探测器得到的投影数据为从与待处理投影数据所源自的探测器相邻且与其在同一列的探测器得到的投影数据。

根据本发明的一个实施例，其中，在所述迭代计算过程中还使得用于从与待处理投影数据所源自的探测器相邻且与其在同一列的前探测器和后探测器得到的投影数据的相应加权系数的比值等于后探测器的行高的一半加上待处理投影数据所源自的探测器的行高的一半所得到的和与前探测器的行高的一半加上待处理投影数据所源自的探测器的行高的一半所得到的和的比值。

根据本发明的一个实施例，其中，所述宽探测器的行高为所述窄探测器的行高的两倍，并且其中：

(1)如果与待处理投影数据所源自的探测器相邻的前探测器行和后探测器行均为窄探测器行，则

(2)如果与待处理投影数据所源自的探测器相邻的前探测器行为宽探测器行，而后探测器行为窄探测器行，则

以及

(3)如果与待处理投影数据所源自的探测器相邻的前探测器行为窄探测器行，而后探测器行为宽探测器行，则

其中，w_-1为用于从位于所述前探测器行中且与待处理投影数据所源自的探测器在同一列中的探测器得到的投影数据的加权系数，w₀为用于所述待处理投影数据的加权系数，w₁为用于从位于所述后探测器行中且与待处理投影数据所源自的探测器在同一列中的探测器得到的投影数据的加权系数。

根据本发明的一个实施例，所述反投影处理包括：通过矩阵法、迭代法、滤波反投影法或二维傅里叶重建法来对插值处理后的投影数据进行反投影处理以得到重建的CT图像。

根据本发明的一个实施例，其中，所述反投影处理包括：

对插值处理后的投影数据进行卷积滤波；以及

反投影经卷积滤波后的投影数据以得到重建的CT图像。

根据本发明的另一方面，提供另一种用于CT设备的螺旋图像的重建方法，其中所述CT设备包括多行探测器，所述多行探测器包括具有较小行高的窄探测器和具有较大行高的宽探测器，所述方法包括：

对从所述多行探测器得到的投影数据直接进行插值处理而忽略窄探测器和宽探测器之间的不同行高；以及

对插值处理后的投影数据进行反投影处理以得到重建的CT图像。

根据本发明的又一方面，提供一种用于CT设备的螺旋图像的重建装置，其中所述CT设备包括多行探测器，所述多行探测器包括具有较小行高的窄探测器和具有较大行高的宽探测器，所述装置包括：

用于对从所述多行探测器得到的投影数据进行预加权处理、以使得从窄探测器行得到的投影数据具有与从宽探测器行得到的投影数据相同的切片方向灵敏度曲线的部件；

用于对预加权处理后的投影数据直接进行插值处理而忽略窄探测器和宽探测器之间的不同行高的部件；以及

用于对插值处理后的投影数据进行反投影处理以得到重建的CT图像的部件。

根据本发明的一个实施例，其中，所述用于进行预加权处理的部件包括：

用于在待处理投影数据是从窄探测器得到的情况下对该待处理投影数据以及从与该窄探测器在同一列但在不同行中的至少一个探测器得到的投影数据进行加权平均、并将加权平均值作为相应的经预加权处理后的投影数据的部件；

用于在待处理投影数据是从宽探测器得到的情况下将该待处理投影数据直接作为相应的经预加权处理后的投影数据的部件。

根据本发明的一个实施例，其中，所述用于对该待处理投影数据以及从所述至少一个探测器得到的投影数据进行加权平均的部件进一步包括：

用于确定分别用于该待处理投影数据以及从所述至少一个探测器得到的投影数据的相应的加权系数的部件；以及

用于将该待处理投影数据以及从所述至少一个探测器得到的投影数据乘以各自的加权系数、然后相加的部件。

根据本发明的一个实施例，其中，所述用于确定分别用于该待处理投影数据以及从所述至少一个探测器得到的投影数据的相应的加权系数的部件确定所述相应的加权系数以使得从窄探测器得到的投影数据与从宽探测器得到的投影数据具有相同的切片曲线。

根据本发明的一个实施例，其中，所述用于确定分别用于该待处理投影数据以及从所述至少一个探测器得到的投影数据的相应的加权系数的部件包括：用于迭代地计算所述相应的加权系数以使得窄探测器的经预加权处理后的投影数据的切片方向灵敏度曲线的半高全宽与宽探测器的相等的部件。

根据本发明的再一方面，提供一种计算机程序产品，其被配置为存储用于在计算机***上执行的程序指令，所述程序指令使计算机***执行如权利要求1-11中任一项所述的方法。

附图说明

为了对本公开内容更透彻的理解，下面参考结合附图所进行的下列描述，在附图中：

图1是示出用于本发明的实施例中的CT设备的探测器阵列的示意图；

图2是示出现有技术的等高插值方法中的合并处理的示意图；

图3是示出根据本发明的图像重建方法中的预加权处理和插值处理步骤的示意图；

图4A-4B是示出如何确定加权系数的示意图；

图5A-5B分别示出利用现有技术的图像重建方法和本发明的图像重建方法得到的CT图像的效果；

图6A-6C分别示出利用进行等高插值的现有技术图像重建方法、根据本发明的预加权和不等高插值相结合的图像重建方法以及根据本发明的只使用不等高插值的图像重建方法得到的CT图像的效果；以及

图7A-7B分别示出通过现有技术的图像重建方法和本发明的图像重建方法得到的头体CT图像的效果。

具体实施方式

下面将详细描述本发明的具体实施例，但本发明并不限于下述具体实施例。

图1示出了用于本发明的实施例中的CT设备的探测器阵列的示意图。在图1中只是示意性地示出了沿着z轴(即螺旋扫描的旋转轴)方向的多列探测器中的一列探测器。实际上，CT设备的探测器可以是多行×多列探测器。在图1中，该列探测器中的每个探测器代表一个探测器行。如图1所示，探测器阵列在z轴方向上具有不同的行高，其中沿着z轴方向中间的探测器行的行高较窄而两端的探测器行的行高较宽。可以将具有较窄的行高Δh1的探测器称为窄探测器，将具有较宽的行高Δh2的探测器称为宽探测器。但是，图1的探测器阵列只是一种示例，实际上探测器阵列的行高可以是另外的情况。

将从探测器采集到的数据进行光电转换、放大、模数转换等处理可以得到CT扫描的原始数据，也称为扫描数据或投影数据。对于螺旋扫描，在进行反投影之前需要对投影数据进行插值。对于探测器的行高不一致的情况，一般地在现有技术中首先进行合并处理，即将来自若干个相邻窄探测器(在同一列中)的投影数据进行合并，以保持与宽探测器的数据的一致性。如图2所示，宽探测器的行高是窄探测器的行高的2倍，则将两个相邻窄探测器(在同一列中)的投影数据合并在一起。经过这样合并后得到的投影数据集合相当于由等高的宽探测器采集到的投影数据集合。然后，可以对投影数据进行行间插值。这样的插值方式称为等高插值，即先对来自不等高探测器的投影数据进行合并处理以得到等高性质的投影数据，然后再进行插值。然后对经过插值处理的投影数据通过矩阵法、迭代法、滤波反投影法或二维傅里叶重建法来进行反投影以得到重建图像。

本发明采用一种不同的图像重建方法，该方法与上述方法相比可以更好地消除图像中的伪影。下面将详细地描述根据本发明的实施例的图像重建方法。但是，在下面的描述中，重点描述了本发明与现有技术的不同之处，对于本领域技术人员已知或公知的一些处理步骤，没有做详细的描述，以免混淆本发明。

本发明在图像重建过程中引入了预加权的处理步骤，并基于该步骤进行不同于等高插值的插值处理。现在参考图3，其示出了根据本发明的一个实施例的图像重建方法中的预加权处理步骤和插值处理步骤。

对于从不等高的探测器阵列(即探测器具有不同的行高)得到的投影数据，具体的预加权处理步骤如下：如果待处理投影数据是从窄探测器得到的，则对该待处理投影数据以及从与该窄探测器在同一列但在不同行中的至少一个探测器得到的投影数据进行加权平均，并将加权平均值作为相应的经预加权处理后的投影数据；如果待处理投影数据是从宽探测器得到的，则将该待处理投影数据直接作为相应的经预加权处理后的投影数据。这样的预加权处理是为了使从窄探测器得到的投影数据具有与从宽探测器得到的投影数据相同的切片方向灵敏度曲线，从而使得在各个方向采集到的数据具有一致性。

在一个优选实施例中，与该窄探测器在同一列但在不同行中的至少一个探测器是指与该窄探测器相邻且在同一列中的探测器。在图3中，投影数据P_m、P_n、P_o源自窄探测器D_m、D_n、D_o，投影数据P_b来自宽探测器D_b。P_b的经预加权处理后的值P′_b仍然等于P_b。P_m的经预加权处理后的相应值P′_m为P_m、P_m-1(来自与D_m相邻且在D_m之前的行中的探测器D_m-1)、P_m+1(来自与D_m相邻且在D_m之后的行中的探测器D_m+1)这三个投影数据的加权平均。同样，P_n的经预加权处理后的相应值P′_n为P_n、P_n-1(来自与D_n相邻且在D_n之前的行中的探测器D_n-1)、P_n+1(来自与D_n相邻且在D_n之后的行中的探测器D_n+1)这三个投影数据的加权平均，P_o的经预加权处理后的相应值P′_o为P_o、P_o-1(来自与D_o相邻且在D_o之前的行中的探测器D_o-1)、P_o+1(来自与D_o相邻且在D_o之后的行中的探测器D_o+1)这三个投影数据的加权平均。

可以用下面的公式来表示预加权处理：

其中，P′_i为对投影数据P_i经过预加权处理后的得到的投影数据，w_k为加权系数，且

i、k为整数，P_i+k(k＝1或-1)为与从与P_i源自的探测器相邻且与其在同一列中的探测器的投影数据。将w_k的值确定为使得从窄探测器得到的投影数据与从宽探测器得到的投影数据具有相同的切片曲线(slice profile)。

参照图4来详细地说明如何确定加权系数w_k的值。在图4中，以对图3中的投影数据P_n进行预加权处理时所需计算的加权系数为例来进行说明。如上所述，P_n的经预加权处理后的相应值P′_n为P_n、P_n-1(来自与D_n相邻且在D_n之前的行中的探测器D_n-1)、P_n+1(来自与D_n相邻且在D_n之后的行中的探测器D_n+1)这三个投影数据的加权平均。从图3可以看出，D_n、D_n+1、D_n+1均为窄探测器。

图4A为对来自窄探测器D_n的P_n和来自宽探测器D_b的P_b的经预加权处理后的值进行卷积运算的示意图。图4A的上半部分中左边的图表示P_n的预加权处理，即P_n、P_n-1、P_n+1这三个投影数据乘以各自的加权系数然后相加，其中三个矩形分别对应于这三个投影数据，矩形的宽为窄探测器的行高，矩形的高表示投影数据的相应加权系数的大小。图4A的上半部分中右边的图表示用于窄探测器的卷积函数。在本实施例中由于采用线性插值，因此该卷积函数为三角波函数。该函数的峰值为1，宽度为窄探测器的行高的2倍。卷积函数与加权平均值之间的符号“*”表示卷积运算符号。图4A的下半部分中左边的图表示P_b的预加权处理，即P_b乘以加权系数1，其中该图中的矩形的宽为宽探测器的行高，矩形的高度表示P_b的加权系数1。图4A的下半部分中右边的图表示用于宽探测器的卷积函数，该函数的峰值为1，宽度为宽探测器的行高的2倍。同样，卷积函数与加权平均值之间的符号“*”表示卷积运算符号。图4B示意性地示出经过卷积运算所得到的切片方向灵敏度曲线。由于各加权系数是用于加权平均运算，因此各加权系数之和应等于1。通过迭代地调整图4A上半部分中左边图中三个矩形的高度，即各个加权系数的大小，可以使得P′_n的切片方向灵敏度曲线与P′_b的切片方向灵敏度曲线具有相等的半高全宽(FWHM)，这时的各个加权系数值即是合适的值。

在一个优选实施例中，在迭代计算过程中还使得用于与待处理投影数据所源自的探测器相邻且与其在同一列的前探测器和后探测器得到的投影数据的相应加权系数的比值等于后探测器的行高的一半加上待处理投影数据所源自的探测器的行高的一半所得到的和与前探测器的行高的一半加上待处理投影数据所源自的探测器的行高的一半所得到的和的比值。也可以说，在迭代计算过程中还使得用于前探测器的投影数据的加权系数与用于后探测器的投影数据的加权系数的比值是上述图中与前、后探测器对应的矩形的中心到中间矩形的中心的距离比的倒数，我们将这种关系称为加权系数与距离的反比关系。

对于P_n的预加权处理来说，以上述方法进行迭代计算，可以得出相应的各个加权系数的值为：

其中w₀为待处理投影数据P_n的加权系数，w_-1为来自位于与投影数据P_n源自的探测器D_n相邻的前一行探测器中且与D_n在同一列的探测器D_n-1的投影数据的加权系数，w₁为来自位于与D_n相邻的后一行探测器中且与D_n在同一列的探测器D_n-1的投影数据的加权系数。由上面的描述可知，各个加权系数与投影数据的值无关，而是与宽、窄探测器的行高关系、待处理投影数据所源自的探测器以及与该探测器相邻且在同一列中的探测器的类型(是宽探测器还是窄探测器)有关。上面确定出的各个加权系数的值适用于来自于以下这样的窄探测器的待处理投影数据：该窄探测器相邻的前后探测器均是窄探测器，且窄探测器的行高是探测器阵列中宽探测器的行高的二分之一。

对于P_m，可以通过相似的方法确定出所需的各个加权系数。与上面的方法的不同之处在于，图4A上半部分中左边图中的三个矩形中最右边的矩形的宽度应该调整为与宽探测器D_m+1的行高对应；图4A上半部分右边图中的三角波的原点右边部分的宽度应该调整为宽探测器的行高的一半加上窄探测器的行高的一半所得到的和，而三角波原点左边部分的宽度不变。通过上面所述的迭代计算方法，可以确定出各个加权系数为：

同样，这些加权系数适用于来自于以下这样的窄探测器的待处理投影数据：该窄探测器相邻的前探测器是窄探测器而后探测器是宽探测器，且窄探测器的行高是探测器阵列中宽探测器的行高的二分之一。

对于P_o，同样可以通过相似的方法确定出所需的各个加权系数。与上面的方法的不同之处在于，图4A上半部分中左边图中的三个矩形中最左边的矩形的宽度应该调整为与宽探测器D_o-1的行高对应；图4A上半部分右边图中的三角波的原点左边部分的宽度应该调整为宽探测器的行高的一半加上窄探测器的行高的一半所得到的和，而三角波原点右边部分的宽度不变。通过上面所述的迭代计算方法，可以确定出各个加权系数为：

同样，这些加权系数适用于来自于以下这样的窄探测器的待处理投影数据：该窄探测器相邻的前探测器是宽探测器而后探测器是窄探测器，且窄探测器的行高是探测器阵列中宽探测器的行高的二分之一。

需要说明的是，上面描述中提到的使两条切片方向灵敏度曲线的半高全宽相等不一定是指使它们绝对相等，也包括使它们近似相等。

在上面的描述中，k＝-1，0，1为例进行了说明。但是，本领域技术人员应当理解的是，k的取值也可以是其它情况，例如，k可以为从-2至2，这时通过迭代计算仍然可以确定出合适的加权系数值以使得窄探测器的预加权处理后的投影数据的切片方向灵敏度曲线与宽探测器的投影数据的切片方向灵敏度曲线具有相等的半高全宽。另外，宽探测器与窄探测器的行高关系也不必是2倍的关系，可以是3倍或其它比例关系。

经过上面的预加权处理，投影内的数据一致性得到确保，并且也减小了投影之间的差异，从而有助于减少重建图像中的伪影。

返回图3，在对投影数据进行了上述预加权处理之后，开始进行插值。如上所述，在现有技术中，对于来自不等高探测器的投影数据，在插值之前需要先进行合并。但是，在本发明中，由于已经进行了预加权处理，因此可以无需进行合并、忽略窄探测器和宽探测器之间的不同行高而直接进行行间插值。

经过插值处理后，即可以通过诸如矩阵法、迭代法、滤波反投影法或二维傅里叶重建法的图像重建算法来对经预加权和插值处理后的投影数据进行反投影以得到重建图像。

在上面描述的实施例中，通过将预加权处理和不等高插值相结合实现了一种新的图像重建方法。实验结果表明，通过这样的图像重建方法得到的CT图像中的伪影得到有效的减少或消除。

图5A和5B分别示出了利用两种图像重建方法得到的CT图像的效果。图5A对应于先合并、然后进行等高插值的现有技术图像重建方法，图5B对应于根据本发明的预加权处理和不等高插值相结合的图像重建方法。可以看出，使用本发明的图像重建方法得到的CT图像的质量更高，其中风车状伪影大大减少。

另外，作为另一实施例，也可以省略掉预加权处理步骤，而只使用不等高插值处理，即对原始数据直接进行不等高插值处理而忽略窄探测器和宽探测器之间的不同行高，之后对插值后的投影数据进行反投影以重建图像。

图6A-C分别示出了通过三种图像重建方法得到的CT图像的效果。图6A对应于先合并、然后进行等高插值的现有技术图像重建方法，图6B对应于单独进行不等高插值的图像重建方法，图6C对应于预加权处理与不等高插值相结合的图像重建方法。可以看出，图6B中的CT图像中的伪影也大大减少了，但由于没有预加权处理，其数据的不一致性带来了新的伪影-黑影。

图7A和7B分别示出通过现有技术的图像重建方法和本发明的图像重建方法得到的头体CT图像的效果。图7A对应于先合并、然后进行等高插值的现有技术图像重建方法，图7B对应于根据本发明的预加权处理和不等高插值相结合的图像重建方法。可以看出，在图7B的图中，阶梯状伪影大大减少。

本发明的预加权处理和不等高插值相结合的图像重建方法采用了物理上更好的采样，不会影响z轴的空间分辨率。这种图像重建方法不仅减少软组织区域中的伪影，还可以减少骨头区域中的伪影。根据本发明的图像重建方法在减少伪影的同时实现了对具有较窄行高的窄探测器的有效利用，提高了分辨率，并且使得可以利用更大的螺距进行扫描，从而可以进一步提高扫描速度。并且，本发明的这种图像重建方法不会引入其它伪影，也不会导致诸如分辨率和噪声等其它图像质量因素的劣化。

虽然上述已经结合附图描述了本发明的具体实施例，但是本领域技术人员在不脱离本发明的精神和范围的情况下，可以对本发明进行各种改变、修改和等效替代。这些改变、修改和等效替代都意为落入随附的权利要求所限定的精神和范围之内。

Claims (17)

1.一种用于CT设备的螺旋扫描图像的重建方法，其中所述CT设备包括多行探测器，所述多行探测器包括具有较小行高的窄探测器和具有较大行高的宽探测器，所述方法包括：

对从所述多行探测器得到的投影数据进行预加权处理，以使得从窄探测器得到的投影数据具有与从宽探测器得到的投影数据相同的切片方向灵敏度曲线；

对预加权处理后的投影数据直接进行插值处理而忽略窄探测器和宽探测器之间的不同行高；以及

对插值处理后的投影数据进行反投影处理以得到重建的CT图像。

2.如权利要求1所述的重建方法，其中，所述预加权处理包括对每个投影数据执行以下操作：

如果待处理投影数据是从窄探测器得到的，则对该待处理投影数据以及从与该窄探测器在同一列但在不同行中的至少一个探测器得到的投影数据进行加权平均，并将加权平均值作为相应的经预加权处理后的投影数据；

如果待处理投影数据是从宽探测器得到的，则将该待处理投影数据直接作为相应的经预加权处理后的投影数据。

3.如权利要求2所述的重建方法，其中，所述对该待处理投影数据以及从所述至少一个探测器得到的投影数据进行加权平均的步骤进一步包括：

确定分别用于该待处理投影数据以及从所述至少一个探测器得到的投影数据的相应的加权系数；以及

将该待处理投影数据以及从所述至少一个探测器得到的投影数据乘以各自的加权系数，然后相加。

4.如权利要求2或3所述的重建方法，其中，确定加权系数的步骤包括确定所述相应的加权系数以使得从窄探测器得到的投影数据与从宽探测器得到的投影数据具有相同的切片曲线。

5.如权利要求2或3所述的重建方法，其中，确定加权系数的步骤包括迭代地计算所述相应的加权系数以使得窄探测器的经预加权处理后的投影数据的切片方向灵敏度曲线的半高全宽与宽探测器的相等。

6.如权利要求5所述的重建方法，其中，从所述至少一个探测器得到的投影数据为从与待处理投影数据所源自的探测器相邻且与其在同一列的探测器得到的投影数据。

7.如权利要求6所述的重建方法，其中，在所述迭代计算过程中还使得用于从与待处理投影数据所源自的探测器相邻且与其在同一列的前探测器和后探测器得到的投影数据的相应加权系数的比值等于后探测器的行高的一半加上待处理投影数据所源自的探测器的行高的一半所得到的和与前探测器的行高的一半加上待处理投影数据所源自的探测器的行高的一半所得到的和的比值。

8.如权利要求6所述的重建方法，其中，所述宽探测器的行高为所述窄探测器的行高的两倍，并且其中：

(1)如果与待处理投影数据所源自的探测器相邻的前探测器行和后探测器行均为窄探测器行，则

(2)如果与待处理投影数据所源自的探测器相邻的前探测器行为宽探测器行，而后探测器行为窄探测器行，则

以及

(3)如果与待处理投影数据所源自的探测器相邻的前探测器行为窄探测器行，而后探测器行为宽探测器行，则

其中，w_-1为用于从位于所述前探测器行中且与待处理投影数据所源自的探测器在同一列中的探测器得到的投影数据的加权系数，w₀为用于所述待处理投影数据的加权系数，w₁为用于从位于所述后探测器行中且与待处理投影数据所源自的探测器在同一列中的探测器得到的投影数据的加权系数。

9.如权利要求1所述的重建方法，其中，所述反投影处理包括：通过矩阵法、迭代法、滤波反投影法或二维傅里叶重建法来对插值处理后的投影数据进行反投影处理以得到重建的CT图像。

10.如权利要求1所述的重建方法，其中，所述反投影处理包括：

对插值处理后的投影数据进行卷积滤波；以及

反投影经卷积滤波后的投影数据以得到重建的CT图像。

11.一种用于CT设备的螺旋图像的重建方法，其中所述CT设备包括多行探测器，所述多行探测器包括具有较小行高的窄探测器和具有较大行高的宽探测器，所述方法包括：

对从所述多行探测器得到的投影数据直接进行插值处理而忽略窄探测器和宽探测器之间的不同行高；以及

对插值处理后的投影数据进行反投影处理以得到重建的CT图像。

12.一种用于CT设备的螺旋图像的重建装置，其中所述CT设备包括多行探测器，所述多行探测器包括具有较小行高的窄探测器和具有较大行高的宽探测器，所述装置包括：

用于对从所述多行探测器得到的投影数据进行预加权处理、以使得从窄探测器行得到的投影数据具有与从宽探测器行得到的投影数据相同的切片方向灵敏度曲线的部件；

用于对预加权处理后的投影数据直接进行插值处理而忽略窄探测器和宽探测器之间的不同行高的部件；以及

用于对插值处理后的投影数据进行反投影处理以得到重建的CT图像的部件。

13.如权利要求12所述的重建装置，其中，所述用于进行预加权处理的部件包括：

用于在待处理投影数据是从窄探测器得到的情况下对该待处理投影数据以及从与该窄探测器在同一列但在不同行中的至少一个探测器得到的投影数据进行加权平均、并将加权平均值作为相应的经预加权处理后的投影数据的部件；

用于在待处理投影数据是从宽探测器得到的情况下将该待处理投影数据直接作为相应的经预加权处理后的投影数据的部件。

14.如权利要求13所述的重建装置，其中，所述用于对该待处理投影数据以及从所述至少一个探测器得到的投影数据进行加权平均的部件进一步包括：

用于确定分别用于该待处理投影数据以及从所述至少一个探测器得到的投影数据的相应的加权系数的部件；以及

用于将该待处理投影数据以及从所述至少一个探测器得到的投影数据乘以各自的加权系数、然后相加的部件。

15.如权利要求13或14所述的重建装置，其中，所述用于确定分别用于该待处理投影数据以及从所述至少一个探测器得到的投影数据的相应的加权系数的部件确定所述相应的加权系数以使得从窄探测器得到的投影数据与从宽探测器得到的投影数据具有相同的切片曲线。

16.如权利要求13或14所述的重建装置，其中，所述用于确定分别用于该待处理投影数据以及从所述至少一个探测器得到的投影数据的相应的加权系数的部件包括：用于迭代地计算所述相应的加权系数以使得窄探测器的经预加权处理后的投影数据的切片方向灵敏度曲线的半高全宽与宽探测器的相等的部件。

17.一种计算机程序产品，其被配置为存储用于在计算机***上执行的程序指令，所述程序指令使计算机***执行如权利要求1-11中任一项所述的方法。

Images