CN109166161B - 一种基于噪声伪影抑制卷积神经网络的低剂量ct图像处理*** - Google Patents
一种基于噪声伪影抑制卷积神经网络的低剂量ct图像处理*** Download PDFInfo
- Publication number
- CN109166161B CN109166161B CN201810722915.1A CN201810722915A CN109166161B CN 109166161 B CN109166161 B CN 109166161B CN 201810722915 A CN201810722915 A CN 201810722915A CN 109166161 B CN109166161 B CN 109166161B
- Authority
- CN
- China
- Prior art keywords
- image
- dose
- low
- noise artifact
- neural network
- 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
- 238000013527 convolutional neural network Methods 0.000 title claims abstract description 61
- 238000012545 processing Methods 0.000 title claims abstract description 44
- 230000001629 suppression Effects 0.000 title claims abstract description 27
- 238000000034 method Methods 0.000 claims abstract description 41
- 238000012549 training Methods 0.000 claims abstract description 27
- 238000013507 mapping Methods 0.000 claims abstract description 8
- 238000013528 artificial neural network Methods 0.000 claims abstract description 7
- 238000007781 pre-processing Methods 0.000 claims abstract description 6
- 230000006870 function Effects 0.000 claims description 13
- 230000004913 activation Effects 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 4
- 239000006185 dispersion Substances 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000000903 blocking effect Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 9
- 238000003745 diagnosis Methods 0.000 abstract description 6
- 238000004458 analytical method Methods 0.000 abstract description 2
- 238000013170 computed tomography imaging Methods 0.000 abstract 1
- 238000002591 computed tomography Methods 0.000 description 5
- 230000005855 radiation Effects 0.000 description 5
- 210000001519 tissue Anatomy 0.000 description 5
- 210000001015 abdomen Anatomy 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 230000000007 visual effect Effects 0.000 description 4
- 206010028980 Neoplasm Diseases 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000013135 deep learning Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 210000003484 anatomy Anatomy 0.000 description 2
- 201000011510 cancer Diseases 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000000638 solvent extraction Methods 0.000 description 2
- 208000021959 Abnormal metabolism Diseases 0.000 description 1
- 206010011732 Cyst Diseases 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 208000031513 cyst Diseases 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000006371 metabolic abnormality Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 210000000952 spleen Anatomy 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
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/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- 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
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- 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/20084—Artificial neural networks [ANN]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Biophysics (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Health & Medical Sciences (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种基于噪声伪影抑制卷积神经网络的低剂量CT图像处理***,包括:图像预处理模块,用于获得多组匹配的低剂量CT图像和常规剂量CT图像并将和相减获得噪声伪影图像Ns;噪声伪影抑制卷积神经网络建立模块用于将作为训练图像,Ns作为标签图像,建立与Ns之间的映射卷积神经网络;网络训练模块,用于通过降低神经网络损失函数,来训练噪声伪影抑制卷积神经网络;网络处理模块,用于将待处理低剂量CT图像输入映射卷积神经网络进行处理,得到预测的噪声伪影图像噪声伪影抑制模块,用于将减去得到噪声伪影抑制后的图像本发明可有效的抑制低剂量CT数据中的噪声伪影,处理后的图像质量可满足临床分析、诊断等要求,提高低剂量CT成像的图像效果。
Description
技术领域
本发明涉及图像处理技术,尤其涉及一种基于噪声伪影抑制卷积神经网络的低剂量CT图像处理***。
背景技术
X射线计算机断层成像(X-ray Computer Tomography,CT)技术是通过对物体进行射线投影测量而获取物体精确而无损的横截面衰减信息的成像技术,是目前常规有效的临床医学诊断工具之一,为临床医生的诊断和预防提供丰富的三维人体器官组织信息,已成为医学成像领域内不可或缺的检查诊断方法。然而,随着CT断层扫描在临床诊断尤其是常规检查中的普及,CT扫描中的辐射剂量问题已经引起了人们越来越多的关注,大量的临床研究表明超过正常范围的CT辐射剂量易诱发人体新陈代谢异常乃至癌症等疾病,增加人体患癌的风险。但另一方面,为了获取更加清晰的CT图片,提高确诊率,传统的做法是增加扫描时的电流值或者电压值,或者投影的角度,但这将显著增加患者所受的X射线辐射剂量。X射线剂量的控制和临床上对CT图像质量的需求一直以来都是不可避免的矛盾。如何在图像质量无明显下降的基础上,以最低的辐射剂量获得最佳的CT诊断影像已成为行业共识。
当前在提高低剂量CT图像质量的方法主要分为两大类:基于投影空间数据处理的和基于图像空间数据处理。基于投影空间数据的方法主要通过对低剂量CT投影数据的校正,恢复复原和去噪来为重建提供更准确、噪声更少的投影数据,以提高其重建的质量,如结构自适应滤波器和双边滤波器是两种有效的滤波器。还有一些学者和公司通过研究投影空间的数据模型并基于此建立不同算法来抑制低剂量投影数据中的噪声,提高成像质量。另一类方法是直接通过图像空间处理技术来提高已重建的低剂量图像质量,具有不依赖原始投影数据和处理速度快的优点,通常使用非线性的处理方法进行,如小波(Wavelet)变换方法通过保持图像边缘信息来去除伪影及噪声,然而此类方法主要基于图像的局部信息,忽略了图像中重要的非局部性质,也难以达到满意的效果。再如基于字典学习的稀疏表示图像处理算法,这种方法通过训练得到一组过完备字典(基),在基于字典表示的低剂量CT图像处理过程中,通过控制参数,可以使伪影及噪声得不到较好的表示,从而达到去除伪影及噪声的目的。基于字典学习的稀疏表示方法已经被证明在低剂量腹部CT图像处理中具有一定的效果,腹部低剂量CT图像通过训练字典处理后,能够获得较好的效果,能够使病人在腹部CT扫描中把所受剂量降低到原来的五分之一。但是这种方法需要训练不同方向的高频细节图像,缺乏断层图像之间的相关性,难以扩展到三维体数据的处理,计算量过大且耗时,难以在实际三维医疗图像处理***中广泛应用。
随着大数据集、大样本的普遍,深度学习在工业界及学术界都受到了广泛的关注,也逐渐应用于CT图像领域中。如Chen等人采用了一种RED-CNN编码网络,大幅度降低了CT图像中的噪声伪影,提高了肿瘤病变组织的识别率。基于深度学习方法对图像空间建模,通过卷积神经网络强大的特征表示能力来分辨伪影及噪声与人体解剖结构。基于深度学习的方法,测试时间短、处理效果好,是数据量充足时应首先考虑的算法。为此,本发明将提出一种基于卷积神经网络的方法,充分利用卷积网络强大的表示能力,实现低剂量CT图像噪声伪影及解剖结构之间的分解。
虽然RED-CNN并与传统的方法相比不管是在主观视觉效果上,还是在峰值信噪比、结构相似性和均方根误差等客观评价指标上,都是最优的,对低剂量CT图像的去噪效果达到了目前较先进的水平,但RED-CNN的网络复杂度较高,运算耗时较长。
发明内容
发明目的:本发明针对现有技术存在的问题,提供一种基于噪声伪影抑制卷积神经网络的低剂量CT图像处理***,该***与RED-CNN相比,可以大幅度降低网络参数数量,同时具有较好去噪去伪影效果。
技术方案:本发明所述的基于噪声伪影抑制卷积神经网络的低剂量CT图像处理***包括:
图像预处理模块,用于获得多组匹配的低剂量CT图像Vs ld和常规剂量CT图像Vs rd,并将低剂量CT图像和常规剂量CT图像相减获得噪声伪影图像Ns=Vs ld-Vs rd;
噪声伪影抑制卷积神经网络建立模块,用于将低剂量CT图像Vs ld作为训练图像,噪声伪影图像Ns作为标签图像,建立低剂量CT图像Vs ld与噪声伪影图像Ns之间的映射卷积神经网络,作为噪声伪影抑制卷积神经网络;
网络训练模块,用于将多个低剂量CT图像Vs ld作为输入,对应的噪声伪影图像作为输出,通过降低神经网络损失函数,来训练和学习噪声伪影抑制卷积神经网络中的参数;
进一步的,所述图像预处理模块具体包括:
低剂量CT图像处理单元,用于将低剂量扫描下的CT投影数据通过解析的FBP重建算法获得低剂量CT图像Vs ld;
常规剂量CT图像处理单元,用于将常规剂量扫描下的CT投影数据经过GDSIR迭代重建算法获得常规剂量CT图像Vs rd。
进一步的,所述常规剂量CT图像处理单元具体用于执行如下步骤:
获取常规剂量扫描下的CT投影数据;
建立重建目标函数为:
式中,G为投影矩阵,u为重建CT图像,W为CT投影数据的统计权重,按照CT投影数据p的方差来计算;||·||w为加权L2范数,||·||0为L0范数,λ为正则化参数,Ei为在重建图像u的i个位置选取三维图像块的操作,I为图像块总数,D为特征字典,αi为第i个图像块的稀疏表示向量,Mi为第i个图像块的均值,T为稀疏度,ε为误差限度;
对上述重建目标函数进行迭代,从而得到常规剂量CT图像数据Vs rd。
进一步的,所述噪声伪影抑制卷积神经网络包含CBR模块、分支模块和残差模块。
进一步的,所述CBR模块用于对输入的提取低剂量CT图像进行卷积、尺度变化和ReLU激活,以提取低剂量CT图像的低层特征信息。
进一步的,所述分支模块用于通过增加网络宽度混合CBR模块提取的特征信息,以提高网络的表示能力。
进一步的,所述残差模块通过卷积、尺度变化和ReLU激活来缩短训练时间,降低同等表示能力下特征卷积核的冗余度,避免训练中梯度弥散。
进一步的,所述损失函数Loss为:
有益效果:本发明与现有技术相比,其显著优点是:本发明可以有效的处理高噪声污染下低剂量CT图像中星条状伪影及噪声和特征结构成分,处理效果优于目前图像空间神经网络方法,其处理后的低剂量CT图像质量能够更好的满足临床分析和诊断的要求,为CT扫描中的降低对扫描者的辐射剂量伤害做出贡献。
附图说明
图1为本发明实施的方法流程图;
图2为本发明实施例中五幅典型的轴向训练图(第一行:常规剂量图;第二行:低剂量图);
图3为本发明实施例中轴向常规剂量CT图像和低剂量CT图像(a:常规剂量;b:低剂量);
图4为本发明实施例中轴向临床低剂量CT图像使用RED-CNN方法和使用本发明方法NAS-CNN处理后的结果(a:RED-CNN;b:NAS-CNN);
图5为本发明实施例中轴向临床低剂量CT图像对比实验中结果的局部放大图(a、b、c、d分别为常规剂量CT图像、低剂量CT图像、RED-CNN处理后的结果和NAS-CNN处理后的结果);
图6为本发明实施例中矢向常规剂量CT图像和冠向低剂量CT图像(a:常规剂量;b:低剂量);
图7为本发明实施例中矢向临床低剂量CT图像使用RED-CNN方法和使用本发明方法NAS-CNN处理后的结果的结果(a:RED-CNN;b:NAS-CNN);
图8为本发明实施例中矢向临床低剂量CT图像对比实验中结果的局部放大图(a、b、c、d分别为常规剂量CT图像、低剂量CT图像、RED-CNN处理后的结果和NAS-CNN处理后的结果)。
具体实施方式
本实施例提供了一种基于噪声伪影抑制卷积神经网络的低剂量CT图像处理***,如图1所示,包括图像预处理模块、噪声伪影抑制卷积神经网络建立模块、网络训练模块、网络处理模块和噪声伪影抑制模块,下面对各模块进行详细介绍。
图像预处理模块具体包括低剂量CT图像处理单元和常规剂量CT图像处理单元,低剂量CT图像处理单元用于将低剂量扫描下的CT投影数据通过解析的FBP重建算法获得低剂量CT图像Vs ld;常规剂量CT图像处理单元用于将常规剂量扫描下的CT投影数据经过GDSIR(Global Dictionary Based Statistical Iterative Reconstruction)迭代重建算法获得常规剂量CT图像Vs rd,具体执行如下步骤:
A、获取常规剂量扫描下的CT投影数据;
B、建立重建目标函数为:
式中,G为投影矩阵,u为重建CT图像,W为CT投影数据的统计权重,按照CT投影数据p的方差来计算;||·||w为加权L2范数,||·||0为L0范数,λ为正则化参数,Ei为在重建图像u的i个位置选取三维图像块的操作,I为图像块总数,D为特征字典,αi为第i个图像块的稀疏表示向量,Mi为第i个图像块的均值,T为稀疏度,ε为误差限度;
C、对上述重建目标函数进行迭代,从而得到常规剂量CT图像数据Vs rd。
噪声伪影抑制卷积神经网络建立模块用于将低剂量CT图像Vs ld作为训练图像,噪声伪影图像Ns作为标签图像,建立低剂量CT图像Vs ld与噪声伪影图像Ns之间的映射卷积神经网络,作为噪声伪影抑制卷积神经网络。所述噪声伪影抑制卷积神经网络包含CBR模块、分支模块和残差模块。CBR模块用于对输入的提取低剂量CT图像进行卷积、尺度变化和ReLU激活,以提取低剂量CT图像的低层特征信息。分支模块用于通过增加网络宽度混合CBR模块提取的特征信息,以提高网络的表示能力。残差模块通过卷积、尺度变化和ReLU激活来缩短训练时间,降低同等表示能力下特征卷积核的冗余度,避免训练中梯度弥散。
网络训练模块用于将多个低剂量CT图像Vs ld作为输入,对应的噪声伪影图像作为输出,通过降低神经网络损失函数,来训练和学习噪声伪影抑制卷积神经网络中的参数。具体的,将训练集中的低剂量CT图像Vs ld和噪声伪影图像Ns按照一定的尺寸n×n×t和像素间隔l1×l2×l3进行分块(如:图块尺寸为65×65×32像素,分块间隔为12×12×12像素),分别得到图像块集Ps ld和Ps N,特别的,这里的两种图像块选取的位置是完全匹配的。将图像块集Ps ld和Ps N放入网络中,通过降低神经网络损失函数Loss,也就是预测的噪声伪影图像块集与实际噪声伪影图像块集Ps N的欧氏距离,来训练和学习NAS-CNN中的参数,最终得到泛化能力较强的神经网络。损失函数Loss定义为:
网络处理模块用于将待处理低剂量CT图像Vt ld输入已训练好的映射卷积神经网络进行处理,得到预测的噪声伪影图像。具体的,首先,将需要处理的低剂量CT图像Vt ld按照尺寸为n×n×t和像素间隔为l1×l2×l3进行分块,得到图像块集Pt ld;然后,将Pt ld输入训练完成的NAS-CNN中,得到预测的噪声伪影图像块集/>接下来,按照像素间隔l1×l2×l3将图像块集/>组合成噪声伪影图像/>
下面对本发明方法进行评估。
1、效果评估准则
首先获得多组腹部数据,实验中所使用Low Dose Challenge比赛公布的数据,来自Somatom Definition AS+CT设备,具体的扫描参数为:管电压100KVp,管电流为360mAs(常规剂量)/85mAs(低剂量),探测器大小为736×64,每个探测器单元尺寸为1.2856×1.0947mm2,射线源到物体中心和探测器中心的距离分别为59.5cm和108.56cm,全角度模式下每圈采集1152个投影数据,螺距为0.6,其它参数采用机器默认值。分别通过FDK(Feldkamp,Davis,Kress Algorithm)和GDSIR重建后得到重建图像,重建图像大小为512×512,像素大小为0.8×0.8mm2,层厚为1mm,三维连续性较好。
选取其中九组扫描数据作为训练数据,其中五幅典型的轴向训练图如图2所示;一组扫描数据作为测试数据,其中所选取的轴向CT图像如图3所示(常规剂量与低剂量)和所选取的矢向CT图像如图6所示(常规剂量与低剂量)。低剂量CT图像、正常剂量CT图像及处理后图像的显示窗宽为300HU(Housfield Units,HU),窗位为50HU。
2、视觉评估
通过观察图3、图4、图5的常规剂量和低剂量的CT图像,以及RED-CNN方法和本发明NAS-CNN方法处理后的图像,可以看到RED-CNN方法虽然能去除低剂量CT图像中的噪声和条状伪影,但在处理过程中,解剖结构成分丢失了部分组织细节,部分区域有一定的模糊,像肝脏、脾静脉血管和血管囊肿区域;而使用本发明方法处理后的CT图像,视觉效果有了明显改善,图像组织能够很好的保留下来,噪声伪影基本上被完全去除掉,处理后的图像,图像视觉纹理更加接近于常规剂量下CT图像。
3、量化评估
为了量化的验证本发明方法的有效性,我们通过计算比较了整张图像(低剂量CT图、RED-CNN处理后的图像和本发明NAS-CNN处理后的图像)与常规剂量CT图的峰值信噪比与结构相似度,这里峰值信噪比PSNR的定义为:
其中此处I代表正常剂量CT图像,K代表含此处代表待计算的图像,LI代表图像I可代表的最大图像像素值,i,j分别为图像的像素索引。结构相似度SSIM的定义为:
其中μI、μK分别是图像I,K的均值,σI、σK分别是I,K的标准差,σIK是图像I与K的协方差,C1和C2为两个常数,其中C1=(0.01×L)2,C2=(0.03×L)2。从下表1可以看到本发明的分解方法能够大幅度降低分解后解剖结构成分中的噪声,提高信噪比,获得更加接近正常剂量的CT图像。
表1
为对比了RED-CNN与本发明所提NAS-CNN网络的计算耗时,我们在GPU为NVIDIAGTX1080的环境下使用TensorFlow平台,对512×512×200大小的测试图进行计算,发现RED-CNN平均耗时3.72秒,而本发明NAS-CNN平均耗时1.08秒。
从上述实验可以看到,采用本发明的方法可以有效快速的处理低剂量CT图像,可获得接近正常剂量水平的CT信息的人体解剖组织衰减图像。而且本发明方法中,网络一旦构建好后,无需重复训练,实际处理时间短,速度快,具有较大的应用范围。
以上所揭露的仅为本发明一种较佳实施例而已,不能以此来限定本发明之权利范围,因此依本发明权利要求所作的等同变化,仍属本发明所涵盖的范围。
Claims (2)
1.一种基于噪声伪影抑制卷积神经网络的低剂量CT图像处理***,其特征在于,包括:
图像预处理模块,用于获得多组匹配的低剂量CT图像Vs ld和常规剂量CT图像Vs rd,并将低剂量CT图像和常规剂量CT图像相减获得噪声伪影图像Ns=Vs ld-Vs rd;具体包括:
低剂量CT图像处理单元,用于将低剂量扫描下的CT投影数据通过解析的FBP重建算法获得低剂量CT图像Vs ld;
常规剂量CT图像处理单元,用于将常规剂量扫描下的CT投影数据经过GDSIR迭代重建算法获得常规剂量CT图像Vs rd,具体用于执行如下步骤:
获取常规剂量扫描下的CT投影数据;
建立重建目标函数为:
式中,G为投影矩阵,u为重建CT图像,W为CT投影数据的统计权重,按照CT投影数据p的方差来计算;||·||w为加权L2范数,||·||0为L0范数,λ为正则化参数,Ei为在重建图像u的i个位置选取三维图像块的操作,I为图像块总数,D为特征字典,αi为第i个图像块的稀疏表示向量,Mi为第i个图像块的均值,T为稀疏度,ε为误差限度;
对上述重建目标函数进行迭代,从而得到常规剂量CT图像数据Vs rd;
噪声伪影抑制卷积神经网络建立模块,用于将低剂量CT图像Vs ld作为训练图像,噪声伪影图像Ns作为标签图像,建立低剂量CT图像Vs ld与噪声伪影图像Ns之间的映射卷积神经网络,作为噪声伪影抑制卷积神经网络,所述噪声伪影抑制卷积神经网络包含CBR模块、分支模块和残差模块,所述CBR模块用于对输入的提取低剂量CT图像进行卷积、尺度变化和ReLU激活,以提取低剂量CT图像的低层特征信息,所述分支模块用于通过增加网络宽度混合CBR模块提取的特征信息,以提高网络的表示能力,所述残差模块通过卷积、尺度变化和ReLU激活来缩短训练时间,降低同等表示能力下特征卷积核的冗余度,避免训练中梯度弥散;
网络训练模块,用于将多个低剂量CT图像Vs ld作为输入,对应的噪声伪影图像作为输出,通过降低神经网络损失函数,来训练和学习噪声伪影抑制卷积神经网络中的参数;
噪声伪影抑制模块,用于将待处理低剂量CT图像Vt ld减去预测的噪声伪影图像Nt p,得到噪声伪影抑制后的图像Vt p。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810722915.1A CN109166161B (zh) | 2018-07-04 | 2018-07-04 | 一种基于噪声伪影抑制卷积神经网络的低剂量ct图像处理*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810722915.1A CN109166161B (zh) | 2018-07-04 | 2018-07-04 | 一种基于噪声伪影抑制卷积神经网络的低剂量ct图像处理*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109166161A CN109166161A (zh) | 2019-01-08 |
CN109166161B true CN109166161B (zh) | 2023-06-30 |
Family
ID=64897375
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810722915.1A Active CN109166161B (zh) | 2018-07-04 | 2018-07-04 | 一种基于噪声伪影抑制卷积神经网络的低剂量ct图像处理*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109166161B (zh) |
Families Citing this family (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7237624B2 (ja) * | 2019-02-07 | 2023-03-13 | 浜松ホトニクス株式会社 | 画像処理装置および画像処理方法 |
CN110070510A (zh) * | 2019-04-26 | 2019-07-30 | 东北大学 | 一种基于vgg-19提取特征的cnn医学图像降噪方法 |
CN110378982B (zh) * | 2019-07-23 | 2023-09-12 | 上海联影医疗科技股份有限公司 | 重建图像处理方法、装置、设备及存储介质 |
JP7245740B2 (ja) * | 2019-07-25 | 2023-03-24 | 富士フイルムヘルスケア株式会社 | 画像処理装置、画像処理方法及びx線ct装置 |
CN110570492B (zh) * | 2019-09-11 | 2021-09-03 | 清华大学 | 一种基于神经网络的ct伪影抑制方法、设备以及介质 |
CN111009019B (zh) * | 2019-09-27 | 2021-07-16 | 北京航空航天大学 | 基于深度学习的微分相衬ct不完备数据重建方法 |
CN112581554B (zh) * | 2019-09-30 | 2024-02-27 | 中国科学院深圳先进技术研究院 | 一种ct成像方法、装置、存储设备及医学成像*** |
CN110677649B (zh) * | 2019-10-16 | 2021-09-28 | 腾讯科技(深圳)有限公司 | 基于机器学习的去伪影方法、去伪影模型训练方法及装置 |
CN111080736B (zh) * | 2019-12-11 | 2023-09-08 | 电子科技大学 | 一种基于稀疏变换的低剂量ct图像重建方法 |
CN111179366B (zh) * | 2019-12-18 | 2023-04-25 | 深圳先进技术研究院 | 基于解剖结构差异先验的低剂量图像重建方法和*** |
CN110992295B (zh) * | 2019-12-20 | 2022-04-19 | 电子科技大学 | 基于小波-red卷积神经网络的低剂量ct重建方法 |
CN111325686B (zh) * | 2020-02-11 | 2021-03-30 | 之江实验室 | 一种基于深度学习的低剂量pet三维重建方法 |
CN111325737B (zh) * | 2020-02-28 | 2024-03-15 | 上海志唐健康科技有限公司 | 低剂量ct图像处理方法、装置和计算机设备 |
CN111325695B (zh) * | 2020-02-29 | 2023-04-07 | 深圳先进技术研究院 | 基于多剂量等级的低剂量图像增强方法、***及存储介质 |
CN111445406B (zh) * | 2020-03-24 | 2023-05-05 | 广东工业大学 | 一种低剂量ct图片质量改善方法、***及设备 |
CN111968195B (zh) * | 2020-08-20 | 2022-09-02 | 太原科技大学 | 用于低剂量ct图像降噪及去伪影的双注意力生成对抗网络 |
CN112116677B (zh) * | 2020-09-23 | 2024-01-23 | 赣南师范大学 | 一种基于低维流形先验的低剂量ct重建方法 |
CN112330575B (zh) * | 2020-12-03 | 2022-10-14 | 华北理工大学 | 一种卷积神经网络医学ct图像去噪方法 |
CN112446840B (zh) * | 2020-12-07 | 2024-01-19 | 明峰医疗***股份有限公司 | 一种基于深度学习的ct图像黑带伪影消除方法及*** |
CN112598759B (zh) * | 2020-12-15 | 2022-09-13 | 太原科技大学 | 抑制低剂量ct图像中伪影噪声的多尺度特征生成对抗网络 |
CN112767273B (zh) * | 2021-01-21 | 2023-10-20 | 中山大学 | 一种应用特征解耦的低剂量ct图像修复方法及*** |
CN112927318B (zh) * | 2021-02-22 | 2022-04-22 | 明峰医疗***股份有限公司 | 低剂量pet图像的降噪重建方法及计算机可读存储介质 |
CN112927230A (zh) * | 2021-04-26 | 2021-06-08 | 常州市第一人民医院 | 一种影像图像的处理分析*** |
CN113298167A (zh) * | 2021-06-01 | 2021-08-24 | 北京思特奇信息技术股份有限公司 | 一种基于轻量级神经网络模型的文字检测方法及*** |
CN113256529B (zh) * | 2021-06-09 | 2021-10-15 | 腾讯科技(深圳)有限公司 | 图像处理方法、装置、计算机设备及存储介质 |
CN113450427B (zh) * | 2021-06-29 | 2023-09-01 | 深圳高性能医疗器械国家研究院有限公司 | 基于联合字典学习和深度网络的pet图像重建方法 |
CN113379868A (zh) * | 2021-07-08 | 2021-09-10 | 安徽工程大学 | 基于卷积稀疏编码网络的低剂量ct图像噪声伪影分解方法 |
CN113570705B (zh) * | 2021-07-28 | 2024-04-30 | 广州瑞多思医疗科技有限公司 | 一种三维剂量重建方法、装置、计算机设备及存储介质 |
CN113436118B (zh) * | 2021-08-10 | 2022-09-27 | 安徽工程大学 | 一种基于多尺度卷积编码网络的低剂量ct图像复原方法 |
CN116433785A (zh) * | 2021-12-31 | 2023-07-14 | 中国石油天然气股份有限公司 | Ct图像环状伪影的处理方法、装置、***及存储介质 |
CN115828396B (zh) * | 2023-01-03 | 2023-04-28 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 基于混合约束条件下的民宿平面自动生成方法及*** |
CN117315063A (zh) * | 2023-09-07 | 2023-12-29 | 先进能源科学与技术广东省实验室 | 基于深度学习的低剂量ct图像重建方法和*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103473745A (zh) * | 2013-09-16 | 2013-12-25 | 东南大学 | 一种基于区别性字典的低剂量ct图像处理方法 |
CN105118066A (zh) * | 2015-09-16 | 2015-12-02 | 东南大学 | 一种基于三维区别性特征表示的低剂量ct图像分解方法 |
CN108122265A (zh) * | 2017-11-13 | 2018-06-05 | 深圳先进技术研究院 | 一种ct重建图像优化方法及*** |
-
2018
- 2018-07-04 CN CN201810722915.1A patent/CN109166161B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103473745A (zh) * | 2013-09-16 | 2013-12-25 | 东南大学 | 一种基于区别性字典的低剂量ct图像处理方法 |
CN105118066A (zh) * | 2015-09-16 | 2015-12-02 | 东南大学 | 一种基于三维区别性特征表示的低剂量ct图像分解方法 |
CN108122265A (zh) * | 2017-11-13 | 2018-06-05 | 深圳先进技术研究院 | 一种ct重建图像优化方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN109166161A (zh) | 2019-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109166161B (zh) | 一种基于噪声伪影抑制卷积神经网络的低剂量ct图像处理*** | |
CN108961237B (zh) | 一种基于卷积神经网络的低剂量ct图像分解方法 | |
Lu et al. | Iterative reconstruction of low-dose CT based on differential sparse | |
CN109102550B (zh) | 基于卷积残差网络的全网络低剂量ct成像方法及装置 | |
Yang et al. | Low-dose CT image denoising using a generative adversarial network with Wasserstein distance and perceptual loss | |
Yin et al. | Domain progressive 3D residual convolution network to improve low-dose CT imaging | |
Kuanar et al. | Low dose abdominal CT image reconstruction: an unsupervised learning based approach | |
CN112348936B (zh) | 一种基于深度学习的低剂量锥束ct图像重建方法 | |
CN105118066B (zh) | 一种基于三维区别性特征表示的低剂量ct图像分解方法 | |
CN106373163B (zh) | 一种基于三维投影图区别性特征表示的低剂量ct成像方法 | |
CN110599530B (zh) | 基于双正则约束的mvct图像纹理增强方法 | |
CN103226815A (zh) | 一种低剂量ct图像滤波方法 | |
CN113706643B (zh) | 一种基于同形适应学习的头部ct金属伪影校正方法 | |
CN108846813A (zh) | 基于mfdf分解框架与nsst的医学ct图像去噪方法 | |
CN111080736B (zh) | 一种基于稀疏变换的低剂量ct图像重建方法 | |
CN116167929A (zh) | 基于残差多尺度特征提取的低剂量ct图像去噪网络 | |
Wu et al. | Unsharp structure guided filtering for self-supervised low-dose CT imaging | |
Liu et al. | Cascade resunet with noise power spectrum loss for low dose ct imaging | |
Ma et al. | Low dose CT reconstruction assisted by an image manifold prior | |
Yim et al. | Limited-angle CT reconstruction via data-driven deep neural network | |
Chen et al. | Low-dose dental CT image enhancement using a multiscale feature sensing network | |
CN115731158A (zh) | 一种基于残差域迭代优化网络的低剂量ct重建方法 | |
CN115456890A (zh) | 基于多尺度双域判别器的生成对抗医学ct图像去噪方法 | |
Slavine et al. | Phantom and preclinical studies for image improvement in clinical CT | |
KR102329938B1 (ko) | 뉴럴 네트워크를 이용한 콘빔 단층촬영 영상 처리 방법 및 그 장치 |
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 |