CN1820710A

CN1820710A - 预测生命体内造影剂流动的方法

Info

Publication number: CN1820710A
Application number: CNA200610007032XA
Authority: CN
Inventors: 厄恩斯特·克洛茨; 安娜贝拉·劳舍尔
Original assignee: Siemens AG
Current assignee: Siemens AG
Priority date: 2005-02-14
Filing date: 2006-02-14
Publication date: 2006-08-23
Also published as: DE102005006657A1; US20060239917A1; JP2006223861A; US7477929B2

Abstract

本发明涉及一种预测造影剂在生命体内、尤其是在患者体内的流动的方法，其中：将包含具有公知注射流动变化过程的造影剂的特定测试药丸注入体内，优选注入血管、特别是静脉中，利用断层造影方法在包含多个测量时刻的有限时间段Z期间观察和确定造影剂在体内至少一个地点的时间浓度变化过程，从获得的关于造影剂分布的测量数据中，利用线性原因/效果依据来预测另一次造影剂注入的造影剂浓度的时间变化过程，为了预测造影剂在体内至少一个事先测量过的地点x的时间浓度变化过程_R(t)而采用以下计算公式(见图下)。

Description

预测生命体内造影剂流动的方法

技术领域

本发明涉及一种预测造影剂在生命体内、尤其是在患者体内的流动的方法，其中将包含具有公知注射流动变化过程的造影剂的确定测试药丸注入体内，优选注入血管、特别是静脉中，利用断层造影方法在包含多个测量时刻的有限时间段Z期间观察和确定造影剂在体内至少一个地点的时间浓度变化过程，并从获得的关于造影剂分布的测量数据中利用线性原因/效果依据预测另一次造影剂注入的造影剂浓度的时间变化过程。

背景技术

对于断层造影方法，尤其是在计算机断层造影或NMR断层造影领域的断层造影方法，由于在对特定身体区域显示中出现的小对比度，因此为了显示该特定身体区域优选施加造影剂，由此获得该身体区域的对比度较强的图像。但造影剂大多具有缺点，即其从生物学观点来看难以吸收，因此其剂量应当保持尽可能地小。但由于被检查体的生物可变性，无法对一种特定造影剂在体内的观察点的浓度变化过程是如何随着时间发展的得出足够精确和适用于普遍情况的结论。因此对于被检查体来说必须基于造影剂或测试药丸的测试注射来观察该测试注射的效果，尤其是注射之后浓度值在被检查体内的感兴趣地点的时间变化过程。

在结合CT检查的应用中，在测试药丸注射后的一段特定时间内，并且在采用尽可能少的辐射剂量的条件下，间接通过出现的HU值的变化来测量浓度。由于只对造影剂的成像作用感兴趣，并且在造影剂的成像效果和浓度之间形成线性关系，因此关于造影剂的绝对浓度的结论还不知道而且也是次要的。这种测试检查的图像的分辨率也是很低的。

公知基于对这种测试药丸的效果的认识利用傅立叶变换预先计算正式的造影剂注入效果，并由此确定必须的造影剂剂量。

发明内容

本发明要解决的技术问题是，找到一种改进的预测方法，可以对检查特定身体区域所必需的造影剂剂量作出更准确的预测。

发明人已经具有如下认识：

为了对应于上述技术问题找到一种更好的方法，该方法从测试药丸注射之后的增强值的时间变化过程中预测在另一次不同的注射之后在相同地点的增强值的时间变化。这应用于优化的造影剂方案规划。在此增强值应理解为造影剂在图像显示中的浓度改变的效果。在CT检查中增强值例如是所确定的HU值。

测试药丸的变化过程可以通过以下线性函数表示：

b_T(t)＝F_TΘ(t-t_0T)Θ(t_FT-t)(1)

其中F_T表示流速，t_0T和t_FT表示测试药丸b_T(t)的开始时间和结束时间。在此Θ(t)表示如下定义的海维希德函数或所谓的阶梯函数：

