当前位置: 首页 > article >正文

【从VAE到LDM】Variational Auto Encoder原理以及关于Latent Diffusion的思考

论文链接:High-Resolution Image Synthesis with Latent Diffusion Models
官方实现:CompVis/latent-diffusion、CompVis/stable-diffusion
视频讲解:一个视频看懂VAE的原理以及关于latent diffusion的思考

前言

目前的扩散模型范式基本上都变成做Latent Diffusion,包括像DIT(diffusion Transformer),实际上也是基于Latent Diffusion的思想,它只是把里面的UNet换成了Transformer。

传统的扩散模型一般是输入一张图像 x x x,然后对它进行增噪,再通过去噪采样得到一个新的图像 x ^ \hat x x^。但这个过程的计算都是在Pixel Space上面进行的,采样空间太大,学习的噪声维度和图像的维度是相同的。当进行高分辨率图像生成时,需要的计算资源会急剧增加,虽然 DDIM 等工作已经对此有所改善,但效果依然有限。

因此就产生了Latent Space上进行扩散操作的想法,经过测试证明是可以的。但压缩到一个比较小的程度,比如说压缩到32*32大小,会有信息的损失。作者认为压缩分为两种:感知压缩和语义压缩。

作者认为感知压缩损失的这些信息实际上是很好学到的,因为它有很强的相关性。但是整体语义压缩损失实际上是比较重要的,在Latent Diffusion模型里,实际上是用扩散模型去学习语义。

实际上的过程:

  • 先把 x x x编码成一个隐藏向量 z z z
  • z z z增噪去噪,学得 z ^ \hat z z^
  • 解码得到图像 x ^ \hat x x^
    在这里插入图片描述

注意:
Encoder和Decoder的过程实际上是在做感知方面的压缩和解码,中间的扩散模型处理语义信息Encoder和Decoder就来自于VAE,Variational Auto Encoder其实是从AutoEncoder来的,用了一些方法把自编码器变成了生成模型,这个方法在2013年提出。实际上VAE的结果是不太好的,那为什么把中间换成扩散模型效果会变好?就这说明编码器和解码器的能力是足够的,问题出在中间的部分。

所以学习Latent Diffusion之前必须了解VAE的原理

问题定义

什么是AE(Auto Encoder)

AutoEncoder自编码器,它实际上是一种自监督学习图像特征的一个方法,它的原理很简单:
在这里插入图片描述
首先有一个 x i x_i xi,把它输入到Encoder中编码成一个 z i z_i zi,再把这个 z i z_i zi输入到Decoder中,给它解码成 x ^ i \hat x_i x^i,这就是重构图像,然后再去跟真实的 x i x_i xi做重构损失的计算,最后来优化编码器和解码器。如果训练好给定一个 z i z_i zi,就可以得到重构的比较好的 x i x_i xi,就可以认为 z i z_i zi实际上已经是学到了 x i x_i xi中比较丰富的细节特征了。

但是AE为什么不能作为生成模型?

假设现在有一个数据集它里边有n个图像 [ x 1 , x 2 . . . . x n ] {[x_1,x_2....x_n]} [x1,x2....xn],通过Encoder可以给它编码成 [ z 1 , z 2 . . . . z n ] {[z_1,z_2....z_n]} [z1,z2....zn]。因为生成模型有一个核心的诉求:生成一个全新的图像,即数据集里没有出现过的。现在如果想要一个全新的 Z Z Z,但肯定不能从 [ z 1 , z 2 . . . . z n ] {[z_1,z_2....z_n]} [z1,z2....zn]中来,但是也不能随便来。如果知道 Z Z Z服从哪一种分布,就可以从这个分布中取到一个新的 z z z,就可以通过解码得到一张新的图片。但训练时没有对 z z z的分布进行约束,导致 z z z的分布可能是不规律的,就不能生成新的图像。

假设 z z z服从正态分布 P ( z ) ∽ N ( 0 , 1 ) P(z)∽N(0,1) P(z)N(0,1)
如果假设 z z z服从正态分布 P ( z ) ∽ N ( 0 , 1 ) P(z)∽N(0,1) P(z)N(0,1)的话,很难确定正态分布N(0,1)中每个点和真实数据中每个点的对应关系,即 x i x_i xi z i z_i zi之间的对应关系

假设已经知道后验概率 P θ ( z ∣ x i ) P_θ(z|x_i) Pθ(zxi)
此时不再去找确切的 z z z了,而是假设已经知道后验概率 P θ ( z ∣ x i ) P_θ(z|x_i) Pθ(zxi),如果知道了这个分布,就可以从中取得一个新的 z i z_i zi,然后 z i z_i zi再输入到Decoder中,输出的 x ^ i \hat x_i x^i就可以跟 x i x_i xi进行重构损失的计算来优化Decoder。

思路的转变
这时思路就从想优化的auto encoder变成一个生成模型变成了需要找到 P θ ( z ∣ x i ) P_θ(z|x_i) Pθ(zxi)这么一个后验概率,接下来就可以通过变分贝叶斯的方法找到这个分布。

