CN101001569B - 评估组织特性和检测病人整体心脏损害的系统 - Google Patents

Abstract

评估组织特性包括识别受损组织或天然组织成分比例的变更，例如心脏中正常心肌和纤维组织数量的变化，识别细胞外成分或液体(例如浮肿或细胞外基质蛋白)数量的增加，或检测肿瘤细胞或炎症介体渗出到患者例如人体的所关注组织中，本发明提供的评估方法是：例如在病人注射造影剂后从第一次采集获得包括所关注区域的组织的第一影像，并例如在病人注射造影剂后在第二次随继的采集中获得包括所关注区域的组织的第二影像。随继的采集可以例如在一段时间之后获得，以便确定在这段时间中是否有损害发生。所关注区域可以包括以下各项中的至少一项，例如，心脏、血液、肌肉、脑、神经、骨骼、骨骼肌、肝、肾、肺、胰腺、内分泌、胃肠和/或泌尿生殖器组织。确定第一影像和第二影像的所关注区域的整体特性，从而可以比较第一影像和第二影像的整体特性，以便确定整体组织特性改变的潜在可能。这种比较可以包括比较均值、平均特性、直方图形状例如倾斜和峰态、或直方图中的强度分布。

Description

评估组织特性和检测病人整体心脏损害的系统

有关申请和优先权要求

本申请与2004年2月6日提交的序列号60/542,547的美国临时专利申请有关并要求其优先权，所述专利的公开内容全部通过引用结合在本文中。

技术领域

本发明涉及诊断，而更具体地说，涉及整体组织特性如整体组织损伤的检测。

背景技术

阿霉素是一种蒽环类抗菌素，由土壤微生物分离得来。它的抗肿瘤效果与和topoisomerase-2酶的相互作用及产生DNA双条断裂有关。此外，这种药物会产生具有高度细胞毒性的细胞内自由基。阿霉素被认为是有最广泛作用的抗肿瘤药物之一。阿霉素不仅在乳腺癌、软组织肉瘤和其他的实体肿瘤的现代治疗方面被认为是典型的重要治疗组成成分之一，而且在治疗急性白血病、何杰金氏病、非何杰金氏淋巴瘤与许多儿童时期癌症方面，它也被认为是联合化疗的重要组成成份。因此，对于许多晚期癌症病人来说，阿霉素用作他们医学疗法的重要部分。

在成人和儿童中，阿霉素疗法的应用通常受它对心脏毒性蓄积量的限制。在那些接受超过500-550mg/m²剂量的病人中，会发生伴有严重充血性心力衰竭的不可逆心肌病变的严重危险。不幸的是，促成充血性心力衰竭的剂量范围变化很广(在一份超过七年之久的1487病例的研究报告中，变化范围为30-880mg/m²)。那些老年或休息时就有左室收缩功能轻度减低(左室射血分数〔LVEF〕≤50％)的病人是最危险的。在西方工业化国家里，通常是患有癌症的老年病人和有某种程度的基础心脏疾病的人们常常最需用阿霉素治疗，然而对于他们，由于潜在的心肌毒性可能被停止医疗。

由阿霉素引发的心肌病变的检测方法之一，是伴随做心导管检查以测量左和右心室压时，做活体组织检查。不幸的是，这种方法涉及侵入性过程，而且随着时间的推移不适合做重复的测量。放射性核素心室造影术也广泛用来筛选那些处于阿霉素引发心肌病变发展危险中的病人。LVEF进展减少10％或更多的病人，或在治疗过程中射血分数跌落低于50％的病人都处于发展为不可逆心肌中毒的最最危险中。虽然这些信息作为潜在的筛选术是有用的，对某些病人而言，但所观察到的LVEF的跌落太晚，不能避免不可逆的心肌病变的发展。正因为这个原因，对于接受化疗的病人阿霉素的总用量可能受到不适当的限制。对许多病人而言，重要的是常常在药物疗法获得最大利益之前阿霉素治疗就停用了。一种能用于精确检测那些会继续发展心肌毒性病人的非侵入性的、广泛可用的方法，会有明显的临床效用。

在过去的7年里，研究人员已确定了磁共振成像检查能用于在持续心肌损害的病人中识别左室内的坏死组织。所述技术包括在静脉内注射钆螯合物后用非选择性预备射频脉冲获得梯度-回声脉冲序列。在心肌坏死的区域中，注射造影剂后20分钟收集的影像上发生升高的信号密度，这与因心肌细胞膜破坏及毛细管渗透性增强而细胞外容积扩大相一致。这种方法已被用来识别患持续的急性或慢性Q波(ST段抬高)的病人中的透壁性心肌细胞坏死，和在患持续无Q波(无ST段抬高)的心肌梗塞的病人中的心内膜下(非透壁性)损伤。由磁共振成像检查发现的坏死量，与在冠状动脉血管再形成后收缩期增厚的恢复，呈现一种相反的关系。在注射造影剂后20分钟没有钆的极度增强，则与在持续的ST段抬高或无ST段抬高的心肌梗塞后心肌的生存能力和左室收缩的随后改善相关联。虽然有些人觉得延迟增强技术可能在急性梗塞方面过高地估计了心肌细胞坏死的区域，最近，一项在动物中作标签的研究指出延迟增强技术确实可以识别心肌梗塞(MI)后早期心肌细胞的坏死。据信在梗塞边缘地带，由于来自邻近有生命区域的束缚(tethering)缘故，死亡细胞可以迁移。

