CN105025794B - 谱ct的结构传播恢复 - Google Patents
谱ct的结构传播恢复 Download PDFInfo
- Publication number
- CN105025794B CN105025794B CN201480009706.2A CN201480009706A CN105025794B CN 105025794 B CN105025794 B CN 105025794B CN 201480009706 A CN201480009706 A CN 201480009706A CN 105025794 B CN105025794 B CN 105025794B
- Authority
- CN
- China
- Prior art keywords
- local
- data set
- projection data
- image data
- 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
- 230000003595 spectral effect Effects 0.000 title claims abstract description 54
- 238000011084 recovery Methods 0.000 title claims description 29
- 238000000034 method Methods 0.000 claims abstract description 40
- 238000000605 extraction Methods 0.000 claims description 6
- 238000003384 imaging method Methods 0.000 claims description 6
- 230000001902 propagating effect Effects 0.000 claims description 4
- 230000001629 suppression Effects 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 claims description 2
- 238000002591 computed tomography Methods 0.000 description 19
- 239000000463 material Substances 0.000 description 19
- 238000000701 chemical imaging Methods 0.000 description 8
- 230000005855 radiation Effects 0.000 description 8
- 230000000644 propagated effect Effects 0.000 description 7
- 238000000354 decomposition reaction Methods 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 4
- 230000002596 correlated effect Effects 0.000 description 3
- 230000000875 corresponding effect Effects 0.000 description 3
- 230000009977 dual effect Effects 0.000 description 3
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 2
- 230000004075 alteration Effects 0.000 description 2
- 230000002146 bilateral effect Effects 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 229910052740 iodine Inorganic materials 0.000 description 2
- 239000011630 iodine Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 239000002355 dual-layer Substances 0.000 description 1
- 238000000295 emission spectrum Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
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/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/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/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
-
- 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
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- 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
-
- 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
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/408—Dual energy
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- High Energy & Nuclear Physics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Pulmonology (AREA)
- Quality & Reliability (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
一种方法包括：获得来自谱扫描的投影数据或根据所述投影数据生成的图像数据中的至少一种，从所述投影数据或所述图像数据中的所述至少一种中选择局部参考数据集，确定针对所选择的参考数据集的噪声模式，基于所述噪声模式根据所述参考数据集来估计潜在的局部结构，并且基于所估计的潜在的局部结构来恢复所述投影数据或所述图像数据中的至少一种。
Description
技术领域
下文总体上涉及谱投影数据和/或谱图像数据处理，并且更具体地涉及谱投影数据和/或谱图像数据的结构传播恢复，并且结合对计算机断层摄影(CT)的具体应用而进行描述。然而，下文还适用于其他成像模态。
背景技术
谱(或多能量)CT已经利用同时以两种不同的光子能量(例如，利用双能量CT)采集的两个衰减值来求解包括被扫描的材料的质量衰减系数的光电和康普顿贡献，并且然后通过其光电和康普顿贡献的值来识别未知的材料。该方法非常适合具有接近诊断能量范围的平均值的k边缘能量的诸如碘的材料。因为两个基函数的任何两个线性无关的和跨越整个衰减系数空间，所以任何材料都能够由两种其他材料的线性组合来表示，所述两种其他材料一般称为基础材料，例如水和碘。
基础材料图像提供诸如单色图像、材料抵消图像、有效原子序数图像和电子密度图像的新应用。随着最近的技术进步，存在若干用于执行诸如双源、快速kVp切换和双层探测器配置的双能量CT采集的方法。另外，定量成像是当前医学成像界的主要趋势之一。谱CT支持该趋势，因为额外的谱信息改进能够测得的关于被扫描的目标及其材料成分的定量信息。
双能量材料分解是用于提取表示两种基础材料的射束硬化无关的线积分的数学方案，所述射束硬化无关的线积分即被扫描的目标内的光电吸收和康普顿散射。然而，谱CT遭受固有的噪声挑战，尤其是在进行了材料分解的情景中。因为材料分解是不适定问题，所以分解显著放大了噪声，其中，所获得的噪声在各个材料之间是高度反相关的。遗憾的是，如果直接根据材料线积分来重建图像，则图像噪声非常多。因此，所获得的图像质量很低，这可能降低图像的临床价值。
WO2011/064683是涉及增强传统CT扫描的图像质量以便减少由被扫描的对象接收到的X射线剂量的在先申请。
WO2013/011418是涉及谱CT中的局部噪声减少的在先申请。其描述了一种按体素的体素降噪方法。
发明内容
本文中描述的各方面解决以上提及的问题和其他问题。
下文描述了一种使用恢复算法来改进谱CT图像质量的方法，所述恢复算法从成像研究(即，投影数据和/或图像数据)移除噪声和/或伪影，同时保留潜在的目标结构和谱信息。所述算法能够被应用于投影域、图像域、或投影域和图像域两者中。一般地，参考数据集，即图像或正弦图，根据谱扫描的投影或图像数据来导出并且结合确定的噪声模式被使用，以估计潜在的局部目标结构，所述潜在的局部目标结构也被利用作为额外的约束。所述估计以及任选地所述额外的约束提供移除噪声和伪影的鲁棒且高质量的恢复。
在一个方面中，一种方法包括：获得来自谱扫描的投影数据或者根据所述投影数据生成的图像数据中的至少一种，从所述投影数据或所述图像数据中的所述至少一种中选择局部参考数据集，确定针对所选择的参考数据集的噪声模式，基于所述噪声模式根据所述参考数据集来估计潜在的局部结构，并且基于所估计的潜在的局部结构来恢复所述投影数据或所述图像数据中的至少一种。
在另一个方面中，一种投影数据和/或图像数据处理器包括：存储器，其存储结构传播算法；以及处理器，其结合来自谱扫描的投影数据或图像数据中的至少一种来运行所述结构传播算法以从所述投影数据或所述图像数据中的所述至少一种移除噪声或伪影中的至少一种，同时保留潜在的目标结构和谱信息。
在另一方面中，一种计算机可读存储介质被编码有计算机程序指令。所述计算机可读指令当由处理器运行时，令所述处理器：获得来自谱扫描的投影数据或根据所述投影数据生成的图像数据中的至少一种，从所述投影数据或所述图像数据中的所述至少一种中选择局部参考数据集，确定针对所选择的参考数据集的噪声模式，基于所述噪声模式根据所述参考数据集来估计潜在的局部结构，并且基于所估计的潜在的局部结构来恢复投影数据或图像数据中的至少一种。
附图说明
本发明可以采取各种部件和部件的布置以及各种步骤和步骤的安排的形式。附图仅出于图示优选的实施例的目的并且不应当被解释为对本发明的限制。
图1示意性地图示了结合采用结构传播恢复算法的投影数据和/或图像数据处理器的范例成像系统；
图2示意性地图示了结构传播恢复算法的范例模型；
图3图示了根据第一投影数据重建的现有技术图像；
图4图示了使用本文中描述的结构传播恢复方法结合第一投影数据的恢复的图像；
图5图示了根据第二投影数据重建的现有技术图像；
图6图示了使用本文中描述的结构传播恢复方法结合第二投影数据的恢复的图像；
图7图示了根据第三投影数据重建的现有技术图像；
图8图示了使用本文中描述的结构传播恢复方法结合第三投影数据的恢复的图像；
图9图示了用于处理投影数据和/或图像数据的范例方法。
具体实施方式
下文描述了一种使用恢复算法来改进谱CT图像质量的方法，所述恢复算法从谱图像研究的投影数据和/或图像数据移除噪声和/或伪影，同时保留潜在的目标结构和谱信息。
首先参考图1，图示了谱成像系统100，例如谱计算机断层摄影(CT)扫描器。所图示的谱成像系统100利用如下面更详细地讨论的kVp切换以产生谱投影数据。谱成像系统100包括大体固定机架102和旋转机架104，旋转机架104由固定机架102可旋转地支撑并且绕检查区域106关于z轴旋转。
诸如X射线管的辐射源108与旋转机架104一起旋转并且发射贯穿检查区域106的辐射。辐射源电压控制器110控制辐射源108的平均或峰值发射电压。在一个实例中，这包括在扫描的视图之间、在扫描的视图内、和/或其他情况下在两个或更多个发射电压(例如，80keV与140keV，100keV与120keV，等等)之间切换发射电压。
探测器阵列112对着相对于辐射源108在检查区域106的对面的角度弧。探测器阵列112探测贯穿检查区域106的辐射并且生成指示所述辐射的投影数据。在所述扫描是多能量扫描并且辐射源电压在用于扫描的至少两个发射电压之间切换的情况下，探测器阵列112生成针对辐射源电压中的每个的投影数据(在本文中也称为正弦图)。
再次地，所图示的谱成像系统100利用kVp切换。在变型中，谱成像系统100包括以两个不同的发射电压发射辐射以产生谱投影数据的至少两个辐射源108，并且/或者探测器阵列112包括产生谱投影数据的能量分辨探测器。在另一变型中，谱成像系统100包括用于产生谱投影数据的以上方法和/或其他方法的组合。
重建器114重建谱投影数据，生成指示对象的被扫描部分或被定位在检查区域106中的目标的体积图像数据。这包括重建在一个或多个发射电压处的谱图像数据和/或在整个发射谱上的常规(非谱)图像数据。卧榻或对象支撑物116将对象或目标支撑在检查区域106中。
操作者控制台118包括诸如监视器的人类可读输出设备和诸如键盘、鼠标等的输入设备。驻存在控制台118上的软件允许操作者经由图形用户接口(GUI)或以其他方式与谱成像系统100交互和/或操作谱成像系统100。这可以包括选择多能量谱成像协议，在所述多能量谱成像协议中，发射电压在两个或更多个电压之间切换。
投影数据和/或图像数据处理器120包括至少一个微处理器122，至少一个微处理器122运行存储在诸如物理存储器124或其他非暂态存储介质的计算机可读存储介质中的至少一个计算机可读指令。微处理器122还可以运行由载波、信号或其他暂态介质承载的一个或多个计算机可读指令。投影数据和/或图像数据处理器120能够是控制台118和/或其他计算系统的部分。
所述至少一个计算机可读指令包括(一个或多个)结构传播恢复指令126，(一个或多个)结构传播恢复指令126能够被应用到谱投影数据(或谱域)和/或重建的图像数据(或图像域)。在(一个或多个)结构传播恢复指令126被应用到两者的情况下，其首先被应用到谱投影数据，生成恢复的投影数据，恢复的投影数据然后被重建以生成图像数据，并且然后(一个或多个)结构传播恢复指令126被应用到重建的图像数据。
如下面更详细地描述的，运行(一个或多个)结构传播恢复指令126包括：从输入的投影数据和/或图像数据中选择参考投影和/或图像数据集；根据所选择的参考投影和/或图像数据集来确定噪声模型；基于所述噪声模型根据参考数据集来估计潜在的局部目标结构；并且将所估计的潜在的局部目标结构传播通过输入的投影数据和/或图像数据；恢复投影数据和/或图像数据。
任选的，所估计的潜在的局部目标结构能够被用作对恢复的约束。另外，移除的纹理和/或噪声中的一些能够被添加回。应当理解，(一个或多个)结构传播恢复指令126通过移除噪声和/或伪影来改进谱数据的质量，同时保留潜在的目标结构和谱信息。
所恢复的投影数据和/或所恢复的图像数据能够被进一步处理、经由显示监视器被显示、被制成胶片、被存储在数据储存库中(例如，影像存储与传输系统或PACS、电子病历或EMR、辐射信息系统或RIS、医院信息系统或HIS、等等)、和/或以其他方式被利用。
图2示意性地图示了(一个或多个)结构传播恢复指令126的各范例模块200。
参考结构选择器202从输入的投影数据和/或图像数据中选择局部参考数据集，所述输入的投影数据和/或图像数据能够由成像系统100(图1)和/或其他系统生成。针对参考数据集选择的标准包括能够根据其导出或估计潜在的局部目标结构的参考数据。
合适的参考图像数据的范例包括但不限于：特定单色图像(例如，70keV)，因为噪声在各个材料之间是反相关的并且存在其中对比噪声比(CNR)为最优的特定单色图像；非谱(或常规)CT图像，组合的谱/非谱图像，例如基于所有能量或全谱的非谱CT图像，其可以允许更准确地提取结构；仅基于部分谱的低能量图像和/或高能量图像；和/或其他参考图像。
合适的参考投影数据的范例包括但不限于：特定单色正弦图(例如，70keV)，因为噪声在各个材料之间是反相关的；以及非谱(或常规)CT正弦图，组合的谱/非谱正弦图，例如，基于所有能量或全谱的没有材料分解的非谱CT正弦图；仅基于部分谱的没有材料分解的低能量正弦图和/或高能量正弦图；和/或其他参考投影数据。
在一个实例中，利用全局方法来选择参考数据。对于该方法，能够针对全部接收到的投影数据和/或图像数据来选择单个参考数据集。所述选择能够基于以下标准中的一个或多个：最小总体变分(total variation)；最小熵；在图像/正弦图的局部标准偏差上的最小中值；局部噪声估计的最小平均值；基于两个ROI的输入的CNR图像，等等。在方程1中示出了单色图像中的总体变分选择标准的范例：
方程1
其中，e是单色能量并且
参考图像R中的i,j像素。
在另一实例中，利用局部方法来选择参考数据。对于该方法，针对每个像素，选择最优参考数据集片(或投影数据和/或图像数据的子集)。所述片能够在来自潜在的参考数据集中的一个中的像素的位置附近选择。该方法利用以下事实：输入的投影数据和/或图像数据中的不同区域可能具有针对局部结构偏差的不同的最优参考数据集。
噪声建模器204对所选择的参数数据集的噪声模型进行建模。合适的建模方法的范例包括但不限于：蒙特卡罗估计、解析法、直接提取和/或其他方法。所获得的噪声模式或模型被利用于估计局部结构参考数据集。范例解析方法被描述在Wunderlich的“ImageCovariance and Lesion Detectability in Direct Fan-Beam X-Ray ComputedTomography”(Phys.Med.Biol.53(2008)，第2472至2493页)中。
直接提取方法的范例被描述在于2009年11月25日提交的并且题为“ENHANCEDIMAGE DATA/DOSE REDUCTION”的申请序列号为61/264340的专利申请中，通过引用将其整体并入本文中；并且被描述在于2009年12月15日提交的并且题为“ENHANCED IMAGE DATA/DOSE REDUCTION”的申请序列号为61/286477的专利申请中，通过引用将其整体并入本文中。在本文中也预见到其他方法。
恢复器206恢复投影数据和/或图像数据。所图示的恢复器206包括结构导出器208和结构传播器210。
结构导出器208根据参考数据集来估计潜在的局部结构。所述估计可以改进参考数据集的局部CNR，这有助于估计数据集并且使得能够进行非常准确的结构估计。结构导出器208利用由噪声建模器204生成的噪声模型作为引导，这有助于在噪声与潜在的目标结构之间进行区分。
对于估计，结构导出器208可以利用各种算法，例如双边滤波、扩散滤波、全变差去噪、均值漂移等。通过非限制性举例的方式，下文描述了一种使用双边算法来估计局部结构的方法。该范例包括任选的尖峰噪声抑制。
对于参考数据集R中的每个体素Ri,j,k，结构导出器208执行以下：
1、提取体素Ri,j,k周围的n个体素的子体积；
2、基于方程2来计算核
方程2：
其中，α是控制权重的竞争力的参数并且
是在以上讨论的噪声建模中估计的对Ri,j,k的局部噪声水平估计；
3、将局部核乘以空间核：
其中，能够基于方程3来确定具有标准偏差的三维(3D)空间高斯核；
方程3：
其中，dx是像素的尺寸(例如，以毫米或mm为单位)，dz是切片宽度(例如，以mm为单位)，并且σ空间是控制权重的竞争力的参数；
4、将wi',j',k'归一化以使和等于一；
5、应用尖峰噪声抑制如下：如果中心权重wi',j',k'>w阈值并且α<α最大，则α＝α*αult并且返回步骤2；并且
6、基于方程4来估计目标结构：
方程4：
结构传播器210将所估计的结构传播到输入的投影数据和/或图像数据，任选地，所述估计被利用作为对恢复的额外的约束。移除的纹理和/或噪声中的一些在恢复期间能够被添加回。这可以有助于控制最终图像外观。
为此，结构传播器210根据结构导出器208的输出来导出局部结构，并且根据噪声建模器204的输出来定义结构水平。该方法的强处在于其利用具有改进的CNR的非常良好地定义的局部结构。由于处理是在投影数据和/或图像数据上进行的，所以投影数据和/或图像数据的强度值非常准确地被保留。
通过非限制性举例的方式，在一个实例中，结构传播器210针对投影数据和/或图像数据中的每个体素Vi,j,k执行以下的N次迭代：
1、提取体素Vi,j,k周围的n个体素的子体积；
2、基于方程4来计算局部核
方程4：
其中，β是输入参数；
3、将局部核乘以空间核
4、将wi′,j′,k′归一化以使和等于一；
5、处理投影数据和/或图像数据，如方程5中所示：
方程5：
其中，
是输入的目标数据集并且Iter是当前迭代的索引。
在步骤5中，图像纹理和/或噪声中的一些能够被添加回所述目标数据集以控制其最终外观。存在针对纹理和/或噪声的两个任选的源：参考投影数据和/或图像数据或者目标投影数据和/或图像数据。
获得最终的恢复，如方程6或方程7中所示：
方程6:
或者
方程7：
其中，δ是输入参数。
图3-8结合在没有恢复的情况下生成的图像示出了恢复的范例结果。
对于第一组投影数据和/或图像数据，图3表示现有技术图像300，其中，未利用(一个或多个)结构传播恢复指令126，并且图4表示其中利用(一个或多个)结构传播恢复指令126的图像400。如所示出的，区域402比对应的区域302噪声更少，而没有图像400中的结构的可见损失。
对于第二组投影数据和/或图像数据，图5表示现有技术图像500，其中，未利用(一个或多个)结构传播恢复指令126，并且图6表示其中利用(一个或多个)结构传播恢复指令126的图像600。如所示出的，区域602比对应的区域502噪声更少，而没有图像600中的结构的可见损失。
对于第三组投影数据和/或图像数据，图7表示现有技术图像700，其中，未利用(一个或多个)结构传播恢复指令126，并且图8表示利用(一个或多个)结构传播恢复指令126的图像800。如所示出的，区域802比对应的区域702噪声更少，而没有图像800中的结构的可见损失。
图9图示了根据本文中的公开内容的范例方法。
应当理解，动作的顺序不是限制性的。因此，本文中预见到其他顺序。另外，可以省略一个或多个动作和/或可以包括一个或多个额外的动作。
在902，获得输入的投影数据和/或图像数据。
在904，从所述投影数据和/或所述图像数据中选择局部参考数据集。
在906，确定针对所选择的参考数据集的噪声模型。
在908，基于所述噪声模型来估计潜在的局部结构。
在910，将所估计的结构传播到所述投影数据和/或所述图像数据，恢复所述投影数据和/或所述图像数据。
在912，所恢复的投影数据和/或所恢复的图像数据能够被进一步处理、经由显示监视器被显示、被制成胶片、被归档在数据储存库中、和/或以其他方式被利用。
以上可以通过被编码或被嵌入在计算机可读存储介质中的计算机可读指令的方式来实施，所述计算机可读指令当由(一个或多个)计算机处理器运行时，令所述(一个或多个)处理器运行所描述的动作。额外地或备选地，所述计算机可读指令中的至少一个由信号、载波或其他暂态介质承载。
己经参考优选实施例描述了本发明。他人在阅读并理解前述详细的描述之后可以进行修改和更改。旨在将本发明解释为包括所有这样的修改和更改，只要它们落入权利要求书或其等价要件的范围内。
Claims (15)
1.一种成像方法，包括：
获得来自谱扫描的投影数据或根据来自谱扫描的投影数据生成的图像数据中的至少一种；
从所述投影数据或所述图像数据中的所述至少一种中选择处于至少一个能量或能量范围的局部参考数据集；
确定针对所选择的局部参考数据集的噪声模式；
通过以下操作基于所述噪声模式根据所述局部参考数据集来估计所述局部参考数据集的潜在的局部结构：提取在所述局部参考数据集中的体素周围的体素的第一子体积，基于所提取的第一子体积计算第一局部核，将所述第一局部核乘以第一空间核来产生第一权重，将所述第一权重归一化以使和等于一，并且基于所提取的第一子体积和归一化的第一权重来估计所述潜在的局部结构；并且
通过以下操作恢复所述投影数据或所述图像数据中的所述至少一种：提取在所述投影数据或所述图像数据中的所述至少一种中的体素周围的体素的第二子体积，基于所估计的潜在的局部结构计算第二局部核，将所述第二局部核乘以第二空间核来产生第二权重，将所述第二权重归一化以使和等于一，并且基于所提取的第二子体积和归一化的第二权重来将所述潜在的局部结构传播到所述投影数据或所述图像数据中的所述至少一种。
2.如权利要求1所述的方法，还包括：
利用所估计的潜在的局部结构作为对所述恢复的约束。
3.如权利要求1-2中的任一项所述的方法，还包括：
将移除的纹理或噪声中的至少一种的子部分添加回所恢复的投影数据或图像数据。
4.如权利要求1-2中的任一项所述的方法，其中，所述局部参考数据集是仅从所述投影数据中选择的。
5.如权利要求1-2中的任一项所述的方法，其中，所述局部参考数据集是仅从所述图像数据中选择的。
6.如权利要求1-2中的任一项所述的方法，其中，所述局部参考数据集是从所述投影数据中选择的并且所述投影数据被恢复，并且所述方法还包括：
从所恢复的投影数据中选择第二局部参考数据集；
确定针对所选择的第二局部参考数据集的第二噪声模式；
基于所述第二噪声模式根据所述第二局部参考数据集来估计第二潜在的局部结构；并且
基于所估计的第二潜在的局部结构来恢复所述图像数据。
7.如权利要求1-2中的任一项所述的方法，其中，所述局部参考数据集包括以下中的一项：特定单色图像、非谱图像、组合的谱图像/组合的非谱图像、低能量图像或高能量图像。
8.如权利要求1-2中的任一项所述的方法，其中，所述局部参考数据集包括以下中的一项：特定单色正弦图、非谱正弦图、组合的正弦图/组合的非谱正弦图、低能量正弦图或高能量正弦图。
9.如权利要求1-2中的任一项所述的方法，还包括：
基于全局方法来选择所述局部参考数据集，在所述全局方法中，针对所述投影数据或所述图像数据中的至少一种的整体来选择单个局部参考数据集。
10.如权利要求9所述的方法，其中，对单个局部参考数据集的所述选择基于以下中的一个或多个被执行：最小总体变分；最小熵；在局部标准偏差上的最小中值；局部噪声估计的最小平均值；以及基于两个感兴趣区域的输入的对比噪声比图像。
11.如权利要求1-2中的任一项所述的方法，还包括：
基于局部方法来选择所述局部参考数据集，在所述局部方法中，针对每个像素，选择最优局部参考数据集片。
12.如权利要求1-2中的任一项所述的方法，其中，所述噪声模式基于以下中的一个或多个被确定：蒙特卡罗估计、解析法、或直接提取。
13.如权利要求1-2中的任一项所述的方法，还包括：
通过以下方式基于所述噪声模式根据所述局部参考数据集来估计所述潜在的局部结构：提取一体素周围的体素的子体积，计算局部核，将所述局部核乘以空间核，将乘积归一化为等于一的和，应用尖峰噪声抑制，并且估计所述结构。
14.一种投影数据和/或图像数据处理器（120），包括：
存储器（124），其被配置为存储结构传播算法；以及
微处理器（122），其结合来自谱扫描的投影数据或图像数据中的至少一种来运行所述结构传播算法以从所述投影数据或所述图像数据中的所述至少一种移除噪声或伪影中的至少一种，同时保留潜在的目标结构和谱信息，
其中，所述微处理器：
从所述投影数据或所述图像数据中的所述至少一种中选择处于至少一个能量或能量范围的局部参考数据集；
确定针对所选择的局部参考数据集的噪声模式；
通过以下操作基于所述噪声模式根据所述局部参考数据集来估计所述局部参考数据集的潜在的局部结构：提取在所述局部参考数据集中的体素周围的体素的第一子体积，基于所提取的第一子体积计算第一局部核，将所述第一局部核乘以第一空间核来产生第一权重，将所述第一权重归一化以使和等于一，并且基于所提取的第一子体积和归一化的第一权重来估计所述潜在的局部结构；并且
通过以下操作恢复投影数据或图像数据中的所述至少一种：提取在所述投影数据或所述图像数据中的所述至少一种中的体素周围的体素的第二子体积，基于所估计的潜在的局部结构计算第二局部核，将所述第二局部核乘以第二空间核来产生第二权重，将所述第二权重归一化以使和等于一，并且基于所提取的第二子体积和归一化的第二权重来将所述潜在的局部结构传播到所述投影数据或所述图像数据中的所述至少一种。
15.一种编码有计算机可读指令的计算机可读存储介质，所述计算机可读指令当由处理器运行时，被配置为令所述处理器：
获得来自谱扫描的投影数据或根据来自谱扫描的投影数据生成的图像数据中的至少一种；
从所述投影数据或所述图像数据中的所述至少一种中来选择处于至少一个能量或能量范围的局部参考数据集；
确定针对所选择的局部参考数据集的噪声模式；
通过以下操作基于所述噪声模式根据所述局部参考数据集来估计所述局部参考数据集的潜在的局部结构：提取在所述局部参考数据集中的体素周围的体素的第一子体积，基于所提取的第一子体积计算第一局部核，将所述第一局部核乘以第一空间核来产生第一权重，将所述第一权重归一化以使和等于一，并且基于所提取的第一子体积和归一化的第一权重来估计所述潜在的局部结构；并且
通过以下操作恢复投影数据或图像数据中的所述至少一种：提取在所述投影数据或所述图像数据中的所述至少一种中的体素周围的体素的第二子体积，基于所估计的潜在的局部结构计算第二局部核，将所述第二局部核乘以第二空间核来产生第二权重，将所述第二权重归一化以使和等于一，并且基于所提取的第二子体积和归一化的第二权重来将所述潜在的局部结构传播到所述投影数据或所述图像数据中的所述至少一种。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US201361767300P | 2013-02-21 | 2013-02-21 | |
| US61/767,300 | 2013-02-21 | ||
| PCT/IB2014/058895 WO2014128595A1 (en) | 2013-02-21 | 2014-02-11 | Structure propagation restoration for spectral ct |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105025794A CN105025794A (zh) | 2015-11-04 |
| CN105025794B true CN105025794B (zh) | 2021-11-02 |
Family
ID=50159317
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201480009706.2A Active CN105025794B (zh) | 2013-02-21 | 2014-02-11 | 谱ct的结构传播恢复 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US10282820B2 (zh) |
| EP (1) | EP2958494B1 (zh) |
| JP (1) | JP6472088B2 (zh) |
| CN (1) | CN105025794B (zh) |
| WO (1) | WO2014128595A1 (zh) |
Families Citing this family (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3077990A1 (en) | 2013-12-06 | 2016-10-12 | Koninklijke Philips N.V. | Bone segmentation from image data |
| WO2016087394A1 (en) | 2014-12-05 | 2016-06-09 | Koninklijke Philips N.V. | Eyi chusystem for generating spectral computed tomography projection data |
| US10657679B2 (en) | 2015-03-09 | 2020-05-19 | Koninklijke Philips N.V. | Multi-energy (spectral) image data processing |
| DE102015206630B4 (de) * | 2015-04-14 | 2022-05-05 | Siemens Healthcare Gmbh | Multispektrale CT-Bildgebung |
| KR101697501B1 (ko) * | 2015-07-23 | 2017-01-18 | 서울대학교산학협력단 | Ct 이미지의 잡음 저감 방법 및 장치 |
| CN108475419B (zh) * | 2015-12-15 | 2022-04-26 | 皇家飞利浦有限公司 | 用于计算机断层摄影的数据处理的方法 |
| EP3270788B1 (en) * | 2016-05-20 | 2021-09-29 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for computed tomography |
| EP3661334B1 (en) | 2016-05-31 | 2021-07-07 | Koninklijke Philips N.V. | Apparatus for generating x-rays |
| EP3475916B1 (en) | 2016-11-30 | 2019-10-23 | Koninklijke Philips N.V. | Bone and hard plaque segmentation in spectral ct |
| JP7278224B2 (ja) * | 2017-06-02 | 2023-05-19 | コーニンクレッカ フィリップス エヌ ヴェ | 医用画像の病変の定量化される態様 |
| EP3467771A1 (en) * | 2017-10-05 | 2019-04-10 | Koninklijke Philips N.V. | Image feature annotation in diagnostic imaging |
| US11417034B2 (en) | 2018-01-16 | 2022-08-16 | Koninklijke Philips N.V. | Spectral imaging with a non-spectral imaging system |
| EP3616620A1 (en) * | 2018-08-28 | 2020-03-04 | Koninklijke Philips N.V. | Spectral dual-layer ct-guided interventions |
| EP3893205A1 (en) | 2020-04-08 | 2021-10-13 | Koninklijke Philips N.V. | Suppression of motion artifacts in computed tomography imaging |
| CN113706419B (zh) * | 2021-09-13 | 2024-07-19 | 上海联影医疗科技股份有限公司 | 一种图像处理方法和系统 |
| CN114972088A (zh) * | 2022-05-06 | 2022-08-30 | 北京师范大学 | 一种基于最小熵的固定模式噪声去除方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050111751A1 (en) * | 2003-11-26 | 2005-05-26 | Avinash Gopal B. | Method for random point and patterned noise reduction in digital images |
| US20090052612A1 (en) * | 2007-08-22 | 2009-02-26 | Xiaoye Wu | System and method of optimizing a monochromatic representation of basis material decomposed ct images |
| CN102655810A (zh) * | 2009-12-15 | 2012-09-05 | 皇家飞利浦电子股份有限公司 | 计算机断层摄影装置 |
| CN102667852A (zh) * | 2009-11-25 | 2012-09-12 | 皇家飞利浦电子股份有限公司 | 增强图像数据/剂量减小 |
| WO2013011418A2 (en) * | 2011-07-15 | 2013-01-24 | Koninklijke Philips Electronics N.V. | Spectral ct |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE69616031T2 (de) | 1995-12-21 | 2002-06-20 | Koninklijke Philips Electronics N.V., Eindhoven | Rauschreduzierung in einem bild |
| US7254209B2 (en) | 2003-11-17 | 2007-08-07 | General Electric Company | Iterative CT reconstruction method using multi-modal edge information |
| US20100061656A1 (en) | 2006-11-08 | 2010-03-11 | Koninklijke Philips Electronics N. V. | Noise reduction of an image signal |
| DE102008063311A1 (de) | 2008-12-30 | 2010-07-08 | Siemens Aktiengesellschaft | Verfahren zur Rauschreduktion von CT-Bilddaten und Bildbearbeitungssystem |
| DE102009010501A1 (de) | 2009-02-25 | 2010-09-09 | Siemens Aktiengesellschaft | Verfahren zur Rauschreduktion von CT-Bilddaten und Bildbearbeitungssystem |
-
2014
- 2014-02-11 WO PCT/IB2014/058895 patent/WO2014128595A1/en active Application Filing
- 2014-02-11 CN CN201480009706.2A patent/CN105025794B/zh active Active
- 2014-02-11 US US14/766,999 patent/US10282820B2/en active Active
- 2014-02-11 EP EP14706342.4A patent/EP2958494B1/en active Active
- 2014-02-11 JP JP2015558575A patent/JP6472088B2/ja active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050111751A1 (en) * | 2003-11-26 | 2005-05-26 | Avinash Gopal B. | Method for random point and patterned noise reduction in digital images |
| US20090052612A1 (en) * | 2007-08-22 | 2009-02-26 | Xiaoye Wu | System and method of optimizing a monochromatic representation of basis material decomposed ct images |
| US7724865B2 (en) * | 2007-08-22 | 2010-05-25 | General Electric Company | System and method of optimizing a monochromatic representation of basis material decomposed CT images |
| CN102667852A (zh) * | 2009-11-25 | 2012-09-12 | 皇家飞利浦电子股份有限公司 | 增强图像数据/剂量减小 |
| CN102655810A (zh) * | 2009-12-15 | 2012-09-05 | 皇家飞利浦电子股份有限公司 | 计算机断层摄影装置 |
| WO2013011418A2 (en) * | 2011-07-15 | 2013-01-24 | Koninklijke Philips Electronics N.V. | Spectral ct |
Also Published As
| Publication number | Publication date |
|---|---|
| EP2958494B1 (en) | 2022-12-21 |
| CN105025794A (zh) | 2015-11-04 |
| JP2016507320A (ja) | 2016-03-10 |
| US20150379694A1 (en) | 2015-12-31 |
| WO2014128595A1 (en) | 2014-08-28 |
| EP2958494A1 (en) | 2015-12-30 |
| JP6472088B2 (ja) | 2019-02-20 |
| US10282820B2 (en) | 2019-05-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105025794B (zh) | 谱ct的结构传播恢复 | |
| US10147168B2 (en) | Spectral CT | |
| US9959640B2 (en) | Iterative image reconstruction with a sharpness driven regularization parameter | |
| EP2923332B1 (en) | Projection data de-noising | |
| CN105793894B (zh) | 根据图像数据来进行骨骼分割 | |
| EP3134867B1 (en) | Restoration of low contrast structure in de-noise image data | |
| CN107430779B (zh) | 多能量（谱）图像数据处理 | |
| CN110574073B (zh) | 能谱计算机断层摄影(ct)成像中的残余碘伪影的探测和/或校正 | |
| Chen et al. | Line integral alternating minimization algorithm for dual-energy x-ray CT image reconstruction | |
| Gomi et al. | Development of a novel algorithm for metal artifact reduction in digital tomosynthesis using projection-based dual-energy material decomposition for arthroplasty: A phantom study | |
| WO2016097981A1 (en) | Penalized maximum likelihood material decomposition | |
| EP2864964B1 (en) | Image reconstruction in interleaved multi-energy imaging | |
| Zhu et al. | Feasibility of PET-enabled dual-energy CT imaging: First physical phantom and patient results | |
| Zhang et al. | Perfusion Computed Tomography Image Reconstruction Using Spatio-Temporal Constrained Edge-preserving Prior |
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 |