一、VAE

1. 定义

变分贝叶斯:如下图,有一个分布集合,其中有一点 P θ ( z ∣ x i ) P_θ(z|x_i) Pθ(zxi),其中有一个已知分布族,这个分布符合高斯分布,现在希望从这个分布族里找到一个最接近 P P P的分布并做近似。如果最开始取的是 x 1 x_1 x1点的分布,然后通过不断的优化,慢慢找到 x 2 x_2 x2这个地方。这个分布就是离 P P P的距离是最近的,所以最终通过优化找到这个分布。
在这里插入图片描述
从上面的高斯分布族取出一个分布 q φ ( z ∣ x i ) q_φ(z|x_i) qφ(zxi),此时需要优化 φ ,就需要找到已知分布中离已知 P 最近的一个分布。自然就想到一个度量方法,那就是KL散度。现在的问题就变成了最小化 K L ( q φ ( z ∣ x i ) ∣ ∣ P θ ( z ∣ x i ) ) KL(q_φ(z|x_i) || P_θ(z|x_i)) KL(qφ(zxi)∣∣Pθ(zxi))。如果以KL散度作为这个距离的衡量标准的话(虽然它不是一个严格意义的距离,它不服从对称性),这个方法就叫做变分贝叶斯法。

2. 变分下界

想要最小化 K L ( q φ ( z ∣ x i ) ∣ ∣ P θ ( z ∣ x i ) ) KL(q_φ(z|x_i) || P_θ(z|x_i)) KL(qφ(zxi)∣∣Pθ(zxi)),但是这里面有一个难点就是这个我们不知道 P θ ( z ∣ x i ) P_θ(z|x_i) Pθ(zxi),所以就很难去通过计算方法求得它,另一种方法就是变分下界方法。

根据KL散度定义展开:
K L ( q φ ( z ∣ x i ) ∣ ∣ P θ ( z ∣ x i ) ) = ∫ q φ ( z ∣ x i ) l o g q φ ( z ∣ x i ) P θ ( z ∣ x i ) d z = ∫ q φ ( z ∣ x i ) l o g ( q φ ( z ∣ x i ) ) d z − ∫ q φ ( z ∣ x i ) l o g ( P θ ( z ∣ x i ) ) d z = E q φ ( z ∣ x i ) [ l o g ( q φ ( z ∣ x i ) ) ] − E q φ ( z ∣ x i ) [ l o g ( P θ ( z , x i ) ) ] + E q φ ( z ∣ x i ) [ l o g ( P θ ( x i ) ) ] = E q φ ( z ∣ x i ) [ l o g ( q φ ( z ∣ x i ) ) ] − E q φ ( z ∣ x i ) [ l o g ( P θ ( z , x i ) ) ] + l o g ( P θ ( x i ) ) = − E L B O + l o g ( P θ ( x i ) ) ⇨ l o g ( P θ ( x i ) ) = K L + E L B O KL(q_φ(z|x_i) || P_θ(z|x_i)) \\ = \int q_φ(z|x_i) log \frac{q_φ(z|x_i)}{P_θ(z|x_i)} dz \\ = \int q_φ(z|x_i) log (q_φ(z|x_i))dz - \int q_φ(z|x_i) log (P_θ(z|x_i))dz \\ = E q_φ(z|x_i) [log (q_φ(z|x_i))] - E q_φ(z|x_i) [log (P_θ(z,x_i))] +E q_φ(z|x_i) [log (P_θ(x_i))] \\ = E q_φ(z|x_i) [log (q_φ(z|x_i))] - E q_φ(z|x_i) [log (P_θ(z,x_i))] +log (P_θ(x_i)) \\ = -ELBO + log (P_θ(x_i)) \\ ⇨ log (P_θ(x_i)) = KL + ELBO KL(qφ(zxi)∣∣Pθ(zxi))=qφ(zxi)logPθ(zxi)qφ(zxi)dz=qφ(zxi)log(qφ(zxi))dzqφ(zxi)log(Pθ(zxi))dz=Eqφ(zxi)[log(qφ(zxi))]Eqφ(zxi)[log(Pθ(z,xi))]+Eqφ(zxi)[log(Pθ(xi))]=Eqφ(zxi)[log(qφ(zxi))]Eqφ(zxi)[log(Pθ(z,xi))]+log(Pθ(xi))=ELBO+log(Pθ(xi))log(Pθ(xi))=KL+ELBO

  • 第三行到第四行根据数学期望展开,后一项的 P θ ( z ∣ x i ) P_θ(z|x_i) Pθ(zxi)展开为 P θ ( z , x i ) / P θ ( x i ) P_θ(z,x_i)/P_θ(x_i) Pθ(z,xi)/Pθ(xi)
  • 第四行的最后一项,是服从 q φ ( z ∣ x i ) q_φ(z|x_i) qφ(zxi)的数学期望,但 E q φ ( z ∣ x i ) [ l o g ( P θ ( x i ) ) ] E q_φ(z|x_i) [log (P_θ(x_i))] Eqφ(zxi)[log(Pθ(xi))]里没有 z z z,所以这项实际上就是 l o g ( P θ ( x i ) ) log (P_θ(x_i)) log(Pθ(xi))
  • 第五行的 E q φ ( z ∣ x i ) [ l o g ( q φ ( z ∣ x i ) ) ] − E q φ ( z ∣ x i ) [ l o g ( P θ ( z , x i ) ) ] E q_φ(z|x_i) [log (q_φ(z|x_i))] - E q_φ(z|x_i) [log (P_θ(z,x_i))] Eqφ(zxi)[log(qφ(zxi))]Eqφ(zxi)[log(Pθ(z,xi))]被称为变分下界即 − E L B O -ELBO ELBO l o g ( P θ ( x i ) ) log (P_θ(x_i)) log(Pθ(xi))被称为边界似然
  • 第七行通过移项得