利用磁共振成像，可以将心脏结构成像，且对左室功能可直接进行评估，具有时间与空间的高分辨率。由于声窗不限制影像的采集，特别对于巨大的或不一般的体质，磁共振成像的应用程度很高。这种影像增高的清晰度允许研究者们进行左室结构与功能的定量测量，比放射性核素和超声检查方法获得的定量测量有更高的精确度。用90％的功率就能检测出左室功能减低的病人中5％的LVEF改变，p值为0.05，在平行研究设计中样本大小为每组5例病人。根据操作者的经验，在相同的研究设计中，同样5％的LVEF改变则要求大于每组100病人的超声心动描记术的评估。同理，用延迟增强磁共振成像技术获得的提高空间分辨率(1mm²像素大小)可以检测出微小梗塞，而这种微小梗塞迄今只能从例如血清样本中检测到心脏酶升高才能评估，但利用放射性核素或超声心动描记术的技术却无法看到。

在延迟增强成像过程中，将造影剂注射进病人，在造影剂已有机会分布到需成像的区域中，使造影剂存留在受损组织中而不存留在健康组织中以后才采集影像。这种延迟增强成像可用来例如识别心肌梗塞，因为梗塞区域的坏死组织会保留造影剂，而造影剂会被健康组织清除掉。这样，梗塞会呈现为较高密度的局限区域。按照惯例，延迟增强成像可用于识别诸如心脏组织、脑组织、神经组织等组织中的组织损害的局限区域。

发明内容

本发明的实施例提供用于评估组织特性的方法、系统和/或计算机程序产品，包括识别受损组织或天然组织成分比例的变更，例如心脏中正常心肌和纤维组织数量的变化，识别细胞外成分或液体(例如浮肿或细胞外基质蛋白)的数量增加，或检测肿瘤细胞或炎症介体渗出到患者例如人体的所关注的组织中，其方法是：例如，在病人注射造影剂后从第一次采集中获得包括所关注区域的组织的第一影像，并例如在病人注射造影剂后在第二次随继的采集中获得包括所关注区域的组织的第二影像。随继的采集可以例如在一段时间之后获得，以便确定在这段时间中是否有损害发生。所关注区域可以包括以下各项中的至少一项，例如，心脏、血液、肌肉、脑、神经、骨骼、骨骼肌、肝、肾、肺、胰腺、内分泌、胃肠和/或泌尿生殖器组织。确定第一影像和第二影像的所关注区域的整体特性，从而可以比较第一影像和第二影像的整体特性，以便确定例如由所关注区域组织的整体损害引起整体组织特性改变的潜在可能。这种比较例如可以包括比较平均特性、直方图形状例如倾斜和峰态、或直方图中的强度分布。

在本发明又一些实施例中，整体特性是所关注区域的像素/体素特性，所述特性基于所关注区域中基本上全部像素/体素。整体特性也可以是所关注区域中像素/体素的平均强度。所关注区域中的组织可以是心脏组织，脑组织和/或神经组织中的至少一种。第一影像和第二影像可以是磁共振成像(MRI)影像。

虽然本发明的某些实施例是针对检测整体组织特性，例如患者如人体的整体损坏来加以说明的，但本发明的附加实施例可以包括检测脊椎动物或无脊椎动物、重建组织和/或合成组织的整体损坏。因此，本发明的某些实施例不应被认为限于检测病人的整体损害。

本发明的特定实施例提供用于检测病人整体心脏损害的方法、系统和/或计算机程序产品。在病人注射造影剂后，获得第一心脏影像。在病人注射造影剂后，还要获得第二心脏影像。确定第一心脏影像的强度和第二心脏影像的强度，并将第一心脏影像的强度和第二心脏影像的强度进行比较，以便确定整体心脏损害的潜在可能性。在本发明的某些实施例中，影像强度的增强表示可能存在整体心脏损害。

在本发明的又一些实施例中，第一心脏影像和第二心脏影像是磁共振成像(MRI)影像和/或X射线计算机体层摄影(CT)影像。而且，第一心脏影像的强度和第二心脏影像的强度可以是各影像的平均强度。

在本发明的附加实施例中，还获得对应于第一心脏影像的心脏外所关注区域的第一影像。利用心脏外所关注区域的第一影像的数据对第一心脏影像上正常心肌组织中像素强度的变化进行校正。同理，获得对应于第二心脏影像的心脏外所关注区域的第二影像并使用心脏外所关注区域的第二影像的数据对第二心脏影像上正常心肌组织中像素强度的变化进行校正。利用校正的第一心脏影像和校正的第二心脏影像获得第一心脏影像的强度和第二心脏影像的强度。例如，可以相对于没有造影剂情况下正常心肌组织来测量由造影剂的存在引起的增加的亮度的量度。在所有对象中可以不把正常心肌抑制到同样的暗度，且这种改变可以计数。

本专业的技术人员会理解，根据本公开内容，可以以方法、系统、和/或计算机程序产品的形式提供本发明的实施例。

附图说明

图1是按照本发明实施例的MRI系统的方框图。

图2是按照本发明实施例的数据处理系统的方框图。

图3是按照本发明实施例的数据处理系统的方框图。

图4A和4B是说明按照本发明某些实施例的操作的流程图。

图5是说明按照本发明某些实施例的操作的流程图。

图6是左心室3个短轴平面的3D描绘。

图7为左心室的中(中间平面)短轴视图的延迟增强MRI影像，带有相应的强度直方图。

图8是所关注区域(ROI)内体素的强度直方图。

图9是用于研究病人的自动相关性测量的图形。

图10是两个单独的病人的影像和平均体素强度。

图11是对一个病人相隔25天采集的中短轴视图。

图12是用于再现断层位置的影像计划软件的屏幕捕获。

具体实施方式

现参阅附图对本发明作更详细的说明，附图中示出了本发明的实施例。但本发明不应被认为限于本文所提出的实施例。而是，提供这些实施例是使本公开内容更透彻和完全，以将本发明的范围充分传达给本专业的技术人员。在全文中相同的编号指相同的元件。在本文中词语“和/或”包括一个或多个关联所列词语的任何以及全部组合。

