CN107886554B - 流数据的重构 - Google Patents

流数据的重构 Download PDF

Info

Publication number
CN107886554B
CN107886554B CN201710906720.8A CN201710906720A CN107886554B CN 107886554 B CN107886554 B CN 107886554B CN 201710906720 A CN201710906720 A CN 201710906720A CN 107886554 B CN107886554 B CN 107886554B
Authority
CN
China
Prior art keywords
dimensional
contrast agent
flowsheet
flow
projection
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
CN201710906720.8A
Other languages
English (en)
Other versions
CN107886554A (zh
Inventor
S.沙弗
S.格里施
M.科瓦施克
C.罗科尔
K.罗亚蒂
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 Healthcare GmbH
Original Assignee
Siemens Healthcare GmbH
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 Healthcare GmbH filed Critical Siemens Healthcare GmbH
Publication of CN107886554A publication Critical patent/CN107886554A/zh
Application granted granted Critical
Publication of CN107886554B publication Critical patent/CN107886554B/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
    • 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
    • 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/037—Emission tomography
    • 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/486—Diagnostic techniques involving generating temporal series of image 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/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
    • 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/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/5288—Devices using data or image processing specially adapted for radiation diagnosis involving retrospective matching to a physiological signal
    • 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/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • 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/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • 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
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/38—Registration of image sequences
    • 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/501—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 head, e.g. neuroimaging or craniography
    • 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/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
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • 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/10016—Video; Image sequence
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10116—X-ray image
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30101—Blood vessel; Artery; Vein; Vascular
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00—Image generation
    • G06T2211/40—Computed tomography
    • G06T2211/404—Angiography
    • 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/424—Iterative
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00—Image generation
    • G06T2211/40—Computed tomography
    • G06T2211/432—Truncation
    • 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)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Algebra (AREA)
  • Pulmonology (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Physiology (AREA)
  • Quality & Reliability (AREA)
  • Vascular Medicine (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本申请涉及流数据的重构，其中具体描述了用于促进流数据的重构的技术。依照一个方面，框架接收四维投影图像数据集并且配准图像数据集中的一对或多对时间上相邻的投影图像。可以基于经配准的对确定二维流图。框架然后可以根据心相对二维流图进行分类，并且基于经分类的二维流图重构三维流图。

Description

流数据的重构
技术领域
本公开一般地涉及流数据的自动化或部分自动化的重构。
背景技术
血管造影术是用于基于诊断成像方法而表示血管的常见方法。为了X射线检查之下的血管的改进的表示，已经开发了数字减影血管造影术（DSA）。DSA是使用在介入放射学中以清楚地可视化骨头或致密软组织环境中的脉管系统的X射线技术。图像通过在已经将对比剂引入到感兴趣的结构或组织中之后从随后的图像减去“强化前图像”或掩膜而产生。这些图像可以用于提供示出感兴趣的结构或组织随时间的移动的时间分辨或时变信息。
从血管造影图像导出流信息在相当长的时间内是感兴趣的对象。当前的技术利用二维（2D）DSA图像，所述图像在跨一系列心相的高帧率时间系列处利用静态C臂位置来采集。大多数方法然后使用光学流方案来逐帧导出指示对比度改变的幅度和方向的向量。虽然2D方法已被广泛地讨论，但是它们受到重叠的脉管系统的挑战，并且典型地在存在植入物的情况下失效。
发明内容
本文描述了用于促进流数据的重构的技术。依照一个方面，框架接收四维投影图像数据集并且配准图像数据集中的一对或多对时间上相邻的投影图像。二维流（2D）图可以基于经配准的对而确定。2D流图的每一个像素可以包括表示流幅度和二维流方向（u’, v’）的第一流向量。框架然后可以根据心相对2D流图分类，并且基于经分类的流图重构三维（3D）流图。3D流图中的每一个体元可以包括表示流幅度和3D流方向（x’, y’, z’）的流向量。
提供该概述来以简化的形式引入以下在随后的详细描述中进一步描述的概念的选择。不意图标识所要求保护的主题的特征或本质特征，也不意图将其用于限制所要求保护的主题的范围。另外，所要求保护的主题不限于解决本公开的任何部分中指出的任何或全部缺点的实现方式。
附图说明
本公开及其许多随附方面的更加完整的领会将容易被获取，因为该领会在结合随附各图来考虑时通过参照以下详细描述变得更好理解。另外，应当指出的是，相同的数字贯穿各图用于索引相同的元素和特征。
图1是图示了示例性成像系统的框图；
图2示出重构流数据的示例性方法；
图3示出如由颈内动脉供给的患者的大脑脉管系统的三维（3D）血管造影图像（或4D DSA）的时间系列；
图4a示出两个连续时间点处的投影图像Pt和Pt+1；
图4b示出在其已经配准到体积V4Dt之后的体积V4Dt+1；
图5图示了心电图（ECG）追踪信号如何用于将2D流图分类到相位仓（bin）中；
图6示出在颈内动脉的根部处可辨认的示例性时间-对比剂浓度（TCC）曲线；
图7是基于四维（4D）DSA数据执行的相位装仓的示意图；以及
图8图示了示例性层析重构技术。
具体实施方式
在以下描述中，阐述众多具体细节，诸如具体组件、设备、方法等的示例，以便提供本发明的实施例的透彻理解。然而，对本领域技术人员将明显的是，这些具体细节不需要用于实践本发明的实施例。在其它实例中，并未详细描述公知的材料或方法，以便避免不必要地使本发明的实施例模糊。在本发明容许各种修改和可替换形式的同时，其具体实施例作为示例在附图中被示出并且将在本文中被详细描述。然而，应当理解到，不存在将本发明限制到所公开的特定形式的意图，而是相反，本发明要覆盖落在本发明的精神和范围内的所有修改、等同物和可替换方案。
如本文所使用的术语“x射线图像”可以意指可见x射线图像（例如显示在视频屏幕上）或x射线图像的数字表示（例如对应于x射线检测器的像素输出的文件）。如本文所使用的术语“治疗中x射线图像”可以是指在任何介入过程的治疗实施阶段期间的任何时间点处捕获的图像，所述时间点可以包括辐射源接通或关断的时间。有时，为了方便描述，锥束CT（CBCT）成像数据可以在本文中用作示例性成像模态。然而，将领会到，来自包括但不限于X射线放射照相术、MRI、CT、PET（正电子发射断层摄影术）、PET-CT、SPECT、SPECT-CT、MR-PET、3D超声图像等的任何类型的成像模态的数据也可以使用在本发明的各种实施例中。
除非另行陈述，否则如从以下讨论明显的，将领会到，诸如“分段”、“生成”、“配准”、“确定”、“对准”、“定位”、“处理”、“计算”、“选择”、“估计”、“检测”、“追踪”等之类的术语可以是指计算机系统或类似的电子计算设备的动作和过程，所述动作和过程将表示为计算机系统的寄存器和存储器内的物理（例如电子）量的数据操纵和变换成类似地表示为计算机系统的存储器或寄存器或其它这样的信息存储、传输或显示设备内的物理量的其它数据。本文所描述的方法的实施例可以使用计算机软件来实现。如果以遵从公认标准的编程语言编写，设计成实现该方法的指令序列可以被编译以用于在各种硬件平台上执行和用于对接到各种操作系统。此外，没有参照任何特定编程语言来描述本发明的实施例。将领会到，各种编程语言可以用于实现本发明的实施例。
如本文所使用的，术语“图像”或“图像数据”是指包括离散图像元素（例如针对2D图像的像素和针对3D图像的体元）的多维数据。图像可以是例如通过（锥束）计算机断层摄影术、磁共振成像、超声或本领域技术人员已知的任何其它医学成像系统收集的主体的医学图像。图像还可以从非医学情境（诸如例如遥感系统、电子显微镜等）提供。尽管图像可以被认为是从R2到R4或R8的函数，但是本发明的方法不限于这样的图像，并且可以应用于任何维度的图像，例如2D图片或3D体积。对于2维或3维图像，图像域典型地为2维或3维矩形阵列，其中每一个像素或体元可以参考一组2或3个相互正交的轴来寻址。如本文所使用的术语“数字”和“数字化的”将酌情是指以经由数字采集系统或经由从模拟图像的转换采集的数字或数字化格式的图像或体积。
为了缓解由传统的2D方法遭遇的问题，本框架以生成血液/对比剂（contrast）混合物流方向的时变3D表示的旋转方式采集DSA数据。以其典型形式，这样的4D-DSA数据仅用作用于视觉检查的定性工具。然而，所合并的对比剂动态性允许评价定量度量，诸如针对给定心相的对比剂/血液流混合物的方向。
用于导出2D流向量的传统方法典型地使用光学流的方案，合并投影视图中的空间对比度改变和时间帧之间的对比度改变二者，从而导致针对每一个投影图像像素的2D流向量。如这样应用，该方法将由于旋转采集中的投影帧之间的角度差异（高达1.5度）而失败。
为了合并来自旋转采集的投影图像的3D流向量可视化，来自2D连续投影帧的流向量可以在相位装仓重构之前首先量化。本框架的一个方面从现有4D-DSA数据集生成心相特定对比剂/血液混合物3D流图。3D流图中的每一个体元可以包括表示3D流方向（x’, y’,z’）和流幅度的流向量。首先登记4D-DSA数据集中的时间上相邻的投影图像以计及角度差异。然后可以执行基于投影的流计算以确定所登记的投影图像之间的流向量。3D流图可以基于流向量而重构。本文将更加详细地描述这些和其它示例性特征和优点。
图1是图示了示例性成像系统100的框图。成像系统100包括用于实现如本文所描述的框架的计算机系统101。计算机系统101可以进一步通过有线或无线网络连接到成像设备102和工作站103。成像设备102可以是放射学扫描仪，诸如磁共振（MR）扫描仪、X射线或CT扫描仪。在一些实现方式中，成像设备102采用锥束CT（或C臂CT、锥束体积CT、平板CT等）成像技术来采集从患者的感兴趣的解剖体或结构的不同视图重构的体积或3D数据集。成像设备102可以包括安装在例如C臂上的扫描仪，C臂绕患者的身体旋转以采集表示多个视图的不同图像。成像设备102的不同实现方式可以包括例如固定房间C臂、移动U臂、移动O臂、移动C臂等等。
计算机系统101可以是台式个人计算机、便携式膝上型计算机、另一便携式设备、迷你计算机、大型计算机、服务器、存储系统、专用数字器具或具有配置成存储一批数字数据项的存储子系统的另一设备。在一个实现方式中，计算机系统101包括经由输入-输出接口121耦合到一个或多个非暂时性计算机可读介质105（例如计算机储存或存储器）、输出设备108（例如监视器、显示器、打印机等）和各种输入设备110（例如鼠标、键盘、触摸垫、语音识别模块等）的处理器或中央处理单元（CPU）104。计算机系统101还可以包括支持电路，诸如高速缓存、电源、时钟电路和通信总线。再进一步地，计算机系统101可以提供有图形控制器芯片，诸如支持高性能图形功能的图形处理单元（GPU）。
要理解的是，本技术可以实现在各种形式的硬件、软件、固件、专用处理器或其组合中。在一个实现方式中，本文所描述的技术可以由重构单元106实现。重构单元106可以包括有形地体现在非暂时性计算机可读介质105中的计算机可读程序代码。非暂时性计算机可读介质105可以包括随机存取存储器（RAM）、只读存储器（ROM）、磁性软盘、闪速存储器和其它类型的存储器或其组合。计算机可读程序代码由CPU 104执行以控制和/或处理来自成像设备102的图像数据。
照此，计算机系统101是通用计算机系统，其在执行计算机可读程序代码时成为专用计算机系统。计算机可读程序代码不意图受限于任何特定的编程语言及其实现。将领会到，各种编程语言及其编码可以用于实现本文所包含的本公开的教导。计算机系统101还可以包括操作系统和微指令代码。本文所描述的各种技术可以实现为微指令代码的部分或实现为应用程序或软件产品的部分或其组合，该技术经由操作系统被执行。各种其它的外围设备，诸如附加的数据存储设备和打印设备，可以连接到计算机系统101。
工作站103可以包括计算机和适当的外设（诸如键盘和显示器），并且可以结合整个系统100来操作。例如，工作站103可以与成像设备102通信，使得由成像设备102收集的图像数据可以在工作站103处再现并且在显示器上查看。工作站103可以包括用户接口，其允许放射学家或任何其他有技能的用户（例如医生、技术人员、操作者、科学家等）操纵图像数据。例如，用户可以标识图像数据中的感兴趣的结构或区，或者经由用户接口使用预先限定的描述符来注释感兴趣的结构或区。另外，工作站103可以直接与计算机系统101通信以显示经处理或重构的图像数据。例如，放射学家可以交互式地操纵经处理的图像数据的所显示的表示，并且从各种视点和在各种读取模式中查看它。
图2示出重构流数据的示例性方法200。示例性方法200可以由之前参照图1描述的计算机系统101中的重构单元106实现。应当指出的是，在以下讨论中，将使用相同的标号来引用图1中描述的特征。
在202处，重构单元106基于时变图像数据执行4D DSA重构以生成感兴趣的对象的4D（或时间分辨的3D）DSA数据集V4D。4D DSA数据集V4D表示在感兴趣的对象的血管状结构中注入的对比剂流或动态性的相位信息。4D DSA数据集V4D可以包括允许在其经过期间的任何时间处和在任何期望的查看（或投影）角度处穿过血管状结构的对比剂剂量的查看的3D血管造影图像体积的时间系列。可以从4D DSA数据集V4D提取一组时间-对比剂浓度（TCC）曲线。每一个TCC曲线提供针对给定图像元素（例如针对2D数据的像素、针对3D体积的体元）的关于对比剂密度和进而心相的信息，并且通过将各个投影图像重投影到体积域中来生成。
感兴趣的对象可以是为了调查或检查所标识的任何生物对象，诸如患者或主体的脑部、心脏、腿部、手臂等等的部分。感兴趣的对象包括一个或多个血管状结构（例如血管、动脉、血管树或网络等）。所述一个或多个血管状结构可以是动态的或时变的结构，其可以填充有对比剂或介质以用于观察其随时间的传播。在一些实现方式中，还重构在感兴趣的对象中植入的设备（例如流分配设备）的静态（即无时态）3D图像数据。
时变图像数据可以是通过使用成像设备102执行旋转扫描或角度采集而采集的一组数字化的2D DSA投影图像。可以利用X射线延迟经由成像设备102采集被减投影图像的序列以捕获尽可能多的具有混浊的心相。更特别地，可以首先经由成像设备102采集掩模图像数据集，使得可以从对应的时变的对比剂填充的投影图像数据集减去它。掩模图像简单地为在配给对比剂（或介质）以填充要调查的感兴趣的被辐照对象的血管状结构之前的相同区域的图像。实际的角度和时变的2D投影数据可以基于在将对比介质注入到血管状结构之前或之后、当对比剂的首次流入变得可见时发起的对比剂增强采集。掩模和填充行进都可以遵循相同的采集轨迹。轨迹可以覆盖3D DSA的整个视场（FOV）范围。
成像设备102可以是具有单个、两个或多个成像平面的扫描仪或C臂系统。例如，成像设备102可以是基于平板的X射线扫描仪，其包括至少一对X射线源和X射线检测器。可替换地，成像设备102可以包括覆盖至少一对X射线源和X射线检测器的旋转CT机架。在其它实现方式中，成像设备102是MR投影扫描仪。在再其它的实现方式中，成像设备102是覆盖至少一对光源和光学检测器的旋转光学CT机架。还可以使用其它类型的成像设备102，诸如角度采样超声。可以在每一个角度位置处采集单个图像。可替换地，在每一个角度位置处采集具有时间延迟的两个图像，以便避免必须执行随后的图像配准。
图3示出如由颈内动脉供给的患者的大脑脉管系统的三维（3D）血管造影图像（或4D DSA）302的时间系列。基于通过单个旋转成像平面采集的时变图像数据而重构图像302。在2014年6月12日提交的美国申请序列号14/302,596（现在为美国公开号2014/0376791）中描述用于执行由单个旋转平面C臂系统采集的时变图像数据的4D-DSA重构的方法，所述美国申请通过引用并入于此。这些方法确定血管状结构的时变体积衰减曲线，从而导致包括时间维度的3D加时间（或4D-DSA）体积数据集。4D-DSA数据集还可以从时变和投影角度变化数据导出。可以在执行时间分辨的3D DSA的内插中使用置信度值或曲线。这样的框架可以被应用一次，或者以迭代的方式应用。4D-DSA数据集还可以基于例如从在多个角度处采集的时变2D投影数据导出的初始时变3D投影数据集而动态地并且迭代地重构。
在2015年12月3日提交的题为“Tomography system and method for generatinga sequence of volume images of a vasculature”的德国申请号102015224176.9（还是2015年12月9日提交的PCT申请号PCT/EP2015/079102）中描述用于执行通过双C臂系统采集的时变图像数据的4D-DSA重构的方法，所述德国申请通过引用并入于此。这些技术是基于包括两个同时旋转的平面的血管造影双平面系统。可以明显改进时间分辨的体积的经重构的系列的精度，因为来自这两个平面的信息可以用于减轻由于血管重叠所致的精度问题。
返回到图2，在204处，重构单元106基于4D DSA图像数据配准时间上相邻（或连续采集）的旋转投影图像对。每一对时间上相邻的投影图像可以通过Pt和Pt+1来指代。为了执行配准，重构单元106可以首先标识对应于每一个投影图像对Pt和Pt+1的V4D体积（V4Dt和V4Dt+1）。接着，重构单元106可以通过确定从与投影图像相关联的已知系统几何结构的投影位移来配准投影图像Pt与投影图像Pt+1。更特别地，对于每一个投影图像，存在唯一地描述系统的采集几何结构G（源和检测器位置）。对于两个相邻的投影图像Pt和Pt+1及其相应的系统几何结构Gt和Gt+1，可以导出变形场以配准体积V4Dt+1与体积V4Dt。该变形场D(t+1►t)然后可以被正向投影到投影图像Pt+1中以生成2D投影位移（或变形向量场）d(t+1►t)。这样的投影位移可以应用于投影图像Pt以生成经配准的投影图像Pt+1►t。
可替换地，重构单元106可以通过使用V4Dt的采集几何结构正向投影V4Dt+1来配准相邻的旋转投影图像Pt和Pt+1以创建经配置的图像Pt+1►t。该示例性方法通过图4a-b来图示。图4a示出具有它们之间的小角度增量的两个连续时间点（t和t+1）处的投影图像Pt和Pt+1。对比剂在采集之间前进，如通过如由体积V4Dt和V4Dt+1捕获的结构404a-b的不同填充状态可见的。图4b示出在其已经配准（或平移）到体积V4Dt之后的体积V4Dt+1。使用V4Dt的采集几何结构的体积V4Dt+1->t的正向投影导致经配准的投影图像对410（Pt和Pt+1->t）。可以通过经由逐帧投影方法合并来自每一个体积的对比剂中的改变来生成T配准矩阵。
在208处，单元106基于相邻配准的投影图像生成2D流图。2D流图的每一个像素表示投影图像的特定2D位置（u, v）处的流向量。可以处理对应于不同投影角度的N个投影图像以生成T个2D流图，其中T = N-1，因为2个投影图像用于计算流图。每一个流向量表示经配准的投影图像Pt和Pt+1之间的对比剂流的2D方向（u’, v’），以及对比剂流的幅度。这样的流向量可以通过使用基于投影的流计算方法（例如光学流）来确定。
在210处，重构单元106根据心相对2D流图进行分类。可以使用心电图（ECG）追踪信号或从4D DSA数据集自身导出的底层时间-对比剂浓度（TCC）曲线来执行分类。更特别地，经处理的投影数据（所建立的流向量）的ECG门控或TCC分布可以用于将2D流图分类到相位仓中。分类可以导致包含M个径向分布的流图的N个相位仓。
图5图示了ECG追踪信号504如何用于将2D流图分类到相位仓502中。每一个竖直阴影区域502表示对应于ECG追踪信号的特定相位的相位仓。图6示出可以用于对2D流图进行分类的颈内动脉604的根部处可分辨的示例性TCC 602。TCC 602表示随时间的浓度/图像强度，并且可以从4D-DSA数据集提取。TCC 602可以单独地或者与从ECG信号监视确定的信息组合地使用以提供用于投影图像分类的另外的输入。图7是基于4D DSA数据执行的相位装仓的示意图701。时间-对比剂浓度（TCC）曲线峰702对应于心脏舒张，而谷704对应于心脏收缩。在心脏收缩期间，非混浊血流入到脉管系统中，稀释对比剂并且导致入射x射线的较低衰减。2D流图可以使用由TCC指示的心相而分类到相位仓706中。
返回到图2，在212处，重构单元106基于经分类的2D流图而生成3D流图。3D流图的每一个体元可以表示特定3D位置（x, y, z）和心相处的一个或多个3D流向量。每一个3D流向量表示特定流幅度和3D流方向（x’, y’, z’）。可以在例如工作站103处显示3D流图。这样的3D流图可以用于呈现流改变的详细模式，以用于评价和/或诊断疾病，诸如动脉狭窄或动脉瘤流入。
为了生成3D流图，可以分离地且通过断层摄影术重构每一个心相仓中的2D流图，以得出3D空间中的流向量的所估计的分布。断层摄影术重构一般牵涉通过将系统几何结构考虑在内而将独立于强度或流数据的2D投影图像反向投影到3D空间中。参见例如Feldkamp等人的Practical cone-beam algorithm，J. Opt. Soc. Am.，第1卷，1984，第1612-619页，其通过引用并入于此。在一些实现方式中，通过使用2D流图中的流向量的简单反向投影和在三维空间中累积流向量来执行断层摄影术重构。
图8图示了示例性断层摄影术重构技术。每一个2D投影图801a-c与（u, v）空间相关联，而3D图像803与（x, y, z）空间相关联。可以从分组在特定心相仓中的2D流图提取2D流向量804a-c。可以通过反向投影并且代数组合流向量804a-c来获取3D流向量802。虽然投影域中的每一个2D流向量804a-c仅具有两个方向的分量（u’, v’）和针对2D空间中的特定像素（u, v）的流幅度值，但是从对应于特定心相仓的不同角度位置反向投影这些2D流向量导致具有三个方向的分量（x’, y’, z’）和针对3D空间中的特定体元（x, y, z）的流幅度值的3D流向量802。
总结来说，系统采集几何结构建立2D投影图（801a-c）与3D体积803之间的关系。经反向投影的信息具有多个方向的分量（u, v, u’, v’）和幅度值。每一个图可以分离地反向投影以实现可以组合以得出最终的方向（x’, y’, z’）和幅度值的三个不同的体积。
为了进一步改进流图重构，可以通过在经重构的3D流图上使用正则化矩阵来强制执行所述流的平滑。可以通过迭代优化框架来实现正则化矩阵和3D流图的重构。更特别地，优化框架可以在优化经重构的流图到所导出的2D流图的拟合与在相邻流相位之间应用平滑正则化矩阵之间交替。可以在3D流图的用户限定的子区上执行优化，这允许更快的计算。
尽管已经以特定于结构特征和/或方法学步骤的语言描述了一个或多个以上描述的实现方式，但是要理解到，可以在没有所描述的特定特征或步骤的情况下实践其它实现方式。更确切地，特定特征和步骤作为一个或多个实现方式的优选形式而被公开。

Claims (18)

1.一种体现指令程序的非暂时性计算机可读介质，所述指令程序可由机器执行以执行包括以下的步骤：
接收现有的四维数字减影血管造影术（DSA）数据集；
配准从相应的四维图像体积中导出的时间上相邻的投影图像对，以生成经配准的时间上相邻的投影图像对，所述相应的四维图像体积是从所述现有的四维数字减影血管造影术（DSA）数据集中标识出来的；
基于所述经配准的时间上相邻的投影图像对确定二维流图，其中二维流图的像素包括表示在至少一个所述经配准的时间上相邻的投影图像对中的两个时间上相邻的投影图像之间的对比剂的第一流幅度和二维流方向的对比剂的第一流向量；
根据心相对所述二维流图进行分类以生成经分类的二维流图；以及
通过对所述经分类的二维流图中的流向量进行反向投影并对三维空间中的所述流向量进行代数组合来重构三维流图，其中所述经重构的三维流图的体元包括表示所述对比剂的第二流幅度和三维流方向的所述对比剂的第二流向量。
2.权利要求1所述的非暂时性计算机可读介质，其中指令可由机器执行以按照所述心相执行对二维流图的分类，以通过使用心电图（ECG）追踪信号生成所述经分类的二维流图。
3.权利要求1所述的非暂时性计算机可读介质，其中指令可由机器执行以按照所述心相执行对二维流图的分类，以通过使用从所述现有的四维DSA数据集导出的时间-对比剂浓度（TCC）曲线生成所述经分类的二维流图。
4.一种重构流数据的方法，包括：
接收现有的四维投影图像数据集；
配准从相应的四维图像体积中导出的时间上相邻的投影图像对，以生成经配准的时间上相邻的投影图像对，所述相应的四维图像体积是从所述现有的四维投影图像数据集中标识出来的；
基于所述经配准的时间上相邻的投影图像对确定对比剂的二维流图；
根据心相对所述对比剂的二维流图进行分类以生成所述对比剂的经分类的二维流图；以及
通过对所述经分类的二维流图中的流向量进行反向投影并对三维空间中的所述流向量进行代数组合来重构所述对比剂的三维流图。
5.权利要求4所述的方法，还包括基于二维时变图像数据重构所述现有的四维投影图像数据集。
6.权利要求5所述的方法，还包括由成像设备通过在每一个角度位置处采集单个图像来采集二维时变图像数据。
7.权利要求5所述的方法，还包括由成像设备通过在每一个角度位置处采集具有时间延迟的两个图像来采集二维时变图像数据。
8.权利要求4所述的方法，其中配准所述时间上相邻的投影图像对包括：
基于所述相应的四维图像体积并基于与至少一个所述时间上相邻的投影图像对相关联的系统几何结构确定投影位移；以及
向所述至少一个所述时间上相邻的投影图像对中的一个图像应用投影位移以生成经配准的投影图像。
9.权利要求4所述的方法，其中配准从所述相应的四维图像体积中导出的所述时间上相邻的投影图像对，所述相应的四维图像体积是从所述现有的四维投影图像数据集中标识出来的包括：
使用所述相应的四维图像体积中的第二图像体积的采集几何结构正向投影所述相应的四维图像体积中的第一图像体积以生成经配准的投影图像。
10.权利要求4所述的方法，其中确定所述对比剂的二维流图包括通过使用基于投影的流计算方法来确定流向量。
11.权利要求4所述的方法，其中根据心相对所述对比剂的二维流图进行分类包括使用心电图（ECG）追踪信号对所述对比剂的二维流图进行分类。
12.权利要求4所述的方法，其中根据心相对所述对比剂的二维流图进行分类包括使用从所述现有的四维投影图像数据集导出的时间-对比剂浓度（TCC）曲线对所述对比剂的二维流图进行分类。
13.权利要求4所述的方法，其中重构所述对比剂的三维流图包括分离地且通过断层摄影术重构至少一个所述心相中的所述对比剂的所述经分类的二维流图。
14.权利要求4所述的方法，还包括在所述对比剂的经重构的三维流图上应用正则化矩阵。
15.权利要求14所述的方法，其中应用正则化矩阵包括在所述对比剂的所述经重构的三维流图的用户限定的子区上应用正则化矩阵。
16.一种重构流数据的系统，包括：
用于存储计算机可读程序代码的非暂时性存储器设备；以及
与存储器设备通信的处理器，处理器利用计算机可读程序代码来操作以执行包括以下的操作
接收现有的四维投影图像数据集；
配准从相应的四维图像体积中导出的时间上相邻的投影图像对，以生成经配准的时间上相邻的投影图像对，所述相应的四维图像体积是从所述现有的四维 投影图像数据集中标识出来的；
基于所述经配准的时间上相邻的投影图像对确定对比剂的二维流图；
根据心相对所述对比剂的二维流图进行分类以生成所述对比剂的经分类的二维流图；以及
通过对所述经分类的二维流图中的流向量进行反向投影并对三维空间中的所述流向量进行代数组合来重构所述对比剂的三维流图。
17.权利要求16所述的系统，其中处理器利用计算机可读程序代码来操作以按照所述心相执行对所述对比剂的二维流图的分类，以通过使用心电图（ECG）追踪信号生成所述对比剂的所述经分类的二维流图。
18.权利要求16所述的系统，其中处理器利用计算机可读程序代码来操作以按照所述心相执行对所述对比剂的二维流图的分类，以通过使用从所述现有的四维投影图像数据集导出的时间-对比剂浓度（TCC）曲线生成所述对比剂的所述经分类的二维流图。
CN201710906720.8A 2016-09-30 2017-09-29 流数据的重构 Active CN107886554B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP16191701.8 2016-09-30
EP16191701.8A EP3300664B1 (en) 2016-09-30 2016-09-30 Reconstruction of flow data

Publications (2)

Publication Number Publication Date
CN107886554A CN107886554A (zh) 2018-04-06
CN107886554B true CN107886554B (zh) 2021-11-09

Family

ID=57083167

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710906720.8A Active CN107886554B (zh) 2016-09-30 2017-09-29 流数据的重构

Country Status (4)

Country Link
US (1) US11317875B2 (zh)
EP (1) EP3300664B1 (zh)
JP (1) JP6479919B2 (zh)
CN (1) CN107886554B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102019200888A1 (de) * 2019-01-24 2020-07-30 Siemens Healthcare Gmbh Bestimmen eines Bilddatensatzes
JP7165600B2 (ja) * 2019-02-28 2022-11-04 キヤノンメディカルシステムズ株式会社 Ｘ線診断装置及び医療情報処理装置
CN116867440A (zh) * 2023-04-21 2023-10-10 武汉迈瑞医疗技术研究院有限公司 一种胎心的超声成像方法和系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1620991A (zh) * 2003-11-26 2005-06-01 通用电气公司 心脏显示方法和装置
US20060140482A1 (en) * 2003-06-18 2006-06-29 Thomas Koehler Motion compensated reconstruction technique
US20080205726A1 (en) * 2007-02-23 2008-08-28 Siemens Aktengesellschaft Method for providing extended possibilities when imaging a patient's heart
CN102696056A (zh) * 2009-08-17 2012-09-26 米斯特雷塔医疗有限公司 用于四维血管造影和荧光透视的系统和方法
CN104252714A (zh) * 2013-06-25 2014-12-31 美国西门子医疗解决公司 时变数据的重建

Family Cites Families (39)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5522393A (en) * 1994-05-24 1996-06-04 Duke University Multi-dimensional real-time ultrasonic blood flow imaging apparatus and method
US6246898B1 (en) * 1995-03-28 2001-06-12 Sonometrics Corporation Method for carrying out a medical procedure using a three-dimensional tracking and imaging system
US6447450B1 (en) * 1999-11-02 2002-09-10 Ge Medical Systems Global Technology Company, Llc ECG gated ultrasonic image compounding
US6503202B1 (en) * 2000-06-29 2003-01-07 Acuson Corp. Medical diagnostic ultrasound system and method for flow analysis
US7372984B2 (en) * 2004-05-05 2008-05-13 California Institute Of Technology Four-dimensional imaging of periodically moving objects via post-acquisition synchronization of nongated slice-sequences
US8340744B2 (en) 2005-12-15 2012-12-25 Koninklijke Philips Electronics N.V. System, apparatus, and method for reproducible and comparable flow acquisitions
WO2008085193A2 (en) * 2006-08-14 2008-07-17 University Of Maryland Quantitative real-time 4d strees test analysis
EP2054856B1 (en) * 2006-08-15 2010-10-06 Koninklijke Philips Electronics N.V. Motion compensation in energy-sensitive computed tomography
US9672651B2 (en) * 2006-10-17 2017-06-06 Koninklijke Philips N.V. Four-dimensional reconstruction of regions exhibiting multiple phases of periodic motion
JP5591440B2 (ja) 2007-01-17 2014-09-17 株式会社東芝 医用画像表示装置
US8542900B2 (en) * 2007-03-08 2013-09-24 Sync-Rx Ltd. Automatic reduction of interfering elements from an image stream of a moving organ
US11064964B2 (en) * 2007-03-08 2021-07-20 Sync-Rx, Ltd Determining a characteristic of a lumen by measuring velocity of a contrast agent
US8073224B2 (en) * 2007-07-09 2011-12-06 Siemens Aktiengesellschaft System and method for two-dimensional visualization of temporal phenomena and three dimensional vessel reconstruction
JP2009279290A (ja) 2008-05-26 2009-12-03 Hitachi Medical Corp 医用画像診断装置
JP5613182B2 (ja) * 2009-03-06 2014-10-22 バイオ−ツリー システムズ， インコーポレイテッド 血管分析方法および装置
US8771189B2 (en) * 2009-03-18 2014-07-08 Siemens Medical Solutions Usa, Inc. Valve assessment from medical diagnostic imaging data
US9019305B2 (en) * 2009-04-10 2015-04-28 Siemens Medical Solutions Usa, Inc. Method of visualization of contrast intensity change over time in a DSA image
US20110235885A1 (en) * 2009-08-31 2011-09-29 Siemens Medical Solutions Usa, Inc. System for Providing Digital Subtraction Angiography (DSA) Medical Images
US8731262B2 (en) * 2010-06-03 2014-05-20 Siemens Medical Solutions Usa, Inc. Medical image and vessel characteristic data processing system
WO2012011014A1 (en) * 2010-07-20 2012-01-26 Koninklijke Philips Electronics N.V. 3d flow visualization
US8768031B2 (en) * 2010-10-01 2014-07-01 Mistretta Medical, Llc Time resolved digital subtraction angiography perfusion measurement method, apparatus and system
US8929632B2 (en) * 2011-09-22 2015-01-06 Siemens Aktiengesellschaft Temporal difference encoding for angiographic image sequences
US8948487B2 (en) * 2011-09-28 2015-02-03 Siemens Aktiengesellschaft Non-rigid 2D/3D registration of coronary artery models with live fluoroscopy images
CN103857340A (zh) 2011-10-10 2014-06-11 莫纳什大学 心脏成像方法
US8463012B2 (en) * 2011-10-14 2013-06-11 Siemens Medical Solutions Usa, Inc. System for comparison of medical images
DE102011086771A1 (de) 2011-11-22 2013-05-23 Siemens Aktiengesellschaft Computertomographieanlage und Verfahren zum Ermitteln von Volumeninformationen zu einem Körper
WO2013192234A1 (en) * 2012-06-18 2013-12-27 The Research Foundation For The State University Of New York Systems and methods for identifying historical vasculature cases
US20140071125A1 (en) * 2012-09-11 2014-03-13 The Johns Hopkins University Patient-Specific Segmentation, Analysis, and Modeling from 3-Dimensional Ultrasound Image Data
US10820872B2 (en) * 2012-12-19 2020-11-03 Koninklijke Philips N.V. X-ray controlled contrast agent injection
JP6382036B2 (ja) * 2013-09-30 2018-08-29 キヤノンメディカルシステムズ株式会社 超音波診断装置及び画像処理装置
US9478048B2 (en) * 2014-07-02 2016-10-25 Siemens Medical Solutions Usa, Inc. Prior image based three dimensional imaging
US9962142B2 (en) * 2014-11-14 2018-05-08 The Regents Of The University Of California Ultrasound-based volumetric particle tracking method
US10460452B2 (en) * 2014-11-14 2019-10-29 The Regents Of The University Of California Ultrasound-based volumetric particle tracking method
US20160135775A1 (en) * 2014-11-17 2016-05-19 Wisconsin Alumni Research Foundation System And Method For Time-Resolved, Three-Dimensional Angiography With Physiological Information
US9508157B2 (en) * 2014-12-12 2016-11-29 Siemens Medical Solutions Usa, Inc. Reconstruction of aneurysm wall motion
JP6513413B2 (ja) 2015-02-13 2019-05-15 キヤノンメディカルシステムズ株式会社 医用画像処理装置及び磁気共鳴イメージング装置
US10818073B2 (en) 2015-03-10 2020-10-27 Wisconsin Alumni Research Foundation System and method for time-resolved, three-dimensional angiography with flow information
US10271905B2 (en) * 2015-04-20 2019-04-30 The Johns Hopkins University Patient-specific virtual intervention laboratory to prevent stroke
US10524756B2 (en) * 2015-08-27 2020-01-07 Varian Medical Systems International Methods and systems for image artifacts reduction

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060140482A1 (en) * 2003-06-18 2006-06-29 Thomas Koehler Motion compensated reconstruction technique
CN1620991A (zh) * 2003-11-26 2005-06-01 通用电气公司 心脏显示方法和装置
US20080205726A1 (en) * 2007-02-23 2008-08-28 Siemens Aktengesellschaft Method for providing extended possibilities when imaging a patient's heart
CN102696056A (zh) * 2009-08-17 2012-09-26 米斯特雷塔医疗有限公司 用于四维血管造影和荧光透视的系统和方法
CN104252714A (zh) * 2013-06-25 2014-12-31 美国西门子医疗解决公司 时变数据的重建

Also Published As

Publication number Publication date
US20180092608A1 (en) 2018-04-05
EP3300664B1 (en) 2019-04-17
JP2018057847A (ja) 2018-04-12
US11317875B2 (en) 2022-05-03
CN107886554A (zh) 2018-04-06
EP3300664A1 (en) 2018-04-04
JP6479919B2 (ja) 2019-03-06

Similar Documents

Publication Publication Date Title
US9754390B2 (en) Reconstruction of time-varying data
US9684980B2 (en) Prior image based three dimensional imaging
US8055050B2 (en) Motion compensation in energy-sensitive computed tomography
CN111540025B (zh) 预测用于图像处理的图像
US9508157B2 (en) Reconstruction of aneurysm wall motion
US10083511B2 (en) Angiographic roadmapping mask
US20150279084A1 (en) Guided Noise Reduction with Streak Removal for High Speed C-Arm CT
US20230097849A1 (en) Creation method of trained model, image generation method, and image processing device
Copeland et al. Spatio-temporal data fusion for 3D+ T image reconstruction in cerebral angiography
CN107886554B (zh) 流数据的重构
Chen et al. Prior image constrained compressed sensing (PICCS) and applications in X-ray computed tomography
CN107170021B (zh) 时变数据的细化重构
JP2019528116A (ja) ４ｄコンピュータ断層撮影から呼吸性血液量変化を決定するシステム及び方法
Zhang et al. Dynamic estimation of three‐dimensional cerebrovascular deformation from rotational angiography
Manhart et al. Iterative denoising algorithms for perfusion C-arm CT with a rapid scanning protocol
Manhart et al. Guided noise reduction with streak removal for high speed flat detector CT perfusion
Keil et al. Dynamic cone-beam reconstruction using a variational level set formulation
Manhart et al. Dynamic reconstruction with statistical ray weighting for C-arm CT perfusion imaging
JPWO2019044983A1 (ja) Ｘ線ｃｔ装置および画像生成方法

Legal Events

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