我们的目标是最小化 K L KL KL散度,因为在上面的推导过程中 z z z是随机变量,而 l o g ( P θ ( x i ) ) log (P_θ(x_i)) log(Pθ(xi))中不含 z z z,说明在这个里面它是一个定值,所以在 l o g ( P θ ( x i ) ) = K L + E L B O log (P_θ(x_i)) = KL + ELBO log(Pθ(xi))=KL+ELBO中最小化KL散度的任务就转换成了最大化 E L B O ELBO ELBO

3. 最大化 ELBO

E L B O = E q φ ( z ∣ x i ) [ l o g ( P θ ( z , x i ) ) ] − E q φ ( z ∣ x i ) [ l o g ( q φ ( z ∣ x i ) ) ] = E q φ ( z ∣ x i ) [ l o g P θ ( z ) P θ ( x i ∣ z ) ] − E q φ ( z ∣ x i ) [ l o g ( q φ ( z ∣ x i ) ) ] = E q φ ( z ∣ x i ) [ l o g P θ ( z ) ] + E q φ ( z ∣ x i ) [ l o g P θ ( x i ∣ z ) ] − E q φ ( z ∣ x i ) [ l o g ( q φ ( z ∣ x i ) ) ] = − ( E q φ ( z ∣ x i ) [ l o g ( q φ ( z ∣ x i ) ) ] − E q φ ( z ∣ x i ) [ l o g P θ ( z ) ] ) + E q φ ( z ∣ x i ) [ l o g P θ ( x i ∣ z ) ] = − E q φ ( z ∣ x i ) [ l o g ( q φ ( z ∣ x i ) ) l o g P θ ( z ) ] + E q φ ( z ∣ x i ) [ l o g P θ ( x i ∣ z ) ] = − K L ( q φ ( z ∣ x i ) ∣ ∣ P θ ( z ) ) + E q φ ( z ∣ x i ) [ l o g P θ ( x i ∣ z ) ] ELBO = E q_φ(z|x_i) [log (P_θ(z,x_i))] - E q_φ(z|x_i) [log (q_φ(z|x_i))] \\ = E q_φ(z|x_i) [logP_θ(z)P_θ(x_i|z)] - E q_φ(z|x_i) [log (q_φ(z|x_i))] \\ = E q_φ(z|x_i) [logP_θ(z)] + E q_φ(z|x_i)[logP_θ(x_i|z)] - E q_φ(z|x_i) [log (q_φ(z|x_i))] \\ = -(E q_φ(z|x_i) [log (q_φ(z|x_i))] - E q_φ(z|x_i) [logP_θ(z)]) + E q_φ(z|x_i)[logP_θ(x_i|z)] \\ = -E q_φ(z|x_i) [\frac{log (q_φ(z|x_i))}{logP_θ(z)}] + E q_φ(z|x_i)[logP_θ(x_i|z)] \\ = -KL(q_φ(z|x_i) || P_θ(z)) + E q_φ(z|x_i)[logP_θ(x_i|z)] ELBO=Eqφ(zxi)[log(Pθ(z,xi))]Eqφ(zxi)[log(qφ(zxi))]=Eqφ(zxi)[logPθ(z)Pθ(xiz)]Eqφ(zxi)[log(qφ(zxi))]=Eqφ(zxi)[logPθ(z)]+Eqφ(zxi)[logPθ(xiz)]Eqφ(zxi)[log(qφ(zxi))]=(Eqφ(zxi)[log(qφ(zxi))]Eqφ(zxi)[logPθ(z)])+Eqφ(zxi)[logPθ(xiz)]=Eqφ(zxi)[logPθ(z)log(qφ(zxi))]+Eqφ(zxi)[logPθ(xiz)]=KL(qφ(zxi)∣∣Pθ(z))+Eqφ(zxi)[logPθ(xiz)]

  • 第二行将 l o g ( P θ ( z , x i ) ) log (P_θ(z,x_i)) log(Pθ(z,xi))展开为 P θ ( z ) P θ ( x i ∣ z ) P_θ(z)P_θ(x_i|z) Pθ(z)Pθ(xiz)
  • 第三行再 l o g log log展开;
  • 第四行移项并合并得到第五行;
  • 第五行中的第一项就是前面的 K L KL KL散度;后一项 l o g P θ ( x i ∣ z ) logP_θ(x_i|z) logPθ(xiz),为给定 z z z之后, x i x_i xi的最大似然,在生成模型中最大似然用模型来学。

