CN103976750B - 用于观察脉管系统以及灌注组织的系统 - Google Patents
用于观察脉管系统以及灌注组织的系统 Download PDFInfo
- Publication number
- CN103976750B CN103976750B CN201310757218.7A CN201310757218A CN103976750B CN 103976750 B CN103976750 B CN 103976750B CN 201310757218 A CN201310757218 A CN 201310757218A CN 103976750 B CN103976750 B CN 103976750B
- Authority
- CN
- China
- Prior art keywords
- image
- ray
- successive
- filtered
- error image
- 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
Links
- 230000002792 vascular Effects 0.000 title abstract description 6
- 238000001914 filtration Methods 0.000 claims abstract description 54
- 239000002872 contrast media Substances 0.000 claims abstract description 34
- 238000000034 method Methods 0.000 claims abstract description 26
- 230000005855 radiation Effects 0.000 description 36
- 238000003384 imaging method Methods 0.000 description 11
- 230000008569 process Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 210000003484 anatomy Anatomy 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 230000010412 perfusion Effects 0.000 description 4
- 210000001367 artery Anatomy 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 238000002583 angiography Methods 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 210000003462 vein Anatomy 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- BUGBHKTXTAQXES-UHFFFAOYSA-N Selenium Chemical compound [Se] BUGBHKTXTAQXES-UHFFFAOYSA-N 0.000 description 1
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000001195 anabolic effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000010253 intravenous injection Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000005622 photoelectricity Effects 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 230000008521 reorganization Effects 0.000 description 1
- 229910052711 selenium Inorganic materials 0.000 description 1
- 239000011669 selenium Substances 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
-
- 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/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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- 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
-
- 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/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20224—Image subtraction
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Vascular Medicine (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
本发明涉及用于观察脉管系统以及灌注组织的系统。一种系统和方法包括接收患者体积的掩模x射线图像，接收包括造影剂的患者体积的多个相继x射线图像，从所述多个相继x射线图像的每一个中减去所述掩模x射线图像以生成多个相继x射线差值图像，基于一个或多个滤波参数对所述多个相继x射线差值图像的每一个进行滤波以生成多个相继已滤波x射线差值图像，基于权重将所述多个相继已滤波x射线差值图像的每一个与所述多个相继x射线差值图像中的对应的一个组合以生成多个组合相继x射线图像，以及相继地显示所述多个组合相继x射线图像。
Description
相关申请的交叉引用
本申请涉及并要求享有2012年11月14日提交的美国专利申请序列号61/726035的优先权，出于所有目的，所述文献的内容由此通过引用被全文引入。技术领域
下述的实施涉及对在造影剂存在于患者体积内时采集的血管造影x射线图像的处理。
背景技术
根据常规血管造影x射线成像，使用造影剂来增强患者解剖结构内的血液载送结构的对比度。例如，将造影剂引入患者某体积内(例如，通过静脉注射)，并且当造影剂位于该体积内时采集该体积的x射线图像。在x射线图像中，包含造影剂的结构显得比它们以其它方式看起来要更暗。
根据DSA(数字减影血管造影)，从如上所述那样采集的x射线图像中减去患者体积的“掩模图像”。掩模图像是在不存在造影剂的情况下采集的并表示背景解剖细节。结果得到的图像旨在仅描绘包括造影剂的患者体积的血管和灌注组织组成部分。
需要对患者体积的不同血液载送部分提供有效图示的系统。
附图说明
通过考虑如随附附图中所述的下列说明，实施例的构造和使用将变得显而易见，其中相同的附图标记指示相似的部件，并且其中：
图1示出了根据一些实施例的系统；
图2是根据一些实施例的过程的流程图；
图3示出了根据一些实施例的x射线图像；
图4示出了根据一些实施例的差值x射线图像；
图5示出了根据一些实施例的用户界面；
图6示出了根据一些实施例的已滤波x射线图像；
图7示出了根据一些实施例的用户界面；
图8示出了根据一些实施例的已滤波x射线图像；
图9示出了根据一些实施例的用户界面；
图10示出了根据一些实施例的组合x射线图像；
图11示出了根据一些实施例的用户界面；
图12示出了根据一些实施例的组合x射线图像；
图13示出了根据一些实施例的用户界面；以及
图14是根据一些实施例的系统的功能框图。
具体实施方式
一些实施例提供未由当前x射线成像系统提供的不包含毛细血管中的灌注的大脉管系统或者不包含脉管系统的毛细血管灌注和／或它们的选择性混合的显不。
提供下列说明书以使本领域任何技术人员能够实现和使用所述实施例，并且阐述了预期用于实施所述实施例的最佳方式。然而，对于本领域技术人员而言，各种修改将仍是显而易见的。
图1示出了根据一些实施例的系统1。系统1包括x射线成像系统10、控制和处理系统20以及操作者终端30。通常，以及根据一些实施例，x射线成像系统10将造影剂引入患者体积中，并且采集该患者体积的x射线图像。控制和处理系统20控制x射线成像系统10，并且从其接收所采集的图像。控制和处理系统20如下所述那样处理图像，并且将已处理的图像提供至终端30，由此以供显示。这样的处理可以基于由终端30接收并且由终端30提供至控制和处理系统20的用户输入。
x射线成像系统10包括C形臂11，C形臂11上安装有辐射源12和辐射检测器13。C形臂11安装在支撑装置14上，并且被配置为相对于支撑装置14而顺时针或逆时针平移。该平移使得辐射源12和辐射检测器13围绕中心体积旋转，同时维持它们之间的物理关系。实施例不局限于基于C形臂的成像系统。
辐射源12可以包括任意合适的辐射源，包括但不局限于Diabo1oTMx射线管。在一些实施例中，辐射源12发射能量范围在50至150keV内的电子、光子或其他类型的辐射。
辐射检测器13可以包括用以基于接收的x射线辐射采集图像的任意系统。在一些实施例中，辐射检测器13是平板成像设备，其使用闪烁体层和布置成二维阵列的固态无定形硅光电二极管。闪烁体层接受光子，并且与接收到的光子的强度成比例地产生光。光电二极管的阵列接收光，并且将接收的光的强度记录为存储的电荷。
在其他实施例中，辐射检测器13将接收的光子转换成电荷，而无需闪烁体层。光子被无定形硒光电导体阵列直接吸收。光电导体将光子直接转换成存储的电荷。辐射检测器13可以包括CCD或者基于管的照相机，包括其内设置有闪烁体、反射镜和照相机的防光外壳。
由辐射检测器13转换和存储的电荷表示由从辐射源12发射的x射线产生的辐射场的每个位置处的辐射强度。辐射场的特定位置处的辐射强度表示沿着辐射源12和辐射场的该特定位置之间的发散线的组织的衰减属性。因而，由辐射检测器13采集的辐射强度集可以表示这些组织的二维投影图像。
造影剂注射器17可以包括任意已知设备或者适于可控制地将造影剂引入患者体积内的设备。如上所述，包含造影剂的结构在x射线图像中显得比它们以其他方式看起来更暗。造影剂注射器17可以包括用于一种或多种造影剂的每一种的储存器，以及诸如以中空针为终端的医用级管状物的患者接口。
系统20可以包括任意通用或专用计算系统。因此，系统20包括被配置为执行处理器可执行的程序代码以使得系统20如本文所述的那样工作的一个或多个处理器21，以及用于存储程序代码的存储设备22。存储设备22可以包括一个或多个硬盘、固态随机存取存储器以及安装在相应接口(例如USB端口)中的可移除介质(例如，指状储存器)。
存储设备22存储系统控制程序23的程序代码。一个或多个处理器21可以执行系统控制程序23以移动C形臂14，使得辐射源12发射辐射，控制检测器13采集图像，使注射器17将造影剂引入患者15的体积内，并且执行任何其他功能。在这一点上，系统20包括用于与系统10通信的x射线系统接口24以及造影剂注射器接口25。
从系统10采集的图像在数据存储设备22中被存储为DICOM或另一数据格式的已采集图像26。每个已采集图像26可进一步与其采集的细节相关联，所述采集的细节包括但不局限于成像平面位置和角度、成像位置、辐射源至检测器的距离、被成像患者的解剖结构、患者位置、造影剂推射注射剖面、x射线管电压、图像分辨率和辐射剂量。
处理器21可以执行系统控制程序23以处理已采集图像26，从而产生已处理图像27。已处理图像27可以通过系统20的UI界面28而被提供至终端30。UI界面28还可以接收来自终端30的输入，这用于如下所述地控制对已采集图像26的处理。
终端30可以仅包括耦合至系统20的显示设备和输入设备。终端30显示从系统20接收的已处理图像27，并且接收指定滤波参数和／或权重的用户输入。滤波参数和／或权重被传送至系统20，并由此被用于产生新的已处理图像27，以便随后由终端30显示。在一些实施例中，终端30是独立的计算设备，诸如但不局限于，桌面型计算机、膝上型计算机、平板计算机以及智能手机。
系统10、系统20和终端30中的每个可以包括它们的工作所必需的其他元件，以及用于提供除本文所述的那些功能之外的功能的附加元件。
根据所图示的实施例，系统20控制系统10的元件。系统20还处理从系统10接收的图像。而且，系统20从终端30接收输入，并且将已处理图像提供至终端30。实施例不局限于执行这些功能中的每一个功能的单一系统。例如，系统10可以由专用控制系统来控制，其中已采集图像通过计算机网络或者经由物理存储介质(例如DVD)被提供至独立的图像处理系统。
在根据一些实施例的操作的一个示例中，图2是根据一些实施例的过程200的流程图。可以使用硬件、软件或手动装置的任意合适组合来执行过程200和本文所述的其他过程。体现这些过程的软件可以由任何非瞬时有形介质来存储，包括硬盘、软盘、CD、DVD、闪盘或磁带。下文将针对系统1的元件来描述这些过程的示例，但是实施例不局限于此。
首先，在S210，采集患者体积的掩模x射线图像。如本文背景技术部分所述的，掩模x射线图像在将造影剂引入患者体积内之前被采集并且包括患者体积的解剖结构。接着，在S220，接收患者体积的多个相继x射线图像，其中患者体积包括造影剂。
根据一些实施例，依次执行S210和S220，同时维持患者处于特定的成像位置。例如，以及参考系统1的元件，患者15定位在平台16上，以将患者15的特定体积布置在辐射源12和辐射检测器13之间。可以调节平台16以根据需要协助定位该患者体积。随后，系统20指示系统10移动C形臂11，使得辐射源12和辐射检测器13将从期望的投影角度生成患者体积的图像。
系统20还指示造影剂注射器17将造影剂引入患者体积中。辐射源11由高能发生器供能，以在造影剂被引入患者体积中之前、其间和之后以连贯的间隔朝辐射检测器13发射x射线辐射。正如本领域中已知的，造影剂引入的参数(例如流速、位置、体积)以及x射线辐射发射的参数(例如，x射线管电压、剂量)可以由系统控制程序23来控制。辐射检测器13接收发射的辐射，并且对于每个间隔产生一组数据(即投影图像)。间隔在持续时间方面无需相同。
S210的掩模图像是在将造影剂引入患者体积中之前形成的投影图像之一，并且S220的多个相继x射线图像中的每一个是在患者体积包括造影剂的同时产生的投影图像。
图3包括在S210和S220接收的x射线图像310-380的图示。图像310-380是时间上连贯的图像，其中与较大参考数字相关联的图像是在与较小参考数字相关联的图像之后采集的。如上所述，x射线图像310-380可能已经由独立的成像系统产生，并且由执行过程200的、独立且可能位于远程的图像处理系统在S210和S220处接收。
根据图3的示例，在将造影剂引入被成像患者体积中之前采集x射线图像310和320。因而，x射线图像310和320仅包括患者解剖结构的图像。图像330和340示出了造影剂到患者体积内的动脉中的开始引入，其中黑线385表示装有造影剂的注射设备(例如，针)。图像350示出了造影剂从动脉到毛细血管灌注中的流动，图像360示出了毛细血管灌注，图像370示出了造影剂从毛细血管灌注到静脉中的流动，而图像380示出了清除造影剂之后的患者体积。
返回过程200，在S230，从多个相继x射线图像的每一个中减去在S210接收的掩模图像。结果得到的图像在本文中称为多个相继x射线差值图像。正如本领域中已知的，可以例如通过将图像320再次与所述多个相继x射线图像的每一个配准并从中减去相关联的已配准图像来执行该减法。图4示出了在S230基于掩模图像320和图3的相继x射线图像330-380而产生的多个相继x射线差值图像430-480。
接着，在S240，起初相继地显示所述多个相继x射线差值图像430-480。图5示出了根据一些实施例的用于显示多个相继x射线差值图像的界面500。界面500可以响应于由处理器21进行的系统控制程序23的程序代码的执行而由终端30的显示设备来显示。
界面500包括用于显示x射线图像的区域510。根据一些实施例，所显示的图像可以是由成像系统10采集的图像和／或基于已采集图像生成的x射线图像。图5的区域510示出了图4的x射线图像450。
界面500还包括用于显示已知或变得已知的任意合适信息的区域520。例如，在界面500在图像采集期间被显示给操作者的情况中，区域520可以呈现一个或多个图像采集参数、患者生物统计数据等。区域520可以另外或者作为替换显示与区域510中显示的图像相关的信息，诸如但不局限于柱状图、滤波参数和基于图像确定的生物信息。
区域520可以另外或替换地包括一个或多个用户界面控件。这些控件可以允许操作者改变图像采集参数，改变界面500显示的信息，控制由系统20执行的图像处理，或者执行任何其他合适的操作。
在这点上，在一些实施例中，操作者操作控件530，以将滤波参数和权重传递至诸如系统20的图像处理系统。控件530包括滑块540和滑条550。操作者可以使用触屏、鼠标、键盘和任意其他输入设备沿着滑条550移动滑块540。正如下文将详细描述的，与滑块540的每个位置相关联的滤波参数和权重被用于生成图像以供在区域510中显示。
实施例不局限于上文针对图5讨论的内容和布置。任意一个或多个用户界面控件可以用于改变本文讨论的滤波参数和权重。这种控件包括但不局限于文本输入框、标度盘、姿势等。而且，在一些实施例中，操作者可以经由输入设备(例如操纵杆、鼠标、触摸板)、而不使用任何可视觉操作的用户界面控件来将滤波参数和权重传递至图像处理系统。在这样的实施例中，如此改变的参数和／或权重的唯一视觉反馈可以是结果得到的图像的显示序列的改变。
返回S240，将假设图像430至480在区域510中被相继显示。根据一些实施例，可以在图像430至480显示之前对其施加已知的处理，以便于增强边缘、调节亮度、校准视野和／或使得图像符合显示设备的显示属性。随着流程通过S250而继续，图像可以以连续循环被显示。
在S250，确定滤波参数和／或权重是否已经改变。根据本示例，可应用的滤波参数和权重起初为“零”，正如滑块540在滑条550上的中心位置所表示的。当滑块540维持在该位置时，随着流程在S250和S260之间循环，图像的相继显示继续。
为了示出一些实施例的目的，现在假设操作者使用输入设备如箭头560所示的那样将滑块540移动至位置570。系统20接收该移动的指示并且确定滤波参数和权重中的至少一个已经改变。因此，流程从S250进行至S270。在一些实施例中，在执行S270和S280期间，当前显示的图像序列继续被显示，直至如下所述在S290新的图像序列准备好被显示。一些实施例允许独立于过程200的操作而停止连续显示。
在S270，基于一个或多个滤波参数而对所述多个相继x射线差值图像的每一个进行滤波。结果得到的图像在本文中被称为多个相继已滤波x射线差值图像，以便于将这些图像与本文讨论的其他图像区分开来。
在一些实施例中，中心位置右侧的所有位置与特定带通滤波器相关联，而中心位置左侧的所有位置与特定的带阻滤波器相关联。特定的带通滤波器旨在通过与血管相关联的图像数据。在一些实施例中，这种带通滤波器保持尺寸在3至40像素之间的图像数据结构。相反地，与中心位置左侧的位置相关联的特定带阻滤波器旨在抑制尺寸在3至40像素之间的图像数据结构，以便于保留和加强灌注组织。
因而，响应于滑块540到位置570的滑动，基于带通滤波器对所述多个相继x射线差值图像的每一个进行滤波。图6示出了相继已滤波x射线差值图像630-680，它们是根据本示例基于图像430-480和带通滤波器而生成的。如所示，已经从图像630-680中过滤出图像430-480的小结构。
在S280，基于权重将S270处生成的每个图像与S230处生成的图像中的相应的一个图像相组合。根据一些实施例，权重指示每个图像对结果得到的组合图像贡献的量。例如，权重0.75可能指示，当在S280处将图像630与图像430组合以生成相应的组合图像时，图像630的权重将是相应的图像430的权重的三倍(即，0.75比0.25)大。在S280处可以利用任意合适的基于权重的组合算法。
在本示例中，位置570指示权重为1。因而，当在S280生成结果得到的组合图像时，图像630-680被100％地加权，而图像430-480被0％地加权。即，组合图像等同于图像630-680。在S290相继显示所述组合图像，并且流程返回至S250以等待滤波参数和／或权重的下一改变。
在S280之后和在S290显示组合图像之前，可以向组合图像应用已知的处理。如上所述，这样的处理可以增强边缘、调节亮度、校准视野和／或使图像符合显示设备的显示属性。
图7示出了根据本示例的在流程在S250和S260之间循环时组合图像之一的显示。如上提及的，相继地且连续地显示的组合图像与图像630-680相同，因而图7示出了图像660的显示。
现在将假设，操作者在箭头710的方向上将滑块540移动至位置720。流程因而从S250进行至S270。本示例的位置720对应于带阻过滤器和权重1。因而，在S270，基于带阻滤波器对S230处生成的多个相继x射线差值图像的每一个进行滤波，以生成新的多个相继已滤波x射线差值图像。
图8示出了根据本示例的基于图像430-480和带阻滤波器生成的相继已滤波x射线差值图像830-880。如所示，已从图像830-880中过滤出图像430-480的较大结构。
由于位置720指示了权重1，因此图像830-880被100％地加权，而图像430-480被0％地加权，因而在S280生成的组合图像等同于图像830-880。这些组合图像在S290被相继显示，并且流程返回S250以等待滤波参数和／或权重的下一改变。
图9示出了根据本示例的在流程返回S250并且在S250和S260之间循环时组合图像中的图像860的显示。如所示，操作者现在沿着箭头910的方向将滑块540移动至位置920。
由于位置920位于滑条550的左侧，因此在S270基于带阳滤波器对在S230处生成的所述多个相继x射线差值图像的每一个进行滤波，以生成多个相继已滤波x射线差值图像。根据本示例，带阻滤波器与前述讨论的相同，因此在S270处生成的图像是图8的图像830-880。根据一些实施例，已滤波图像830-880和630-680可以在他们的最初生成之后被高速缓存，以便于消除重新生成这些图像的需要。
流程随后进行至S280。与上述示例不同，位置920对应于0至1之间的权重。因而，在S280处生成的每个组合图像将包括其对应的已滤波和未滤波图像的部分。例如，假设位置920对应于权重0.5。因此，在S280，将把在S270处生成的已滤波图像以及S230的差值图像均等地组合以生成多个组合的相继x射线图像。
图10示出了在S280处基于图像430-480和图像830-880生成的组合相继已滤波x射线图像1030-1080，而图11示出了在S290处且在流程返回S250并在S250和S260之间循环时的组合图像的显示。
图11示出了滑块540沿着箭头1110的方向到位置1120的移动。位置1120位于滑条550的右侧，因而在S270基于上述的带通滤波器对在S230处生成的所述多个相继x射线差值图像(例如图像430-480)的的每一个进行滤波，以生成多个相继已滤波x射线差值图像。替换地，从高速缓存中取回这些已经生成的已滤波x射线差值图像(即图像630-680)。
接着，在S280，使用对应于位置1120的权重(例如，0.75)来将S270处生成的已滤波图像与S230的差值图像中的对应的差值图像相组合。在本示例中，在组合过程中，已滤波图像将被赋予是差值图像的权重的3倍(即，0.75比0.25)的权重。
图12示出了组合的相继已滤波x射线图像1230-1280，根据一些实施例，它们是在S280基于图像430-480和图像630-680被生成的。在S290相继显示这些组合图像。图13示出了该显示，随着流程返回S250并且在S250和S260之间循环，该显示在连续的循环中继续。
一些实施例因而有效地提供未由当前x射线成像系统提供的对不包含毛细血管中的灌注的大脉管、不包含脉管的毛细血管灌注和／或它们的选择性混合的逐步改变的显示。
实施例不局限于带通滤波器和带阻滤波器。而且，虽然本示例将相同的带通滤波器与滑块540的所有中心右侧位置相关联并将相同的带阻滤波器与所有中心左侧位置相关联，但是一些实施例还允许操作者修改滤波器通过或阻止的波段。可以在不修改S280处使用的权重的情况下进行这样的修改。即，操作者可以使用合适的用户界面控件仅调节所通过／所阻止的波段的尺寸和／或位置(在不改变权重的情况下)，由此使得生成和显示如上所述的新的组合相继x射线图像。
图14是根据一些实施例的系统1400的功能框图。根据一些实施例，系统1400的硬件和／或软件实施可以执行处理200。
例如，在S210接收掩模图像1405，并且在S220接收N个相继的“造影”图像1410。在S230从每个造影图像1410中减去(元件1415)掩模图像1405，并且在S270对每个结果得到的差值图像应用滤波器1420。滤波器1420从控制系统1430接收滤波参数(例如，带通或带阻)。控制系统1430还提供权重α，在S280基于该权重α组合已滤波图像(元件1440)和差值图像(元件1450)。使组合图像经受处理器1460的额外处理，并且随后在S290由显示器1470来显不。
一些实施例针对彩色x射线图像进行操作。在这点上，在S220接收单个彩色x射线图像，其中造影剂被色彩编码以指示在不同时间点其存在到达最高点之处。过程200如上所述那样工作，以在S280生成单个图像，其是减去掩模的原始彩色图像和减去掩模的原始彩色图像的带通／带阻已滤波版本的加权组合。
本领域技术人员将认识到，可以配置对上述实施例的各种改编和修改，而不脱离权利要求的范围和精神。因而，应当理解，可以以与本文具体描述的不同方式来实施权利要求。
Claims (26)
1.一种系统，包括：
界面，其用于：
接收患者体积的掩模x射线图像；以及
接收包括造影剂的患者体积的多个相继x射线图像；以及
图像数据处理器，其用于：
从所述多个相继x射线图像中的每一个中减去所述掩模x射线图像，以生成多个相继x射线差值图像；
基于一个或多个滤波参数对所述多个相继x射线差值图像中的每一个进行滤波，以生成多个相继已滤波x射线差值图像；以及
基于第一权重将所述多个相继已滤波x射线差值图像中的每一个与所述多个相继x射线差值图像中的对应的一个组合，以生成第一多个组合相继x射线图像；以及
显示器，其用于相继地显示所述第一多个组合相继x射线图像。
2.根据权利要求1所述的系统，所述图像数据处理器还用于：
基于一个或多个第二滤波参数对所述多个相继x射线差值图像中的每一个进行滤波，以生成第二多个相继已滤波x射线差值图像；以及
基于所述第一权重将所述第二多个相继已滤波x射线差值图像中的每一个与所述多个相继x射线差值图像中的对应的一个组合，以生成第二多个组合相继x射线图像；以及
所述显示器用于相继地显示所述第二多个组合相继x射线图像。
3.根据权利要求2所述的系统，其中所述图像数据处理器在显示所述第一多个组合相继x射线图像期间基于所述一个或多个第二滤波参数对所述多个相继x射线差值图像中的每一个进行滤波并且基于所述第一权重将所述第二多个相继已滤波x射线差值图像中的每一个与所述多个相继x射线差值图像中的对应的一个组合。
4.根据权利要求3所述的系统，其中所述一个或多个滤波参数限定带阻滤波器，并且所述一个或多个第二滤波参数限定带通滤波器。
5.根据权利要求3所述的系统，其中所述一个或多个滤波参数限定第一带阻滤波器，并且所述一个或多个第二滤波参数限定与第一带阻滤波器不同的第二带阻滤波器。
6.根据权利要求3所述的系统，其中所述一个或多个滤波参数限定第一带通滤波器，并且所述一个或多个第二滤波参数限定与第一带通滤波器不同的第二带通滤波器。
7.根据权利要求2所述的系统，其中所述一个或多个滤波参数限定带阻滤波器，并且所述一个或多个第二滤波参数限定带通滤波器。
8.根据权利要求2所述的系统，所述图像数据处理器还用于：
基于第二权重将所述第二多个相继已滤波x射线差值图像中的每一个与所述多个相继x射线差值图像中的对应的一个组合，以生成第二多个组合相继x射线图像；以及
所述显示器用于相继地显示所述第二多个组合相继x射线图像。
9.根据权利要求8所述的系统，其中所述图像数据处理器在显示所述第一多个组合相继x射线图像期间基于所述第二权重将所述第二多个相继已滤波x射线差值图像中的每一个与所述多个相继x射线差值图像中的对应的一个组合。
10.根据权利要求1所述的系统，所述图像数据处理器还用于：
基于一个或多个第二滤波参数对所述多个相继x射线差值图像中的每一个进行滤波，以生成第二多个相继已滤波x射线差值图像；以及
基于第二权重将所述第二多个相继已滤波x射线差值图像中的每一个与所述多个相继x射线差值图像中的对应的一个组合，以生成第二多个组合相继x射线图像；以及
所述显示器用于相继地显示所述第二多个组合相继x射线图像。
11.根据权利要求10所述的系统，其中所述图像数据处理器在显示所述第一多个组合相继x射线图像期间基于所述一个或多个第二滤波参数对所述多个相继x射线差值图像中的每一个进行滤波并且基于所述第二权重将所述第二多个相继已滤波x射线差值图像中的每一个与所述多个相继x射线差值图像中的对应的一个组合。
12.根据权利要求10所述的系统，其中所述一个或多个滤波参数限定带阻滤波器，并且所述一个或多个第二滤波参数限定带通滤波器。
13.根据权利要求12所述的系统，其中所述图像数据处理器在显示所述第一多个组合相继x射线图像期间基于所述一个或多个第二滤波参数对所述多个相继x射线差值图像中的每一个进行滤波并且基于所述第二权重将所述第二多个相继已滤波x射线差值图像中的每一个与所述多个相继x射线差值图像中的对应的一个组合。
14.一种方法，包括：
接收患者体积的掩模x射线图像；
接收包括造影剂的患者体积的多个相继x射线图像；
从所述多个相继x射线图像中的每一个中减去所述掩模x射线图像，以生成多个相继x射线差值图像；
基于一个或多个滤波参数对所述多个相继x射线差值图像中的每一个进行滤波，以生成多个相继已滤波x射线差值图像；以及
基于第一权重将所述多个相继已滤波x射线差值图像中的每一个与所述多个相继x射线差值图像中的对应的一个组合，以生成第一多个组合相继x射线图像；以及
相继地显示所述第一多个组合相继x射线图像。
15.根据权利要求14所述的方法，还包括：
基于一个或多个第二滤波参数对所述多个相继x射线差值图像中的每一个进行滤波，以生成第二多个相继已滤波x射线差值图像；
基于所述第一权重将所述第二多个相继已滤波x射线差值图像中的每一个与所述多个相继x射线差值图像中的对应的一个组合，以生成第二多个组合相继x射线图像；以及
相继地显示所述第二多个组合相继x射线图像。
16.根据权利要求15所述的方法，其中在显示所述第一多个组合相继x射线图像期间，基于所述一个或多个第二滤波参数对所述多个相继x射线差值图像中的每一个进行滤波，并且基于所述第一权重将所述第二多个相继已滤波x射线差值图像中的每一个与所述多个相继x射线差值图像中的对应的一个组合。
17.根据权利要求16所述的方法，其中所述一个或多个滤波参数限定带阻滤波器，并且所述一个或多个第二滤波参数限定带通滤波器。
18.根据权利要求16所述的方法，其中所述一个或多个滤波参数限定第一带阻滤波器，并且所述一个或多个第二滤波参数限定与第一带阻滤波器不同的第二带阻滤波器。
19.根据权利要求16所述的方法，其中所述一个或多个滤波参数限定第一带通滤波器，并且所述一个或多个第二滤波参数限定与第一带通滤波器不同的第二带通滤波器。
20.根据权利要求15所述的方法，其中所述一个或多个滤波参数限定带阻滤波器，并且所述一个或多个第二滤波参数限定带通滤波器。
21.根据权利要求15所述的方法，还包括：
基于第二权重将所述第二多个相继已滤波x射线差值图像中的每一个与所述多个相继x射线差值图像中的对应的一个组合，以生成第二多个组合相继x射线图像；以及
相继地显示所述第二多个组合相继x射线图像。
22.根据权利要求21所述的方法，其中在显示所述第一多个组合相继x射线图像期间，基于所述第二权重将所述第二多个相继已滤波x射线差值图像与所述多个相继x射线差值图像中的对应的一个组合。
23.根据权利要求14所述的方法，还包括：
基于一个或多个第二滤波参数对所述多个相继x射线差值图像中的每一个进行滤波，以生成第二多个相继已滤波x射线差值图像；
基于第二权重将所述第二多个相继已滤波x射线差值图像中的每一个与所述多个相继x射线差值图像中的对应的一个组合，以生成第二多个组合相继x射线图像；以及
相继地显示所述第二多个组合相继x射线图像。
24.根据权利要求23所述的方法，其中在显示所述第一多个组合相继x射线图像期间，基于所述一个或多个第二滤波参数对所述多个相继x射线差值图像中的每一个进行滤波，并且基于所述第二权重将所述第二多个相继已滤波x射线差值图像中的每一个与所述多个相继x射线差值图像中的对应的一个组合。
25.根据权利要求23所述的方法，其中所述一个或多个滤波参数限定带阻滤波器，并且所述一个或多个第二滤波参数限定带通滤波器。
26.根据权利要求25所述的方法，其中在显示所述第一多个组合相继x射线图像期间，基于所述一个或多个第二滤波参数对所述多个相继x射线差值图像中的每一个进行滤波，并且基于所述第二权重将所述第二多个相继已滤波x射线差值图像中的每一个与所述多个相继x射线差值图像中的对应的一个组合。
Applications Claiming Priority (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US201261726035P | 2012-11-14 | 2012-11-14 | |
| US61/726035 | 2012-11-14 | ||
| US14/064470 | 2013-10-28 | ||
| US14/064,470 US9320486B2 (en) | 2012-11-14 | 2013-10-28 | System for viewing vasculature and perfuse tissue |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103976750A CN103976750A (zh) | 2014-08-13 |
| CN103976750B true CN103976750B (zh) | 2017-06-23 |
Family
ID=50681740
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310757218.7A Active CN103976750B (zh) | 2012-11-14 | 2013-11-14 | 用于观察脉管系统以及灌注组织的系统 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US9320486B2 (zh) |
| CN (1) | CN103976750B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR102153345B1 (ko) * | 2013-05-30 | 2020-09-09 | 삼성전자주식회사 | 방사선 촬영 장치 및 방사선 촬영 장치의 제어 방법 |
| CN104287764B (zh) * | 2014-09-11 | 2017-05-31 | 沈阳东软医疗系统有限公司 | 一种ct灌注成像方法和设备 |
| US10013779B2 (en) * | 2015-06-22 | 2018-07-03 | Toshiba Medical Systems Corporation | Metal artifact reduction for 3D-digtial subtraction angiography |
| US10147171B2 (en) * | 2016-09-21 | 2018-12-04 | General Electric Company | Systems and methods for generating subtracted images |
| US10297038B2 (en) | 2017-03-30 | 2019-05-21 | Siemens Healthcare Gmbh | Determination and presentation of flow transit curves |
| JP6945345B2 (ja) * | 2017-05-22 | 2021-10-06 | シャープ株式会社 | 表示装置、表示方法及びプログラム |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1521695A (zh) * | 2003-02-13 | 2004-08-18 | ��ʽ���綫֥ | 从图像中减少噪声的图像处理设备 |
| CN102665558A (zh) * | 2009-12-18 | 2012-09-12 | 佳能株式会社 | X射线图像处理设备、x射线图像处理方法和计算机程序的存储介质 |
Family Cites Families (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4504908A (en) * | 1982-03-15 | 1985-03-12 | General Electric Company | Matched filter for X-ray temporal subtraction |
| US4444196A (en) * | 1982-06-21 | 1984-04-24 | Diagnostic Technology, Inc. | Digital intravenous subtraction angiography |
| US4551800A (en) * | 1982-11-26 | 1985-11-05 | General Electric Company | Integrated hybrid image remasking in a subtraction angiography method |
| DE69724906T2 (de) * | 1996-10-16 | 2004-07-22 | Koninklijke Philips Electronics N.V. | Numerisches Bildverarbeitungsverfahren zur automatischen Extraktion bandförmiger Objekte |
| US7826884B2 (en) * | 2006-07-31 | 2010-11-02 | Siemens Medical Solutions Usa, Inc. | Live fluoroscopic roadmapping including targeted automatic pixel shift for misregistration correction |
| US7907217B2 (en) * | 2006-10-31 | 2011-03-15 | Siemens Medical Solutions Usa, Inc. | Systems and methods of subtraction angiography utilizing motion prediction |
| CN102172330B (zh) * | 2007-07-10 | 2013-03-27 | 株式会社东芝 | X射线摄影装置以及图像处理显示装置 |
| US8437519B2 (en) | 2007-10-19 | 2013-05-07 | Siemens Medical Solutions Usa, Inc. | Automated image data subtraction system suitable for use in angiography |
| US8050474B2 (en) * | 2008-04-11 | 2011-11-01 | Siemens Medical Solutions Usa, Inc. | System for generation of a composite medical image of vessel structure |
| US8150127B2 (en) | 2008-05-28 | 2012-04-03 | Siemens Medical Solutions Usa, Inc. | Method for automatically synchronizing the review of two DSA scenes |
| US8498463B2 (en) * | 2008-08-13 | 2013-07-30 | Koninklijke Philips N.V. | Mask construction for cardiac subtraction |
| US8355557B2 (en) * | 2008-10-10 | 2013-01-15 | Siemens Aktiengesellschaft | System and method for decomposed temporal filtering for X-ray guided intervention application |
| US8299413B2 (en) * | 2009-01-08 | 2012-10-30 | Siemens Aktiengesellschaft | Method for pixel shift calculation in digital subtraction angiography and X-ray diagnostic imaging system for generating images in digital subtraction angiography |
| 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 |
| US20130241552A1 (en) * | 2010-12-07 | 2013-09-19 | Hitachi Medical Corporation | Magnetic resonance imaging apparatus and contrast-enhanced image acquisition method |
| US8553963B2 (en) * | 2011-02-09 | 2013-10-08 | Siemens Medical Solutions Usa, Inc. | Digital subtraction angiography (DSA) motion compensated imaging system |
| DE102011005777B4 (de) * | 2011-03-18 | 2018-04-26 | Siemens Healthcare Gmbh | Angiographiesystem zur angiographischen Untersuchung eines Patienten und angiographisches Untersuchungsverfahren |
| US8463012B2 (en) | 2011-10-14 | 2013-06-11 | Siemens Medical Solutions Usa, Inc. | System for comparison of medical images |
| US8848996B2 (en) * | 2012-02-17 | 2014-09-30 | Siemens Medical Solutions Usa, Inc. | System for suppressing vascular structure in medical images |
-
2013
- 2013-10-28 US US14/064,470 patent/US9320486B2/en not_active Expired - Fee Related
- 2013-11-14 CN CN201310757218.7A patent/CN103976750B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1521695A (zh) * | 2003-02-13 | 2004-08-18 | ��ʽ���綫֥ | 从图像中减少噪声的图像处理设备 |
| CN102665558A (zh) * | 2009-12-18 | 2012-09-12 | 佳能株式会社 | X射线图像处理设备、x射线图像处理方法和计算机程序的存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| US9320486B2 (en) | 2016-04-26 |
| US20140133731A1 (en) | 2014-05-15 |
| CN103976750A (zh) | 2014-08-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103976750B (zh) | 用于观察脉管系统以及灌注组织的系统 | |
| US10470732B2 (en) | System and method for generating a time-encoded blood flow image from an arbitrary projection | |
| EP3545837B1 (en) | Medical image providing apparatus and medical image processing method of the same | |
| JP6104525B2 (ja) | 画像処理装置及びｘ線診断装置 | |
| CN102670225B (zh) | 自动成像造影剂确定系统 | |
| CN102988071B (zh) | 造影剂灌注自适应成像系统 | |
| US12070346B2 (en) | Medical image diagnosis apparatus and medical information display controlling device | |
| CN102542584B (zh) | 图像处理装置以及图像处理方法 | |
| JP6214646B2 (ja) | 血管造影図における時間的な解剖学的標的タグ付け | |
| US20080317195A1 (en) | Medical-diagnosis assisting apparatus, medical-diagnosis assisting method, and radiodiagnosis apparatus | |
| US10032293B2 (en) | Computed tomography (CT) apparatus and method of reconstructing CT image | |
| EP3061073A2 (en) | Image visualization | |
| CN103860185B (zh) | 多能量图像的确定 | |
| US10052078B2 (en) | Segmentation of moving structure in image data | |
| KR20200086919A (ko) | 단층 영상 촬영 장치 및 단층 영상 생성 방법 | |
| US11179122B2 (en) | Bolus imaging | |
| JP6479919B2 (ja) | 流動データの再構築 | |
| US20150164450A1 (en) | System and Method for Real Time 4D Quantification | |
| KR101590874B1 (ko) | 의료 영상 장치의 영상 재구성 방법 및 그 의료 영상 장치 | |
| US10297038B2 (en) | Determination and presentation of flow transit curves | |
| US10165997B2 (en) | System for acquiring a three-dimensional image of arteries and veins | |
| KR102273022B1 (ko) | 단층 촬영 장치 및 그에 따른 단층 영상 복원 방법 | |
| Badea et al. | A LabVIEW Platform for Preclinical Imaging Using Digital Subtraction Angiography and Micro‐CT | |
| JP6498516B2 (ja) | 医用画像処理装置及びｘ線ｃｔ装置 | |
| CA3066504A1 (en) | Systems and methods for medical imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20170210 Address after: Erlangen, Germany Applicant after: SIEMENS HEALTHCARE GmbH Address before: American Pennsylvania Applicant before: SIEMENS MEDICAL SOLUTIONS USA, Inc. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240902 Address after: German Phu F Haim Patentee after: Siemens Medical AG Country or region after: Germany Address before: Erlangen, Germany Patentee before: SIEMENS HEALTHCARE GmbH Country or region before: Germany |