CN104424659A - 图像数据的基于图像的运动补偿 - Google Patents
图像数据的基于图像的运动补偿 Download PDFInfo
- Publication number
- CN104424659A CN104424659A CN201410429770.8A CN201410429770A CN104424659A CN 104424659 A CN104424659 A CN 104424659A CN 201410429770 A CN201410429770 A CN 201410429770A CN 104424659 A CN104424659 A CN 104424659A
- Authority
- CN
- China
- Prior art keywords
- image
- data
- image data
- view data
- sub
- 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
Links
- 230000033001 locomotion Effects 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 claims abstract description 61
- 230000005855 radiation Effects 0.000 claims abstract description 34
- 238000005457 optimization Methods 0.000 claims abstract description 9
- 238000005259 measurement Methods 0.000 claims description 38
- 238000004590 computer program Methods 0.000 claims description 29
- 238000001228 spectrum Methods 0.000 claims description 12
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 4
- 210000004351 coronary vessel Anatomy 0.000 claims description 4
- 238000003860 storage Methods 0.000 claims description 4
- 238000011478 gradient descent method Methods 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 2
- 230000002123 temporal effect Effects 0.000 abstract description 2
- 230000001131 transforming effect Effects 0.000 abstract 1
- 238000002591 computed tomography Methods 0.000 description 15
- 230000006870 function Effects 0.000 description 8
- 239000002872 contrast media Substances 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 238000002360 preparation method Methods 0.000 description 3
- 230000011218 segmentation Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000000747 cardiac effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000010412 perfusion Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000000541 pulsatile effect Effects 0.000 description 2
- 108010022579 ATP dependent 26S protease Proteins 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000002586 coronary angiography Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000008014 freezing Effects 0.000 description 1
- 238000007710 freezing Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- PNDPGZBMCMUPRI-UHFFFAOYSA-N iodine Chemical compound II PNDPGZBMCMUPRI-UHFFFAOYSA-N 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000000007 visual effect Effects 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/503—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 the heart
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- 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]
-
- 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]
-
- 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/5205—Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce 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/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
- A61B6/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
- G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
-
- 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
-
- 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/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4064—Arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
- A61B6/4078—Fan-beams
-
- 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/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4064—Arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
- A61B6/4085—Cone-beams
-
- 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/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
- A61B6/4435—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
- A61B6/4441—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure the rigid structure being a C-arm or U-arm
-
- 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/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy 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/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/507—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 determination of haemodynamic parameters, e.g. perfusion CT
-
- 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/5288—Devices using data or image processing specially adapted for radiation diagnosis involving retrospective matching to a physiological signal
-
- 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/54—Control of apparatus or devices for radiation diagnosis
- A61B6/541—Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal
-
- 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/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- 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/20048—Transform domain processing
-
- 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/20172—Image enhancement details
- G06T2207/20201—Motion blur correction
-
- 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
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/412—Dynamic
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/416—Exact reconstruction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/421—Filtered back projection [FBP]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/436—Limited angle
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- High Energy & Nuclear Physics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Pulmonology (AREA)
- Cardiology (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Architecture (AREA)
- Computer Graphics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Quality & Reliability (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及图像数据运动补偿的方法。图像数据是由测量数据的空间三维重建,在辐射源和检查对象之间相对旋转运动下在角度范围内录取测量数据。角度范围包括多个子角度范围,每个子角度范围对应一个时间控制点。在迭代方法第一步骤期间确定图像数据的图像度量的至少一个值,图像度量是图像数据中运动伪影的量度。基于图像数据根据图像度量能确定图像数据运动补偿的运动场。在第二步骤中确定子图像数据,子图像数据分别对应来自子角度范围的测量数据的空间三维重建。在第三步骤中借助优化方法根据图像度量在控制点上确定图像数据的运动场,在第四步骤中将子图像数据按照运动场描述的运动进行变换。在第五步骤中通过拼合子图像数据产生新的图像数据。
Description
技术领域
本发明涉及一种用于图像数据的运动补偿的方法以及计算机程序、计算机程序产品和计算机断层造影***。
背景技术
计算机断层造影是主要用于医学诊断以及用于材料检查的成像方法。在计算机断层造影中为了录取空间上三维的图像数据,辐射源以及与该辐射源一起工作的辐射探测器围绕检查对象旋转。在旋转运动期间在一个角度范围内录取测量数据。测量数据是多个投影,所述投影包含关于辐射通过检查对象的衰减的信息。测量数据也称为原始数据。从这些测量数据中然后可以重建图像数据,例如借助所谓的滤波反投影或借助迭代的重建方法。如果检查对象在录取期间运动,则在图像数据的重建中会由于运动而形成模糊和伪影。
因此一般公知的是,在断层造影录取期间的心脏运动会导致不一致和通过运动伪影损害图像数据。由此会严重限制图像数据的临床可用性。为了避免运动伪影,在现代的断层造影心脏成像中使用借助回顾性或前瞻性触发的、测量数据的相位相关的录取。因此在前瞻性触发中仅录取在大约心脏的静止阶段的一定的时间间隔中的数据并且用于图像重建。该方案的共同目标是“冻结”心脏运动,以减小数据不一致性和由此提高图像质量。
在部分和周期运动的对象的情况中用于提高图像质量的一种方案从文献DE 10 2011 083 643 A1公知。在那里公开了一种用于确定CT图像数据组的图像空间中的运动场(Bewegungsfeld)的方法。在此对于预先给出的运动阶段和投影角度范围采集CT***的投影数据组。此外通过将利用第一分析重建算法在使用运动场的条件下运动补偿地重建的、第一图像分辨率的CT图像数据组的、基于原始数据的或变换到图像空间中的成本函数最小化,进行运动场的迭代确定。
发明内容
本发明要解决的技术问题是,如何能够基于图像数据对运动的检查对象执行图像数据的运动补偿。
以下描述关于要求保护的装置和关于要求保护的方法的技术问题的按照本发明的解决方案。在此提到的特征、优点或替换实施方式同样也适用于其他要求保护的内容,反之亦然。换言之,例如针对***的具体权利要求也可以利用关于方法描述或要求保护的特征来扩展。方法的相应功能特征在此通过相应的具体模块构造。
本发明涉及一种用于对运动的检查对象的图像数据进行运动补偿的方法。图像数据是由测量数据的空间三维重建,其中在辐射源和检查对象之间相对旋转运动的情况下在角度范围内录取测量数据。该角度范围包括多个子角度范围,其中每个子角度范围对应于一个时间控制点。特别地,时间控制点可以是通过时间间隔定义的时刻,在该时间间隔期间在各自的子角度范围内录取测量数据。本发明的核心是以下的迭代方法。在按照本发明的方法的第一步骤期间确定对于图像数据的图像度量L的至少一个值,其中图像度量是对于图像数据中运动伪影的量度。发明人现在已经认识到,基于图像数据地根据图像度量可以确定用于图像数据的运动补偿的运动场s,其中运动场是按照对于图像数据的运动的量度。为此在第二步骤中确定子图像数据,其中子图像数据分别对应于来自于子角度范围的测量数据的空间三维重建。在第三步骤中借助优化方法根据图像度量在控制点上确定图像数据的运动场s,从而在第四步骤中可以将子图像数据按照通过运动场描述的运动进行变换。在第五步骤中通过将子图像数据拼合产生新的图像数据。迭代方法的这些步骤一直进行,直到达到中断标准。新的图像数据通过被变换的子图像数据的变换被运动补偿并且相对于初始的图像数据具有改善的图像质量。通过按照本发明基于图像地进行该方法,与其他基于原始数据的方法相比,需要更小的计算开销并且可以特别地在所谓的工作站上和独立于重建计算机进行。由此证明按照本发明的方法是对于临床实践特别实用的。
按照本发明的另一方面,角度范围是对于断层造影完整的重建所需的最小角度范围。在这样的情况中,角度范围通过给出,其中说明扇形角度,也就是在径向方向上通过辐射源发射的辐射扇形的开口角度。
按照本发明的另一方面,通过首先进行图像数据到频率域的傅里叶变换,然后借助滤波函数掩模图像数据并且借助傅里叶变换又将掩模的频谱转换到空间域,确定子图像数据。本发明的该方面的优点是,子图像数据的确定可以利用小的计算开销实现并且无需访问原始数据。
按照本发明的另一方面,分别借助图像数据的正投影和对应于子角度范围的一部分投影后的图像数据的然后的反投影来确定子图像数据。本发明的该方面也可以利用小的计算开销进行并且无需访问原始数据。
按照本发明的另一方面,在迭代方法执行期间和/或之后,借助运动场进行图像数据的正投影以及然后的运动补偿的反投影。在该方面中尽可能无失真地进行图像数据的运动补偿,从而实现特别高的图像质量。
按照本发明的另一方面,图像数据是分割的检查对象的图像数据。由此可以仅分割在图像数据中的如下区域,该区域的运动应当被校正,从而图像度量是对于待校正的运动的特别有说服力的量度。此外图像数据然后包括更少的像素或体素,从而待处理的数据量被减小并且按照本发明的方法可以更快并且利用更小的计算开销进行。
按照本发明的另一方面,检查对象包括心脏的至少一个部分和/或至少一个冠状动脉。在心脏或冠状动脉的成像中运动补偿是特别有意义的,因为涉及快速运动的检查对象并且为了解释心脏或冠状动脉的图像数据而期望特别高的空间和时间分辨率。
按照另一方面进行分别对于来自于图像数据的各个截面图的图像度量的值的确定。因为图像度量对于在其录取区域中实际上发生运动的那个截面图是特别有说服力的。由此对于各个截面图的图像度量的值的确定通常比对于整组图像数据的图像度量的单个值的确定对运动伪影更敏感。
按照本发明的另一方面,图像度量是图像数据的熵、正性()、可压缩性或总变化。因为这些参数分别是对于运动伪影的良好量度。
按照本发明的另一方面为了确定运动场的优化方法是梯度下降方法,其中梯度是基于图像地被确定。
按照本发明的另一方面,将用于子图像数据的变换的运动场在控制点之间插值。由此可以对于大数量时间点确定运动场,从而可以以特别精确的方式实现图像数据或子图像数据的变换。
本发明还涉及一种具有程序代码的计算机程序,用于当计算机程序在计算机上执行时,执行所有按照本发明的方法步骤。借助按照本发明的计算机程序可以快速、可靠和可重复地执行按照本发明的方法。
此外本发明涉及一种以机器可读的载体形式的计算机程序产品,在其上可调用地存储了按照本发明的计算机程序。计算机程序在机器可读的载体上的存储允许以技术简单的方式和可重复地在不同的计算机上执行按照本发明的方法。
本发明同样涉及一种用于图像处理的计算机,包括用于存储计算机程序的存储器以及用于执行存储的计算机程序的处理器,其中在存储器上存储了计算机程序,当计算机程序在计算机上执行时,其执行按照本发明的方法。
本发明还涉及一种计算机断层造影***,包括辐射源以及与辐射源共同工作的辐射探测器,构造为用于在辐射源和检查对象之间相对旋转运动的情况下在角度范围内采集测量数据,还包括重建单元,构造为用于从测量数据重建图像数据,还包括用于图像处理的按照本发明的计算单元。
附图说明
以下借助在参考附图中示出的实施例详细描述和解释本发明。其中:
图1示出按照本发明的计算机断层造影***,包括计算机断层造影装置,
图2示出按照本发明的计算机断层造影***,包括C形臂X射线设备,
图3示出按照本发明的方法的框图,和
图4示出借助傅里叶变换确定子图像数据的图例。
具体实施方式
图1和图2分别示出按照本发明的计算机断层造影***1。在此处示出的例子中患者P在录取测量数据MEAS时卧于患者卧榻8上。患者卧榻8构造为将患者P在断层造影录取期间沿着***轴9移动。在录取测量数据MEAS时辐射源2和与辐射源2共同工作的辐射探测器5围绕***轴9旋转。辐射源9在此这样发出辐射,使得该辐射对于辐射探测器5是基本上可探测的。辐射通过各自的检查对象衰减,特别是通过辐射的吸收和反射。测量数据MEAS在此处示出的例子中是患者P的身体部位的多个投影,其中投影分别说明辐射通过患者P的身体部位的衰减。
在图1中示出的具有计算机断层造影装置的实施例中X射线探测器是具有多个行和列的探测器,而在图2中示出的C形臂X射线设备具有以平板探测器形式的X射线探测器。X射线探测器可以分别既作为闪烁晶体计数器又作为直接转换的X射线探测器构造。此外其还可以作为计数的X射线探测器构造,其构造为探测并且计数单个光子。此外计算机断层造影装置在图1中示出的例子中还具有两对互相共同工作的以X射线管形式的辐射源2、4和以X射线探测器形式的辐射探测器5、4。由此在此示出的计算机断层造影装置特别适合于多重能量录取,在所述多重能量录取中两个X射线管分别发出具有不同的能量谱的X射线辐射。在此处未示出的另一个实施方式中按照本发明的计算机断层造影***1的计算机断层造影装置仅具有分别一个X射线管作为辐射源和一个X射线探测器作为辐射探测器。在计算机断层造影装置中X射线管和X射线探测器集成在机架6中。在C形臂X射线设备中X射线管和X射线探测器通过C形臂7相连,所述C形臂又固定在机架6上。计算机断层造影装置的机架6可以这样构造,使得其可以围绕垂直于***轴9的至少一个轴翻转。在图2中示出的C形臂X射线设备的C形臂7分别可以沿着两个箭头中的一个摆动或旋转。
在此示出的计算机断层造影***1附加地还具有用于将造影剂注射到患者P的血液循环中的造影剂注射器11。由此可以借助造影剂这样录取测量数据MEAS,使得例如患者P的血管,特别是跳动的心脏的心室,在图像数据F中可以以提高的对比度被示出。此外利用造影剂注射器11还存在执行灌注测量的可能性,建议的方法同样适合于所述灌注测量。造影剂一般地理解为这样的试剂,其在成像方法中改善身体的结构和功能的显示。在此处呈现的申请的范围内造影剂既理解为传统的造影剂,诸如碘,也理解为示踪物,诸如18F,11C或13N。
此外按照本发明的计算机断层造影***还包括计算机10,其也称为工作站。在此示出的计算机10构造为用于控制计算机断层造影***的各个单元,例如用于控制患者卧榻8、造影剂注射器11和X射线管以及X射线探测器。计算机10与输出单元以及输入单元相连。输出单元例如是一个(或多个)LCD、等离子体或OLED显示器。输出单元上的输出包括例如图形用户界面或图像数据的输出。输入单元构造为用于输入数据诸如患者数据以及用于输入和选择对于计算机断层造影***的各个单元的参数。输入单元例如是键盘、鼠标、所谓的触摸屏或用于语音输入的麦克风。
在此处示出的例子中计算机10还构造为借助计算机断层造影装置或C形臂X射线设备的数据连接获得测量数据MEAS和借助重建单元从测量数据MEAS中重建图像数据F。在本发明的替换实施方式中计算机10与以重建计算机形式的计算***相连,通过另一个数据连接可以向其传输测量数据MEAS,由此计算***可以从测量数据MEAS中借助重建单元重建图像数据F。特别地在此计算机10可以作为客户机并且计算***可以作为服务器起作用。在该替换实施方式的变形中计算机10不具有重建单元并且不是构造为用于将测量数据MEAS重建为图像数据F。不取决于,重建单元是作为本地计算机10的部分还是作为必要时作为服务器起作用的、单独的计算***的部分,重建单元既可以作为软件也可以作为硬件实现。
此外计算机10还包括用于存储计算机程序Prg1-Prgn的存储器以及用于执行存储的计算机程序Prg1-Prgn的处理器。在此处示出的实施方式中在存储器上存储至少一个计算机程序,当计算机程序在计算机10上运行时,其执行按照本发明的方法的全部方法步骤。用于执行按照本发明的方法的方法步骤的计算机程序包括程序代码。此外计算机程序还可以作为可执行的文件构造和/或存储在计算机10之外的其他计算***上。例如计算机断层造影***1可以这样构造,使得计算机10将用于执行按照本发明的方法的计算机程序经过内联网或经过因特网加载到其内部工作存储器中。此外计算机程序可以可调用地存储在以机器可读的载体形式的计算机程序产品上。特别地,机器可读的载体可以是CD、DVD、蓝光盘、记忆棒或硬盘。
此外计算机10在图1中示出的实施方式中构造为借助EKG数据连接12获得并且处理患者P的EKG信号。借助获得的EKG信号,计算机10特别地构造为在所谓的前瞻性触发的范围内确定录取测量数据MEAS的时刻并且在定义的时刻开始相应的录取。此外计算机10构造为在所谓的回顾性触发的范围内将所录取的投影的仅一部分用于重建图像数据F以及提供关于心脏阶段的信息,其简化对已经重建的图像数据的进一步处理。本发明的典型应用范围是心脏成像和/或冠状血管成像。由此在以下描述的按照本发明的方法的检查对象特别地可以包括患者P的跳动的心脏和/或至少一个冠状动脉。
图3示出了按照本发明的方法的框图。在按照本发明的方法之前,在辐射源2和检查对象之间相对旋转运动的情况下在角度范围Ω内进行测量数据MEAS的录取。角度范围Ω在图1中示出的本发明的实施方式中垂直于***轴9延伸。一般公知,在的角度范围Ω下必须录取用于图像数据的空间三维重建的一组测量数据MEAS,以满足完整性标准。
在此说明了扇形角度,也就是扇形在径向方向上的开口角度。在C形臂X射线设备的情况下,人们不是提到扇形,而是通常提到X射线辐射的圆锥体,但是该圆锥体也具有定义的开口角度原则上,本发明可以基于完整的也可以基于不完整的测量数据MEAS的重建REKON。特别地,角度范围可以大于、小于或等于典型地,角度范围Ω是连续的角度范围。然而也可以考虑,按照本发明的方法利用不连续的角度范围Ω进行,例如,因为角度范围Ω的单个部分分别相应于作为检查对象的心脏的静止阶段。
角度范围Ω可以被划分为多个子角度范围ω_j。如果测量数据MEAS在角度范围Ω=360°中被录取,则其可以被划分为20个子角度范围,其中j=1,2…20,其中成立ω_j=18°。每个子角度范围ω_j现在可以对应于一个时间控制点t_j。时间控制点t_j例如是时间重点(Schwerpunkt)或时间间隔的开始或结束,在所述时间间隔中已经录取属于各自的子角度范围ω_j的投影。如果测量数据MEAS的录取例如持续一秒,则在不同的时间控制点之间的间隔在20个等距的子角度范围ω_j和保持相同的录取速度的情况下为50ms。当然也可以选择其它数量的子角度范围ω_j或时间控制点t_j。此外在本发明的其他实施方式中子角度范围ω_j或时间控制点t_j可以不是等距的。
在此处描述的按照本发明的方法的初始化INIT中,从测量数据MEAS中进行图像数据F的空间三维重建RECON。此外在初始化INIT的范围中可以确定图像数据F中的所谓的感兴趣区域ROI。通过选择感兴趣区域ROI可以减小图像数据F的范围,从而可以更快和以更小的计算开销执行按照本发明的方法。此外可以更灵敏地确定在选择的图像区域中的运动伪影。感兴趣区域ROI例如可以通过图像数据F的分割来确定。在这样的情况中,图像数据F在以下是指分割的图像数据。特别地可以将检查区域,也就是例如患者P的心脏和/或冠状动脉,在图像数据F中作为感兴趣区域ROI分割。
优化OPT是迭代的并且对于每个迭代包括多个步骤,所述步骤的目的在于基于图像地确定图像数据F的运动场s,并且具体来说借助图像度量L来确定,其中图像度量L是对于图像数据F中的运动伪影的量度。运动场s说明在图像空间内部单个图像元素(例如体素或像素)的运动,并且也可以作为矢量场描述。因为运动场s应当在控制点t_j被确定,所以运动场s特别地可以是具有三个空间的和一个时间的分量的四维矢量场。以下省略时间分量的依赖关系,只要这对于理解是不需要的。图像度量例如是图像数据F的熵、正性、可压缩性或总变化。图像的熵例如通过
给出,其中P(h,s)是对于一个体素在图像数据F内部具有例如以豪恩斯菲尔德单位h形式的强度的概率。概率分布P(h,s)可以借助不同的、专业人员常用的方法来计算。在优化OPT的第一步骤(i)中确定对于图像数据F的图像度量L的值。图像度量的值可以关于图像数据F的整个重建的体积来确定。替换地,图像度量L的多个值可以分别来对于自于图像数据F的子体积,特别是以单个截面图的形式的子体积被确定。如果确定图像度量L的多个值,则它们可以对于优化OPT的目的被综合,例如通过值的相加。
在优化OPT的另一个步骤(ii)中,确定子图像数据f_j,其中子图像数据分别对应于来自于子角度范围ω_j的测量数据的空间三维重建。在辐射源2和检查对象之间相对旋转运动期间在子角度范围ω_j内部已经录取的测量数据,在以下也称为子测量数据MEAS_j。在此分别对于一个时间控制点t_j存在一个对应的子角度范围ω_j和对应的子测量数据MEAS_j。然而在步骤(ii)中按照本发明不进行原始数据的重建,而是仅确定子图像数据f_j,其相应于由子测量数据MEAS_j的重建。由此可以基于图像地并且由此特别快速和以相对小的计算开销实现按照本发明的方法。子图像数据f_j的这样的确定例如可以借助傅里叶变换进行,如在对图4的描述中详细解释的那样。子图像数据f_j的确定也可以借助图像数据F的正投影和对应于子角度范围ω_j的一部分投影后的图像数据F的然后的反投影来进行。
在优化OPT的下一步骤(iii)中借助优化方法根据图像度量L确定图像数据F在控制点t_j上的运动场s。特别地可以这样进行优化OPT,使得图像度量的值被优化。换言之,优化OPT包括,图像度量的值力求趋于(entgegen strebt)极值,例如趋于最小值或最大值,但是不一定达到极值。原则上为了实现本发明可以使用用于该目的的不同的合适优化方法。在此详细地以所谓的梯度下降方法描述本发明。在每个迭代k中运动场s通过以下公式(1)描述,其中τ是待确定的参数:
作为对于第一迭代的开始值,原则上可以选择任意的场,特别是其项全部为0的场。梯度被基于图像地确定。因为运动场应当在控制点t_j处被确定,所以需要梯度为了简化,在以下省略迭代角标k。新的、迭代地产生的图像数据F可以作为子图像数据f_j的叠加、特别是作为和被显示。利用三维位置矢量x成立:
梯度可以特别地作为差分商计算,从而成立 其中e是三维的单位矢量。
在优化OPT的下一步骤(iv)中,将子图像数据f_j按照通过运动场s描述的运动进行转换。通过在第k+1个迭代步骤中将第k个迭代步骤的各个子图像数据f_j的图像元素(例如像素或体素)相应于运动场s移动,进行这样的转换。即转换也可以包括三维图像空间中子图像数据f_j的强度值的移动。此外子图像数据f_j的转换可以包括子图像数据f_j的插值。在本发明的实施方式中也插值运动场s。因为按照至此描述的方法,运动场仅在离散的时间上的控制点t_j处定义。
在优化OPT的下一步骤(v)中,通过转换的子图像数据f_j的拼合,特别是通过累加产生新的图像数据F。在本发明的实施方式中在步骤(v)之后迭代k结束并且下一个迭代k+1以步骤(i)开始,其中新拼合的图像数据代替初始的图像数据。此外,优化OPT包括中断标准的查询STOP。原则上,在此处描述的每个步骤(i)至(v)之间或期间都可以进行查询STOP。例如查询STOP可以包括,是否已经进行了一定的最大数量k_max的迭代,或者在两个先后跟随的迭代之间度量的值是否发生小于极限值的改变。查询STOP的步骤也可以包括多个中断标准的查询。在图3所示的例子中在步骤(i)之后查询中断标准。在满足中断标准Y的情况下以步骤(vi)继续,在不满足中断标准N的情况下以步骤(ii)继续。
步骤(vi)包括图像数据F的正投影以及借助运动场s然后进行的运动补偿的反投影。通过运动补偿的反投影可以弃用子图像数据f_j的迭代,从而借助步骤(vi)产生的图像数据F的质量特别高。由此特别有利地在迭代优化OPT的结束执行该步骤。然而步骤(vi)也可以在优化OPT期间被多次执行,例如用于输出优化OPT的可视的中间结果。
图4示出了借助傅里叶变换对子图像数据的确定的图例。断层造影重建基于投影的累加。按照所谓的傅里叶切片定理,一个投影的傅里叶变换与在选择的投影方向上衰减分布的一维傅里叶谱相同。图4a)示出了二维图像数据F的二维傅里叶谱。此处示出的傅里叶谱是空间频率的谱,其中强度的亮度是在频率域中。图4b)示出了掩模,该掩模在此处讨论的例子中被用于借助滤波函数对图像数据F的频谱进行掩模。因为单个投影的录取的时刻是已知的,所以也可以确定按照傅里叶切片定理计算的一维傅里叶谱的时序。掩模或滤波函数满足从傅里叶谱中选择扇区的目的,其中一个扇区相应于一个子角度范围ω_j。由此掩模或滤波函数也被称作定向的。图4c)示出掩模后的或滤波后的傅里叶谱。即按照在此示出的实施方式从一维的傅里叶谱的扇区中确定子图像数据f_j。由此也可以将扇区的时间控制点t_j根据与其对应的投影数据一一对应地确定。由此掩模或滤波函数可以确定特别是以时间重点形式的时间控制点t_j。并且图4d)最后根据掩模后的频谱到空间域中的傅里叶变换示出了子图像数据f_j。
按照本发明的方法的各个步骤不一定必须按照在此描述的顺序执行。而是专业人员理解,可以改变步骤的顺序,只要在技术上是有意义的。
Claims (15)
1.一种用于对运动的检查对象的图像数据(F)进行运动补偿的方法,
其中,所述图像数据(F)是由测量数据(MEAS)的空间三维重建(RECON),其中,在辐射源(2)和检查对象之间相对旋转运动的情况下在角度范围(Ω)内录取所述测量数据(MEAS),其中,所述角度范围(Ω)包括多个子角度范围,其中每个子角度范围(ω_j)对应于一个时间控制点(t_j),
其中,所述方法包括以下步骤:
(i)确定对于图像数据(F)的图像度量(L)的至少一个值,其中,所述图像度量(L)是对于图像数据(F)中运动伪影的量度,
(ii)确定子图像数据(f_j),其中,所述子图像数据(f_j)分别对应于来自于子角度范围(ω_j)的测量数据的空间三维重建,
(iii)借助优化方法根据所述图像度量(L)在控制点(t_j)上确定图像数据(F)的运动场(s),其中,所述运动场(s)是对于图像数据(F)的运动的量度,
(iv)将所述子图像数据(t_j)按照通过所述运动场(s)描述的运动进行变换,
(v)通过将所述子图像数据(f_j)拼合产生新的图像数据(F),
其中,迭代地一直进行步骤(i)至(v),直到达到中断标准。
2.根据权利要求1所述的方法,
其中,所述角度范围是对于断层造影地完整的重建所需的最小角度范围。
3.根据权利要求1或2所述的方法,
其中,借助以下步骤确定所述子图像数据(f_j):
-图像数据(F)到频率域的傅里叶变换,
-借助滤波函数来掩模图像数据(F)的频谱,
-将掩模后的频谱转换到空间域的傅里叶变换。
4.根据权利要求1或2所述的方法,
其中,分别借助图像数据(F)的正投影以及对应于子角度范围(ω_j)的一部分投影后的图像数据(F)的然后的反投影,来确定所述子图像数据(f_j)。
5.根据权利要求1至4中任一项所述的方法,还包括:
(vi)借助所述运动场(s)进行图像数据(F)的正投影以及然后的运动补偿的反投影。
6.根据权利要求1至5中任一项所述的方法,
其中,所述图像数据(F)是被分割的检查对象的图像数据。
7.根据权利要求1至6中任一项所述的方法,
其中,所述检查对象包括心脏的至少一个部分和/或至少一个冠状动脉。
8.根据权利要求1至7中任一项所述的方法,
其中,进行分别对于来自于图像数据(F)的各个截面图的图像度量(L)的值的确定。
9.根据权利要求1至8中任一项所述的方法,
其中,所述图像度量(L)是以下图像参数之一:
-熵,
-正性,
-可压缩性,
-总变化。
10.根据权利要求1至9中任一项所述的方法,
其中,用于确定所述运动场(s)的优化方法是梯度下降方法,其中,基于图像来确定梯度
11.根据权利要求1至10中任一项所述的方法,
其中,对用于子图像数据(f_j)的变换的运动场(s)在控制点(t_j)之间进行插值。
12.一种具有程序代码的计算机程序,用于当所述计算机程序在计算机(10)上执行时,执行按照权利要求1至11中任一项所述的方法的所有方法步骤。
13.一种以机器可读的载体形式的计算机程序产品,在其上可调用地存储了按照权利要求12所述的计算机程序。
14.一种用于图像处理的计算机(10),包括用于存储计算机程序的存储器以及用于执行存储的计算机程序的处理器,其中,在所述存储器上存储了计算机程序,当计算机程序在所述计算机(10)上执行时,其执行按照权利要求1至11中任一项所述的方法的所有方法步骤。
15.一种计算机断层造影***(1),包括辐射源(2)以及与辐射源(2)共同工作的辐射探测器(3),所述辐射探测器被构造为用于在辐射源(2)和检查对象之间相对旋转运动的情况下在角度范围内采集测量数据,还包括重建单元,所述重建单元被构造为用于从测量数据(MEAS)中重建图像数据(F),以及还包括按照权利要求14的用于图像处理的计算机(10)。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102013217351.2 | 2013-08-30 | ||
DE102013217351.2A DE102013217351B4 (de) | 2013-08-30 | 2013-08-30 | Bildbasierte Bewegungskompensation von Bilddaten |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104424659A true CN104424659A (zh) | 2015-03-18 |
CN104424659B CN104424659B (zh) | 2020-06-19 |
Family
ID=52470355
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410429770.8A Active CN104424659B (zh) | 2013-08-30 | 2014-08-28 | 图像数据的基于图像的运动补偿 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9619869B2 (zh) |
CN (1) | CN104424659B (zh) |
DE (1) | DE102013217351B4 (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110473269A (zh) * | 2019-08-08 | 2019-11-19 | 上海联影医疗科技有限公司 | 一种图像重建方法、***、设备及存储介质 |
CN110870776A (zh) * | 2018-08-31 | 2020-03-10 | 西门子医疗有限公司 | 用于三维数字减影血管造影图像的方法和装置 |
CN111093506A (zh) * | 2017-07-27 | 2020-05-01 | 皇家飞利浦有限公司 | 运动补偿心脏瓣膜重建 |
CN113616223A (zh) * | 2020-05-07 | 2021-11-09 | 西门子医疗有限公司 | 提供掩模图像 |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015112816A1 (en) * | 2014-01-24 | 2015-07-30 | E-Image Data Corporation | Multi-mode image capture systems and methods |
CN104318595A (zh) * | 2014-10-21 | 2015-01-28 | 沈阳东软医疗***有限公司 | Ct图像的运动向量场的计算方法及装置 |
KR101783964B1 (ko) * | 2015-04-09 | 2017-10-10 | 삼성전자주식회사 | 단층 촬영 장치 및 그에 따른 단층 영상 복원 방법 |
KR101728046B1 (ko) | 2015-08-27 | 2017-04-18 | 삼성전자주식회사 | 단층 영상 복원 장치 및 그에 따른 단층 영상 복원 방법 |
JP6437163B1 (ja) * | 2015-10-28 | 2018-12-12 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | コンピュータ断層撮影画像生成装置 |
FI20155856A (fi) * | 2015-11-18 | 2017-05-19 | Teknologian Tutkimuskeskus Vtt Oy | Menetelmiä ja laitteita läpitunkevaa kuvantamista varten |
US9931095B2 (en) | 2016-03-30 | 2018-04-03 | General Electric Company | Method for segmenting small features in an image volume |
US10565744B2 (en) | 2016-06-30 | 2020-02-18 | Samsung Electronics Co., Ltd. | Method and apparatus for processing a medical image to reduce motion artifacts |
KR20180041007A (ko) * | 2016-10-13 | 2018-04-23 | 삼성전자주식회사 | 의료 영상 처리 장치 및 방법 |
US10915603B2 (en) * | 2016-12-09 | 2021-02-09 | Korea Advanced Institute Of Science And Technology | Method for estimating suitability as multi-screen projecting type theatre system |
JP6882135B2 (ja) * | 2017-10-06 | 2021-06-02 | キヤノン株式会社 | 画像処理装置、画像処理方法及びプログラム |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1809841A (zh) * | 2003-06-18 | 2006-07-26 | 皇家飞利浦电子股份有限公司 | 运动补偿的重建技术 |
US20120275656A1 (en) * | 2011-04-28 | 2012-11-01 | Jan Boese | Method for supplying a 3d x-ray image data record for a moving object with highly absorbent material |
CN102982582A (zh) * | 2011-05-05 | 2013-03-20 | 西门子公司 | 对周期性运动的对象获得三维x射线图像数据组的方法 |
CN103198497A (zh) * | 2011-09-28 | 2013-07-10 | 西门子公司 | 确定运动场和利用运动场进行运动补偿重建的方法和*** |
CN103190927A (zh) * | 2011-09-28 | 2013-07-10 | 西门子公司 | 确定运动场和产生运动补偿的ct 图像数据组的方法和*** |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8224056B2 (en) | 2009-12-15 | 2012-07-17 | General Electronic Company | Method for computed tomography motion estimation and compensation |
DE102010024684B8 (de) | 2010-06-23 | 2019-08-22 | Siemens Healthcare Gmbh | Verfahren zur Rekonstruktion von Bilddaten eines bewegten Untersuchungsobjektes, Steuer- und Recheneinheit, Computertomographiesystem und Computerprogramm |
DE102011083641A1 (de) | 2011-09-28 | 2013-03-28 | Siemens Ag | Verfahren, Rechensystem und CT-System zur Bestimmung eines Bewegungsfeldes und zur Erzeugung eines bewegungskompensierten Cardio-CT-Bilddatensatzes |
-
2013
- 2013-08-30 DE DE102013217351.2A patent/DE102013217351B4/de active Active
-
2014
- 2014-08-12 US US14/457,215 patent/US9619869B2/en active Active
- 2014-08-28 CN CN201410429770.8A patent/CN104424659B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1809841A (zh) * | 2003-06-18 | 2006-07-26 | 皇家飞利浦电子股份有限公司 | 运动补偿的重建技术 |
US20120275656A1 (en) * | 2011-04-28 | 2012-11-01 | Jan Boese | Method for supplying a 3d x-ray image data record for a moving object with highly absorbent material |
CN102982582A (zh) * | 2011-05-05 | 2013-03-20 | 西门子公司 | 对周期性运动的对象获得三维x射线图像数据组的方法 |
CN103198497A (zh) * | 2011-09-28 | 2013-07-10 | 西门子公司 | 确定运动场和利用运动场进行运动补偿重建的方法和*** |
CN103190927A (zh) * | 2011-09-28 | 2013-07-10 | 西门子公司 | 确定运动场和产生运动补偿的ct 图像数据组的方法和*** |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111093506A (zh) * | 2017-07-27 | 2020-05-01 | 皇家飞利浦有限公司 | 运动补偿心脏瓣膜重建 |
CN111093506B (zh) * | 2017-07-27 | 2023-08-01 | 皇家飞利浦有限公司 | 运动补偿心脏瓣膜重建 |
CN110870776A (zh) * | 2018-08-31 | 2020-03-10 | 西门子医疗有限公司 | 用于三维数字减影血管造影图像的方法和装置 |
CN110870776B (zh) * | 2018-08-31 | 2023-08-25 | 西门子医疗有限公司 | 用于三维数字减影血管造影图像的方法和装置 |
CN110473269A (zh) * | 2019-08-08 | 2019-11-19 | 上海联影医疗科技有限公司 | 一种图像重建方法、***、设备及存储介质 |
CN110473269B (zh) * | 2019-08-08 | 2023-05-26 | 上海联影医疗科技股份有限公司 | 一种图像重建方法、***、设备及存储介质 |
CN113616223A (zh) * | 2020-05-07 | 2021-11-09 | 西门子医疗有限公司 | 提供掩模图像 |
Also Published As
Publication number | Publication date |
---|---|
US9619869B2 (en) | 2017-04-11 |
DE102013217351A1 (de) | 2015-03-05 |
CN104424659B (zh) | 2020-06-19 |
DE102013217351B4 (de) | 2023-05-11 |
US20150063534A1 (en) | 2015-03-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104424659A (zh) | 图像数据的基于图像的运动补偿 | |
US10748293B2 (en) | Tomography apparatus and method for reconstructing tomography image thereof | |
Tang et al. | Temporal resolution improvement in cardiac CT using PICCS (TRI‐PICCS): Performance studies | |
CN109389655B (zh) | 时变数据的重建 | |
US9576391B2 (en) | Tomography apparatus and method of reconstructing a tomography image by the tomography apparatus | |
CN107427274B (zh) | 断层扫描设备及其用于重构断层扫描图像的方法 | |
US8055050B2 (en) | Motion compensation in energy-sensitive computed tomography | |
JP2021502836A (ja) | 機械学習を使用した画像生成 | |
CN103198497A (zh) | 确定运动场和利用运动场进行运动补偿重建的方法和*** | |
JP2017189612A (ja) | 放射線画像診断装置及び医用画像処理装置 | |
EP3107457B1 (en) | Tomography apparatus and method of reconstructing a tomography image by the tomography apparatus | |
US8768030B2 (en) | CT measurement with multiple X-ray sources | |
Rohkohl et al. | Interventional 4D motion estimation and reconstruction of cardiac vasculature without motion periodicity assumption | |
CN103027705B (zh) | 产生运动补偿的ct 图像数据组的方法和*** | |
CN102013089A (zh) | 用于噪声减少的迭代ct图像滤波器 | |
CN102376084B (zh) | 使用各向异性噪声模型对ct图像的迭代图像滤波 | |
CN101991428A (zh) | 用于在心脏ct中改进时间分辨率的ct图像重建 | |
JP2016152916A (ja) | X線コンピュータ断層撮像装置及び医用画像処理装置 | |
KR20170017156A (ko) | 단층 영상 장치 및 그에 따른 단층 영상 복원 방법 | |
Van Stevendaal et al. | A motion‐compensated scheme for helical cone‐beam reconstruction in cardiac CT angiography | |
Pack et al. | Four‐dimensional computed tomography of the left ventricle, Part I: motion artifact reduction | |
JP6479919B2 (ja) | 流動データの再構築 | |
JPWO2019044983A1 (ja) | X線ct装置および画像生成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220207 Address after: Erlangen Patentee after: Siemens Healthineers AG Address before: Munich, Germany Patentee before: SIEMENS AG |
|
TR01 | Transfer of patent right |