推理步骤:
此时最大化 E L B O ELBO ELBO的任务就变成了最小化 K L ( q φ ( z ∣ x i ) ∣ ∣ P θ ( z ) ) KL(q_φ(z|x_i) || P_θ(z)) KL(qφ(zxi)∣∣Pθ(z)),同时最大化 E q φ ( z ∣ x i ) [ P θ ( x i ∣ z ) ] E q_φ(z|x_i)[P_θ(x_i|z)] Eqφ(zxi)[Pθ(xiz)]

q φ ( z ∣ x i ) ∽ N ( μ , σ 2 ) q_φ(z|x_i)∽N(\mu,\sigma^2) qφ(zxi)N(μ,σ2)是从高斯分布中取的随便一个分布

当时我们希望 q φ ( z ∣ x i ) ∽ N ( μ , σ 2 ) q_φ(z|x_i)∽N(\mu,\sigma^2) qφ(zxi)N(μ,σ2)和真实的后验概率越近越好

通过推导发现, q φ ( z ∣ x i ) ∽ N ( μ , σ 2 ) q_φ(z|x_i)∽N(\mu,\sigma^2) qφ(zxi)N(μ,σ2)需要和先验概率 P θ ( z ) P_θ(z) Pθ(z)越近越好,而最开始我们已经假设了 P θ ( z ) ∽ N ( 0 , 1 ) P_θ(z)∽N(0,1) Pθ(z)N(0,1),这就是为什么我们希望 q φ ( z ∣ x i ) ∽ N ( μ , σ 2 ) q_φ(z|x_i)∽N(\mu,\sigma^2) qφ(zxi)N(μ,σ2)逼近标准正态分布

4. 思路整理

  1. 先希望得到 P θ ( z ∣ x i ) P_θ(z|x_i) Pθ(zxi),目的是得到与 x i x_i xi对应的 z z z,这样好优化decode
  2. 用变分贝叶斯的方法得到 q φ ( z ∣ x i ) q_φ(z|x_i) qφ(zxi) q φ ( z ∣ x i ) q_φ(z|x_i) qφ(zxi)就是离 P θ ( z ∣ x i ) P_θ(z|x_i) Pθ(zxi)最近的
  3. K L KL KL散度优化这两个之间的距离
  4. 通过推导转化为 E L B O ELBO ELBO
  5. E L B O ELBO ELBO又转化为 − K L ( q φ ( z ∣ x i ) ∣ ∣ P θ ( z ) ) + E q φ ( z ∣ x i ) [ P θ ( x i ∣ z ) ] -KL(q_φ(z|x_i) || P_θ(z)) + E q_φ(z|x_i)[P_θ(x_i|z)] KL(qφ(zxi)∣∣Pθ(z))+Eqφ(zxi)[Pθ(xiz)]。第一项就是前面的 K L KL KL散度;后一项 l o g P θ ( x i ∣ z ) logP_θ(x_i|z) logPθ(xiz),为给定 z z z之后, x i x_i xi的最大似然,在生成模型中最大似然用模型来学。所以 P θ ( x i ∣ z ) P_θ(x_i|z) Pθ(xiz)就是decode; l o g P θ ( x i ∣ z ) logP_θ(x_i|z) logPθ(xiz)前面有一个期望值 E E E,这个期望是服从 q φ ( z ∣ x i ) q_φ(z|x_i) qφ(zxi)的。
    换句话说,如果在每次编码中得到的 q φ ( z ∣ x i ) ∽ N ( μ , σ 2 ) q_φ(z|x_i)∽N(\mu,\sigma^2) qφ(zxi)N(μ,σ2)分布后,要从这个里面重新取一个 z z z,通过蒙特拉洛的方法,让每次取得的 z z z都能最大的还原 x i x_i xi,这就是为什么VAE的训练过程是这样的(下一节解释)。

5. VAE的训练过程

VAE的训练过程如下图:
在这里插入图片描述

  • 先输入一个 x x x到encode
  • 得到一个 μ , σ \mu,\sigma μ,σ,通过这两个值得到高斯分布 N ( μ , σ 2 ) N(\mu,\sigma^2) N(μ,σ2)
  • 从高斯分布中采样一个 z z z,再把这个 z z z输入到decode中得到 x ^ \hat x x^,然后和 x x x做重构损失计算
  • 期间通过 K L KL KL散度尽量让 N ( μ , σ 2 ) N(\mu,\sigma^2) N(μ,σ2)和N(0,1)做近似 K L ( N ( μ , σ 2 ) ∣ N ( 0 , 1 ) ) KL(N(\mu,\sigma^2) | N(0,1)) KL(N(μ,σ2)N(0,1))

