CN115937038A - 一种基于wgan-gp模型的网格形变数据增强方法 - Google Patents

一种基于wgan-gp模型的网格形变数据增强方法 Download PDF

Info

Publication number
CN115937038A
CN115937038A CN202211682574.2A CN202211682574A CN115937038A CN 115937038 A CN115937038 A CN 115937038A CN 202211682574 A CN202211682574 A CN 202211682574A CN 115937038 A CN115937038 A CN 115937038A
Authority
CN
China
Prior art keywords
deformation
net
train
model
layer
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.)
Pending
Application number
CN202211682574.2A
Other languages
English (en)
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.)
University of Shanghai for Science and Technology
Original Assignee
University of Shanghai for Science and Technology
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 University of Shanghai for Science and Technology filed Critical University of Shanghai for Science and Technology
Priority to CN202211682574.2A priority Critical patent/CN115937038A/zh
Publication of CN115937038A publication Critical patent/CN115937038A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Abstract

本发明公开了一种基于WGAN‑GP模型的网格形变数据增强方法,涉及图像处理、数据增强技术领域,包括:构建训练数据集和测试数据集;构建WGAN‑GP模型;基于训练数据集训练WGAN‑GP模型;使用测试数据集评估模型性能,确定网格形变数据增强模型;将随机噪声分别与医疗内窥镜图像作为输入,输出对应的形变网格,将医疗内窥镜图像分别和形变网格进行warp操作得到数据增强后的医疗内窥镜图像。本发明构建WGAN‑GP模型自动生成形变网格,通过形变网格对医疗内窥镜图像进行形变增强,增强后的医疗内窥镜图像真实且具有一定多样性,有效解决数据样本稀缺和数据样本不均衡的问题,进而提高基于医疗内窥镜图像的人工智能算法模型精度及泛化性。

Description

