CN112801886B - 一种基于图小波变换的动态pet图像去噪方法和系统 - Google Patents

一种基于图小波变换的动态pet图像去噪方法和系统 Download PDF

Info

Publication number
CN112801886B
CN112801886B CN202011535207.0A CN202011535207A CN112801886B CN 112801886 B CN112801886 B CN 112801886B CN 202011535207 A CN202011535207 A CN 202011535207A CN 112801886 B CN112801886 B CN 112801886B
Authority
CN
China
Prior art keywords
image
wavelet
graph
scale
matrix
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011535207.0A
Other languages
English (en)
Other versions
CN112801886A (zh
Inventor
盛玉霞
易利群
柴利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University of Science and Engineering WUSE
Original Assignee
Wuhan University of Science and Engineering WUSE
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuhan University of Science and Engineering WUSE filed Critical Wuhan University of Science and Engineering WUSE
Priority to CN202011535207.0A priority Critical patent/CN112801886B/zh
Publication of CN112801886A publication Critical patent/CN112801886A/zh
Application granted granted Critical
Publication of CN112801886B publication Critical patent/CN112801886B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10104Positron emission tomography [PET]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20048Transform domain processing
    • G06T2207/20064Wavelet transform [DWT]

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)

Abstract

本发明提供了一种基于图小波变换的动态PET图像去噪方法和系统,通过充分利用合成图像的结构信息和病灶信息,以及图小波变换的多尺度分解特性,克服了现有技术缺陷,实现了简单、高效、高去噪性能地对动态PET图像去噪的功能。本发明利用合成图像构图保证了在滤除噪声的同时保留了PET图像更多的结构信息和更多的病灶信息;本发明利用图小波变换的多尺度分解特性,将噪声PET图像中的信号和大部分噪声分离开,大大提高了去噪性能;本发明的图小波变换可以通过切比雪夫多项式逼近算法快速实现,很大程度提高了计算效率;与现有去噪方法的实验比较表明,本发明能获得较好的去噪效果。

Description