注意:

  • N ( μ , σ 2 ) N(\mu,\sigma^2) N(μ,σ2)就是前面讲的 q φ ( z ∣ x i ) q_φ(z|x_i) qφ(zxi),对 P θ ( z ∣ x i ) P_θ(z|x_i) Pθ(zxi)做最佳近似, N ( 0 , 1 ) N(0,1) N(0,1)是先验分布 P θ ( z ) P_θ(z) Pθ(z)
  • 找到了真实的后验概率 P θ ( z ∣ x i ) P_θ(z|x_i) Pθ(zxi)的最佳近似 q φ ( z ∣ x i ) q_φ(z|x_i) qφ(zxi)之后,就可以从 q φ ( z ∣ x i ) q_φ(z|x_i) qφ(zxi)中去一个 z z z,我们认为这个 z z z是从后验分布离拿到的,所以就和 x i x_i xi是对应的;
  • z z z通过decode后得到的 x ^ \hat x x^,就可以和真实的 x x x做重构损失计算,这一步就是前面推导出的 E q φ ( z ∣ x i ) [ l o g P θ ( x i ∣ z ) ] E q_φ(z|x_i)[logP_θ(x_i|z)] Eqφ(zxi)[logPθ(xiz)]
  • 有的论文认为不断的从高斯分布中采样 z 1 , z 2 z_1,z_2 z1,z2,让它们都去重构 x i x_i xi,这也是一种数据增强的方法,这就使得decode更具有泛化性。

二、Latent Space

Latent Space是什么样的呢?
在这里插入图片描述
这是在MNIST数据集上做了一个训练,图里一共有10种颜色,代表是0到9,做了对比试验。这个就是隐空间的分布,这是一个二维的分布,实际上隐空间的维度可能是很高的,这里只是做了一个投影。

  1. 只有重构损失,即AE:有区域,无规则,不服从正态分布,所以很难取出一个 z z z
  2. 只有 K L KL KL散度:无区域,有规则,即正态分布;
  3. VAE:有区域,有规则;

从这个图可知:

  • 重构损失会导致最终的特征是否有区域性;
  • K L KL KL散度对应最终的特征是否有规则;
  • 虽然上述分布是在隐空间中的,但实际上也是一种特征图,AE就是在学特征,其中的每个向量都是特征向量,而相似的特征必然在一个维度上,将其投影到平面上是比较相近的;
  • 图3可以看到这是一个从左上角到右下角连续的过程,最后生成了数字2;
  • K L KL KL使得隐空间符合设好的已知分布,重构损失会学到一定的特征,而学特征就是Decode的过程,如果给定了一个 z z z,它能很好的还原一个 x x x,那么我们认为 z z z已经包含了 x x x比较丰富的特征。
  • 图2到3加了重构损失后,就是对图2做了聚类的操作,把图像特征更像的聚集到一起,但是并不影响整体的分布。

三、LDM vs VAE

这两个模型到底差在哪?VAE的效果很一般,虽然LDM也用了VAE编码器解码器的结构,但效果好很多。

LDM训练过程简单示意图:
在这里插入图片描述
VAE训练过程简单示意图:
在这里插入图片描述

问题一(为什么LDM更好):

  • 我们认为LDM和VAE的先验概率都是服从 P ( z ) ∽ N ( 0 , 1 ) P(z)∽N(0,1) P(z)N(0,1)
  • 在训练VAE时,用了最小化 K L ( q ∣ ∣ P ( z ) ) KL(q || P(z)) KL(q∣∣P(z))的方法,让数据尽量接近 N ( 0 , 1 ) N(0,1) N(0,1) q q q是后验概率 P ( z ∣ x i ) P(z|x_i) P(zxi)的近似值。这里就会产生一个问题, P ( z ∣ x i ) P(z|x_i) P(zxi)不一定是真实的高斯分布,即后验概率,所以在做 q q q P ( z ∣ x i ) P(z|x_i) P(zxi)的近似值时就会产生偏差,生成效果就不是很好;
  • Diffusion的训练过程时不断的增噪,最后得到 N ( 0 , 1 ) N(0,1) N(0,1),采样时直接选取一个 x ^ i \hat x_i x^i是可以直接生成图像的,所以扩散模型通过不断的增噪去噪学得了标准正态分布 N ( 0 , 1 ) N(0,1) N(0,1)到真实空间 P ( x ) P(x) P(x)的映射关系;
  • 扩散模型可以学任意分布到标准正态之间的映射,VAE是假设的 P ( z ∣ x i ) P(z|x_i) P(zxi),用符合高斯分布的 q q q逼近它,因为 P ( z ∣ x i ) P(z|x_i) P(zxi)具体是什么分布是不知道,但现在可以直接用扩散模型做 P ( z ∣ x i ) P(z|x_i) P(zxi)到标准正态分布之间的映射,生成的时候就可以取标准正态分布中的 ε ε ε,得到一个服从 P ( z ∣ x i ) P(z|x_i) P(zxi)分布的 z z z,有了这个 z z z,就可以通过Decode解出 x ^ i \hat x_i x^i