本文中所用的术语仅是为了说明特定的实施例，而不是旨在限制本发明。本文中用到的单数形式“一”、“一个”和“所述”应也包括多数形式，除非上下文另有明确说明。还应理解，在本说明书中所用的词语“包括”和/或“包括”说明存在所述的特征、总体、步骤、操作、元件、和/或部件，但并不排除附加一个或多个其它特征、总体、步骤、操作、元件、部件、和/或它们的组合。

应理解，虽然在本文中用词语第一，第二等来说明各种元件、部件、区域、层和/或段，但这些元件、部件、区域、层和/或段不应受这些词语的限制。这些词语仅是用来使元件、部件、区域、层或段与另一区域、层或段相区别。因此，以下讨论的第一元件、部件、区域、层或段也可以称之为第二元件、部件、区域、层或段，这都不背离本发明的内容。

除非另有定义，本文所用的所有术语(包括科技术语)都具有本发明所属专业一般技术人员通常理解的相同含意。还应理解，有些术语，例如在常用词典中所定义的那些术语，应解释为具有与相关技术和本公开内容上下文中的意义相一致的意义，而不应解释为理想化或过于正式的意义，除非在本文中作这样明确的定义。

如本专业的技术人员所理解，本发明可以作为方法、系统、和/或计算机程序产品来实施。因此，本发明可采取以下形式：全部硬件实施例，全部软件实施例，或将软件和硬件方面加以组合的实施例，在本文中统称为“电路”或“模块”。而且，本发明可以采取在计算机可用存储介质上的计算机程序产品的形式，计算机可用程序代码在介质中实施。可以使用任何适用的计算机可读介质，包括硬盘、CD-ROM、光存储装置、传输介质，例如那些支持因特网或内部网的传输介质、或磁存储装置。

实现本发明操作的计算机程序代码可以用面向对象编程语言来写入，例如

Smalltalk或C++。但是，实现本发明操作的计算机程序代码也可以用传统的，过程式编程语言来写入，例如“C”编程语言。程序代码可以完全在用户的计算机上执行、部分在用户的计算机上执行，作为单独的软件包，部分在用户的计算机上且部分在远程计算机上或全部在在远程计算机上执行。在后者情况下，远程计算机可以通过局域网(LAN)或广域网(WAN)连接到用户的计算机上，或此连接可以连接到外部计算机上(例如使用因特网业务提供商通过因特网来连接)。而且，用户的计算机、远程计算机、或二者，可结合到其它系统中，例如MRI和/或X射线计算机体层摄影(CT)系统中。

以下参阅按照本发明实施例的方法、设备(系统)以及计算机程序产品的流程图图示和/或方框图对本发明加以说明。应理解，流程图图示和/或方框图的每个框，以及流程图图示和/或方框图中的框组合可以用计算机程序指令来实现。可以将这些程序指令提供到通用计算机、专用计算机的处理器上或其它可编程数据处理设备上以产生一种机器，这样通过计算机处理器或其它可编程数据处理设备执行的指令创建了用于实现在流程图和/或方框图各框所规定的功能/动作的构件。

这些计算机程序指令也可以存储在能引导计算机或其它可编程数据处理设备以特定方式工作的计算机可读存储器中，使得这样存储在计算机可读存储器中的指令产生了一种制造品，包括实现在流程图和/或方框图各框所规定的功能/动作的指令工具。

还可以把这些计算机程序指令装入计算机或其它可编程数据处理设备上，使一系列操作步骤能在计算机或其它可编程数据处理设备上进行，以产生由计算机实现的过程，这样在计算机或其它可编程数据处理设备上执行的指令提供用于实现在流程图和/或方框图各框所规定的功能/动作的步骤。

MRI过程已可以确定地用于识别患继发于冠状动脉硬化的缺血性心肌病的病人的心肌细胞损害和LVEF。这些过程可以识别局部心脏损害。但据信这种非侵入性成像尚未用于识别患继发于化疗的心肌病病人的整体心脏损害。心肌细胞损害的早期检测可以提供一个机会来调节药物剂量并降低和/或减少与化疗相关联的心脏毒性的效果。这样，可以在没有心肌细胞损害的情况下给病人作最大剂量的化疗，且化疗药物的所需疗效能更充分实现。虽然本发明的实施例在阿霉素治疗方面特别有用，本发明的实施例也可用在其它化疗或治疗方案中，和/或欲检测整体心脏损害的诊断环境中。

本发明的实施例利用在所关注区域影像中所关注区域整体特性的比较来检测例如可能因损害而引起的组织特性的改变。所关注区域整体特性是所关注区域的一种特性，它基于所关注区域中全部或基本上全部的像素/体素的一种或多种特性。因此，在本发明的某些实施例中，整体特性可以与所关注区域中像素的位置基本上无关。整体特性的实例可以包括(但不限于)所关注区域中的像素/体素特性的统计分析，例如平均强度、强度值直方图或其它统计分析。使用影像整体特性的比较可以检测损害图形是随机的和/或是以所比较影像的分辨率不能检测的损害。本发明的实施例还可以使用整体特性不仅检测对一个区域的损害，也可检测在天然组织内的正常比例中未被发现的材料异常累积。本发明的实施例还可以用于分子成像策略，例如将具有分子识别部位的造影剂引入组织区域并量化目标或分子过程的存在。因此，本发明的特定实施例可以用于检测癌变、炎症、肿大或浮肿、斑痕组织等。而且，本发明的实施例可用来定义在有机系统中的组织内起作用的代谢途径。本发明的特定实施例提供在注射造影剂后使用非侵入性成像检测整体心脏损害。适用于本发明实施例的非侵入性技术包括磁共振成像(MRI)、超声、X射线计算机体层摄影(CT)、单光子发射计算机体层摄影(SPECT)和/或正电子发射体层摄影(PET)。可以在第一或基线影像和第二影像之间进行比较，对影像的对比度进行分析以便检测整体心脏损害的存在。在本文中使用的术语影像是指一种空间信号，可以对其进行评估以获得信号强度的所需量度。