Θ (t) = \{\begin{matrix} 1 & t > 0 \\ 0 & t < 0 \end{matrix} \cdot \cdot \cdot (2) .

已消耗的造影剂的量作为F_T(t_FT-t_0T)给出。增强值c_T(t)基于测试药丸的注入在预定位置的变化过程通过下式给出

c_{T} (t) = c_{0} + {\tilde{c}}_{T} (t) \cdot \cdot \cdot (3)

其中，在CT检查的情况下c₀表示没有注入造影剂时的HU值，

是由造影剂引起的分量。由这些数据可以预则正式的、也就是用于检查的药丸注入在相同位置的增强值的变化过程c_R(t)：

c_{R} (t) = c_{0} + {\tilde{c}}_{R} (t) \cdot \cdot \cdot (4)

这在改变了的造影剂药丸b_R(t)以下式被应用或者说注射时成立：

b_R(t)＝F_RΘ(t-t_0R)Θ(t_FR-t) (5)

公知通过简单的和相应于延长的药丸注射时间而错开累加测得的测试药丸曲线可以推导出正式的药丸注入和由此产生的造影剂浓度。该方法的功效很可靠，但其缺点是将造影剂流动的持续时间限制为测试药丸的持续时间的整数倍。

同样还公知采用对所测量的测试药丸曲线的傅立叶变换来推导出正式的药丸注入效果。通过对测试药丸曲线进行傅立叶变换的变通方式，可以避开持续时间是测试药丸持续时间整数倍的限制。但由于测试药丸曲线一般由少数强噪声数据点组成，因此傅立叶变换的计算在数值上非常不稳定，并由此不适用于临床上日常记录的数据。

此外这两种方法还在测试药丸数据不包括整个增强曲线、也就是扫描得太早或太晚时出现问题。由于对于稳定的预测造影剂再循环也起着很重要的作用，因此对这两种方法都必须对测试药丸实际进行很长时间的测量。

根据本发明上述问题是通过另一种计算方法来解决的，该方法也是基于线性原因/效果关系。此外对于在测试运行期间没有测量的检查时间段采用一个生理模型进行预测。

采用简单线性法则进行的预测：

该方法的重要假定是原因(＝造影剂注入)和作用效果(在CT检查的情况下是HU值的提高)之间的线性性。该假定可以通过以下数学关系表达：

\tilde{c} (t) = {&Integral;}_{- \infty}^{\infty} {dt}^{'} k (t - t^{'}) b (t^{'}) \cdot \cdot \cdot (6)

其中b在以下表示药丸曲线，c表示造影剂曲线。k在此是特定于患者的并描述身体对注入造影剂的反应。

如果考察该公式的傅立叶变换，则得出

\tilde{C} (ξ) = K (ξ) B (ξ) \cdot \cdot \cdot (7)

其中函数f(t)的傅立叶变换在此通过下式给出

F (ξ) = {&Integral;}_{- \infty}^{\infty} dtf (t) \exp (iξt) . \cdot \cdot \cdot (8)

关系式(7)既适用于测试药丸也适用于待预测的药丸。这样利用测试药丸来确定K(ξ)，其中ξ是傅立叶空间中随时间变化的变量：

K (ξ) = \frac{{\tilde{C}}_{T} (ξ)}{B_{T} (ξ)} \cdot \cdot \cdot (9)

从而能随后用下面的公式确定浓度变化的傅立叶变换

{\tilde{C}}_{R} (ξ) = \frac{B_{R} (ξ)}{B_{T} (ξ)} {\tilde{C}}_{T} (ξ) . \cdot \cdot \cdot (10)

然后通过反向傅立叶变换

f (t) = \frac{1}{2 π} {&Integral;}_{- \infty}^{\infty} dξF (ξ) \exp (- iξt) \cdot \cdot \cdot (11)

给出所寻找的函数该过程一般用于通过傅立叶变换确定造影剂的变化过程。但由于测试药丸的傅立叶变换已知，因此可以针对一次药丸注射的任意另一个矩形来预测造影剂变化过程，还可以对完全任意的造影剂注射来预测造影剂变化过程。在此“矩形”应理解为在注射期间造影剂在特定时间段内的恒定流量，其被图形地和理想地表示为矩形。对于任意另一个矩形，也遵守上述详细的公式推导，任意的注射曲线也遵守该最终结果。

两个药丸(1)和(5)的傅立叶变换通过下式给出：

B_{I} (ξ) = F_{I} \frac{1}{iξ} (\exp (iξ t_{FI}) - \exp (i {ξt}_{0 I})) \cdot \cdot \cdot (12)

从而所寻找的增强曲线的傅立叶变换是

{\tilde{C}}_{R} (ξ) = \frac{F_{R}}{F_{T}} \frac{\exp (iξ t_{FR}) - \exp (i {ξt}_{0 R})}{\exp (iξ t_{FT}) - \exp (i {ξt}_{0 T})} {\tilde{C}}_{T} (ξ) . \cdot \cdot \cdot (13)

所寻找的增强曲线则通过下式给出

{\tilde{c}}_{R} (t) = \frac{1}{2 π} \frac{F_{R}}{F_{T}} {&Integral;}_{- \infty}^{\infty} dξexp (- iξt) \frac{\exp (i {ξt}_{FR}) - \exp (iξ t_{0 R})}{\exp (i {ξt}_{FT}) - \exp (iξ t_{0 T})} {\tilde{C}}_{T} (ξ) . \cdot \cdot \cdot (14)

该表达式可以在不知道傅立叶变换情况下积分，其中展开所出现的相位因子

\frac{\exp (iξ t_{FR}) - \exp (i {ξt}_{0 R})}{\exp (iξ t_{FT}) - \exp (iξ t_{0 T})} = \frac{\exp (i {ξt}_{0 R}) 1 - \exp (i {ξΔ}_{R})}{\exp (iξ t_{0 T}) 1 - \exp (iξ Δ_{T})}

= \frac{\exp (iξ t_{0 R})}{\exp (i {ξt}_{0 T})} (1 - \exp (iξ Δ_{R})) Σ_{n = 0}^{\infty} {(\exp (iξ Δ_{T}))}^{n}

= \frac{\exp (iξ t_{0 R})}{\exp (i {ξt}_{0 T})} (1 - \exp (i {ξΔ}_{R})) Σ_{n = 0}^{\infty} (\exp (inξ Δ_{T}) \cdot \cdot \cdot (15)

其中采用了简化式Δ_I＝t_FI-t_0I。由此给出

{\tilde{c}}_{R} (t) = \frac{1}{2 π} \frac{F_{R}}{F_{T}} {&Integral;}_{- \infty}^{\infty} dξexp (- iξ (t - t_{0 R} + t_{0 T})) (1 - \exp (iξ Δ_{R})) Σ_{n = 0}^{\infty} \exp (inξ Δ_{T}) {\tilde{C}}_{T} (ξ)

= \frac{1}{2 π} \frac{F_{R}}{F_{T}} Σ_{n = 0}^{\infty} {&Integral;}_{- \infty}^{\infty} dξ (1 - \exp (iξ Δ_{R})) \exp (- iξ (t - t_{0 R} + t_{0 T} + n Δ_{T}))) {\tilde{C}}_{T} (ξ) . \cdot \cdot \cdot (16)

= \frac{F_{R}}{F_{T}} Σ_{n = 0}^{\infty} ({\tilde{c}}_{T} (t - t_{0 R} + t_{0 T} - n Δ_{T}) - {\tilde{c}}_{T} (t - t_{FR} + t_{0 T} - n Δ_{T}))

所出现的无穷和表示对于具体的计算是没有问题的，因为所出现的测试药丸的响应函数

分别在特定时刻之前和之后被假定为消失，也就是说，只要还没有药丸流动或者在注射之后过了一段较长时间还没有药丸流动，则可以不观察增强。因此在实际情况下只需累加有限数目的和项。

此外当待预测的造影剂变化过程的持续时间不是测试药丸持续时间的整数倍时，该公式还表示对简单累加的概括，对于整数倍来说该公式减小为一个简单和，其中测试药丸曲线随时间的推移而累加。

对于一般的造影剂变化过程b_R(t)来说，类似地给出

{\tilde{c}}_{R} (t) = \frac{1}{2 π} \frac{i}{F_{T}} {&Integral;}_{- \infty}^{\infty} dξξexp (- iξt) \frac{B_{R} (ξ)}{\exp (i {ξt}_{FT}) - \exp (iξ t_{0 T})} {\tilde{C}}_{T} (ξ) . \cdot \cdot \cdot (17)

通过展开分母并利用以下认知

FT(b′_R)(ξ)＝-iξB_R(ξ) (18)

得到

{\tilde{c}}_{R} (t) = \frac{1}{2 π} \frac{1}{F_{T}} Σ_{n = 0}^{\infty} {&Integral;}_{- \infty}^{\infty} dξ \exp (- iξ (t + t_{0 T} - n Δ_{T})) F T ({b^{'}}_{R}) {(ξ) \tilde{C}}_{T} (ξ) \cdot \cdot \cdot (19)

= \frac{1}{F_{T}} Σ_{n = 0}^{\infty} {&Integral;}_{- \infty}^{\infty} {dt}^{'} {\tilde{c}}_{T} (t + t_{0 T} - n Δ_{T} - t^{'}) {b^{'}}_{R} (t^{'}) \cdot \cdot \cdot (20)

从而对于任意造影剂变化过程也只需要数据而不需要其傅立叶变换。

利用该方法基本上覆盖了事先通过基于测试药丸注入的测量值而被测量过的预测时间段。如果还要对超过测试测量的中断时刻进行预测，则发明人建议，还采用简单的生理模型并从中获得预测值。

为此采用生理系统的以下特定差分方程：

\frac{&PartialD;}{&PartialD; t} b (x, t) + v \frac{&PartialD;}{&PartialD; x} b (x, t) - D \frac{{&PartialD;}^{2}}{&PartialD; x^{2}} b (x, t) = F δ^{(1)} (x) Θ (t - t_{0}) Θ (t_{F} - t) \cdot \cdot \cdot (21)

其中方程式左端是一维的导热方程，其通过右端的一个源项得到补充。

该右侧的源项描述测试药丸，其中F表示流速，t₀和t_F表示观察的开始和结束时间。一维的Delta函数δ⁽¹⁾(x)在此描述在开始位置0的注射穿刺点。左侧基本上由两个过程确定：造影剂以漂移速度v运动，同时刚开始为矩形的测试药丸以扩散常量D扩散。在此通过扩散常量D来基本上描述穿越时间差，并模拟血流速度的波动。

该差分方程的解通过确定该差分算子的Green函数来给定：

(\frac{&PartialD;}{&PartialD; t} + v \frac{&PartialD;}{&PartialD; x} - D \frac{{&PartialD;}^{2}}{{&PartialD; x}^{2}}) G (x - x^{'}, t - t^{'}) = δ^{(1)} (x - x^{'}) δ^{(1)} (t - t^{'}) \cdot \cdot \cdot (22)

由此接着可以用下式直接给出差分方程的解

b (x, t) = {&Integral;}_{- \infty}^{\infty} {Fδ}^{(1)} (x^{'}) Θ (t^{'} - t_{0}) Θ (t_{F} - t^{'}) G (x - x^{'}, t - t^{'}) d t^{'} {dx}^{'} . \cdot \cdot \cdot (23)

一维导热方程的Green函数具有下列形式：

G (x - x^{'}, t - t^{'}) = Θ (t - t^{'}) \frac{1}{\sqrt{4 D (t - t^{'}) π}} \exp (- \frac{{(x - x^{'} - v (t - t^{'}))}^{2}}{4 D (t - t^{'})}) \cdot \cdot \cdot (24)

由此给出了差分方程的解，如下所示：

b (x, t) = {&Integral;}_{- \infty}^{\infty} F δ^{(1)} (x^{'}) Θ (t^{'} - t_{0}) Θ (t_{F} - t^{'}) G (x, t; x^{'}, t^{'}) {dt}^{'} {dx}^{'}

= F {&Integral;}_{- \infty}^{\infty} {dt}^{'} Θ (t^{'} - t_{0}) Θ (t_{F} - t^{'}) Θ (t - t^{'}) \frac{1}{\sqrt{4 D (t - t^{'}) π}} \exp (- \frac{{(x - v (t - t^{'}))}^{2}}{4 D (t - t^{'})})

= FΘ (t - t_{0}) (Θ (t_{F} - t) {&Integral;}_{t_{0}}^{t} {dt}^{'} + Θ (t - t_{F}) {&Integral;}_{t_{0}}^{t_{F}} {dt}^{'}) \frac{1}{\sqrt{4 D (t - t^{'}) π}} \exp (- \frac{{(x - v (t - t^{'}))}^{2}}{4 D (t - t^{'})})

= \frac{F}{2 v} Θ (t - t_{0}) (Θ (t_{F} - t) (1 - \exp (\frac{xv}{D}) + \exp (\frac{xv}{D}) Erf (\frac{x + v (t - t_{0})}{\sqrt{4 D (t - t_{0})}}) - Erf (\frac{x - v (t - t_{0})}{\sqrt{4 D (t - t_{0})}}))

+ Θ (t - t_{F}) (\exp (\frac{xv}{D}) Erf (\frac{x + v (t - t_{0})}{\sqrt{4 D (t - t_{0})}}) - Erf (\frac{x - v (t - t_{0})}{\sqrt{4 D (t - t_{0})}})

- \exp (\frac{xv}{D}) Erf (\frac{x + v (t - t_{F})}{\sqrt{4 D (t - t_{F})}}) + Erf (\frac{x - v (t - t_{F})}{\sqrt{4 D (t - t_{F})}}))) \cdot \cdot \cdot (25)

误差函数Erf通过下式给出：

Erf (x) = \frac{2}{\sqrt{π}} {&Integral;}_{- \infty}^{x} \exp ({- t}^{2}) dt . \cdot \cdot \cdot (26)

由于误差函数只对其自变数的零点附近的值才具有不等于±1的值，因此可如下近似对x＞0和v＞0时的解：

b (x, t) = \frac{F}{2 v} Θ (t - t_{0}) (Θ (t_{F} - t) (1 - Erf (\frac{x - v (t - t_{0})}{\sqrt{4 D (t - t_{0})}}))

+ Θ (t - t_{F}) (Erf (\frac{x - v (t - t_{F})}{\sqrt{4 D (t - t_{F})}})) - Erf (\frac{x - v (t - t_{0})}{\sqrt{4 D (t - t_{0})}}))

因此如果通过下式描述测试药丸曲线

b (x, t) = c_{0} + C Θ (t - t_{0 T}) (Θ (t_{FT} - t) (1 - Erf (A \frac{B - (t - t_{0 T})}{\sqrt{(t - t_{0 T})}}))

+ Θ (t - t_{FT}) (Erf (A \frac{B - (t - t_{FT})}{\sqrt{(t - t_{FT})}})) - Erf (A \frac{B - (t - t_{0 T})}{\sqrt{(t - t_{0 T})}})) \cdot \cdot \cdot (27)

则通过下式给出所寻找的造影剂曲线：

b (x, t) = c_{0} + C \frac{F_{R}}{F_{T}} Θ (t - t_{0 R}) (Θ (t_{FR} - t) (1 - Erf (A \frac{B - (t - t_{0 R})}{\sqrt{(t - t_{0 R})}}))

+ Θ (t - t_{FR}) (Erf (A \frac{B - (t - t_{FR})}{\sqrt{(t - t_{FR})}})) - Erf (A \frac{B - (t - t_{0 R})}{\sqrt{(t - t_{0 R})}})) . \cdot \cdot \cdot (28)

因此，通过对测试药丸曲线的拟和确定差分方程的自由参数，从而能随后用这些参数进行预测。

或者，在本发明的范围中还可以用Gamma变量分布来匹配该生理函数，但由此也可以非完全地描述测试药丸曲线。

根据上述展示的基本思想，发明人建议对于本身公知的用于预测生命体内、尤其是患者体内造影剂流动的方法，其中将包含具有公知注射流动变化过程的造影剂的特定测试药丸注入体内，优选注入血管、特别是静脉中；利用断层造影方法在包含多个测量时刻的有限时间段z期间观察和确定造影剂在体内至少一个地点的时间浓度变化过程；并从获得的关于造影剂分布的测量数据中借助线性原因/效果依据预测另一次造影剂注入的造影剂浓度的时间变化过程，这样进行改进，其中为了预测造影剂在体内至少一个事先测量过的地点x的时间上的浓度变化过程

而采用以下计算公式：

{\tilde{c}}_{R} (t) = \frac{1}{F_{T}} Σ_{n = - \infty}^{\infty} {&Integral;}_{- \infty}^{\infty} d t^{'} {\tilde{c}}_{R} (t + t_{0 T} - n Δ_{T} - t^{'}) {b^{'}}_{R} (t^{'})

其中，

为对应于从时刻t推移到时刻t+t_0T-nΔ_T-t′的浓度，

F_T为测试药丸的造影剂的流速，

b′_R(t′)为所注入的造影剂药丸的变化过程的时间导数，

t为预测时刻

t’对应于积分变量。

在此要指出，所给定的从-∞到+∞的界限当然是理论上的，实际中替换成测量前后相应距离较远的时刻。

在采用在注射时间期间具有恒定流量的造影剂注射来作为测试药丸和要计算的正式药丸的情况下，建议用以下公式计算浓度变化过程

的预测：

{\tilde{c}}_{R} (t) = \frac{1}{2 π} \frac{F_{R}}{F_{T}} {&Integral;}_{- \infty}^{\infty} dξexp (- iξ (t - t_{0 R} + t_{oT})) (1 - \exp (iξ Δ_{R})) Σ_{n = 0}^{\infty} \exp (inξ Δ_{T}) {\tilde{C}}_{T} (ξ)

= \frac{1}{2 π} \frac{F_{R}}{F_{T}} Σ_{n = 0}^{\infty} {&Integral;}_{- \infty}^{\infty} dξ (1 - \exp (iξ Δ_{R})) \exp (- iξ (t - t_{0 R} + t_{0 T} + n Δ_{T}))) {\tilde{C}}_{T} (ξ)

= \frac{F_{R}}{F_{T}} Σ_{n = 0}^{\infty} ({\tilde{c}}_{T} (t - t_{0 R} + t_{0 T} - n Δ_{T}) - {\tilde{c}}_{T} (t - t_{FR} + t_{0 T} - n Δ_{T}))

其中采用以下标记：

F_T为测试药丸的造影剂的流速

F_R为正式药丸的造影剂的流速

ξ为积分变量

t_0R为正式药丸的起始时刻

t_0T为测试药丸的起始时刻。

如果另外在不存在来自测试药丸注射的测量值的时刻根据生理计算模型预测造影剂浓度，则本发明是很有利的。

在此作为生理计算模型，优选采用以下差分方程：

\frac{&PartialD;}{&PartialD; t} b (x, t) + v \frac{&PartialD;}{&PartialD; x} b (x, t) - D \frac{{&PartialD;}^{2}}{&PartialD; x^{2}} b (x, t) = F δ^{(1)} (x) Θ (t - t_{0}) Θ (t_{F} - t)

其中采用以下标记：

b(x，t)为药丸在时刻t位置x的浓度变化

F为造影剂的流量

δ⁽¹⁾为Delta函数，

Θ(t-t₀)为用于描述药丸注射的开始的海维赛德阶梯函数

Θ(t_F-t)为用于描述药丸注射的结束的海维赛德阶梯函数。

根据该生理模型，为了预测没有来自测试药丸注射的测量值的时刻，首先用以下函数近似测试药丸注射之后的测量的浓度变化过程：

b (x, t) = c_{0} + C Θ (t - t_{0 T}) (Θ (t_{FT} - t) (1 - Erf (A \frac{B - (t - t_{0 T})}{\sqrt{(t - t_{0 T})}}))

+ Θ (t - t_{FT}) (Erf (A \frac{B - (t - t_{FT})}{\sqrt{(t - t_{FT})}})) - Erf (A \frac{B - (t - t_{0 T})}{\sqrt{(t - t_{0 T})}}))

接着用这样计算出的参数A、B、C和c₀根据以下公式估计期待的浓度变化过程：

b (x, t) = c_{0} + C \frac{F_{R}}{F_{T}} Θ (t - t_{0 R}) (Θ (t_{FR} - t) (1 - Erf (A \frac{B - (t - t_{0 R})}{\sqrt{(t - t_{0 R})}}))

+ Θ (t - t_{FR}) (Erf (A \frac{B - (t - t_{FR})}{\sqrt{(t - t_{FR})}})) - Erf (A \frac{B - (t - t_{0 R})}{\sqrt{(t - t_{0 R})}}))

其中采用以下标记：

A为第一函数常量，基本上与测试药丸曲线的宽度间接成正比，

B为第二函数常量，基本上与测试药丸曲线的峰值成正比，

b(x，t)为药丸在时刻t位置x的浓度变化

C为第三函数常量，与测试药丸曲线下包含的面积成正比，

c₀为在注射造影剂药丸之前的增强值，

Erf()为误差函数，

F_R为正式药丸的造影剂的流速，

F_T为测试药丸的造影剂的流速，

t_FR为正式药丸的结束时刻，

t_FT为测试药丸的结束时刻，

t_0R为正式药丸的起始时刻，

t_0T为测试药丸的起始时刻，

x为被观察地点

Θ为用于描述药丸注射的开始和结束的海维赛德阶梯函数。

如果对于预测时刻既根据线性模型预测又根据生理模型预测，则可以采用两个预测值的简单平均值或两个预测值的加权平均值。

附图说明

下面借助附图中结合优选实施例详细解释本发明，其中采用以下附图标记：1：计算机断层造影系统；2：X射线管；3：检测器；4：系统轴；5：外壳；6：可移动患者卧榻；7：患者；8：自动注射器的控制导线；9：计算单元；10：数据和控制导线；11：受控的注射单元；12：静脉入口；13：对造影剂浓度的预定目标输入；14-17：增强曲线；18：预测曲线；18+：图18的预测范围的统计上限；18-：曲线18的预测范围的统计下限；19，20：观察测试药丸注入的时限；A：第一函数常量，基本上间接与测试药丸曲线的宽度成正比；B：第二函数常量，基本上与测试药丸曲线的峰值成正比；b(x，t)：药丸在时刻t在地点x的浓度变化；C：第三函数常量，与测试药丸曲线下包含的面积成正比；c₀：在注射造影剂药丸之前的增强值；Erf()：误差函数；F_R：正式药丸的造影剂流速；F_T：测试药丸的造影剂流速；t_FR：FR的预测时刻；t_FT：FT的测试时刻；t_0R：FR预测的起始时刻；t_0T：FT的起始时刻；t₁-p_N：程序和程序模块；x：被观察地点。

具体示出：

图1：示出CT系统的示意图；

图2：示出测试药丸注射的流程；

图3：示出注射测试药丸之后的浓度变化过程；

图4：示出正式的造影剂注射的期望的、预测的浓度变化过程；

图5：示出图4的正式造影剂注射的流程；

图6：示出测试药丸注射之后的增强数据的理论变化过程；

图7：示出正式的造影剂注射之后增强数据的预测变化过程；

图8：示出具有附加显示的造影剂注射之后测量值的图7。

具体实施方式

图1示出优选采用本发明方法的计算机断层造影系统1，其以不同的形式采用。在所示例子中该计算机断层造影系统1具有X射线管2和与其相对设置的检测器3，它们可旋转地设置在一个支架上。在X射线管2和检测器3旋转期间，患者7在患者卧榻6上沿着系统轴4移动地经过X射线管2和检测器3，从而相对于患者进行螺旋形扫描。X射线管2和检测器3与支架一起位于外壳5中，外壳5通过数据和控制导线10与计算单元9连接。为了注射造影剂由计算单元9通过控制导线8控制注射器11，该注射器11通过静脉入口12以期望的流速并在期望的时刻向患者注入造影剂。

如果注射所谓的测试药丸，则注射器上的造影剂流速就如图2所示。图2示出时间轴t上的流速F_T。测试药丸的注射在时刻t_0T开始，并具有以阴影示出的矩形变化过程。

由于这种测试药丸注射和测试药丸注射期间进行的测试扫描，可以确定在测试药丸注射之后发生的体内造影剂分布。

图3示出增强曲线，也就是造影剂在患者的被观察位置(如心脏的心房)的图像响应。在纵坐标上绘出浓度值c_T(x)，其与从图像显示中测量的增强值相关。所显示的造影剂在被检查地点随时间变化的浓度曲线展示出典型的陡峭上升和大致对应于测试药丸注射的持续时间的短暂的平稳，以及接下来的陡峭下降和后面缓慢停下的低平稳。

这种造影剂注入的目的是，在利用计算机断层造影设备进行断层造影检查期间获得被观察区域内的足够造影剂浓度，以保证例如对心脏动脉的很好的显示。因此需要造影剂具有允许相应显示的特定浓度。但同时浓度不应当过高，因为应当将造影剂的负生物作用保持得尽可能地小。

图4用阴影矩形示出在预定时间间隔期间造影剂浓度的这种期望范围13，该时间间隔对应于检查时间间隔。由于测试药丸数据已知，现在应当进行正式的造影剂注射，其最后会导致造影剂浓度在扫描期间的足够的、但不是太大的变化。

图5示例性示出流速FR的这种变化过程，这可以通过本发明的方法找出。

图6示出一系列测量点，表示为关于时间t的HU值，如对于测试药丸(对应于图2)来说在心房区域内找出的。时间t₀在此对应于测试药丸注射的开始t_0T。除了基于测试药丸注射确定的HU值之外还显示出点式的三条曲线，其中第一陡峭曲线14对应于造影剂从注射位置到测量位置的直接传送。下一条较小的曲线15示出惯性的较慢的造影剂，最后曲线16示出由于造影剂在血液循环中的再循环而造成的影响。曲线14至16的累加对应于体内所有三种效果的和，并对应于在观察点x的实际发现的浓度变化过程17。

如果基于图6中显示的测试药丸的测量值和先前示出的计算方法预先计算预测的浓度值，或者说在此是按照HU单位的增强值，则给出如图7所示的曲线18。在曲线18上面的曲线18+和下面的曲线18-分别限定了期待的统计置信区间。在优选的将线性依据与生理依据组合的计算方法中，可以超过测试药丸的测量时间边界(该边界在此通过垂直的虚线19和20来表示)来预测期待的浓度值，从而在图7中通过生理模型近似越过时间边界20的预测值。

图8中示出以空正方形表示的实际确定的值和理论预测曲线18之间的比较。

可以理解，本发明的上述特征不仅能以给定的组合，还能在不偏离本发明范围的情况下以其它组合或单独使用。

此外还要指出，本发明的方法不仅可用于计算机断层造影系统，还能用于NMR断层造影系统或与C型X射线系统一起使用。

总之，通过本发明的方法改善了对造影剂注射之后增强值的预测，并由此在推断时表现出一种改善了的可能性，即如何基于给定的或期望的患者体内预定地点的造影剂浓度来预先计算对应的造影剂浓度，尤其是还有其时间变化过程。

Claims

1.一种预测生命体内、尤其是在患者体内造影剂流动的方法，其中：

1.1.将包含具有公知注射流动变化过程的造影剂的特定测试药丸注入体内，优选注入血管、特别是静脉中，

1.2.利用断层造影方法在包含多个测量时刻的有限时间段Z期间观察和确定造影剂在体内至少一个地点的时间浓度变化过程，

1.3.从获得的关于造影剂分布的测量数据中，利用线性原因/效果依据来预测另一次造影剂注入的造影剂浓度的时间变化过程，

其特征在于，

1.4.为了预测造影剂在体内至少一个事先测量过的地点的时间浓度变化过程

采用以下计算公式：

{\tilde{c}}_{R} (t) = \frac{1}{F_{T}} Σ_{n = - \infty}^{\infty} {&Integral;}_{- \infty}^{\infty} d t^{'} {\tilde{c}}_{R} (t + t_{0 T} - n Δ_{T} - t^{'}) {b^{'}}_{R} (t^{'})

其中，

为对应于从时刻t推移到时刻t+t_0T-nΔ_T-t′的浓度，

F_T为测试药丸的造影剂的流速，

b′_R(t′)为所注入的造影剂药丸的变化过程的时间导数，

t为预测时刻

t’对应于积分变量。

2.根据权利要求1所述的方法，其特征在于，采用在注射时间期间具有恒定流动的造影剂注射来作为测试药丸和最终要计算的药丸，并用以下公式计算浓度变化过程

的预测：

{\tilde{c}}_{R} (t) = \frac{1}{2 π} \frac{F_{R}}{F_{T}} {&Integral;}_{- \infty}^{\infty} dξexp (- iξ (t - t_{0 R} + t_{0 T})) (1 - \exp (iξ Δ_{R})) Σ_{n = 0}^{\infty} \exp (inξ Δ_{T}) {\tilde{C}}_{T} (ξ)

= \frac{1}{2 π} \frac{F_{R}}{F_{T}} Σ_{n = 0}^{\infty} {&Integral;}_{- \infty}^{\infty} dξ (1 - \exp (iξ Δ_{R})) \exp (- iξ (t - t_{0 R} + t_{0 T} + n Δ_{T})) {\tilde{C}}_{T} (ξ)

= \frac{F_{R}}{F_{T}} Σ_{n = 0}^{\infty} ({\tilde{c}}_{T} (t - t_{0 R} + t_{0 T} - n Δ_{T}) - {\tilde{c}}_{T} (t - t_{FR} + t_{0 T} - n Δ_{T}))

其中采用以下标记：

F_T为测试药丸的造影剂流速，

F_R为正式药丸的造影剂流速，

ξ为积分变量，

t_0R为正式药丸的起始时刻，

t_0T为测试药丸的起始时刻。

3.根据权利要求1或2所述的方法，其特征在于，为了预测不存在来自测试药丸注射的测量值的时刻，根据生理计算模型预测所述造影剂浓度。

4.根据权利要求3所述的方法，其特征在于，作为生理计算模型采用以下差分方程：

\frac{&PartialD;}{&PartialD; t} b (x, t) + v \frac{&PartialD;}{&PartialD; x} b (x, t) - D \frac{{&PartialD;}^{2}}{{&PartialD; x}^{2}} b (x, t) = {Fδ}^{(1)} (x) Θ (t - t_{0}) Θ (t_{F} - t)

其中采用以下标记：

b(x，t)为药丸在时刻t在地点x的浓度变化，

F为造影剂的流动，

δ⁽¹⁾为Delta函数，

Θ(t-t₀)为用于描述药丸注射的开始的海维赛德阶梯函数，

Θ(t_F-t)为用于描述药丸注射的结束的海维赛德阶梯函数。

5.根据权利要求3或4所述的方法，其特征在于，为了预测不存在来自测试药丸注射的测量值的时刻，首先用以下函数近似测试药丸注射之后的测量的浓度变化过程：

b (x, t) = c_{0} + CΘ (t - t_{0 T}) (Θ (t_{FT} - t) (1 - Erf (A \frac{B - (t - t_{0 T})}{\sqrt{(t - t_{0 T})}}))

+ Θ (t - t_{FT}) (Erf (A \frac{B - (t - t_{FT})}{\sqrt{(t - t_{FT})}})) - Erf (A \frac{B - (t - t_{0 T})}{\sqrt{(t - t_{0 T})}}))

接着用这样计算出的参数A，B，C，c₀根据以下公式估计期待的浓度变化过程：

b (x, t) = c_{0} + C \frac{F_{R}}{F_{T}} Θ (t - t_{0 R}) (Θ (t_{FR} - t) (1 - Erf (A \frac{B - (t - t_{0 R})}{\sqrt{(t - t_{0 R})}}))

+ Θ (t - t_{FR}) (Erf (A \frac{B - (t - t_{FR})}{\sqrt{(t - t_{FR})}})) - Erf (A \frac{B - (t - t_{0 R})}{\sqrt{(t - t_{0 R})}}))

其中采用以下标记：

B为第二函数常量，基本上与测试药丸曲线的峰值成正比，

b(x，t)为药丸在时刻t在地点x的浓度变化，

C为第三函数常量，与测试药丸曲线下包含的面积成正比，

c₀为在注射造影剂药丸之前的增强值，

Erf()为误差函数，

F_R为正式药丸的造影剂的流速，

F_T为测试药丸的造影剂的流速，

t_FR为正式药丸的结束时刻，

t_FT为测试药丸的结束时刻，

t_0R为正式药丸的起始时刻，

t_0T为测试药丸的起始时刻，

x为被观察地点，

Θ为用于描述药丸注射的开始和结束的海维赛德阶梯函数。

6.根据权利要求3至5中任一项所述的方法，其特征在于，对于根据线性模型和生理模型进行预测的时刻，采用两个预测值的平均值。

7.根据权利要求3至5中任一项所述的方法，其特征在于，对于根据线性模型和生理模型进行预测的时刻，采用两个预测值的加权平均值。