CN104424659B - 图像数据的基于图像的运动补偿 - Google Patents

图像数据的基于图像的运动补偿 Download PDF

Info

Publication number
CN104424659B
CN104424659B CN201410429770.8A CN201410429770A CN104424659B CN 104424659 B CN104424659 B CN 104424659B CN 201410429770 A CN201410429770 A CN 201410429770A CN 104424659 B CN104424659 B CN 104424659B
Authority
CN
China
Prior art keywords
image data
image
sub
motion
data
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.)
Active
Application number
CN201410429770.8A
Other languages
English (en)
Other versions
CN104424659A (zh
Inventor
T.奥尔门丁格
H.布鲁德
T.弗洛尔
C.罗科尔
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 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 AG filed Critical Siemens AG
Publication of CN104424659A publication Critical patent/CN104424659A/zh
Application granted granted Critical
Publication of CN104424659B publication Critical patent/CN104424659B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/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
    • 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)
  • Quality & Reliability (AREA)
  • Architecture (AREA)
  • Computer Graphics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明涉及图像数据运动补偿的方法。图像数据是由测量数据的空间三维重建，在辐射源和检查对象之间相对旋转运动下在角度范围内录取测量数据。角度范围包括多个子角度范围，每个子角度范围对应一个时间控制点。在迭代方法第一步骤期间确定图像数据的图像度量的至少一个值，图像度量是图像数据中运动伪影的量度。基于图像数据根据图像度量能确定图像数据运动补偿的运动场。在第二步骤中确定子图像数据，子图像数据分别对应来自子角度范围的测量数据的空间三维重建。在第三步骤中借助优化方法根据图像度量在控制点上确定图像数据的运动场，在第四步骤中将子图像数据按照运动场描述的运动进行变换。在第五步骤中通过拼合子图像数据产生新的图像数据。

Description