在本文中使用的术语“整体损害”是指组织组成和/或功能的改变，这种改变具有基本上是随机分布的图形的形式和/或具有不能以对其进行分析以检测损害的影像的分辨率检测的图形的形式。因此，例如，“整体心脏损害”可以是指心脏损害和/或用纤维组织例如斑痕组织来代替天然心肌组织，这导致的坏死和/或纤维化，所述坏死和/或纤维化具有基本上是随机分布的图形的形式和/或具有不能以对其进行分析以检测损害的影像的分辨率检测的图形的形式。可以通过按照本发明实施例的强度分析来检测的整体心脏损害包括，例如：病毒性心肌病、酒精性心肌病、产后心肌病和/或原发性扩张性心肌病。整体损害还可以包括不成比例数量的其它异常，例如水肿(过量液体)、纤维化(斑痕组织)等等。因此，本发明的实施例可以提供对整体异常组织的检测。

适用于本发明实施例的造影剂可以包括顺磁性镧系元素螯合物，和/或连接到高分子的顺磁性镧系元素，例如钆DPTA。用于灌注成像的MR造影剂的其它实例包括使用含有氧化铁或镝的敏感剂(susceptibility agent)，通过在血液和细胞内室之间引起磁矩的大幅波动而将局部不均一性引入磁场。也可在引入会引发心肌病的其它药物例如可卡因和/或酒精后进行成像。这些波动导致邻近氢核子的T2-星缩短，引起信号强度的损失。在本发明的特定实施例中，对每个影像使用同一种造影剂。

此外，本发明的某些实施例可以提供不注射造影剂的对比度/强度分析。例如，灌注成像的另一实例是通过T2预备的真实FISP，或3D-T2-加权的序列策略使用取决于血氧水平(BOLD)的心脏成像作不注射造影剂的心肌灌注或损害的评估。其它技术使用内生的对比度，包括自旋标记和磁化转移对比度。因此，在本发明的某些实施例中，不注射造影剂也可检测所关注区域的整体特性。

按照本发明实施例的示范系统10示于图1。如图1所示，强度分析/MRI系统10包括MRI采集系统11，它可以包括MRI控制系统电路12、MRI脉冲激励系统电路14以及MRI信号测量系统电路16。MRI控制系统电路12控制MRI采集系统11的工作，以便在病人的心脏周期或其一部分期间获得并提供MRI影像。MRI控制系统电路12还可以将所采集的影像加以组合并传送到工作站20或其它这种数据处理系统上供进一步分析和/或显示。工作站20可以在MRI室内，或者可远离MRI室。控制MRI脉冲激励系统电路14和MRI信号测量系统电路16以便采集能提供病人心脏的MRI影像的MRI信号。

传统的MRI系统，例如由通用电气医疗系统、西门子、菲立普、Varian、Bruker、马可尼、日立和东芝提供的那些MRI系统可用来提供在注射造影剂后收集的所需MRI影像帧。

虽然示范性强度分析/MRI系统示于图1并以特定的功能和/或操作分工进行说明，但本专业的技术人员会理解，也可以使用其它功能和/或操作分工而仍然可从本发明的内容获益。例如，MRI控制系统电路12可以和MRI脉冲激励系统电路14或MRI信号测量系统电路16组合。因此，本发明不应被认为限于特定的体系结构或MRI功能/操作的分工，而是应涵盖能够进行本文所述操作的任何体系结构或MRI功能/操作的分工。

图2示出数据处理系统230的示范性实施例，适用于按照本发明的实施例提供工作站20和/或MRI控制系统电路12。数据处理系统230通常包括输入装置232，例如键盘板或键盘、显示器234、以及存储器236，它们都与处理器238通信。数据处理系统230还可以包括扬声器244以及I/O数据端口246，它们也与处理器238通信。I/O数据端口246可用来在数据处理系统230和另一计算机系统或网络之间传送信息。这些部件可以是例如在许多传统数据处理系统中所使用的传统部件，可以将它们配置成按照本文中所述的那样工作。

图3是数据处理系统的方框图，示出按照本发明实施例的系统、方法以及计算机程序产品。处理器238通过地址/数据总线348与存储器236通信。处理器238可以是任何市售或定制的微处理器。存储器236代表存储器装置的整体结构，它含有用于实现数据处理系统230的功能的软件和数据。存储器236可以包括(但不限于)以下类型的装置：缓存、ROM、PROM、EPROM、EEPROM、闪存、SRAM以及DRAM。

如图3所示，存储器236可以包括在数据处理系统230中使用的数种类别的软件和/或数据：操作系统352、应用程序354、输入/输出(I/O)装置驱动器358、以及数据356。如本专业的技术人员会理解的，操作系统352可以是适用于数据处理系统的任何操作系统，例如International Business Machines Corporation，Armonk，NY提供的OS/2、AIX或System390，或Mictrosoft Corporation，Redmond，WA提供的Windows95、Windows98、Windows2000、WindowsNT或WindowsXP，Unix或Linux。可以将操作系统配置成支持基于TCP/IP或其它网络通信协议的连接。I/O装置驱动器358通常包括软件程序，由应用程序354通过操作系统352存取，以便与诸如I/O数据端口246和某些存储器236部件等装置通信。应用程序354说明实现数据处理系统230各种特性的程序，并且最好包括能支持按照本发明实施例的操作的至少一种应用程序。最后，数据356代表由应用程序354、操作系统352、I/O装置驱动器358以及驻留在存储器236中的其它软件程序所使用的静态和动态数据。