一种基于图小波变换的动态PET图像去噪方法和系统
技术领域
本发明属于PET图像去噪技术领域,具体涉及一种基于图小波变换的动态PET图像去噪方法和系统。
背景技术
正电子发射计算机断层扫描(Positron Emission Tomography,简称PET)是一种强大的功能成像模式,PET通过使用特定的示踪剂来检测分子水平的活性组织,如新陈代谢,血液流动和神经受体浓度。PET图像可以显示分子水平的示踪剂分布,从而可以对生理和药理过程进行定量测量。它在肿瘤学,心脏病学和神经病学中具有广泛的应用。例如,在放射肿瘤学应用中,PET有助于肿瘤的定位和分期。因此,PET具有巨大的临床应用价值。但是,与MRI和CT相比,PET图像具有高对比度和较低的空间分辨率。另外,受注射剂的剂量和扫描时间等物理因素的限制,最后重建的PET图像信噪比(SNR)低,这会影响其检测和定量准确性,也会影响医生诊断的正确性。这推动了近年来针对PET图像去噪方法的研究。
临床上,高斯滤波器被用于PET图像降噪。这是一种简单的降噪算法,但是它在去噪过程中会平滑重要的图像结构。其他后处理方法,例如双边滤波、各向异性传播和非局部均值滤波等也被用于PET图像去噪。但是由于PET图像中的噪声较复杂,一般是包括高斯噪声和泊松噪声的混合噪声,所以这些自然图像的去噪方法不能很好地直接应用于PET图像。
近年来,图信号处理理论(GSP)已成为图像处理中的有用工具,尤其是在自然图像去噪中。Guo等人(S.Guo,Y.Sheng,L.Chai,J.Zhang,Graph filtering approach to PETimage denoising,International Conference on Industrial ArtificialIntelligence pp.1-6,2019)设计了一种低通图滤波器,对噪声PET图像在图傅里叶域进行低通滤波。该方法得到的去噪PET图像质量相比于高斯滤波和非局部均值滤波都有所提高,但是计算量太大,效率低。
发明内容
本发明要解决的技术问题是:提供一种基于图小波变换的动态PET图像去噪方法和系统,用于简单、高效地对动态PET图像去噪。
本发明为解决上述技术问题所采取的技术方案为:一种基于图小波变换的动态PET图像去噪方法,包括以下步骤:
S1:通过基于合成图像的构图方法构建图加权邻接矩阵A将图像信号转换为图信号;
S2:通过图小波变换对噪声PET图像进行第一次图小波多尺度分解,同时通过切比雪夫多项式逼近快速实现图小波变换;
S3:对第二个尺度的图小波高频系数进行第二次图小波多尺度分解;
S4:对两次分解得到的第一尺度图小波低频系数进行重建得到去噪后的PET图像。
按上述方案,所述的步骤S1中,具体步骤为:
S11:将单个时间帧表示为It,t=1,2,...,T,将不同时间帧相加得到先验的合成图像
Figure GDA0003798074680000021
S12:对大小为N×N的合成图像Ic构造大小为N2×N2的图加权邻接矩阵A,通过阈值高斯核权重函数定义合成图像中的像素点i和像素点j之间的边缘权重:
Figure GDA0003798074680000022
式(1)中:
Figure GDA0003798074680000023
表示以像素点i为中心的图像块,维度为n×n,n为正整数,
Figure GDA0003798074680000024
表示以像素点j为中心的图像块,维度为n×n,n为正整数,
pi表示图像块中心像素点的横坐标,
pj表示图像块中心像素点的纵坐标,
η表示控制图像块相似度参数,η为非负数,
θ表示控制整体相似度的缩放指标,
Figure GDA0003798074680000025
代表距离参数为k的关于图像块xi的最近K邻居;
通过图加权邻接矩阵A构造图拉普拉斯矩阵L:
L=D-A (2),
式(2)中:D表示对角矩阵,大小为N2×N2,且D的每个对角元素dii是图加权邻接矩阵A第i行元素之和。
进一步的,所述的步骤S2中,具体步骤为:
对大小为N×N的噪声PET图像进行第一次图小波变换,图小波尺度为R=[R1,R2,...,RJ],并通过切比雪夫多项式逼近得到大小为N2×(J+1)的图小波系数B1,J+1=[B1,0,B1,1,...,B1,J]:
Figure GDA0003798074680000031
Figure GDA0003798074680000032
Figure GDA0003798074680000033
Figure GDA0003798074680000034
Figure GDA0003798074680000035
Figure GDA0003798074680000036
Figure GDA0003798074680000037
Figure GDA0003798074680000038
Figure GDA0003798074680000039
Figure GDA00037980746800000310
Figure GDA00037980746800000311
式(3)~(13)中:f表示空间域的图像信号,
j表示图小波第j+1个尺度,j=1,2,...,J,
M表示切比雪夫多项式阶数,M为大于0的正整数,
L表示拉普拉斯矩阵,
h(·)表示图谱尺度核函数,
g(·)表示图小波核函数,
γ控制h(0)与图小波核函数g的最大值保持一致,
c0,k表示图谱尺度函数的切比雪夫多项式逼近系数,
cj,k表示图小波函数的切比雪夫多项式逼近系数,
Figure GDA0003798074680000041
代表切比雪夫多项式,
λmax表示拉普拉斯矩阵L进行特征值分解后的最大特征值。
进一步的,所述的步骤S3中,具体步骤为:
重复步骤S2,对第一次分解得到的第二尺度图小波系数B1,j=1进行第二次图小波多尺度分解,得到第二次分解的图小波系数B2,J+1=[B2,0,B2,1,...,B2,J]。
进一步的,所述的步骤S4中,具体步骤为:
S41:构造重建的图小波系数
Figure GDA0003798074680000042
其中b0=B1,0,b1=B2,0,[b2,…,bJ]=0;
S42:进行图小波逆变换,通过共轭梯度方法求解式(14),得到去噪后的PET图像:
Figure GDA0003798074680000043
式(14)中:
Figure GDA0003798074680000044
表示去噪后的PET图像,W表示算子矩阵,且
Figure GDA0003798074680000045
W*表示矩阵W的共轭。
一种基于图小波变换的动态PET图像去噪系统,包括依次连接的邻接矩阵构建模块、图小波变换模块和图小波逆变换模块;邻接矩阵构建模块用于接收动态PET噪声图像并将图像信号转换为图信号;图小波变换模块用于对图信号进行第一次图小波变换,并对第一次图小波变换得到的第二尺度图小波系数进行第二次图小波变换,得到两次分解的图小波系数;图小波逆变换模块用于重建系数得到动态PET去噪图像。
进一步的,还包括合成图像构建模块,合成图像构建模块的输出端连接邻接矩阵构建模块的输入端,用于构建先验的合成图像。
进一步的,还包括系数矩阵构造模块,系数矩阵构造模块串联在图小波变换模块和图小波逆变换模块之间,用于构造图小波系数矩阵。
一种计算机存储介质,其内存储有可被计算机处理器执行的计算机程序,该计算机程序执行一种基于图小波变换的动态PET图像去噪方法。
本发明的有益效果为:
1.本发明的一种基于图小波变换的动态PET图像去噪方法和系统,通过充分利用合成图像的结构信息和病灶信息,以及图小波变换的多尺度分解特性,克服了现有技术缺陷,实现了简单、高效、高去噪性能地对动态PET图像去噪的功能。
2.本发明利用合成图像构图保证了在滤除噪声的同时保留了PET图像更多的结构信息和更多的病灶信息。
3.本发明利用图小波变换的多尺度分解特性,将噪声PET图像中的信号和大部分噪声分离开,大大提高了去噪性能。
4.本发明的图小波变换可以通过切比雪夫多项式逼近算法快速实现,很大程度提高了计算效率;与现有去噪方法的实验比较表明,本发明能获得较好的去噪效果。
附图说明
图1是本发明实施例的流程图。
图2是本发明实施例的功能框图。
图3是本发明实施例的三种不同去噪算法的去噪结果对比图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细的说明。
参见图1,本发明的实施例的一种基于图小波变换的动态PET图像去噪方法,包括以下步骤:
S1:通过基于合成图像的构图方法构建图加权邻接矩阵A将图像信号转换为图信号:
S11:将单个时间帧表示为It,t=1,2,...,24,将不同时间帧相加得到先验的合成图像
Figure GDA0003798074680000061
S12:对大小为175×175的合成图像Ic构造大小为1752×1752的图加权邻接矩阵A,通过阈值高斯核权重函数定义合成图像中的像素点i和像素点j之间的边缘权重:
Figure GDA0003798074680000062
式(1)中:
Figure GDA0003798074680000063
表示以像素点i为中心的图像块,维度为n×n,n=5,
Figure GDA0003798074680000064
表示以像素点j为中心的图像块,维度为n×n,n=5,
pi表示图像块中心像素点的横坐标,
pj表示图像块中心像素点的纵坐标,
η表示控制图像块相似度参数,η=0.0008,
θ表示控制整体相似度的缩放指标,θ=20,
Figure GDA0003798074680000067
代表距离参数为k的关于图像块xi的最近K邻居,k=11;
通过图加权邻接矩阵A构造图拉普拉斯矩阵L:
L=D-A (2),
式(2)中:D表示对角矩阵,大小为1752×1752,且D的每个对角元素dii是图加权邻接矩阵A第i行元素之和;
S2:通过图小波变换对噪声PET图像进行第一次图小波多尺度分解,同时通过切比雪夫多项式逼近快速实现图小波变换:
对大小为175×175的噪声PET图像fin进行第一次图小波变换,图小波尺度为R=[1.2719,0.3719,0.1087,0.0318],并通过切比雪夫多项式逼近得到大小为1752×(4+1)的图小波系数B1,5=[B1,0,B1,1,...,B1,4]:
Figure GDA0003798074680000065
Figure GDA0003798074680000066
Figure GDA0003798074680000071
Figure GDA0003798074680000072
Figure GDA0003798074680000073
Figure GDA0003798074680000074
Figure GDA0003798074680000075
Figure GDA0003798074680000076
Figure GDA0003798074680000077
Figure GDA0003798074680000078
Figure GDA0003798074680000079
式(3)~(13)中:fin表示单帧PET噪声图像信号,
j表示图小波第j+1个尺度,j=1,2,...,4,
M表示切比雪夫多项式阶数,M为大于0的正整数,M=51,
L表示拉普拉斯矩阵,
h(·)表示图谱尺度核函数,
g(·)表示图小波核函数,
γ控制h(0)与图小波核函数g的最大值保持一致,γ=0.9999,
c0,k表示图谱尺度函数的切比雪夫多项式逼近系数,
cj,k表示图小波函数的切比雪夫多项式逼近系数,
Figure GDA00037980746800000710
代表切比雪夫多项式,
λmax表示拉普拉斯矩阵L进行特征值分解后的最大特征值,λmax=31.4489。
S3:对第二个尺度的图小波高频系数进行第二次图小波多尺度分解:
重复步骤S2,对第一次分解得到的第二尺度图小波系数B1,j=1进行第二次图小波多尺度分解,得到第二次分解的图小波系数B2,5=[B2,0,B2,1,...,B2,4];
S4:对两次分解得到的第一尺度图小波低频系数进行重建得到去噪后的PET图像。
S41:构造重建的图小波系数
Figure GDA0003798074680000081
其中b0=B1,0,b1=B2,0,[b2,…,b4]=0;
S42:进行图小波逆变换,通过共轭梯度方法求解式(14),得到去噪后的PET图像:
Figure GDA0003798074680000082
式(14)中:
Figure GDA0003798074680000083
表示去噪后的PET图像,W表示算子矩阵,且
Figure GDA0003798074680000084
W*表示矩阵W的共轭。
本实施例中,使用动态PET图像中的第24帧PET噪声图像进行去噪实验,如图3所示。从图3可以看出,与高斯滤波算法和非局部均值滤波算法相比,基于图小波变换的去噪算法处理后的PET图像保留了更多结构信息,图像细节和病灶更突出,PET图像的质量更高。
以上实施例仅用于说明本发明的设计思想和特点,其目的在于使本领域内的技术人员能够了解本发明的内容并据以实施,本发明的保护范围不限于上述实施例。所以,凡依据本发明所揭示的原理、设计思路所作的等同变化或修饰,均在本发明的保护范围之内。

