CN117617997A - 时间分辨的血管造影数据集的时空融合 - Google Patents
时间分辨的血管造影数据集的时空融合 Download PDFInfo
- Publication number
- CN117617997A CN117617997A CN202311059982.7A CN202311059982A CN117617997A CN 117617997 A CN117617997 A CN 117617997A CN 202311059982 A CN202311059982 A CN 202311059982A CN 117617997 A CN117617997 A CN 117617997A
- Authority
- CN
- China
- Prior art keywords
- time
- reconstruction
- contrast agent
- resolved
- angiographic
- 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
Links
- 230000004927 fusion Effects 0.000 title claims abstract description 23
- 239000002872 contrast media Substances 0.000 claims abstract description 125
- 238000000034 method Methods 0.000 claims abstract description 35
- 238000002583 angiography Methods 0.000 claims abstract description 30
- 230000001360 synchronised effect Effects 0.000 claims abstract description 20
- 230000002792 vascular Effects 0.000 claims description 73
- 230000005484 gravity Effects 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 6
- 208000022211 Arteriovenous Malformations Diseases 0.000 description 11
- 230000005744 arteriovenous malformation Effects 0.000 description 11
- 210000001367 artery Anatomy 0.000 description 10
- 210000004204 blood vessel Anatomy 0.000 description 10
- 238000002347 injection Methods 0.000 description 7
- 239000007924 injection Substances 0.000 description 7
- 230000002123 temporal effect Effects 0.000 description 7
- 238000010521 absorption reaction Methods 0.000 description 4
- 239000008280 blood Substances 0.000 description 4
- 210000004369 blood Anatomy 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 210000004556 brain Anatomy 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 206010003226 Arteriovenous fistula Diseases 0.000 description 2
- 210000003484 anatomy Anatomy 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007170 pathology Effects 0.000 description 2
- 208000034710 Cerebral arteriovenous malformation Diseases 0.000 description 1
- 208000002263 Intracranial Arteriovenous Malformations Diseases 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 210000000709 aorta Anatomy 0.000 description 1
- 201000000034 arteriovenous malformations of the brain Diseases 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000003073 embolic effect Effects 0.000 description 1
- 238000010230 functional analysis Methods 0.000 description 1
- 230000000004 hemodynamic effect Effects 0.000 description 1
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000002673 radiosurgery Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
- 208000019553 vascular disease Diseases 0.000 description 1
- 230000006439 vascular pathology Effects 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/0035—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
- A61B5/0042—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the brain
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/481—Diagnostic techniques involving the use of contrast agents
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
- G06T7/0014—Biomedical image inspection using an image reference approach
- G06T7/0016—Biomedical image inspection using an image reference approach involving temporal comparison
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/60—Type of objects
- G06V20/69—Microscopic objects, e.g. biological cells or cellular parts
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/20—ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/04—Indexing scheme for image data processing or generation, in general involving 3D image data
-
- 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/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- 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
- G06T2207/30048—Heart; Cardiac
-
- 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
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
- G06T2207/30104—Vascular flow; Blood flow; Perfusion
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Optics & Photonics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Vascular Medicine (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Quality & Reliability (AREA)
- Computer Graphics (AREA)
- Physiology (AREA)
- Multimedia (AREA)
- Neurology (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
血管造影记录应该更有信服力。为此,提出了用于时间分辨的血管造影数据集的时空融合的方法。从造影剂在不同位置记录的血管造影的三维数据集中获得各自的4D重建(1,2)。在两个4D重建(1,2)中确定共同血管区域(11)。对于每个造影剂团注,在共同血管区域(11)中确定相应的时间点或时间曲线(23)。最后,两个4D重建(1,2)被同步和融合。
Description
技术领域
本发明涉及一种用于时间分辨的血管造影数据集的时空融合的方法。此外,本发明涉及一种相应的图像处理装置,以及一种血管造影装置。此外,本发明涉及一种计算机程序和一种电子可读的数据载体。
背景技术
所谓的AVM(动静脉畸形)通常是患者大脑中的血管团,并且可以通过几条大动脉供血。AVF(动静脉瘘)应理解为动脉和静脉之间的“短路”,往往很难定位和发现。此外,由于所谓的分路或分流,静脉出流可以以组合的方式进行,这意味着不是只有一个出流,而是多个出流。
为了更好地了解具有两个入流的AVM中的流动条件,目前正在分别研究入流及其对整个血管系统的影响。例如,获取两个单独的3D-DSA图像或4D-DSA图像(DSA:数字减影血管造影)或类似图像。对于第一图像,在第一入流中添加造影剂,对于第二图像,在第二入流中添加造影剂。然后两张照片在空间上融合。
3D或4D-DSA图像的这种融合在K.-K.chen的文章“,,Application of timeresolved 3D digital subtraction angiography to plan cerebral arteriovenousmalformation radio surgery”,26.Januar 2017,10.3174/ajnr.A5074.”中示出。
发明内容
因此,本发明要解决的技术问题在于改进复杂血管系统的显示。
根据本发明,该问题通过一种用于时间分辨的血管造影数据集的时空融合的方法、图像处理装置和血管造影装置来解决。此外,提供了相应的计算机程序和数据载体。本发明的有益的改进设计也由本发明给出。
为了评估例如AVM的整个病理,以组合形式包含所有相关4D流信息(移动3D图像或3D视频)的图像将有助于理解AVM架构,从而促进治疗规划。例如,可以使用基于X射线的医学成像模式,例如X射线血管造影系统或CT扫描仪来获取4D流数据。核磁共振成像(磁共振成像)也允许记录这样的4D图像。此外,在治疗期间,所有3D和/或4D图像的所有相关结构的组合可视化是有帮助的。
根据本发明,提供了一种用于时间分辨的血管造影数据集的时空融合的方法。因此,形成两个或两个以上的时间分辨的血管造影数据集,这些数据集不仅在空间上而且在时间上被融合。
首先,提供在第一位置接受第一造影剂的患者的第一时间分辨的3D血管造影数据集。所述提供可以这样实现,其可以通过在计算机的存储单元中、通过云或互联网或通过其他数据载体提供。或者,提供也可以由血管造影单元直接完成,该血管造影单元获取或获得患者各自的时间分辨三维血管造影数据集。该初始数据集应代表患者在初始位置接受造影剂后的记录,例如AVM的第一入流。然后,造影剂从第一入流开始分布在AVM或血管团中。
同样,提供在与第一位置不同的第二位置接受第二造影剂的该患者的第二时间分辨血管造影3D数据集。与提供第一数据集一样,提供第二数据集也可以以上述各种方式进行。第二造影剂的给药发生在不同于第一位置的第二位置,例如对AVM的第二入流。这使得分析这两种入流的影响成为可能。
然后从第一时间分辨的三维血管造影数据集获得第一4D重建。从连续拍摄的三维图像(三个空间维度)中获得三个空间维度和一个时间维度的4D重建。例如,这种4D重建是空间视频记录,显示造影剂在被检查的血管系统中的扩散。类似地,从第二时间分辨的3D血管造影数据集获得第二4D重建。例如,这提供了两个关于两种不同的造影剂进入血管系统的视频记录。
在进一步的步骤中,确定由两个4D重建所代表的共同血管区域,其中血管区域来自患者或是其血管系统的一部分。因此,在4D重建中确定了这个共同的血管区域,因此在两个图像中都显示了这个区域。例如,如果AVM具有两个入流,并且造影剂是在两个不同的入流中进行的,则共同血管区域代表由两个入流共同供应的血管部位。
然后,从第一4D重建确定与第一造影剂的第一造影剂团注在共同血管区域中的出现有关的第一时间点或第一时间曲线(或称为时间历程)。在第一造影剂给药时注入的造影剂量,即造影剂团注(简称团注),从第一位置输送到共同的血管区域。这种运输以一定的速度进行,其中,在血管中流动的血液稀释了造影剂。因此,与注射一起给药的造影剂以一定的时间曲线(时间-对比度曲线)出现在下游血管区域。造影剂浓度在注射后逐渐增加,在断层成像中更久之后再次减少。造影剂浓度在血管点或血管区域的这种随时间的变化被称为造影剂团注在(共同)血管区域的出现。关于这种“出现”,可以确定造影剂浓度的一个时间点或整个时间曲线。例如,确定的时间可以是团注到达共同血管区域的时间。
类似地,从第二4D重建中确定与第二造影剂给药的第二造影剂团注在共同血管区域中的出现有关的第二时间点或第二时间曲线。例如,在这里可以确定第二团注到达共同血管区域的时间。
在此之后,可以根据第一和第二时间点或时间曲线来同步第一4D重建和第二4D重建。例如,如果第一和第二时间点分别是各自团注到达共同血管区域的时间点,则这两个到达时间点可用于同步4D重建或视频。因此,例如,在同步之后,两个造影剂团注在同一时间到达共同的血管区域。
最后,将两个同步的4D重建进行融合。融合意味着4D重建在空间上被正确地叠加。因为4D重建也彼此同步,所以融合相当于4D融合。这意味着在融合的4D重建或融合的视频中,不仅发生空间融合,而且还发生单个的重建或视频的时间融合。因此,换句话说,提出了一种在空间和时间上融合两个(或更多)血管造影4D数据集以产生单个时间上同步的组合的4D数据集的方法。
例如,由于所依据的血流动力学,可以通过识别在两个(所有)4D数据集中对比度增强的血管部位(即单个体素或体素的部位)来在时间上配准4D血管造影图像。在两个(所有)4D图像中分别识别这些共同区域被对比度增强的时间点,实现数据集的时间同步。
在根据本发明的方法的一个实施例中规定,第一和第二时间点分别对应于各造影剂团注在共同血管区域中的到达时间点。血管区域可以选择较大或较小。例如,血管区域可以包括一团血管或单个血管。各自的到达地点由图像处理监控。例如,当造影剂到达血管区域时,在血管区域中可以检测到黑化。例如,就黑化而言,可以监测单个血管,或者在更大的血管区域内形成平均值。如上所述,造影剂在血管区域的浓度最初持续增加并达到最大值。当造影剂流出时,浓度再次连续下降。例如,如果浓度或黑化超过检测阈值,这可以定义为造影剂团注的到达时间。这样,对于第一造影剂的团注和第二造影剂的团注,可以分别确定到达时间。然后可以将由此获得的第一时间点和第二时间点等同以进行同步。
在本发明方法的替代设计中规定,第一和第二时间点分别对应于各造影剂团注在共同血管区域中的最大浓度的时间点。由于造影剂浓度具有特定的时间曲线,因此可以从该时间曲线中选择具有该过程代表性的定义时间点。在本例中,使用浓度最大值来确定浓度历程的代表性时间。这样的最大值可以通过图像处理可靠地确定。
在本发明方法的另一备选方案中规定,第一和第二时间点分别对应于各自造影剂团注的各自重心在共同血管区域中的驻留时间点。例如,可以从相对较短的造影剂团注的时间曲线中确定相应的造影剂团注的重心。例如,相应的造影剂团注的重心可以定义为造影剂团注在特定血管区域中的相应曲线函数的重心。这种重心也可以通过对造影剂的功能分析来自动确定。
在另一个实施例中,这样的时间点被用于同步,在该时间点第一和第二时间曲线的互相关(Kreuzkorrelation)(分别必要时与造影剂入流和出流或仅入流相关)具有最大值。因此,不是从每个单独的时间曲线中直接确定一个时间点,而是两种造影剂的时间曲线直接受到互相关或类似的影响。因此,使用每个造影剂团注的整个时间曲线来执行同步。形象地说,这两个时间序列是以这样一种方式叠加的,它们具有尽可能大的一致性。因此,对于每个团注,不仅使用单个时间点进行同步,而且使用整个时间曲线。
根据另一个实施例,从4D重建中确定第一和第二造影剂团注都流入的多个血管部位,并且两个造影剂给药的各自造影剂团注首先到达的血管部位被自动地确定为共同血管区域。第一造影剂和第二造影剂的造影剂经常流入几个共同的血管部位。由于造影剂随着时间的推移而变得越来越稀释,因此选择尽可能早的时间进行同步是有利的。因此有利的是,使用第一共同的血管部位用于同步,第一造影剂给药和第二造影剂给药的造影剂都流入所述第一共同的血管部位。在这个第一共同的血管部位中能分别最好地识别、尤其分辨造影剂。
在另一个实施例中规定,至少一个另外的4D重建与第一和第二4D重建同步并融合。例如,不仅两个视频记录被同步和融合,而且三个或更多的视频记录被同步和融合。相关数据集可以通过不同的成像模式(血管造影系统,CT扫描仪等)获得。这样可以更好地分析一个非常复杂的血管系统。
此外,在一个实施例中可以规定,第一和第二4D重建与3D血管造影数据集融合。通过这种方式,例如,来自较小血管区域的两个视频记录可以被投影或嵌入到静态的整体血管结构中。因此,例如,如果主治医生在较大的血管树或整个血管树中看到要治疗的区域,他可以更好地计划治疗。
此外,在一个实施例中可以规定,第一和第二4D重建与三维断层扫描数据集融合。例如,三维断层扫描数据集详细反映了患者的解剖结构。因此,由两个4D重建显示的血管区域可以精确地插入到患者的空间解剖结构中。这种类型的融合成像也有助于规划或治疗医生更好地评估血管情况。
根据本发明,还通过用于时间分辨的血管造影数据集的时空融合的图像处理装置解决上述问题,图像处理装置具有存储单元,其被设计成提供患者的第一时间分辨的3D血管造影数据集,患者在第一位置接受第一造影剂给药,以及提供患者的第二时间分辨的3D血管造影数据集,患者在与第一位置不同的第二位置接受了第二造影剂给药,以及图像处理单元,其设计用于从第一时间分辨的3D血管造影数据集获得第一4D重建和从第二时间分辨的3D血管造影数据集获得第二4D重建,确定由两个4D重建所代表的患者的共同血管区域,从第一4D重建确定与第一造影剂给药的第一造影剂团注在共同血管区域中的出现有关的第一时间点或时间曲线,从第二4D重建确定与第二造影剂给药的第二造影剂团注在共同血管区域中的出现有关的第二时间点或时间曲线,根据第一和第二时间点或时间曲线同步第一4D重建和第二4D重建,并融合两个同步的4D重建。
因此,图像处理设备能够根据上述方法同步和融合4D血管造影重建。上述与本发明方法相关的优点和改进也比照适用于图像处理设备。其中,各个方法步骤可被视为图像处理装置的相应单元的功能特征。
此外,根据本发明提供了一种用于时间分辨的血管造影数据集的时空融合的血管造影装置,具有造影单元,其设计用于在第一位置获得第一造影剂的患者中检测第一时间分辨的3D血管造影数据集,并在与第一位置不同的第二位置获得第二造影剂的患者中检测第二时间分辨的3D血管造影数据集,以及图像处理单元,其设计用于其设计用于从第一时间分辨的3D血管造影数据集获得第一4D重建和从第二时间分辨的3D血管造影数据集获得第二4D重建,确定由两个4D重建所代表的患者的共同血管区域,从第一4D重建确定与第一造影剂给药的第一造影剂团注在共同血管区域中的出现有关的第一时间点或时间曲线,从第二4D重建确定与第二造影剂给药的第二造影剂团注在共同血管区域中的出现有关的第二时间点或时间曲线,根据第一和第二时间点或时间曲线同步第一4D重建和第二4D重建,并融合两个同步的4D重建。
这里的血管造影装置能够自行生成三维血管造影数据集或者说3D血管造影数据集。当然,所述血管造影装置还可以具有用于缓存所获得的三维血管造影数据集的存储装置。但是,此存储单元不必必然存在。
关于血管造影装置的改进和优点,上述关于图像处理装置的内容是适用的。
根据本发明的另一个方面,提供具有指令的计算机程序。当指令由至少一个上述类型的图像处理装置执行时,这些指令使图像处理装置执行本发明的方法。这些指令可以作为程序代码提供。程序代码可以作为二进制代码或汇编程序和/或作为编程语言的源代码,例如C语言,和/或作为程序脚本,例如Python。
此外,这里还提供了存储所述计算机程序的电子可读的数据载体。数据载体可以是DVD,内存条等。
附图说明
现在根据所附附图对本发明进行更详细的说明,其中显示:
图1示出两个4D血管造影数据集的时空融合的示意图;
图2示出本发明方法的实施例的示意框图;
图3示出本发明的血管造影装置的示意图;和
图4示出造影剂在血管(区域)中的时间曲线或时间对比曲线。
具体实施方式
下面更详细地描述的实施例代表了本发明的优选实施例。在附图中,相同或功能相同的元件可以被标记为相同的参考符号。在适当的情况下,不一定就不同的图形重复描述相同或功能相同的元件。
如果在本公开的范围内提及,本发明的图像处理装置或血管造影装置或其图像处理单元的部件被设置、构造、设计或类似地用于执行或实现特定功能、实现特定效果或服务于特定目的,这可以被理解为,除了该部件的原理或理论上的可用性或用于该功能、效果或目的的适用性之外,通过适当的调整、编程、物理构造等,该部件能够具体而实际地执行或实现该功能、实现效果或服务于该目的。
对于一个具体的实施例,可以执行以下步骤:
a)记录具有不同造影剂注射部位(通常是病变血管区域的不同动脉供应,例如在病人的大脑中)的患者的两个或多个对比度增强的三维数据集。
b)血管造影数据集的4D重建(例如,使用基于流体限制的方法,如DE 10 201 7200 489A1中所述)。
c)确定重建4D图像中每个血管的体素(或血管的体素的区域)的时间动态(例如团注到达时间,最大黑化时间等)。
d)识别在两个(或多个)血管造影4D图像中对比度增强的血管体素。
e)基于这些共同的血管体素在时间上同步4D图像(这里也称为4D重建)。团注到达时间(或其他时间度量,例如最大黑化时间)定义用于时间同步的“锚时间点间”。
图1示出了用于分别说明两个血管造影4D数据集或4D重建1和2的时空融合的工作流程。第一4D重建是从(在第一位置接受第一造影剂的)患者的第一时间分辨的3D血管造影数据集中获得的。第二4D重建是从在与第一位置不同的第二位置接受第二造影剂的该患者的第二时间分辨的3D血管造影数据集中获得的。
图1以二维方式显示4D重建1和2。在第一4D重建1所依据的第一造影剂给药中,将造影剂注入第一动脉3。造影剂从该第一动脉3流入血管团5的第一血管部段4。
在第一4D重建1的选定的视图中,可以看到血管结构内造影剂的时间动态。因此,造影剂在注射后约150毫秒到达第一动脉3。大约200-300毫秒后,达到更细的血管结构6。大约400毫秒后,造影剂在静脉出流7中排出。
类似地,第二4D重建2基于向第二动脉8注入造影剂。血液或造影剂从这第二动脉8流入更细的血管结构9,然后再流入静脉出流10。时间刻度指示造影剂何时到达相应的区域。
由于血管团5具有许多动静脉畸形和分流,来自第一动脉3的造影剂与来自第二动脉8的造影剂流入同一血管部位。该血管部位可称为共同血管区域11。因此,在这个共同血管区域11中,造影剂在第一造影剂给药中进入第一动脉3,在第二造影剂给药中进入第二动脉8。由于这个原因,在4D重建1和2中都显示了这个共同血管区域11。
由于共同区域11在两个4D重建1和2中都被再现,因此两个4D重建1和2可以在空间上彼此配准。这意味着两个4D重建1和2在空间上对齐并相互表示。这样的空间配准导致以不同颜色表示第一4D重建1和第二4D重建2的血管部段的空间配准图像12。受两种造影剂影响的血管部位呈现混合着色(例如,如果第一造影剂的血管结构变为红色且第二造影剂的血管结构变为蓝色,则为紫色)。
最后,对两个4D重建1和2进行时间配准。由此得到时间配准图像13,其简化地再现了两个同步的4D重建1和2的融合。在时间配准图像13中,与在空间配准图像12中一样,可以识别两个造影剂剂量或两个4D重建1和2的所有血管部段。然而,对于融合的4D数据集,造影剂流现在是同步的,即在时间上是同步的。在本实施例中,时间点被设置为造影剂到达共同的区域11的同步时间。因此,融合导致了一个三维视频记录,其基于两个彼此在空间和时间上被配准的三维视频记录。
图2在流程图中示出了根据本发明的用于时间分辨血管造影数据集的时空融合的方法的示例的更详细的流程。在第一步骤S1中,提供或检测或记录患者的第一时间分辨的3D血管造影数据集。为此,患者在第一位置接受了第一造影剂。在进一步的步骤S2中,提供或记录该患者的第二时间分辨的3D血管造影数据集,其中患者在不同于第一位置的第二位置接受了第二造影剂。
在步骤S3中,从第一时间分辨的三维血管造影数据集获得第一4D重建。同样,在步骤S4中,从第二时间分辨的3D血管造影数据集获得第二4D重建或第二4D数据集。
在步骤S5中,确定由两个4D重建或数据集分别表示的共同血管区域。患者的这个共同血管区域对于随后在空间和时间上相互登记4D数据集是很重要的。为此,在步骤S6中,从第一4D重建中确定与共同血管区域中的第一造影剂的第一造影剂团注的出现有关的第一时间点或第一时间曲线。以同样的方式,从第二4D重建中确定与第二造影剂给药的第二造影剂团注在共同血管区域中出现有关的第二时间点或第二时间曲线。
在另一步骤S8中,基于第一时间点和第二时间点或第一时间曲线第二时间曲线同步两个4D数据集,即第一4D重建和第二4D重建。最后,在步骤S9中,两个同步的4D重建被融合。通过这种方式,例如,在一个共同的视图中获得两个时间上彼此一致的3D视频记录。
图3示意性地示出具有血管造影单元15和图像处理单元16的本发明的血管造影装置14的实施例。所述血管造影单元可以具有C形臂17,其带有X射线发射器18和相对的X射线检测器19。此外,血管造影单元15可以具有控制单元20,用于控制记录操作,以便在必要时分别获取和检测相应的时间分辨的3D血管造影数据集。可选地,所述血管造影装置具有存储单元21,以存储所检测的时间分辨的3D血管造影数据集,并且能够使其可用于随后的图像处理。
血管造影装置14的图像处理单元16在某些情况下也可实现为具有存储单元21和/或数据接口的独立设备,通过该存储单元21和/或数据接口可提供相应的3D数据集。此外,图像处理单元16可以具有计算单元22,利用该计算单元22分别从3D数据集中获得4D重建和4D数据集。此外,利用计算单元22可以确定共同血管区域以及在该共同血管区域中第一造影剂团注的第一时间点或时间曲线和第二造影剂团注的第二时间点或时间曲线。最后,计算单元22还被设计成分别基于时间点和时间曲线同步两个4D重建,并且融合两个同步的4D重建,并且在适当的情况下在屏幕27上表示。
下面描述了一些细节,这些细节可以用来进一步改进所描述的方法和设备,或者在适当的情况下应该考虑这些细节。共同的血管区域可以由一个血管体素或一个血管体素区域来描述,并且可以位于静脉侧或动脉侧,这取决于两个注射部位和患者的血管病理。理想情况下,早期增强的共同体素(或共同血管体素的区域)用于时间上的配准以将不准确性最小化,因为在后期增强的体素中,由于非增强血液的交叉流动以及造影剂团注的分散,造影剂增强发生减弱。
通常,两个(或多个)4D血管造影数据集已经在空间上配准或者说登记,因为它们通常是在同一过程中检测的,并且患者在检测之间没有移动或几乎没有移动。但是,如有必要,必须在时间配准或同步之前实施一次或多次空间配准。空间配准的方法是从现有技术中已知的。
在一个特殊实施例中—如上所述—可以相应地融合两个以上的4D血管造影数据集。在另一实施例中,至少两个4D血管造影数据集与至少一个3D血管造影数据集(例如静态CTA或MRA图像)融合。在另一实施例中,至少两个4D血管造影数据集与至少一个3D断层数据集(例如静态CT或MR图像)融合。
根据本发明的方法可用于脑外和脑内的血管病理。例如,该方法还可以用于分析四肢或面部的末梢AVM。
在进一步的实施例中,用于两个或多个4D血管造影数据集的时间同步的时间参数可以变化。这些变化的可能性从图4中可见,图4显示了血管中造影剂浓度的时间曲线,即时间-对比度曲线。如果需要,时间-对比度曲线的一部分也用于同步。X射线介质浓度通常与X射线衰减或吸收系数μ成正比。例如,在时间t=0时进行造影剂给药。然后,在血管树中的指定位置发生洪水阶段。黑化或吸收系数在时间BAT(团注到达时间)达到阈值24。例如,该阈值24是检测阈值。这意味着只有从检测阈值24开始,造影剂才能被明确地检测出来。在团注到达时间BAT之后,造影剂浓度增加,因此吸收系数μ根据时间曲线23增加到最大值或峰值25。在峰值时间PT达到该峰值25。随后,特定血管部位的吸收系数再次连续下降。造影剂浓度的整个时间曲线23具有重心26,其通过重心时间CGT表征。
因此,给出了一些时间参数,可用于两个或多个4D血管造影数据集的时间同步,即相对于两个或多个造影剂剂量:
a)相对应的团注到达时间
b)最大造影剂浓度下的相对应的峰值时间PT
c)相应时间-对比度曲线的重心的相对应的重心时间CGT 23
d)同步时间点的选择,时间-对比度曲线23的归一化的互相关(Kreuzkorrelation)在所述同步时间点达到最大值
e)以及由时间-对比度曲线23产生的其他时间参数也可用于同步。
上述实施例基于根据具有不同造影剂注入位置的多个4D图像生成组合4D血管图的方法。这使得能够更完整地绘制供应器官(例如大脑)或血管疾病(例如AVM)的整个血管系统的图像。在AVM手术的情况下,这种时空融合的4D图像可用于基于共同数据集规划所有可能的动静脉栓塞途径。
Claims (13)
1.一种用于时间分辨的血管造影数据集的时空融合的方法,
其特征在于,
-提供(S1)患者的第一时间分辨的3D血管造影数据集,所述患者在第一位置接受第一造影剂给药,
-提供(S2)所述患者的第二时间分辨的3D血管造影数据集,所述患者在与第一位置不同的第二位置接受第二造影剂给药,
-从第一时间分辨的3D血管造影数据集获得(S3)第一4D重建(1),
-从第二时间分辨的3D血管造影数据集获得(S4)第二4D重建(2)的,
-确定(S5)两个4D重建所代表的患者的共同血管区域(11),
-从第一4D重建(1)确定(S6)与第一造影剂给药的第一造影剂团注在共同血管区域(11)中的出现有关的第一时间点或第一时间曲线(23),
-从第二4D重建(2)确定(S7)与第二造影剂给药的第二造影剂团注在共同血管区域(11)中的出现有关的第二时间点或第二时间曲线(23),
-根据第一时间点和第二时间点或第一时间曲线和第二时间曲线(23)同步(S8)第一4D重建(1)和第二4D重建(2);和
-融合两个同步的4D重建(1,2)。
2.根据权利要求1所述的方法,其特征在于,第一时间点和第二时间点分别对应于各造影剂团注在共同血管区域(11)中的到达时间点(BAT)。
3.根据权利要求1所述的方法,其特征在于,第一时间点和第二时间点分别对应于各造影剂团注在共同血管区域(11)中的最大浓度(25)的时间点(PT)。
4.根据权利要求1所述的方法,其特征在于,第一时间点和第二时间点分别对应于各造影剂团注的各自重心(26)在共同血管区域(11)中的驻留时间点(CGT)。
5.根据权利要求1所述的方法,其特征在于,将第一时间曲线和第二时间曲线的互相关具有最大值的时间点用于同步。
6.根据前述权利要求之一所述的方法,其特征在于,从4D重建(1,2)中确定第一造影剂团注和第二造影剂团注都流入其中的多个血管部位,并且两个造影剂给药的各自造影剂团注首先到达的血管部位被自动地确定为共同血管区域(11)。
7.根据前述权利要求之一所述的方法,其特征在于,至少另一个4D重建与第一4D重建和第二4D重建同步并融合。
8.根据前述权利要求之一所述的方法,其特征在于,第一4D重建和第二4D重建(1,2)与3D血管造影数据集融合。
9.根据前述权利要求之一的方法,其特征在于,第一4D重建和第二4D重建(1,2)与三维断层扫描数据集融合。
10.一种用于时间分辨的血管造影数据集的时空融合的图像处理装置,
其特征在于:
-存储单元(21),其被设计用于
o提供患者的第一时间分辨的3D血管造影数据集,所述患者在第一位置接受第一造影剂给药,以及
○提供所述患者的第二时间分辨的3D血管造影数据集,所述患者在与第一位置不同的第二位置接受第二造影剂给药,和
-图像处理单元(16),其设计用于
○从第一时间分辨的3D血管造影数据集获得第一4D重建和从第二时间分辨的3D血管造影数据集获得第二4D重建,
○确定由两个4D重建所代表的患者的共同血管区域,
○从第一4D重建(1)确定与第一造影剂给药的第一造影剂团注在共同血管区域中的出现有关的第一时间点或第一时间曲线,
○从第二4D重建(2)确定与第二造影剂给药的第二造影剂团注在共同血管区域中的出现有关的第二时间点或第二时间曲线,
○根据第一和第二时间点或时间曲线同步第一4D重建(1)和第二4D重建(2),并且
○融合两个同步的4D重建(1,2)。
11.一种用于时间分辨的血管造影数据集的时空融合的血管造影装置,
其特征在于:
-造影单元(15),其设计用于
○检测患者的第一时间分辨的3D造影数据集,所述患者在第一位置接受第一造影剂给药,并且
o检测所述患者的第二时间分辨的3D造影数据集,所述患者在与所述第一位置不同的第二位置接受第二造影剂给药,和
-图像处理单元(16),其设计用于
○从第一时间分辨的3D血管造影数据集获得第一4D重建(1)和从第二时间分辨的3D血管造影数据集获得第二4D重建(2),
○确定由两个4D重建所代表的患者的共同血管区域(11),
○从第一4D重建(1)确定与第一造影剂给药的第一造影剂团注在共同血管区域中的出现有关的第一时间点或第一时间曲线(23),
o从第二4D重建(2)确定与第二造影剂给药的第二造影剂团注在共同血管区域中的出现有关的第二时间点或第二时间曲线(23),
o根据第一时间点和第二时间点或第一时间曲线和第二时间曲线(23)同步第一4D重建(1)和第二4D重建(2),并且
o融合两个同步的4D重建(1,2)。
12.一种计算机程序,当计算机程序在权利要求10的图像处理装置上执行时,执行根据权利要求1至9中的任一项所述的方法的步骤。
13.一种电子可读的数据载体,其中存储了根据权利要求12所述的计算机程序。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102022208655.4A DE102022208655B3 (de) | 2022-08-22 | 2022-08-22 | Räumlich-zeitliche Fusion von zeitaufgelösten angiographischen Datensätzen |
DE102022208655.4 | 2022-08-22 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117617997A true CN117617997A (zh) | 2024-03-01 |
Family
ID=88094436
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311059982.7A Pending CN117617997A (zh) | 2022-08-22 | 2023-08-21 | 时间分辨的血管造影数据集的时空融合 |
Country Status (3)
Country | Link |
---|---|
US (1) | US20240062340A1 (zh) |
CN (1) | CN117617997A (zh) |
DE (1) | DE102022208655B3 (zh) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102223843B (zh) | 2008-11-26 | 2013-10-23 | 皇家飞利浦电子股份有限公司 | 冠状动脉树的可视化 |
WO2015044433A1 (en) | 2013-09-27 | 2015-04-02 | Koninklijke Philips N.V. | Merging vessel maps |
DE102017200489A1 (de) | 2017-01-13 | 2018-07-19 | Siemens Healthcare Gmbh | Verfahren zur Reduzierung von aufgrund von Gefäßüberlagerungen auftretenden Artefakten in einem vierdimensionalen Angiographiedatensatz, Bildverarbeitungseinrichtung, Computerprogramm und elektronisch lesbarer Datenträger |
DE102021210879A1 (de) | 2021-09-29 | 2023-03-30 | Siemens Healthcare Gmbh | Bereitstellen eines Vergleichsdatensatzes |
-
2022
- 2022-08-22 DE DE102022208655.4A patent/DE102022208655B3/de active Active
-
2023
- 2023-08-07 US US18/230,790 patent/US20240062340A1/en active Pending
- 2023-08-21 CN CN202311059982.7A patent/CN117617997A/zh active Pending
Also Published As
Publication number | Publication date |
---|---|
US20240062340A1 (en) | 2024-02-22 |
DE102022208655B3 (de) | 2023-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6695478B2 (ja) | 血管セグメント化 | |
US6073042A (en) | Display of three-dimensional MRA images in which arteries can be distinguished from veins | |
Fain et al. | High-spatial-resolution contrast-enhanced MR angiography of the renal arteries: a prospective comparison with digital subtraction angiography | |
Rofsky et al. | MR angiography in the evaluation of atherosclerotic peripheral vascular disease | |
Heautot et al. | Analysis of cerebrovascular diseases by a new 3-dimensional computerised X-ray angiography system | |
Hyodoh et al. | Usefulness of preoperative detection of artery of Adamkiewicz with dynamic contrast-enhanced MR angiography | |
Mallouhi et al. | Detection and characterization of intracranial aneurysms with MR angiography: comparison of volume-rendering and maximum-intensity-projection algorithms | |
JP6214646B2 (ja) | 血管造影図における時間的な解剖学的標的タグ付け | |
Pitton et al. | MRI versus helical CT for endoleak detection after endovascular aneurysm repair | |
US20190117179A1 (en) | Systems And Methods For Deciding Management Strategy in Acute Ischemic Strokes Using Rotational Angiography | |
US9125616B2 (en) | Collateral blood flow assessment | |
Struffert et al. | Monitoring of balloon test occlusion of the internal carotid artery by parametric color coding and perfusion imaging within the angio suite: first results | |
Sekine et al. | 4D flow MRI assessment of extracranial-intracranial bypass: qualitative and quantitative evaluation of the hemodynamics | |
Ohno et al. | Contrast-enhanced MR perfusion imaging and MR angiography: utility for management of pulmonary arteriovenous malformations for embolotherapy | |
Vogt et al. | Contrast-enhanced time-resolved 4D MRA of congenital heart and vessel anomalies: image quality and diagnostic value compared with 3D MRA | |
Carr et al. | High-resolution breath-hold contrast-enhanced MR angiography of the entire carotid circulation | |
Zhang et al. | Diagnostic accuracy of time-resolved 2D projection MR angiography for symptomatic infrapopliteal arterial occlusive disease | |
CN114502069A (zh) | 加速mrt检查 | |
CN117617997A (zh) | 时间分辨的血管造影数据集的时空融合 | |
Kubo et al. | Thoracoabdominal-aortoiliac MDCT angiography using reduced dose of contrast material | |
Hansch et al. | Quantitative evaluation of MR perfusion imaging using blood pool contrast agent in subjects without pulmonary diseases and in patients with pulmonary embolism | |
Crean | Technical considerations for ACHD imaging | |
US8660333B2 (en) | Functional imaging | |
KR101744424B1 (ko) | Mra 영상을 이용하여 관류 특성 확인용 추가 영상을 얻기 위한 방법 및 시스템 | |
Castillo | Digital subtraction angiography (DSA): basic principles |
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 |