图像数据的基于图像的运动补偿
技术领域
本发明涉及一种用于图像数据的运动补偿的方法以及计算机程序、计算机程序产品和计算机断层造影系统。
背景技术
计算机断层造影是主要用于医学诊断以及用于材料检查的成像方法。在计算机断层造影中为了录取空间上三维的图像数据，辐射源以及与该辐射源一起工作的辐射探测器围绕检查对象旋转。在旋转运动期间在一个角度范围内录取测量数据。测量数据是多个投影，所述投影包含关于辐射通过检查对象的衰减的信息。测量数据也称为原始数据。从这些测量数据中然后可以重建图像数据，例如借助所谓的滤波反投影或借助迭代的重建方法。如果检查对象在录取期间运动，则在图像数据的重建中会由于运动而形成模糊和伪影。
因此一般公知的是，在断层造影录取期间的心脏运动会导致不一致和通过运动伪影损害图像数据。由此会严重限制图像数据的临床可用性。为了避免运动伪影，在现代的断层造影心脏成像中使用借助回顾性或前瞻性触发的、测量数据的相位相关的录取。因此在前瞻性触发中仅录取在大约心脏的静止阶段的一定的时间间隔中的数据并且用于图像重建。该方案的共同目标是“冻结”心脏运动，以减小数据不一致性和由此提高图像质量。
在部分和周期运动的对象的情况中用于提高图像质量的一种方案从文献DE 102011 083 643 A1公知。在那里公开了一种用于确定CT图像数据组的图像空间中的运动场(Bewegungsfeld)的方法。在此对于预先给出的运动阶段和投影角度范围采集CT系统的投影数据组。此外通过将利用第一分析重建算法在使用运动场的条件下运动补偿地重建的、第一图像分辨率的CT图像数据组的、基于原始数据的或变换到图像空间中的成本函数最小化，进行运动场的迭代确定。
发明内容
本发明要解决的技术问题是，如何能够基于图像数据对运动的检查对象执行图像数据的运动补偿。
以下描述关于要求保护的装置和关于要求保护的方法的技术问题的按照本发明的解决方案。在此提到的特征、优点或替换实施方式同样也适用于其他要求保护的内容，反之亦然。换言之，例如针对系统的具体权利要求也可以利用关于方法描述或要求保护的特征来扩展。方法的相应功能特征在此通过相应的具体模块构造。
本发明涉及一种用于对运动的检查对象的图像数据进行运动补偿的方法。图像数据是由测量数据的空间三维重建，其中在辐射源和检查对象之间相对旋转运动的情况下在角度范围内录取测量数据。该角度范围包括多个子角度范围，其中每个子角度范围对应于一个时间控制点。特别地，时间控制点可以是通过时间间隔定义的时刻，在该时间间隔期间在各自的子角度范围内录取测量数据。本发明的核心是以下的迭代方法。在按照本发明的方法的第一步骤期间确定对于图像数据的图像度量L的至少一个值，其中图像度量是对于图像数据中运动伪影的量度。发明人现在已经认识到，基于图像数据地根据图像度量可以确定用于图像数据的运动补偿的运动场s，其中运动场是按照对于图像数据的运动的量度。为此在第二步骤中确定子图像数据，其中子图像数据分别对应于来自于子角度范围的测量数据的空间三维重建。在第三步骤中借助优化方法根据图像度量在控制点上确定图像数据的运动场s，从而在第四步骤中可以将子图像数据按照通过运动场描述的运动进行变换。在第五步骤中通过将子图像数据拼合产生新的图像数据。迭代方法的这些步骤一直进行，直到达到中断标准。新的图像数据通过被变换的子图像数据的变换被运动补偿并且相对于初始的图像数据具有改善的图像质量。通过按照本发明基于图像地进行该方法，与其他基于原始数据的方法相比，需要更小的计算开销并且可以特别地在所谓的工作站上和独立于重建计算机进行。由此证明按照本发明的方法是对于临床实践特别实用的。
按照本发明的另一方面，角度范围是对于断层造影完整的重建所需的最小角度范围。在这样的情况中，角度范围通过
Figure BDA0000561344140000021
给出，其中
Figure BDA0000561344140000022
说明扇形角度，也就是在径向方向上通过辐射源发射的辐射扇形的开口角度。
按照本发明的另一方面，通过首先进行图像数据到频率域的傅里叶变换，然后借助滤波函数掩模图像数据并且借助傅里叶变换又将掩模的频谱转换到空间域，确定子图像数据。本发明的该方面的优点是，子图像数据的确定可以利用小的计算开销实现并且无需访问原始数据。
按照本发明的另一方面，分别借助图像数据的正投影和对应于子角度范围的一部分投影后的图像数据的然后的反投影来确定子图像数据。本发明的该方面也可以利用小的计算开销进行并且无需访问原始数据。
按照本发明的另一方面，在迭代方法执行期间和/或之后，借助运动场进行图像数据的正投影以及然后的运动补偿的反投影。在该方面中尽可能无失真地进行图像数据的运动补偿，从而实现特别高的图像质量。
按照本发明的另一方面，图像数据是分割的检查对象的图像数据。由此可以仅分割在图像数据中的如下区域，该区域的运动应当被校正，从而图像度量是对于待校正的运动的特别有说服力的量度。此外图像数据然后包括更少的像素或体素，从而待处理的数据量被减小并且按照本发明的方法可以更快并且利用更小的计算开销进行。
按照本发明的另一方面，检查对象包括心脏的至少一个部分和/或至少一个冠状动脉。在心脏或冠状动脉的成像中运动补偿是特别有意义的，因为涉及快速运动的检查对象并且为了解释心脏或冠状动脉的图像数据而期望特别高的空间和时间分辨率。
按照另一方面进行分别对于来自于图像数据的各个截面图的图像度量的值的确定。因为图像度量对于在其录取区域中实际上发生运动的那个截面图是特别有说服力的。由此对于各个截面图的图像度量的值的确定通常比对于整组图像数据的图像度量的单个值的确定对运动伪影更敏感。
按照本发明的另一方面，图像度量是图像数据的熵、正性(
Figure BDA0000561344140000031
)、可压缩性或总变化。因为这些参数分别是对于运动伪影的良好量度。
按照本发明的另一方面为了确定运动场的优化方法是梯度下降方法，其中梯度
Figure BDA0000561344140000032
是基于图像地被确定。
按照本发明的另一方面，将用于子图像数据的变换的运动场在控制点之间插值。由此可以对于大数量时间点确定运动场，从而可以以特别精确的方式实现图像数据或子图像数据的变换。
本发明还涉及一种具有程序代码的计算机程序，用于当计算机程序在计算机上执行时，执行所有按照本发明的方法步骤。借助按照本发明的计算机程序可以快速、可靠和可重复地执行按照本发明的方法。
此外本发明涉及一种以机器可读的载体形式的计算机程序产品，在其上可调用地存储了按照本发明的计算机程序。计算机程序在机器可读的载体上的存储允许以技术简单的方式和可重复地在不同的计算机上执行按照本发明的方法。
本发明同样涉及一种用于图像处理的计算机，包括用于存储计算机程序的存储器以及用于执行存储的计算机程序的处理器，其中在存储器上存储了计算机程序，当计算机程序在计算机上执行时，其执行按照本发明的方法。
本发明还涉及一种计算机断层造影系统，包括辐射源以及与辐射源共同工作的辐射探测器，构造为用于在辐射源和检查对象之间相对旋转运动的情况下在角度范围内采集测量数据，还包括重建单元，构造为用于从测量数据重建图像数据，还包括用于图像处理的按照本发明的计算单元。
附图说明
以下借助在参考附图中示出的实施例详细描述和解释本发明。其中：
图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延伸。一般公知，在
Figure BDA0000561344140000061
的角度范围Ω下必须录取用于图像数据的空间三维重建的一组测量数据MEAS，以满足完整性标准。
在此
Figure BDA0000561344140000062
说明了扇形角度，也就是扇形在径向方向上的开口角度。在C形臂X射线设备的情况下，人们不是提到扇形，而是通常提到X射线辐射的圆锥体，但是该圆锥体也具有定义的开口角度
Figure BDA0000561344140000071
原则上，本发明可以基于完整的也可以基于不完整的测量数据MEAS的重建REKON。特别地，角度范围可以大于、小于或等于
Figure BDA0000561344140000072
典型地，角度范围Ω是连续的角度范围。然而也可以考虑，按照本发明的方法利用不连续的角度范围Ω进行，例如，因为角度范围Ω的单个部分分别相应于作为检查对象的心脏的静止阶段。
角度范围Ω可以被划分为多个子角度范围ω_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的熵、正性、可压缩性或总变化。图像的熵例如通过
Figure BDA0000561344140000081
给出，其中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)描述，其中τ是待确定的参数：
Figure BDA0000561344140000082
作为对于第一迭代的开始值，原则上可以选择任意的场，特别是其项全部为0的场。梯度
Figure BDA0000561344140000083
被基于图像地确定。因为运动场应当在控制点t_j处被确定，所以需要梯度
Figure BDA0000561344140000091
为了简化，在以下省略迭代角标k。新的、迭代地产生的图像数据F可以作为子图像数据f_j的叠加、特别是作为和被显示。利用三维位置矢量x成立：
Figure BDA0000561344140000092
梯度可以特别地作为差分商计算，从而成立
Figure BDA0000561344140000094
其中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 (14)

