InfoGAN:一种无监督生成方法

过去几年发表于各大 AI 顶会论文提出的 400 多种算法中,公开算法代码的仅占 6%,其中三分之一的论文作者分享了测试数据,约 54% 的分享包含“伪代码”。这是今年 AAAI 会议上一个严峻的报告。人工智能这个蓬勃发展的领域正面临着实验重现的危机,就像实验重现问题过去十年来一直困扰着心理学、医学以及其他领域一样。最根本的问题是研究人员通常不共享他们的源代码。

可验证的知识是科学的基础,它事关理解。随着人工智能领域的发展,打破不可复现性将是必要的。为此,PaperWeekly 联手百度 PaddlePaddle 共同发起了本次论文有奖复现,我们希望和来自学界、工业界的研究者一起接力,为 AI 行业带来良性循环。

作者丨黄泽宇

学校丨中山大学数学学院18级本科生

研究方向丨图像识别、VQA、生成模型和自编码器

论文复现代码:

http://aistudio.baidu.com/#/projectdetail/23600

GAN

生成对抗网络(Generative Adversarial Nets)是一类新兴的生成模型,由两部分组成:一部分是判别模型(discriminator)D(·),用来判别输入数据是真实数据还是生成出来的数据;另一部分是是生成模型(generator)G(·),由输入的噪声生成目标数据。GAN 的优化问题可以表示为:

其中 Pdata 是生成样本,noise 是随机噪声。而对于带标签的数据,通常用潜码(latent code)c 来表示这一标签,作为生成模型的一个输入,这样我们有:

然而当我们遇到存在潜在的类别差别而没有标签数据,要使 GAN 能够在这类数据上拥有更好表现,我们就需要一类能够无监督地辨别出这类潜在标签的数据,InfoGAN 就给出了一个较好的解决方案。

互信息(Mutual Information)

互信息是两个随机变量依赖程度的量度,可以表示为:

要去直接优化 I(c;G(z,c)) 是极其困难的,因为这意味着我们要能够计算后验概率(posterior probability)P(c|x),但是我们可以用一个辅助分布(auxiliary distribution)Q(c|x),来近似这一后验概率。这样我们能够给出互信息的一个下界(lower bounding):

InfoGAN

在 InfoGAN 中,为了能够增加潜码和生成数据间的依赖程度,我们可以增大潜码和生成数据间的互信息,使生成数据变得与潜码更相关:

▲图1. InfoGAN的整体结构图

由上面的,对于一个极大化互信息的问题转化为一个极大化互信息下界的问题,我们接下来就可以定义:

在论文的附录中,作者证明了:

于是:

故 LI (G, Q) 是互信息的一个下界。作者指出,用蒙特卡罗模拟(Monte Carlo simulation)去逼近 LI (G, Q) 是较为方便的,这样我们的优化问题就可以表示为:

实现

在实现中,D(x)、G(z, c) 和 Q(x) 分别用一个 CNN (Convolutional Neural Networks)、CNN、DCNN (DeConv Neural Networks) 来实现。同时,潜码 c 也包含两部分:一部分是类别,服从 Cat(K = N,p = 1/N),其中 N 为类别数量;另一部分是连续的与生成数据有关的参数,服从 Unif(−1,1) 的分布。

在此应指出,Q(c|x) 可以表示为一个神经网络 Q(x) 的输出。对于输入随机变量 z 和类别潜码 c,实际的 LI(G, Q) 可以表示为:

其中 · 表示内积(inner product),c 是一个选择计算哪个 log 的参数,例如 ci = 1 而 cj = 0(∀j = 1,2,···,i − 1,i + 1,···,n),那么 z 这时候计算出的 LI(G,Q) 就等于 log(Q(z,c)i)。这里我们可以消去 H(c),因为 c 的分布是固定的,即优化目标与 H(c) 无关:

而对于参数潜码,我们假设它符合正态分布,神经网络 Q(x) 则输出其预测出的该潜码的均值和标准差, 我们知道,对于均值 μ,标准差 σ 的随机变量,其概率密度函数为:

要计算参数潜码 c 的,就是要计算 log p(c),即:

设 Q(x) 输出的参数潜码 c 的均值 μ,标准差 σ 分别为 Q(x)μ 和 Q(x)σ,那么对于参数潜码 c:

同样的,我们可以消去 H(c),因为 c 的分布是固定的,那么:

实验

首先,通过和普通的 GAN 比较 LI ,作者证明了 InfoGAN 确实能够优化这一互信息的下界 2。

作者在 MNIST 手写数字数据集(3)、3D 面部数据集(4)、3D 椅子数据集(5)、SVHN 街景房号数据集(6)以及 CelebA 人脸数据集(7)上进行了模型的相关测试。

▲图2. MNIST手写字符数据集上的结果

▲图3. 3D面部数据集上的结果

▲图4. 3D椅子数据集上的结果

▲图5. SVHN街景房号数据集上的结果

▲图6. CelebA人脸数据集上的结果

作者展示了这些数据集上学习到的类别潜码(从上至下变化)和参数潜码(从左至右变化,由 -2 到 2),我们可以看出,InfoGAN 不仅能够很好地学习数据之间的类型差别,也能够很好地学习到数据本身的一些易于区分的特点,而且生成模型对这些特点的泛化能力还是很好的。

再论InfoGAN的LI

读完论文,我们发现,对于类别潜码,这个 LI 本质上是 x 与 G(z, c) 之间的 KL 散度:

也就是说:

而 min DKL(c||Q(G(z, c))) 意味着减小 c 与 Q(G(z, c)) 的差别。

▲图7. 普通GAN和InfoGAN的LI在训练过程中的比较

如果我们不考虑 Q(x)σ 的影响,LI 的优化过程:

而也意味着减小 c 与 Q(G(z, c))μ 的差。

再纵观整个模型,我们会发现这一对 LI 优化的过程,实质上是以 G 为编码器(Encoder), Q 为解码器(Decoder),生成的图像作为我们要编码的码(code),训练一个自编码器(Autoencoder),也就是说,作者口中的信息论优化问题,本质上是无监督训练问题。

关于PaddlePaddle

在 PaddlePaddle中,一个极为重要的概念即是 fluid.Program(),在官方文档里常见的 exe.run(program= fluid.default_startup_program())的 fluid.default_startup_program() 就是其中一个例子。

在这一使用中可以了解到,我们要用 exe.run() 中的 program 参数运行指定的 fluid.Program(),而官方文档指出,当该参数未指定时,会运行 fluid.default_main_program(),而 fluid.default_main_program() 代表的是未指定 fluid.Program() 的所有操作

注意,这里说的是“所有”,由于PaddlePaddle没有计算依赖检测机制,即使在计算 fetch_list 中的值的时候不会用到操作也会被计算,这一点与 TensorFlow 极其不同,作者本人在使用过程中踩了很大的坑,还望各位注意。在执行多种任务的时候不要一股脑全部写在 fluid.default_main_program() 之中, 这样极其浪费资源,也容易造成一些问题。

一个新的 fluid.Program() 被创建之后,可以在 fluid.program_guard() 中指定该 fluid.Program() 中的操作与变量:

PaddlePaddle中还需要注意的一点是,fluid.Variable 的命名空间是全局的,也就是说在同一或者不同 fluid. Program() 间,同名(fluid.Variable 的 name 属性相同)的 fluid.Variable 所指向的变量是相同的,所以同一名称在同一或者不同 fluid.Program () 中可以被使用多次,而不用担心 TensorFlow 中会出现的 reuse 问题。

要对一个操作的中的权值的名称进行定义(权值命名为 W1,偏置命名为 b1):

要在之后使用这些 fluid.Variable,例如在 Optimizer 中使用:

#可以直接用名称指代对应的fluid.Variable

parameter_list = ["W1","b1"]

#构建optimizer

optimizer = fluid.optimizer.AdamOptimizer()

#指定optimizer优化的目标和对象

optimizer.minimize(loss, parameter_list=parameter_list)

在构建完基本的运算操作后,便可以开始初始化操作了:

#初始化fluid.Executor(指定执行程序位置)