一种基于WGAN-GP模型的网格形变数据增强方法
技术领域
本发明涉及图像处理、数据增强技术领域,更具体的说是涉及一种基于WGAN-GP模型的网格形变数据增强方法。
背景技术
近年来,AI技术在医学影像分析、医疗决策等场景领域得到广泛应用,驱动传统医疗研究迈进了全新的阶段。作为“AI+医学”的一个重要研究方向,AI驱动的医学影像分析技术日趋成熟。然而,由于医学影像数据比较敏感和隐私,且高质量的数据标注工作非常困难,大部分医学影像数据集规模都较小,这也限制了AI模型在医学影像分析领域的应用。数据增强技术通过对已有医学影像数据(尤其是数量有限的类别)进行数据扩充,使得数据样本更加丰富、类别更加均衡,有效提升AI模型(尤其是深度神经网络模型)在医学影像分析领域中的任务性能。
目前在医疗内窥镜图像数据增强技术领域中,平移、翻转、旋转、缩放、裁剪、仿射变换等几何变换技术常被用于对已有的医疗内窥镜图像进行数据增强;这类传统方法在一定程度上缓解了数据样本不足的问题,但这类方法依赖于人为预设参数,且无法针对医疗内窥镜图像进行专门设计,导致增强后的图像数据真实性较差,进而导致增强后的数据对后续的病灶诊断、病灶分割等任务精度提升有限。近年来,生成对抗网络(GAN)因其强大的图像生成能力,逐渐被应用在医学影像数据增强领域,GAN通过输入随机噪声直接生成逼真的医学影像,以达到扩充医学影像数据集的目的;然而通过GAN模型直接生成医学影像,这类模型的训练需要大量的数据基础,医疗内窥镜图像数据这类医学影像数据因其敏感性和隐私性,数据集规模往往较小,所以这类方法并不适用。
因此,提出一种数据样本量需求小,且能够自动生成形变网格对医疗内窥镜图像进行形变增强以达到扩充数据集的目的的网格形变数据增强方法,使增强后的医疗内窥镜图像真实且具有一定多样性是本领域技术人员亟需解决的问题。
发明内容
有鉴于此,本发明提供了一种基于WGAN-GP模型的网格形变数据增强方法,该方法针对医疗内窥镜图像的特点,通过构建WGAN-GP模型自动生成形变网格,进而通过形变网格对医疗内窥镜图像进行形变增强以达到扩充数据集的目的,增强后的医疗内窥镜图像真实且具有一定多样性,为了实现上述目的,本发明采用如下技术方案:
本发明通过WGAN-GP模型直接生成形变网格,进而通过形变网格与医疗内窥镜图像进行warp操作后得到数据增强后的医疗内窥镜图像,并非通过模型直接生成图像,降低了模型训练过程的数据样本量需求;同时在模型训练的损失函数中设计了平滑形变项和病变区域LA大形变项,其中平滑形变项约束模型产生平滑形变,保证了形变的真实性和合理性,病变区域LA大形变项鼓励病变区域产生较大形变,保证了病变区域形变的显著性和多样性。
一种基于WGAN-GP模型的网格形变数据增强方法,包括以下步骤:
S1:构建训练数据集和测试数据集,将归一化后的医疗内窥镜图像数据按照7:3的比例划分为训练集Xtrain和测试集Xtest,将训练集图像及其相对应的病变区域LA标注打包成训练数据集;
S2:构建WGAN-GP模型,所述WGAN-GP模型包括基于残差结构的生成器网络NetG和基于残差结构的判别器网络NetD
S3:基于训练数据集迭代训练WGAN-GP模型,设置WGAN-GP模型中判别器网络NetD和生成器网络NetG的损失函数,分别通过损失函数交替迭代优化判别器网络NetD和生成器网络NetG,直至模型收敛,保存模型参数;
S4:使用医疗内窥镜图像测试集Xtest评估训练完成后的生成器网络NetG,评估性能符合需求的生成器网络NetG即为网格形变数据增强模型modelaug
S5:采样满足标准正态分布的随机噪声{z1,z2,…,zk},分别将随机噪声{z1,z2,…,zk}与医疗内窥镜图像X作为网格形变数据增强模型modelaug的模型输入,输出对应的形变网格{g1,g2,…,gk},将医疗内窥镜图像X分别和形变网格{g1,g2,…,gk}进行warp操作后得到数据增强后的医疗内窥镜图像{Xaug_1,Xaug_2,…,Xaug_k}。
可选的,所述步骤S2中所述生成器网络NetG为基于残差结构的卷积神经网络,包含第一残差结构块、第二残差结构块、上采样块、卷积层、批归一化层、最大值池化层、线性层、reshape层以及Tanh激活函数,其中第一残差结构块和第二残差结构块由卷积层、批归一化层、LeakyReLU激活函数交替串联组成,并通过残差连接将低层输出与高层输出进行连接;上采样块由卷积层、批归一化层、LeakyReLU激活函数、上采样层交替串联组成。
可选的,生成器网络NetG的输入为医疗内窥镜图像训练集Xtrain和随机噪声z,生成器网络NetG的输出为形变网格g。
可选的,所述步骤S2中所述判别器网络NetD为基于残差结构的卷积神经网络,包含第三残差结构块、第四残差结构块、卷积层、层归一化层、最大值池化层、自适应平均池化层、reshape层以及线性层,其中第三残差结构块和第四残差结构块由卷积层、层归一化层、ReLU激活函数交替串联组成,并通过残差连接将低层输出与高层输出进行连接。
可选的,判别器网络NetD的输入为医疗内窥镜图像训练集Xtrain和形变图像Xtrain_aug,判别器网络NetD的输出为对应的值NetD(Xtrain)和NetD(Xtrain_aug)。
可选的,所述步骤S3具体包含以下步骤:
S3-1:初始化生成器网络NetG和判别器网络NetD
S3-2:将医疗内窥镜图像训练集Xtrain和随机噪声z作为NetG的输入,输出形变网格g;
S3-3:将医疗内窥镜图像训练集Xtrain和形变网格g进行warp操作得到形变图像Xtrain_aug,将医疗内窥镜图像训练集Xtrain和形变图像Xtrain_aug输入判别器网络NetD,输出对应的值NetD(Xtrain)和NetD(Xtrain_aug);
S3-4:构造判别器网络NetD的损失函数为LD=Lm1×LGrad_penalty,其中LD1=-NetD(Xtrain)+NetD(Xtrain_aug),LGrad_penalty为梯度惩罚项,λ1为权重参数;
S3—5:构造生成器网络NetG的损失函数为LG=LG12×Lsmooth3×LLA_deform,其中LG1=-NetD(Xtrain_aug),Lsmooth为平滑形变项,LLA_deform为病变区域LA大形变项,λ2和λ3为权重参数;
S3-6:构造训练过程中监测指标项Lmonitor=NetD(Xtrain)-NetD(Xtrain_aug),用于观察WGAN-GP模型的训练趋势;
S3-7:分别通过损失函数LD和LG交替迭代优化判别器网络NetD和生成器网络NetG,观察监测指标项Lmonitor的收敛情况,当监测指标项Lmonitor曲线收敛时停止模型训练,保存训练后的判别器网络NetD和生成器网络NetG的模型参数。
可选的,所述步骤S3-5中所述平滑形变项Lsmooth为约束生成器网络NetG产生平滑形变的损失函数项,
Figure BDA0004019066910000041
Figure BDA0004019066910000042
其中N为△g内像素点个数,p表示△g内任意一点,px和py分别为p在X方向和Y方向的位置,△g=g-g0,g0为原图-原图的形变网格。
可选的,所述步骤S3-5中所述病变区域LA大形变项LLA_deform为鼓励病变区域LA产生较大形变的损失函数项,
Figure BDA0004019066910000051
Figure BDA0004019066910000052
其中N为△g内像素点个数,R为△g内属于病变区域LA位置的任意一点,R′为△g内不属于病变区域LA位置的任意一点,a为常数(a>1),△g=g-g0,g0为原图-原图的形变网格。
经由上述的技术方案可知,与现有技术相比,本发明公开提供了一种基于WGAN-GP模型的网格形变数据增强方法,具有以下优点:
1.本发明针对医疗内窥镜图像实际采集过程中病人的呼吸或者移动会导致器官组织产生局部位移或变形等特点,该方法通过构建WGAN-GP模型自动生成形变网格,进而通过形变网格对医疗内窥镜图像进行形变增强,增强后的医疗内窥镜图像真实且具有一定多样性。
2.本发明方法通过WGAN-GP模型直接生成形变网格,进而通过形变网格对医疗内窥镜图像进行形变增强,并非通过模型直接生成图像,降低了模型训练过程的数据样本量需求。
3.本发明方法在模型训练的损失函数中设计了平滑形变项和病变区域LA大形变项,其中平滑形变项约束模型产生平滑形变,保证了形变的真实性和合理性,病变区域LA大形变项鼓励病变区域产生较大形变,保证了病变区域形变的显著性和多样性。
4.本发明方法不仅可以对医疗内窥镜图像进行数据增强,还可以对相应的病变区域标签进行同步的数据增强,有效扩充医疗内窥镜图像及其病变区域标注数据集,进而解决数据样本稀缺和数据样本不均衡的问题;利用扩充后的数据集训练人工智能算法模型(尤其是深度神经网络模型)进行病变诊断、病变分割等医学图像任务,可以有效提高最终模型精度及泛化性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为一种基于WGAN-GP模型的网格形变数据增强方法的示意图;
图2为本发明中生成器网络NetG的模型结构图;
图3为第一残差结构块的结构图;
图4为第二残差结构块的结构图;
图5为上采样块的结构图;
图6为本发明中医疗内窥镜图像训练集Xtrain和形变网格g进行warp操作的示意图;
图7为本发明中判别器网络NetD的模型结构图;
图8为第三残差结构块的结构图;
图9为第四残差结构块的结构图;
图10为一种基于WGAN-GP模型的网格形变数据增强方法在胃肠道内窥镜数据集上的应用示例。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示为本发明提供的一种基于WGAN-GP模型的网格形变数据增强方法,该方法具体包括以下步骤:
S1:构建训练数据集和测试数据集,包括将归一化后的医疗内窥镜图像数据按照7∶3的比例划分为训练集Xtrain和测试集Xtest,将训练集图像及其相对应的病变区域LA标注打包成训练数据集。
医疗内窥镜图像数据以RGB格式存储,将图像各颜色通道数值都除以255,使其数值归一化到0-1,将归一化后的医疗内窥镜图像数据按照7:3的比例划分为训练集Xtrain和测试集Xtest,医疗内窥镜图像训练集Xtrain用于训练WGAN-GP模型,医疗内窥镜图像测试集Xtest用于评估训练完成后的WGAN-GP模型。
所构建的训练数据集还包括与医疗内窥镜图像训练集Xtrain相对应的病变区域LA标注,病变区域LA标注以二值图方式存储,病变区域内像素用1表示,病变区域外像素用0表示。
S2:构建WGAN-GP模型,所述WGAN-GP模型包括基于残差结构的生成器网络NetG和基于残差结构的判别器网络NetD
生成器网络NetG的模型结构如图2所示,该网络为基于残差结构的卷积神经网络,由第一残差结构块(如图3所示)、第二残差结构块(如图4所示)、上采样块(如图5所示)、卷积层、批归一化层、最大值池化层、线性层、reshape层以及Tanh激活函数组成;第一残差结构块和第二残差结构块由卷积层、批归一化层、LeakyReLU激活函数交替串联组成,并通过残差连接将低层输出与高层输出进行连接;输入特征图经过第一残差结构块后输出特征图尺寸不变,而输入特征图经过第二残差结构块后输出特征图尺寸缩小为原来的
Figure BDA0004019066910000071
上采样块由卷积层、批归一化层、LeakyReLU激活函数、上采样层交替串联组成。
生成器网络NetG的输入为医疗内窥镜图像训练集Xtrain和随机噪声z(z分布满足标准正态分布,即均数为0、标准差为1的正态分布N(0,1)),生成器网络NetG的输出为形变网格g;医疗内窥镜图像训练集Xtrain(size=3×512×512为例)经过对应的卷积层、批归一化层、最大值池化层、第一残差结构块和第二残差结构块后,输出特征图FX的size=128×64×64;随机噪声z(以维度=100为例)经过对应的线性层、reshape层和第一残差结构块后,输出特征图Fz的size=128×64×64;将特征图FX和Fz沿通道维度进行拼接,拼接后特征图FX+z的size=256×64×64;将FX+z经过3个上采样块、卷积层和Tanh激活函数,调整维度顺序后输出size=512×512×2的形变网格g。
图6为医疗内窥镜图像训练集Xtrain和形变网格g进行warp操作的示意图,如图所示,形变网格g是大小为(H,W,2)的数组,其中H,W分别为医疗内窥镜图像训练集Xtrain长和宽方向的尺寸;以H=W=512为例,即形变网格g为size=(512,512,2)的数组。
形变网格g在(h,w)处对应的值(g1(h,w),g2(h,w))指示形变图像Xtrain_aug在(h,w)处对应的像素值由医疗内窥镜图像训练集Xtrain在(g1(h,w),g2(h,w))处对应的像素值决定,即
Figure BDA0004019066910000081
形变网格g的数值范围为[-1,1],对应于医疗内窥镜图像训练集Xtrain最左上、最左下、最右上、最右下的位置分别用(-1,-1)、(-1,1)、(1,-1)、(1,1)来表示。
形变网格g指示的坐标(g1(h,w),g2(h,w))的数值范围在[-1,1]之间,需要根据一个浮点型的坐标值在医疗内窥镜图像训练集Xtrain中对像素值进行采样,本方法采用双线性插值的方法进行填充,即医疗内窥镜图像训练集Xtrain在(g1(h,w),g2(h,w))处对应的像素值
Figure BDA0004019066910000082
由(g1(h,w),g2(h,w))周围的四个像素值采用双线性插值方法计算得到。
判别器网络NetD的模型结构如图7所示,该网络为基于残差结构的卷积神经网络,由第三残差结构块(如图8所示)、第四残差结构块(如图9所示)、卷积层、层归一化层、最大值池化层、自适应平均池化层、reshape层以及线性层组成;第三残差结构块和第四残差结构块由卷积层、层归一化层、ReLU激活函数交替串联组成,并通过残差连接将低层输出与高层输出进行连接;输入特征图经过第三残差结构块后输出特征图尺寸不变,而输入特征图经过第四残差结构块后输出特征图尺寸缩小为原来的
Figure BDA0004019066910000091
判别器网络NetD的输入为医疗内窥镜图像训练集Xtrain和形变图像Xtrain_aug,判别器网络NetD的输出为对应的值NetD(Xtrain)和NetD(Xtrain_aug);医疗内窥镜图像训练集Xtrain和形变图像Xtrain_aug(size=3×512×512为例)输入判别器网络NetD后,经过卷积层、层归一化层、最大值池化层、第三残差结构块和2个第四残差结构块后,特征图尺寸逐渐减小,通道数逐渐增加,特征图size变为256×32×32;接着经过自适应平均池化层后特征图size变为256×1×1,经reshape层降维后经过2个线性层,最终输出对应的值NetD(Xtrain)和NetD(Xtrain_aug)。
S3:基于训练数据集迭代训练WGAN-GP模型,设置WGAN-GP模型中判别器网络NetD和生成器网络NetG的损失函数,分别通过损失函数交替迭代优化判别器网络NetD和生成器网络NetG,直至模型收敛,保存模型参数。
S3-1:使用kaiming_normal参数初始化方式分别初始化生成器网络NetG和判别器网络NetD的参数。
S3-2:将医疗内窥镜图像训练集Xtrain和随机噪声z作为生成器网络NetG的输入,输出形变网格g。
S3-3:将医疗内窥镜图像训练集Xtrain和形变网格g进行warp操作得到形变图像Xtrain_aug,将医疗内窥镜图像训练集Xtrain和形变图像Xtrain_aug输入判别器网络NetD,输出对应的值NetD(Xtrain)和NetD(Xtrain_aug)。
S3-4:构造判别器网络NetD的损失函数为LD=LD11×LGrad_penalty,其中LD1=-NetD(Xtrain)+NetD(Xtrain_aug),LGrad_penalty为梯度惩罚项,λ1为权重参数。
WGAN-GP模型的迭代优化过程是一个极大极小博弈问题,即寻找生成器网络NetG和判别器网络NetD之间的平衡点;生成器网络NetG的目标是使输出的形变网格g与医疗内窥镜图像训练集Xtrain进行warp操作后得到的形变图像Xtrain_aug的分布尽可能接近于医疗内窥镜图像训练集Xtrain的分布;而判别器网络NetD的目标则是尽可能区分开医疗内窥镜图像训练集Xtrain和形变图像Xtrain_aug;不断交替迭代训练生成器网络NetG和判别器网络NetD,最终达到平衡点。
在WGAN-GP模型的训练过程中,使用W=NetD(Xtrain)-NetD(Xtrain_aug)来近似表示医疗内窥镜图像训练集Xtrain分布与形变图像Xtrain_aug分布之间的Wasserstein距离,即W越小,形变图像Xtrain_aug的分布越接近于医疗内窥镜图像训练集Xtrain的分布。
判别器网络NetD的目标是尽可能区分开医疗内窥镜图像训练集Xtrain和形变图像Xtrain_aug,即使得W=NetD(Xtrain)-NetD(Xtrain_aug)尽可能大,则判别器网络NetD损失函数的第一部分可表示为LD1=-NetD(Xtrain)+NetD(Xtrain_aug)。
使用W=NetD(Xtrain)-NetD(Xtrain_aug)来近似表示医疗内窥镜图像训练集Xtrain和形变图像Xtrain_aug分布之间的Wasserstein距离,需要满足Lipschitz连续条件,因此在判别器网络NetD的损失函数中引入梯度惩罚项
Figure BDA0004019066910000101
Figure BDA0004019066910000102
其中
Figure BDA0004019066910000103
随机数β∈[0,1)。
因此,将判别器网络NetD的损失函数定义为:LD=LD11×LGrad_penalty,其中λ1为权重参数。
S3-5:构造生成器网络NetG的损失函数为LG=LG12×Lsmooth3×LLA_deform,其中LG1=-NetD(Xtrain_aug),Lsmooth为平滑形变项,LLA_deform为病变区域LA大形变项,λ2和λ3为权重参数。
生成器网络NetG的目标是使形变图像Xtrain_aug的分布尽可能接近于医疗内窥镜图像训练集Xtrain的分布,即使得W=NetD(Xtrain)-NetD(Xtrain_aug)尽可能小,因为NetD(Xtrain)这一项与生成器网络NetG无关,所以生成器网络NetG损失函数的第一部分可表示为LG1=-NetD(Xtrain_aug)。
为了使得生成器网络NetG输出的形变网格g所对应的形变更加真实且平滑,在生成器网络NetG的损失函数中引入了平滑形变项以约束生成器网络NetG产生平滑形变,
Figure BDA0004019066910000111
Figure BDA0004019066910000112
其中N为△g内像素点个数,p表示△g内任意一点,px和py分别为p在X方向和Y方向的位置,△g=g-g0,g0为原图-原图的形变网格。
为了使得生成器网络NetG输出的形变网格g所对应的形变在病变区域LA处产生较大形变,在生成器网络NetG的损失函数中引入了病变区域LA大形变项;在所构造的病变区域LA大形变项损失函数中,对病变区域LA标注区域和非病变区域LA标注区域的形变量赋予不同的权重,即
Figure BDA0004019066910000113
Figure BDA0004019066910000114
其中N为△g内像素点个数,R为△g内属于病变区域LA位置的任意一点,R′为△g内不属于病变区域LA位置的任意一点,a为常数(a>1),△g=g-g0,g0为原图-原图的形变网格。
因此,将生成器网络NetG的损失函数定义为:LG=LG12×Lsmooth3×LLA_deform,其中λ2和λ3为权重参数。
S3-6:构造训练过程中监测指标项Lmonitor=NetD(Xtrain)-NetD(Xtrain_aug),用于观察WGAN-GP模型的训练趋势。
S3-7:分别通过损失函数LD和LG交替迭代优化判别器网络NetD和生成器网络NetG,观察监测指标项Lmonitor的收敛情况,当监测指标项Lmonitor曲线收敛时停止模型训练,保存训练后的判别器网络NetD和生成器网络NetG的模型参数。
每一个训练周期内,通过损失函数LD迭代优化判别器网络NetD,通过损失函数LG迭代优化生成器网络NetG,相关训练参数如下表所示:
Figure BDA0004019066910000121
S4:使用医疗内窥镜图像测试集Xtest评估训练完成后的生成器网络NetG,评估性能符合需求的生成器网络NetG即为网格形变数据增强模型modelaug
S5:采样满足标准正态分布的随机噪声{z1,z2,…,zk},分别将随机噪声{z1,z2,…,zk)与医疗内窥镜图像X作为网格形变数据增强模型modelaug的模型输入,输出对应的形变网格{g1,g2,…,gk),将医疗内窥镜图像X分别和形变网格{g1,g2,…,gk}进行warp操作后得到数据增强后的医疗内窥镜图像{Xaug_1,Xaug_2,…,Xaug_k}。
如图10所示为本发明提供的基于WGAN-GP模型的网格形变数据增强方法在胃肠道内窥镜息肉分割数据集上的应用示例,可见本方法不仅可以对内窥镜图像进行数据增强,还可以对相应的病变区域标签进行同步的数据增强。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (8)

