CN115880216A - 提供比较数据集 - Google Patents

提供比较数据集 Download PDF

Info

Publication number
CN115880216A
CN115880216A CN202211175577.7A CN202211175577A CN115880216A CN 115880216 A CN115880216 A CN 115880216A CN 202211175577 A CN202211175577 A CN 202211175577A CN 115880216 A CN115880216 A CN 115880216A
Authority
CN
China
Prior art keywords
data set
image
registered
reg
contrast agent
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.)
Pending
Application number
CN202211175577.7A
Other languages
English (en)
Inventor
A.比克霍尔德
M.科瓦谢克
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.)
Siemens Healthineers AG
Original Assignee
Siemens Healthineers AG
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 Siemens Healthineers AG filed Critical Siemens Healthineers AG
Publication of CN115880216A publication Critical patent/CN115880216A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74Image or video pattern matching; Proximity measures in feature spaces
    • G06V10/75Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
    • G06V10/751Comparing pixel values or logical combinations thereof, or feature values having positional relevance, e.g. template matching
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/20Analysis of motion
    • G06T7/254Analysis of motion involving subtraction of images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/38Registration of image sequences
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/25Determination of region of interest [ROI] or a volume of interest [VOI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10132Ultrasound image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20224Image subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular
    • G06T2207/30104Vascular flow; Blood flow; Perfusion

Abstract

本发明涉及一种用于提供比较数据集的方法,包括提供映射第一造影剂流的第一数据集和映射第二造影剂流的第二数据集;在空间上配准第一数据集和第二数据集;借助第一和/或第二造影剂流的映射的灌注的方向,识别第一数据集和第二数据集中感兴趣区域的至少一个血管区段的映像;在时间上配准第一数据集和第二数据集,以便最小化至少一个血管区段中第一和第二造影剂流的映射的灌注的时间差;通过比较配准的第一数据集和配准的第二数据集识别第一和第二造影剂流之间的偏差;基于配准的第一数据集和配准的第二数据集提供比较数据集;其中,比较数据集具有至少一个表征偏差的参数。本发明还涉及一种提供单元、一种医学成像设备和一种计算机程序产品。

Description

提供比较数据集
技术领域
本发明涉及一种用于提供比较数据集的方法、一种提供单元、一种医学成像设备和一种计算机程序产品。
背景技术
为了检测检查对象、例如人类和/或动物患者的感兴趣区域(英文:region-of-interest;缩写:ROI)的变化,通常使用尤其基于X射线的差分成像方法,例如数字减影血管造影(DSA)。待检测的变化例如可以包括造影剂在检查对象的血管区段中的扩散运动和/或流动运动、尤其造影剂流。在DSA中,通常从一个或多个填充图像中减去掩模图像并提供差分图像,其中,掩模图像不带布置在其中的造影剂的感兴趣区域,而至少一个填充图像映射带有布置在其中的造影剂的感兴趣区域。为了映射造影剂流的时间动态,例如造影剂流的灌注(Anflutung),在感兴趣区域中通常按时间顺序记录多个填充图像并与掩模图像组合成时间分辨的差分图像。
通常实施术前DSA,以便尤其在术前计划对血管畸形和/或狭窄和/或动脉瘤的尤其脑血管治疗。为了还能够在血管畸形和/或狭窄和/或动脉瘤的治疗期间或之后检测造影剂流的时间动态变化,可以在术前或术中实施另外的DSA。通常,医疗操作者在视觉上比较第一DSA和另外的DSA的差分图像。不利地,术前和术中或术后的造影剂施用之间的差异、例如注射时间点和/或注射速率和/或造影剂浓度和/或导管位置可能导致差分图像的错误比较。因此,对造影剂流的时间动态的变化的定量分析可能会被阻碍或篡改。
发明内容
因此,本发明所要解决的技术问题是,能够实现对检查对象的感兴趣区域中造影剂流的变化的时间动态的定量分析。
所述技术问题按照本发明通过一种用于提供比较数据集的方法解决,所述方法包括:
-提供时间分辨的第一数据集,所述第一数据集映射在第一时间段内检查对象的感兴趣区域中的第一造影剂流;
-提供时间分辨的第二数据集,所述第二数据集映射在第一时间段之后的第二时间段内检查对象的感兴趣区域中的第二造影剂流;
-在空间上配准第一数据集和第二数据集;
-借助第一和/或第二造影剂流的映射的灌注的方向,识别第一数据集和第二数据集中感兴趣区域的至少一个尤其传入的血管区段的映像;
-在时间上配准第一数据集和第二数据集,以便最小化所述至少一个血管区段中第一和第二造影剂流的映射的灌注的时间差;
-通过比较配准的第一数据集和配准的第二数据集识别第一和第二造影剂流之间的偏差;
-基于配准的第一数据集和配准的第二数据集提供所述比较数据集;
其中,所述比较数据集具有至少一个表征所述偏差的参数。
所述技术问题按照本发明还通过一种提供单元解决,所述提供单元设计用于实施前述类型的方法。
所述技术问题按照本发明还通过一种医学成像设备解决,所述医学成像设备具有前述类型的提供单元,其中,所述医学成像设备设计用于记录和/或提供第一数据集和第二数据集。
所述技术问题按照本发明还通过一种计算机程序产品解决,所述计算机程序产品具有能够直接加载到提供单元的存储器中的计算机程序,所述计算机程序具有多个程序段,以便当这些程序段由所述提供单元执行时实施前述类型的方法的所有步骤。
在第一方面,本发明涉及一种尤其计算机实现的用于提供比较数据集的方法。在该方法中,提供时间分辨的第一数据集,该第一数据集映射在第一时间段内检查对象的感兴趣区域中的第一造影剂流。此外,提供时间分辨的第二数据集,该第二数据集映射在第一时间段之后的第二时间段内检查对象的感兴趣区域中的第二造影剂流。此外,在另外的步骤中,在空间上配准第一和第二数据集。此外,借助第一和/或第二造影剂流的映射的灌注方向识别第一数据集和第二数据集中感兴趣区域的至少一个尤其传入的血管区段的映像。此外,在时间上配准第一数据集和第二数据集,以便最小化所述至少一个血管区段中第一和第二造影剂流的映射的灌注的时间差。此外,通过比较配准的第一数据集和配准的第二数据集识别第一和第二造影剂流之间的偏差。此后,基于配准的第一数据集和配准的第二数据集提供所述比较数据集。其中,所述比较数据集具有至少一个表征所述偏差的参数。
在此,所建议的用于提供比较数据集的方法的上述步骤可以相继和/或至少部分地同时实施。此外,所建议的方法的步骤可以至少部分地、尤其完全地由计算机实现。
第一数据集的提供可以包括第一数据集的接收和/或生成、尤其重建。此外,第二数据集的提供可以包括第二数据集的接收和/或生成、尤其重建。
第一和/或第二数据集的接收尤其可以包括获取和/或读出计算机可读数据存储器和/或从数据存储单元、例如数据库接收。此外,第一和/或第二数据集可以由一个或多个尤其具有相同或不同成像模态的医学成像设备的提供单元提供。至少一个医学成像设备例如可以设计为磁共振成像系统(MRT)和/或计算机断层成像系统(CT)和/或医用X射线设备和/或正电子发射断层成像系统(PET)和/或超声设备。
备选或附加地,第一数据集可以由预先获取的第一图像数据生成、尤其重建。与此类似,第二数据集可以由预先获取的第二图像数据生成、尤其重建。
时间分辨的第一数据集可以映射检查对象的二维(2D)和/或三维(3D)空间分辨的感兴趣区域。检查对象例如可以是人类和/或动物患者。此外,感兴趣区域可以包括检查对象的空间区域,该空间区域具有解剖对象、例如器官、尤其中空器官,和/或肿瘤,和/或解剖对象的预定义区段。有利地,感兴趣区域可以包括至少一个血管区段、例如动脉和/或静脉。此外,第一数据集可以映射在第一时间段内的感兴趣区域,其中,第一时间段包括多个第一记录时间点。在这种情况下,第一数据集可以映射在第一时间段期间造影剂、尤其造影剂团在感兴趣区域中的第一造影剂流、尤其扩散运动、尤其灌注和/或流动运动。有利地,第一数据集可以具有多个像点,这些像点分别具有时间强度曲线,这些时间强度曲线将感兴趣区域中的第一造影剂流映射为至少一个血管区段中的时间上的强度变化。
时间分辨的第二数据集可以映射检查对象的二维(2D)和/或三维(3D)空间分辨的感兴趣区域。此外,第二数据集可以映射在第二时间段内的感兴趣区域,其中,第二时间段包括多个第二记录时间点。在这种情况下,第二数据集可以映射在第二时间段期间造影剂、尤其造影剂团在感兴趣区域中的第二造影剂流、尤其扩散运动、尤其灌注和/或流动运动。有利地,第二数据集可以具有多个像点,这些像点分别具有时间强度曲线,这些时间强度曲线将感兴趣区域中的第二造影剂流映射为至少一个血管区段中的时间上的强度变化。
第一和第二数据集的空间配准可以例如包括第一和/或第二数据集、尤其各自像点的刚性和/或非刚性变换、尤其平移和/或旋转和/或缩放。有利地,空间配准可以基于映射在第一和第二数据集中的解剖和/或几何特征。解剖特征例如可以包括器官、尤其中空器官,和/或解剖标志、尤其口。此外,几何特征例如可以包括轮廓和/或标记对象。在此,空间配准可以包括确定变换和/或将变换应用于第一和/或第二数据集,使得第一和第二数据集中的相应的解剖特征和/或几何特征的图像之间的空间偏差最小化。
对感兴趣区域的至少一个血管区段的映像的识别可以包括识别、尤其分割第一数据集的像点和第二数据集的像点,这些像点映射至少一个血管区段。在此可以借助感兴趣区域的第一和/或第二造影剂流的映射的灌注方向识别至少一个血管区段、尤其至少一个相对于感兴趣区域和/或感兴趣区域的局部传入、尤其用于供应的血管区段。第一和/或第二造影剂流的灌注可以描述在至少一个血管区段、尤其至少一个血管区段在第一和/或第二数据集中的映像中的造影剂浓度的尤其局部的上升、尤其预定义的强度变化。在这种情况下,灌注的方向可以借助造影剂浓度的上升、尤其预定义的强度变化的时间分辨的映像沿着映射至少一个血管区段的像点、尤其沿着至少一个血管区段的中心线来确定。如果至少一个血管区段中的灌注方向指向感兴趣区域,则传入的、尤其用于供应的血管区段可以被识别为所述至少一个血管区段。
有利地可以借助第一和/或第二造影剂流的映射的灌注方向识别尤其共同血管、例如动脉和/或静脉的一个或多个尤其连续的血管区段。
根据第一变型方案可以借助第一和第二造影剂流的映射的灌注方向识别第一和第二数据集中感兴趣区域的至少一个血管区段的映像。
根据第二变型方案可以借助第一或第二造影剂流的映射的灌注方向识别第一和第二数据集中感兴趣区域的至少一个血管区段的映像。在此,至少一个血管区段在第一和第二数据集中的映像另外可以由于第一和第二数据集的空间配准来识别。
第一和第二数据集的在时间上的配准可以包括基于在至少一个血管区段中第一和第二造影剂流的映射的灌注同步第一和第二数据集的在空间上彼此对应的像点、尤其第一和第二数据集的相应像点的时间强度曲线。在这种情况下,第一和第二数据集的在空间上彼此对应的像点可以映射至少一个血管区段中的共同的空间位置并且尤其也用于同步第一和第二数据集的其它像点的时间强度曲线。此外,至少一个血管中第一和第二造影剂流的映射的灌注的时间差可以通过比较达到或超出第一和第二数据集的在空间上彼此对应的像点的时间强度曲线的预定义强度阈值的时间点来确定。有利地,通过第一和第二数据集的时间配准、尤其同步第一和第二数据集的在空间上彼此对应的像点的时间强度曲线可以使至少一个血管区段中第一和第二造影剂的映射的灌注的时间差最小化。
比较尤其在空间和时间上配准的第一和第二数据集以便识别第一和第二造影剂流之间的偏差可以包括比较第一和第二数据集的在空间上彼此对应的像点的配准的时间强度曲线。通过上述比较可以识别第一和第二造影剂流之间的偏差。尤其可以识别配准的第一和/或第二数据集的映射第一和第二造影剂流之间偏差的像点。有利地可以首先通过比较配准的第一和第二数据集来定性地识别偏差。
此外,比较数据集的提供可以包括配准的第一和第二数据集的定量比较,尤其包括配准的第一和第二数据集的映射第一和第二造影剂流之间偏差的像点之间的定量比较。比较数据集可以至少一个具有表征所述偏差的参数。在此,表征偏差的至少一个参数可以尤其针对整个感兴趣区域全局地、尤其逐像点地局部地或者尤其针对至少一个血管区段区域性地描述、尤其量化该偏差。此外,比较数据集可以具有多个表征偏差的参数,这些参数局部或区域性地描述偏差。此外,比较数据集可以具有检查对象的感兴趣区域的2D和/或3D空间分辨的比较图像,其中,该比较图像的像点映射第一和第二造影剂流之间的偏差、例如差异或比率。此外,比较图像可以是时间分辨的。在这种情况下,比较图像的像点的图像值可以至少部分对应于所述多个表征偏差的参数。
比较数据集的提供可以例如包括存储在计算机可读存储介质上和/或将比较数据集的图形表示显示在显示单元上和/或传输给提供单元。尤其地,比较数据集的提供可以包括通过显示装置显示比较数据集的图形表示、尤其至少一个表征偏差的参数。此外,比较数据集的提供可以包括显示、尤其叠加或并排地(英文:side-by-side)显示配准的第一和第二数据集的图形表示。此外,用于显示图形表示的显示装置的至少一个显示参数可以根据比较数据集、尤其至少一个表征偏差的参数例如通过颜色编码来调整。
所建议的实施方式可以有利地尤其通过提供具有至少一个表征偏差的参数的比较数据集使得能够定量分析第二造影剂流相对于第一造影剂流在检查对象的感兴趣区域中的时间动态的变化。
在所建议的方法的另外的有利实施方式中,第一数据集的提供可以包括接收预先获取的第一图像数据。其中,第一数据集可以由第一图像数据重建。此外,第二数据集的提供可以包括接收预先获取的第二图像数据。其中,第二数据集可以由第二图像数据重建。
第一和/或第二图像数据的接收尤其可以包括获取和/或读出计算机可读数据存储器和/或从数据存储单元、例如数据库接收。此外,第一和/或第二图像数据可以由一个或多个尤其具有相同或不同成像模态的医学成像设备的提供单元提供。
第一图像数据可以映射在第一时间段内、尤其在多个第一记录时间点检查对象的感兴趣区域。第一图像数据可以分别至少部分地、尤其完全地映射感兴趣区域。此外,第一图像数据可以以至少部分不同的映射几何形状、尤其映射方向和/或在不同的第一记录时间点映射感兴趣区域。第一图像数据可以例如具有第一时间段内感兴趣区域的层图像和/或投影图像。
第二图像数据可以映射在第二时间段内、尤其在多个第二记录时间点检查对象的感兴趣区域。第二图像数据可以分别至少部分地、尤其完全地映射感兴趣区域。此外,第二图像数据可以以至少部分不同的映射几何形状、尤其映射方向和/或在不同的第二记录时间点映射感兴趣区域。第二图像数据可以例如具有第二时间段内感兴趣区域的层图像和/或投影图像。
有利地,第一数据集可以由第一图像数据重建。此外,第二数据集可以由第二图像数据重建。在这种情况下,第一数据集的重建可以包括第一图像数据相对彼此的空间和/或时间配准。此外,第一数据集的重建可以包括对第一图像数据应用2D或3D重建、例如逆Radon变换和/或滤波反投影和/或逆傅里叶变换。与此类似,第二数据集的重建可以包括第二图像数据相对彼此的空间和/或时间配准。此外,第二数据集的重建可以包括对第二图像数据应用2D或3D重建、例如逆Radon变换和/或滤波反投影和/或逆傅里叶变换。在此,第一和第二数据集的重建可以相同或不同。此外,第一和/或第二数据集的重建可以包括对相应的图像数据应用图像伪影校正、尤其运动校正和/或金属伪影校正。
所建议的实施方式可以实现检查对象的感兴趣区域中第一和第二造影剂流的改进映射。
在所建议的方法的另外的有利实施方式中,第一图像数据可以具有检查对象的从至少部分不同的、尤其不共线的第一投影方向的多个第一投影图像。其中,时间分辨的第一数据集可以由第一投影图像重建。此外,第二图像数据可以具有检查对象的从至少部分不同的、尤其不共线的第二投影方向的多个第二投影图像。其中,时间分辨的第二数据集可以由第二投影图像重建。
用于记录第一和第二图像数据的至少一个医学成像设备可以有利地具有源和探测器,所述源和探测器能够相对于检查对象、尤其感兴趣区域定位在定义的布置中。在医学成像设备设计为医用X射线设备、尤其医用C形臂X射线设备和/或计算机断层成像系统(CT)的设计方案中,所述源可以是X射线源,并且所述探测器可以是X射线探测器。
至少部分不同的第一和第二投影方向可以分别描述在相应图像数据的记录时间点、尤其在第一或第二记录时间点至少一个医学成像设备的源和探测器、尤其探测器中心点之间的射线、尤其中央射线和/或中心射线的走向。尤其地,第一和第二投影方向可以分别描述至少一个医学成像设备相对于检查对象、尤其感兴趣区域的角度和/或源和探测器的定义的布置的等中心、尤其旋转中心的角度。在这种情况下,等中心点可以描述这样的空间点,源和探测器的定义的布置可以尤其在第一和第二图像数据的记录期间围绕该空间点移动、尤其旋转。有利地,至少部分不同的第一和第二投影方向可以分别延伸穿过尤其共同的等中心。在这种情况下,在第一图像数据的记录期间的等中心可以与在第二图像数据的记录期间的等中心相同或不同地相对于检查对象、尤其感兴趣区域定位。
有利地,多个第一和第二投影图像可以映射检查对象的分别2D空间分辨的感兴趣区域。有利地,第一数据集可以如此由多个第一投影图像重建,使得第一数据集映射3D空间分辨且时间分辨的感兴趣区域。此外,第二数据集可以如此由多个第二投影图像重建,使得第二数据集映射3D空间分辨且时间分辨的感兴趣区域。尤其地,第一和第二数据集可以分别按照4D DSA由相应的投影图像重建。
所建议的实施方式可以实现检查对象的感兴趣区域中第一和第二造影剂流的改进的、尤其少遮盖的映射。
在所建议的方法的另外的有利实施方式中,第一图像数据和/或第二图像数据可以具有至少一个掩模图像。此外,第一图像数据可以具有多个第一填充图像。其中,第一数据集的重建可以包括从多个第一填充图像中减去至少一个掩模图像。此外,第二图像数据可以具有多个第二填充图像。其中,第二数据集的重建可以包括从多个第二填充图像中减去至少一个掩模图像。
有利地,第一和/或第二图像数据可以具有至少一个掩模图像,尤其分别具有至少一个掩模图像。如果第一图像数据具有至少一个掩模图像,则至少一个掩模图像可以映射时间上在第一造影剂流的灌注、尤其第一造影剂流在至少一个血管区段中的灌注之前的感兴趣区域。如果第二图像数据具有至少一个掩模图像,则至少一个掩模图像可以映射时间上在第二造影剂流的灌注、尤其第二造影剂流在至少一个血管区段中的灌注之前的感兴趣区域。此外,至少一个掩模图像可以具有第一和/或第二图像数据、尤其第一和/或第二投影图像的所有特征和特性。有利地,第一和/或第二图像数据可以具有多个掩模图像、尤其针对相对于检查对象的每个映射几何形状分别具有掩模图像。
第一图像数据可以具有多个第一填充图像。在这种情况下,第一填充图像可以映射在第一时间段内、尤其在多个第一记录时间点检查对象的感兴趣区域中的第一造影剂流。所述多个第一填充图像可以具有第一图像数据、尤其第一投影图像的所有特征和特性。
此外,第二图像数据可以具有多个第二填充图像。在这种情况下,第二填充图像可以映射在第二时间段内、尤其在多个第二记录时间点检查对象的感兴趣区域中的第二造影剂流。所述多个第二填充图像可以具有第二图像数据、尤其第二投影图像的所有特征和特性。
根据第一变型,第一数据集可以通过从多个第一填充图像中减去至少一个掩模图像来确定。尤其地,第一数据集可以包括从多个第一填充图像中减去至少一个掩模图像的第一差分图像。此外,第二数据集可以通过从多个第二填充图像中减去至少一个掩模图像来确定。尤其地,第二数据集可以包括从多个第二填充图像中减去至少一个掩模图像的第二差分图像。如果仅第一图像数据或者仅第二图像数据具有至少一个掩模图像,则该至少一个掩模图像可以有利地在相应被减去之前与第一和/或第二填充图像尤其在空间上配准。备选地,第一和第二图像数据可以分别具有至少一个掩模图像、尤其至少一个第一掩模图像和至少一个第二掩模图像。在这种情况下,第一数据集可以通过从多个第一填充图像中减去至少一个第一掩模图像来确定。此外,第二数据集可以通过从多个第二填充图像中减去至少一个第二掩模图像来确定。
根据第二变型,第一和/或第二图像数据可以具有多个掩模图像,这些掩模图像以至少部分不同的映射几何形状、尤其投影方向映射感兴趣区域。在这种情况下可以有利地由多个掩模图像重建至少一个掩模数据集。此外,第一数据集可以通过从多个第一填充图像中减去至少一个掩模数据集来确定。此外,第二数据集可以通过从多个第二填充图像中减去至少一个掩模数据集来确定。如果第一图像数据和第二图像数据分别具有多个掩模图像,则可以由掩模图像重建相应的掩模数据集、尤其第一掩模数据集和第二掩模数据集。在这种情况下,第一数据集可以通过从多个第一填充图像中减去第一掩模数据集来确定。此外,第二数据集可以通过从多个第二填充图像中减去第二掩模数据集来确定。
通过分别从第一和第二填充图像中减去至少一个掩模图像、尤其至少一个掩模数据集,可以有利地去除、尤其掩蔽映射感兴趣区域的在时间和/或空间上不可改变的区域的图像区域。由此可以有利地实现,第一数据集基本上仅映射第一时间段内的第一造影剂流。由此还可以实现,第二数据集基本上仅映射第二时间段内的第二造影剂流动。
在所建议的方法的另外的有利实施方式中,配准的第一和第二数据集可以具有多个像点,这些像点分别具有时间强度曲线。其中,对第一和第二造影剂流之间的偏差的识别可以包括对配准的第一和第二数据集的在空间上彼此对应的像点的时间强度曲线的比较,尤其是对时间强度曲线的斜率、方差、平均强度值、最大强度值和/或累积强度值的比较。此外,可以借助配准的第一和第二数据集的在空间上彼此对应的像点的时间强度曲线的比较来确定所述至少一个表征偏差的参数,所述参数描述填充延迟和/或流量比率。
配准的第一和第二数据集的在空间上彼此对应的像点的时间强度曲线的比较可以例如在第一和第二时间段内的相应灌注开始之后的至少一个时间点进行,例如通过比较时间强度曲线在至少一个时间点的强度值。备选或附加地,配准的第一和第二数据集的在空间上彼此对应的像点的时间强度曲线的比较可以在第一和第二时段内的相应灌注开始之后的预定义时间段上进行,例如通过比较时间强度曲线在预定义时间段内的平均强度值和/或累积强度值和/或强度值变化。备选或附加地,配准的第一和第二数据集的在空间上彼此对应的像点的时间强度曲线的比较可以在预定义的强度值范围内进行,例如通过比较从配准的第一和第二数据集的在空间上彼此对应的像点的时间强度曲线的强度值范围的下限开始直至达到或超出强度值范围的上限的时长。
配准的第一和第二数据集的在空间上彼此对应的像点的时间强度曲线的斜率可以映射第一和第二造影剂流的相应灌注的速度。此外,配准的第一和第二数据集的在空间上彼此对应的像点的时间强度曲线的方差可以映射感兴趣区域的灌注。此外,平均强度值可以映射第一和第二造影剂流的尤其平均的体积流量。此外,最大强度值可以映射第一和第二造影剂流的尤其最大的体积流量。此外,累积强度值可以映射在预定义时间段期间第一和第二造影剂流的造影剂体积。通过比较配准的第一和第二数据集的在空间上彼此对应的像点的时间强度曲线的斜率、方差、平均强度值和/或累积强度值可以有利地识别、尤其量化第一和第二造影剂流之间的偏差。
有利地,确定至少一个表征偏差的参数可以包括通过比较配准的第一和第二数据集的在空间上彼此对应的像点的时间强度曲线、尤其时间强度曲线的斜率、方差、平均强度值、最大强度值和/或累积强度值来尤其局部地、区域性地或全局地量化偏差。尤其地,至少一个表征偏差的参数可以映射配准的第一和第二数据集的在空间上彼此对应的像点的时间强度曲线的方差、平均强度值、最大强度值和/或累积强度值的差异或比率。此外,确定至少一个表征偏差的参数可以包括基于配准的第一和第二数据集的在空间上彼此对应的像点的时间强度曲线的斜率、方差、平均强度值和/或累积强度值的偏差来量化在第一和第二时间段之间血流动力学参数、尤其血流动力学参数的变化。在这种情况下,至少一个表征偏差的参数例如可以描述填充延迟和/或流量比率。在此,填充延迟可以描述在第一和第二造影剂流的灌注直至达到预定义的、尤其标准化的强度值之间的时长的时间差。此外,流量比率可以描述第一和第二造影剂流的尤其最大或平均的体积流量的尤其标准化的比率。有利地,填充延迟可以借助配准的第一和第二数据集的在空间上彼此对应的像点的时间强度曲线的斜率和/或累积强度值的比较来确定。此外,流量比率可以借助平均强度值、最大强度值和/或累积强度值来确定。此外,至少一个表征偏差的参数可以描述血管壁的血流动力学压力和/或剪切应力和/或至少一个血管区段中的流速。为此可以另外基于流动模型、例如计算流体动力学(英文:computational fluid dynamics;缩写:CFD)和/或通过借助配准的第一和第二数据集确定例如用于灌注速度的边界条件来确定至少一个表征偏差的参数。
所建议的实施方式可以有利地使得能够通过生物物理量、例如填充延迟和/或流量比率量化第一和第二造影剂流的不同的时间动态。
在所建议的方法的另外的有利实施方式中,可以在配准的第一和第二数据集中分别确定至少一个第一像点和至少一个第二像点,这些像点映射所述至少一个血管区段中的两个不同的空间位置。此外,由至少一个第二像点映射的空间位置可以相对于由至少一个第一像点映射的空间位置布置在下游。其中,可以将造影剂团到达时间(英文:bolusarrival time;缩写:BAT)分别确定为第一和第二造影剂流在由至少一个第二像点映射的空间位置处灌注的时间点。此外可以基于至少一个第二像点的造影剂团到达时间和配准的第一和第二数据集的至少一个第一像点的时间强度曲线的比较、尤其时间强度曲线的斜率、方差、平均强度值、最大强度值和/或累积强度值的比较识别所述偏差。
在配准的第一和第二数据集中分别确定至少一个第一像点和至少一个第二像点可以手动、半自动或全自动进行。例如,至少一个第一像点和/或至少一个第二像点可以借助医疗操作者通过输入单元的输入、尤其借助配准的第一和/或第二数据集的图形表示预先确定。尤其可以借助医疗操作者的输入仅预先确定至少一个第一像点或者仅预先确定至少一个第二像点。在这种情况下,可以例如借助解剖图谱和/或至少一个血管区段的分割自动确定分别另外的至少一个像点。备选地可以借助医疗操作者的输入既预先确定至少一个第一像点又预先确定至少一个第二像点。此外,至少一个第一像点和至少一个第二像点可以全自动地、例如借助在配准的第一和第二数据集中感兴趣区域和至少一个血管区段的分割来确定。
所述两个不同的空间位置可以分别由配准的第一和第二数据集的一个或多个像点映射,尤其由至少一个第一像点和至少一个第二像点映射。此外,由至少一个第二像点映射的空间位置可以尤其沿至少一个血管区段中的流动方向和/或灌注方向相对于由至少一个第一像点映射的空间位置布置在下游。由此,由至少一个第一像点映射的空间位置可以在时间上分别在由至少一个第二像点之前被第一和第二造影剂流灌注(angeflutet)。
有利地,第一和第二造影剂流在由至少一个第二像点映射的空间位置处的造影剂团到达时间可以通过将相应的时间强度曲线与预定强度阈值进行比较来确定。达到或超出预定强度阈值可以表示第一或第二造影剂流在由至少一个第二像点映射的空间位置处的灌注的开始。有利地,第一和第二造影剂流在由至少一个第二像点映射的空间位置处的造影剂团到达时间可以识别配准的第一和第二数据集的至少一个第一像点的时间强度曲线的具有最大斜率的点和/或区段。有利地可以基于配准的第一和第二数据集的至少一个第一像点的时间强度曲线在识别到的点处的比较来识别所述偏差。
所建议的实施方式可以实现第一和第二造影剂流的时间动态的精确比较。
在所建议的方法的另外的有利实施方式中,可以在配准的第一和第二数据集中分别确定至少一个第一像点和至少一个第二像点,这些像点映射所述至少一个血管区段中的两个不同的空间位置。在这种情况下,对第一和第二造影剂流之间的偏差的识别可以包括确定配准的第一和第二数据集中的至少一个第一像点的时间强度曲线、尤其时间强度曲线的斜率、方差、平均强度值、最大强度值和/或累积强度值之间的第一比率。此外,对所述偏差的识别可以包括确定配准的第一和第二数据集中的至少一个第二像点的时间强度曲线、尤其时间强度曲线的斜率、方差、平均强度值、最大强度值和/或累积强度值之间的第二比率。此外,对所述偏差的识别可以包括确定第一比率和第二比率之间的第三比率,其中,结合第三比率识别偏差。
由于所述两个不同空间位置沿至少一个血管区段的不同空间布置,第一和第二造影剂流在这两个不同空间位置处灌注的时间点可以彼此不同。尤其地,由至少一个第一像点映射的空间位置可以在时间上在由至少一个第二像点映射的空间位置之前分别被第一和第二造影剂流灌注。
配准的第一和第二数据集中的至少一个第一像点的时间强度曲线之间的第一比率的确定可以包括确定配准的第一数据集中的至少一个第一像点的时间强度曲线和配准的第二数据集中的至少一个第一像点的时间强度曲线之间的第一商或第一差。因此,第一比率可以映射第一和第二造影剂流之间的在由至少一个第一像点映射的空间位置处的偏差。
此外,配准的第一和第二数据集中的至少一个第二像点的时间强度曲线之间的第二比率的确定可以包括确定配准的第一数据集中的至少一个第二像点的时间强度曲线和配准的第二数据集中的至少一个第二像点的时间强度曲线之间的第二商或第二差。因此,第二比率可以映射第一和第二造影剂流之间的在由至少一个第二像点映射的空间位置处的偏差。
第三比率的确定可以包括确定第一比率和第二比率之间的商或差,尤其第一商和第二商之间或第一差和第二差之间的商或差。在此,除了第一和第二造影剂流之间的在由至少一个第一像点和至少一个第二像点映射的空间位置处的偏差之外,第三比率还可以映射第一和第二造影剂流在两个不同空间位置之间的空间偏差。
所建议的实施方式可以有利地实现第一和第二造影剂流之间的时空比较。
在所建议的方法的另外的有利实施方式中,感兴趣区域可能具有血管畸形、尤其动静脉畸形和/或静脉畸形和/或毛细血管畸形,和/或狭窄和/或动脉瘤。其中,所述至少一个血管区段可以相对于血管畸形和/或狭窄和/或动脉瘤是传入的、尤其用于供应的和/或传出的、尤其用于导出的。在此,分别相对于血管畸形和/或狭窄和/或动脉瘤,至少一个第一像点映射内部或近端的空间位置,并且至少一个第二像点映射远端的空间位置。备选地,分别相对于血管畸形和/或狭窄和/或动脉瘤,至少一个第一像点映射近端的空间位置,并且至少一个第二像点映射内部的空间位置。
相对于血管畸形和/或狭窄和/或动脉瘤的近端空间位置可以表示至少一个血管区段中的这样的空间位置,该空间位置在时间上在血管畸形和/或狭窄和/或动脉瘤之前被第一和第二造影剂流灌注。尤其地,近端空间位置可以布置在至少一个血管区段的传入部分中。在这种情况下,相对于血管畸形和/或狭窄和/或动脉瘤近端布置的空间位置可以在至少一个血管区段中布置在血管畸形和/或狭窄和/或动脉瘤的尤其直接的空间邻近处,尤其与血管畸形和/或狭窄和/或动脉瘤邻接或者在空间上间隔布置。
此外,相对于血管畸形和/或狭窄和/或动脉瘤的远端空间位置可以表示至少一个血管区段中的这样的空间位置,该空间位置在时间上在血管畸形和/或狭窄和/或动脉瘤之后被第一和第二造影剂流灌注。尤其地,远端空间位置可以布置在至少一个血管区段的传出部分中。在这种情况下,相对于血管畸形和/或狭窄和/或动脉瘤远端布置的空间位置可以在至少一个血管区段中布置在血管畸形和/或狭窄和/或动脉瘤的尤其直接的空间邻近处,尤其与血管畸形和/或狭窄和/或动脉瘤邻接或者在空间上间隔布置。
有利地,可以布置在血管畸形和/或狭窄和/或动脉瘤的近端、远端或内部的两个不同的空间位置可以通过至少一个血管区段在血流动力学上连接。
由于所述两个不同空间位置在至少一个血管区段的不同空间布置,第一和第二造影剂流在这两个不同空间位置处灌注的时间点可以彼此不同。尤其地,由至少一个第一像点映射的空间位置可以在时间上在由至少一个第二像点映射的空间位置之前分别被第一和第二造影剂流灌注。
有利地,在第一和第二时间段之间、尤其在所述方法开始之前可以已经发生血管畸形、狭窄和/或动脉瘤的尤其治疗和/或手术和/或介入性的变化。因此,在血管畸形和/或狭窄和/或动脉瘤的内部或远端,流动动态、尤其至少一个血流动力学参数可能已经发生变化。有利地,改变的流动动态可以由第三比率映射,尤其地,注射动态从改变的流动动态中被清除。
所建议的实施方式可以有利地实现第一和第二造影剂流的定量比较,以便评估血管畸形、狭窄和/或动脉瘤的变化。
在所建议的方法的另外的有利实施方式中,可以分别接收第一和第二造影剂流的注射参数。此外,多个时间强度曲线可以基于注射参数被标准化。其中,对第一和第二造影剂流之间的偏差的识别可以包括配准的第一和第二数据集的在空间彼此对应的像点的标准化的时间强度曲线的比较。
第一和第二造影剂流的注射参数的接收尤其可以包括获取和/或读出计算机可读取的数据存储器和/或从数据存储单元、例如数据库接收。此外,注射参数可以由注射装置、尤其自动注射装置的提供单元提供。备选或附加地可以借助医疗操作者通过输入单元的另外的输入接收注射参数。
第一和第二造影剂流的注射参数可以分别描述注射速率、尤其体积流量和/或注射压力和/或造影剂浓度和/或注射时间点和/或空间导管位置。通过基于注射参数对时间强度曲线进行标准化可以识别第一和第二造影剂流之间的偏差,第一和第二造影剂流的注射动态从该偏差中被清除。
在第二方面,本发明涉及一种提供单元,该提供单元设计用于实施所建议的用于提供比较数据集的方法。
所建议的提供单元的优点基本上与所建议的用于提供比较数据集的方法的优点相符。在此提到的特征、优点或备选实施方式也可以转用于其它要求保护的技术方案,反之亦然。
提供单元可以有利地包括计算单元、存储单元和/或接口。在此,提供单元、尤其提供单元的部件可以设计用于实施所建议的用于提供比较数据集的方法的各个单独的步骤。尤其地,接口可以设计用于提供时间分辨的第一和第二数据集并且用于提供比较数据集。此外,计算单元和/或存储单元可以设计用于在空间和时间上配准第一和第二数据集以及用于识别第一和第二造影剂流之间的偏差。
在第三方面,本发明涉及一种医学成像设备,其具有所建议的提供单元,其中,该医学成像设备设计用于记录和/或提供第一数据集和第二数据集。
医学成像设备例如可以设计为磁共振成像系统(MRT)和/或计算机断层成像系统(CT)和/或医用X射线设备和/或正电子发射断层成像系统(PET)和/或超声设备。此外,医学成像设备可以设计用于记录包括感兴趣区域在内的检查对象的第一和第二数据集。备选地,医学成像设备可以设计用于记录包括感兴趣区域在内的检查对象的第一和第二图像数据、尤其第一和第二投影图像。在这种情况下,医学成像设备还可以设计用于由第一图像数据重建和提供第一数据集。此外,医学成像设备可以设计用于由第二图像数据重建和提供第二数据集。
所建议的医学成像设备的优点基本上与所建议的用于提供比较数据集的方法和/或所建议的提供单元的优点相符。在此提到的特征、优点或备选实施方式也可以应用于其它要求保护的技术方案,反之亦然。
在第四方面,本发明涉及一种计算机程序产品,其具有能够直接加载到提供单元的存储器中的计算机程序,所述计算机程序具有多个程序段,以便当这些程序段由提供单元执行时实施用于提供比较数据集的方法的所有步骤。
本发明还可以涉及一种计算机可读存储介质,在所述计算机可读取的存储介质上存储有能够由提供单元读取和执行的多个程序段,以便当这些程序段由提供单元执行时实施用于提供比较数据集的方法的所有步骤。
很大程度上基于软件的实现具有以下优点:当前已经使用的提供单元和/或训练单元也可以以简单的方式通过软件更新改装,以便以根据本发明的方式工作。除了计算机程序之外,这样的计算机程序产品必要时可以包括附加的组成部分、例如文档和/或附加部件以及硬件部件、例如用于软件使用的硬件密钥(加密狗等)。
附图说明
在附图中示出并且以下更详细说明本发明的实施例。不同图中的相同特征使用相同的附图标记。在附图中:
图1示意性地示出所建议的用于提供比较数据集的方法的有利的实施方式;
图2示意性地示出配准的第一和第二数据集的在空间上彼此对应的时间强度曲线;
图3至图5示意性示出所建议的用于提供比较数据集的方法的另外有利实施方式;
图6示意性地示出具有动脉瘤和至少一个血管区段的感兴趣区域;
图7和图8示意性地示出沿血管区段的不同空间位置的时间强度曲线;
图9示意性地示出所建议的提供单元;
图10示意性地示出所建议的医学成像设备。
具体实施方式
图1示意性地示出所建议的用于提供PROV-CD比较数据集CD的方法的一种有利实施方式。在该方法中可以提供PROV-D1时间分辨的第一数据集D1,该第一数据集映射在第一时间段内检查对象的感兴趣区域中的第一造影剂流。此外可以提供PROV-D2时间分辨的第二数据集D2,该第二数据集映射在第一时间段之后的第二时间段内检查对象的感兴趣区域中的第二造影剂流。此外可以将第一数据集D1和第二数据集D2在空间上彼此配准SREG-D1-D2。此外可以借助第一和/或第二造影剂流的映射的灌注方向识别ID-V第一数据集D1和第二数据集D2中感兴趣区域的至少一个尤其传入的血管区段的映像。此外可以在时间上配准第一数据集D1和第二数据集D2,以便最小化至少一个血管区段中第一和第二造影剂流的映射的灌注的时间差。此外可以通过比较配准的第一数据集D1.REG和配准的第二数据集D2.REG识别ID-DIFF第一和第二造影剂流之间的偏差。此外可以基于配准的第一数据集D1.REG和配准的第二数据集D2.REG提供PROV-CD比较数据集CD,其中,比较数据集CD具有至少一个表征偏差的参数P.DIFF。
图2示出配准的第一数据集D1.REG和配准的第二数据集D2.REG的在空间上彼此对应的时间强度曲线的示意图。有利地,配准的第一数据集D1.REG和配准的第二数据集D2.REG可以具有多个像点,这些像点分别具有时间强度曲线ID1.REG和ID2.REG。在此,对第一和第二造影剂流之间的偏差的识别ID-DIFF可以包括配准的第一数据集D1.REG和配准的第二数据集D2.REG的在空间上彼此对应的像点的时间强度曲线ID1.REG和ID2.REG的比较,尤其时间强度曲线ID1.REG和ID2.REG的斜率、方差、平均强度值、最大强度值和/或累积强度值的比较。此外,所述至少一个表征偏差的参数P.DIFF可以借助配准的第一数据集D1.REG和配准的第二数据集D2.REG的在空间上彼此对应的像点的时间强度曲线ID1.REG和ID2.REG的比较来确定,所述参数描述填充延迟和/或流量比率。
通过第一数据集D1.REG和第二数据集D2.REG的空间和时间配准,在空间上彼此对应的时间强度曲线ID1.REG和ID2.REG可以相对于第一和第二造影剂流的分别映射的灌注、尤其相对于造影剂团到达时间tBAT在时间上配准。有利地可以比较在空间上彼此对应的时间强度曲线ID1.REG和ID2.REG,以便确定表征偏差的参数P.DIFF。为此,例如可以确定造影剂团到达时间tBAT和直至达到或超出时间强度曲线ID1.REG和ID2.REG的预定强度阈值的时间点之间的时间差。例如可以为强度阈值预先确定最大强度值、尤其填充强度值IFill的33%的强度值。在这种情况下可以确定造影剂团到达时间tBAT和相应的直至达到或超出预定强度阈值、尤其33%IFill的时间点t33%,D1.REG、t33%,D2.REG之间的时间差。备选地可以确定造影剂团到达时间tBAT和直至达到或超出时间强度曲线下方面积的预定面积阈值、例如从造影剂团到达时间tBAT开始的时间强度曲线下方面积的5%的时间点之间的时间差。通过比较时间差(t33%D1.REG-tBAT)和(t33%,D2.REG-tBAT)可以识别在第一和第二时间段之间至少一个血管区段的填充率的变化。
图3示出所建议的用于提供PROV-CD比较数据集CD的方法的另外的有利实施方式的示意图。其中,第一数据集D1的提供PROV-D1包括接收REC-ID1预先获取的第一图像数据ID1。此外可以由第一图像数据ID1重建RECO-D1第一数据集D1。此外,第二数据集D2的提供PROV-D2可以包括接收REC-ID2预先获取的第二图像数据D2。其中,第二数据集D2可以由第二图像数据ID2重建RECO-D2。
有利地,第一图像数据ID可以具有检查对象的从至少部分不同的第一投影方向的多个第一投影图像。其中,时间分辨的第一数据集D1可以由第一投影图像重建RECO-D1。此外,第二图像数据ID2可以具有检查对象的从至少部分不同的第二投影方向的多个第二投影图像。其中,时间分辨的第二数据集D2可以由第二投影图像重建RECO-D2。
图4示意性地示出所建议的用于提供PROV-CD比较数据集CD的方法的另外的有利实施方式。在此,第一图像数据ID1和/或第二图像数据ID2可以具有至少一个掩模图像,尤其分别具有掩模图像MI1和MI2。此外,第一图像数据ID1可以具有多个第一填充图像FI1。其中,第一数据集D1的重建RECO-D1可以包括从多个第一填充图像FI1中减去DIFF-MI1-FI1至少一个掩模图像MI1。此外,第二图像数据ID2可以具有多个第二填充图像FI2。其中,第二数据集D2的重建RECO-ID2可以包括从多个第二填充图像FI2中减去DIFF-MI2-FI2至少一个掩模图像MI2。
图5示出所建议的用于提供PROV-CD比较数据集CD的方法的另外的有利实施方式的示意图。在此可以分别接收REC-IP第一和第二造影剂流的注射参数IP。此外,时间强度曲线可以基于注入参数IP被标准化。此外,对第一和第二造影剂流之间的偏差DIFF的识别ID-DIFF可以包括配准的第一数据集D1.REG和配准的第二数据集D2.REG的在空间彼此对应的像点的标准化时间强度曲线的比较。
在图6中示意性地示出具有动脉瘤AV和至少一个血管区段V的示例性的感兴趣区域。至少一个血管区段V可以相对于动脉瘤AV是传入的、尤其用于供应的和/或传出的、尤其用于导出的。在此,在图6中通过箭头示出至少一个血管区段V中的流动方向FD、尤其第一和第二造影剂流的灌注方向。
有利地可以在配准的第一数据集D1.REG和配准的第二数据集D2.REG中分别确定至少一个第一像点和至少一个第二像点,这些像点映射至少一个血管区段V中的两个不同的空间位置。其中,分别相对于动脉瘤AV,至少一个第一像点可以映射内部P2或近端P1的空间位置,并且至少一个第二像点可以映射远端的空间位置P3。备选地,分别相对于动脉瘤AV,至少一个第一像点可以映射近端的空间位置P1,并且至少一个第二像点可以映射内部的空间位置P2。在下文中,这三个不同的空间位置P1、P2和P3也称为第一空间位置P1、第二空间位置P2和第三空间位置P3。
有利地,对第一和第二造影剂流之间的偏差的识别ID-DIFF可以包括确定配准的第一数据集D1.REG和第二数据集D2.REG中的至少一个第一像点的时间强度曲线、尤其时间强度曲线的斜率、方差、平均强度值、最大强度值和/或累积强度值之间的第一比率。此外,对所述偏差的识别ID-DIFF可以包括确定配准的第一数据集D1.REG和第二数据集D2.REG中的至少一个第二像点的时间强度曲线、尤其时间强度曲线的斜率、方差、平均强度值、最大强度值和/或累积强度值之间的第二比率。此外可以确定第一比率和第二比率之间的第三比率,其中,借助第三比率识别ID-DIFF偏差。
在图7中示意性地示出配准的第一数据集D1.REG的像点的时间强度曲线ID1.REG,P1、ID1.REG,P2和ID1.REG,P3,这些像点映射至少一个血管区段V中在图6中示出的不同空间位置P1、P2和P3。以下描述的对配准的第一数据集D1.REG的像点的时间强度曲线ID1.REG,P1、ID1.REG,P2和ID1.REG,P3的比较和/或分析可以类似地转用于配准的第二数据集D2.REG的像点的时间强度曲线ID2.REG,P1、ID2.REG,P2和ID2.REG,P3(在此未示出)。
有利地,时间强度曲线ID1.REG,P1、ID1.REG,P2和ID1.REG,P3可以在时间上相对于第一造影剂流的分别映射的灌注、尤其相对于造影剂团到达时间tBAT在时间上配准。
为了确定第一和第二比率,可以确定造影剂团到达时间tBAT和直至达到或超出相应的时间强度曲线的预定强度阈值的时间点之间的时间差。例如可以为强度阈值预先确定最大强度值、尤其填充强度值IFill的33%的强度值。
例如,分别相对于动脉瘤AV,至少一个第一像点可以映射近端的空间位置P1,并且至少一个第二像点可以映射远端的空间位置P3。在这种情况下可以确定造影剂团到达时间tBAT和相应的直至达到或超出至少一个第一像点的时间强度曲线ID1.REG,P1和ID2.REG,P1的预定强度阈值、尤其33%IFill的时间点t33%,D1.REG,P1、t33%,D2.REG,P1之间的时间差:
RP1=(t33%,D2.REC,P1-t33%,D1.REG,P1) (1).此外,第二比率可以确定为造影剂团到达时间tBAT和相应的直至达到或超出至少一个第二像点的时间强度曲线ID1.REG,P3和ID2.REG,P3的预定强度阈值、尤其33%IFill的时间点t33%,D1.REG,P3、t33%,D2.REG,P3之间的时间差:
RP3=(t33%,D2,REG,P3-t33%,D1,REG,P3) (2).
此外,第三比率DRP3,P1可以确定为第二比率RP3和第一比率RP1之差:
DRP3,P1=RP3-RP1
=(t33%,D2,REG,P3-t33%,D1,REG,P3)-(t33%,D2,REG,P1-t33%,D1,REG,P1) (3).
在这种情况下,第三比率DRP3,P1可以描述第三空间位置P3相对于第一空间位置P1在第一和第二时间段之间改变的填充延迟。至少一个表征偏差的参数P.DIFF可以例如包括第三比率DRP3,P1
备选地,分别相对于动脉瘤AV,至少一个第一像点可以映射内部的空间位置P2,并且至少一个第二像点可以映射远端的空间位置P3。在这种情况下可以使第一比率RP2的计算相应地适应由至少一个第一像点映射的第二空间位置P2。此外,第三比率DRP3,P2可以描述第三空间位置P3相对于第二空间位置P2在第一和第二时间段之间改变的填充延迟:
DRP3,P2=(RP3-RPZ)
=(t33%,D2,REG,P3-t33%D1,REG,P3)-(t33%,D2,REG,P2-t33%,D1,REG,P2) (4).
图8示出配准的第一数据集D1.REG的像点的时间强度曲线ID1.REG,P1、ID1.REG,P2和ID1.REG,P3的示意图,这些像点映射至少一个血管区段V中在图6中示出的不同空间位置P1、P2和P3。在此,时间强度曲线ID1.REG,P1、ID1.REG,P2和ID1.REG,P3仅为了图解以相同的斜率示出。此外,时间强度曲线ID1.REG,P1、ID1.REG,P2和ID1.REG,P3在该实施方式中未在时间上彼此配准。
为了确定空间上的第一感兴趣位置IP1、尤其第一空间位置P1和/或第二空间位置P2和在至少一个血管区段V中相对于第一感兴趣位置IP1布置在下游的第一参考位置RP1、尤其第三空间位置P3之间的相对的流变化,可以应用菲克定律、尤其质量守恒:
Figure BDA0003864187150000211
Figure BDA0003864187150000212
其中,QIP1可以表示造影剂的累积质量,VIP1表示体积流量,CIP1和CRP1表示相应空间位置IP1或RP1处的造影剂浓度。
当空间上的第一感兴趣位置IP1和空间上的第一参考位置RP1处的造影剂浓度之差变为最大时,造影剂的累积速度可以达到局部最大值。该局部最大值可以在第一和/或第二造影剂流在下游的第一参考位置RP1处灌注的时间点、尤其到达空间上的第一参考位置RP1的造影剂团到达时间tBAT,RP1达到。
因此,由等式(6)得到:
Figure BDA0003864187150000213
Figure BDA0003864187150000214
由于CIP1∝IIP1,因此
CIP1(tBAT,RP1)=k·IIP1(tBAT,RP1) (9).
由等式(8)和(9)得出:
Figure BDA0003864187150000221
其中,Vrel,IP1=K·VIP1
此外,由等式(7)和(9)得出:
Figure BDA0003864187150000222
在配准的第一数据集D1.REG和第二数据集D2.REG中映射的相对体积流之间的比率根据等式(10)得出:
Figure BDA0003864187150000223
在这种情况下,
Figure BDA0003864187150000224
描述至少一个血管区段V中的第一感兴趣位置IP1处在第一和第二时间段内的体积流量的比率。
假设第一和第二造影剂流的注射参数是恒定的并且在第一和第二时间段之间至少一个血管区段中的血管体积保持基本上相同,可以假设
Figure BDA0003864187150000225
由等式(11)至(13)得出:
Figure BDA0003864187150000231
/>
通过将用于空间上的第一感兴趣位置IP1的相对体积流量
Figure BDA0003864187150000232
与空间上的第二感兴趣位置IP2的相对体积流量/>
Figure BDA0003864187150000233
进行比较可以使用等式(14)确定流量FRIP1,IP2
Figure BDA0003864187150000234
在此,空间上的第二感兴趣位置IP2可以表示至少一个血管区段V中的不同于空间上的第一感兴趣位置IP1的另外的空间位置。此外,空间上的第二参考位置RP2可以在至少一个血管区段V中布置在空间上的第二感兴趣位置IP2的下游。尤其地,第一参考位置RP1和第二参考位置RP2可以相同或不同。
对于第二空间位置P2作为第一感兴趣位置IP1、第三空间位置P3作为第一参考位置RP1并且作为第二感兴趣位置IP2,得到相对流量FRP2,P3
Figure BDA0003864187150000241
其中,RP2布置在第三空间位置P3的下游。
图9示出所建议的提供单元PRVS的有利实施方式的示意图。提供单元PRVS可以有利地包括计算单元CU、存储单元MU和/或接口IF。在这种情况下,提供单元PRVS、尤其提供单元PRVS的部件可以设计用于实施所建议的用于提供PROV-CD比较数据集CD的方法的各个单独的步骤。尤其地,接口IF可以设计用于提供PROV-D1、PROV-D2时间分辨的第一和第二数据集并且用于提供PROV-CD比较数据集。此外,计算单元CU和/或存储单元MU可以设计用于在空间上配准SREG-D1-D2并且在时间上配准TREG-D1-D2第一数据集D1和第二数据集D2以及用于识别ID-DIFF第一和第二造影剂流之间的偏差。
图10示例性地为所建议的医学成像设备示出医用C形臂X射线设备37的示意图。在此,医用C形臂X射线设备37有利地包括探测器34、尤其X射线探测器和X射线源33。医用C形臂X射线设备37可以有利地设计用于记录和/或提供第一数据集D1和第二数据集D2。尤其地,医用C形臂X射线设备可以设计用于记录包括检查对象31的多个第一投影图像的第一图像数据ID1并且用于记录包括检查对象31的多个第二投影图像的第二图像数据ID2。
为了记录第一和第二投影图像,C形臂X射线设备37的臂38围绕一个或多个轴线可移动地支承。此外,医用C形臂X射线设备37可以包括能够实现C形臂X射线设备37在空间中移动的移动装置39。此外,用于记录布置在患者支承装置32上的检查对象31的第一和第二投影图像的提供单元PRVS可以向X射线源33发送信号24。随后X射线源33可以发射X射线束。当X射线束在与检查对象31相互作用之后撞击在探测器34的表面上时,探测器34可以向提供单元PRVS发送信号21。提供单元PRVS例如可以借助信号21来接收第一和第二投影图像。
此外,该系统可以具有输入单元42、例如键盘和显示装置41、例如监视器和/或显示器。输入单元42可以优选地集成到显示装置41中,例如在电容和/或电阻输入显示器的情况下。在这种情况下,通过医疗操作者在输入单元42处的输入可以实现对医用C形臂X射线设备37的控制、尤其可以实现所建议的用于提供PROV-CD比较数据集CD的方法。为此,输入单元42例如可以向提供单元PRVS发送信号26。
显示装置41可以有利地设计用于尤其叠加或并排地(英文:side-by-side)显示配准的第一数据集D1.REG和/或配准的第二数据集D2.REG和/或比较数据集CD、尤其至少一个表征偏差的参数P.DIFF的图形表示。为此,提供单元PRVS可以向显示装置41发送信号25。此外,用于显示图形表示的显示装置41的至少一个显示参数可以根据比较数据集CD、尤其至少一个表征偏差的参数P.DIFF例如通过颜色编码来调整。
在所描述的图中包含的示意图不以任何方式描绘比例尺或尺寸比例。
最后还应再一次说明的是,以上详细描述的方法和装置仅是实施例,本领域技术人员能够以不同方式修改这些实施例,只要不脱离本发明的范围即可。此外,不定冠词“一个”不排除相关的特征也可以存在多个。同样地,术语“单元”和“元件”不排除相关的部件由多个共同起作用的子部件构成,这些子部件必要时也可以在空间上分散。

Claims (12)

1.一种用于提供(PROV-CD)比较数据集(CD)的方法,所述方法包括:
-提供(PROV-D1)时间分辨的第一数据集(D1),所述第一数据集映射在第一时间段内检查对象(31)的感兴趣区域中的第一造影剂流;
-提供(PROV-D2)时间分辨的第二数据集(D2),所述第二数据集映射在第一时间段之后的第二时间段内检查对象(31)的感兴趣区域中的第二造影剂流;
-在空间上配准(SREG-D1-D2)第一数据集(D1)和第二数据集(D2);
-借助第一和/或第二造影剂流的映射的灌注的方向,识别(ID-V)第一数据集(D1)和第二数据集(D2)中感兴趣区域的至少一个尤其传入的血管区段(V)的映像;
-在时间上配准(TREG-D1-D2)第一数据集(D1)和第二数据集(D2),以便最小化所述至少一个血管区段(V)中第一和第二造影剂流的映射的灌注的时间差;
-通过比较配准的第一数据集(D1.REG)和配准的第二数据集(D2.REG)识别(ID-DIFF)第一和第二造影剂流之间的偏差;
-基于配准的第一数据集(D1.REG)和配准的第二数据集(D2.REG)提供(PROV-CD)所述比较数据集(CD);
其中,所述比较数据集(CD)具有至少一个表征所述偏差的参数(P.DIFF)。
2.根据权利要求1所述的方法,其中,所述第一数据集(D1)的提供(PROV-D1)包括接收(REC-ID1)预先获取的第一图像数据(ID1),
其中,所述第一数据集(D1)由第一图像数据(ID1)重建(RECO-D1),
其中,所述第二数据集(D2)的提供(PROV-D2)包括接收(REC-ID2)预先获取的第二图像数据(ID2),
其中,所述第二数据集(D2)由第二图像数据(ID2)重建(RECO-D2)。
3.根据权利要求2所述的方法,
其中,所述第一图像数据(ID)具有检查对象(31)的从至少部分不同的第一投影方向的多个第一投影图像,
其中,时间分辨的第一数据集(D1)由所述第一投影图像重建,
其中,所述第二图像数据(ID2)具有检查对象(31)的从至少部分不同的第二投影方向的多个第二投影图像,
其中,时间分辨的第二数据集(D2)由所述第二投影图像重建。
4.根据权利要求2或3所述的方法,
其中,所述第一图像数据(ID1)和/或第二图像数据(ID2)具有至少一个掩模图像(MI1,MI2),
其中,所述第一图像数据(ID1)具有多个第一填充图像(FI1),
其中,所述第一数据集(D1)的重建(RECO-D1)包括从多个第一填充图像(FI1)中减去(DIFF-MI1-FI1)至少一个掩模图像(MI1,MI2),
其中,所述第二图像数据(ID2)具有多个第二填充图像(FI2),
其中,所述第二数据集(D2)的重建(RECO-D2)包括从多个第二填充图像(FI2)中减去(DIFF-MI2-FI2)至少一个掩模图像(MI1,MI2)。
5.根据前述权利要求中任一项所述的方法,
其中,配准的第一数据集(D1.REG)和配准的第二数据集(D2.REG)具有多个像点,这些像点分别具有时间强度曲线,
其中,对第一和第二造影剂流之间的偏差的识别包括对配准的第一数据集(D1.REG)和配准的第二数据集(D2.REG)的在空间上彼此对应的像点的时间强度曲线的比较,尤其是对时间强度曲线的斜率、方差、平均强度值、最大强度值和/或累积强度值的比较,
其中,借助对配准的第一数据集(D1.REG)和配准的第二数据集(D2.REG)的在空间上彼此对应的像点的时间强度曲线的比较来确定所述至少一个表征偏差的参数(P.DIFF),所述参数描述填充延迟和/或流量比率。
6.根据权利要求5所述的方法,
其中,在配准的第一数据集(D1.REG)和配准的第二数据集(D2.REG)中分别确定至少一个第一像点和至少一个第二像点,这些第一像点和第二像点映射所述至少一个血管区段(V)中的两个不同的空间位置,
其中,由至少一个第二像点映射的空间位置相对于由至少一个第一像点映射的空间位置布置在下游,
其中,将造影剂团到达时间分别确定为第一和第二造影剂流在由至少一个第二像点映射的空间位置处灌注的时间点,
其中,基于至少一个第二像点的造影剂团到达时间以及对配准的第一数据集(D1.REG)和配准的第二数据集(D2.REG)的至少一个第一像点的时间强度曲线的比较来识别所述偏差(ID-DIFF)。
7.根据权利要求5所述的方法,
其中,在配准的第一数据集(D1.REG)和配准的第二数据集(D2.REG)中分别确定至少一个第一像点和至少一个第二像点,这些像点映射所述至少一个血管区段(V)中的两个不同的空间位置,
其中,对第一和第二造影剂流之间的偏差(ID-DIFF)的识别包括:
-确定配准的第一数据集(D1.REG)和配准的第二数据集(D2.REG)中的至少一个第一像点的时间强度曲线之间的第一比率;
-确定配准的第一数据集(D1.REG)和配准的第二数据集(D2.REG)中的至少一个第二像点的时间强度曲线之间的第二比率;
-确定所述第一比率和第二比率之间的第三比率,
其中,结合所述第三比率(ID-DIFF)识别偏差(DIFF)。
8.根据权利要求6或7所述的方法,
其中,所述感兴趣区域具有血管畸形和/或狭窄和/或动脉瘤(AV),
其中,所述至少一个血管区段(V)相对于血管畸形和/或狭窄和/或动脉瘤(AV)是传入的和/或传出的,
其中,分别相对于血管畸形和/或狭窄和/或动脉瘤(AV),所述至少一个第一像点映射内部或近端的空间位置,并且所述至少一个第二像点映射远端的空间位置。
或者其中,分别相对于血管畸形和/或狭窄和/或动脉瘤(AV),所述至少一个第一像点映射近端的空间位置,并且所述至少一个第二像点映射内部的空间位置。
9.根据权利要求5至8中任一项所述的方法,
其中,分别接收(REC-IP)第一和第二造影剂流的注射参数(IP),
其中,多个时间强度曲线基于注射参数被标准化(NORM-IC),
其中,对第一和第二造影剂流之间的偏差(DIFF)的识别(ID-DIFF)包括对配准的第一数据集(D1.REG)和配准的第二数据集(D2.REG)的在空间上彼此对应的像点的标准化的时间强度曲线的比较。
10.一种提供单元(PRVS),所述提供单元设计用于实施根据前述权利要求中任一项所述的方法。
11.一种医学成像设备,具有根据权利要求10所述的提供单元(PRVS),其中,所述医学成像设备设计用于记录和/或提供(PROV-D1,PROV-D2)第一数据集(D1)和第二数据集(D2)。
12.一种计算机程序产品,具有能够直接加载到提供单元(PRVS)的存储器中的计算机程序,所述计算机程序具有多个程序段,以便当这些程序段由所述提供单元(PRVS)执行时实施根据权利要求1至9中任一项所述的方法的所有步骤。
CN202211175577.7A 2021-09-29 2022-09-26 提供比较数据集 Pending CN115880216A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102021210879.2A DE102021210879A1 (de) 2021-09-29 2021-09-29 Bereitstellen eines Vergleichsdatensatzes
DE102021210879.2 2021-09-29

Publications (1)

Publication Number Publication Date
CN115880216A true CN115880216A (zh) 2023-03-31

Family

ID=85477069

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211175577.7A Pending CN115880216A (zh) 2021-09-29 2022-09-26 提供比较数据集

Country Status (3)

Country Link
US (1) US20230102646A1 (zh)
CN (1) CN115880216A (zh)
DE (1) DE102021210879A1 (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102022208655B3 (de) 2022-08-22 2023-10-12 Siemens Healthcare Gmbh Räumlich-zeitliche Fusion von zeitaufgelösten angiographischen Datensätzen

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5543976B2 (ja) 2008-11-26 2014-07-09 コーニンクレッカ フィリップス エヌ ヴェ 冠状動脈ツリーの視覚化
DE102011081730B4 (de) 2011-08-29 2017-03-23 Siemens Healthcare Gmbh Verfahren zum Vergleich von Flussgeschwindigkeiten

Also Published As

Publication number Publication date
DE102021210879A1 (de) 2023-03-30
US20230102646A1 (en) 2023-03-30

Similar Documents

Publication Publication Date Title
US9867584B2 (en) Method of determining the blood flow through coronary arteries
US9633454B2 (en) Fluid-dynamic analysis of a vascular tree using angiography
US9754390B2 (en) Reconstruction of time-varying data
US9788807B2 (en) Processing apparatus for processing cardiac data
US7853061B2 (en) System and method to improve visibility of an object in an imaged subject
US8285021B2 (en) Three-dimensional (3D) reconstruction of the left atrium and pulmonary veins
US8611989B2 (en) Multi-planar reconstruction lumen imaging method and apparatus
US20110150309A1 (en) Method and system for managing imaging data, and associated devices and compounds
US20090281418A1 (en) Determining tissue surrounding an object being inserted into a patient
US10052032B2 (en) Stenosis therapy planning
WO2016145010A1 (en) System and method for time-resolved, three-dimensional angiography with flow information
JP2006051359A (ja) 医療的侵襲処置計画、並びに侵襲型器具の位置決定及び進路誘導のための方法並びに装置
US8755583B2 (en) Method and apparatus for generating a functional data set of a perfused region of the human or animal body
JP2019512295A (ja) 撮像方法およびシステム
CN114126491B (zh) 血管造影图像中的冠状动脉钙化的评估
US20220051401A1 (en) Providing a scene with synthetic contrast
CN115880216A (zh) 提供比较数据集
US20190076105A1 (en) Hemodynamic parameters for co-registration
Zhang et al. Dynamic estimation of three‐dimensional cerebrovascular deformation from rotational angiography
US20230172571A1 (en) Providing a result data set
US20220262000A1 (en) Vessel registration using functional information
US20230410307A1 (en) Method and system for visualization
US20190216417A1 (en) System and Method for Non-Invasive, Quantitative Measurements of Blood Flow Parameters in Vascular Networks
Zhang Recovery of cerebrovascular morphodynamics from time-resolved rotational angiography

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