又如图3所示，应用程序354可以包括强度分析应用程序360。强度分析应用程序360可以进行本文所述的用于评估影像的操作，以便检测可能与整体心脏损害相关联的强度改变。存储器236的数据部分356，如图3的实施例所示，可以包括影像数据362，例如MRI影像数据，它包括供作比较的所关注区域组织的第一和第二影像。

虽然本发明是例如将强度分析应用程序360作为图3中的应用程序进行说明的，但是如本专业的技术人员会理解的，其它配置也可使用而仍然可从本发明的内容获益。例如，可以将强度分析应用程序360结合到操作系统352、I/O装置驱动器358或数据处理系统230的其它这种逻辑分支中。因此，本发明不应被认为限于图3的配置，而应涵盖能够进行本文所述操作的任何配置。

图4A示出按照本发明特定实施例的操作。如图4A所示，获得病人组织的所关注区域的第一影像(框400)。影像的获得可以例如通过从成像系统，例如上述成像系统来采集影像，和/或从影像数据的数据库、文件或其它存储装置中获得影像。例如，病人的影像可以作为用于比较的第一影像被保持在历史数据库中供随后再调用。成像的病人组织的所关注区域可以例如包括心脏、血液、肌肉、脑、神经、骨骼、骨骼肌、肝、肾、肺、胰腺、内分泌、胃肠和/或泌尿生殖器组织。在本发明的特定实施例中，组织可以是人体组织。在其它实施例中，组织可以是动物的组织。

又如图4A所示，在一段时间之后，例如数小时、数日、数周、数月或甚至数年之后，获得所关注区域中组织的第二影像(框402)，供与第一影像作比较。供比较的第二影像反映所关注区域中组织特性的任何改变。第二个作比较的影像可以利用相应的第一影像来采集和登记(在相同的断层位置拍摄)。也可以参阅第一影像如上面所述那样获得第二影像。因此，例如，比较影像可以是历史影像，也可以是最近采集的影像。

对第一影像和第二影像进行评估，以便确定影像的一个或多个整体特性(框404)。影像的整体特性例如可以是所关注区域中像素/体素的平均强度。整体特性也可以是所关注区域中像素/体素的统计分析。例如，可以确定所关注区域中像素/体素的标准偏差、均值、或其它统计分析。所关注区域中像素/体素的特性直方图也可作为整体特性提供。为提供整体特性而评估的像素/体素的特性可以包括强度、颜色、饱和度和/或各个像素/体素的其它特性，以及多个像素/体素的相对特性，例如对比度比例等。

将所述评估结果提供给用户，或者可以为进一步分析而提供所述评估结果(框406)。例如可以对第一影像和第二影像进行比较，并可以将平均强度之差作为结果提供给用户。而且，可以确定特性的直方图和/或基线影像和比较影像之间的特性之差并将其作为结果来提供。此外，直方图可以与代表特定损害、疾病和/或状态的直方图轮廓库作图形匹配。例如可以把确定的结果以图形用户界面的一部分的形式提供。

所关注区域中的组织影像整体特性的评估结果可用于对例如因组织损害或上述其它情况导致的组织特性的改变的检测(也许是早期检测)。这种整体特性评估可以适用于检测导致所关注区域中不同组织特性的随机图形的组织特性或者是以不能用于检测组织特性的图形的分辨率成像的组织特性。

图4B示出利用注射造影剂的按照本发明特定实施例的操作。如图4B所示，获得病人组织所关注区域的基线影像(框450)。影像的获得可以例如通过从成像系统，例如图1所示的MRI系统来采集影像，和/或从影像数据的数据库、文件或其它存储装置中获得影像。例如，病人的影像可以被保持在历史数据库中，供随后调用作为用于比较的基线影像。基线影像可以是没有注射造影剂、注射造影剂后、和/或注射造影剂后一段时间，例如20分钟之后所拍摄的影像。所成像的病人组织的所关注区域可以例如包括心脏、血液、肌肉、脑、神经、骨骼、骨骼肌、肝、肾、肺、胰腺、内分泌、胃肠和/或泌尿生殖器组织。在本发明的特定实施例中，组织可以是人体组织。在其它实施例中，组织可以是动物的组织。

又如图4B所示，在注射造影剂后，获得所关注区域中组织的影像(框452)，供与基线影像作比较。供比较的影像反映了造影剂对所关注区域中组织的影响。在本发明的特定实施例中，所述影像可以是心肌延迟增强(MDE)影像。比较影像可以用相应的基线影像来采集和登记(在相同的断层位置拍摄)。比较影像也可如上述参阅基线影像来获得。因此，例如，比较影像可以是历史影像，也可以是最近采集的影像。

对基线影像和比较影像进行评估，以便确定影像的一个或多个整体特性(框454)。影像的整体特性例如可以是所关注区域中像素/体素的平均强度。整体特性也可以是所关注区域中像素/体素的统计分析。例如，可以确定所关注区域中像素/体素的标准偏差、均值、或其它统计分析。所关注区域中像素/体素的特性直方图也可作为整体特性提供。为提供整体特性而评估的像素/体素的特性可以包括强度、颜色、饱和度和/或各个像素/体素的其它特性，以及多个像素/体素的相对特性，例如对比度比例等。

将所述评估的结果提供给用户，或者可以为进一步分析而提供所述评估结果(框456)。例如可以将基线影像和比较影像进行比较，并且可以将平均强度之差作为结果提供给用户。而且，可以确定特性的直方图和/或基线和比较影像之间的特性之差并将其作为结果来提供。此外，将直方图与代表特定损害、疾病和/或状态的直方图轮廓库作图形匹配。例如可以把确定的结果作为图形用户界面的一部分来提供。

