CN110136218B - 基于噪声生成机制与数据驱动紧框架的ct投影去噪重建方法及装置 - Google Patents
基于噪声生成机制与数据驱动紧框架的ct投影去噪重建方法及装置 Download PDFInfo
- Publication number
- CN110136218B CN110136218B CN201910243267.6A CN201910243267A CN110136218B CN 110136218 B CN110136218 B CN 110136218B CN 201910243267 A CN201910243267 A CN 201910243267A CN 110136218 B CN110136218 B CN 110136218B
- Authority
- CN
- China
- Prior art keywords
- data
- noise
- projection
- projection data
- sinogram
- 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
- 238000000034 method Methods 0.000 title claims abstract description 66
- 230000007246 mechanism Effects 0.000 title claims abstract description 42
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
- 238000005065 mining Methods 0.000 claims abstract description 12
- 238000005457 optimization Methods 0.000 claims description 27
- 230000006870 function Effects 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 16
- 238000012545 processing Methods 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 7
- 238000000354 decomposition reaction Methods 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 7
- 230000003190 augmentative effect Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 125000004122 cyclic group Chemical group 0.000 claims description 2
- 238000001914 filtration Methods 0.000 claims description 2
- 238000002591 computed tomography Methods 0.000 description 36
- 238000010586 diagram Methods 0.000 description 7
- 230000000694 effects Effects 0.000 description 4
- 230000005855 radiation Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 201000010099 disease Diseases 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 230000014759 maintenance of location Effects 0.000 description 2
- 230000005012 migration Effects 0.000 description 2
- 238000013508 migration Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 238000012885 constant function Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002599 functional magnetic resonance imaging Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 208000032839 leukemia Diseases 0.000 description 1
- 238000011551 log transformation method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000013526 transfer learning Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- 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
- 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
- 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]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明属于CT图像重建技术领域，特别涉及一种基于噪声生成机制与数据驱动紧框架的CT投影去噪重建方法及装置，该方法中：首先通过刻画投影数据噪声统计特性，并挖掘投影数据的稀疏先验信息；结合两者数据建立最大后验概率模型；然后对模型迭代求解，获取去噪后的正弦图数据；通过重建获取CT图像。本发明克服传统低剂量CT重建算法因采用近似的方式刻画投影数据的噪声导致重建图像易残留严重的噪声和伪影的问题，能够在尽可能保留投影数据有用信息的同时去除噪声，进而有效抑制了重建图像的噪声和伪影，在提升低剂量CT重建图像质量方面具有较高的应用价值。
Description
技术领域
本发明属于CT图像重建技术领域，特别涉及一种基于噪声生成机制与数据驱动紧框架 的CT投影去噪重建方法及装置。
背景技术
计算机断层成像(Computed Tomography,CT)利用X射线穿透物体的透射信息，通过重 建算法计算物体内部衰减特性分布。作为20世纪后期最伟大的科技成果之一，CT技术可以通 过非接触、非破坏的方式高效率地获取物体内部结构信息。1972年，第一台CT设备由工程师 Hounsfield成功研发，随后在医学领域获取了快速发展，掀起了临床疾病诊断的革命。CT可 以提供病人清晰准确的结构图像，有利于医师高效准确地辨别病变，目前已成为临床疾病诊 断不可或缺的工具之一。
随着CT技术在临床上的广泛应用，过多辐射剂量所造成的危害逐渐成为人们关注的焦 点。当辐射剂量超过安全范围时，会影响人体的身体健康，增加癌症、白血病等疾病的患病 风险。1990年，Naidich等人首次提出低剂量CT(Low-Dose Computed Tomography,LDCT)的 概念，引起了CT领域研究者的广泛关注。目前，LDCT的实现主要有两种方式，一是通过降 低X射线源的管电流来降低单角度投影采集时的X射线强度；二是减少投影数量，即稀疏角度 问题。由于实际临床扫描速度较快，稀疏采样不易实现，而降低管电流可操作性强，简便有 效，目前已成为临床上实现LDCT的主要手段。但降低管电流时，由于探测器接收到的光子数 大量减少，采集得到的投影往往被光子统计噪声和探测系统热噪声严重污染，是CT图像质量 严重下降的本质原因。因此，亟需一种保证去噪以及细节保留能力的低剂量CT投影域重建方 法，以解决CT临床中过多辐射剂量的危害。
发明内容
为此，本发明提供一种基于噪声生成机制与数据驱动紧框架的CT投影去噪重建方法及 装置，克服低剂量CT中噪声的严重影响，提高CT图像重建质量。
按照本发明所提供的设计方案，一种基于噪声生成机制与数据驱动紧框架的CT投影去 噪重建方法，包含如下内容：
刻画投影数据噪声统计特性，并挖掘投影数据的稀疏先验信息；
结合噪声统计特性和稀疏先验信息，建立最大后验概率模型；
对最大后验概率模型进行迭代求解，获取去噪后的正弦图数据；
对去噪后的正弦图数据进行重建，得到重建后的CT图像。
上述的，刻画投影数据噪声统计特性中，通过分析投影数据噪声本质来源，考虑光子统 计噪声和系统热噪声统计特性，建立噪声生成机制，通过该噪声生成机制对投影数据的统计 特性进行刻画。
优选的，挖掘投影数据的稀疏先验信息中，将稀疏先验信息嵌入噪声生成机制，并通过 数据驱动紧框架挖掘投影数据的稀疏先验信息。
上述的，建立最大后验概率模型，包含如下内容：首先，采用基于小波变换的数据驱动 紧框架对投影数据进行稀疏编码，获取投影数据和先验信息的联合分布；然后，结合噪声生 成机制，获取后验概率分布；基于该后延概率分布，采用MAP框架估计正弦图，并得到最 大后验概率模型。
上述的，对最大后验概率模型进行迭代求解中，将模型求解等价为带约束的优化问题； 并利用增广拉格朗日函数将带约束的优化问题转化为无约束的优化问题；采用交替方向承子 ADMM法将该无约束的优化问题进行问题分解；根据预先设定的循环迭代条件，对分解后的 问题进行迭代求解，得到去噪后的正弦图数据。
优选的，用交替方向承子ADMM法将该无约束的优化问题进行问题分解中，依据接收 射线量子、正弦图及尺度参数进行问题分解，获取若干子问题；然后，根据预先设定的循环 迭代条件，对分解后的该若干子问题进行迭代求解，得到去噪后的正弦图数据。
优选的，对分解后的若干子问题迭代求解过程，包含如下内容：初始化求解参数，获取 投影数据和正弦图；依次对子问题、求解参数和小波基函数进行参数更新求解，并判断是否 达到预设最大迭代次数，若是，则判定当前更新求解过程的正弦图数据即为去噪后的正弦图 数据，否则，返回参数更新求解步骤继续执行。
上述的，采用经典解析FBP算法对去噪后的正弦图数据进行重建，得到重建后的CT图 像。
一种基于噪声生成机制与数据驱动紧框架的CT投影去噪重建装置，包含：数据处理模 块、模型构建模块、模型求解模块和图像重建模块，其中，
数据处理模块，用于刻画投影数据噪声统计特性，并挖掘投影数据的稀疏先验信息；
模型构建模块，用于结合噪声统计特性和稀疏先验信息，建立最大后验概率模型；
模型求解模块，用于对最大后验概率模型进行迭代求解，获取去噪后的正弦图数据；
图像重建模块，用于对去噪后的正弦图数据进行重建，得到重建后的CT图像。
上述的装置中，模型构建模块包含：稀疏编码单元、概率分布计算单元和模型获取单元， 其中，
稀疏编码单元，用于采用基于小波变换的数据驱动紧框架对投影数据进行稀疏编码，获 取投影数据和先验信息的联合分布；
概率分布计算单元，用于结合噪声生成机制，获取后验概率分布；
模型获取单元，用于基于后延概率分布，采用MAP框架估计正弦图，并得到最大后验 概率模型。
本发明的有益效果：
本发明中，首先通过刻画投影数据噪声统计特性，并挖掘投影数据的稀疏先验信息；结 合两者数据建立最大后验概率模型；然后对模型迭代求解，获取去噪后的正弦图数据；通过 重建获取CT图像；克服了传统低剂量CT重建算法因采用近似的方式刻画投影数据的噪声导 致重建图像易残留严重的噪声和伪影的问题；并进一步引入噪声生成机制来充分刻画投影数 据固有的“泊松+高斯”统计特性，同时利用基于学习策略的数据驱动紧框架对投影数据进行稀 疏信息挖掘，利用ADMM算法对模型进行高效求解，能够在尽可能保留投影数据有用信息 的同时去除噪声，进而有效抑制了重建图像的噪声和伪影，在提升低剂量CT重建图像质量 方面具有较高的应用价值。
附图说明：
图1为实施例中方法流程示意图；
图2为实施例中计算过程示意图；
图3为实施例中装置示意图；
图4为实施例中模型构建模块示意图；
图5为实施例中仿真实验重建图像对比图。
具体实施方式：
为使本发明的目的、技术方案和优点更加清楚、明白，下面结合附图和技术方案对本发 明作进一步详细的说明。
基于投影预处理的LDCT重建算法的基本思路是首先对LDCT系统直接采集的投影数据 进行预处理以抑制投影域噪声，而后采用FBP算法对去噪后的投影数据进行重建，从而改善 LDCT图像的质量。现阶段国内外学者对该类算法的研究主要分为两种方式：一种是通过改 进滤波器的性能以实现去噪效果的改善；另一种是基于投影数据的统计特性建立统计迭代模 型对投影数据进行噪声抑制。该种方法将投影数据去噪过程看作优化函数问题，在投影数据 统计特性的基础上，建立统计迭代模型进行去噪，由于这种类方法更具有针对性，因此受到 了国内外学者的广泛关注并在实际应用中取得了较好的效果。通过对LDCT多次扫描所得投 影数据的统计特性进行探究，发现对采集的投影数据进行校正和对数变换后，其噪声服从一 种非平稳高斯分布，且方差与均值满足非线性关系。为此，本发明实施例提出基于迁移学习 的fMRI视觉编码模型。所谓迁移学习，通俗来说就是将从源领域学习到的东西应用到新的 领域。迁移学习可以解决小数据问题，能够将用于大数据的模型迁移到小数据上。本发明实 施例中，参见图1所示，提供一种基于噪声生成机制与数据驱动紧框架的CT投影去噪重建 方法，包含如下内容：
S101、刻画投影数据噪声统计特性，并挖掘投影数据的稀疏先验信息；
S102、结合噪声统计特性和稀疏先验信息，建立最大后验概率模型；
S103、对最大后验概率模型进行迭代求解，获取去噪后的正弦图数据；
S104、对去噪后的正弦图数据进行重建，得到重建后的CT图像。
首先通过刻画投影数据噪声统计特性，并挖掘投影数据的稀疏先验信息；结合两者数据 建立最大后验概率模型；然后对模型迭代求解，获取去噪后的正弦图数据；通过重建获取CT 图像；克服了传统低剂量CT重建算法因采用近似的方式刻画投影数据的噪声导致重建图像 易残留严重的噪声和伪影的问题，提升图像重建质量。
本发明实施例中，刻画投影数据噪声统计特性中，通过分析投影数据噪声本质来源，考 虑光子统计噪声和系统热噪声统计特性，建立噪声生成机制，通过该噪声生成机制对投影数 据的统计特性进行刻画。优选的，挖掘投影数据的稀疏先验信息中，将稀疏先验信息嵌入噪 声生成机制，并通过数据驱动紧框架挖掘投影数据的稀疏先验信息。通过引入噪声生成机制 以充分刻画投影数据的统计特性并设计正则化项充分挖掘先验信息进而提高模型的去噪以及 细节保留能力。
本发明另一个实施例，建立最大后验概率模型中，采用基于小波变换的数据驱动紧框架 对投影数据进行稀疏编码，获取投影数据和先验信息的联合分布；然后，结合噪声生成机制， 获取后验概率分布；基于该后延概率分布，采用MAP框架估计正弦图，并得到最大后验概 率模型。
上述的，对最大后验概率模型进行迭代求解中，将模型求解等价为带约束的优化问题； 并利用增广拉格朗日函数将带约束的优化问题转化为无约束的优化问题；采用交替方向承子 ADMM法将该无约束的优化问题进行问题分解；根据预先设定的循环迭代条件，对分解后的 问题进行迭代求解，得到去噪后的正弦图数据。用交替方向承子ADMM法将该无约束的优 化问题进行问题分解中，依据接收射线量子、正弦图及尺度参数进行问题分解，获取若干子 问题；然后，根据预先设定的循环迭代条件，对分解后的该若干子问题进行迭代求解，得到 去噪后的正弦图数据。优选的，对分解后的若干子问题迭代求解过程，包含如下内容：初始 化求解参数，获取投影数据和正弦图；依次对子问题、求解参数和小波基函数进行参数更新 求解，并判断是否达到预设最大迭代次数，若是，则判定当前更新求解过程的正弦图数据即 为去噪后的正弦图数据，否则，返回参数更新求解步骤继续执行。
投影重建技术中，FBP算法是在傅立叶变换理论基础之上的一种空域处理技术。它的特 点是在反投影前将每一个采集投影角度下的投影进行卷积处理，从而改善点扩散函数引起的 形状伪影，重建的图像质量较好。本发明另一个实施例中，采用经典解析滤波反投影FBP算 法对去噪后的正弦图数据进行重建，得到重建后的CT图像。
本发明实施例中具体计算过程可设计为如下内容：
参见图2所示，步骤1中：建立噪声生成机制，充分刻画投影数据的噪声统计特性：
在理想情况下，由Beer-Lambert定律可知，投影测量数据的表达式为：
其中
为一维向量，表示正弦图，即测量数据对数变换后的数据，N表示投影路径数目。 I0i(i＝1,…N)表示沿投影路径i的入射X射线的强度。
由上述噪声来源分析可知，测量数据F不可避免的存在噪声，可采用如下模型进行描述：
F＝Q+ε, (2)
其中
表示探测器接收到的射线的量子，
表示系统热噪声，
分别表示实数集 合、自然数集合，一般假设该噪声服从高斯分布，可用下式表示：
其中σ2为噪声的方差。
由于X射线在发射过程中会不可避免地出现量子波动，因此探测器接收到的量子可看作 理想射线被扰动的结果，泊松分布可以较好刻画该过程：
其中O为理想情况下探测器接收到的射线量子，即射线量子的平均数量。令Y为理想的正弦 图，由Beer-Lambert定律可得：
结合式(4)和(5)可得：
由式(3)可知：
结合式(6)和(7)可得：
通过分析投影数据噪声的本质来源，式(8)综合考虑了光子统计噪声和系统热噪声的统 计特性，充分利用了投影数据的统计特性，相比于传统刻画方式，能更充分地描述隐含在正 弦图中的噪声生成机制。将在该噪声机制框架的基础上，设计针对性去噪算法对低剂量投影 数据进行去噪。
步骤2中：采用Hessian惩罚项将稀疏先验模型嵌入到噪声生成机制模型中。并选用基 于学习策略的数据驱动紧框架来挖掘投影数据的稀疏先验信息，建立基于噪声生成机制与数 据驱动紧框架的最大后验概率模型：
对于稀疏先验信息的引入，本文采用Hessian惩罚项将稀疏先验模型嵌入到统计模型中。 通常来讲，可以采用如下公式对数据进行稀疏编码：
其中，
为指数幂分布(也称为超拉普拉斯分布)，b为尺度参数，f:
为 线性算子，cf为与算子f相关的常数。通过采用基于小波变换的数据驱动紧框架对投影数据 进行稀疏编码。令W表示小波算子，则f(Y)＝WY。假设参数b的先验分布为p(b)(通常为 常数函数)，则由式(9)可得数据Y和b的联合分布p(Y,b)：
结合式(8)和式(10)可得到如下后验概率分布：
其中p(F)为与观测值F相关的常数。基于该后验概率分布，我们可以采用MAP框架来估计 正弦图Y，即解决如下优化问题：
步骤3中：采用交替方向乘子法设计高效求解算法：令Z＝WY，则式(12)可改写为：
s.t.WY＝Z.
对于该模型，可等价为求解如下带约束的优化问题：
s.t.WY＝Z.
利用增广Lagrangian函数法将带约束的优化问题(14)转化为如下无约束的优化问题：
其中，
为Lagrangian乘子，μ为非负常数。
采用ADMM方法将Lμ(Q,Y,Z,b,Λ)原问题分解为四个子问题，实现对上述优化模型的求 解，分解得到的四个子问题分别为L(Q)、L(Y)、L(Z)、L(b)。
(1)子问题L(Q)可表示为：
对于每一个Qi可采用下式进行求解：
(2)子问题L(Y)可表示为：
对于该优化问题，通过采用加速一阶方法对其进行求解：给定
采用二次逼近的方法对式 (18)的第二项进行近似，可得：
其中⊙表示逐元素相乘，参数τ可设置为I0的最大值。对式(19)求导并令其导函数为0， 则可得到如下公式：
由数据驱动紧框架的特性可得WTW＝I，此时Y的求解公式为：
(3)子问题L(Z)可表示为：
该式的局部最优解的闭式解为：
其中Sλ(·)为门限算子，定义如下：
其中
(4)子问题L(b)可表示为：
对上式求导并令其导函数为0，可得：
步骤4：实现并运行算法，获得较低噪声的重建图像：现给出基于噪声生成机制与数据 驱动紧框架的投影去噪算法，获得低剂量CT重建图像的具体方法如下：
输入：
I0i，初始μ,b,ρ,Λ＝0和Y＝max{ln(I0)-ln(F),0}，执行如下迭代过程： (1)更新Q子问题，对于每一个Qi，公式如下：
(2)更新Y子问题，公式如下：
(3)更新Z子问题，公式如下：
(4)更新b子问题，公式如下：
(5)参数更新，公式如下：
(6)更新小波基函数，公式如下：
重组操作：
{YY,ZZ,H}←{Yk+1,Zk+1,Wl}, (32)
更新小波基函数：
MΣNT＝ZZYY T,Hl+1←MNT, (33)
更新变量Z：
ZZ←Γλ((Hl+1)TYY), (34)
变量替换：
{Z,W}←{ZZ,Hl+1}, (35)
(7)迭代轮数增加：k←k+1。若“未达到迭代设定的最大次数”则返回第(1)步，否则停 止，迭代达到最大轮数退出时，Y(N)即输出去噪后的正弦图数据，然后采用经典解析算法FBP 对去噪后正弦图进行重建即可得到重建后的CT图像。
基于上述的重建方法，本发明实施例还提供一种基于噪声生成机制与数据驱动紧框架的 CT投影去噪重建装置，参见图3所示，包含：数据处理模块101、模型构建模块102、模型 求解模块103和图像重建模块104，其中，
数据处理模块101，用于刻画投影数据噪声统计特性，并挖掘投影数据的稀疏先验信息；
模型构建模块102，用于结合噪声统计特性和稀疏先验信息，建立最大后验概率模型；
模型求解模块103，用于对最大后验概率模型进行迭代求解，获取去噪后的正弦图数据；
图像重建模块104，用于对去噪后的正弦图数据进行重建，得到重建后的CT图像。
上述的装置中，模型构建模块102包含：稀疏编码单元2001、概率分布计算单元2002 和模型获取单元2003，其中，
稀疏编码单元2001，用于采用基于小波变换的数据驱动紧框架对投影数据进行稀疏编码， 获取投影数据和先验信息的联合分布；
概率分布计算单元2002，用于结合噪声生成机制，获取后验概率分布；
模型获取单元2003，用于基于后延概率分布，采用MAP框架估计正弦图，并得到最大 后验概率模型。
为进一步说明本发明的有效性，本发明实施例中还通过仿真实验对本发明的有效性做进 一步验证，参见图5所示，(a)为参考图像，(b)为低剂量CT重建图像，(c)为本发明实施例技术方案的重建图像；通过图示也可看出，本发明中方案能够在尽可能不破坏投影数据 有用信息的基础上有效抑制噪声，对重建图像的噪声和伪影具有较好的去除效果，在提升低 剂量CT重建图像质量方面具有较高的应用价值
除非另外具体说明，否则在这些实施例中阐述的部件和步骤的相对步骤、数字表达式和 数值并不限制本发明的范围。
基于上述的方法，本发明实施例还提供一种服务器，包括：一个或多个处理器；存储装 置，用于存储一个或多个程序，当所述一个或多个程序被所述一个或多个处理器执行，使得 所述一个或多个处理器实现上述的方法。
基于上述的方法，本发明实施例还提供一种计算机可读介质，其上存储有计算机程序， 其中，该程序被处理器执行时实现上述的方法。
本发明实施例所提供的装置，其实现原理及产生的技术效果和前述方法实施例相同，为 简要描述，装置实施例部分未提及之处，可参考前述方法实施例中相应内容。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统和装置 的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在这里示出和描述的所有示例中，任何具体值应被解释为仅仅是示例性的，而不是作为 限制，因此，示例性实施例的其他示例可以具有不同的值。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附 图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。
附图中的流程图和框图显示了根据本发明的多个实施例的系统、方法和计算机程序产品 的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个 模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现 规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现中，方框中所标注的功 能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地 执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/ 或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动 作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
在本申请所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过 其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅 为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结 合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的 相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或 通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部 件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元 上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个 单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在 一个处理器可执行的非易失的计算机可读取存储介质中。基于这样的理解，本发明的技术方 案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体 现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备 (可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或 部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、 随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码 的介质。
最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技 术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进 行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发 明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变 化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术 方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。 因此，本发明的保护范围应所述以权利要求的保护范围为准。
Claims (6)
1.一种基于噪声生成机制与数据驱动紧框架的CT投影去噪重建方法，其特征在于，包含如下内容：
刻画投影数据噪声统计特性，并挖掘投影数据的稀疏先验信息；
结合噪声统计特性和稀疏先验信息，建立最大后验概率模型；
对最大后验概率模型进行迭代求解，获取去噪后的正弦图数据；
对去噪后的正弦图数据进行重建，得到重建后的CT图像；
刻画投影数据噪声统计特性中，通过分析投影数据噪声本质来源，考虑光子统计噪声和系统热噪声统计特性，建立噪声生成机制，通过该噪声生成机制对投影数据的统计特性进行刻画；
挖掘投影数据的稀疏先验信息中，将稀疏先验信息嵌入噪声生成机制，并通过数据驱动紧框架挖掘投影数据的稀疏先验信息；
建立最大后验概率模型，包含如下内容：首先，采用基于小波变换的数据驱动紧框架对投影数据进行稀疏编码，获取投影数据和先验信息的联合分布；然后，结合噪声生成机制，获取后验概率分布；基于该后延概率分布，采用MAP框架估计正弦图，并得到最大后验概率模型；
对最大后验概率模型进行迭代求解中，将模型求解等价为带约束的优化问题；并利用增广拉格朗日函数将带约束的优化问题转化为无约束的优化问题；采用交替方向承子ADMM法将该无约束的优化问题进行问题分解；根据预先设定的循环迭代条件，对分解后的问题进行迭代求解，得到去噪后的正弦图数据。
2.根据权利要求1所述的基于噪声生成机制与数据驱动紧框架的CT投影去噪重建方法，其特征在于，用交替方向承子ADMM法将该无约束的优化问题进行问题分解中，依据接收射线量子、正弦图及尺度参数进行问题分解，获取若干子问题；然后，根据预先设定的循环迭代条件，对分解后的该若干子问题进行迭代求解，得到去噪后的正弦图数据。
3.根据权利要求1所述的基于噪声生成机制与数据驱动紧框架的CT投影去噪重建方法，其特征在于，对分解后的若干子问题迭代求解过程，包含如下内容：初始化求解参数，获取投影数据和正弦图；依次对子问题、求解参数和小波基函数进行参数更新求解，并判断是否达到预设最大迭代次数，若是，则判定当前更新求解过程的正弦图数据即为去噪后的正弦图数据，否则，返回参数更新求解步骤继续执行。
4.根据权利要求1所述的基于噪声生成机制与数据驱动紧框架的CT投影去噪重建方法，其特征在于，采用滤波反投影FBP算法对去噪后的正弦图数据进行重建，得到重建后的CT图像。
5.一种基于噪声生成机制与数据驱动紧框架的CT投影去噪重建装置，其特征在于，包含：数据处理模块、模型构建模块、模型求解模块和图像重建模块，其中，
数据处理模块，用于刻画投影数据噪声统计特性，并挖掘投影数据的稀疏先验信息；
模型构建模块，用于结合噪声统计特性和稀疏先验信息，建立最大后验概率模型；
模型求解模块，用于对最大后验概率模型进行迭代求解，获取去噪后的正弦图数据；
图像重建模块，用于对去噪后的正弦图数据进行重建，得到重建后的CT图像；
刻画投影数据噪声统计特性中，通过分析投影数据噪声本质来源，考虑光子统计噪声和系统热噪声统计特性，建立噪声生成机制，通过该噪声生成机制对投影数据的统计特性进行刻画；
挖掘投影数据的稀疏先验信息中，将稀疏先验信息嵌入噪声生成机制，并通过数据驱动紧框架挖掘投影数据的稀疏先验信息；
建立最大后验概率模型，包含如下内容：首先，采用基于小波变换的数据驱动紧框架对投影数据进行稀疏编码，获取投影数据和先验信息的联合分布；然后，结合噪声生成机制，获取后验概率分布；基于该后延概率分布，采用MAP框架估计正弦图，并得到最大后验概率模型；
对最大后验概率模型进行迭代求解中，将模型求解等价为带约束的优化问题；并利用增广拉格朗日函数将带约束的优化问题转化为无约束的优化问题；采用交替方向承子ADMM法将该无约束的优化问题进行问题分解；根据预先设定的循环迭代条件，对分解后的问题进行迭代求解，得到去噪后的正弦图数据。
6.根据权利要求5所述的基于噪声生成机制与数据驱动紧框架的CT投影去噪重建装置，其特征在于，模型构建模块包含：稀疏编码单元、概率分布计算单元和模型获取单元，其中，
稀疏编码单元，用于采用基于小波变换的数据驱动紧框架对投影数据进行稀疏编码，获取投影数据和先验信息的联合分布；
概率分布计算单元，用于结合噪声生成机制，获取后验概率分布；
模型获取单元，用于基于后延概率分布，采用MAP框架估计正弦图，并得到最大后验概率模型。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910243267.6A CN110136218B (zh) | 2019-03-28 | 2019-03-28 | 基于噪声生成机制与数据驱动紧框架的ct投影去噪重建方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910243267.6A CN110136218B (zh) | 2019-03-28 | 2019-03-28 | 基于噪声生成机制与数据驱动紧框架的ct投影去噪重建方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110136218A CN110136218A (zh) | 2019-08-16 |
| CN110136218B true CN110136218B (zh) | 2022-10-28 |
Family
ID=67568673
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910243267.6A Active CN110136218B (zh) | 2019-03-28 | 2019-03-28 | 基于噪声生成机制与数据驱动紧框架的ct投影去噪重建方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110136218B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110488350B (zh) * | 2019-09-20 | 2021-10-29 | 西南石油大学 | 基于卷积神经网络的地震反演大数据生成方法 |
| CN111652951B (zh) * | 2020-05-07 | 2023-06-06 | 中国工程物理研究院材料研究所 | 一种稀疏角度快中子ct成像方法 |
| CN111724452B (zh) * | 2020-06-17 | 2024-05-28 | 深圳高性能医疗器械国家研究院有限公司 | 一种低剂量ct图像重建方法 |
| CN112330565A (zh) * | 2020-11-12 | 2021-02-05 | 中国人民解放军战略支援部队信息工程大学 | 基于改善的U-net的低剂量CT投影域中图像去噪方法 |
| CN113011593B (zh) * | 2021-03-15 | 2021-11-02 | 北京百度网讯科技有限公司 | 消除量子测量噪声的方法及系统、电子设备和介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103810733A (zh) * | 2014-02-28 | 2014-05-21 | 南方医科大学 | 一种稀疏角度x射线ct图像的统计迭代重建方法 |
| CN106780641A (zh) * | 2016-11-14 | 2017-05-31 | 西安交通大学 | 一种低剂量x射线ct图像重建方法 |
| JP2017196404A (ja) * | 2016-04-27 | 2017-11-02 | 東芝メディカルシステムズ株式会社 | 画像処理装置、ｘ線ｃｔ装置及び画像処理方法 |
-
2019
- 2019-03-28 CN CN201910243267.6A patent/CN110136218B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103810733A (zh) * | 2014-02-28 | 2014-05-21 | 南方医科大学 | 一种稀疏角度x射线ct图像的统计迭代重建方法 |
| JP2017196404A (ja) * | 2016-04-27 | 2017-11-02 | 東芝メディカルシステムズ株式会社 | 画像処理装置、ｘ線ｃｔ装置及び画像処理方法 |
| CN106780641A (zh) * | 2016-11-14 | 2017-05-31 | 西安交通大学 | 一种低剂量x射线ct图像重建方法 |
Non-Patent Citations (2)
| Title |
|---|
| 低剂量CT的线性Bregman迭代重建算法;王丽艳等;《电子与信息学报》;20131015(第10期);全文 * |
| 基于离散剪切波正则化的低剂量CT图像统计重建算法;张海燕等;《计算机工程与科学》;20180115(第01期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110136218A (zh) | 2019-08-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110136218B (zh) | 基于噪声生成机制与数据驱动紧框架的ct投影去噪重建方法及装置 | |
| EP3195265B1 (en) | Iterative image reconstruction with a sharpness driven regularization parameter | |
| Xu et al. | A practical cone-beam CT scatter correction method with optimized Monte Carlo simulations for image-guided radiation therapy | |
| Jia et al. | GPU-based fast low-dose cone beam CT reconstruction via total variation | |
| EP2958494B1 (en) | Structure propagation restoration for spectral ct | |
| EP2881039A1 (en) | X-ray computer tomography image pick-up device and image reconstruction method | |
| Gong et al. | Adaptive iterative reconstruction based on relative total variation for low-intensity computed tomography | |
| Xie et al. | Deep efficient end-to-end reconstruction (DEER) network for few-view breast CT image reconstruction | |
| EP2831844B1 (en) | Iterative image reconstruction with regularization | |
| Gong et al. | Image reconstruction model for limited-angle CT based on prior image induced relative total variation | |
| Zhang et al. | Limited angle CT reconstruction by simultaneous spatial and Radon domain regularization based on TV and data-driven tight frame | |
| Xu et al. | Quantifying the importance of the statistical assumption in statistical x-ray CT image reconstruction | |
| CN113870138A (zh) | 基于三维U-net的低剂量CT图像去噪方法及系统 | |
| Schmidt et al. | Addressing CT metal artifacts using photon‐counting detectors and one‐step spectral CT image reconstruction | |
| Jin et al. | Bone-induced streak artifact suppression in sparse-view CT image reconstruction | |
| CN106233327B (zh) | 去噪图像数据中对低对比度结构的恢复 | |
| Yu et al. | Weighted adaptive non-local dictionary for low-dose CT reconstruction | |
| Qi et al. | Iterative image reconstruction using modified non-local means filtering for limited-angle computed tomography | |
| Schöndube et al. | Evaluation of a novel CT image reconstruction algorithm with enhanced temporal resolution | |
| Chen et al. | Joint-MAP tomographic reconstruction with patch similarity based mixture prior model | |
| Cierniak | An analytical iterative statistical algorithm for image reconstruction from projections | |
| Li et al. | Adaptive non-local means filtering based on local noise level for CT denoising | |
| Ma et al. | Low dose CT reconstruction assisted by an image manifold prior | |
| Hsieh et al. | An algorithm to estimate the object support in truncated images | |
| Arya et al. | Regularization based modified SART iterative method for CT image reconstruction |
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 |