专利名称：用于截断补偿的方法和设备的制作方法
技术领域：
本发明一般涉及计算机断层扫描(CT)图像重建中的方法和设备，尤其涉及用于截断补偿的方法和设备。
背景技术：
在某些扫描条件下，病人部分身体可能会超出检测器所能测量的区域之外，这会导致图像伪影和物体成像不完整。已经公布了一些已知方法，其解决减少伪影，而不是超出视野(FOV)的病人部分身体成像。然而，有必要对超出FOV的病人部分身体成像，例如在2003年6月24日提交的、具有代理人案卷号为129993的、且要求2002年10月4日所提交的、临时申请序号为60/416,072的申请的利益的、题为“用于截断补偿的方法和设备”的共同待审申请序号##/###,###中所描述的。另外，有必要改进超出视野的病人部分身体的图像质量。
发明内容
一方面，提供了一种确定投影是否截断的方法。该方法包括计算扫描物体的每个投影视图中的所有采样的和；确定所计算的和的最大值；当索引k的所有采样的和小于最大值的预先确定的百分比时，求投影视图索引k的多个采样m的均值；将均值和门限t进行比较，当均值大于t时，确定投影为截断，当均值不大于t时，确定投影为未截断。
另一方面，提供了一种成像设备。这种成像设备包括射线源；检测器，其对射线响应，且置于接收该源所放射射线的位置上；以及计算机，其操作地耦合到该射线源和该检测器。配置该计算机以计算扫描物体的每个投影视图中的所有采样的和；确定所计算的和的最大值；当投影视图索引k的所有采样的和小于该最大值的预先确定的百分比时，求索引k中多个采样m的均值；将均值和门限t进行比较，当均值大于t时，确定投影为截断，当均值不大于t时，确定投影为未截断。
再一个方面，提供了一种计算机可读介质，包含配置的编码程序，以指示计算机计算扫描物体的每个投影视图中的所有采样的和；确定所计算的和的最大值；当索引k的所有采样的和小于所述大值的预先确定的百分比时，求投影视图索引k的多个采样m的均值；将均值和门限t进行比较，当均值大于t时，确定投影为截断，当均值不大于t时，确定投影为未截断；并按τ(k)=k2-kk2-k1ξ(k1)+k-k1k2-k1ξ(k2)]]>估计总衰减τ(k)，其中k1和k2是与含至少一个确定为截断投影的截断区域相邻的未截断的视图的多个视图位置的均值；按ξ(k)=Σi=1Np(i,k)]]>计算ξ(k)。进一步配置该程序以指示计算机按λ(k)＝τ(k)-ξ(k)确定衰减差额λ(k)；按η(k)=π2Rl2(k)-xl(k)pl(k)-Rl2(k)arcsin(xl(k)Rl(k))+π2Rr2(k)-xr(k)pr(k)-Rr2(k)arcsin(xr(k)Rr(k))]]>计算要增加的衰减总量η(k)，其中pl(k)和pr(k)分别是左右投影边界上的数个视图的平均采样值幅值，xl(k)、xr(k)、Rl(k)和Rr(k)分别是左右拟合柱面的位置和半径；通过计算比率ϵ(k)=η(k)μwλ(k),]]>将η(k)和λ(k)进行比较，其中μw是水的衰减系数；将ε(k)和门限q进行比较；当ε(k)不大于q时，至少使用η(k)和λ(k)中的一个校正图像。进一步配置该程序以指示计算机当ε(k)大于q时，基于和k1-n视图以及k2+n视图有关的数据，计算ηn(k)，其中n是整数；以及使用ηn(k)校正图像。
图1是CT成像系统实施例图；图2是图1所示系统的方框示意图；图3是仅出现在部分投影角度上的投影截断图；图4是不进行截断校正时的重建图像；图5是根据代理人案卷号为129993描述的方法作截断校正得到的图像；图6是按照投影视图的函数(function of projection view)绘制的所测量的每个投影的衰减总量图；图7是拟合(fitting)过程的图示；图8说明多个图像，且包括图像，其仅用单个视图来估计左截断投影的幅值和斜率pl，sl，且也用单个视图来估计右截断投影的幅值和斜率pr，sr；图9说明如何基于所估计的柱面计算缺少的投影采样；图10是用这里所述方法和设备对图4和图5所示的同样的体模(phantom)扫描进行校正得到的图像。
具体实施例方式
于此提供的截断补偿方法和设备用于扩展CT系统的视野。参考附图，对该设备和方法进行详细的描述，所有附图中的相同的附图标记都指相同的元件。附图的用意在于例证而不用作限制，且附于此的目的是有助于对本发明中示例性方法和设备的实施例进行解释。
在一些已知的CT成像系统配置中，射线源投射的扇形光束被校准在位于笛卡儿坐标系的X-Y平面内，该平面通常称为成像平面。射线光束穿过被成像的物体，如病人。经物体衰减之后，光束照射射线检测器阵列。检测器阵列接收到的衰减后的射线光束的强度取决于物体对射线光束的衰减程度。阵列中的每个检测器元件产生单独的电信号，它是该检测器处的光束衰减测量值。从所有检测器独立地获得衰减测量值以产生透射轮廓(transmissionprofile)。
在第三代CT系统中，射线源和检测器阵列随台架在成像平面内围绕被成像的物体旋转，从而射线光束和物体交叉的角度以常量变化。检测器阵列在一个台架角度上获得的一组射线衰减测量(即，投影数据)被称为一个视图。对物体的扫描包含射线源和检测器旋转一周期间在不同台架角度或角度上获得的一组视图。
在轴向扫描中，处理投影数据以重建相应于穿过物体的二维层(twodimensional slice)图像。从一组投影数据重建图像的方法技术上称之为过滤背景投影技术。该处理将扫描所获得的衰减测量值转换成称为“CT数”或“Hounsfield单位”的整数，用以控制显示设备上相应像素的亮度。
为减少总扫描时间，可作螺旋状扫描。为进行螺旋状扫描，在获取规定量的层数据的时候，要移动病人。这种系统从一次扇形光束螺旋状扫描中产生单个螺旋。由扇形光束映射的螺旋产生投影数据，根据该数据重建每个规定层中的图像。
如这里的用法，在所引述的单数形式元素或步骤之前的“a”或“an”应理解成不排除复数形式的所述元素或步骤，除非明确说明。另外，本发明所称的一个实施例也不要解释成排除包含所述特性的其它实施例的存在。
如这里的用法，短语“重建图像”也不要解释成排除本发明中只产生表示图像的数据而不产生可视图像的实施例。因此，如这里的用法，术语“图像”泛指可视图像和表示可视图像的数据。然而，许多实施例产生(或配置来产生)至少一种可视图像。
图1是CT成像系统(10)的图。图2是图1所示系统(10)的方框示意图。在这个示例性实施例中，CT成像系统10包括台架12，表示第三代CT成像系统。台架12有射线源14，将锥形X射线光束16投射到台架12另一端的检测器阵列18的方向上。
检测器组18由多个检测器行(图中未显示)组成，包括多个检测器元件20，都用于感知穿过物体(如病人22)的投射出的X射线光束。每个检测器元件20产生电信号，以表示照射的射线光束的强度，也即光束穿过物体或病人22时的衰减。成像系统10有多层检测器18，能提供多个表示物体22的图像。这些多个图像中的每一幅图像相应于其中的单独层。层的厚度或缝隙取决于检测器行的厚度。
在扫描获取射线投影数据期间，台架12和其上安装的组件围绕旋转中心24旋转。图2仅显示单行检测器元件20(即检测器行)。然而，多层检测器阵列18包括由检测器元件20构成的多个平行的检测器行，因此相应于多个准平行或平行层的投影数据可以在一次扫描中同时获得。
台架12的旋转和射线源14的操作由CT系统10中的控制机械装置26管理。控制机械装置26包括射线控制器28，其为射线源14提供电源和定时信号；以及台架电动机控制器30，其控制台架12的旋转速度和位置。控制机械装置26中的数据采集系统DAS 32从检测器元件20采样模拟数据，并把这些数据转换成数字信号，用于后续处理。图像重建器34从DAS 32接收采样的数字化的射线数据，并进行高速图像重建。重建图像作为计算机36的输入，将图像存储于海量存储设备38中。
计算机36也通过具有键盘的控制台40接收来自操作员的命令和扫描参数。相关的阴极射线管显示器42使操作员能观察来自计算机36的重建图像和其它数据。计算机36使用操作员提供的命令和参数为DAS 32、射线控制器28和台架电动机控制器30提供控制信号和信息。此外，计算机36操纵工作台电动机控制器44，以控制机动化的工作台46，把病人22放到台架12中。尤其是，工作台46将病人22的部分身体移动通过台架开口48。
在一个实施例中，计算机36包括设备50，如磁盘驱动器或CD-ROM驱动器，以读取计算机可读介质52(如磁盘或CD-ROM)中的指令或数据。在另一个实施例中，计算机36执行存储在固件(未示出)中的指令。通常，对图2的DAS 32、重建器34和计算机36中至少一个的所含的处理器进行编程，以执行以下描述的程序。当然，这种方法不仅限于在CT系统10中使用，并可用于其它类型或种类的成像系统中。在一个实施例中，对计算机36编程，以执行上述功能，相应地，如这里的用法，术语“计算机”不限于技术上所指的集成电路意义上的计算机，而泛指计算机、处理器、微控制器、微机、可编程控制器、ASIC和其它可编程电路。尽管这里描述的方法针对医疗场合，可预期本发明的益处对非医用成像系统，如工业场合或运输场合中常使用的系统，如但不局限于飞机场或其它运输中心的行李扫描CT系统，也是十分关键的。
已经公布了一些已知方法，其解决减少伪影，而不是超出视野的病人部分身体成像。然而，有必要对超出视野的病人部分身体成像。这在许多领域内如肿瘤学、旋转血管造影术、融合成像系统以及商用CT扫描器十分有用。至少有一种已知的多层CT扫描器，其现有硬件将重建视野限制在大约50cm范围内。尽管这对大多数医疗应用是足够的，也有必要对视野进行扩展，以对该视野之外的物体成像。这对诸如肿瘤学或CT/PET等应用非常有益。对肿瘤学应用来说，需要一个较大的视野。这主要是因为对射线诊疗方案而言，病人的四肢常常放在扫描视野之外，以更好的定位瘤的位置。已知的CT重建算法不考虑截断投影，且产生的图像有严重的伪影。从治疗方案讲，这些伪影可能影响精确估计衰减路径。这里所描述的是一种算法途径，以将重建视野(FOV)扩大到由检测器硬件所限制的FOV之外。校正算法可应用于不同的重建算法，包括但不受限于全扫描、半扫描/分割、螺旋状扫描、以及基于心脏部分(cardiac sector)的算法。另外，系统10的配置采用了这里描述的算法。
图3是出现在部分投影角度上的投影截断图。在某些扫描条件下，病人22的部分身体超出检测器18所能测量的区域，这会导致图像伪影和成像物体不完整。X射线管和检测器18被牢牢固定在绕病人22旋转的框架上。在旋转过程中，位于完全采样视野60内，会获得不间断的测量。穿过物体22的完全采样视野60之外的任何区域的射线衰减只能在有限的旋转角范围内测量到，且该区域被称为部分采样视野区域。换句话说，位于完全采样视野60内的部分被定位在扇形16之内，因此从所有台架角度都能测量到，且所测得的数据被定义成完全采样视野数据。然而，一些部分在某些角度上位于扇形16之内，而在另一些角度上却位于扇形16之外，与这些部分有关的采集数据被定义成部分采样视野数据。在这个例子中，时针3所指位置上的投影不存在截断，但时针12所指位置上的投影严重截断。因此，可依赖未截断的投影(如图3中大约时针3所指位置上)估计截断视图的截断总量(如图3中大约时针12所指位置上)。校正处理初期步骤是在处理投影之前用软件完成扇形光束到平行光束的转换。在一个实施例中，该初期步骤是第一步。该处理技术众所周知，且不需要专门的数据收集。当转换完成之后，将所有检测器通道上的投影集中起来，以获得下面所详细描述的总衰减曲线。
图4是不进行截断校正时的重建图像80。容易观察到，图像80包含截断伪影82。
代理人案卷号为129993描述了一种针对大多数情况的截断补偿算法。然而，在某些情况下，代理人案卷号为129993中描述的补偿算法不能获得最优结果。图5是根据代理人案卷号为129993中的方法作截断校正后的图像。尽管图像84比图像80有很大改进，但校正区域86(恢复手臂)内仍存在失真。在该例子中，胸模置于最外端区域距离相同中心29cm的位置上(在58cm FOV边缘)。尽管上述校正消除了50cm视野以内的伪影，但所恢复手臂86的形状和强度从某种程度上而言存在失真。
相信失真主要由所估计的缺少衰减不精确引起。为了说明起见，图6是根据投影视图的函数绘制的所测得的每个投影的总衰减90的图。为清楚起见，对该图进行了放大。在衰减测量值90中，下降区92相应于截断视图。从该图中可清楚看到，对多个未截断投影(衰减测量值90中的近似平坦区域)而言，为每个投影所测得的投影测量值的和在所有投影角度上不是常数。这是由非理想化校准和扇形到平行转换处理(parallel rebinning process)引起的。尽管变化较小(小于2％)，它也能表示整个截断投影采样中的主要部分。换句话说，和所估计的缺少投影相比，该变化百分比要大得多。因此，把其他的信息用于截断估计，而不是用于代理人案卷号为129993中所描述的补偿，将会带来更少的、甚至没有失真。
在多个体模试验期间，发现仅投影边界上的幅值和斜率有时不足以估计拟合投影的柱面物体的大小和位置，因为有时边界采样的斜率偏离所扫描物体的内部结构。例如，如果一个气穴碰巧位于截断投影边界采样的切线上，在外插投影采样中会出现明显的不连续(在投影角度上)。这会引起阴影和拖尾伪影。
当投影采样有噪声时(由于有限的x射线光子统计)，发生另一种复杂情况。在这种情况下，投影边界的幅值和斜率都变得不可靠。高频变化(从一个视图到另一个视图)可导致严重的拖尾以及物体形状高度失真。临床应用上，这常发生在扫描体大的病人时，或当使用低流量(low-flux)获取协议时。这里描述的方法和设备至少部分解决上述所描述的难题。
更明确地讲，投影被声明为截断或未截断。要声明为截断投影，必须满足一个条件，即投影边界采样的任一端(在几个通道上的均值)必须小于预先确定的门限，除先前代理人案卷号为129993描述的条件之外。多个采样的均值是指多个边缘采样的均值。即如果投影包含N个采样，索引从0到N-1，使用从0到m-1左面m个投影采样的均值以及从N-m到N-1右面m个采样的均值。除上述测试之外，仍使用代理人案卷号为129993中的测试。即计算所有视图的所有N个采样的和。计算所有投影视图中的最大值，且当一个特定视图的和小于预先确定的最大值的百分比时，该视图被声明为截断的。门限值非零以考虑非理想校准。从算式上讲，如果下列条件存在，投影被声明为截断pl=1mΣi=1mp(i,k)>t]]>或pr=1mΣi=1mp(N-i,k)>t---(1)]]>其中N是检测器通道的数量，k是投影视图索引，m是采样数，t是门限值。在示例性的实施例中，t约等于0.42时，从经验上讲，性能最好。在其它实施例中，t约在0.25和0.58之间，t约在0.33和0.5之间，t约在0.375和0.46之间。
用上面描述的确定方法，识别的未截断视图在图6中以双箭头显示。为了估计缺少投影的量，相邻未截断投影的总衰减用于估计。尤其，截断视图的总衰减通过对位于截断视图区域两端的两个未截断投影线性内插获得。对图6中所示的例子来说，第一组截断视图位于投影视图号290(投射角106°)和445(投射角163°)之间。从视图号289和446获得的总衰减用于估计截断视图的总衰减。通过对两个未截断视图抽样线性内插进行估计，如图6中的点线所示。从算式上而言，所估计的总衰减τ(k)是τ(k)=k2-kk2-k1ξ(k1)+k-k1k2-k1ξ(k2)---(2)]]>其中k1和k2是和截断区域相邻的未截断视图的视图位置，ξ(k)是通过将视图k的投影测量值求总和而获得的总衰减。从算式上而言，它由下面的等式描述ξ(k)=Σi=1Np(i,k)---(3)]]>为改进算法的稳健性，可以使用几个未截断投影的投影视图的均值。例如，可以使用视图285到289的平均总衰减来代替只根据视图289获得的总衰减，且k1为视图285到289的平均值。在示例性的实施例中，取五个视图的平均值，但在其它实施例中，取非五个视图的平均值(如3-7个视图或2-8个视图)。另外，k2可指单个视图或与k1无关的几个视图的平均值。例如，参照图6中的例子，若k1指视图285到289的平均值，则在一个实施例中，k2指视图446到450的平均值，在另一个实施例中，k2指446。另外，k1和k2无关，且k1可指单个视图或与k2无关的几个视图。
用于截断补偿的投影增加量λ(k)是所估计的总衰减和所测量的总衰减之间的差值λ(k)＝τ(k)-ξ(k)(4)下一步，估计可拟合每个截断投影边界的水柱面物体的大小和位置以推导待估计的投影。基于m个投影视图的截断投影边界采样的平均值的幅值和斜率计算拟合柱面的位置和半径xl(k)=-sl(k)pl(k)4μw2,]]>和Rl(k)=pl2(k)4μw2+xl2(k)---(5)]]>xr(k)=-sr(k)pr(k)4μw2,]]>和Rr(k)=pr2(k)4μw2+xr2(k)---(6)]]>其中pl，sl和pr，sr分别是左右数个视图的投影边界的采样平均值幅值和斜率，xl、xr、Rl和Rr分别是左右拟合柱面的位置和半径，μw是水的衰减系数。拟合过程如图7所示。
此前，表明在缺少x射线光子的条件下，参数估计过程会变得不可靠，因为所估计的参数从一个视图到另一个视图会有极大波动。图8说明包含图像102的多个图像100，其只使用单个视图估计左截断投影的幅值和斜率pl，sl，也只使用单个视图估计右截断投影的幅值和斜率pr，sr。图像102中清楚可见由参数波动引起的拖尾伪影。为克服这种可能的不足，在一个实施例中，使用几个视图的投影采样的均值，以估计截断投影的幅值和斜率pl，slpr，sr。在示例性的实施例中，使用7个视图的投影采样估计幅值和斜率。多个图像100中的平均图像104使用7个视图估计幅值和斜率。在另一个实施例中，5到9个视图用于估计幅值和斜率。之后，用平均值确定拟合柱面的位置和半径。基于拟合柱面，对于投影k增加的柱面的总衰减η(k)(在最初投影之外的柱面部分)按下式计算η(k)=π2Rl2(k)-xl(k)pl(k)-Rl2(k)arcsin(xl(k)Rl(k))+π2Rr2(k)-xr(k)pr(k)-Rr2(k)arcsin(xr(k)Rr(k))---(7)]]>在一个实施例中，将η(k)与所估计的增量衰减测量值λ(k)(如等式4描述)进行核对。在理想情况下，这两个量应该相等。但在实际情况下，因为物体的缺少部分极少为一个理想柱面，所以这两个量之间总是有一些差别。因此，计算这两个量的比值ε(k)ϵ(k)=η(k)μwλ(k)---(8)]]>当ε(k)完全不等于1时，说明用拟合的柱面参数有问题。这可能是由被扫描物体的不规则形状引起的，正如上面所描述的。因此，当ε(k)＞q时，不能使用所估计的参数作投影补偿。q是一个预先确定的门限，且在示例性实施例中设为2。在另一个示例性实施例中，q在1.5和2.5之间。在另一个实施例中，q在1.75和2.25之间。在另一个示例性实施例中，q在1.9和2.1之间。
对ε(k)＞q的投影来说，可使用参数拟合较好的两个投影(即ε(k)≤q)的柱面参数的内插。即如果视图k1到k2的ε(k)值超过门限，可使用视图k1-n和k2+n的拟合柱面参数来估计这些视图的柱面参数，其中n是参数，n＝1表示第一个较好拟合的视图。在示例性的实施例中，选择n＝5，采用线性内插足够好。即xl(k)、xr(k)、Rl(k)和Rr(k)是基于xl(k1-n)、xl(k2+n)、xr(k1-n)、xr(k2+n)、Rl(k1-n)、Rl(k2+n)、Rr(k1-n)和Rr(k2+n)的线性内插，且使用内插按等式7计算ηn(k)。比值ε(k)也从同样的视图组作线性内插。当然，也可使用其它内插方法，如高阶插值。另外，沿视图k的额外低通过滤可被应用于所估计的参数以改进参数的连续性。在其它实施例中，n在2和8之间，或n在3和7之间。
最后，基于所估计的柱面计算缺少的投影采样。图9图示这个过程。因为这些柱面的总衰减不能匹配按一致性条件(等式4)估计的整个缺少衰减，要作进一步调整。例如，x轴(检测器通道刻度i)按比率缩放，使初始估计投影(图9中点线所示)在x上进一步扩展(如图9中实粗线所示)或压缩。这可通过改变检测器通道i的刻度距离，或对初始计算出的投影进行内插来实现。
用这里描述的方法和设备，代理人案卷号为129993所描述的方法和设备的性能和稳健性大大提高。图10是用这里所述方法和设备对图4和图5所示同样的体模扫描进行校正得到的图像。位于扫描视野之外的物体106形状和密度的精确性明显好转，这是本发明的技术效果和贡献。当然，其它技术效果也存在。
尽管上面所描述的系统和方法只使用了总衰减、边界采样的幅值和斜率来估计缺少投影的分布，其他信息也能用于估计。例如，对X线断层摄影术而言，可使用Helgason-Ludwig条件(HL条件)进一步提炼上述技术。另外，设置不同的门限以确保算法在错误测量条件下正常工作。例如，可设置图9中拉伸率的上下限，以防止由于不可靠的测量造成误差增加。另外，为所计算的斜率sl，sr设置一个合理的范围。如果被扫描物体的材料特征和水完全不同，也可使用已知材料的衰减系数来计算等式3和4中的大小和位置。
因为内插数据不可能有和完全采样视野内的数据一样的图像质量，所以对FOV外推的图像加以标记就十分有用。在一个实施例中，提供了在重建图像中描绘代表完全采样视野数据和部分采样视野数据区域的方法。也可以在CT号中用色彩代码和变换的CT号来标记。因为标记可能影响图像数据的可视性，提供了打开和关闭标记的简易方式。因此，允许系统10的用户打开和关闭该标记。
尽管本发明参照不同的特定实施例进行描述，本领域的技术人员理解本发明能够在权利要求的实质和范围内加以修改。
权利要求
1.一种成像设备(10)，包括射线源(14)；检测器(18)，其对射线响应，并置于接收从所述源放射的射线的位置；以及计算机(36)，其操作地耦合到所述射线源和所述检测器，且所述计算机配置为计算扫描物体(22)的每个投影视图的所有采样的和；确定所计算的和的最大值；当索引k的所有采样的和小于所述最大值的预先确定的百分比时，求投影视图索引k的多个采样m的均值；将均值和门限t进行比较；当均值大于t时，确定投影为截断；以及当均值不大于t时，确定投影为未截断。
2.根据权利要求1所述的设备(10)，其中进一步配置所述计算机(36)以将均值和门限t进行比较，其中t大约在0.25和0.58之间。
3.根据权利要求1所述的设备(10)，其中进一步配置所述计算机(36)按下式估计总衰减τ(k)=k2-kk2-k1ξ(k1)+k-k1k2-k1ξ(k2)]]>其中k1和k2是与含至少一个确定为截断的投影的截断区域相邻的未截断视图的视图位置，且按ξ(k)=Σi=1Np(i,k)]]>计算ξ(k)。
4.根据权利要求1所述的设备(10)，其中进一步配置所述计算机(36)以按下式估计总衰减τ(k)=k2-kk2-k1ξ(k1)+k-k1k2-k1ξ(k2)]]>其中k1和k2是与含至少一个确定为截断的投影的截断区域相邻的未截断的视图的多个视图位置的均值，且按ξ(k)=Σi=1Np(i,k)]]>计算ξ(k)。
5.根据权利要求4所述的设备(10)，其中进一步配置所述计算机(36)为按λ(k)＝τ(k)-ξ(k)确定衰减差额λ(k)；按η(k)=π2Rl2(k)-xl(k)pl(k)-Rl2(k)arcsin(xl(k)Rl(k))+π2Rr2(k)-xr(k)pr(k)-Rr2(k)arcsin(xr(k)Rr(k))]]>计算要增加的衰减总量η(k)，其中pl(k)和pr(k)分别是左右投影边界上的数个视图的平均采样值幅值，xl(k)、xr(k)、Rl(k)和Rr(k)分别是左右拟合柱面的位置和半径；通过计算比率ϵ(k)=η(k)μwλ(k),]]>将η(k)和λ(k)进行比较，其中μw是水的衰减系数；将ε(k)和门限q进行比较；当ε(k)不大于q时，至少使用η(k)和λ(k)中的一个校正图像；以及当ε(k)大于q时基于和k1-n视图以及k2+n视图有关的数据，计算ηn(k)，其中n是整数；以及使用ηn(k)校正图像。
全文摘要
一种方法包括计算扫描物体的每个投影视图中的所有采样的和；确定所计算的和的最大值；当索引k的所有采样的和小于最大值的预先确定的百分比时，求投影视图索引k的多个采样m的均值；将均值和门限t进行比较；当均值大于t时，确定投影为截断，当均值不大于t时，确定投影为未截断。
文档编号A61B6/03GK1575760SQ200410061719
公开日2005年2月9日 申请日期2004年6月30日 优先权日2003年7月3日
发明者谢强, 布赖恩·格雷科维茨, 爱德华·H·赵, 斯科特·M·麦科拉什, 埃米·L·霍斯特 申请人:Ge医药系统环球科技公司