1.一种基于WGAN-GP模型的网格形变数据增强方法,其特征在于,包括以下步骤:
S1:构建训练数据集和测试数据集,将归一化后的医疗内窥镜图像数据按照7:3的比例划分为训练集Xtrain和测试集Xtest,将训练集图像及其相对应的病变区域LA标注打包成训练数据集;
S2:构建WGAN-GP模型,所述WGAN-GP模型包括基于残差结构的生成器网络NetG和基于残差结构的判别器网络NetD
S3:基于训练数据集迭代训练WGAN-GP模型,设置WGAN-GP模型中判别器网络NetD和生成器网络NetG的损失函数,分别通过损失函数交替迭代优化判别器网络NetD和生成器网络NetG,直至模型收敛,保存模型参数;
S4:使用医疗内窥镜图像测试集Xtest评估训练完成后的生成器网络NetG,评估性能符合需求的生成器网络NetG即为网格形变数据增强模型modelaug
S5:采样满足标准正态分布的随机噪声{z1,z2,…,zk},分别将随机噪声{z1,z2,…,zk}与医疗内窥镜图像X作为网格形变数据增强模型modelaug的模型输入,输出对应的形变网格{g1,g2,…,gk},将医疗内窥镜图像X分别和形变网格{g1,g2,…,gk}进行warp操作后得到数据增强后的医疗内窥镜图像Xaug_1,Xaug_2,…,Xaug_k}。
2.根据权利要求1所述的一种基于WGAN-GP模型的网格形变数据增强方法,其特征在于,步骤S2中所述生成器网络NetG为基于残差结构的卷积神经网络,包含第一残差结构块、第二残差结构块、上采样块、卷积层、批归一化层、最大值池化层、线性层、reshape层以及Tanh激活函数,其中第一残差结构块和第二残差结构块由卷积层、批归一化层、LeakyReLU激活函数交替串联组成,并通过残差连接将低层输出与高层输出进行连接;上采样块由卷积层、批归一化层、LeakyReLU激活函数、上采样层交替串联组成。
3.根据权利要求2所述的一种基于WGAN-GP模型的网格形变数据增强方法,其特征在于,所述生成器网络NetG的输入为医疗内窥镜图像训练集Xtrain和随机噪声z,生成器网络NetG的输出为形变网格g。
4.根据权利要求1所述的一种基于WGAN-GP模型的网格形变数据增强方法,其特征在于,步骤S2中所述判别器网络NetD为基于残差结构的卷积神经网络,包含第三残差结构块、第四残差结构块、卷积层、层归一化层、最大值池化层、自适应平均池化层、reshape层以及线性层,其中第三残差结构块和第四残差结构块由卷积层、层归一化层、ReLU激活函数交替串联组成,并通过残差连接将低层输出与高层输出进行连接。
5.根据权利要求4所述的一种基于WGAN-GP模型的网格形变数据增强方法,其特征在于,所述判别器网络NetD的输入为医疗内窥镜图像训练集Xtrain和形变图像Xtrain_aug,判别器网络NetD的输出为对应的值NetD(Xtrain)和NetD(Xtrain_aug)。
6.根据权利要求1所述的一种基于WGAN-GP模型的网格形变数据增强方法,其特征在于,步骤S3具体包含以下步骤:
S3-1:初始化生成器网络NetG和判别器网络NetD
S3-2:将医疗内窥镜图像训练集Xtrain和随机噪声z作为NetG的输入,输出形变网格g;
S3-3:将医疗内窥镜图像训练集Xtrain和形变网格g进行warp操作得到形变图像Xtrain_aug,将医疗内窥镜图像训练集Xtrain和形变图像Xtrain_aug输入判别器网络NetD,输出对应的值NetD(Xtrain)和NetD(Xtrain_aug);
S3-4:构造判别器网络NetD的损失函数为LD=LD11×LGrad_penalty,其中LD1=-NetD(Xtrain)+NetD(Xtrain_aug),LGrad_penalty为梯度惩罚项,λ1为权重参数;
S3-5:构造生成器网络NetG的损失函数为LG=LG12×Lsmooth3×LLA_deform,其中LG1=-NetD(Xtrain_aug),Lsmooth为平滑形变项,LLA_deform为病变区域LA大形变项,λ2和λ3为权重参数;
S3-6:构造训练过程中监测指标项Lmonitor=NetD(Xtrain)-NetD(Xtrain_aug),用于观察WGAN-GP模型的训练趋势;
S3-7:分别通过损失函数LD和LG交替迭代优化判别器网络NetD和生成器网络NetG,观察监测指标项Lmonitor的收敛情况,当监测指标项Lmonitor曲线收敛时停止模型训练,保存训练后的判别器网络NetD和生成器网络NetG的模型参数。
7.根据权利要求6所述的一种基于WGAN-GP模型的网格形变数据增强方法,其特征在于,步骤S3-5中所述平滑形变项Lsmooth为约束生成器网络NetG产生平滑形变的损失函数项,
Figure FDA0004019066900000031
Figure FDA0004019066900000032
其中N为Δg内像素点个数,p表示Δg内任意一点,px和py分别为p在X方向和Y方向的位置,Δg=g-g0,g0为原图-原图的形变网格。
8.根据权利要求6所述的一种基于WGAN-GP模型的网格形变数据增强方法,其特征在于,步骤S3-5中所述病变区域LA大形变项LLA_deform为鼓励病变区域LA产生较大形变的损失函数项,
Figure FDA0004019066900000033
Figure FDA0004019066900000034
其中N为Δg内像素点个数,R为Δg内属于病变区域LA位置的任意一点,R'为Δg内不属于病变区域LA位置的任意一点,α为常数(α>1),Δg=g-g0,g0为原图-原图的形变网格。
CN202211682574.2A 2022-12-27 2022-12-27 一种基于wgan-gp模型的网格形变数据增强方法 Pending CN115937038A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211682574.2A CN115937038A (zh) 2022-12-27 2022-12-27 一种基于wgan-gp模型的网格形变数据增强方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211682574.2A CN115937038A (zh) 2022-12-27 2022-12-27 一种基于wgan-gp模型的网格形变数据增强方法

