用于误差补偿的方法和系统的制作方法

专利名称：用于误差补偿的方法和系统的制作方法
用于误差补偿的方法和系统
本发明涉及用于生成核函数集(a set of kernel)的方法，用于误差补偿 的方法和系统，计算机可读介质和程序单元，尤其涉及X射线散射的基于 巻积的误差补偿方法。
计算机断层摄影术(CT)是使用数字处理从围绕单一旋转轴获得的一 系列二维x射线图像生成研究对象(感兴趣对象，检査对象)的内部的三 维图像的过程。CT图像的重建可以通过应用适当的算法进行。
CT成像的基本原理是通过CT系统的探测器获得检查对象的投影数 据。投影数据代表辐射束所经过的对象的信息。为了从投影数据当中生成 图像，这些投影数据(线积分)可以被反投影，从而产生二维图像，艮P， 代表盘。从多个这样的二维图像当中可以重建所谓的体素表示，即，三维 像素的表示。在探测器已经以平面的形式布置的情况下，获得二维投影数 据并且反投影的结果是三维体素。也就是说，在现代的、更复杂的所谓的 "锥束"CT和重建方法中，二维探测器即具有以矩阵的形式布置的多个探 测元件的探测器的投影数据在单一重建步骤中被直接反投影到体素的三维 分布中。
散射辐射是锥束X射线计算机断层摄影术中的伪影的主要来源。通过 在重建图像中导致伪影例如噪声、条纹和低频不均匀性，所谓的杯状伪影， 散射阻碍了柔和对比度即具有低对比度的部分的可见性。尤其在使用防散 射滤线栅效率低的介入性X射线系统的体积成像中，需要用于散射补偿的 可靠和准确的回顾性方法。 一种用于校正的方法是所谓的基于巻积的方法， 该方法经常用于估计射线摄影图像的散射背景。例如，这样的基于巻积的 方 去在"Computerized scatter correction in diagnostic radiology", K.RMaher 和J.RMalone， Contemporary Physics 38(2)， 131-148， 1997中得以描述。
尽管这些基于巻积的校正方法的确增加了重建图像的质量，但是重建 图像仍然呈现出伪影，特别是在体积图像中。希望提供用于生成核函数集的备选方法，用于误差补偿的方法和系统， 计算机可读介质和程序单元，其可以在误差补偿中表现出更大的准确度或 者可以在重建图像中较少地倾向于产生伪影。
根据独立权利要求的用于生成核函数集的备选方法，用于误差补偿的 方法和系统，计算机可读介质和程序单元可以满足该需要。
根据用于针对由成像系统记录的物理对象的投影图像的巻积误差补偿 生成核函数集的方法的一个典型实施例，所述方法包括以这样的方式计算 核函数集使得对于投影图像的每个像素计算用于误差补偿的通常不对称的 散射分布，从而表示在X射线源和像素之间由射束定义的体积中发生的X 射线散射。
根据一个典型实施例， 一种用于对物理对象的图像进行误差补偿的方 法包括接收被成像物理对象的原始投影图像，将原始投影图像转换成水当
量图像(water equivalent image)(水当量厚度的图像)，从所述水当量图像 并且可能从它的梯度提取许多标量参数，根据用于生成核函数集的方法的 一个典型实施例通过使所提取参数与预先计算的核函数的参数相关来确定 至少一个预先计算的核函数，和通过使用所确定的所述至少一个预先计算 的核函数来补偿原始投影图像的误差。
根据一个典型实施例， 一种用于对物理对象的图像进行误差补偿的系
统包括接收单元，其适于接收被成像物理对象的原始投影图像；计算单 元，其适于将原始投影图像转换成水当量图像，特别是计算相应的梯度图 像，并且从水当量厚度的图像并且特别从所述梯度提取许多参数；确定单 元，其适于根据用于生成核函数集的方法的一个典型实施例通过使所提取 参数与预先计算的核函数的参数相关来确定至少一个预先计算的核函数； 和补偿单元，其适于通过使用所确定的所述至少一个预先计算的核函数补 偿原始投影图像的误差。
根据一个典型实施例， 一种断层摄影装置包括辐射源，辐射探测器，
和用于根据一个典型实施例进行误差补偿的系统，其中，所述辐射探测器 适于记录代表被成像物理对象的原始投影图像的信息的数据。
根据一个典型实施例，提供了一种计算机可读介质，用于对物理对象的图像进行误差补偿的程序存储在所述计算机可读介质中，当由处理器执 行时，所述程序适于控制一种方法，所述方法包括接收被成像物理对象 的原始投影图像，将原始投影图像转换成水当量图像，根据用于生成核函 数集的方法的一个典型实施例通过使所提取参数与预先计算的核函数的参 数相关来确定至少一个预先计算的核函数，和通过使用所确定的所述至少 一个预先计算的核函数补偿原始投影图像的误差。
根据一个典型实施例，提供了一种用于对物理对象的图像进行误差补
偿的程序单元，当由处理器执行时，所述程序适于控制一种方法，所述方 法包括接收被成像物理对象的原始投影图像，将原始投影图像转换成水 当量图像，特别是计算相应的梯度图像，从所述水当量图像并且特别从所 述梯度图像提取许多标量参数，根据用于生成核函数集的方法的一个典型 实施例通过使所提取参数与预先计算的核函数的参数相关来确定至少一个 预先计算的核函数，和通过使用根据用于生成核函数集的方法的一个典型 实施例确定的所述至少一个预先计算的核函数补偿原始投影图像的误差。
提供一种用于预先计算核函数的方法可以被看作是本发明的一个典型 实施例的要旨，取决于射线穿透被成像对象的位置，所述核函数足够解决 沿着射线生成的散射分布的不对称。可以看到本发明的一个典型方面在于 本发明精确地解决了如下事实，即，被探测散射X射线量子的大部分可以 源自靠近被成像物理对象的边界的区域并且沿着到这样的位置的路径生成 的散射分布可能是高度不对称的。典型实施例可以提供校正方案，所述校 正方案可以提供可能性以更加定量地估计和校正射线摄影图像中的散射和 锥束计算机断层摄影(CT)采集的投影。由此，与基于巻积的方法相比， 可能减小图像伪影并且因此可能增强低对比度可见性，所述基于巻积的方 法不会考虑对射线穿透被成像对象的位置的依赖性，例如，被考虑像素与 物理对象的中心相关还是与边界区域相关。优选地，核函数集的计算以这 样的方式进行，即使得对于投影图像的每个像素计算用于误差补偿的不对 称的散射分布，从而表示沿着从X射线源到像素的射线发生的X射线散射， 其中不对称可以表示不存在对称轴。特别地，该不对称甚至可以在不使用 防散射滤线栅的情况下存在。
提供的巻积相关散射估计方案(不基于严格数学意义上的巻积)使用预先计算散射核函数，所述散射核函数取决于在所述像素处的对象衰减和 从投影图像导出的进一步性质(例如总对象大小或它的最大深度的估计值， 或在所述像素处的水当量图像中的衰减梯度的估计值)，确定从X射线源到 探测元件的射线的散射分布。可以通过总计所有这样的射线的贡献获得总 散射图像。可以用实验方法或数值地生成所述核函数。这些核函数是可使 用的以便对投影图像,)(x,力进行误差补偿，所述投影图像由初级部分P和 散射部分S组成，艮口， "(x,力^(;c,;;) + S0c，力。
该重建方法可以用于断层摄影装置的领域，例如计算机断层摄影装置， 特别是X射线计算机断层摄影术。
在下面，将描述用于生成核函数集的方法的进一步典型实施例。然而， 这些实施例也适用于用于误差补偿的方法和系统，断层摄影装置，计算机 可读介质和程序单元。
根据用于生成核函数集的方法的另一典型实施例，通过使用X射线模
体(phantom)作为模型用实验方法确定核函数集。特别地，在核函数集的 计算中可以使用实验测量的结果。
根据用于生成核函数集的方法的另一典型实施例，通过使用对几何模 型的散射模拟，优选地假设类水散射特性或其他材料的散射特性，计算核 函数集。优选地，核函数集的每个核函数是几何模型的参数的函数。
也就是说，为了生成预先计算的核函数，可以使用几何模型的笔形束 蒙特卡罗(Monte-Carlo)散射模拟离线生成归一化散射分布，所述散 射分布可以以这样的方式被参数化，即使得它考虑校正系统几何形状，例 如断层摄影系统的几何形状，校正束波谱，例如断层摄影系统的相应辐射 源的能量谱，和正确的防散射滤线栅，例如防散射滤线栅是否用在断层摄 影系统中和哪个特定防散射滤线栅用在断层摄影系统中。通过总计撞击在 投影图像中的单个探测器像素上的射线的预先计算的贡献，从这些散射分 布可获得散射图像W(u)的估计值。
根据用于生成核函数集的方法的另一典型实施例，至少一个参数是几 何模型的半径。优选地，所述核函数进一步是所述几何模型的投影中心与 穿透笔形束撞击在探测器上的位置之间的移位的函数。
也就是说，对于给定的系统配置，例如断层摄影系统配置，可以作为
ii模型参数M例如至少一个半径的函数，和作为模型相对于用于模拟的笔形 束的位置移位(r,(D)的函数离线生成独立核函数^^。(;c，力，其中(r,O)是表示 在平行于断层摄影系统的探测器平面的平面中的移位的极坐标。通过作为 位置移位(。①)的函数计算核函数，有可能取决于像素位置，例如像素是边 界像素还是中心像素，而解决散射变化。
根据用于生成核函数集的方法的另一典型实施例，所述几何模型是椭 球体模型。优选地，核函数集的每个核函数是几何模型的^2和,*3以及模型 的中心与笔形束穿透模型的位置之间的移位r,①的函数，所述移位可能导致 几何模型的投影中心与穿透笔形束撞击在探测器上的位置之间的移位，其 中和r3可以是椭球体模型的半轴。
对于这些模型参数M例如^^和^可以计算预先计算的核函数，艮口， 可以作为这些模型参数M的函数和作为模型相对于笔形束的位置移位(r,O)) 的函数计算预先计算的核函数^c，力。对于模型参数M^^^)的每个组合， 可以在笔形束与平行于探测器的平面中的椭球体模型之间的相对位置的变 化下生成散射核函数A^，^Oc,力，其中椭球体模型相对于笔形束的位置移位
由极坐标o,①)表示。
根据用于生成核函数集的方法的另一典型实施例，所述几何模型是球 体模型。优选地，核函数集的每个核函数是球体模型的半径i 的和模型的中 心与笔形束穿透模型的位置之间的移位r,①的函数，所述移位可能导致几何 模型的投影中心与穿透笔形束撞击在探测器上的位置之间的移位。
在下面，将描述用于误差补偿的方法的进一步典型实施例。然而，这 些实施例也适用于用于生成核函数集的方法，用于误差补偿的系统，断层 摄影装置，计算机可读介质和程序单元。
根据用于误差补偿的方法的另一典型实施例，原始投影图像被归一化。
根据用于误差补偿的方法的另一典型实施例，核函数集的每个核函数 是成像系统的几何形状、成像系统的射束波谱和/或成像系统的防散射滤线 栅参数的函数。
在该上下文中"归一化"表示这样的事实，即，数量？表示由空气值归
一化的初级辐射的探测强度，使得/> = 1对应于直接辐射并且^ = 0对应于完全
吸收。通过归一化投影图像和将它转换成水当量图像，有可能通过使用预先计算的核函数提供误差补偿图像的高效方式。
根据用于误差补偿的方法的另一典型实施例，根据7^,力=—ln(P(°)(W》
将原始投影图像转换成水当量图像，其中,)代表原始投影图像，roc,力代 表水当量厚度r的图像，并且;/代表水的适当衰减值。
特别地，生成水当量图像可能是合适的，原因是水在人体中占多数， 因此导致简单但仍然足够正确的模型。
根据另一典型实施例，用于误差补偿的方法还包括通过总计与所有像 素对应的所有核函数的贡献计算像素阵列中给定像素处的总散射。
对所有核函数的所有贡献进行的这种总计是用于计算总散射贡献对在 给定像素处测量的强度，例如在计算机断层摄影装置的探测器元件处测量 的强度的影响的高效方式。该总散射贡献以后可以用于补偿通过散射引入 投影图像的误差。
根据用于误差补偿的方法的另一典型实施例，在给定像素处的总散射
由下式定义
其中W(/,刀是在像素(/,y)处的总散射，w表示像素的面积，并且
-w-/)是核函数，其指示从在像素(/,力的位置处撞击像素ou)
的射线引入的散射并且取决于代表几何模型的参数的M ，和代表几何模型
相对于像素阵列的中心的位置移位的(。①)。
根据用于误差补偿的方法的另一典型实施例，使用椭球体模型计算核
函数，其中似代表椭球体模型的半轴^,2,3，其中/"^r^wrt04/兀)并且,'^S，
」=物理对象的最大横截面积，并且5=物理对象的最大厚度。
相应的参数/!^和^可以从物理图像的投影图像提取并且可以代表椭球
体模型的半轴。优选地，根据r(:c,力=—ln(pW(^))将投影图像转换成水当量
厚度r的图像，其中卢代表原始投影图像，其中r(u)代表水当量厚度r的 图像，并且//表示水的适当衰减值。参数^可以被指定为在投影图像上的物 理对象的阴影，例如衰减超过某个阈值的投影区域的面积，除以系统的几何放大因数的平方。参数5可以被指定为低通滤波之后的roc,力的最大值或 来自r的直方图的百分位数，两者都可以使强衰减变化的影响最小化。在一
个备选实施例中，可以从模型的正投影到采集投影的最小平方拟合确定模 型参数。
在该上下文中应当注意的是对于给定模型M，不同移位值相当于在笔
形束位置的水当量厚度的不同值，从在零移位的模型的最大厚度向下变化 到在几乎等于模型的空间范围的移位处的几乎零厚度。转而，对于简单几
何模型和固定移位角①，在考虑范围内的水厚度r的给定值明确地确定r的 相应值，使得在区间(o,rmax]中可以假设Kr)为唯一方程。
对于该实施例，考虑撞击在带有标记(W)的探测器像素上的射线的散射
贡献。在另一像素(/,力的位置，该射线产生由表达式&,^,))，，"(/",卜/)大 致描述的散射贡献，其中对于所利用核函数，位移半径r由在像素(w)处的 水厚度指定，并且移位角o(/U)可以在以"衰减质量的中心"(q,o作为原
点的坐标系中被选择作为像素(W)的极角，(c,，c2)被指定为
^4^.2^(w)ffl。然后可以通过总计所有射线(w)的贡献获得在像
〉j7W) A/ 、/乂
素o)处的总散射，产生^(/，力=^5;&禅，/))，,/)(卜^-/)，其中求和在探测
器的所有像素(w)上进行，并且w表示像素面积。
根据用于误差补偿的方法的另一典型实施例，使用球体模型计算核函
数，其中在给定像素处的总散射由下式定义
其中是在像素(/，y)处的总散射，w表示像素面积，并且
&卵,^ ,，稀力)婦力(''-0是核函数，其指示从在像素(!',力的位置处撞 击像素OU)的射线引入的散射并且取决于代表球体几何模型的半径的J ，代 表水当量厚度r的相应图像的梯度的g，和代表椭球体几何模型相对于像素
阵列的中心的位置移位的(r,cD)。
r
根据用于误差补偿的方法的另一典型实施例根据/^t.V^7和
14^丄.g和d^arg(gradr)选择参数/ 和r，其中T =物理对象的水当量厚度，并 4
根据该典型实施例，可以使用球体几何模型，其可以具有的显著优点 在于它并不需要对每个投影估计全局模型参数，而是基于对每个单一射线 局部估计这样的参数。该变型使用球体几何模型(模体)，并且如先前所述 的变型一样，也通过相对于笔形束的模体偏移工作。
应用于投影P ，所述方法首先可以要求计算水当量厚度r:-(ln。/^的 相应图像的梯度，对于每个探测器元件所述梯度具有某个幅度值g叫gmdrl 和方向①-arg(gradr)。为了估计给定源射线的散射贡献，水当量厚度r的局 部值，梯度幅度g和方向①，然后可以唯一地确定所利用的球体模体的参数 (i ,。cD)，其中i 可以表示球体的半径，并且(r,d))可以是它在平行于探测器 的平面中的位置偏移。以这样的方式进行映射(r，g)-(i 力，即使得所利用
的球体的平行投影将具有水当量厚度r和在笔形束的位置处的厚度梯度g。
这通过以下公式实现i^二.V^7， "二.g。应当注意的是以该方式，位4 4
置偏移在平图像区域中将接近于零，同时它在陡梯度，例如在对象边界附 近变成接近于球体半径。使用该方法，对于给定系统几何形状和射束质量， 与在基于椭球体模型(核函数)的方法的情况下的四个参数相比，取决于 三个参数(凡,',①)预先计算巻积核函数。
当通过截断影响投影图像时，例如在物理对象大于可能的成像对象的 情况下，这样的球体模型可能是特别高效的。尽管椭球体模型可以由于该 截断受到模型参数;!-^的错误估计的影响，基于球体核函数的方法由于它 的模型参数的局部估计可以受到截断影响。
根据另一典型实施例，所述方法进一步包括通过使用总散射以乘法方
式计算第一误差补偿图像。优选地，根据卢"+1)=^^执行乘法校正，其
中W表示从投影图像估计的散射图像。
乘法方式可以是特别有利的，原因是它可以表现出增加的收敛稳定性 并且可以具有的附加优点在于避免了负投影值。使用后一校正方案，假设估计的散射相同，在带有高衰减的区域中与带有低衰减的区域相比更少的 散射量可以被校正。与最大被减散射量的预定阈值可以被指定以避免负投 影值的减法校正相比，可以使用乘法校正自动避免这样的影响。与减法校 正相比，乘法校正可能需要在全分辨率投影图像上执行，并且因此在每次 迭代中在应用校正步骤之前可以对估计的粗略散射分布再次至少部分上采 样。
根据另一典型实施例，所述方法进一步包括通过使用总散射以减法方 式计算第一误差补偿图像。优选地，根据卢"+1)=,°)-S(")执行减法校正，其 中表示从投影图像估计的散射图像。
根据另一典型实施例，所述方法进一步包括通过将第一误差补偿图像 用作投影图像计算第二误差补偿图像。也就是说，可以以迭代方式，例如 在4一5次重复中执行校正。也就是说，在散射图像，)(x,力的估计之后，然 后该图像用于校正原始采集的投影图像,V,力(包含初级和散射辐射两者 的贡献)，从而产生真实初级图像的估计值卢"Oc,力。由于初始散射恶化
(scatter deteriorated)的投影图像")导致稍稍失真的厚度图像r ，估计和校 正步骤优选地以迭代方式重复多次，直到达到估计的初级图像的收敛(这 通常在大约4一5次迭代中获得)。由于散射分布是平滑的，可以使用强下 采样探测器像素栅格执行散射估计以便减小计算工作量。
可以看出本发明的一个典型方面在于在核函数生成期间引入了所利用 模体的可变偏移。基于椭球体核函数和基于球体核函数的方案利用这样的 偏移，并且因此潜在地能够适当地解决在对象边界附近产生的散射分布的 不对称。两个估计方案可以具有应用于X射线体积成像的高潜力。尤其基 于预先计算的球体核函数的方案可以产生不同身体区域(例如头，胸和骨 盆区域)的准确结果，并且它的性能可以不受截断的存在影响。最重要地， 这些身体区域的最佳校正因数几乎可以是相同的。关于计算成本，球体方 法可能比椭球体方法稍稍要求高些，原因是所有可能球体配置的散射核函 数优选地被读取并且同时被存储在存储器中。对于该方法的最高效使用， 该数据可以保持在存储器中而不是当该方法应用于旋转采集的投影序列时 重复被读取。
为了改善基于椭球体核函数的方法，当应用于胸和骨盆的投影时所述
16方法可能受到截断的影响，有可能使用模型的正向投影通过优化算法更稳 健地估计模型参数，所述优化算法也可独立地应用于每个采集投影。这是 由于这样的事实，g卩，该方法至少依赖于每个投影的两个全局参数的近似 估计，其中的一个全局参数在截断的情况下难以估计。
此外，根据本发明的一个典型方面考虑了散射补偿的校正步骤的两个 不伺方案，目卩，减法和乘法校正。每个方案可以与根据典型实施例的每个 散射估计算法组合。尽管减法校正可能偶尔产生削波相关条纹伪影并且可 能遭受迭代估计校正过程的令人不满意的稳定性，它可以有较少的计算时 间消耗。备选地，乘法校正可以在所有情况下产生有利结果。由于乘法校 正可能需要在较高分辨率投影图像上被执行，在每次迭代中在应用校正步 骤之前可以对估计的粗略散射分布再次上采样。
可以通过计算机程序，即通过软件，或者通过使用一个或多个专用电 子优化电路，即在硬件中，或者以混合形式，即通过软件部件和硬件部件， 实现对物理对象的投影图像的误差补偿。计算机程序可以以任何合适的编
程语言例如0++写入并且可以存储在计算机可读介质例如CD-ROM上。而
且，计算机程序可从网络获得，例如万维网，它可以从所述网络下载到图 像处理单元或处理器，或任何合适的计算机。
在该上下文中应当注意的是本发明并不限于计算机断层摄影，而是可
以包括使用基于C形臂的3D旋转X射线成像，正电子发射断层摄影等。 也应当注意的是该技术特别可用于对患者的不同身体区域例如头、胸或骨 盆区域进行医学成像。
从和参考下文描述的实施例将显而易见和阐明本发明的这些和其他方 面。在该申请的任何地方描述的公开实施例和方面可以彼此混合和/或组合。
在下面将参考下图描述本发明的典型实施例。


图1显示了计算机断层摄影系统的简化示意图2显示了用于生成椭球体核函数的几何图形的示意草图3显示了根据本发明的典型实施例的误差补偿方法的示意性流程图;
图4显示了一些典型散射图像。图中的图示是示例性的。在不同的图中，类似或相同的元件带有类似 或相同的参考符号。
图1显示了计算机断层摄影系统的一个典型实施例，其中投影数据可 以由根据本发明的一个实施例的校正方法处理。
图1中描绘的计算机断层摄影系统100是锥束CT扫描器。图1中描绘
的CT扫描器包括可围绕旋转轴102旋转的扫描架101。扫描架101借助于 电机103驱动。参考数字105指示辐射源例如X射线源，其发射多色或单 色辐射。
参考数字106指示孔径系统，其将从辐射源单元发射的辐射束成形为 锥形辐射束107。锥束107被引导使得它穿透布置在扫描架101的中心，即 在CT扫描器的检查区中的感兴趣对象110，并且撞击在探测器115 (探测 单元)上。从图1可以看出，探测器115与辐射源单元105相对地布置在 扫描架101上，使得探测器115的表面由锥束107覆盖。图1中描绘的探 测器115包括多个探测元件115a，每个能够探测由感兴趣对象110散射、 衰减或穿过其中的X射线。图1中示意性显示的探测器115是二维探测器， 即，单个探测器元件被布置在平面中，这样的探测器被用在所谓的锥束断 层摄影术中。
在扫描感兴趣对象110期间，辐射源单元105，孔径系统106和探测器 115沿着扫描架101沿着由箭头117指示的方向旋转。为了使带有辐射源单 元105，孔径系统106和探测器115的扫描架101旋转，电机103连接到电 机控制单元120，所述电机控制单元连接到控制单元125 (其也可以被指示 并且用作计算、重建或确定单元)。
在图1中，感兴趣对象IIO是设置在操作台112上的人。在对头110a， 胸或人110的任何其他部分的扫描期间，当扫描架101围绕人110旋转时， 操作台112可以沿着平行于扫描架101的旋转轴102的方向移动人110。这 可以使用电机113进行。通过这样，头沿着螺旋扫描路径被扫描。操作台 112也可以在扫描期间停止，由此测量信号切片。
探测器115连接到控制单元125。控制单元125接收探测结果，即，来 自探测器115的探测元件115a的读数，并且在这些读数的基础上确定扫描 结果。此外，控制单元125与电机控制单元120通信以使带有电机103和113的扫描架101的运动与操作台112协调。
控制单元125适于从探测器115的读数重建图像。由控制单元125生 成的重建图像可以通过接口输出到显示器(未在图1中显示)。
控制单元125可以由数据处理器或计算机实现以处理来自探测器115 的探测器元件115a的读数。
图1中所示的计算机断层摄影装置可以采集患者的头或胸的计算机断 层摄影数据。换句话说，当扫描架101旋转并且当操作台112线性移位时， 然后相对于患者由X射线源105和探测器115执行螺旋扫描。在采集这些 数据之后，将数据转送到控制单元125，并且对测量数据进行回顾性分析。
图2显示了用于生成椭球体核函数的几何图形的示意草图。通过参考 该草图，将描述带有可变偏移的椭球体核函数的典型实施例。该方法解决 了如下事实，即，源自靠近被成像对象的边界的区域(例如对象中心相比 于边界区域)的散射分布是高度不对称的，而在已知方法中在散射射线和 模型之间不使用偏移，所以产生的散射分布的不对称通常未被准确地解决。
根据椭球体模型，投影图像POc，力首先被归一化并且然后根据公式 r = -(lnp)/一被转换成水当量厚度r(x,力的等效图像，其中//表示水的近似衰 减值。
然后，从水当量厚度r的图像提取两个标量，用类水衰减和散射特性指 定被成像对象的椭球体模型的参数。特别地，假设均匀椭球体具有在平行
于探测器表面的平面中的半轴/",^-sqrt04Z;r)，和垂直于探测器的半轴
r3=rmax/2。在这里，^是平行于探测器表面的被成像对象的横截面积的测量
值并且被指定为对象阴影(其被定义为带有的水当量厚度在某个阈值例如
10mm之上的投影中的区域)的面积除以系统的几何放大因数的平方。数量 T^是最大水当量厚度的近似测量值。为了计算解决了对像素位置的重要依 赖性的散射核函数，已知的水厚片由椭球体模型替换，并且附加地考虑模 型相对于模拟笔形束的位置偏移。对于模型参数M^n, "的每个组合，在 平行于探测器的平面中在笔形束与椭球体模型之间的相对位置的变化下生 成散射核函数^^，。(^)，其中椭球体模型相对于笔形束的位置移位由极坐 标(/",①)表示。在该上下文中应当注意的是对于给定模型M，不同移位值相 当于在笔形束位置的水当量厚度的不同值，从在零移位处的模型的最大厚
19度向下变化到在几乎等于模型的空间范围的移位处的几乎零厚度。转而，
对于固定移位角f，在考虑范围内的水厚度r的给定值明确地确定r的相应 值，使得在区间(o，7;j中Kr)为唯一方程。在图2中示出了用于生成椭球体 核函数的几何图形。
图2显示了包括列x和行y的平板探测器201。在探测器上探测到的穿 透椭球体模型的射线的散射分布示意性地由平板探测器201上的白区描绘。 可以看出探测到的散射是高度不对称的，即，散射对平板探测器201的左 半部的影响比对平板探测器201的右半部的影响远远更高，这导致留下离 中心的更亮像素。进一步地，图2在示意性形状中显示了用于生成椭球体 核函数^^(x,力的水椭球体202。椭球体202由几个参数特别是半轴/",203， 。204表征，而^未在图2中被描绘，原因是它垂直于图2中所示的平面延 伸。此外，在图2中标示了非零位置移位r205， g口，从焦斑207延伸到平 板探测器201的中心的焦线206和水椭球体202的中心之间的非零移位。 与^相同，未在图2中显示移位角0，原因是它在平行于探测器表面的平面 中被定义。水厚度r在图2中被描绘成208，而线209示意性地显示了不同 散射射线。
使用模型和位置相关核函数^^，。(x,力，在另一像素(/,力的位置处撞击 在带有标记(W)的探测器像素上的射线的散射贡献由表达式 - W - /)给出，其中所利用的核函数的位置移位的长度卩通过在 像素(W)处的水厚度指定，并且移位角C^,/)可以在以"衰减质量的中心"
(…&)作为原点的坐标系中被选择作为像素(W)的极角，(c,^)被指定为
A ) i______ "、
、、
xnw)H。这可以在单一、近似椭球体形状的对象的情况下
2/(W)" 、"
提供不对称散射核函数分布的合适定向。然后通过总计所有射线(w)的贡献
获得在像素a力处的总散射，从而产生,(/,力=^2>乾仰,,)脾,,々'-^-/)。
4.乂
除了r,^^^-i 之外，即代替水椭球体使用水球体，用于生成球体核
函数的几何图形与图2中所描绘的相同。然而，可以不同地计算偏移。图3显示了根据本发明的典型实施例的误差补偿方法的示意性流程图。 该实施例特别涉及椭球体几何模型。该方法独立地处理每个采集的投影图 像并且可以包括以下序列
1.根据7" = -(lnP(。)(u))/z/将由初级部分P和散射部分S组成的采集的 归一化投影图像")(W)^(x,力+ S(;c,力转换成水当量厚度r的图像，其中一 表示水的适当衰减值(步骤S301)。
2.从图像r提取许多标量参数，指定被成像对象的简单几何模型的 参数。例如，可以根据例如^和S为标量的/","2-sqrt(^/7r)和r^5计算带有 平行于探测器平面的圆形横截面和类水衰减的均匀椭球体对象模型的半轴 r.。在该特定实施例中，即在椭球体情况下，^是平行于探测器表面的被成 像对象的最大横截面积的近似测量值，并且S是被成像对象的最大水当量厚 度的适当测量值。^可以被指定为对象阴影，例如衰减在某个阈值之上的投 影区域的面积除以系统的几何放大因数的平方。为了最小化局部强衰减变
化的影响，s可以被指定为低通滤波之后r(x,力的最大值或来自r的直方图 的百分位数。在一个备选实施例中，从模型的正向投影到采集投影的最小 平方拟合确定模型参数(步骤302)。
3.通过总计撞击在单个探测器像素上的射线的预先计算的贡献获得 散射图像W(x,力的估计值。为了该目的，考虑正确的系统几何形状、射束波 谱和防散射滤线栅参数，使用参数对象模型的笔形束蒙特卡罗
(Monte-Cario)散射模拟离线生成归一化散射分布X(x,力。对于给定系统配 置，作为模型参数M的函数，和为了解决像素位置相对于投影对象中心的 重要依赖性作为模型相对于笔形束的位置移位(。o)的函数离线生成独立核 函数4,,。"力，其中"O)是表示在平行于探测器的平面中的移位的极坐标。 重要的是应当注意对于给定模型M，不同移位值相当于在笔形束位置的水 当量厚度的不同值，从在零移位处的模型的最大厚度向下变化到在几乎等 于模型的空间范围的移位处的几乎零厚度。转而，对于简单几何模型和固
定移位角cD ，在考虑范围内的水厚度r的给定值明确地确定了 r的相应值，
使得在区间(o,rmaj中，可以假没k:o为唯一方程。
现在，考虑了撞击在标记为(w)的探测器像素上的射线的散射贡献。在
21另一像素(/,_/)的位置，该射线产生了由表达式/^(歌,),錄,)(/"，卜/)大致描述 的散射贡献，其中对于所利用的核函数，移位半径r由在像素(W)处的水厚
度指定，并且移位角O(W)可以在以"衰减质量的中心"(C,2)作为原点的
坐标系中被选择作为像素(W)的极角，(c',c2)被指定为
^~一.5>(w).f"。然后可以通过总计所有射线ou)的贡献获得在像
素(/,_/)处的总散射，从而产生S。(/,力=w. - W -/)，其中求和在
探测器的所有像素(W)上进行，并且w表示像素面积(步骤303)。
4.使用估计的散射f)(x，力，然后校正所采集投影图像")(x,力，从 而产生真实初级图像的估计值,V,力(步骤304)。由于初始散射恶化的投 影图像,)导致稍稍失真的厚度图像r，当以迭代方式重复1至4 (步骤301 至304)大约四次直到达到估计的初级图像的收敛时获得最佳结果，其中重 复2次(步骤302)是最佳的。由于散射分布是平滑的，可以使用强下采样 探测器像素栅格执行以上步骤以便减小计算工作量。
可以以减法或乘法方式执行校正。根据公式卢"+"-,)-，执行n次迭
代减法校正(n^1)。然而，发现根据= :、 的乘法校正增加了收
敛的稳定性并且具有的附加优点是避免了负投影值。
图4显示了一些典型散射图像。在图2的上图中特别地在椭球体模型 中显示了本发明的一个典型实施例的误差校正方法的结果，而图2的下图 显示了基于预先计算散射核函数的已知方法的结果，其并不使用模型的位 置偏移并且因此并未准确地解决尤其靠近对象边界的射线的不对称散射贡 献。具体地，图4a在上部分中显示了针对具有行y和列;c的二维探测器描 绘的估计的散射图像。在图4a的下部分中显示了沿着通过图像的中心水平 横截面的相应轮廓。图4b在上部分中显示了针对具有行y和列;c的二维探 测器描绘的估计的散射图像。在图4b的下部分中显示了沿着通过图像的中 心水平横截面的相应轮廓。图4c在上部分中显示了针对具有行y和列x的二维探测器描绘的模拟
背景实况(ground truth)。在图4a的下部分中显示了沿着通过图像的中心 水平横截面的相应轮廓。图4d描绘了已知方法的相同图像。图4c和4d是 相同的，原因是将相同的背景实况用于比较。
图4e在上部分中显示了在行y和列;c的二维绘图中估计图像与背景实 况的比率。可以清楚地看出平均比率大约为1，但是仍然显示了大约5%的 稍稍高估，而散射形状被良好地近似，这可以从灰度级值的相对均匀分布 看出。相反地，图4f显示了高度不均匀分布。特别在图像中心，对应于对 象的最大厚度的区域，散射被大大高估，而在边界附近高估要小得多。高 估平均大约为44%。
在测试的实现方式中，根据本发明的一个典型实施例的校正方法的计 算工作量仅仅是己知的基于巻积的方法的大约两倍高，其结果在图4的下 部分中被描绘。 一般而言根据本发明的一个典型实施例的校正方法可以允 许潜在地远远更准确地估计散射分布，这是由于这样的事实，即，散射对 进一步参数而不仅是水当量厚度的依赖性得以解决。应当注意的是词语"包 括"并不排除其他元件或步骤并且"一"并不排除多个。而且结合不同实 施例描述的元件可以组合。应当注意的是权利要求中的参考符号不应当被 理解成限制权利要求的范围。
权利要求
1. 一种用于生成核函数集的方法，所述核函数集用于对由成像系统记录的物理对象的投影图像进行卷积误差补偿，所述方法包括以如下方式计算所述核函数集，即使得对于所述投影图像中的每个像素计算用于误差补偿的不对称散射分布，从而表示在由X射线源和所述像素之间的射束定义的体积中产生的X射线散射。
2. 根据权利要求1所述的方法，其中，使用X射线模体作为模型用实验方法确定所述核函数集。
3. 根据权利要求1所述的方法，其中，使用几何模型的散射模拟计算所述核函数集。
4. 根据权利要求2或3所述的方法，其中，所述核函数集的每个核函数是所述几何模型的参数的函数。
5. 根据权利要求4所述的方法，其中，所述参数中的至少一个是所述几何模型的半径。
6. 根据权利要求4或5所述的方法，其中，所述核函数进一步是所述几何模型的投影中心与穿透笔形束撞 击在探测器上的位置之间的移位的函数。
7. 根据权利要求2—6中任一项所述的方法， 其中，所述几何模型是椭球体模型。
8. 根据权利要求7所述的方法，其中，所述核函数集中的每个核函数是所述几何模型的。^和6的以及 所述模型的中心与所述笔形束穿透所述模型的位置之间的移位r,(D的函数。
9. 根据权利要求2—6中任一项所述的方法，其中，所述几何模型是球体模型。
10. 根据权利要求9所述的方法，其中，所述核函数集中的每个核函数是所述球体模型的半径i 的和所述 模型的中心与所述笔形束穿透所述模型的位置之间的移位;",O的函数。
11. 根据前述权利要求中任一项所述的方法，其中，所述核函数集中的每个核函数是所述成像系统的几何形状、所 述成像系统的射束波谱和/或所述成像系统的防散射滤线栅参数的函数。
12. —种用于对物理对象的图像进行误差补偿的方法，所述方法包括接收被成像物理对象的原始投影图像；将所述原始投影图像转换成水当量图像，特别是计算相应的梯度图像； 从所述水当量厚度的图像并且特别从所述梯度图像提取许多参数； 根据权利要求l一ll的任何一个通过使所提取的参数与预先计算的核函数的参数相关来确定至少一个预先计算的核函数；和通过使用所述确定的至少一个预先计算的核函数补偿所述原始投影图像的误差。
13. 根据权利要求12所述的方法，其中，归一化所述原始投影图像。
14. 根据权利要求12或13所述的方法，其中，根据r(x,力=—ln(iJ(。w))将所述原始投影图像转换成水当量图像，其中，,)代表所述原始投影图像， r(;c,力代表所述水当量厚度r的图像；并且A代表水的适当衰减值。
15. 根据权利要求12—14中任一项所述的方法，进一步包括 通过总计对应于所有像素的所有核函数的贡献计算在像素阵列中的给定像素处的总散射。
16. 根据权利要求15所述的方法，其中，在所述给定像素处的总散射由下式定义 "("力=w. $ ^^亂，(w) ('■",_/-,)，其中W(!',力是在像素处的总散射， ^表示像素的面积，并且^W雄^U-/)是核函数，其指示从在所述像素(U)的位置处撞击 像素(W)的射线引入的散射并且取决于 代表所述几何模型的参数的M;和代表所述椭球体几何模型相对于所述像素阵列的中心的位置移位的 C 。
17. 根据权利要求16所述的方法，其中，使用椭球体模型计算所述核函数；并且 其中，M代表所述椭球体模型的半轴。^,^ 其中，。-r2-s^^04/;r)并且^-S，4 =所述物理对象的最大横截面积，并且 ^ =所述物理对象的最大厚度。
18. 根据权利要求15所述的方法，其中，使用球体模型计算所述核函数；并且 其中，在给定像素处的总散射由下式定义<formula>formula see original document page 4</formula> 其中w(/，刀是在像素(/,力处的总散射，w表示所述像素的面积，并且脾，,々'-W -0是核函数，其指示从在所述像素的 位置处撞击像素(/M)的射线引入的散射并且取决于-代表所述球体几何模型的半径的i ; 代表所述水当量厚度r的相应图像的梯度的g，和 代表所述椭球体几何模型相对于像素阵列的中心的位置移位的(r①)。
19. 根据权利要求18所述的方法，其中，J4 + g2禾口r^二.g禾口0^arg(gradr)， 4 4r二所述物理对象的水当量厚度，并且g = |gradr| 。
20. 根据权利要求16—19中任一项所述的方法，进一步包括 通过使用所述总散射以乘法方式计算第一误差补偿图像。
21. 根据权利要求20所述的方法，进一步包括根据户("+') = V:执行所述乘法校正，其中 ^表示从所述投影图像,)估计的散射图像。
22. 根据权利要求16—19中任一项所述的方法，进一步包括 通过使用所述总散射以减法方式计算第一误差补偿图像。
23. 根据权利要求22所述的方法，进一步包括 根据尸("+1) = sw执行所述乘法校正，其中表示从所述投影图像卢")估计的散射图像。
24. 根据权利要求20—23中任一项所述的方法，进一步包括-通过将所述第一误差补偿图像用作所述投影图像来计算第二误差补偿图像。
25. —种用于对物理对象的图像进行误差补偿的系统，所述系统包括 接收单元，其适于接收被成像物理对象的原始投影图像； 计算单元，其适于将所述原始投影图像转换成水当量图像，特别是计算相应的梯度图像，并且从所述水当量厚度的图像并且特别从所述梯度图 像提取许多参数；确定单元，其适于根据权利要求l一ll的任何一个通过使所提取的参数与预先计算的核函数的参数相关来确定至少一个预先计算的核函数；和补偿单元，其适于通过使用所述确定的至少一个预先计算的核函数补 偿所述原始投影图像的误差。
26. —种断层摄影装置，其包括 辐射源；辐射探测器；和用于根据权利要求25进行误差补偿的系统；其中，所述辐射探测器适于记录代表所述被成像物理对象的所述原始 投影图像的数据。
27. —种计算机可读介质，用于对物理对象的图像进行误差补偿的程 序存储在所述计算机可读介质中，当由处理器执行时，所述程序适于控制 一种方法，所述方法包括接收被成像物理对象的原始投影图像；将所述原始投影图像转换成水当量图像，特别是计算相应的梯度图像;从所述水当量厚度的图像并且特别从所述梯度图像提取许多参数；根据权利要求l一ll的任何一个通过使所提取的参数与预先计算的核 函数的参数相关来确定至少一个预先计算的核函数；和通过使用所述确定的至少一个预先计算的核函数补偿所述原始投影图 像的误差。
28. —种用于对物理对象的图像进行误差补偿的程序单元，当由处理 器执行吋，所述程序适于控制一种方法，所述方法包括接收被成像物理对象的原始投影图像；将所述原始投影图像转换成水当量图像，特别是计算相应的梯度图像;从所述水当量厚度的图像并且特别从所述梯度图像提取许多参数；根据权利要求l一ll的任何一个通过使所提取的参数与预先计算的核 函数的参数相关来确定至少一个预先计算的核函数；和通过使用所述确定的至少一个预先计算的核函数补偿所述原始投影图 像的误差。
全文摘要
一种用于生成核函数集的方法，所述核函数集用于对由成像系统记录的物理对象的投影图像进行卷积误差补偿，所述方法包括以如下方式计算所述核函数集，即使得对于所述投影图像中的每个像素计算用于误差补偿的不对称散射分布，从而表示由X射线源到所述像素的射束产生的X射线散射。
文档编号G06T5/00GK101473348SQ200780022897
公开日2009年7月1日 申请日期2007年6月13日 优先权日2006年6月22日
发明者J·维格特, J·蒂默, M·贝尔特拉姆, N·J·努尔德霍尔克 申请人:皇家飞利浦电子股份有限公司