专利名称：医学数据的时间配准的制作方法
技术领域：
本发明涉及计算机辅助诊断，其更具体但非排他地涉及用于(使用运动特征的)
医学图像数据的时间配准的装置和方法。
背景技术：
在工业化国家中，心血管疾病是一个非常主要的致死原因。其早期诊断和治疗对 于降低死亡率和改善患者的生活质量而言是极其重要的。医学成像和计算机辅助诊断通过 提供例如关于患者内部器官的有用信息在辅助医生、放射学家方面发挥着越来越重要的作 用。当前，诸如计算机断层摄影(CT)或磁共振成像(MR)等的无创医学成像程序不仅允许 描述例如心脏解剖结构的3D图像的数据采集，还允许描述心脏解剖结构和功能的4D图像 序列的数据采集(即，还包含时间分量)。 为了实现相同或不同患者的不同医学图像数据集之间的客观比较，需要图像数据 配准。同样，为了使用来自许多不同患者的医学图像建立有效的平均心脏运动模型，不仅必 须在空间上，而且必须在时间上对数据进行配准。例如，为了使独立跳动的各心脏(每个都 在其本身的患者坐标系中)在空间上彼此共配准，可以执行Procrustes分析，其将每个患 者坐标系变换到公共模型坐标系中，并且可以通过针对样本中所有患者在给定心脏相点上 对心脏表面上的可识别标志物位置的运动进行平均而在模型空间中计算平均模型。
通常，将心动周期描述为起始于心电图(ECG)的R峰并具有1/r的绝对持续时间 的100%完整周期，其中，r为心率。因此，每个相点具有位于周期的0%到100%之间的专 有时间位置。然而，这并不一定意味着表示诸如末期收縮(收縮期的末期)等的可识别事 件的点在心动周期内具有这种专有位置。这种生理相点的时间失准的原因可能是由于采集 参数的不同(例如，触发点从R峰偏移以及在连续帧的采集中的不同间隔)、心动周期的长 度的不同或心脏动态特性的不同。例如，一位患者的心脏可能较另一位患者的心脏具有较 长的收縮期而具有较短的舒张期。还公知的是，当心率上升时，收縮期的持续时间可能不如 舒张期的持续时间降低得大，这可能是可识别事件并非线性地对应于心动周期的R-R间隔 中的点的原因。 图1示出了表示四位不同患者的不同相点的标志物(即，心脏上的可识别点)的 平均位移的运动特征的示例。使用由Vembar等人在Med. Phys. 200330 (7)p. 1683ff中公开 的方法，针对不同心率的影响来补偿运动特征。对四位不同患者中的每一位来说，舒张期的 峰和收縮期的峰清晰可见。尽管对象之间的舒张峰的时间位置很好地吻合在经校正的R-R 间隔的约20%处，但是对象之间的收縮峰的时间位置显著变化。 图1的示例示出了可识别生理事件的相应相点的时间对准是无法保证的，当将两 个独立的心脏运动模式定量地彼此比较时，这可能是严重的问题。
发明内容
本发明的优选实施例试图克服现有技术的以上缺点中的一个或多个。
根据本发明的一方面，提供了一种供计算机系统使用以对时变医学数据进行比较 的数据结构，所述数据结构包括可执行以完成以下步骤的计算机代码 接收至少一个第一数据集，其包括表示多个第一时间处的医学参数的第一数据；
接收至少一个第二数据集，其包括表示多个第二时间处的所述医学参数的第二数 据；以及 对所述第一和/或第二数据中的至少一些进行处理以增加表示相应的多个可识 别事件的多个所述第一和第二时间之间的相关程度。 通过处理表示单一医学参数的第一和第二数据，这提供的优势为降低了所需的处 理能力，使增加多个所述第一和第二时间之间的相关程度的过程显著更快。例如，根据每 位患者的预先确定的标志物生成心脏的运动特征，并且彼此配准，以便估计表示可识别事 件的不同特征中的数据之间的时间偏差，因此允许不同患者之间的医学图像数据的时间对 准。因此，配准是一维过程，其更加易于执行并且因而更快。本领域技术人员将会意识到在 本发明的背景下，"相关"意味着以下两者之间的相似程度(i)多个第一时间，其表示一个 数据集中的一序列可识别或可检测事件，和(ii)多个第二时间，其表示另一个数据集中的 对应的事件序列。 所述医学参数涉及患者体内的至少一个预先确定的可识别位置。 这提供的优势是对医学数据的分割已经存在，因此简化了所述第一和第二数据集
的生成。 所述医学参数可表示多个可识别位置的平均位移值。 这提供的优势是诸如收縮或舒张等的重要生理事件在所述第一和第二数据集内 是可识别的，因此允许不同数据集之间的时间配准。 通过使所述第一和第二数据集之间的全局相似性度量最大化可以增加多个所述 第一和第二时间之间的相关程度。 全局相似性度量可包括所述第一和/或第二数据点的互相关和/或距离平方和。
处理所述第一和/或第二数据的步骤可包括调整至少一个第一和/或第二时间以 增加所述相关程度。 所述计算机代码是可执行的以限制对至少一个所述第一和/或第二时间进行调 整的量。 同样，所述第一和第二数据集为心脏特征数据。 根据本发明的另一方面，提供了一种在其上存储有如上所定义的数据结构的计算 机可读介质。 根据本发明的又一方面，提供了一种用于处理时变医学数据的医学数据处理装
置，所述装置包括至少一个适于处理如上所定义的数据结构的处理器。 根据本发明的又一方面，提供了一种医学成像装置，其包括 至少一个成像设备，其用于形成医学数据； 如上所定义的医学数据处理装置；以及 至少一个显示设备，其用于在处理所述第一和/或第二数据之后显示所述第一和 第二数据集的表示。 根据本发明又一方面，提供了一种比较时变医学数据的方法，所述方法包括
接收至少一个第一数据集，其包括表示多个第一时间处的医学参数的第一数据；
接收至少一个第二数据集，其包括表示多个第二时间处的所述医学参数的第二数 据；以及 对所述第一和/或第二数据集中的至少一些进行处理以增加表示相应的多个可 识别事件的多个所述第一和第二时间之间的相关程度。
通过参考附图，仅以示例且不具有任何限制意义的方式描述本发明的优选实施 例，其中 图1示出了来自4位不同患者的运动特征的示例，其中纵轴表示以[mm]为单位的 平均顶点位移，横轴表示以[% ]为单位的R-R峰之间的心动周期；
图2是体现本发明的医学成像数据处理装置的各部件的示意性表示；以及
图3是示出体现本发明的处理医学图像数据的方法的流程图。
具体实施例方式
参考图2，医学成像数据处理装置具有计算机断层摄影(CT)成像装置l，其包括围 绕圆形框架4相对地成对布置的X射线源2和探测器3。处理器5a将患者6的CT图像数 据集处理成4D图像序列7a，并且将其与来自不同时间的相同患者、不同患者或代表性参考 平均模型的参考图像数据7b进行比较。处理器5a生成3D图像8a和8b，其包括序列中的 每个图像中的标志物，并且使不同数据集7a、7b在空间上彼此共配准。另外，处理器5a根 据3D图像序列8a和8b生成运动特征9a和9b，并且执行程序10，程序10产生空间和/或 时间上对准的特征，其被显示在显示设备11上，从而使得医师或放射学家可直接地比较运 动特征和/或4D图像数据序列。 参考图3，处理器5a在步骤S10获得患者6的图像数据并且获得参考图像数据。 之后处理器5a在步骤S20生成来自患者6和参考数据的4D数据集。在步骤S30，使用例 如形状跟踪方法识别标志物，并且针对患者6和参考数据产生包括标志物的3D图像数据序 列8a、8b。在步骤S40，处理器5a使3D患者图像数据序列8a与3D参考图像数据序列8b 空间上对准，并且在步骤S50为这两个3D图像数据序列提供运动特征。运动特征表示所有 被识别标志物的平均位移幅度。在步骤S60，处理器5a对所述患者6和参考图像数据7a、 7b的运动特征9a和9b在时间上进行对准。为了使运动特征彼此对准，可在时域上使用 Procrustes分析。 程序10允许运动特征9a、9b中的一个或两个的给定相点的位移，以便使应用例如 互相关和/或距离平方和的运动特征9a、9b之间的全局相似性度量最大化。程序10也可 应用正则化约束，以避免运动特征9a、9b的相点的过度局部位移。
在步骤S70，在显示单元11上显示由处理器5a产生的输出。 上述方法尤其允许相对于R-R周期中的生理相点而非等距相点或仅使用心率进 行补偿的那些相点来比较两个不同心脏的跳动模式。这改善了两个跳动心脏的时间特性的 可比较性，而且还允许在时域中建立更准确的模型，其可得到具有更高预测值的模型。
除了运动特征，可以使用反映显著相点的任何其他总运动模式特征。用于心脏图像而言，可以代替地使用左心室体积曲线，这是因为其指示左心室的收縮和舒张。 相同的方法也可应用于其他周期性运动的4D图像数据集，诸如呼吸运动以及复
杂的非周期性运动，例如在不同时间的利用不同肌肉的关节弯曲。另外，除了以上使用的单
值运动特征，也可以使用多值运动特征。 本领域技术人员将会意识到仅以示例且不具有任何限制意义的方式描述了上述 实施例，并且可能存在各种变化和修改而不偏离由随附的权利要求书所限定的本发明的范 围。
权利要求
一种供计算机系统使用以对时变医学数据进行比较的数据结构，所述数据包括可执行以完成以下步骤的计算机代码接收至少一个第一数据集，其包括表示多个第一时间处的医学参数的第一数据；接收至少一个第二数据集，其包括表示多个第二时间处的所述医学参数的第二数据；以及对所述第一和/或第二数据中的至少一些进行处理以增加表示相应的多个可识别事件的多个所述第一和第二时间之间的相关程度。
2. 根据权利要求1所述的数据结构，其中，所述医学参数涉及患者体内的至少一个可 识别位置。
3. 根据权利要求2所述的数据结构，其中，所述医学参数表示多个所述可识别位置的 平均位移值。
4. 根据权利要求1所述的数据结构，其中，通过使所述第一和第二数据集之间的全局 相似性度量最大化而增加多个所述第一和第二时间之间的所述相关程度。
5. 根据权利要求4所述的数据结构，其中，所述全局相似性度量包括表示所述第一和 第二数据的点之间的互相关和/或距离平方和。
6. 根据权利要求1所述的数据结构，其中，处理所述第一和/或第二数据的步骤包括调 整至少一个第一和/或第二时间以增加所述相关程度。
7. 根据权利要求6所述的数据结构，其中，所述计算机代码是可执行的以限制对至少 一个所述第一和/或第二时间进行调整的量。
8. 根据权利要求1所述的数据结构，其中，所述第一和第二数据集为心脏特征数据。
9. 一种计算机可读介质，在其上存储有根据权利要求1所述的数据结构。
10. —种用于处理时变医学数据的医学数据处理装置，所述装置包括至少一个适于处 理如权利要求1所述的数据结构的处理器。
11. 一种医学成像装置，包括 至少一个成像设备，其用于形成医学数据； 根据权利要求10所述的医学数据处理装置，以及至少一个显示设备，其用于在处理所述第一和/或第二数据后显示所述第一和第二数 据集的表示。
12. —种比较时变医学数据的方法，所述方法包括接收至少一个第一数据集，其包括表示多个第一时间处的医学参数的第一数据； 接收至少一个第二数据集，其包括表示多个第二时间处的所述医学参数的第二数据；以及对所述第一和/或第二数据中的至少一些进行处理以增加表示相应的多个可识别事 件的多个所述第一和第二时间之间的相关程度。
全文摘要
公开了一种供计算机系统使用以对时变医学数据(9a，9b)进行比较的数据结构。该数据结构执行以下步骤接收第一数据集(9a)，其包括表示多个第一时间处的医学参数的第一数据；接收第二数据集(9b)，其包括表示多个第二时间处的所述医学参数的第二数据；以及对所述第一和/或第二数据集进行处理以增加表示可识别事件的多个所述第一和第二时间之间的相关或相似程度。
文档编号G06T7/00GK101765863SQ200780100843
公开日2010年6月30日 申请日期2007年12月12日 优先权日2006年12月19日
发明者C·洛伦茨, J·冯贝格 申请人:皇家飞利浦电子股份有限公司