1.一种用于对运动的检查对象的图像数据进行运动补偿的方法，
其中，所述图像数据获自测量数据的空间三维重建，其中，在辐射源和检查对象之间相对旋转运动的情况下在角度范围内录取所述测量数据，其中，所述角度范围包括多个子角度范围，其中每个子角度范围对应于一个时间控制点，
其中，所述方法包括以下步骤：
(i)确定对于图像数据的图像度量的至少一个值，其中，所述图像度量是对于图像数据中运动伪影的量度，
(ii)确定子图像数据，其中，所述子图像数据分别对应于来自于子角度范围的测量数据的空间三维重建，
(iii)借助优化方法根据所述图像度量在时间控制点上确定图像数据的运动场，其中，所述运动场是对于图像数据的运动的量度，
(iv)将所述子图像数据按照通过所述运动场描述的运动进行转换，
(v)通过将转换的子图像数据拼合产生新的图像数据，
其中，迭代地一直进行步骤(i)至(v)，直到达到中断标准,
其中，进行分别对于来自于图像数据的各个截面图的图像度量的值的确定。
2.根据权利要求1所述的方法，
其中，所述角度范围是对于断层造影地完整的重建所需的最小角度范围。
3.根据权利要求1或2所述的方法，
其中，借助以下步骤确定所述子图像数据：
-图像数据到频率域的傅里叶变换，
-借助滤波函数来掩模图像数据的频谱，
-将掩模后的频谱转换到空间域的傅里叶变换。
4.根据权利要求1或2所述的方法，
其中，分别借助图像数据的正投影以及对应于子角度范围的一部分投影后的图像数据的然后的反投影，来确定所述子图像数据。
5.根据权利要求1或2所述的方法，还包括：
(vi)借助所述运动场进行图像数据的正投影以及然后的运动补偿的反投影。
6.根据权利要求1或2所述的方法，
其中，所述图像数据是被分割的检查对象的图像数据。
7.根据权利要求1或2所述的方法，
其中，所述检查对象包括心脏的至少一个部分。
8.根据权利要求1或2所述的方法，
其中，所述检查对象包括至少一个冠状动脉。
9.根据权利要求1或2所述的方法，
其中，所述图像度量是以下图像参数之一：
-熵，
-正性，
-可压缩性，
-总变化。
10.根据权利要求1或2所述的方法，
其中，用于确定所述运动场的优化方法是梯度下降方法，其中，基于图像来确定梯度
Figure FDA0002264219700000021
11.根据权利要求1或2所述的方法，
其中，对用于子图像数据的变换的运动场在控制点之间进行插值。
12.一种机器可读的载体，在其上可调用地存储了计算机程序，具有程序代码，用于当所述计算机程序在计算机上执行时，执行按照权利要求1至11中任一项所述的方法的所有方法步骤。
13.一种用于图像处理的计算机，包括用于存储计算机程序的存储器以及用于执行存储的计算机程序的处理器，其中，在所述存储器上存储了计算机程序，当计算机程序在所述计算机上执行时，其执行按照权利要求1至11中任一项所述的方法的所有方法步骤。
14.一种计算机断层造影系统，包括辐射源以及与辐射源共同工作的辐射探测器，所述辐射探测器被构造为用于在辐射源和检查对象之间相对旋转运动的情况下在角度范围内采集测量数据，还包括重建单元，所述重建单元被构造为用于从测量数据中重建图像数据，以及还包括按照权利要求13的用于图像处理的计算机。
CN201410429770.8A 2013-08-30 2014-08-28 图像数据的基于图像的运动补偿 Active CN104424659B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102013217351.2A DE102013217351B4 (de) 2013-08-30 2013-08-30 Bildbasierte Bewegungskompensation von Bilddaten
DE102013217351.2 2013-08-30