所关注区域中的组织影像整体特性的评估结果可用于检测组织的损害，也许是早期检测组织的损害。这种检测可以为导致在受损和健康组织中存在不同浓度的造影剂的损伤。这种整体特性评估可以适合于检测导致所关注区域中受损组织的随机图形的损害和/或以不能用于检测受损组织的图形的分辨率成像的损害。因此，例如，在1.5Tesla MRI成像系统的情况下，典型的心肌梗塞就不会被认为是整体影像，并且对梗塞位置上影像中增加的强度的检测结果及其位置不会被认为是受损组织的随机图形或不能以MRI成像系统的分辨率检测的受损组织的图形。

图5示出按照本发明特定实施例的操作。如图5所示，对病人注射造影剂(框400)，并采集病人心脏至少一部分的影像(框402)。在本发明的特定实施例中，所采集的灌注影像可以是心肌延迟增强(MDE)影像。在MDE中，在注射造影剂例如钆DPTA后20分钟，一些造影剂已渗入到坏死(死亡)组织中，并会呈现为很亮(故为延迟增强)。这些影像可以用相应的基线灌注影像来采集和登记(在相同的断层位置拍摄)。

对所采集的影像进行评估，并将影像的强度与基线影像进行比较(框404)。基线影像是病人心脏的影像，可以是以前采集的、也是在注射造影剂后采集的影像。基线影像可以是在施行治疗方案之前已采集，或是在较早评估时所采集的影像。影像的比较可以是影像平均强度的比较，见以下详述。如果与基线影像相比影像的强度没有增加(框406)，则提示不存在整体心脏损害(框408)。如果与基线影像相比影像的强度已增加(框406)，则提示可能存在整体心脏损害(框410)。

在本发明又一些实施例中，可以利用影像处理技术自动地或部分自动地进行整体影像特性，例如心脏影像强度的评估。自动比较例如还可以包括登记相互有差异的影像。可以利用传统图形识别和/或对准技术来提供这种登记，以便使影像的相应像素或影像的各部分各自与病人体内大致相同的物理位置相关联。

在本发明的特定实施例中，可以将病人带到MRI室中，使他们仰卧在MRI台上，加上ECG导线和呼吸选通风箱(gating bellows)。例如可在1.5 Tesla GE CV_i扫描器上进行MRI扫描，在胸部周围加上相控阵表面线圈，以优化信噪或其它MRI扫描器。可以利用梯度回声技术采集影像，重复时间(TR)和回声时间(TE)均基于对象的R-R时间间隔。多断层冠状梯度回声序列可用于获得胸部的搜索影像并定位左心室。为对象静脉注射钆造影剂(0.2mmole/kg钆(Prohance，BraccoDiagnostics，Princeton，NJ)。可以将注射时间记录下来。

注射造影剂时间后大约20分钟，先加非选择性饱和脉冲再用快速梯度回声采集三个短轴视图(基底部，中间部，和心尖部)的延迟增强影像。可以区划用于采集这些影像的界标，从房室沟内的冠状窦水平延伸到横过二尖瓣环。可以使用38cm视野，每段24个视图，8mm断层厚度，2NEX，256×256成像矩阵，以及0.75矩形视野来采集这些影像。可以从140msec到160msec调节用于延迟增强影像的倒序时间(TI)，以便提供均匀的暗背景。此外，在这三个短轴断层位置，快速-梯度-重新调用的回声脉冲序列可以与相位编码排序一起使用。可以对这些影像进行相位敏感重构，以减少当TI改变时在重要视图中观察到的视在对比强度的变化。此外，相位敏感重构可以降低对随着钆造影剂注射后增加的延迟而在组织T₁中的改变的敏感度。

影像采集完成后，其位置、测量值和代表性影像可以以电子方式传送到数据库。这种信息可供MRI技术人员在每次扫描时通过PC工作站使用，便于在随后的研究中重新定位(登录)断层位置。

图12示出用于计划影像断层的软件的屏幕捕获。这种软件可提供影像计划断层以及定位坐标的电子复制件，它们均被保存，供研究中随后访问时检索。这对提高MRI技术人员从以前的访问中再现断层位置的能力很有效果。在图12所示实例中，示出心脏的长轴视图，以及结果延迟增强短轴视图。

在延迟增强采集时，可以在全部多断层采集时确定包括LV心肌的所关注区域(ROI)。可以避免与LV腔内血池相关联的高信号强度。可以从基线和延迟增强的影像中记录ROI中每个(或所选)体素的信号强度和位置(x，y和z坐标)。也可以通过使用体外空气/空间中的单独ROI，从强度中减去例如没有造影剂的单独ROI的平均强度来导出数值。ROI可以用来确定在两个影像之间强度的改变，如在以下”实例”中所述。

虽然以上就心脏的特定视图、区域、面积和/或断层对本发明实施例作了说明，但也可使用心脏的其它视图、区域、面积和/或断层。而且，可以使用少于或多于三个断层。此外，可以沿心脏的长轴或短轴拍摄影像。因此，本发明的某些实施例不应被认为限于心脏的特定视图，而是可以包括允许进行强度分析的心脏的任何视图和/或视图数量，以检测整体心脏损害。

通常，在治疗前或治疗早期，或作为诊断心脏情况改变的初始参考点，来获得第一基线影像。供比较的后继影像可以每天、每周、或按照其它固定的或可变的时间间隔拍摄，或在计划治疗，如细胞毒素治疗前或后采集。

现在，在以下非限制性实例中对本发明进行更详细的说明。

实例

如以上简述的，按照传统方法，通过在非增强LV心肌内确定具有高于背景强度＞2标准偏差的信号强度的体素的位置来识别患缺血性心肌病的病人的心肌细胞坏死。通过确定超增强的透壁程度(表示为从心内到心外膜表面线性延伸的高强度像素的数量相对于从心内膜到心外膜的总距离的比例)来量化坏死的数量。由于在冠状动脉血流减少的背景下，心肌坏死在从心内到心外膜表面的波阵面中继续进行，所述方法对于心肌梗塞后评估坏死数量是很有用的。

但是，所述方法可能并不同样适用于会引起在随机分布图形中整个LV心肌敏感组织坏死(例如，整体损害)的过程。为克服这种限制，可对LV中三个短轴断层位置(心尖部，中间部和基底部)内的体素，且在一些实施例中为全部体素，进行抽样，并且识别3维空间中每个体素的强度以及x、y和z坐标(图6))。图6为左心室3个短轴平面(基底部，中间部，和心尖部)的三维图示。在每个平面中，每个断层面上的小方形网格分界体素。分析时，记录每个体素的影像强度以及x，y和z坐标。这样，用延迟增强技术识别的、与引起心肌细胞坏死的随机分布过程关联的高强度像素就可特显出来(影像上的白色斑点)。