问题二(VAE自己的问题):
有两张图的分布情况如下图,这两张图都采样了同一个点,这一个点都通过一个decode去学习
在这里插入图片描述
z z z重构出来的 x x x在和左边学习还是右边学习的情况时有分歧,所以模型在训练过程中就会不断的取average,所以它学习高频特征的能力很不好。

四、Latent Diffusion Models

扩散模型(Diffusion Models, DMs)通常在像素空间中操作,因此训练这些模型需要消耗大量的GPU资源,而推理过程基于序列化生成也会产生较高的计算成本。为在有限的计算资源下实现高质量的生成效果和灵活性,Latent Diffusion Models (LDM) 提出了将扩散过程转移到预训练自编码器的潜在空间中。相比于基于像素的 DMs,这种方法显著降低了计算需求,并且在空间维度上表现出更好的可扩展性。

此外,LDM 通过引入跨注意力(cross-attention)层,使模型能够更灵活地根据不同类型的条件输入进行图像生成。这一创新方法大幅提升了扩散模型的训练和推理效率,使得高质量图像生成在消费级显卡上成为可能。当前广受欢迎的 DALL-E 2 和 Stable Diffusion 就应用了 LDM 技术,进一步推动了扩散模型在图像生成领域的应用。

1. 方法介绍

本方法的整体结构如下图所示,主要分为三部分:最左侧的红框对应于感知图像压缩,中间的绿框对应 Latent Diffusion Models,右侧的白框表示生成条件,下面将分别介绍这三个部分。
在这里插入图片描述

1.1 感知图像压缩

LDM 把图像生成过程从原始的图像像素空间转换到了一个隐空间,具体来说,对于一个维度为 x ∈ R H × W × 3 \mathbf{x}\in\mathbb{R}^{H\times W\times 3} xRH×W×3的 RGB 图像,可以使用一个 encoder E \mathcal{E} E 将其转换为隐变量 z = E ( x ) \mathbf{z}=\mathcal{E}(\mathbf{x}) z=E(x),也可以用一个 decoder D \mathcal{D} D 将其从隐变量转换回像素空间 x ~ = D ( E ( x ) ) \tilde{\mathbf{x}}=\mathcal{D}(\mathcal{E}(\mathbf{x})) x~=D(E(x))在转换时会将图像下采样,作者测试了一系列下采样倍数 f ∈ { 1 , 2 , 4 , 8 , 16 , 32 } f\in\{1, 2, 4, 8, 16, 32\} f{1,2,4,8,16,32},发现下采样 4-16 倍的时候可以比较好地权衡效率和质量。
在这里插入图片描述

在进行图像压缩时,为了防止压缩后的空间是某个高方差的空间,需要进行正则化。作者使用了两种正则化,第一种是 KL-正则化,也就是将隐变量和标准高斯分布使用一个 KL 惩罚项进行正则化;第二种是 VQ-正则化,也就是使用一个 vector quantization 层进行正则化。

1.2 Latent Diffusion Models

实际上 latent diffusion models 和普通的扩散模型没有太大区别,只是因为从像素空间变到了隐空间,所以维度降低了。训练的优化目标也没有太大变化,普通的扩散模型优化目标为:
L D M = E x , ϵ ∼ N ( 0 , 1 ) , t [ ∣ ∣ ϵ − ϵ θ ( x t , t ) ∣ ∣ 2 2 ] L_{DM}=E_{x,ϵ∼N(0,1),t} [∣∣ϵ−ϵ_θ(x_t,t)∣∣_2^2] LDM=Ex,ϵN(0,1),t[∣∣ϵϵθ(xt,t)22]
而 Latent Diffusion Models 的优化目标只是套了一层 autoencoder:
L L D M = E E ( x ) , ϵ ∼ N ( 0 , 1 ) , t [ ∣ ∣ ϵ − ϵ θ ( x t , t ) ∣ ∣ 2 2 ] L_{LDM}=E_{\mathcal{E}_(x),ϵ∼N(0,1),t} [∣∣ϵ−ϵ_θ(x_t,t)∣∣_2^2] LLDM=EE(x),ϵN(0,1),t[∣∣ϵϵθ(xt,t)22]
在采样时,首先从隐空间随机采样噪声,在去噪后再用 decoder 转换到像素空间即可。

1.3 条件生成

为了进行条件生成,需要学习 ϵ θ ( x t , t , y ) \epsilon_\theta(\mathbf{x}_t,t,y) ϵθ(xt,t,y),这里使用的方法是在去噪网络中加入 cross attention 层,条件通过交叉注意力注入。