Publications (2)

Publication Number Publication Date
CN104424659A CN104424659A (zh) 2015-03-18
CN104424659B true 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)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
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 삼성전자주식회사 단층 영상 복원 장치 및 그에 따른 단층 영상 복원 방법
EP3331449B1 (en) * 2015-10-28 2018-12-19 Koninklijke Philips N.V. Computed tomography image generation apparatus
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
EP3434187A1 (en) * 2017-07-27 2019-01-30 Koninklijke Philips N.V. Motion compensated cardiac valve reconstruction
JP6882135B2 (ja) * 2017-10-06 2021-06-02 キヤノン株式会社 画像処理装置、画像処理方法及びプログラム
EP3616621B1 (de) * 2018-08-31 2023-03-15 Siemens Healthcare GmbH Verfahren zur 3d dsa und vorrichtung
CN110473269B (zh) * 2019-08-08 2023-05-26 上海联影医疗科技股份有限公司 一种图像重建方法、系统、设备及存储介质
DE102020205762A1 (de) * 2020-05-07 2021-11-11 Siemens Healthcare Gmbh Bereitstellen eines Maskenbildes

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004111946A1 (en) * 2003-06-18 2004-12-23 Philips Intellectual Property & Standards Gmbh Motion compensated reconstruction technique
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
DE102011017710B4 (de) 2011-04-28 2021-03-04 Siemens Healthcare Gmbh Verfahren zum Bereitstellen eines 3D-Röntgenbilddatensatzes zu einem sich bewegenden Objekt mit stark absorbierendem Material
DE102011075287B4 (de) 2011-05-05 2016-07-21 Siemens Healthcare Gmbh Verfahren zum Gewinnen eines 3D-Röntgenbilddatensatzes zu einem sich periodisch bewegenden Bildobjekt
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
DE102011083646A1 (de) 2011-09-28 2013-03-28 Siemens Aktiengesellschaft Verfahren, Rechensystem und CT-System zur Bestimmung eines Bewegungsfeldes und zur Erzeugung eines bewegungskompensierten CT-Bilddatensatzes eines sich teilweise bewegenden Objektes
DE102011083643A1 (de) 2011-09-28 2013-03-28 Siemens Aktiengesellschaft Verfahren, Rechensystem und CT-System zur Bestimmung eines Bewegungsfeldes und zur bewegungskompensierenden Rekonstruktion mit diesem Bewegungsfeld