对影像中体素强度变化的校正也可通过确定目标区域中，通常是心脏外1cm直径的圆形所关注区域(ROI)中的体素强度来识别。对于每个心尖部、中间部和基底部断层，可以确定在给定强度下的像素数，并将心脏外的ROI的强度从像素中减去。在某些实施例中，对于每个断层，全部体素的平均强度以及在最高40％的分布中的峰值体素强度都可被确定(图6)。这样，就可识别高强度像素区域相对于它们在左心室中的位置。

图7为在LV中间短轴视图中的示范性延迟增强MR影像(上部各小图)。心肌为灰色，血池为白色。造影剂注射后20分钟ROI(红线)内体素的数量(y轴)和强度(x轴)示于下部各小图。造影剂被所有心肌细胞吸收，但注射后20分钟，它没有从坏死细胞中被清除。如图中所示，在健康正常病人中造影剂摄取的平均强度很低(最左图)，在患梗塞的病人中最高(从左数第三图)。中间平均强度显示在与阿霉素心肌病人关联的直方图上(从左数第二图)。

为确定静脉注射后20分钟MRI评估钆造影剂摄取的位置和数量的效用，对进行了年龄(从35到50岁)和性别匹配的4组参与者进行了截面研究。这些包括：

a)(组I)：4个对象(1男，3女)，没有疾病，不用心脏药物，且通过MRI具有正常的心脏收缩和舒张功能，

b)(组II)：3个病人(3女)，在造影剂冠状血管造影上没有冠状动脉腔狭窄，但具有较差的LV射血分数(＜35％)和继发于阿霉素注射的充血性心力衰竭，

c)(组III)：3个病人(2男，1女)，在造影剂冠状血管造影上没有冠状动脉腔狭窄，但具有较差的LV射血分数(＜35％)和继发于原发性扩张性心肌病的充血性心力衰竭，

d)(组IV)：3个病人(2男，1女)，有继发于缺血性心肌病的LV功能紊乱以及ST前段抬高的心肌梗塞。

每组中一个对象的中间短轴影像和影像中体素的强度分布示于图7，所有参与者的全部断层中体素强度的分布示于图8。

图8中示出在注射造影剂后20分钟对象的截面抽样中全部参与者的ROI中体素的百分比(y轴)和强度(x轴)。如图7所示，与正常年龄匹配控制组相比，在因化疗用药而患心肌病的病人中显示有15到30范围的强度增加百分比。所述强度图形看来不同于在患缺血性心肌病的病人中所见的图形。

为确定在左心室的每个断层中高强度像素图形之间的关系，使用了自动校正统计。系列自动校正量度(I)定义如下。设δ_ij为像素i和j之间距离的加权函数，n为像素数，且x_i为第i个像素的强度。则定义：

$I=n\frac{\underset{\mathrm{ij}}{\mathrm{\Σ}}{\mathrm{\δ}}_{\mathrm{ij}}(x_i-\stackrel{\‾}{x})(x_j-\stackrel{\‾}{x})}{\left(\underset{\mathrm{ij}}{\mathrm{\Σ}}{\mathrm{\δ}}_{\mathrm{ij}}\right)\left(\underset{i}{\mathrm{\Σ}}{(x_i-\stackrel{\‾}{x})}^2\right)}.$

I是系列自动校正的量度，当邻近像素都高于或低于均值时所述量度较高(Ripley，1981)。实际上，一直使用表达式

式中d(x_i，x_j)是点x_i和x_j之间的欧几里得距离。

使用这种分析形式，高数量表示在ROI内的图形群集，而低数量更能代表随机关联性。如图9所示，与MI关联的增高信号强度紧密群集在梗塞区域；而与阿霉素毒性关联的那些信号强度则分散在整个LV。患继发于阿霉素注射的心肌病的病人其造影剂摄取的图形是随机的，且显著不同于(p＜0.001)与继发于心肌梗塞的心肌细胞坏死相关联的高信号强度体素的图形。

为了确定在接受化疗的病人中造影剂增强是否与LVEF的跌落相关联，在开始化疗前对病人进行基线MRI检查，然后按照研究计划进行附加MRI检查。还进行了超声波心动描记检查，以监控在各次MRI检查之间病人的左心室功能。一个对象发生了呼吸困难并作了超声波心动描记检查以便确定LVEF。所述对象的LVEF从55％跌落到48％。所述病人进行了MRI测试和影像分析。将所述对象的影像分析与在化疗过程中没有发生LVEF跌落的另一对象作比较。病人的中间短轴视图中的影像和体素强度示于图10。