exe = fluid.Executor(fluid.CPUPlace())

#执行fluid.default_startup_program(),在fluid.program_guard()中

#若没有指定初始化program,则默认为此program

exe.run(program=fluid.default_startup_program())

初始化完成后,可以开始训练啦:

GAN实现

生成对抗网络(Generative Adversarial Nets)是一类新兴的生成模型,由两部分组成:一部分是判别模型(discriminator)D(·),用来判别输入数据是真实数据还是生成出来的数据;另一部分是是生成模型(generator)G(·),由输入的噪声生成目标数据。GAN 的优化问题可以表示为:

其中 Pdata 是生成样本,noise 是随机噪声。我们用一个双层的 MLP 来演示:

通常,一个 GAN 的训练由两部分组成,第一部分是对 D(·) 进行训练,极大化目标函数:

第二部分是对 G(·) 进行训练,极小化目标函数:

以下是两部分优化的定义:

#参考Todd的LSGAN的实现,使用函数获取模型所有变量

defget_params(program, prefix):

all_params = program.global_block().all_parameters()

return[t.namefortinall_paramsift.name.startswith(prefix)]

#G优化程序

G_program = fluid.Program()

withfluid.program_guard(main_program = G_program):

#定义输入数据

Z = fluid.layers.data(name='Z', shape=[Z_dim], dtype='float32')

#执行相关模型的计算

G_sample = generator(Z)

D_fake = discriminator(G_sample)

#计算损失函数

G_loss =0.0- fluid.layers.reduce_mean(fluid.layers.log(D_fake +1e-8))

#定义optimizer优化的变量的范围

theta_G = get_params(G_program,"G")

G_optimizer = fluid.optimizer.AdamOptimizer()

G_optimizer.minimize(G_loss, parameter_list=theta_G)

#D优化程序

D_program = fluid.Program()

withfluid.program_guard(main_program = D_program):

Z = fluid.layers.data(name='Z', shape=[Z_dim], dtype='float32')

X = fluid.layers.data(name='X', shape=[784], dtype='float32')

#在使用数据集时,要注意相应接口传入数据的值的范围

#paddle.dataset.mnist中的数据,范围在[-1, 1]

#要将其转换到sigmoid函数的值域内

X = X *0.5+0.5

G_sample = generator(Z)

D_real = discriminator(X)

D_fake = discriminator(G_sample)

D_loss =0.0- fluid.layers.reduce_mean(fluid.layers.log(D_real +1e-8)

+ fluid.layers.log(1.0- D_fake +1e-8))

theta_D = get_params(G_program,"D")

D_optimizer = fluid.optimizer.AdamOptimizer()

D_optimizer.minimize(D_loss, parameter_list=theta_D)

在定义好这些之后,是时候开训练了:

#定义传入的数据

feeding_withx= {"X": np.array(X_mb).astype('float32'),

"Z": np.array(Z_noise).astype('float32')}

feeding = {"Z": np.array(Z_noise).astype('float32')}

#执行训练操作并获取当前损失函数的值

D_loss_curr = exe.run(feed = feeding_withx, program = D_program,

fetch_list = [D_loss])

G_loss_curr = exe.run(feed = feeding, program = G_program,

fetch_list = [G_loss])

若欲测试模型效果,可再定义一个 Inference:

然后再这样获取 samples:

feeding = {"Z": np.array(Z_noise).astype('float32')}

samples = exe.run(feed = feeding, program = Infer_program,

fetch_list = [G_sample])

后记

本文先前于今年 8 月完成,共享于 PaddlePaddle 论文复现群内,在 10 月 LSGAN 的复现公开之 后,参考该复现更改了模型参数命名和参数列表的实现方法,在此感谢 Todd 同学的复现对本文的帮助。

  • 发表于:
  • 原文链接:https://kuaibao.qq.com/s/20181026B10AE600?refer=cp_1026
  • 腾讯「云+社区」是腾讯内容开放平台帐号(企鹅号)传播渠道之一,根据《腾讯内容开放平台服务协议》转载发布内容。

扫码关注云+社区

领取腾讯云代金券