Also Published As

Publication number Publication date
DE102013217351B4 (de) 2023-05-11
DE102013217351A1 (de) 2015-03-05
US20150063534A1 (en) 2015-03-05
US9619869B2 (en) 2017-04-11
CN104424659A (zh) 2015-03-18

Similar Documents

Publication Publication Date Title
CN104424659B (zh) 图像数据的基于图像的运动补偿
US10748293B2 (en) Tomography apparatus and method for reconstructing tomography image thereof
US9754390B2 (en) Reconstruction of time-varying data
US9576391B2 (en) Tomography apparatus and method of reconstructing a tomography image by the tomography apparatus
JP5860607B2 (ja) トモグラフィデータ収集及び画像再構成のためのシステム及び方法
EP2567359B1 (en) Image data registration for dynamic perfusion ct
US8055050B2 (en) Motion compensation in energy-sensitive computed tomography
US8600137B2 (en) Iterative CT image reconstruction with a four-dimensional noise filter
US8768030B2 (en) CT measurement with multiple X-ray sources
EP3107457B1 (en) Tomography apparatus and method of reconstructing a tomography image by the tomography apparatus
US20080267455A1 (en) Method for Movement Compensation of Image Data
RU2598159C2 (ru) Способ реконструкции изображений для отфильтрованной обратной проекции в томографии с ограниченным углом обзора объекта
CN110751702A (zh) 图像重建方法、系统、装置及存储介质
KR20150099381A (ko) 단층 촬영 장치 및 그에 따른 단층 영상 복원 방법
JP2016152916A (ja) Ｘ線コンピュータ断層撮像装置及び医用画像処理装置
EP3084726B1 (en) Moving structure motion compensation in imaging
KR20170017156A (ko) 단층 영상 장치 및 그에 따른 단층 영상 복원 방법
JP7122301B2 (ja) 画像データ内のノイズ低減
US20160331341A1 (en) Segmentation of moving structure in image data
JP6479919B2 (ja) 流動データの再構築
Manhart et al. Fast dynamic reconstruction algorithm with joint bilateral filtering for perfusion C-arm CT
Lindström Söraas Respiratory Motion Correction in PET Imaging: Comparative Analysis of External Device and Data-driven Gating Approaches

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