CN110866888A - 基于潜在信息表征gan的多模态mri合成方法 - Google Patents

基于潜在信息表征gan的多模态mri合成方法 Download PDF

Info

Publication number
CN110866888A
CN110866888A CN201911114218.9A CN201911114218A CN110866888A CN 110866888 A CN110866888 A CN 110866888A CN 201911114218 A CN201911114218 A CN 201911114218A CN 110866888 A CN110866888 A CN 110866888A
Authority
CN
China
Prior art keywords
image
potential information
network
modality
potential
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.)
Granted
Application number
CN201911114218.9A
Other languages
English (en)
Other versions
CN110866888B (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.)
Sichuan University
Original Assignee
Sichuan University
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 Sichuan University filed Critical Sichuan University
Priority to CN201911114218.9A priority Critical patent/CN110866888B/zh
Publication of CN110866888A publication Critical patent/CN110866888A/zh
Application granted granted Critical
Publication of CN110866888B publication Critical patent/CN110866888B/zh
Expired - Fee Related 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/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/25Fusion techniques
    • G06F18/253Fusion techniques of extracted features
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Software Systems (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Mathematical Physics (AREA)
  • Computational Linguistics (AREA)
  • Health & Medical Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明公开基于潜在信息表征GAN的多模态MRI合成方法,包括步骤:将不同的MRI模态的集体信息输入生成网络;生成网络通过使用不同的编码器分别提取多种MRI模态的潜在信息表征;提取的潜在信息表征进一步传输到潜在空间处理网络进行集成处理;由解码器获得相应的合成目标模态,作为合成图像;将合成图像和真实图像共同输入至鉴别网络;由鉴别网络将真实图像与生成的图像区分开。本发明能够灵活的接收多个输入模态,并将多输入模态进行合成,能够有效避免信息丢失,有效提高合成图像的保真度,获得真实体现检测部位的高品质图像;适用范围宽,计算效率高,且实际运用效果好。

Description

基于潜在信息表征GAN的多模态MRI合成方法
技术领域
本发明属于图像处理技术领域,特别是涉及基于潜在信息表征GAN的多模态MRI合成方法。
背景技术
磁共振成像(MRI)作为一种无创性成像技术,已成为研究神经解剖学的主要成像方式。可以通过使用不同的脉冲序列和参数来产生不同的组织对比图像,从而形成同一解剖结构的不同模态。但是,在实践中,由于时间限制,会造成某些序列发生丢失;某些模态会受到随机噪声和非预期的图像伪影的影响,从而导致图像质量较差。上述因素使得对于同一被检测人员可用的对比图像的数量总是受到限制。因此,需要利用已经成功获得的其他模态合成丢失或损坏的模态;合成后的模态不仅可以替代丢失或损坏的模态,而且具有改善其他图像分析任务的潜在价值。
目前,有许多用于MRI图像的单模态合成方法。需要对不同的病理特征采用其特定的方式,使用各种单模态合成方式提供的有效信息可以从理论上提高合成图像的质量并提高诊断的实用性,但是适用范围窄,实际运用效果差。而多模态合成的方向比使用单模态合成方法获得了更好的合成结果。目前对于多模态合成方法研究较少,且仅有的方式仍然存在很多缺点,会导致信息丢失和综合性能欠佳,造成合成图像不能够真实有效体现检测部位的形态。
作为最流行的深度学习技术之一,生成对抗网络(GAN)已应用于图像合成领域。但是,现有方式在将GAN扩展到多模式转换时,采用的方式是将不同的模态堆叠为输入,以单模态形式优化网络。由于不同的模式呈现出不同的物理状态,因此这些方法大大降低了从每种方式中提取信息的效率,导致合成结果不理想。
发明内容
为了解决上述问题,本发明提出了基于潜在信息表征GAN的多模态MRI合成方法,能够灵活的接收多个输入模态,并将多输入模态进行合成,能够有效避免信息丢失,有效提高合成图像的保真度,获得真实体现检测部位的高品质图像;适用范围宽,计算效率高,且实际运用效果好。
为达到上述目的,本发明采用的技术方案是:基于潜在信息表征GAN的多模态MRI合成方法,包括步骤:
S100,将不同的MRI模态的集体信息输入生成网络;
S200,生成网络通过使用不同的编码器分别提取多种MRI模态的潜在信息表征;提取的潜在信息表征进一步传输到潜在空间处理网络进行集成处理;由解码器获得相应的合成目标模态,作为合成图像;
S300,将合成图像和真实图像共同输入至鉴别网络;
S400,由鉴别网络将真实图像与生成的图像区分开。
进一步的是,为了从多个模态中提取特定信息,在生成网络中包括有多个相互独立的编码器,每个编码器对应输入一个MRI模态并提取该模态的潜在信息表征;编码器的数量取决于输入模态的数量。编码器在保留模态特定信息的潜在信息表征的同时保证了模态不变性;与使用相同的卷积核来提取不同模态的特征相比,使用特定的编码器分别提取每个模态的潜在表示可以提取更有效的多种模态的潜在信息表征进行合成。
进一步的是,所述编码器包括3个卷积块,每个卷积块结构相同,所述卷积块依次进行填充、归一化、激活和卷积处理;
在第一个卷积块中采用3×3的零填充,进行填充;在其余2个卷积块中应用1×1的零填充,进行填充;
在所述归一化处理过程中对数据进行实例归一化处理,即对单个输入图像的全局信息进行归一化调整;此方式能够有效消除由于混洗导致的平均值归一化和批量归一化的标准偏差,并减少了噪声的引入;所述激活函数采用LeakyReLU激活函数,使网络更易于训练;将所述卷积层的步幅大小设置为2,以将特征图的尺寸减半,这样的处理方式不仅可以压缩图像大小,还可以减少最大池化层导致的详细信息丢失。
进一步的是,在所述生成网络中,将编码器所获得的潜在信息表征表示为LR(·)=E*(·|θ);
其中,E*(·|θ)为具有可学习参数θ的编码器函数;以独立地获取输入模式下的潜在信息表征。
进一步的是,所述潜在空间处理网络通过残差网络将每个编码器生成的潜在信息表征集成到单个潜在信息表征中:首先,直接连接不同模态下的潜在信息表征;然后,利用残差块的特性,通过特征映射的方式整合所述潜在信息表征,构成合成图像的潜在信息表征;借助残差网络,能够在不丢失信息的情况下成功地在各个卷积层中传输此类潜在表示,并最终完成多个潜在表示的融合;能够实现灵活地接收多种输入模态并保留其所有潜在信息表征;
所述残差块包括四个残差网络,每个所述残差网络包括镜像填充、批量归一化、LeakyReLU激活和卷积处理。
进一步的是,所述合成图像的潜在信息表征,表示为:
Figure BDA0002273583980000031
其中,R*(·|ψ)为潜在空间处理中具有可学习参数ψ的残差集成函数;LR(T1)为T1输入模态的潜在信息表征,LR(T2)为T2输入模态的潜在信息表征,n为输入模态的个数,FLAIR为目标模态。
进一步的是,将所述合成图像的潜在信息表征通过解码器进行解码,获得相应的合成目标模态作为合成图像;
解码器采用合成的多通道潜在信息表征作为输入,并输出所需的单通道目标模态图像作为合成图像;在所述解码器中首先使用两个连续转置卷积依次恢复图像的尺寸大小,然后通过两个卷积块获得相匹配的输出,最后,通过1×1卷积层转换为单通道的输出图像;其中两个卷积块将填充大小设置为2且卷积核大小设置为5。从而解决了多通道特性下输入输出通道大小不匹配问题;连续的转置卷积层可以稳定地集成潜在表示并完善目标模态的结构特征。
进一步的是,将所述合成图像的潜在信息表征通过解码器进行解码,获得相应的合目标模态作为合成图像;
在解码器中根据合成图像的潜在信息表征通过解码获得合成FLAIR目标模态,表示为:
Figure BDA0002273583980000032
其中,D*(·|η)为具有可学习参数η的解码器函数。
进一步的是,设定所述鉴别网络的输入大小与生成网络的合成图像相同;所述鉴别网络包括5个卷积层,前四个卷积层的步长均为2且采用都4×4卷积核;最后一个卷积层包含sigmoid激活函数,以确定输入是真实图像还是合成图像。
进一步的是,在所述鉴别网络中将真实图像与生成的图像区分开,通过目标函数进行区分,所述目标函数为:
Figure BDA0002273583980000041
其中,鉴别网络中将真实图像与生成的图像进行区别,表示为:
Figure BDA0002273583980000042
其中,X1为T1输入模态,X2为T2输入模态,Y为真实的目标模态;λ1和λ2为权重因子;上式中D是鉴别网络,G是生成网络,E是输入输出的期望值,
Figure BDA0002273583980000045
表示鉴别网络的损失函数;
采用正则化度量标准并使用L1反馈生成器,选择L1损失以减少图像的模糊度,其表示如下:
Figure BDA0002273583980000043
为了处理L1损失函数固有的模糊预测,将嵌入图像生成网络训练中的梯度差损失函数:
Figure BDA0002273583980000044
其中
Figure BDA0002273583980000046
是网络合成的图像,下标x和y表示梯度分别沿横坐标和纵坐标的移动方向;通过损失函数使合成图像和真实图像之间的梯度大小最小化;使解码值维持较大的梯度区域,以有效补偿L1反馈生成器。
本发明通过将L1反馈生成器损失和图像梯度差损失(GDL)共同作为优化LR-cGAN模型的目标函数,以确保合成图像不会严重偏离真实图像。
采用本技术方案的有益效果:
本发明利用来自不同MRI模态的集体信息,由N个端到一个端的多对一多模态MRI合成网络(称为LR-cGAN模型)包含了生成网络和鉴别网络。所提出的多模态图像合成网络是通过基于GAN模型从多个MRI模态中提取潜在信息表征(LR)来执行的;它的生成网络使用N个编码器独立提取N种不同模式的固有潜在表征;然后通过采用残差结构将潜在表征集成到潜在空间处理网络中,利用解码器生成目标模态;最后,使用鉴别网络来区分真实图像和合成图像。能够灵活的接收多个输入模态,并将多输入模态进行合成,能够有效避免信息丢失,有效提高合成图像的保真度,获得真实体现检测部位的高品质图像。适用范围宽,且实际运用效果好。
在本发明中通过GAN网络的添加使得能够呈现合成图像的高频信息,保证合成图像的真实完整。本发明从多种不同MRI模态中生成高质量的合成图像,能够提高GAN在多模态合成时的效率,提高合成结果的准确和真实性,合成图像能够真实有效体现检测部位的形态。
本发明无需使用最大化或平均化来自不同模态的潜在表示,而是直接将它们连接起来,并通过残差网络通过潜在空间处理对其进行融合,从而有效防止信息丢失,提高图像的保真度。
附图说明
图1为本发明的基于潜在信息表征GAN的多模态MRI合成方法流程示意图;
图2为本发明的实施例中LR-cGAN模型的结构示意图;
图3为本发明的实施例中生成网络的结构示意图;
图4为本发明的实施例中生成T1c模态的多模式输入模型与单模式输入模型的合成图像结果对比图;
图5为本发明的实施例中生成FLAIR模态的多模式输入模型与单模式输入模型的合成图像结果对比图;
图6为本发明的实施例中验证模型中关键组成部分的合成图像结果对比图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面结合附图对本发明作进一步阐述。
在本实施例中,参见图1-图3所示,本发明提出了基于潜在信息表征GAN的多模态MRI合成方法,包括步骤:
S100,将不同的MRI模态的集体信息输入生成网络;
S200,生成网络通过使用不同的编码器分别提取多种MRI模态的潜在信息表征;提取的潜在信息表征进一步传输到潜在空间处理网络进行集成处理;由解码器获得相应的合成目标模态,作为合成图像;
S300,将合成图像和真实图像共同输入至鉴别网络;
S400,由鉴别网络将真实图像与生成的图像区分开。
作为上述实施例的优化方案,如图2和图3所示,为了从多个模态中提取特定信息,在生成网络中包括有多个相互独立的编码器,每个编码器对应输入一个MRI模态并提取该模态的潜在信息表征;编码器的数量取决于输入模态的数量。编码器在保留模态特定信息的潜在信息表征的同时保证了模态不变性;与使用相同的卷积核来提取不同模态的特征相比,使用特定的编码器分别提取每个模态的潜在表示可以提取更有效的多种模态的潜在信息表征进行合成。
所述编码器包括3个卷积块,每个卷积块结构相同,所述卷积块依次进行填充、归一化、激活和卷积处理;
在第一个卷积块中采用3×3的零填充,进行填充;在其余2个卷积块中应用1×1的零填充,进行填充;
在所述归一化处理过程中对数据进行实例归一化处理,即对单个输入图像的全局信息进行归一化调整;此方式能够有效消除由于混洗导致的平均值归一化和批量归一化的标准偏差,并减少了噪声的引入;所述激活函数采用LeakyReLU激活函数,使网络更易于训练;将所述卷积层的步幅大小设置为2,以将特征图的尺寸减半,这样的处理方式不仅可以压缩图像大小,还可以减少最大池化层导致的详细信息丢失。
在所述生成网络中,将编码器所获得的潜在信息表征表示为LR(·)=E*(·|θ);
其中,E*(·|θ)为具有可学习参数θ的编码器函数;以独立地获取输入模式下的潜在信息表征。
作为上述实施例的优化方案,如图2和图3所示,所述潜在空间处理网络通过残差网络将每个编码器生成的潜在信息表征集成到单个潜在信息表征中:首先,直接连接不同模态下的潜在信息表征;然后,利用残差块的特性,通过特征映射的方式整合所述潜在信息表征,构成合成图像的潜在信息表征;借助残差网络,能够在不丢失信息的情况下成功地在各个卷积层中传输此类潜在表示,并最终完成多个潜在表示的融合;能够实现灵活地接收多种输入模态并保留其所有潜在信息表征;
所述残差块包括四个残差网络,每个所述残差网络包括镜像填充、批量归一化、LeakyReLU激活和卷积处理。
所述合成图像的潜在信息表征,表示为:
Figure BDA0002273583980000071
其中,R*(·|ψ)为潜在空间处理中具有可学习参数ψ的残差集成函数;LR(T1)为T1输入模态的潜在信息表征,LR(T2)为T2输入模态的潜在信息表征,n为输入模态的个数,FLAIR为目标模态。
将所述合成图像的潜在信息表征通过解码器进行解码,获得相应的合成目标模态作为合成图像;
解码器采用合成的多通道潜在信息表征作为输入,并输出所需的单通道目标模态图像作为合成图像;在所述解码器中首先使用两个连续转置卷积依次恢复图像的尺寸大小,然后通过两个卷积块获得相匹配的输出,最后,通过1×1卷积层转换为单通道的输出图像;其中两个卷积块将填充大小设置为2且卷积核大小设置为5。从而解决了多通道特性下输入输出通道大小不匹配问题;连续的转置卷积层可以稳定地集成潜在表示并完善目标模态的结构特征。
将所述合成图像的潜在信息表征通过解码器进行解码,获得相应的合目标模态作为合成图像;
在解码器中根据合成图像的潜在信息表征通过解码获得合成FLAIR目标模态,表示为:
Figure BDA0002273583980000072
其中,D*(·|η)为具有可学习参数η的解码器函数。
作为上述实施例的优化方案,如图2和图3所示,设定所述鉴别网络的输入大小与生成网络的合成图像相同;所述鉴别网络包括5个卷积层,前四个卷积层的步长均为2且采用都4×4卷积核;最后一个卷积层包含sigmoid激活函数,以确定输入是真实图像还是合成图像。
在所述鉴别网络中将真实图像与生成的图像区分开,通过目标函数进行区分,所述目标函数为:
Figure BDA0002273583980000081
其中,鉴别网络中将真实图像与生成的图像进行区别,表示为:
Figure BDA0002273583980000082
其中,X1为T1输入模态,X2为T2输入模态,Y为真实的目标模态;λ1和λ2为权重因子;上式中D是鉴别网络,G是生成网络,E是输入输出的期望值,
Figure BDA0002273583980000085
表示鉴别网络的损失函数;
采用正则化度量标准并使用L1反馈生成器,选择L1损失以减少图像的模糊度,其表示如下:
Figure BDA0002273583980000083
为了处理L1损失函数固有的模糊预测,将嵌入图像生成网络训练中的梯度差损失函数:
Figure BDA0002273583980000084
其中
Figure BDA0002273583980000086
是网络合成的图像,下标x和y表示梯度分别沿横坐标和纵坐标的移动方向;通过损失函数使合成图像和真实图像之间的梯度大小最小化;使解码值维持较大的梯度区域,以有效补偿L1反馈生成器。
本发明通过将L1反馈生成器损失和图像梯度差损失(GDL)共同作为优化本发明提出LR-cGAN模型的目标函数,以确保合成图像不会严重偏离真实图像。
为了评估多模式输入模型与单模式输入模型相比的影响,并证明我们的模型可以灵活地接收多个输入,我们通过使用不同的模式输入来比较合成图像。具体来说,我们分别使用T2,T1+T2,T1+T2+FLAIR作为生成T1c模态的输入,并使用T2,T1+T2,T1+T2+T1c作为生成FLAIR模态的输入。表1、图4和图5分别定量和定性比较了实验结果,验证本发明所提出的合成方法。
表1中表示了所提出网络模型在不同输入下图像合成的性能:
表1
Figure BDA0002273583980000091
首先,观察到T1c模态的综合结果。如表1所示,由T1+T2+FLAIR合成的T1c的平均PSNR高于仅由T1+T2和T2合成的PSNR。与使用单一模式的T2相比,采用两种模式(T1+T2)可获得更好的结果,PSRN从27.73改善到29.36。这是因为T1模态包含丰富的解剖结构信息,可以更好地合成T1c模态。当将FLAIR模态(T1+T2+FLAIR)合并时,该模型会进一步产生改进的结果,PSNR从29.36增加到30.77。表1中的NRMSE和SSIM值表明相同的结论。可视化结果如图4所示,由三种模态产生的合成图像提供了最佳的图像质量,并同时保留了图像对比度和详细的组织(如箭头所示)。此外,如果仅使用T2方式进行合成,则方框中的病变将大大丢失;尽管通过额外使用T1可以略微改善合成图像的质量,但仍比采用三种模态进行合成的图像质量差,三模态输入模型的差异最小。
然后,分析了FLAIR模态的综合结果。为了进行定量比较,如表1所示,三模态输入模型的合成图像实现了最高的PSNR和SSIM,以及最低的NRMSE。可视化结果如图5所示,,依次为使用T2,T1+T2,T1+T2+T1c作为生成FLAIR模态的输入下的合成结果图,无论从对比度还是图像纹理来看,三模态输入模型的综合结果都显着改善了FLAIR图像的细节(如箭头所示),与真实的FLAIR图像极为相似,三模态输入模型的合成图像与真实图像之间的差异显然最小。
根据以上定性和定量结果,我们的模型不仅可以灵活地接收各种输入,而且可以融合所有输入信息以合成比单模态模型更高质量的图像。
为了研究该方法的关键组成部分的贡献,我们使用T1+T2作为输入来生成FLAIR,以评估三个关键组成部分:对抗网络(GAN),图像梯度差损失(GDL)和潜在空间处理网络(LSPN)。表2和图6分别定量和定性比较了实验结果,验证本发明所提出的合成方法。
表2
Figure BDA0002273583980000101
为了评估本发明提出的LR-cGAN模型的对抗网络的贡献,我们将提出的模型与去除了鉴别器的模型进行了比较。在表2中显示了有关PSNR,SSIM和NRMSE的详细定量比较。如观察到的那样,在具有对抗网络的模型中,尽管SSIM略有降低,但PSNR从28.01升高到28.23,NRMSE从0.178降低到0.170。从这些定量结果中,我们可以清楚地看到对抗训练方法的使用有助于提高合成图像的质量。图像结果如图6所示,合成图像类似于真实图像,因为没有错误地添加任何额外的结构,但是图像损失了可见的高频信息,使整个图像缺乏精细的结构信息;没有对抗网络的模型的误差几乎均匀地大于完整提出的模型。换句话说,对抗网络系统地减少了误差并提供了精细的结构信息。
通过图像梯度差损失来收集图像的边缘信息并改善图像的清晰度。为了评估梯度差损失函数的效果,将其从提出的模型中删除,并保留其他网络模块。定量结果总结在表2中。与没有GDL损失函数的模型相比,本发明模型的PSNR从27.64增加到28.23,并且SSIM和NRMSE也呈现接近良好方向的趋势。这说明所提出的模型明显优于没有GDL损失函数的模型。图6中示出了所提出的方法的合成图像,从没有GDL损失函数的合成图像中,可以看出部分灰质不是很正确;无GDL损失的模型与完整模型之间的综合质量差异。因此,GDL损失的增加不仅纠正了模型的某些错误纹理合成,而且使图像的对比度更接近于真实情况。
为了评估模型的潜在空间处理网络(LSPN)的效果,在提取潜在表示之后,删除潜在空间处理网络并直接对其进行解码以生成合成图像。结果列在表2中,从中我们可以看到LSPN网络在所有三个指标中均显着提高了合成图像的质量,并且NRMSE的改进明显。合成的结果在图6中,合成的图像呈现一定程度的对比度失真,即,图像总体上比真实图像明亮,而暗区比真实图像明显暗。因此,LSPN是整合从不同模态中提取的特征的关键步骤,并且对所提出模型的性能做出了重要贡献。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (10)

1.基于潜在信息表征GAN的多模态MRI合成方法,其特征在于,包括步骤:
S100,将不同的MRI模态的集体信息输入生成网络;
S200,生成网络通过使用不同的编码器分别提取多种MRI模态的潜在信息表征;提取的潜在信息表征进一步传输到潜在空间处理网络进行集成处理;由解码器获得相应的合成目标模态,作为合成图像;
S300,将合成图像和真实图像共同输入至鉴别网络;
S400,由鉴别网络将真实图像与生成的图像区分开。
2.根据权利要求1所述的基于潜在信息表征GAN的多模态MRI合成方法,其特征在于,在生成网络中包括有多个相互独立的编码器,每个编码器对应输入一个MRI模态并提取该模态的潜在信息表征;编码器的数量取决于输入模态的数量。
3.根据权利要求2所述的基于潜在信息表征GAN的多模态MRI合成方法,其特征在于,所述编码器包括3个卷积块,每个卷积块结构相同,所述卷积块依次进行填充、归一化、激活和卷积处理;
在第一个卷积块中采用3×3的零填充,进行填充;在其余2个卷积块中应用1×1的零填充,进行填充;
在所述归一化处理过程中对数据进行实例归一化处理,即对单个输入图像的全局信息进行归一化调整;所述激活函数采用LeakyReLU激活函数;将所述卷积层的步幅大小设置为2,以将特征图的尺寸减半。
4.根据权利要求3所述的基于潜在信息表征GAN的多模态MRI合成方法,其特征在于,在所述生成网络中,将编码器所获得的潜在信息表征表示为LR(·)=E*(·|θ);
其中,E*(·|θ)为具有可学习参数θ的编码器函数;以独立地获取输入模式下的潜在信息表征。
5.根据权利要求4所述的基于潜在信息表征GAN的多模态MRI合成方法,其特征在于,所述潜在空间处理网络通过残差网络将每个编码器生成的潜在信息表征集成到单个潜在信息表征中:首先,直接连接不同模态下的潜在信息表征;然后,利用残差块的特性,通过特征映射的方式整合所述潜在信息表征,构成合成图像的潜在信息表征;
所述残差块包括四个残差网络,每个所述残差网络包括镜像填充、批量归一化、LeakyReLU激活和卷积处理。
6.根据权利要求5所述的基于潜在信息表征GAN的多模态MRI合成方法,其特征在于,所述合成图像的潜在信息表征,表示为:
Figure FDA0002273583970000021
其中,R*(·|ψ)为潜在空间处理中具有可学习参数ψ的残差集成函数;LR(T1)为T1输入模态的潜在信息表征,LR(T2)为T2输入模态的潜在信息表征,n为输入模态的个数,FLAIR为目标模态。
7.根据权利要求6所述的基于潜在信息表征GAN的多模态MRI合成方法,其特征在于,将所述合成图像的潜在信息表征通过解码器进行解码,获得相应的合成目标模态作为合成图像;
解码器采用合成的多通道潜在信息表征作为输入,并输出所需的单通道目标模态图像作为合成图像;在所述解码器中首先使用两个连续转置卷积依次恢复图像的尺寸大小,然后通过两个卷积块获得相匹配的输出,最后,通过1×1卷积层转换为单通道的输出图像;其中两个卷积块将填充大小设置为2且卷积核大小设置为5。
8.根据权利要求7所述的基于潜在信息表征GAN的多模态MRI合成方法,其特征在于,将所述合成图像的潜在信息表征通过解码器进行解码,获得相应的合目标模态作为合成图像;
在解码器中根据合成图像的潜在信息表征通过解码获得合成FLAIR目标模态,表示为:
Figure FDA0002273583970000022
其中,D*(·|η)为具有可学习参数η的解码器函数。
9.根据权利要求8所述的基于潜在信息表征GAN的多模态MRI合成方法,其特征在于,设定所述鉴别网络的输入大小与生成网络的合成图像相同;所述鉴别网络包括5个卷积层,前四个卷积层的步长均为2且采用都4×4卷积核;最后一个卷积层包含sigmoid激活函数,以确定输入是真实图像还是合成图像。
10.根据权利要求9所述的基于潜在信息表征GAN的多模态MRI合成方法,其特征在于,在所述鉴别网络中将真实图像与生成的图像区分开,通过目标函数进行区分,所述目标函数为:
Figure FDA0002273583970000031
其中,鉴别网络中将真实图像与生成的图像进行区别,表示为:
Figure FDA0002273583970000032
其中,X1为T1输入模态,X2为T2输入模态,Y为真实的目标模态;λ1和λ2为权重因子;上式中D是鉴别网络,G是生成网络,E是输入输出的期望值,
Figure FDA0002273583970000033
表示鉴别网络的损失函数;
采用正则化度量标准并使用L1反馈生成器,选择L1损失以减少图像的模糊度,其表示如下:
Figure FDA0002273583970000034
为了处理L1损失函数固有的模糊预测,将嵌入图像生成网络训练中的梯度差损失函数:
Figure FDA0002273583970000035
其中
Figure FDA0002273583970000036
是网络合成的图像,下标x和y表示梯度分别沿横坐标和纵坐标的移动方向;通过损失函数使合成图像和真实图像之间的梯度大小最小化;使解码值维持较大的梯度区域,以有效补偿L1反馈生成器。
CN201911114218.9A 2019-11-14 2019-11-14 基于潜在信息表征gan的多模态mri合成方法 Expired - Fee Related CN110866888B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911114218.9A CN110866888B (zh) 2019-11-14 2019-11-14 基于潜在信息表征gan的多模态mri合成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911114218.9A CN110866888B (zh) 2019-11-14 2019-11-14 基于潜在信息表征gan的多模态mri合成方法

Publications (2)

Publication Number Publication Date
CN110866888A true CN110866888A (zh) 2020-03-06
CN110866888B CN110866888B (zh) 2022-04-26

Family

ID=69654198

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911114218.9A Expired - Fee Related CN110866888B (zh) 2019-11-14 2019-11-14 基于潜在信息表征gan的多模态mri合成方法

Country Status (1)

Country Link
CN (1) CN110866888B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113012086A (zh) * 2021-03-22 2021-06-22 上海应用技术大学 一种跨模态图像的合成方法
CN113674383A (zh) * 2020-05-15 2021-11-19 华为技术有限公司 生成文本图像的方法及装置

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009045579A2 (en) * 2007-06-14 2009-04-09 The Regents Of The University Of California Multimodal imaging probes for in vivo targeted and non-targeted imaging and therapeutics
CN109035356A (zh) * 2018-07-05 2018-12-18 四川大学 一种基于pet图形成像的系统及方法
CN109472817A (zh) * 2018-09-27 2019-03-15 浙江工业大学 一种基于循环生成对抗网络的多序列磁共振图像配准方法
CN109741410A (zh) * 2018-12-07 2019-05-10 天津大学 基于深度学习的荧光编码微球图像生成及标注方法
CN110097512A (zh) * 2019-04-16 2019-08-06 四川大学 基于Wasserstein生成对抗网络的三维MRI图像去噪模型的构建方法及应用
CN110163897A (zh) * 2019-04-24 2019-08-23 艾瑞迈迪科技石家庄有限公司 一种基于合成超声图像的多模态图像配准方法
CN110276736A (zh) * 2019-04-01 2019-09-24 厦门大学 一种基于权值预测网络的磁共振图像融合方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009045579A2 (en) * 2007-06-14 2009-04-09 The Regents Of The University Of California Multimodal imaging probes for in vivo targeted and non-targeted imaging and therapeutics
CN109035356A (zh) * 2018-07-05 2018-12-18 四川大学 一种基于pet图形成像的系统及方法
CN109472817A (zh) * 2018-09-27 2019-03-15 浙江工业大学 一种基于循环生成对抗网络的多序列磁共振图像配准方法
CN109741410A (zh) * 2018-12-07 2019-05-10 天津大学 基于深度学习的荧光编码微球图像生成及标注方法
CN110276736A (zh) * 2019-04-01 2019-09-24 厦门大学 一种基于权值预测网络的磁共振图像融合方法
CN110097512A (zh) * 2019-04-16 2019-08-06 四川大学 基于Wasserstein生成对抗网络的三维MRI图像去噪模型的构建方法及应用
CN110163897A (zh) * 2019-04-24 2019-08-23 艾瑞迈迪科技石家庄有限公司 一种基于合成超声图像的多模态图像配准方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
C. ZHANG: "MS-GAN: GAN-Based Semantic Segmentation of Multiple Sclerosis Lesions in Brain Magnetic Resonance Imaging", 《2018 DIGITAL IMAGE COMPUTING: TECHNIQUES AND 》 *
陈润: "基于人工免疫的疾病预测诊断模型研究", 《数理医药学杂志》 *
陈锟: "生成对抗网络在医学图像处理中的应用", 《生命科学仪器》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113674383A (zh) * 2020-05-15 2021-11-19 华为技术有限公司 生成文本图像的方法及装置
CN113012086A (zh) * 2021-03-22 2021-06-22 上海应用技术大学 一种跨模态图像的合成方法
CN113012086B (zh) * 2021-03-22 2024-04-16 上海应用技术大学 一种跨模态图像的合成方法

Also Published As

Publication number Publication date
CN110866888B (zh) 2022-04-26

Similar Documents

Publication Publication Date Title
CN112507997B (zh) 一种基于多尺度卷积和感受野特征融合的人脸超分辨系统
CN110047138A (zh) 一种磁共振薄层图像重建方法
CN110490082B (zh) 一种有效融合神经网络特征的道路场景语义分割方法
CN112541864A (zh) 一种基于多尺度生成式对抗网络模型的图像修复方法
Tuzel et al. Global-local face upsampling network
KR102359474B1 (ko) 뉴럴 네트워크를 이용한 결측 영상 데이터 대체 방법 및 그 장치
CN111080567A (zh) 基于多尺度动态卷积神经网络的遥感图像融合方法及系统
Du et al. Accelerated super-resolution MR image reconstruction via a 3D densely connected deep convolutional neural network
CN109389585A (zh) 一种基于全卷积神经网络的脑组织提取方法
Upadhyay et al. Robust super-resolution GAN, with manifold-based and perception loss
CN113255571B (zh) 抗jpeg压缩的伪造图像检测方法
CN110866888B (zh) 基于潜在信息表征gan的多模态mri合成方法
CN113781517A (zh) 用于运动估计的系统和方法
CN115375711A (zh) 基于多尺度融合的全局上下文关注网络的图像分割方法
CN117333750A (zh) 空间配准与局部全局多尺度的多模态医学图像融合方法
CN116739899A (zh) 基于saugan网络的图像超分辨率重建方法
CN117079105A (zh) 遥感图像空谱融合方法、装置、电子设备及存储介质
CN114943656A (zh) 一种人脸图像修复方法及系统
CN116664446A (zh) 基于残差密集块的轻量级暗光图像增强方法
CN117151990A (zh) 一种基于自注意力编码解码的图像去雾方法
Wang et al. Brain MR image super-resolution using 3D feature attention network
Sander et al. Autoencoding low-resolution MRI for semantically smooth interpolation of anisotropic MRI
CN112785540B (zh) 一种扩散加权图像的生成系统及方法
CN114663310A (zh) 一种基于多注意力融合的超声图像去噪方法
CN116862765A (zh) 一种医学影像超分辨率重建方法及系统

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220426