Publications (1)

Publication Number Publication Date
CN115937038A true CN115937038A (zh) 2023-04-07

Family

ID=86650735

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211682574.2A Pending CN115937038A (zh) 2022-12-27 2022-12-27 一种基于wgan-gp模型的网格形变数据增强方法

Country Status (1)

Country Link
CN (1) CN115937038A (zh)

Similar Documents

Publication Publication Date Title
CN111627019B (zh) 一种基于卷积神经网络的肝脏肿瘤分割方法及系统
CN113077471B (zh) 一种基于u型网络的医学图像分割方法
CN110599528B (zh) 一种基于神经网络的无监督三维医学图像配准方法及系统
CN111401480B (zh) 一种基于融合注意力机制的乳腺mri自动辅助诊断新方法
CN113674281B (zh) 一种基于深度形状学习的肝脏ct自动分割方法
CN113674253B (zh) 基于U-Transformer的直肠癌CT影像自动分割方法
CN112258530A (zh) 基于神经网络的计算机辅助肺结节自动分割方法
CN112150425A (zh) 一种基于神经网络的无监督血管内超声图像配准方法
CN112489061A (zh) 一种基于多尺度信息与并行注意力机制的深度学习肠道息肉分割方法
CN109492668B (zh) 基于多通道卷积神经网络的mri不同期多模图像表征方法
CN115578404B (zh) 一种基于深度学习的肝脏肿瘤图像增强和分割的方法
CN112381164B (zh) 一种基于多分支注意力机制的超声图像分类方法和装置
CN115512110A (zh) 一种涉及跨模态注意力机制的医学图像肿瘤分割方法
CN113782184A (zh) 一种基于面部关键点与特征预学习的脑卒中辅助评估系统
Wang et al. Anomaly segmentation in retinal images with poisson-blending data augmentation
CN112270667A (zh) 一种基于ti-rads的集成深度学习多标签识别方法
CN113643261B (zh) 一种基于频率注意网络的胸肺部疾病诊断方法
CN112381822B (zh) 一种用于处理肺部病灶区图像的方法和相关产品
CN112750110A (zh) 基于神经网络对肺部病灶区进行评估的评估系统和相关产品
CN115937038A (zh) 一种基于wgan-gp模型的网格形变数据增强方法
CN113362360B (zh) 基于流体速度场的超声颈动脉斑块分割方法
CN114565626A (zh) 基于PSPNet改进的肺部CT图像分割算法
CN113902738A (zh) 一种心脏mri分割方法及系统
CN111798500B (zh) 一种基于层次邻域谱特征的微分同胚非刚性配准算法
CN110570417B (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