如果说不带控制信号的LDM是对 p ( z ) p(z) p(z)建模的话,那么加上控制信号的建模可以表述成 p ( z ∣ y ) p(z∣y) p(zy)。控制信号 y y y可以是文本、图片布局、轮廓图等。当有有一些更为细粒度的生成存在多个控制信号,即对 p ( z ∣ y 1 , ⋯ , y N ) p(z∣y_1,⋯,y_N) p(zy1,,yN)。为了使得生成过程考虑控制信息,作者在原有的Unet backbone上引入了一个交叉注意力机制,来融入控制信号。首先通过一个将 y y y投影到为一个中间表征 τ θ ( y ) ∈ R M × d τ τ_θ(y)∈R^{M×dτ} τθ(y)RM×dτ,随后和DM某层的输出进行融合。

假定 φ i ( z t ) φ_i(z_t) φi(zt)是隐变量 z z z在Unet在层 i i i时间步为 t t t时的输出,下面需要将 φ i ( z t ) φ_i(z_t) φi(zt) τ θ ( y ) τ_θ(y) τθ(y)用交叉注意力机制融合
A t t e n t i o n ( Q , K , V ) = s o f t m a x ( Q K T d ) V Q = W Q ( i ) φ i ( z t ) ; K = W K ( i ) τ θ ( y ) ; V = W V ( i ) τ θ ( y ) Attention(Q,K,V) = softmax(\frac{QK^T}{\sqrt d}) V \\ Q=W_Q^{(i)} φ_i(z_t); K=W_K^{(i)} τ_θ(y); V=W_V^{(i)} τ_θ(y) Attention(Q,K,V)=softmax(d QKT)VQ=WQ(i)φi(zt);K=WK(i)τθ(y);V=WV(i)τθ(y)
此时模型的优化目标为
L L D M = E E ( x ) , ϵ ∼ N ( 0 , 1 ) , t [ ∣ ∣ ϵ − ϵ θ ( x t , t , τ θ ( y ) ) ∣ ∣ 2 2 ] L_{LDM}=E_{\mathcal{E}_(x),ϵ∼N(0,1),t} [∣∣ϵ−ϵ_θ(x_t,t,τ_θ(y))∣∣_2^2] LLDM=EE(x),ϵN(0,1),t[∣∣ϵϵθ(xt,t,τθ(y))22]

2. LDM 使用示例

huggingface Diffusers 将各种 Diffusion Model Pipeline 都包装好了,使用 Diffusion model 就和使用 Transformers 一样地方便:

# !pip install diffusers transformers
from diffusers import DiffusionPipeline

model_id = "CompVis/ldm-text2im-large-256"

# load model and scheduler
ldm = DiffusionPipeline.from_pretrained(model_id)

# run pipeline in inference (sample random noise and denoise)
prompt = "A painting of a squirrel eating a burger"
images = ldm([prompt], num_inference_steps=50, eta=0.3, guidance_scale=6).images

# save images
for idx, image in enumerate(images):
    image.save(f"squirrel-{idx}.png")

3. LDM Pipeline

LDM 的 pipeline 可以简化表示为:Pipeline(prompt, num_inference_steps, latents)。我们暂时考虑没有 negative prompt 和 初始 latent 的输入,那么整个采样过程大致可以表示为:

  1. 首先采用了 BERT 架构模型对 prompt 进行处理,生成 text_hidden_state;同时生成随机噪声 latents
text_hidden_state = LDMBERT(prompt) # shape=[bs, len_seq, d_model] = [1, 77, 1280] 
latents = randn_tensor(latents_shape) 

对于 "CompVis/ldm-text2im-large-256",其中使用了 LDMBert, 参考 huggignface 的 LDMBert 实现,LDMBert 与传统 BERT 架构相似,规模不同,LDMBert 采用 32 层, hidden_size 为 1280,属实比 bert-base 大上不少。同时文本被 padding 到了固定的 77 长度,以此来保证文字的 hidden state 格式为 [batch_size, 77, 1280]

  1. 之后进行传统的扩散模型 backward process:
for t in self.progress_bar(self.scheduler.timesteps):
    noise_pred = self.unet(latents_input, t, encoder_hidden_states=context).sample
    # compute the previous noisy sample x_t -> x_t-1
    latents = self.scheduler.step(noise_pred, t, latents, **extra_kwargs).prev_sample

其中 UNET 为 UNet2DConditionModel,与传统 Unet 不同在于其应用了 Cross Attention 对文字以及图片信息进行综合处理,下文会对改模块做梳理。scheduler 可以选 DDIM 或者其他算法。

  1. 最后对 latent hidden state 进行 decode,生成图片:
latents = 1 / self.vqvae.config.scaling_factor * latents
image = self.vqvae.decode(latents).sample

4. LDM 中的 UNET

backward process 中的 self.unet(...),即 UNET2DCondition(sample, timestep, encoder_hidden_state) 前向推导可以看成五部分,(以下以 CompVis/ldm-text2im-large-256 为例介绍):

  • 准备 time steps :Timesteps 编码信息是 diffusion 中 predict noise residual 模型的标配:
# 经过两次映射得到 timesteps 对应的 embedding
t_emb = self.time_proj(timesteps)
emb = self.time_embedding(t_emb, timestep_cond)
  • pre-process: LDM 只用了一个 2D 卷积对输入的 hidden state 进行处理