Claims (9)

1.一种基于图小波变换的动态PET图像去噪方法,其特征在于:包括以下步骤:
步骤S1:通过基于合成图像的构图方法构建图加权邻接矩阵A将噪声PET图像信号转换为图信号,通过图加权邻接矩阵A构造图拉普拉斯矩阵L;
步骤S2:根据步骤S1得到的图拉普拉斯矩阵L,通过图小波变换对噪声PET图像进行第一次图小波多尺度分解,同时通过切比雪夫多项式逼近快速实现图小波变换得到多个尺度系数;
步骤S3:对步骤S2得到的多个尺度系数中的第二尺度图小波高频系数进行第二次图小波多尺度分解得到多个尺度系数;
步骤S4:对两次分解得到的多个尺度系数中的第一尺度图小波低频系数进行重建得到去噪后的PET图像。
2.根据权利要求1所述的一种基于图小波变换的动态PET图像去噪方法,其特征在于:所述的步骤S1中,具体步骤为:
S11:将单个时间帧的噪声PET图像信号表示为It,t=1,2,...,T,将不同时间帧的噪声PET图像信号相加得到合成图像
Figure FDA0003798074670000011
S12:对大小为N×N的合成图像Ic构造大小为N2×N2的图加权邻接矩阵A,通过阈值高斯核权重函数定义合成图像中的像素点i和像素点j之间的边缘权重:
Figure FDA0003798074670000012
式(1)中:
Figure FDA0003798074670000013
表示以像素点i为中心的图像块,维度为n×n,n为正整数,
Figure FDA0003798074670000014
表示以像素点j为中心的图像块,维度为n×n,n为正整数,
pi表示图像块中心像素点的横坐标,
pj表示图像块中心像素点的纵坐标,
η表示控制图像块相似度参数,η为非负数,
θ表示控制整体相似度的缩放指标,
kNN(xi)代表距离参数为k的关于图像块xi的最近K邻居;
通过图加权邻接矩阵A构造图拉普拉斯矩阵L:
L=D-A (2),
式(2)中:D表示对角矩阵,大小为N2×N2,且D的每个对角元素dii是图加权邻接矩阵A第i行元素之和。
3.根据权利要求2所述的一种基于图小波变换的动态PET图像去噪方法,其特征在于:所述的步骤S2中,具体步骤为:
对大小为N×N的噪声PET图像进行第一次图小波变换,图小波尺度为R=[R1,R2,...,RJ],并通过切比雪夫多项式逼近得到大小为N2×(J+1)的图小波系数B1,J+1=[B1,0,B1,1,...,B1,J]:
Figure FDA0003798074670000021
Figure FDA0003798074670000022
Figure FDA0003798074670000023
Figure FDA0003798074670000024
Figure FDA0003798074670000025
Figure FDA0003798074670000026
Figure FDA0003798074670000027
Figure FDA0003798074670000028
Figure FDA0003798074670000029
Figure FDA00037980746700000210
Figure FDA00037980746700000211
式(3)~(13)中:f表示空间域的图像信号,
j表示图小波第j+1个尺度,j=1,2,...,J,
M表示切比雪夫多项式阶数,M为大于0的正整数,
L表示拉普拉斯矩阵,
h(·)表示图谱尺度核函数,
g(·)表示图小波核函数,
γ控制h(0)与图小波核函数g的最大值保持一致,
c0,k表示图谱尺度函数的切比雪夫多项式逼近系数,
cj,k表示图小波函数的切比雪夫多项式逼近系数,
Figure FDA0003798074670000031
代表切比雪夫多项式,
λmax表示拉普拉斯矩阵L进行特征值分解后的最大特征值。
4.根据权利要求3所述的一种基于图小波变换的动态PET图像去噪方法,其特征在于:所述的步骤S3中,具体步骤为:
重复步骤S2,对第一次分解得到的第二尺度图小波系数B1,j=1进行第二次图小波多尺度分解,得到第二次分解的图小波系数B2,J+1=[B2,0,B2,1,...,B2,J]。
5.根据权利要求4所述的一种基于图小波变换的动态PET图像去噪方法,其特征在于:所述的步骤S4中,具体步骤为:
S41:构造重建的图小波系数
Figure FDA0003798074670000032
其中b0=B1,0,b1=B2,0,[b2,…,bJ]=0;
S42:进行图小波逆变换,通过共轭梯度方法求解式(14),得到去噪后的PET图像:
Figure FDA0003798074670000033
式(14)中:
Figure FDA0003798074670000034
表示去噪后的PET图像,W表示算子矩阵,且
Figure FDA0003798074670000035
W*表示矩阵W的共轭。
6.一种基于图小波变换的动态PET图像去噪系统,用于权利要求1至5中任意一项所述的一种基于图小波变换的动态PET图像去噪方法,其特征在于:包括依次连接的邻接矩阵构建模块、图小波变换模块和图小波逆变换模块;
邻接矩阵构建模块用于接收动态PET噪声图像并将图像信号转换为图信号;
图小波变换模块用于对图信号进行第一次图小波变换,并对第一次图小波变换得到的第二尺度图小波系数进行第二次图小波变换,得到两次分解的图小波系数;图小波逆变换模块用于重建系数得到动态PET去噪图像。
7.根据权利要求6所述的去噪系统,其特征在于:还包括合成图像构建模块,合成图像构建模块的输出端连接邻接矩阵构建模块的输入端,用于构建先验的合成图像。
8.根据权利要求6所述的去噪系统,其特征在于:还包括系数矩阵构造模块,系数矩阵构造模块串联在图小波变换模块和图小波逆变换模块之间,用于构造图小波系数矩阵。
9.一种计算机存储介质,其特征在于:其内存储有可被计算机处理器执行的计算机程序,该计算机程序执行如权利要求1至权利要求5中任意一项所述的一种基于图小波变换的动态PET图像去噪方法。
CN202011535207.0A 2020-12-23 2020-12-23 一种基于图小波变换的动态pet图像去噪方法和系统 Active CN112801886B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011535207.0A CN112801886B (zh) 2020-12-23 2020-12-23 一种基于图小波变换的动态pet图像去噪方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011535207.0A CN112801886B (zh) 2020-12-23 2020-12-23 一种基于图小波变换的动态pet图像去噪方法和系统