图10示出在两个单独的病人接受化疗时在两个时间点的影像和平均体素强度，其中一人在化疗过程中发生呼吸困难。两个病人治疗前的影像示于左边，治疗后的影像示于右边。影像中ROI的平均体素强度示于影像的下面。发生LVEF跌落的病人1(上图)，在接受用于治疗乳腺癌的400mg/m²的蒽环霉素后的第二次检查中发生增高的造影剂摄取和信号强度。第二个病人(下图)，没有发生LVEF跌落且摄取图形显示没有显著改变。如图所示，有LVEF跌落的病人，与第一次检查相比，第二次检查中LV内的体素强度有显著增加，而没有LVEF跌落的病人，第二次检查没有明显的改变。

为了确定在医疗情况没有实质改变的情况下参与者的MRI延迟增强体素强度随时间的变异，对4个人在注射造影剂后两周期间作了两次研究。一个参与者的影像示于图11，全部4个人的两个抽样点的数据示于表1。

表1：在4个人中，MRI强度(均值±标准偏差)和LVEF。

图11示出在病人情况没有改变的情况下隔21天所采集的中间部左心室短轴视图。应当指出，使用本文其它地方所述的软件在第二次采集影像上断层位置近乎精确地复制。造影剂注射后20分钟，ROI内的信号强度没有显著不同，即5.8对6.1(p＝NS)。用这种技术的MRI检查可随时间的推移作再现性采集。

在第一次和第二次检查之间对象的造影剂摄取图形的改变很少，且对于在两个时间点测量的4个人，2次测量之间的校正极好(y＝0.87x+1.2，R²＝0.96)。

基于上述数据，与年龄和性别匹配的控制对象相比，显然患继发于化疗引发心肌毒性的心肌病的病人其造影剂的延迟增强MRI吸收图形被升高。这种造影剂摄取图形是扩散的并随机分布在整个左心室中，其形式明显不同于对持续心肌梗塞病人所观察的心肌细胞损害。在涉及两个接受化疗的病人的项目中，一个病人发生增高的造影剂摄取，与LVEF的跌落一致，但未发生LVEF跌落的另一个人就不是这样。这种方法论和分析方法可高度再现，并呈现很低的观察者内在的可变性。

以上所述是说明本发明，而不应被认为是限制本发明。虽然已对本发明的几个示范性实施例作了说明，但是，本专业的技术人员会很容易理解在实质上不背离本发明的新颖内容和优点的前提下可以作许多改动。因此，所有这些改动应被包括在权利要求书中所定义的本发明的范围之内。

Claims (9)

1.一种用于检测病人整体心脏损害的系统，包括：

用于在所述病人注射造影剂后获得第一心脏影像的装置；

用于在所述病人注射造影剂后获得第二心脏影像的装置；以及

用于将所述第一心脏影像的、所关注区域的像素/体素特性直方图的倾斜、峰态或标准偏差的至少其中之一和/或与所关注区域相关联的至少一个像素/体素强度直方图的形状或分布的至少其中之一和所述第二心脏影像相应的、所关注区域的像素/体素特性直方图的倾斜、峰态或像素/体素特性的标准偏差的至少其中之一和/或与所关注区域相关联的至少一个像素/体素强度直方图的形状或分布的至少其中之一进行比较，以便确定整体心脏损害的潜在可能性的装置。

2.如权利要求1所述的系统，其中所述第一心脏影像和所述第二心脏影像是磁共振成像MRI影像。

3.如权利要求1所述的系统，其中用于比较的装置被配置成比较多个像素/体素直方图的形状或分布。

4.如权利要求1所述的系统，其中还包括：

用于获得所述第一心脏影像的所述心脏外所关注区域的影像数据的装置；

用于利用所述心脏外所关注区域的第一影像的数据对所述第一心脏影像中像素/体素强度的变化进行校正的装置；

用于获得所述第二心脏影像的所述心脏外所关注区域的影像数据的装置；

用于利用所述心脏外所关注区域的第二影像的数据对所述第二心脏影像中像素/体素强度的变化进行校正的装置；以及

用于利用校正的第一心脏影像和校正的第二心脏影像获得第一心脏影像的强度和第二心脏影像的强度的装置。

5.一种用于评估病人组织特性的系统，包括：

用于在第一次采集期间获得包括所关注区域的组织的第一影像的装置；

用于在第二次随继的采集期间获得包括所关注区域的组织的第二影像的装置；以及

用于确定所述第一影像和所述第二影像的所述所关注区域的整体强度特性，以便能比较所述第一影像和所述第二影像的所述整体特性，以便确定所述所关注区域的所述组织的整体损害的潜在可能性的装置，其中所述整体强度特性包括所关注区域中像素/体素特性的标准偏差、像素/体素特性直方图的倾斜或峰态的至少其中之一和/或至少一个像素/体素强度直方图的形状或分布的至少其中之一。

6.如权利要求5所述的系统，其中所述整体强度特性包括所述所关注区域的基于所述所关注区域中基本上全部像素/体素的像素/体素特性。

7.如权利要求5所述的系统，其中用于确定所关注区域的整体强度特性的装置包括用于评估所关注区域中像素/体素特性的标准偏差、像素/体素特性直方图的倾斜和峰态的至少其中两个和/或与所述第一影像和第二影像相关联的多个像素/体素直方图的形状或分布。

8.如权利要求5所述的系统，其中所述所关注区域中的所述组织包括以下各项中的至少一项，即：心脏、血液、肌肉、脑、神经、骨骼、骨骼肌、肝、肾、肺、胰腺、内分泌、胃肠和/或泌尿生殖器组织。

9.如权利要求5所述的系统，其中所述第一影像和所述第二影像包括磁共振成像MRI影像、X射线计算机体层摄影CT影像、超声波影像、单光子发射计算机体层摄影SPECT影像和/或正电子发射体层摄影PET影像的至少其中之一。

Images