sample = nn.Conv2d(
            in_channels, block_out_channels[0], kernel_size=conv_in_kernel, padding=conv_in_padding
        )(sample)
  • down sampling :down sampling 包括了 3 个 CrossAttnDownBlock2D, 和 1 个 DownBlock2D
# down sampling 大致前向推导
down_block_res_samples = (sample,)
for downsample_block in self.down_blocks:
    sample, res_samples = downsample_block(hidden_states=sample, temb=emb, scale=lora_scale)
    # 用于 UNET 的残差链接
    down_block_res_samples += res_samples

其中每个 CrossAttnDownBlock2D 大概前向过程为:

# CrossAttnDownBlock2D
def forward(self, hidden_states, temb, encoder_hidden_states=None)
	output_states = ()
    for resnet, attn in zip(self.resnets, self.attentions):
        hidden_states = resnet(hidden_states, temb)
        hidden_states = attn(
            hidden_states,
            encoder_hidden_states=encoder_hidden_states,
            cross_attention_kwargs=cross_attention_kwargs,
        ).sample
        output_states += (hidden_states,)

    # downsampler = Conv2D 
    hidden_states = downsampler(hidden_states)
    output_states += (hidden_states,)

    return hidden_states, output_states

CompVis/ldm-text2im-large-256 中,每个 CrossAttnDownBlock2D 包含了 2 个 attn(Transformer2DModel)以及 2 个 resnet (ResnetBlock2D)

文字与图像的交互就发生在 Transformer2DModel 当中。每个 Transformer2DModel 先对输入的图像数据进行预处理,将图片格式从如 (batch_size, channel, width, height)(batch_size, num_image_vectors) 转换为 (batch_size, len_seq, hidden_size),而后将 hidden_states 传入 1 层传统 Transformer layer(非 bert 或 GPT 类型),先对图像hidden_states进行 self-attention,而后结合 encoder_hidden_states 进行 cross attention 处理。

  • mid processing:
sample = MidBlock2DCrossAttn()(sample, 
                              emb,
                           encoder_hidden_states)

CompVis/ldm-text2im-large-256 中,upsampling 和 down sampling 之间采用 MidBlock2DCrossAttn 连接,MidBlock2DCrossAttn 包括了 1 个 1 层的 Transformer2DModel 以及 1 个 resnet ResnetBlock2D

  • upsampling :upsampling 采用的模块 UpBlocks 包括了 ("UpBlock2D", "CrossAttnUpBlock2D", "CrossAttnUpBlock2D", "CrossAttnUpBlock2D"),各个模块的架构与 down sampling 中的模块相似。
# upsample_block
for i, upsample_block in enumerate(self.up_blocks):
    sample = upsample_block(
                    hidden_states=sample,
                    temb=emb,
                    res_hidden_states_tuple=res_samples,
                    upsample_size=upsample_size,
                    scale=lora_scale,
                )
  • post-process
# GroupNorm
sample = self.conv_norm_out(sample)
# Silu
sample = self.conv_act(sample)
# Conv2d(320, 4, kernel=(3,3), s=(1,1), padding=(1,1))
sample = self.conv_out(sample)

总结起来,down sampling,midprocess,upsampling 三个步骤中都涉及到了 Transformer2DModel ,实现多模态的信息交互。

参考资料
DIFFUSION 系列笔记| Latent Diffusion Model
CompVis/ldm-text2im-large-256


http://www.kler.cn/a/383199.html

相关文章:

  • Python中 yield的作用
  • 梧桐数据库之以识别优质的移动服务套餐为例讲解SQL实现分享
  • 进程与线程+多线程优势
  • PCL 点云配准 精度评价指标均方根误差(RMSE)
  • 【redis】延迟双删策略
  • 一周模电速成(3) 超详细!入门小白速成!!!
  • 深入解析 Transformer 模型:编码器和解码器的完整实现
  • 易考八股文之谈谈对sentinel的理解和作用?
  • 关于路由笔记
  • Qt C++,Link custom data to the QCustomPlot data
  • 如何用Neo4j做知识图谱
  • yelp数据集上识别潜在的热门商家
  • 机器学习—前向传播的一般实现
  • 爬虫入门urllib 和 request(二)
  • Backend - Python 爬取网页数据并保存在Excel文件中
  • 【React】JSX规则
  • 产品经理必备秘籍:打造有效的产品 Roadmap
  • 原生鸿蒙应用市场:赋能开发者全生命周期服务体验
  • GAN的基本原理
  • Linux学习笔记之vim入门
  • 【数据结构】二叉树——层序遍历
  • HTML5+css3(伪类,动态伪类,结构伪类,否定伪类,UI伪类,语言伪类,link,hover,active,visited,focus)
  • 网络优化如何利用改IP软件解除地域限制
  • VBA02-初识宏——EXCEL录像机
  • Windows核心编程笔记——DLL基础
  • 【AI视频换脸整合包及教程】AI换脸新星:Rope——让换脸变得如此简单