Publications (2)

Publication Number Publication Date
CN112801886A CN112801886A (zh) 2021-05-14
CN112801886B true CN112801886B (zh) 2022-09-30

Family

ID=75807313

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011535207.0A Active CN112801886B (zh) 2020-12-23 2020-12-23 一种基于图小波变换的动态pet图像去噪方法和系统

Country Status (1)

Country Link
CN (1) CN112801886B (zh)

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103854264B (zh) * 2014-03-28 2016-05-11 中国石油大学(华东) 一种基于改进型阈值函数的小波变换图像去噪方法
CN104182944A (zh) * 2014-08-28 2014-12-03 中国石油大学(华东) 基于曲波与小波变换相串联的光学图像去噪方法

Also Published As

Publication number Publication date
CN112801886A (zh) 2021-05-14

Similar Documents

Publication Publication Date Title
Xu et al. 200x low-dose PET reconstruction using deep learning
Nishio et al. Convolutional auto-encoder for image denoising of ultra-low-dose CT
CN111325686B (zh) 一种基于深度学习的低剂量pet三维重建方法
KR20190137913A (ko) 딥 컨볼루션 신경망을 사용하는 의료 이미징을 위한 선량 감소
CN109584324B (zh) 一种基于自动编码器网络的正电子发射型计算机断层显像(pet)重建方法
CN107705261B (zh) 一种图像重建方法和装置
US11049230B2 (en) Feature-based image processing using feature images extracted from different iterations
CN108287324B (zh) 磁共振多对比度图像的重建方法和装置
Wang et al. Nonlinear filtering based on 3D wavelet transform for MRI denoising
US8989462B2 (en) Systems, methods and computer readable storage mediums storing instructions for applying multiscale bilateral filtering to magnetic resonance (RI) images
Gui et al. PML algorithm for positron emission tomography combined with nonlocal fuzzy anisotropic diffusion filtering
Zhu et al. STEDNet: Swin transformer‐based encoder–decoder network for noise reduction in low‐dose CT
CN112488952A (zh) Pet图像的重建方法及重建终端、计算机可读存储介质
CN112801886B (zh) 一种基于图小波变换的动态pet图像去噪方法和系统
CN102890817A (zh) 医学图像的处理方法及系统
CN109584322B (zh) 基于频域方向平滑Shearlet医学PET图像去噪方法
CN111311531A (zh) 图像增强方法、装置、控制台设备及医学成像系统
CN116485925A (zh) 一种ct图像环状伪影抑制方法、装置、设备和存储介质
CN114913262B (zh) 采样模式与重建算法联合优化的核磁共振成像方法及系统
CN116245969A (zh) 一种基于深度神经网络的低剂量pet图像重建方法
CN116342444A (zh) 一种双通道多模态图像融合方法及融合成像终端机
Mahmoud et al. Variant Wasserstein Generative Adversarial Network Applied on Low Dose CT Image Denoising.
Tripathi et al. A dilated convolution-based denoising network for magnetic resonance images
Angelini et al. Fusion of brushlet and wavelet denoising methods for nuclear images
Xin et al. Guided block matching and 4-D transform domain filter projection denoising method for dynamic PET 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