具体实施方式
现在将参考附图描述本发明的具体实施方式。然而,本发明可被实现为各种不同形式,而不应当被解释为限于本文所阐明的实施方式;相反,这些实施方式被提供以便本公开完全并完整,并且将把本发明的范围完全地传递给本领域的技术人员。更详细地,下面的描述中给出的一些实施方式聚焦于适用于肝功能性评估并尤其适用于用于肝功能性评估的基于MRI的图像分析的方法和系统。然而,应当理解,本发明不限于本应用,而是可被应用于其它医学领域、程序和/或分泌或排泄器官,例如包括下面进一步提到的那些。
在一个实施方式中,T1-加权动态肝的特异性对比度增强(DHCE)MRI提供了3D图像数据。在随后的时间所获得的多个3D数据组提供了4D图像数据组。处理该4D图像数据组以在肝的节段或子节段水平上进行肝功能的评估。该节段可以小至该4D图像数据的体素(voxel)水平。根据处理该4D数据组确定肝功能的参数值。
在节段或子节段水平上,相对于肝的输入血流(输入相对血流,irBF)确定肝中的血流。
在节段或子节段水平上相对于肝中的静脉血流确定肝中的血流。
可以在整个肝上确定该肝中的血流。替代或附加地,可以在节段或子节段水平上确定该血流。可以相对于该肝中的动脉血流确定该血流。
在一些实施方式中,肝提取分数(HEF)是在节段或子节段水平(直至体素水平)上确定的。HEF先前只是在整个器官的水平上被确定的。当在节段或子节段水平(直至体素水平)上确定HEF时,新的并且更有效的诊断和治疗机会出现。
肝的功能可以根据其的容积单元来确定。肝或其节段的容积可以从由图像模态(modality)(例如,MRI模态)所提供的3D图像数据来确定。用这种方式可以将局部HEF与肝的具体容积相关,即,HEF/容积被确定。
例如,这允许进行虚拟的手术过程规划。在切除肝的病变部位之前,进行对器官的一部分的手术切除的基于计算机的虚拟规划可允许计算手术之后的残留肝功能。这从临床医生的视角和病人的安全方面来看都是很大的优点。
在一些实施方式中,肝的肝提取分数(HEF)和/或输入相对血流(irBF)是在节段或子节段器官水平上,根据截断奇异值分解(TSVD)计算来确定的。例如这在计算方面是有优势的,其允许临床上可接受的计算时间。
在一些实施方式中,TSVD被有益地用来确定参数映射。该参数映射提供了对器官功能的有效而快速地诊断。
在一些实施方式中,HEF和irBF被以参数映射的形式可视化。
在一些实施方式中,HEF和irBF被以叠加在解剖图像上的参数映射的形式可视化。
计算结果、对手术过程或其它治疗的虚拟规划可以例如被呈现在医学工作站的显示器上。可以使得基于实施方式的计算结果的过程或治疗的规划在例如下文参考图12所描述的系统的医学工作站的显示器上,以由用户输入所操纵的互动的方式看到。
在一个实施方式中,用磁共振成像(MRI)并将用钆塞酸二钠(Gadoliniumethoxybenzyl diethylenetriaminepentaacetic acid,Gd-EOB-DTPA(德国柏林的Schering AG公司的))用作针对肝的造影剂的动态肝功能测试被使用。Gd-EOB-DTPA具有通过肝和肾通道,等同消除的独特性质。
对于Gd-EOB-DTPA,造影剂的吸收和排泄的双路径(50%肝细胞吸收、胆汁排泄)(50%通过肾小球过滤而肾消除),允许本模型用于同时监控和/或确定肝和肾功能两者,这是本方法与这种造影剂相结合的独特性质。因此,动态肝肾的特异性对比度增强(DHRCE)磁共振成像(MRI)被提供并被用于在各实施方式中或为各实施方式获取3D或4D图像。
能用于临床使用并适于本方法的一些方面的其它基于钆的造影剂包含:Bayer Schering Pharma公司的Gadopentate dimeglumine;GE Healthcare公司的Gadodiamide;Gothia/Guerbet公司的Gd-DOTA;Initios Medical AB/Bracco公司的Gadoteridol;BayerSchering Pharma公司的Gadobutrol。
能用于临床使用并适于本方法的其它方面的其它组织特异性造影剂包含:Gothia/Guerbet公司的(SPIO)(80-150nm)、Bayer ScheringPharma公司的(SPIO)(60nm);GE Healthcare公司的Mangafodipir;Initios Medical AB/Bracco公司的Gadobenatedimeglumine、Bayer Schering Pharma公司的Gadofosvesettrinatrium。然而,后面的造影剂可能不适于肝细胞特异性对比度增强,但是适于其它分泌或排泄器官(比如肾)的对比度增强。
本方法和系统的实施方式不限于使用Gd-EOB-DTPA作为造影剂。其它未来或当前可获得的肝或器官特性性造影剂也是可适用的。
当造影剂被以静脉注射的方式给予时,肝中的浓度受血液循环的影响并且随着时间分散。因此,响应函数可被描述为脉冲响应和输入函数之间的卷积。图2A显示了理想情况,此时所关注的器官由短脉冲函数在施加脉冲响应的情况下呈现。如同提到过的,事实上感兴趣的器官是由输入函数呈现的,该输入函数随时间变化并将因此影响该响应函数,如图2B中所示。
为了克服示踪剂(tracer)循环的影响,使用传入血管的相对增强曲线作为输入函数并使用肝相对增强曲线作为响应函数来应用解卷积分析(DA)。为了DA,执行矩阵求逆以及奇异值分解(SVD)。
测量到的肝造影剂增强可能更依赖于肝灌注而不是实际的肝细胞功能,为输入函数规则采取措施(making provision for input functionimperative)。就具有高肝提取率的示踪剂而言,这更加真实。
在理想情况下,为了克服循环的影响,示踪剂的给予应当作为短血管内丸剂被直接提供到肝的传入血液供应(即,门静脉或肝动脉)中。在第一遍通过(first pass)过程中,周边静脉给药(如临床实践中使用的)将在肝内呈现出只有很小百分比的注入示踪剂,等于由该肝接所接收的心输出量小部分(fraction)。然后,肝将恒定地呈现出示踪剂由于循环和同时进行的提取和排泄引起的变化的浓度。当执行肝成像时,将示踪剂直接给予到门静脉或肝动脉内不是在临床情形下执行的。
然而,可以使用解卷积分析(DA)模拟该原理。DA针对正出现在该器官的造影剂的不断变化的浓度,修正器官的时间-活性曲线。该方法已被证实在使用基于傅里叶变换(FT)的解卷积的动物研究中有效。FT是闪烁扫描(scintigraphic)实践中最广泛使用的解卷积模型。
解卷积分析至今还没有用于确定器官的每容积单元的功能。对于那些不是仅仅由血液所灌注,而是另外具有分泌或排泄功能的分泌或排泄器官的情况,尤其如此。
使用奇异值分解(SVD)的矩阵求逆是用于DA的更有利的数学模型,但是至今仍没有用于从图像数据确定肝功能。
对肝细胞特异性造影剂的使用,比如钆塞酸二钠(柏林的BayerSchering Pharma AG公司的(Gd-EOB-DTPA)),已经被证明当其被用在T1-权重MRI中时,改善了所关注的肝损害的检测和特征化。Gd-EOB-DTPA的药效特性类似于99mTc-IDA-族的特性,其中肝细胞通过有机阴离子传送系统(OATS)吸收并在随后经由谷胱甘肽-S-转移酶进行胆汁排泄。药物动力研究显示,Gd-EOB-DTPA的给予剂量的约50%被肝提取并被透过肝-胆通道分泌。剩余的50%被肾排泄消除。因此,Gd-EOB-DTPA的肝吸收和后续的T1松弛缩短(T1-relaxation shortening)取决于肝细胞块的完整性。动态Gd-EOB-DTPA MRI以前已用在动物模型中,以便在各种试验性设定中或使用简要的(summary)参数或使用DA估计肝功能和机能障碍。
Gd-EOB-DTPA的药效特性结合MRI中获得的高分辨率展示了将DHCE-MRI与Gd-EOB-DTPA一起用作为基于成像的肝功能试验的优点,提供在区域和/或甚至节段水平上的功能方面的可辨别的差异。这是以前人社会没有提供的。
图1是示意性说明图像1的示意图,图像1可视化了经由显示通过腹部100的切片(slice)的MRI模态所获得的三维(3D)数据。示出了肝110,其包括软组织(parenchyma)112(肝的功能部位)和包括胆管、门静脉支流和肝动脉支流111的入口进入蒂(portalpeddicle)。还示意了下腔静脉(IVC,inferior vena cava)130(在其它节段,即使从肝泄漏到IVC的肝静脉也可以被可视化)和主动脉120。
图2A是示意性说明了与其脉冲响应卷积的脉冲函数的示意图,而图2B是示意性说明了与脉冲响应卷积的非理想输入函数的示意图;如同上文所述并且在下文详细解释的那样。
图3是示意性说明了解卷积的肝提取(HE)曲线以及肝滞留曲线(HRC)的图解。图3中显示了解卷积的肝提取(HE)曲线(即,脉冲响应)和肝滞留曲线(HRC,其是对在420和1800秒之间的HE曲线中的时间点的单指数拟合)。该HE曲线的峰值和该HRC的Y轴截断之间的比例被定义为肝提取分数(HEF)。在这种情况下,傅里叶分析(FA+尾)被用于DA,而图3中显示的相关示例中的HEF是约17%。图3将在下文被详细阐明。
图4是示意性说明获得肝提取曲线的示意图。
门静脉中的输入函数的相对随时间增强的曲线、以及来自一个试验对象的肝节段V的实质(parenchymal)响应函数被示出。符号(●和x)指示取样点。该实质响应曲线被示出,其中三个ROI的平均的95%的置信区间被放置在肝节段V中。已经以7个点的滑动窗口函数对两个曲线进行了平滑。
图5是示意了包含实施方式的方法2的流程图。在步骤200中,病人被放在磁共振成像器中。在步骤210中,使用T-1加权序列扫描该病人的肝,从而提供包含对该肝及相连的结构和器官的解剖数据的3D病人数据。
然后,在步骤220中,肝特异性造影剂被注射到病人的血流中。如在步骤230中所示,利用磁共振成像器在连续的时间(约10到90分钟期间)扫描该病人的肝。在每次扫描过程中获得新的3D数据组,这提供了四维(4D)数据组,即,提供了用于3D容积(volume)中的时域变化的数据。该4D数据组也被称为动态4D像容积。图6中描绘了这一点。
在所示的步骤240中,从该动态4D像容积中对于肝血液输入和肝实质(parenchyma)提取数据,例如使用下面更详细地描述的方法。该方法可以是由计算机实现的。
在步骤250中,计算将该血液输入传递到该肝实质响应的脉冲响应函数。例如这是通过合适的计算机程序实现的。该计算可如图7中所示的那样完成,从而提供了脉冲响应函数,其也被称为矩阵形式的肝提取曲线。
在步骤260中,肝提取分数和输入相对血流被从计算出的脉冲响应函数中局部地提取出来,以提供了用于进一步处理或分析的数据。
在步骤270中,来自步骤260的数据被用于提供肝提取分数和输入相对血流图像映射和/或节段水平上、或子节段水平上、直到体素水平的表格式结果
例如在图10A到10D中,显示了基于通过DHCE-MRI和根据不同的计算方法的后续的图像处理获得的数据的图像。图10A和10C显示了用TSVD DA分别计算的HEF和irBF(110内部)的参数映射。图10B和10C显示了用FA DA分别计算的HEF和irBF的参数映射。该参数映射是根据颜色代码表300用颜色编码的。背景(这里是在图1中的腹部100的内部)中的解剖情况是黑白显示的,以便清晰地区分以器官的节段(例如,在特定体素处)的形式所示的器官功能(HEF和irBF)。
图8是分段的肝功能性评估的示意图。
肝被分成八个节段(如I到VIII-Segment1到Segmetn8一SI到SVII所示的),所有节段都作为具有它自己的静脉血供应和胆汁(billary)排泄路径的单独的器官发挥功能。从而可以对整个肝内的每个体素(x,y,z)计算HEF。可以使用基于计算机的分段和/或目标识别(例如,基于图像强度或Hounsfield灰度值)来获得该肝容积。可以使用基于肝解剖学界标(landmark)的半自动计算机软件将肝体进一步分成解剖学肝节段。可以通过将HEF与其相应的容积相乘获得在节段或子节段水平上对整个肝的虚拟功能进行测量。
VSegment 1xHEFSegment 1=FSegment 1
VSegment 2xHEFSegment 2=FSegment 2
VSegment 3xHEFSegment 3=FSegment 3
VSegment 4xHEFSegment 4=FSegment 4
VSegment 5xHEFSegment 5=FSegment 5
VSegment 6xHEFSegment 6=FSegment 6
VSegment 7xHEFSegment 7=FSegment 7
VSegment 8xHEFSegment 8=FSegment 8
总函数=求和,FTotal
如果例如通过手术或药物治疗改变该肝功能和/或该肝容积,新的机能的测量可以使用这种技术来获得。此变化成为分数(fraction)。
Fratio=Fpre-surgery/Fpost-surgery
比值Fratio是在治疗之前和之后的函数的比并且可被应用于药物治疗和手术两者。
附加地或替代地,提供对治疗的虚拟规划。例如,可以虚拟规划出对器官的至少一个节段的至少一部分的手术切除。切除后的器官功能可以通过剩余节段的总函数确定。因此向外科医生提供了在手术切除其一部分之后所评估的器官功能是否仍将足够的有价值的信息。依据虚拟手术规划的现实手术因此与其结果适应。
所提供的对血流的测量可被用来进行虚拟规划。可以确定进入不同的肝部位的每个容积节段的输入相对血流。例如这在临床上有关于严重血管化的肿瘤。检查相关4D容积中的子容积或整个容积中的血流是我们感兴趣的。虚拟规划包括对血流的均匀考量,例如,当规划用坏死病诱导药学物质(比如)进行的治疗时。在虚拟规划过程中,可以对于血管化的肿瘤的一个区域限定这种制剂的效果。因此该虚拟规划可提供在药品治疗之后对总的肝功能的测量。
在另一个示例中,用化疗治疗可以在药品治疗的虚拟规划的计算机实现方法中被虚拟规划。因此,在实际治疗之前,提供了对中断或改变该治疗的测量,这对于病人是有利的,也有利于节省成本。
下面给出评估使用动态Gd-EOB-DTPA增强的MRI计算HEF作为节段水平上肝细胞功能的标记的可行性的一个示例。将基于傅里叶的计算方法与截断SVD(TSVD)对比,以进行解卷积分析。
进一步,下面给出使用动态Gd-EOB-DTPA增强的MRI评估具有原发性胆汁性肝硬变(PBC)和原发性硬化性胆管炎(PSC)的病人的HEF、irBF、药动学传递常数和半定量(semi-quantitative)动态参数的示例。
解卷积分析
在数学上,在肝的实施方式中,器官的响应函数可以被描述为脉冲响应和输入函数之间的卷积
其中y(t)是响应函数,h(t)是脉冲函数而x(t)是输入函数。真实肝功能通过脉冲函数来表征。图2A显示出响应y(t)等于脉冲函数x(t),如果输入函数理想的话。我们的输入函数由注射的示踪剂组成,该示踪剂由于血液循环而随着时间逐渐散布开来。因此,我们的输入函数不是理想的,并且将极大地影响响应函数y(t),如图2B中所示。响应函数y(t)和输入函数x(t)可以被测量,但是h(t)是未知的。然而,有了输入和响应函数的信息,可以估计脉冲函数,通过傅里叶分析(FA)或矩阵求逆。FA,被描述为
其中FT是傅里叶变换而FT-1是逆傅里叶变换,有简单的优点,但是会经受x(t)和y(t)的意外终止点导致的高频伪象。为了避免数据的意外终止,可以向x(t)和y(t)的终点(end)添加平滑的附加曲线以使得这些曲线回到零。这通常是附加具有为x(t)和y(t)的最后点的初始高度的从0到π/2的余弦函数完成的。应该注意,必须对每个体素执行两个傅里叶变换。这在计算上是极为需要的,特别是当具有大规模数据组的时候,比如当前病人的4D数据组。历史上,成像模态的分辨率随着其新的发展而增加,从而体素的尺寸被最小化而体素的数量增加。这种趋势进一步带来了未来更多的计算负担,使得在具有足够精确的结果的临床上可接受的计算时间在未来更不可行。
然而,通过将等式1中的卷积格式化为矩阵形式,根据一个实施方式,该等式反之可以通过矩阵求逆来求解,使用如下所示的以及如图4中所示意说明的SVD:
[3.1].
因为A是方阵,它将按照下式被分为SVD:
A=U·W·VT=U·[diag(wi)]·VT [3.2],
其中U和V是正交的(即,它们的逆等于它们的转置)而W是斜交的,其中元素wi比如
w1≥w2≥...wN≥0 [3.3]
通过矩阵求逆求解h(t):
这比上面提到的每个体素所需的两个傅里叶变换具有少的多的计算上的负担。只计算一次整体矩阵,然后对所有体素可用。
如果wi中的一个或更多个是零或接近于零,那么该矩阵求逆变得状况很坏(ill-conditioned)。因此数据中的噪声在最小二乘解(即,等式3.4)被放大,并使得结果没有实际价值。对这个问题的一个解决方案是规则化的原理,或更具体地说是截断SVD(TSVD)。
在TSVD中,阈值(范围为从0到1的c)被限定为n(1-c),其中n是奇异值的总数而c是阈值。对于小于这个切割(cut-off)的奇异值,不计算1/wi,而是用0取代。
有关解卷积的重要的注意事项是计算效率和所需的数据的量(即,4D时间分辨的(resolved)图像数据组的长度)。当将FA与基于TSVD的解卷积对比时,它们对于单个体素或ROI计算差不多快。但是对于多个体素解卷积(即,参数映射的计算),SVD在效率方面更优良,因为图7中所示的矩阵求逆只需要被计算一次然后被应用到所关注的全部体素。相反,必须使用FA对所关注的全部体素执行完全的DA。
而且,获得的数据量被该4D成像协议的长度严格限制。使用从来自20个健康志愿者的平均输入和响应函数所构造的理想输入和响应函数对协议长度的模拟表明,基于SVD的解卷积计算了同样的HEF值,即使该协议长度被缩短。像25分钟那样短的扫描协议被用于成功地计算HEF。图11中描绘了这种模拟的结果。如图11所示,当该协议变得更短时,FA DA过高估计了HEF。
图9是示意性说明对这种不同的模拟计算方法的平均误差和误差直方图的图示。
解卷积模拟
我们执行数字模拟,将具有附加尾的FA(FA+尾)和TSVD进行对比度。理想的输入和脉冲函数是从两个伽玛变量函数构造的。曲线形状被构造为尽可能类似于在活体内测量的那些。然后将这两个曲线卷积以找出响应函数,如等式1中所示。不同数量的正态分布噪声被分别应用于该响应和输入函数,以模拟不同的SNR水平。然后使用两个不同的技术应用DA。该FA+尾技术中的所附加的尾被设定为是该模拟数据的长度的三倍。该TSVD技术中的截断阈值被固定在0,07。对每个SNR水平执行1000次模拟并将使用FA+尾的结果的标准偏差与使用方差比(variance ratio)测试的TSVD获得的结果相比较。
肝提取分数(HEF)和相对血流(RBF)
解卷积肝响应曲线被关于HEF和RBF分析。HEF由Brown等人使用Tc-99-disofenin闪烁扫描术描述为对肝提取效率的测量,并且可以被理解为在示踪剂被直接注射到该肝的传入(afferent)血液供给而没有后续的循环的情况下,已经被提取的示踪剂的百分比。图3显示了来自使用Gd-EOB-DTPA的肝实质(肝提取(HE)曲线)的典型的脉冲响应。该HE曲线可以被分成血管阶段和肝细胞维持阶段(其描述了肝提取)。我们使用在注射示踪剂的时刻之后从420到1800秒的HE曲线数据点的单指数拟合来计算HEF。在420秒的起始点在肉眼观察(visual inspection)该解卷积HE曲线之后被选择。该拟合曲线、该肝滞留曲线(HRC、单指数衰退拟合曲线)然后被向后推断到血管峰值的时间,而HEF被定义为推断的HRC曲线和该HE曲线(也被显示在图3中)的血管峰之间的比值,
RBF(给出了肝中的血流的相对度量)被描述为HE曲线的初始峰值。RBF值被标准化为具有最高RBF的节段,即,具有最高RBF的节段被设置为100%。
图像分析
该输入函数由在该门静脉的种脐(hilar)所放置的所关注的区域(ROI)定义。由于病人的移动,在每个动态获取中调整输入函数ROI,从而该体素代表门静脉血。通过在每个肝节段(I到VIII,其中节段IV被分成IVa和IVb)中放置三个ROI来限定肝响应函数曲线。该ROI中的体素随着时间的相对增强被认为是那个ROI的实质响应函数。在90分钟时间段内使用相等间距(60秒)来内插数据点。图3显示了具有内插数据点的典型输入函数和实质响应函数。在放置ROI时尽可能小心地排除主血脉管和可视的胆管。节段是如同Strasberg SM.肝解剖和肝切除的术语(其内容全部通过参考并入本文):掌握肝(coming to grips with hepatic)Babel.J Am CollSurg 1997;184(4):413-434提议的那样定义和命名的。然而,在其它实施方式中可以使用其它分段方式。
使用用(美国密歇根州的Mathworks)写成的具体的内部软件,用TSVD和FA+尾两者针对每个ROI计算HEF和HBF。因此,每个ROI分别对HEF和RBF产生两个值。对于RSVD,静态截断阈值被设定在c=0.07。增加从0到π/2、初始高度在x(t)和y(t)的最后一个点的余弦函数以便用FA执行DA,且该尾的长度被设置为90分钟的总采样周期的长度的三倍。
使用相同的输入函数作为被用于节段(但其中每个肝体素代表一个响应函数)的ROI来计算HEF和RBF的参数映射。RBF总是被标准化到每个对象的最大RBF值并作为百分比呈现。为了最小化噪声的影响(主要来自病人的移动),通过在输入和响应函数曲线两者中施加七个点的滑动窗口滤波器而使用数据的低通滤波。
将输入和响应函数中的相对造影剂浓度计算为对数比,
其中c(t,ρ)是体素ρ中在时间t的相对示踪剂浓度。S0(ρ)是来自预对比图像(即,基线信号强度)的体素ρ的平均图像强度。S(t,ρ)是在时间t所测量的体素中的图像强度。
药动学分成隔室的(compartmental)建模
在隔室(compartment)模型中,随着时间在不同分区间穿过的底物(substrate)的分布被建模。该模型是基于一阶运动学的(即,浓度的时间导数与该底物的本身的浓度成反比)。如果该模型只由一个隔室组成的话,描述该系统的等式是一维一阶微分等式。
在药动学建模中,一个可以选择包括任意数量的相关隔室,底物在这些隔室之间流动。当超过一个隔室被建模时,该系统等式变成微分等式的系统。在我们的研究中,我们已使用三隔室药动学模型,其由Gambhir等人(J Nucl35 Med 1989;30(9):1507-1518)在使用Tc-IDA作为示踪剂的闪烁扫描术研究中所描述。该模型被显示在图16中,并且可以以数学方式描述为
y(t)=[C]v(t)+fSblood(t)
[C]=[1 1] [9]
在这种系统中,v(t)=(v1(t),v2(t))是代表肝实质和胆汁中的信号的矢量,y(t)是响应函数,而u(t)是到每个容积的流入量。术语f.Sblood(t)代表来自血液池的信号Sblood(t)的分数f,其从肝实质添加到该信号。参数f和{k12,k21,k3}(后面用kij表示)是该模型中的未知量。如图16中所示,k21表示从血液到肝隔室的流速常数,k12表示从肝隔室流到血池中的回流,从肝细胞流到胆小管的胆汁的肝内流动用k32描述,而从肝内流到肝外的胆汁隔室的流动用k3隔室描述。在数学上通过假定k32等于k3,并假定没有从胆小管到肝实质的回流而简化该模型。
在特定情况下,当输入函数是纯推注(bolus)或Dirac脉冲,响应函数y(t)就等同于输入函数x(t)。然而纯推注剂量假设是一种理想情况,但是响应函数y(t)可以作为脉冲响应h(t)和输入函数x(t)之间的卷积而计算,如同关于HEF的计算讨论的那样。在该隔室模型中,该脉冲响应函数可以以分析方式被描述为包含kij参数的两个指数函数的和:
为了达到模型参数的预估,使用迭代方法。在向kij和f赋予起始值之后,通过Eq 10来估计脉冲响应函数h(t),并使用来自门静脉的输入函数x(t)作为等式7中对血池信号的估计,通过计算输出函数yout(t)。另一方面,响应函数y(t)是在肝实质中的ROI中被测量的,而参数kij和f是通过迭代最小化均方差diff=(y(t)-yout(t))2确定的,参看图17。
为了增加找到整体最小值的可能性,使用一组10个随机化的起始值;只有当该方法在10次中有超过6次汇集于同样的最小值时,kij和f的值才被接受。
该算法因此产生了五个参数,k12、k21、k3、f和diff。在图16中限定三个传递常数kij,该因素f限定了源自于血池的ROI中的信号的分数(因此描述了该ROI中的灌注),而diff描述了相对于在该ROI中测量的那些,对于该汇集的响应曲线的良好拟合。
半定量分析
直接从实质的时间-强度曲线获得的半定量参数是最大相对信号强度(Cmax),到最大强度的时间(Tmax)、从Tmax到相对信号强度的百分之五和十的衰减的时间(分别是T5和T10)以及从0到5400秒的AUC。对于一些响应曲线,或者T10或者T5和T10两者,超出了最后的测量时间点,而且没有值被设定。Tmax、T5和T10是按秒测量的。因为在Gd-EOB-DTPA的情况下信号强度半时(TE)比在这种研究中所使用的90分钟总扫描时间长得多,通过双指数拟合来估算TE,该双指数拟合由下式确定
其中g(t)是拟合的曲线而拟合参数c2和Tu描述了对比度吸收(contrast uptake),而C1和TE描述了肝对比度排泄。TE和TU两者是按分钟计算的。如果包括整个响应曲线的话,该双指数拟合不总是汇集,因此以根据经验选择t=240秒作为该拟合的起始点。
统计分析
这三个节段ROI的平均HEF和RBF被认为是那个特定节段的产生的HEF和RBF。分别用这两种DA方法对HEF和RBF计算描述性统计(平均、标准偏差(SD)、差异系数(CV)、中间、最大、最小和范围)。研究产生180个HEF和RBF的成对观测值(20个对象,每个具有9个节段,且用TSVD和FA+尾分析每个对象)。使用非参数Wilcoxon匹配成对测试来比较这两种DA方法的中间HEF和RBF,并使用变量比测试(也被称为F-测试)比较这两种方法的SD。小于0.05的双侧(two-sided)p值被认为是重要的。该Mann-Whitney U-测试被用于比较非成对数据。
图12是一个实施方式的系统1900的示意图。系统1900适于基于计算机确定具有分泌或排泄功能的至少一个器官(比如肝和/或肾)的功能性评估。该系统包含用于处理所述人的四维(4D)图像数据组的单元,该数据组包含用于评估所述至少一个器官的所述功能的数据,其中所述4D图像数据是通过图像模态处理所述人的四维(4D)图像数据组而获得的,所述数据组包含用于对所述肝功能进行评估的数据,其中所述4D图像数据是通过图像模态获得的。在一个实施方式中,用于处理所述4D图像数据的所述单元被布置为执行解卷积分析(DA),该解卷积分析包含基于所述4D图像数据使用奇异值分解(SVD)的矩阵求逆。
在一个实施方式中,系统1900是一种基于计算机的系统,适于确定人的至少一个器官随时间变化的函数。该器官是具有分泌或排泄功能的器官,比如肝和/或肾。该系统包含处理单元,该处理单元被配置为处理通过图像模态获得的一组四维(4D)图像数据,并被配置为基于该组四维(4D)图像数据确定与该至少一个器官的该至少一个器官的每个容积单位的功能有关的参数的值。
通过将该参数的确定的值与以前确定的健康人群的参数的在值进行比较,促进对该器官的机能失调的诊断。
医学工作站1910包含常用计算机部件,比如中央处理单元(CPU)1920、存储器、界面等等。而且,它被装备有用于处理来自数据输入源的数据的适当软件,比如从MRI扫描获得数据。例如,软件可以被存储在能由医学工作站1910访问的计算机可读介质1930上。计算机可读介质1930可包含按照计算机程序1940形式的软件,该计算机程序1940包含合适的代码段190。医学工作站1910进一步包含例如用于显示渲染的显影(rendered visualization)的监控器,以及例如用于人工微调自动规划(否则由软件所提供)的合适的人机接口装置(类似键盘、鼠标等等)。该医学工作站可以是系统1900的一部分。
计算机程序1940能存储在计算机可读介质上,用于由计算装置(比如医学工作站1910的CPU 1920)进行处理,以确定至少一个分泌或排泄器官(比如人的肝和/或肾)随着时间的功能。计算机程序1930包含多个代码段,该代码段包含第一代码段190,其用于根据对通过图像模态获得的人的一组四维(4D)图像数据的处理,来确定与该至少一个器官的该至少一个器官的每容积单元的函数有关的参数的值。
因此,根据将确定的该参数的值与以前确定的健康人群的参数值的对比,使得在该器官的节段中对该器官的机能失调的诊断成为可能。例如,该参数是肝提取分数或输入相对血流。
图18到25以及25A和26B中分别显示了与来自健康人群的值的这种比较的示例。
上面描述的计算或虚拟规划的结果可以在医学工作站1910上的图形用户界面中提供给用户。
图13是一个实施方式的计算机程序的示意图。该计算机程序被布置为通过对具有分泌或排泄功能的至少一个器官(比如肝和/或肾)进行由计算装置进行的评估而进行处理,以便由所提供的计算机进行处理。该计算机程序可以嵌入计算机可读介质上并包含代码段190,该代码段190处理所述人的四维(4D)图像数据组,该数据组包含用于评估所述至少一个器官的所述功能的数据,其中所述4D图像数据是通过图像模态获得的,其包含执行解卷积分析(DA),该解卷积分析包含使用基于所述4D图像数据的奇异值分解(SVD)的矩阵求逆。
在图27中,显示了具有n个体素的ROI 400,所有体素具有不同的肝细胞和血管比例。对于每个体素计算HEF和irBF并绘图。通过线性回归,直线410被拟合到获得的数据点。用这种方式,计算并提供局部HEF和局部irBF(具有用于部分容积效应的补偿)。
其中所描述的用于节段或子节段肝功能、肝灌注和胆汁排泄功能的估计的方法和/或系统具有对于诊断、疾病进展的监控、治疗效力或治疗副作用的评估的益处的一些疾病、医学领域、过程和/或器官诊断包括:
肝病:
·急性肝炎
·慢性肝炎
·原发性硬化性胆管炎
·原发性胆硬化
·囊性纤维化
·对硬化/纤维化分级并监控疾病的进展
·利胆药对肝内胆汁郁积中的胆汁流入的效力的估计
-对肝的其它形式的医学或免疫治疗的影响的估计
-具有NAFLD和NASH的肥胖症
-具有肝功能衰弱的代谢综合症(metabolic syndrome)
-在做肝硬化手术的病人之间对肝细胞癌的检查过程中的肝功能监控
手术:
·肝内胆结石病
·为直肠癌(colorectal cancer)而进行的节段性肝手术的手术前和手术后肝功能的预测
肝次生肿瘤及其他肝的原发性和次发性肿瘤
·对阻塞性黄疸中胆汁流入的扩张效力或EST(括约肌切开取石术)
·肝内和肝外胆树的恶性和良性肿瘤中的胆汁流动的估计
·所有形式的肝肠支路的开放和效能的估计
·肝移植病人的移植状态的监控
肿瘤:
·化疗导致的实质性损伤(NASH、NAFLD、SOS)
示例1
对象
对20个健康的志愿者,10男10女,年龄从22到45岁,进行T1-加权的Gd-EOB-DTPA增强DHCE-MRI。在该研究中在内涵物(inclusion)进行常规血清肝功能试验。试验对象没有肝胆疾病、以前的肝胆手术或嗜酒症病史。
协议
使用具有Philips四通道感测本体线圈的Philips Intera 1.5T扫描仪(荷兰,Best)收集数据。使用T1-加权的3D扰相梯度回波(spoiled-gradient-echo)脉冲序列(循环时间/回波时间/翻转角度4.1ms/2.0ms/10度,可视区域=415毫米,矩阵分辨率256x192,40个切片,切片厚度为10毫米且感测因子R=2)。在单次呼吸中在41个不同的时间点(每个获得的容积有12秒的扫描时间)对该容积成像。三个容积是在对比之前获得的,用于基线计算,接下来38个容积是采样间隔逐步增加的。根据对象的体力、数据获取限制和测试物质的动态来选择采样密度。0.1ml/kg Gd-EOB-DTPA的计量0.25mmo/ml的被注射到右侧前肘脉中,与第四个获得的容积的起点相一致。该对比度是使用强力注射器(power injector)(匹兹堡Medrad的SpectrisMR注射器系统)、以2ml每秒的注射速率注射的,然后立即以同样的注射速率注射20ml盐水(NaCl 0.9%)的丸剂(bolus)。
结果
所有的对象有正常的血清肝功能试验并且没有肾机能不全的迹象。图9将模拟的结果显示为TSVD和FA+尾技术间的SD对比。在较高的SNR值处TSVD比FA+尾执行得更好。然而,当数据包含更多噪声时,TSVD更稳定,具有好得多的标准偏差。
表1显示了用两种DA方法的HEF和RBF结果的概略统计。
表1
1)Wilcoxon匹配对试验 2)变异率试验
从这20个试验对象产生的HEF和RBF结果以图形方式呈现在图14A和14B中,而HEF和RBF在节段水平上的分布显示在图15A和15B中。平均的ROI尺寸是31.9(SD 21.6)个体素。
使用这两种方法,有关HEF或RBF的总结果没有明显的差异(对于HEFp=0.524,而对于RBFp=0.331),但是TSVD产生了更小的SD和更小的CV,尽管SD的差异并不明显(对于HEF,p=0.152,而对于RBF是0.458)。对于左半边肝和右半边肝有中值HEF的差异(使用TSVD对于左边是0.196而对于右边是0.218,而将FA+尾使用DA是0.194vs.0.224),但是只有当使用FA+尾技术时,该差异才明显(使用TSVD时p=0.14vs.使用FA+尾时p=0.011)。在左半边肝和右半边肝之间还有RBF的差异,其中左半边的RBF明显更低,其中中值RBF在使用TSVD时是79.1%而在使用FA+尾时是81.2%。右半边的对应值分别是94.0%和88.4%。这种差异对于两种DA方法都非常明显(使用Mann-Whitney U-测试时p<0.001)。
图10显示了对于一个试验对象中的肝的水平分段平面之间上方的切片的HEF的参数映射(图10A、B)和RBF的参数映射(图10C、D)。通过肉眼观察,HEF看上去在遍及该切片是同源的。该参数HEF映射中的高值(即,接近于100%)是包含血管的高水平的体素的结果,不被认为是反应了该肝的功能。超过100%的值被认为是伪像并被排除。所有的RBF值分别根据每个对象的最高流值计算比例。
在此示例1中,发现使用用Gd-EOB-DTPA做示踪剂的DHCE-MRI来在节段水平上评估HEF和RBF是可行的。还发现对于活体内的解卷积分析,TSVD执行的比FA+尾更好。对于本应用领域来说,TSVD在计算方面需求较小。计算机模拟还显示出,使用TSVD进行DA对噪声数据更不敏感,其中在较低SNR水平具有明显更低的SD,因此TSVD是DA的优选选择。
在对健康对象的闪烁扫描研究中,当使用具有接近总肝清晰度的IDA-类似物,HEF是约100%。在此示例1中,略高于20的平均HEF可以很好地反映已知的事实,即Gd-EOB-DTPA比IDA-化合物具有更低的肝亲合力,具有约50%的肝清晰度。因为Gd-EOB-DTPA具有不同的肝针对性,HEF可能不是描述使用Gd-EOB-DTPA的肝细胞吸收的最佳参数。
一个有趣的发现是左半边肝和右半边肝的肝节段之间的HEF和RBF的观察到的差别,参看图15A和15B。
对象内变化可部分由90分钟的获得期内的移动伪像来解释。这与ROI的部分容积效应结合起来可能导致噪声数据,其中肝ROI未必反映全部动态容积中的肝实质。高分辨率肝功能试验中的移动伪像应当被最小化以提高数据质量。另一方面,HEF中的对象内变化,可以是用以前的技术不可能检测到的真实现象。
在利用DA的每个研究中,输入函数对于获得的结果是至关重要的。肝具有双血管供应,其具有来自门静脉的静脉输入和来自肝动脉的动脉血。我们选择将来自该门静脉中的ROI的随着时间增加的曲线用作输入函数。这样做的理由主要是生理学上的,因为大约75%的传入该肝的血流源自该门静脉。另一个原因是动脉输入函数具有极短的峰值,而且使用在此示例1中的时域分辨率,我们根据经验发现,我们经常错过该动脉峰值,导致我们的对象之间的动脉输入函数中的最大峰值的恼人的差别。该门峰值随着时间稍微散布得更开,而观察到的峰值的差别小得多。在第一个三分钟期间每分钟获得三个容积。
在T1-加权的对比度增强DHCE-MRI中,信号强度依赖于T1-张弛(relaxation)时间,更高的Gd-DTPA浓度减少了T1-张弛时间并增加了图像信号强度。已经显示,图像强度和Gd-DTPA浓度之间的关系对于稳态MRI脉冲序列来说是非线性的,比如此研究中使用的扰相梯度回波。然而,当T1-张弛在40ms到2600ms范围内时,显示该MRI信号以随着T1-张弛的缩短大致地指数增加。
我们所有的测量在这个范围内被估计,使得等式5a对相对造影剂浓度有良好的逼近。
示例2
对具有原发性胆汁性肝硬变(PBC)的病人的研究
对象
对20个健康的志愿者(10男10女),并对具有确诊的PBC的20个病人执行T1-加权的Gd-EOB-DTPA-增强DHCE-MRI。对健康的志愿者在该研究中在内涵物执行常规血清肝功能试验,而对病人来说从他们的临床图表中所记载过的最近的就诊中记录。健康的志愿者没有肝胆病、以前的肝胆手术或嗜酒病史。所有的对象都被要求在检验前禁食至少四个小时。对每个病人,相关的临床资料被记录并且它们与来自肝功能试验的结果一起被用于计算CPS、Mayo危险分数和MELD分数。
MR过程
使用Philips Intera 1.5T扫描仪(荷兰的Best)执行T1-加权的Gd-EOB-DTPA-增强DHCE-MRI,该扫描仪具有根据示例1中的协议的飞利浦四通道感测本体线圈。使用截断奇异值分解(TSVD)进行解卷积分析。像上面描述的那样计算HEF和irBF。通过评估从峰值到2700秒的肝提取曲线以下的面积来定量计算AUC。像上面描述的那样计算半定量参数(SQP)和药动学传递常数,并且还按照从0到5400秒的实质响应曲线以下的面积半定量计算AUC。将Mann-Whitney U-测试用于明显性试验并将明显性水平设定为α=0.5。每个病人和对照(control)中的所有节段产生观测并且在统计分析中,所有观测被认为是独立的观测,即便它们源自一个个体。因此,该研究从对照(control)中对上面文提到的参数中的每一个产生180个观测,并根据PBC病人产生108个观测。
结果
该研究包括(规划的总共20个病人中的)12个病人(1男11女)。表2中呈现了病人的特征、血清肝功能试验(LFT)的结果和相关的临床信息。
表2
*=比例测试,**=学生t-测试,
***=在MR图像上判断为硬化的实质
有PBC的病人通常比对照组(control)年龄大,而且性别分布不同,正如预期。这两个组之间在PK或胆红素水平方面没有明显差异,但是在PBC病人中白蛋白水平明显更低。在病人中AST、ALAT和碱性磷酸酶都明显更高。表3显示了定量参数的结果,而表4显示了半定量参数的结果。
*Mann-Whitney U-试验
在PBC病人中HEF明显更低而irBF明显更高,但是与对照组相比吸收转移k21不同。传送率常数k12和k3在病人中比在对照组中更高,正如标示ROI中的血液的分数的因子f那样。各组之间在拟合良好度参数diff方面没有明显差异。关于半定量参数,在最大强度(Cmax)、排泄半时(TE)或曲线下面积(AUC)方面没有明显差异。在PBC病人中到达最大强度的时间(Tmax)明显更长,但是排泄参数T5和T10更短。HEF和AUC(定量计算的)随着儿童得分的增加而更低。
在这个研究中,我们在有PBC的病人中发现明显更低的HEF,正如预期的一样,而该差异看上去随着疾病的严重程度的增加而增加,如同上面描述的。已经知道,肝硬化导致动脉血流的增加和门血流的减少。或许硬化的肝实质中的增加的动脉峰值可以解释在此研究中注意到的irBF中的差异。因为PBC导致微胆管的闭塞,人们期望到达最大增强的时间更长,因为钆示踪剂在肝细胞中随着时间的增长而累积。当我们只考虑具有肝硬化的形态学证据的病人时,我们找到相对于健康对照组的甚至更大的差异。在图18-21中,在节段水平上显示了这一点,其中在健康对照组和在该研究中获得的MR图像上有肝硬化迹象的5个病人之间比较定量参数。人们期望k3参数在PBC病人中比在对照组中更低,然而情况似乎不是这样。所使用的研究参数似乎能够检测有关吸收差异的实质函数中的恶化,然而不能量化胆排泄中的差异。
示例3
对具有原发性硬化性胆管炎(PSC)的病人的研究
对象
对20个健康的志愿者(10男10女)和具有确诊的PSC的病人执行T1-加权的Gd-EOB-DTPA-增强MRI。对健康的志愿者,在研究中在内涵物执行常规血清肝功能试验,并且对病人来说根据他们的临床图表被记载的最近的就诊来记录它。健康的志愿者没有肝胆疾病、以前的肝胆手术或嗜酒病史。要求所有对象在检查之前禁食至少4个小时。对每个病人,相关的临床数据以及与来自肝功能试验的结果一起被记录。它们被用于计算CPS、Mayo危险得分和MELD得分。
MR过程
使用Philips Intera 1.5T扫描仪(荷兰的Best)执行T1-加权的Gd-EOB-DTPA-增强的MRI,该扫描仪具有根据示例中概述的协议的Philips四通道感测本体线圈。使用截断奇异值分解(TSVD)执行解卷积分析。像上面描述的那样计算HEF和irBF。通过评估从峰值到2700秒的肝提取曲线以下的面积定量计算AUC。像上面描述的那样计算半定量参数(SQP)和药动学传递常数,并且按照计算从0到5400秒的实质响应曲线以下的面积以半定量计算AUC。使用Mann-Whitney U-测试以进行明显性测试并将明显性水平设置为α=0.5。每个病人和对照组中的所有节段产生观测,而在统计分析中,所有的观测被认为是独立的观测,即使它们源自一个个体。因此,该研究对上面提到的参数中的每一个从对照组中产生180个观测,而从PSC病人中产生108个观测。
结果
该研究中包括(规划的20个病人中的)12个病人。表5中总结了所包括的病人和对照组的人口和临床参数。
*=比例试验,**=学生t-试验。
***=在MR图像上判断为硬化的实质
有PSC的病人通常比对照组更老,而性别分布不同。两个组之间在PK或胆红素水平方面没有明显差异,但是在PSC病人中白蛋白水平明显更低。在病人中AST、ALAT和碱性磷酸酶都明显更高。表6和7分别显示了定量和半定量参数的结果。
*Mann-Whitney U-试验
在PSC病人中HEF明显更低而irBF明显更高,且在病人中定量计算的AUC明显更小。吸收传递常数k21在各组之间没有不同。传递率常数k12和k3在病人中比在对照组中更高,但是指示ROI中的血液中的分数的因子f没有不同。在各组之间在拟合良好度参数diff方面有明显差异,其中在病人中有大体更好的拟合。在半定量参数方面,在最大强度(Cmax)、排泄半时(TE)或曲线下面积(AUC)方面没有明显差异。到最大强度(Tmax)的时间在PSC病人中明显更长,但是排泄参数T5和T10更短。
这个研究中的病人人群具有有低MELD和Mayo得分的相对温和的疾病。只有一个病人是儿童B。尽管如此,在指示实质功能的差别的示踪剂的肝摄取中的明显差异可以使用具有Gd-EOB-DTPA的T1-加权的DHCE-MRI来检测。如果针对健康志愿者的结果绘制HEF和AUC的话,看上去如果疾病有更高得分的话AUC参数倾向于更少的吸收,但是这对HEF不明显。如果我们在节段水平上绘制该结果,并且只绘制具有异常对比度增强图案的节段,在正常和异常实质之间,在HEF和AUC中有值得注意的差异(图22和23)。该传递率常数k21和k3看上去在被更严重影响的实质中更高,但是对此的解释是难以理解的(图24和25)。一个有趣的发现是PSC病人中肝实质的小但明显的高灌注,这通过irBF的增加可以理解。理论上,这可能是肝实质中正在发生的发炎过程的结果,或者有可能是硬化或纤维化实质的动脉化(arterialisation)。至于Cmax,在各组之间没有不同,对此可以有一些解释。当目视检查来自具有异常外观的实质的一部分的实质响应曲线时,很显然它不同于健康志愿者的实质响应曲线(图26)。Tmax在病人中明显更高,但是排泄似乎更快(表7),具有更低的T5和T10。或许对此的一个解释是熊脱氧胆酸(ursodeoxycholicacid)治疗(当时2/3的病人正在在内涵物的时间上)的促胆汁分泌效应。
进一步的示例是在胆管中注入扩张剂(stent)以打开被阻塞的胆管。通过注入前和注入后的器官功能或胆汁流动的对比,提供了对扩张剂效力的估计。
总之,公开一种用于在节段水平上对肝细胞功能的每个容积的估计的新颖的方法和系统。在实施方式中,在健康的人类志愿者中使用DHCE-MRI,比如具有动态Gd-EOB-DTPA增强的MRI。呈现了用FA+尾和TSVD两者应用DA的数学模型,而不是使用概略参数。TSVD(与基于傅里叶的DA相比,对噪声数据稍微更不敏感)是用于在肝功能试验中用DHCE MRI获得的数据的解卷积的优选方法。
该方法和/或系统还可用于实现、提供或执行治疗的虚拟规划,如同上文所进一步描述的。
该方法和/或系统还可被应用于具有分泌或排泄功能的其它器官,比如例如胎盘(placenta)、消化系统或胰腺。
该方法和/或系统可被应用于同时确定若干个器官的功能。这些器官之间的函数的分布可以被计算并进一步处理。
本领域的技术人员将会理解,方法明可以被实现为装置、系统、方法或计算机程序产品。相应地,本发明可以采用完全硬件实施方式、软件实施方式或软件和硬件方面相结合的实施方式。此外,本发明可以采用计算机可用存储介质上的计算机程序产品的形式,其中在该计算机可用介质上嵌有计算机可用程序代码。可以使用任何合适的计算机可读介质,包括硬盘、光存储装置、传输介质,比如支持因特网或内部网的那些、或磁存储装置。
上面参考具体实施方式描述了本发明。然而,不同于上述实施方式的其它实施方式同样有可能在本发明的范围内。本发明的范围只由所附权利要求限制。