CN107004270A - 用于计算感兴趣对象的位移的方法和系统 - Google Patents

用于计算感兴趣对象的位移的方法和系统 Download PDF

Info

Publication number
CN107004270A
CN107004270A CN201580068805.2A CN201580068805A CN107004270A CN 107004270 A CN107004270 A CN 107004270A CN 201580068805 A CN201580068805 A CN 201580068805A CN 107004270 A CN107004270 A CN 107004270A
Authority
CN
China
Prior art keywords
image
subset
interest
displacement
freshly harvested
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201580068805.2A
Other languages
English (en)
Other versions
CN107004270B (zh
Inventor
徐泾平
K·卢
W·T·史
C·S·霍尔
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN107004270A publication Critical patent/CN107004270A/zh
Application granted granted Critical
Publication of CN107004270B publication Critical patent/CN107004270B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/20—Analysis of motion
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/20—Analysis of motion
    • G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
    • G06T7/248—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments involving reference images or patches
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/38—Registration of image sequences
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10016—Video; Image sequence
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10132—Ultrasound image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20081—Training; Learning
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

本发明涉及一种计算感兴趣对象的位移的方法，包括根据所述感兴趣对象的预先采集的图像的集合中相邻图像来计算(101)所述感兴趣对象的位移模型的步骤，所述位移模型反映所述感兴趣对象沿时间的位置。所述方法特征在于所述方法还包括以下步骤。基于所述位移模型来确定(102)在所述预先采集的图像的集合的一个周期性时间循环内的来自所述预先采集的图像的集合图像的第一子集(S1)的步骤。根据新采集的图像来识别(103)图像的第二子集(S2)的第一步骤，其中，所述图像的第二子集(S2)中的图像是相继的并且在所述图像的第一子集(S1)中具有相同的最相似图像，其中，通过比较所述新采集的图像中的给定图像与所述图像的第一子集(S1)中的每幅图像来确定相似性水平的第一集合，并且其中，所述最相似图像具有在所述相似性水平的第一集合中的最大相似性水平。在所述图像的第二子集(S2)中选择(104)给定图像作为第一参考图像(I1)的第一步骤。从所述新采集的图像识别(105)图像的第三子集(S3)的第二步骤，其中，所述图像的第三子集(S3)中的图像是相继的并且在所述图像的第一子集(S1)中具有相同的最相似图像，其中，通过比较所述新采集的图像中的给定图像与图像的所述第一子集(S1)中的每幅图像来确定相似性水平的集合，并且其中，所述最相似图像具有在所述相似性水平的集合中的最大相似性水平。在图像的所述第三子集(S3)中选择(106)给定图像作为第二参考图像(I2)的第二步骤。计算(107)所述第二参考图像(I2)与所述第一参考图像(I1)之间的位移的步骤。本发明还涉及位移计算的对应系统。

Description

用于计算感兴趣对象的位移的方法和系统
技术领域
本发明涉及对象移动测量的领域，尤其涉及用于计算感兴趣对象的位移的方法。
背景技术
辐射治疗广泛用在癌症处置或治疗中。辐射治疗的目的是杀死患者身体中的感兴趣对象中发生的病理结构，例如，肿瘤细胞，而不损害感兴趣对象的周围组织。这通过将较高辐射剂量递送到患者身体的感兴趣对象，并且将可能较低辐射剂量递送到感兴趣对象周围的健康器官组织来完成。然而，在整个处置过程期间，患者身体的呼吸运动可以产生关于辐射被应用的位置的不准确。实际上，患者身体的感兴趣对象的位置可以随着呼吸运动而改变，并且甚至患者身体的感兴趣对象的尺寸和形状有时随着通过内部身体肌肉力的变形而改变。
对象跟踪当感兴趣对象在患者身体中时由于身体移动(诸如，由呼吸引起的身体移动)而是具有挑战的。跟踪由呼吸引起的患者身体中的这样的感兴趣对象的移动中的困难可以是进行辐射治疗时的不准确性的重要原因。
超声成像被用于在处置期间提供图像。用于运动跟踪的已知方法是标记被放置在患者的皮肤上。标记然后使用诸如超声的成像设备来跟踪。然而，处理图像的大的集合在对象跟踪中是具有挑战的。当前，还没有可用于跟踪患者身体中的对象的位置的解决方案，因为需要高的计算复杂性。测量对应的表面移动还不足够准确以反映身体中的对象移动。
公开物，speeding-up image registration for repetitive motionscenarios,ISBI(1355-1358)，公开了一种用于针对经历呼吸的器官的图像序列的实时图像配准的方法。在训练阶段期间，图像被配准，并且然后通过采用对图像的维度降低和对应位移的存储来学习图像外观与空间变换之间的关系。对于应用阶段中的每幅图像，训练设置中最相似的图像被计算并且进行位移预测。公开物，respiration induced fiducialmotion tracking in ultrasound using an extended SFA approach,Proc.Of SPIEVol.9419，公开了一种用于实时将在训练阶段中从超声图像序列提取的预估计的目标运动转移到在线数据的方法。所述方法基于提取目标对象的特征点，通过慢特征分析利用特征运动的低维描述，并且找到来自训练数据的最相似的图像特征，用于估计当前或在线对象位置。对于每幅采集的图像，计算最相似的图像，并且基于针对每幅采集的图像的最相似的图像来进行映射。因此，计算复杂性是高的。
发明内容
本发明的目的是提出一种计算感兴趣对象的位移的方法。
独立权利要求限定了本发明。从属权利要求定义了有利实施例。
根据本发明的方法旨在计算感兴趣对象的位移。所述方法包括根据感兴趣对象的预先采集的图像的集合中相邻的图像来计算感兴趣对象的位移模型的步骤，所述位移模型反映感兴趣对象沿时间的位置。所述方法还包括：
-基于所述位移模型来确定所述预先采集的图像的集合的一个周期性时间循环内的来自所述预先采集的图像的集合的图像的第一子集的步骤；
-从新采集的图像识别图像的第二子集的第一步骤，图像的所述第二子集(S2)中的图像是相继的并且在图像的第一子集(S1)中具有相同的最相似图像，其中，通过将所述新采集的图像中的给定图像与图像的所述第一子集(S1)中的每幅图像进行比较来确定相似性水平的集合，并且其中，所述最相似图像具有在相似性水平的所述集合中的最大相似性水平；
-在图像的所述第二子集中选择给定图像作为第一参考图像的第一步骤；
-从所述新采集的图像识别图像的第三子集的第二步骤，其中，图像的所述第三子集(S3)中的图像是相继的并且在图像的所述第一子集(S1)中具有相同的最相似图像，其中，通过将所述新采集的图像中的给定图像与图像的所述第一子集(S1)中的每幅图像进行比较来确定相似性水平的集合，并且其中，所述最相似图像具有在相似性水平的所述集合中的最大相似性水平；
-在图像的所述第三子集中选择给定图像作为第二参考图像的第二步骤；
-计算所述第二参考图像与所述第一参考图像之间的位移的步骤。
利用该方法，能够利用有限的计算资源来计算由感兴趣对象的移动引起的位移。事实上，该方法降低了计算复杂性，这得到在图像采集正在进行时跟踪感兴趣对象的位置的现实的解决方案。其还改进检测身体中的对象的位置改变的准确性。
通过使用该方法，位移计算基于来自新采集的图像的两个相继子集的选定图像。相比于计算两幅相继的新采集的图像之间的位移，其降低了计算复杂性。
本发明还涉及一种用于计算位移的系统。
所述系统包括用于执行根据本发明的上述方法的各个步骤的处理单元。
以下将给出对本发明的详细解释和其他各方面。
附图说明
现在将参考下文描述并且结合附图考虑的实施例来解释本发明的特定方面，其中，相同部分或子步骤以相同方式来指定：
图1描绘了根据本发明的方法的流程图，
图2图示了根据本发明的感兴趣对象的位移模型的范例，
图3图示了根据本发明的图像的第一子集与预先采集的图像的关系，
图4图示了根据本发明的新采集的图像的相似性的识别，
图5图示了根据本发明的第一参考图像和第二参考图像的选择，
图6描绘了根据本发明的系统的示意图，
图7描绘了根据本发明的系统的示意图。
具体实施方式
例如，辐射治疗是在癌症处置或者更一般地在医学治疗中广泛使用的处置。辐射治疗旨在通过将较高辐射剂量递送到肿瘤并且将可能较低剂量递送到健康器官组织来杀死患者身体中的肿瘤细胞，而不损害周围组织。在下文中，接收辐射的患者身体的部分被描述为“感兴趣对象”。例如，在对腹部的肿瘤的处置期间，感兴趣对象可以具有由患者的呼吸引起的特定移动。通过计算由该移动引起的位移，位移信息被发送到用于处置的辐射治疗装备的控制模块。控制模块使用该位移信息来补偿患者身体的感兴趣对象的位移，以便使辐射射束跟随感兴趣对象，即使感兴趣对象具有特定移动。
在处置的开始之前，利用医学成像系统来扫描患者身体的感兴趣对象以采集图像，例如，超声成像、X射线、CT和MRI等。例如，在处置之前几小时或几分钟进行扫描。
在扫描期间，患者例如水平地躺在床上。成像扫描的持续时间等于患者的若干呼吸循环的持续时间。由医学成像系统生成的这些图像在本发明中被称为“预先采集的图像”。在图像扫描期间获得图像允许导出沿时间的感兴趣对象的位移模型。
根据第一方法，通过依次计算在两幅相继的预先采集的图像之间的感兴趣对象的位移来获得位移模型。通过随时间进行此，能够测量患者身体的感兴趣对象的整体患者位移模型。
根据第二方法，通过在扫描期间使用被附接在患者皮肤上的标记来获得所述位移模型。基于被提供到一些传感器的标记的位置信息，能够计算感兴趣对象的位移信息。因此，能够通过沿着时间进行此来测量患者身体的感兴趣对象的整体患者位移模型。
在处置期间，辐射射束被引导到患者身体的感兴趣对象。利用医学成像系统来扫描患者以采集图像，例如，超声成像。扫描将针对处置的整个持续时间持续。由医学成像系统生成的这些图像在本发明中被称为“新采集的图像”。
根据本发明，当在处置期间正利用医学成像系统扫描患者时，基于新采集的图像来计算患者身体的感兴趣对象的位移信息，如下面详述的。
图1描绘了根据本发明的方法的流程图，所述方法用于计算感兴趣对象的位移。
所述方法包括根据感兴趣对象的预先采集的图像的集合中相邻的图像来计算101感兴趣对象的位移模型的步骤，其中，所述位移模型反映所述感兴趣对象沿时间的位置。所述位移模型例如能够通过第一或第二上述方法来确定。
感兴趣对象可以包括各种类型的对象，其中，位置信息需要被确定。例如，在辐射治疗的领域中，感兴趣对象对应于患者身体中的组织或器官。所述感兴趣对象可以是体模，在工业领域中有兴趣对其进行监测。
图2描绘了根据本发明的感兴趣对象的位移模型的范例。
水平轴对应于时间并且垂直轴对应于感兴趣对象的位移。例如，位移模型反映患者身体的沿时间变化的感兴趣对象的位置。对于患者身体中的感兴趣对象，所述位移模型通常是周期性曲线，因为移动是由呼吸移动引起的。一个周期性时间循环是T。例如，部分A是位移模型的一个周期。与位移改变较慢的区域P2和P5相比，P1、P3、P4和P6表示位移改变较快的区域。
方法还包括基于所述位移模型来确定102所述预先采集的图像的集合的一个周期性时间循环内的来自所述预先采集的图像的集合的图像的第一子集S1的步骤。
位移模型的一个周期内的时间点被选择。如图2中图示的，在部分A中选择了一些时间点。存在一些区域，例如P1、P2、P3、P4、P5和P6。
在图2中，t1P1、t2P1、t3P1是区域P1中的选定的时间点，而t2P2是区域P2中的选定的时间点。应当有在其他区域中，例如在P3、P4、P5和P6中，选择的一些其他时间点。
在这些选定的时间点处采集的图像形成在下文中称为“图像的第一子集”的图像的集合。
图3图示了图像的第一子集S1与预先采集的图像S之间的关系。图像的第一子集S1是预先采集的图像S的子集。
如图示的，预先采集的图像S中没有选择时间点的图像不被包括在图像的第一子集S1中。
所述方法还包括从新采集的图像识别103图像的第二子集S2的第一步骤，其中，图像的所述第二子集(S2)中的图像是相继的并且在图像的第一子集(S1)中具有相同的最相似图像，其中，通过将新采集的图像中的给定图像与图像的所述第一子集(S1)中的每幅图像进行比较来确定相似性水平的集合，并且其中，最相似图像具有在相似性水平的所述集合中的最大相似性水平。
所述新采集的图像是在患者的处置期间采集的图像。所述最相似图像是针对新采集的图像的在图像的所述第一子集S1中的代表性图像。
对于新采集的图像，通过与图像的所述第一子集S1中的图像进行比较，导出相似性水平。通过与图像的第一子集S1中的所有图像进行比较，导出相似性水平的集合。每个相似性水平表示新采集的图像与图像的所述第一子集S1中的对应的图像之间的相似性。所述最相似图像是从图像的所述第一子集S1选择的。所述最相似图像与新采集的图像之间的相似性水平在相似性水平的所述集合中是最大的。
例如，在图4中，图像In1是新采集的图像。将该新采集的图像In1与图像的所述第一子集S1中的所有图像进行比较。通过将新采集的图像In1与图像的比较第一子集S1中的所有图像进行比较来导出相似性水平的集合。图像Imt在图像的第一子集S1之中具有最大相似性水平。然后，图像Imt被识别为S1中的新采集的图像In1的最相似图像。
对于每幅新采集的图像，以相同方式来识别最相似图像。
可能发生的是，一些新采集的图像在图像的所述第一子集S1中具有相同的最相似图像。
图像的所述第一子集S1的沿时间的采集帧率低于所述预先采集的图像。针对患者身体的相同感兴趣对象，预先采集的图像的采集帧率与新采集的图像的采集帧率相同。因此，图像的第一子集S1的采集帧率低于新采集的图像。例如，在3维体积超声成像中，在使用具有机械扫掠线性阵列的超声探头时，图像采集帧率在[4Hz-15Hz]的范围中，并且在使用具有电操纵矩阵阵列换能器时，在[20Hz-50Hz]。
在下文中，“3D”用于表示3维；“2D”用于表示2维。
例如，如图5图示的，来自图像[I11…I1n]的新采集的图像在图像的第一子集S1中具有相同最相似图像Imt1。因此，图像[I11…I1n]的集合定义图像的第二子集，在图5中被称为S2。
在新采集的图像的图像的第一子集S1中的最相似图像不同于图像Imt1时图像的第二子集S2被完全地识别。例如，当Imt2被识别为新采集的图像I21的图像的第一子集S1中的最相似图像，图像的第二子集S2完成。
所述方法还包括在图像的第二子集S2中的选择104给定图像I1作为第一参考图像I1的第一步骤。
如图5所示，从图像的第二子集S2中选择图像I1作为第一参考图像。
有利地，第一参考图像I1是在图像的第二子集S2的中心时隙中选择的图像。例如，如果图像的第二子集S2的数量是n，如果n是奇数，则第一参考图像I1的次序被选择为图像的第二子集S2中的第(n+1)/2；如果n是偶数，则第一参考图像I1的次序被选择为图像的第二子集S2中的第n/2-1或者第n/2+1。
有利地，第一参考图像I1是在图像的第二子集S2中随机地或任意地选择的图像。
有利地，第一参考图像I1是图像的第二子集S2中的第一图像。
所述方法还包括从新采集的图像识别105图像的第三子集S3的第二步骤，其中，图像的所述第三子集(S3)中的图像是相继的并且在图像的第一子集(S1)中具有相同的最相似图像，其中，通过将新采集的图像中的给定图像与图像的第一子集(S1)中的每幅图像进行比较来确定相似性水平的集合，并且其中，所述最相似图像具有在相似性水平的集合中的最大相似性水平。
所述新采集的图像是在患者的处置期间采集的图像。
所述最相似图像是针对新采集的图像的在图像的所述第一子集S1中的代表性图像。
对于新采集的图像，通过与图像的所述第一子集S1中的图像进行比较，导出相似性水平。通过与图像的第一子集S1中的所有图像进行比较，导出相似性水平的集合。每个相似性水平表示新采集的图像与图像的所述第一子集S1中的对应的图像之间的相似性。所述最相似图像是从图像的所述第一子集S1选择的。所述最相似图像与新采集的图像之间的相似性水平在相似性水平的所述集合中是最大的。
例如，在图5中，图像I21是新采集的图像。将该新采集的图像I21与图像的第一子集S1中的所有图像进行比较。通过比较新采集的图像I21与图像的第一子集S1中的所有图像来导出相似性水平的集合。图像Imt2在图像的第一子集S1之中具有最大相似性水平。然后，图像Imt2被识别为S1中的新采集的图像I21的最相似图像。
对于每幅新采集的图像，以相同方式来识别最相似图像。
可能发生的是，一些新采集的图像在图像的所述第一子集S1中具有相同的最相似图像。
图像的所述第一子集S1的沿时间的采集帧率低于所述预先采集的图像。针对患者身体的相同感兴趣对象，预先采集的图像的采集帧率与新采集的图像的采集帧率相同。因此，图像的第一子集S1的采集帧率低于新采集的图像。例如，在3D体积超声成像中，在使用具有机械扫掠线性阵列的超声探头时，图像采集帧率在[4Hz-15Hz]的范围中，并且在使用具有电操纵矩阵阵列换能器时，在[20Hz-50Hz]。
例如，如图5图示的，来自图像[I21…I2m]的新采集的图像在图像的第一子集S1中具有相同最相似图像Imt2。因此，图像[I21…I2m]的集合定义图像的第三子集，在图5中被称为S3。
在新采集的图像的图像的第一子集S1中的最相似图像不同于图像Imt3时图像的第二子集S3被完全地识别。例如，当Imt3被识别为新采集的图像I31的图像的第一子集S1中的最相似图像时，图像的第三子集S3完成。
所述方法还包括在图像的第三子集S3中的选择106给定图像I2作为第二参考图像I2的第二步骤。
如图5所示，从图像的第三子集S3选择图像I2作为第二参考图像。
有利地，第二参考图像I2是在图像的第三子集S3的中心时隙中的选择的图像。例如，如果图像的第三子集S3的数量是n，如果n是奇数，则第二参考图像I2的次序被选择为图像的第三子集S3中的第(n+1)/2；如果n是偶数，则第二参考图像I2的次序被选择为图像的第三子集S3中的第n/2-1或者第n/2+1。
有利地，第二参考图像I2是在图像的第三子集S3中随机或任意选择的图像。
有利地，第二参考图像I2是图像的第三子集S3中的第一图像。
所述方法还包括计算107第二参考图像I2与第一参考图像I1之间的感兴趣对象的位移的步骤。
所述位移表示第二参考图像I2中的对象的位置与第一参考图像I1中的对象的位置之间的位置改变。
所述位移是通过计算图像I2上的患者身体的感兴趣对象的位置信息与图像I1上的患者身体的感兴趣对象的位置信息之间的距离来测量的。
例如，在超声成像中，两幅图像之间的位移是通过使用本发明中的方法来计算的。在2D超声成像中，其分别在坐标轴(x、y)的2个方向上被计算。在3D超声成像中，其分别在坐标轴(x、y、z)的3个方向上被计算。
上面计算的位移反映平移的像素。通过与坐标轴的每个方向上的已知空间分辨率进行组合，实际移动距离被确定。
有利地，确定102的步骤适于确定从集合[12、14、16、18、20、22、24]中选择的具有若干图像的图像的所述第一子集S1。
例如，如果患者呼吸频率为约20/min并且患者的呼吸循环的持续时间为约3秒。在3D体积超声成像中，在使用具有机械扫掠线性阵列的超声探头时，图像采集帧率在[4Hz-15Hz]的范围中，并且在使用具有电操纵矩阵阵列换能器时，在[20Hz-50Hz]。因此，能够在一个呼吸周期性时间循环内采集的超声图像的数量在使用具有机械扫掠线性阵列的超声探头时在范围[12；45]中，并且在使用具有电操纵矩阵阵列换能器的超声探头时在[60；150]中。
作为范例，实验示出了从集合[12、14、16、18、20、22、24]选择的图像的数量是现实的。
有利地，确定102的步骤包括在位移模型的周期性时间循环内选择第一区域和第二区域，其中，在第一区域中选择的图像的数量大于在第二区域中选择的图像的数量，并且其中，第一区域内的位移模型的斜率大于第二区域中的位移模型的斜率。
如图2图示的，选择一些区域，例如P1、P2、P3、P4、P5和P6。位移模型的斜率指示位移变化。区域P1的斜率和区域P3的斜率大于区域P2的斜率。区域P4的斜率和区域P6的斜率大于区域P5的斜率。
例如，t1P1、t2P1、t3P1是区域P1中的选定的时间点，而t2P2是区域P2中的选定的时间点。在区域P1中选择的时间点的数量大于在区域P2中选择的时间点的数量。同时，区域P1的斜率大于区域P2的斜率。
应当有在其他区域中，例如在P3、P4、P5和P6中，选择的一些其他时间点。在区域P3中选择的时间点的数量大于在区域P2中选择的时间点的数量。同时，区域P3的斜率大于区域P2的斜率。在区域P4、P6中选择的时间点的数量大于在区域P5中选择的时间点的数量。同时，区域P4、P6的斜率大于区域P5的斜率。
有利地，在以上的确定102的步骤之后，所述方法还包括将标记附接到图像的第一集合S1的步骤，所述标记包括感兴趣对象的呼吸呼气信息和呼吸吸气信息。
例如，患者的呼吸时间循环包括呼气阶段和吸气阶段。图像的第一集合S1中的每幅图像被给予标记以反映不同状态。
有利地，识别103的第一步骤中和识别105的第二步骤中的最相似图像是基于使用在体积到体积互相关准则和平方差的和的最小值准则之中选择的准则而被识别的。
体积到体积互相关是已知的图像处理方法(Hermosillo,G.,Chefd’Hotel,D.,Faugeras,O.:Variational methods for multimodal imagematching.Int.J.Comput.Vision 50,329-343(2002))。应用体积到体积互相关是要找到给定图像与图像的第一子集S1中的图像之间的最大相关系数。
平方差的和的最小值准则也是已知的图像处理方法(Yin,Y.,Hoffman,E.,Lin,C.:Mass preserving nonrigid registration of CT lung images using cubic B-spline.Med.Phys.36(9),4213–4222(2009))。
有利地，识别103的第一步骤包括以下子步骤：
-基于与相继的图像进行比较来将标记附接到新采集的图像，所述标记包括呼吸呼气信息和呼吸吸气信息；
-从图像的所述第一子集S1识别图像的第四子集S4，图像的所述第四子集S4具有相同的标记；
-从新采集的图像识别图像的第二子集S2，图像的所述第二子集S2中的图像在图像的第四子集S4中具有相同最相似图像。
所述标记包括如上所述的对象的呼吸呼气信息和呼吸吸气信息。
有利地，识别105的第二步骤包括以下子步骤：
-基于与相邻的图像进行比较来将标记附接到新采集的图像，所述标记包括呼吸呼气信息和呼吸吸气信息；
-从图像的所述第一子集S1识别图像的第四子集S4，其中，图像的所述第四子集S4具有相同的标记；
-从新采集的图像识别图像的第三子集S3，其中，图像的所述第三集合S3中的图像在图像的第四子集S4中具有相同最相似图像。
所述标记包括如上所述的对象的呼吸呼气信息和呼吸吸气信息。
有利地，所述方法还包括以下：
-基于所述位移来补偿108辐射射束旨在被引导到感兴趣对象的的位置的子步骤。
相应地调节辐射射束方向，以便跟随所述感兴趣对象的移动。因此，患者身体呼吸移动被补偿。
图6描绘了计算位移的根据本发明的实施例的系统200的示意图。系统200包括处理单元230，其用于：
-第一单元(231)，其用于根据感兴趣对象的预先采集的图像的集合的相邻的图像来计算感兴趣对象的位移模型，其中，所述位移模型反映所述感兴趣对象沿时间的位置；
-第二单元(232)，其用于基于所述位移模型来确定所预先采集的图像的集合的一个周期性时间循环内的来自所述预先采集的图像的集合的图像的第一子集S1；
-第三单元(233)，其用于从新采集的图像识别图像的第二子集S2，其中，图像的所述第二子集(S2)中的图像在图像的第一子集(S1)中具有相同的最相似图像，其中，通过将新采集的图像中的给定图像与图像的第一子集(S1)中的每幅图像进行比较来确定相似性水平的集合，并且其中，所述最相似图像具有在相似性水平的集合中的最大相似性水平；
-第四单元(234)，其用于在图像的第二子集S2中选择给定图像作为第一参考图像I1；
-第五单元(235)，其用于从新采集的图像识别图像的第三子集S3，其中，图像的识别第三子集(S3)中的图像在图像的所述第一子集(S1)中具有相同的最相似图像，其中，通过将新采集的图像中的给定图像与图像的所述第一子集(S1)中的每幅图像进行比较来确定相似性水平的集合，并且其中，所述最相似图像具有在相似性水平的集合中的最大相似性水平；
第六单元(236)，其用于在图像的所述第三子集S3中选择给定图像作为第二参考图像I2；
第七单元(237)，其用于计算所述第二参考图像I2与所述第一参考图像I1之间的位移。
根据本发明的系统的各种单元适于执行先前描述的各个步骤。例如，一个或多个单元被实施在一个或多个处理器中。
图7示描绘了根据本发明的实施例的系统。
有利地，系统200适于连接到图像采集模块210。
该图像采集模块210旨在采集预先采集的图像和新采集的图像。
有利地，系统200适于连接到控制模块220。
由处理单元230计算的位移通过计算处理单元230递送，并且被发送到控制模块220。控制模块220被激活以递送辐射射束，以便跟随感兴趣对象的移动。
有利地，图像采集模块210是超声成像采集模块。
尽管已经在附图和前面的描述中详细例示和描述了本发明，但这样的例示和描述应当被认为是例示性或示范性的，而非限制性的。本发明不限于所公开的实施例。本领域技术人员通过研究附图、公开内容以及从属权利要求，在实践请求保护的本发明时能够理解并且实现对所公开的实施例的其他变型。
在权利要求中，词语“包括”不排除其他元件或步骤，并且词语“一”或“一个”不排除多个。单个处理器或其他单元可以完成权利要求书中所记载的若干个项目的功能。尽管在互不相同的从属权利要求中记载了特定措施，但是这并不指示不能有利地使用这些措施的组合。权利要求书中的任何附图标记均不应被解释为对范围的限制。

Claims (14)

1.一种计算感兴趣对象的位移的方法，包括根据所述感兴趣对象的预先采集的图像的集合中相邻的图像来计算(101)所述感兴趣对象的位移模型的步骤，所述位移模型反映所述感兴趣对象沿时间的位置，其特征在于，所述方法还包括：
-基于所述位移模型来确定(102)在所述预先采集的图像的集合的一个周期性时间循环内的来自所述预先采集的图像的集合图像的第一子集(S1)的步骤；
-从新采集的图像识别(103)图像的第二子集(S2)的第一步骤，其中，图像的所述第二子集(S2)中的图像是相继的并且在图像的所述第一子集(S1)中具有相同的最相似图像，其中，通过将所述新采集的图像中的给定图像与图像的所述第一子集(S1)中的每幅图像进行比较来确定相似性水平的集合，并且其中，所述最相似图像具有在相似性水平的所述集合中的最大相似性水平；
-在图像的所述第二子集(S2)中选择(104)给定图像作为第一参考图像(I1)的第一步骤；
-从所述新采集的图像识别(105)图像的第三子集(S3)的第二步骤，其中，图像的所述第三子集(S3)中的图像是相继的并且在图像的所述第一子集(S1)中具有相同的最相似图像，其中，通过将所述新采集的图像中的给定图像与图像的所述第一子集(S1)中的每幅图像进行比较来确定相似性水平的集合，并且其中，所述最相似图像具有在相似性水平的所述集合中的最大相似性水平；
-在图像的所述第三子集(S3)中选择(106)给定图像作为第二参考图像(I2)的第二步骤；
-计算(107)所述第二参考图像(I2)与所述第一参考图像(I1)之间的位移的步骤。
2.根据权利要求1所述的方法，其中，确定(102)的所述步骤适于确定从集合[12、14、16、18、20、22、24]中选择的具有若干图像的图像的所述第一子集(S1)。
3.根据权利要求1所述的方法，其中，确定(102)的所述步骤包括在所述位移模型的周期性时间循环内选择第一区域和第二区域，其中，在所述第一区域中选择的图像的数量大于在所述第二区域中选择的图像的数量，并且其中，所述位移模型在所述第一区域内的斜率大于所述位移模型在所述第二区域中的斜率。
4.根据权利要求1所述的方法，其中，在确定(102)的所述步骤之后，还包括将标记附接到图像的所述第一集合(S1)的步骤，所述标记包括所述感兴趣对象的呼吸呼气信息和呼吸吸气信息。
5.根据前述权利要求中的任一项所述的方法，其中，识别(103)的所述第一步骤中和识别(105)的所述第二步骤中的所述最相似图像是基于使用从体积到体积互相关准则和平方差的和的最小值准则中选择的准则而被识别的。
6.根据权利要求1至5所述的方法，其中，选择(104)的所述第一步骤包括在所述新采集的图像的中心时隙中选择所述第一参考图像(I1)。
7.根据权利要求1至5所述的方法，其中，选择(106)的所述第二步骤包括在所述新采集的图像的中心时隙中选择所述第二参考图像(I2)。
8.根据权利要求4所述的方法，其中，识别(103)的所述第一步骤包括以下子步骤：
-基于与相继的图像进行比较来将标记附接到新采集的图像，其中，所述标记包括呼吸呼气信息和呼吸吸气信息；
-从图像的所述第一子集(S1)中识别图像的第四子集(S4)，其中，图像的所述第四子集(S4)具有相同的标记；
-从新采集的图像中识别图像的所述第二子集(S2)，其中，图像的所述第二子集(S2)中的图像在图像的所述第四子集(S4)中具有相同的最相似图像。
9.根据权利要求4所述的方法，其中，识别(105)的所述第二步骤包括以下子步骤：
-基于与相邻的图像进行比较来将标记附接到新采集的图像，其中，所述标记包括呼吸呼气信息和呼吸吸气信息；
-从图像的所述第一子集(S1)中识别图像的第四子集(S4)，其中，图像的所述第四子集(S4)具有相同的标记；
-从新采集的图像中识别图像的所述第三子集(S3)，其中，图像的所述第三子集(S3)中的图像在图像的所述第四子集(S4)中具有相同的最相似图像。
10.根据前述权利要求中的任一项所述的方法，还包括：
-基于所述位移来补偿(108)旨在被引导到所述感兴趣对象的辐射射束的位置的子步骤。
11.一种用于计算位移的系统(200)，所述系统(200)包括处理单元(230)，所述处理单元包括第一单元(231)，所述第一单元用于根据感兴趣对象的预先采集的图像的集合中相邻的图像来计算所述感兴趣对象的位移模型，所述位移模型反映所述感兴趣对象沿时间的位置，其特征在于，所述处理单元(230)还包括：
-第二单元(232)，其用于基于所述位移模型来确定所预述先采集的图像的集合的一个周期性时间循环内的来自所述预先采集的图像的集合的图像的第一子集(S1)；
-第三单元(233)，其用于从新采集的图像中识别图像的第二子集(S2)，其中，图像的所述第二子集(S2)中的图像是相继的并且在图像的所述第一子集(S1)中具有相同的最相似图像，其中，通过将所述新采集的图像中的给定图像与图像的所述第一子集(S1)中的每幅图像进行比较来确定相似性水平的集合，并且其中，所述最相似图像具有在相似性水平的所述集合中的最大相似性水平；
-第四单元(234)，其用于在图像的所述第二子集(S2)中选择给定图像作为第一参考图像(I1)；
第五单元(235)，其用于从所述新采集的图像中识别图像的第三子集(S3)，其中，图像的所述第三子集(S3)中的图像是相继的并且在图像的所述第一子集(S1)中具有相同的最相似图像，其中，通过将所述新采集的图像中的给定图像与图像的所述第一子集(S1)中的每幅图像进行比较来确定相似性水平的集合，并且其中，所述最相似图像具有在相似性水平的所述集合中的最大相似性水平；
-第六单元(236)，其用于在图像的所述第三子集(S3)中选择给定图像作为第二参考图像(I2)；
-第七单元(237)，其用于计算所述第二参考图像(I2)与所述第一参考图像(I1)之间的位移。
12.根据权利要求11所述的系统，其适于连接到图像采集模块(210)。
13.根据权利要求11所述的系统，其适于连接到控制模块(220)。
14.根据权利要求12所述的系统，其中，所述图像采集模块(210)是超声成像采集模块。
CN201580068805.2A 2014-12-17 2015-12-11 用于计算感兴趣对象的位移的方法和系统 Active CN107004270B (zh)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
CNPCT/CN2014/094132 2014-12-17
CN2014094132 2014-12-17
EP15161267 2015-03-27
EP15161267.8 2015-03-27
PCT/EP2015/079359 WO2016096635A1 (en) 2014-12-17 2015-12-11 Method and system for calculating a displacement of an object of interest

Publications (2)

Publication Number Publication Date
CN107004270A true CN107004270A (zh) 2017-08-01
CN107004270B CN107004270B (zh) 2021-02-19

Family

ID=54848569

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201580068805.2A Active CN107004270B (zh) 2014-12-17 2015-12-11 用于计算感兴趣对象的位移的方法和系统

Country Status (5)

Country Link
US (1) US10346989B2 (zh)
EP (1) EP3234917B1 (zh)
JP (1) JP6692817B2 (zh)
CN (1) CN107004270B (zh)
WO (1) WO2016096635A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111374708A (zh) * 2018-12-28 2020-07-07 深圳迈瑞生物医疗电子股份有限公司 一种胎儿心率检测方法及超声成像装置、存储介质

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016096635A1 (en) * 2014-12-17 2016-06-23 Koninklijke Philips N.V. Method and system for calculating a displacement of an object of interest
CN110443794B (zh) * 2019-05-29 2020-12-18 腾讯科技（深圳）有限公司 一种基于病理图像的图像状态确定方法、装置以及系统

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10118061A (ja) * 1996-10-16 1998-05-12 Aloka Co Ltd 超音波診断装置
US6442476B1 (en) * 1998-04-15 2002-08-27 Research Organisation Method of tracking and sensing position of objects
JP2005160688A (ja) * 2003-12-02 2005-06-23 Hitachi Medical Corp 画像診断装置
CN101013503A (zh) * 2007-01-26 2007-08-08 清华大学 一种医学图像中腹部器官分割方法
CN101027677A (zh) * 2004-06-24 2007-08-29 维森盖特有限公司 用于校正连续视图之间的对象－探测器相对运动的方法
CN101107628A (zh) * 2005-01-19 2008-01-16 皇家飞利浦电子股份有限公司 用于图像对准的图像处理系统和方法
WO2008044572A1 (fr) * 2006-10-04 2008-04-17 Hitachi Medical Corporation dispositif de diagnostic d'images médicales
CN102132322A (zh) * 2008-08-28 2011-07-20 皇家飞利浦电子股份有限公司 用于确定对象的尺寸改变的设备
CN102812494A (zh) * 2010-03-12 2012-12-05 皇家飞利浦电子股份有限公司 血管造影图像中的运动可视化

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6563964B1 (en) 1999-02-08 2003-05-13 Sharp Laboratories Of America, Inc. Image downsampling using redundant pixel removal
WO2002037934A2 (en) 2000-06-05 2002-05-16 Mentor Corporation Automated implantation system for radioisotope seeds
US8379794B2 (en) 2008-09-05 2013-02-19 The Board Of Trustees Of The Leland Stanford Junior University Method to estimate position, motion and trajectory of a target with a single x-ray imager
CN102858229B (zh) 2010-02-18 2015-11-25 皇家飞利浦电子股份有限公司 采用受到跟踪的支气管镜检查进行肿瘤运动模拟和运动补偿的系统和方法
WO2012019162A1 (en) 2010-08-06 2012-02-09 Accuray, Inc. Systems and methods for real-time tumor tracking during radiation treatment using ultrasound imaging
JP5906428B2 (ja) 2012-07-11 2016-04-20 パナソニックＩｐマネジメント株式会社 Ｌｅｄモジュールおよびこれを備えた照明装置
WO2014010073A1 (ja) 2012-07-13 2014-01-16 三菱電機株式会社 Ｘ線位置決め装置、ｘ線位置決め方法及び注目画像撮影方法
US20140107390A1 (en) 2012-10-12 2014-04-17 Elekta Ab (Publ) Implementation and experimental results of real-time 4d tumor tracking using multi-leaf collimator (mlc), and/or mlc-carriage (mlc-bank), and/or treatment table (couch)
JP6266887B2 (ja) * 2013-03-04 2018-01-24 株式会社デンソー 推定装置
KR20140126815A (ko) 2013-04-22 2014-11-03 삼성전자주식회사 호흡 주기 동안 체내 장기의 변화를 추적하는 방법, 장치 및 시스템.
US9601625B2 (en) 2013-07-15 2017-03-21 Taiwan Semiconductor Manufacturing Company Limited Guard ring structure of semiconductor arrangement
KR102185362B1 (ko) * 2013-10-08 2020-12-02 삼성전자주식회사 초음파 프로브 및 이를 포함한 의료 장치
JP6182488B2 (ja) * 2014-03-25 2017-08-16 株式会社日立製作所 位置決め制御装置
WO2016096635A1 (en) * 2014-12-17 2016-06-23 Koninklijke Philips N.V. Method and system for calculating a displacement of an object of interest
US10242281B2 (en) * 2015-07-05 2019-03-26 RideOn Ltd. Hybrid orientation system

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10118061A (ja) * 1996-10-16 1998-05-12 Aloka Co Ltd 超音波診断装置
US6442476B1 (en) * 1998-04-15 2002-08-27 Research Organisation Method of tracking and sensing position of objects
JP2005160688A (ja) * 2003-12-02 2005-06-23 Hitachi Medical Corp 画像診断装置
CN101027677A (zh) * 2004-06-24 2007-08-29 维森盖特有限公司 用于校正连续视图之间的对象－探测器相对运动的方法
CN101107628A (zh) * 2005-01-19 2008-01-16 皇家飞利浦电子股份有限公司 用于图像对准的图像处理系统和方法
WO2008044572A1 (fr) * 2006-10-04 2008-04-17 Hitachi Medical Corporation dispositif de diagnostic d'images médicales
CN101013503A (zh) * 2007-01-26 2007-08-08 清华大学 一种医学图像中腹部器官分割方法
CN102132322A (zh) * 2008-08-28 2011-07-20 皇家飞利浦电子股份有限公司 用于确定对象的尺寸改变的设备
CN102812494A (zh) * 2010-03-12 2012-12-05 皇家飞利浦电子股份有限公司 血管造影图像中的运动可视化

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JUNG-BAE KIM 等: "Real-time moving organ tracking in ultrasound video based on a 3D organ model", 《ULTRASONICS SYMPOSIUM (IUS) , 2012 IEEE INTERNATIONAL, IEEE》 *
VALERIA DE LUCA 等: "Speeding-up image registration for repetitive motion scenarios", 《BIOMEDICAL IMAGING (ISBI), 2012 9TH IEEE INTERNATIONAL SYMPOSIUM ON, IEEE》 *
邓梁 等: "基于时变医学先验信息的约束成像及图像配准方法", 《电子与信息学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111374708A (zh) * 2018-12-28 2020-07-07 深圳迈瑞生物医疗电子股份有限公司 一种胎儿心率检测方法及超声成像装置、存储介质
CN111374708B (zh) * 2018-12-28 2024-02-20 深圳迈瑞生物医疗电子股份有限公司 一种胎儿心率检测方法及超声成像装置、存储介质

Also Published As

Publication number Publication date
EP3234917B1 (en) 2019-08-14
JP2018506322A (ja) 2018-03-08
WO2016096635A1 (en) 2016-06-23
EP3234917A1 (en) 2017-10-25
JP6692817B2 (ja) 2020-05-13
US20170345158A1 (en) 2017-11-30
CN107004270B (zh) 2021-02-19
US10346989B2 (en) 2019-07-09

Similar Documents

Publication Publication Date Title
JP5797352B1 (ja) ３次元物体を追跡するための方法
CN102781331B (zh) 使用飞行时间信息来探测和校正成像扫描中的运动的方法和设备
EP3804624B1 (en) Synthetic parameterized computed tomography from surface data in medical imaging
CN101076282B (zh) 移动目标的动态追踪
KR101982149B1 (ko) 의료 영상의 일부 정보를 활용한 장기 영상 생성 방법 및 장치
Wu et al. Detection of the respiratory rate of standing cows by combining the Deeplab V3+ semantic segmentation model with the phase-based video magnification algorithm
US8655040B2 (en) Integrated image registration and motion estimation for medical imaging applications
US8594769B2 (en) Systems and methods for associating physiological data with image data
CN107334486A (zh) 确定造影剂辅助成像方法的至少一个协议参数
US20140286556A1 (en) Heart imaging method
KR20140096919A (ko) 의료 영상 정합 방법 및 장치
CN102908144A (zh) 用于治疗计划的磁共振成像
CN107133549A (zh) Ect运动门控信号获取方法及ect图像重建方法
US11269036B2 (en) System and method for phase unwrapping for automatic cine DENSE strain analysis using phase predictions and region growing
CN104145295B (zh) 使用不同指示器来确定放射治疗目标的位置变化
CN103445799A (zh) 用于确定患者特定的造影剂脉冲响应函数的方法和装置
CN107004270A (zh) 用于计算感兴趣对象的位移的方法和系统
US9717441B2 (en) Automatic method of predictive determination of the position of the skin
CN107249465A (zh) 用于稀疏角度采样的断层摄影成像装置和方法
KR20150129506A (ko) 의료 영상 정합 방법 및 그 장치
US9636076B2 (en) X-ray CT apparatus and image processing method
CN105578963B (zh) 用于组织剂量估计的图像数据z轴覆盖范围延伸
CN104000618A (zh) 一种用环真光子数门控方法进行呼吸运动门控校正技术
KR20190039498A (ko) 자기 공명 용적 영상화
US11837352B2 (en) Body representations

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant