← 返回 2026-03-10

尺度空间扩散:将尺度空间理论融入扩散模型的统一框架 Scale Space Diffusion

Soumik Mukhopadhyay, Prateksha Udhayanan, Abhinav Shrivastava 📅 2026-03-09 👍 16 2026-07-13 08:35
图像合成 多分辨率生成 尺度空间理论 扩散模型 高效推理

揭示扩散过程与尺度空间的信息层级对应关系,通过广义线性退化实现多分辨率高效生成

前置知识

扩散模型(Diffusion Models)

扩散模型是一类生成模型,通过逐步向数据添加高斯噪声来构建前向过程,然后学习反向去噪过程来生成新样本。在DDPM中,前向过程定义为 $x_t = \sqrt{\alpha_t} x_{t-1} + \sqrt{1-\alpha_t}\epsilon$,其中 $\epsilon \sim \mathcal{N}(0, I)$。通过累积乘积 $\bar{\alpha}_t = \prod_{i=1}^t \alpha_i$,可以将任意时刻的噪声状态表示为 $x_t = \sqrt{\bar{\alpha}_t} x_0 + \sqrt{1-\bar{\alpha}_t}\epsilon$。模型通过预测噪声 $\epsilon_\theta(x_t, t)$ 或预测干净图像 $x_0$ 来学习反向过程。

本文的核心是扩展扩散模型框架,理解标准DDPM的前向、边际和后验分布是理解本文广义线性扩散过程的基础

尺度空间理论(Scale Space Theory)

尺度空间理论是计算机视觉中的基础框架,通过多尺度表示来分析图像结构。其核心思想是通过高斯金字塔(Gaussian Pyramid)对图像进行逐级低通滤波和下采样,从粗糙到精细地表示图像信息。随着尺度减小(分辨率降低),图像的细节逐渐丢失,只保留粗略结构。这种信息层级特性与扩散过程中噪声逐步去除信息的方式有惊人的相似性。

本文发现扩散时间步和尺度空间分辨率之间存在信息层级的对应关系,这是提出Scale Space Diffusion的核心动机

FID(Fréchet Inception Distance)

FID是评估生成模型质量的标准指标,通过计算生成图像和真实图像在Inception网络特征空间中的Fréchet距离来衡量。FID越低表示生成质量越好,计算公式为 $\text{FID} = ||\mu_r - \mu_g||^2 + \text{Tr}(\Sigma_r + \Sigma_g - 2(\Sigma_r\Sigma_g)^{1/2})$,其中 $\mu, \Sigma$ 分别是真实和生成分布的均值和协方差。

本文所有实验结果都通过FID进行评估,是衡量方法性能的核心指标

非各向同性高斯采样

当扩散过程中的退化算子 $M_t$ 是非标量的线性算子(如下采样)时,后验协方差 $\Sigma_{t \to t-1}$ 不再是各向同性的(即 $\sigma^2 I$ 的形式)。采样这种分布需要先用Lanczos算法数值计算协方差矩阵的平方根 $\Sigma_{t \to t-1}^{1/2}$,然后对标准高斯噪声 $\epsilon$ 进行变换:$\eta = \Sigma_{t \to t-1}^{1/2} \epsilon$。

这是Scale Space Diffusion实现的关键技术细节,忽略非各向同性会导致生成图像出现颜色平坦和饱和的伪影

UNet架构

UNet是一种编码器-解码器架构,通过跳跃连接(skip connections)将编码器的特征传递到解码器。在扩散模型中,UNet被广泛用作去噪网络,输入噪声图像 $x_t$ 和时间步 $t$,输出预测的噪声或干净图像。标准UNet通过下采样和上采样块在不同空间尺度上处理特征。

本文提出的Flexi-UNet是对标准UNet的重要改进,使其能够处理多分辨率输入和输出

研究动机

现有扩散模型在处理高分辨率图像生成时存在严重的计算效率问题。当我们观察扩散过程的中间状态时会发现一个有趣的现象:在高噪声水平下,扩散状态实际上只包含相当于低分辨率图像的信息量。例如,当扩散步 $t$ 较大时,信号项 $\sqrt{\bar{\alpha}_t} x_0$ 被噪声项 $\sqrt{1-\bar{\alpha}_t} \epsilon$ 主导,这意味着大部分像素的信号已经被噪声淹没。通过信息分析,我们定义信息量为信号主导像素的期望比例 $\text{Info}(t) = \mathbb{E}_{x_0}[1 - 2\Phi(-s(t)|x_0|)]$,其中 $s(t) = \sqrt{\bar{\alpha}_t/(1-\bar{\alpha}_t)}$ 是信噪比的平方根。然而,现有的扩散模型仍然以全分辨率处理这些高噪声状态,导致大量不必要的计算。以CelebA-256为例,DDPM基线需要87.31小时训练,497.03 GFLOPs的计算量。

本文的目标是本文的目标是正式建立扩散过程和尺度空间理论之间的数学联系,提出一个能够将尺度空间结构直接嵌入扩散过程的统一框架。具体而言,作者希望:(1)量化扩散时间步和尺度空间分辨率之间的信息对应关系;(2)开发一个广义线性扩散过程框架,支持任意线性退化算子;(3)设计能够处理多分辨率输入输出的神经网络架构;(4)在保持生成质量的同时显著提升训练和推理效率。最终目标是让扩散模型能够在适当的分辨率下处理不同噪声水平的状态,避免在高噪声状态下进行不必要的高分辨率计算。

与已有工作不同的是,尽管已有多项工作尝试结合尺度空间和扩散模型,但存在根本性的方法论缺陷。Cascaded Diffusion Model和Matryoshka Diffusion等方法仍然只在最高分辨率下操作或依赖多个独立模型,计算效率提升有限。Pyramidal Flow Matching和PixelFlow等方法在推理时近似尺度变化,但没有从数学上建模分辨率转换过程。UDPM尝试将模糊和子采样融入扩散,但假设后验协方差是各向同性的($\Sigma = \sigma^2 I$),这在实际的下采样退化中并不成立,因为模糊核通常会重叠。本文的独特切入角度是:首先建立扩散过程在广义线性退化下的完整数学理论,然后选择下采样作为退化算子来实现尺度空间的融合,从而得到真正的Scale Space Diffusion框架。

核心方法

本文的方法可以这样理解:想象你在看一张逐渐模糊的图片,先是失去细节(如发丝纹理),然后失去中等特征(如眼睛轮廓),最后只剩下大致的脸部轮廓。扩散过程的噪声添加也遵循类似的信息层级。作者的洞察是:既然高度噪声的扩散状态信息量等同于低分辨率图像,为什么不直接在低分辨率下处理它们呢?技术路线是:将标准扩散过程中的标量系数 $\sqrt{\alpha_t}$ 替换为通用线性算子 $M_t$(如下采样算子),建立广义线性扩散过程的完整数学框架,包括前向转移分布、边际分布和后验分布。然后设计Flexi-UNet架构,让网络能够根据输入分辨率动态激活不同的层,实现分辨率保持和分辨率提升的去噪操作。

本文的核心创新是将扩散过程从标量退化扩展到线性退化。在标准DDPM中,前向过程为 $x_t = \sqrt{\alpha_t} x_{t-1} + \sqrt{1-\alpha_t}\epsilon$,退化系数 $\sqrt{\alpha_t}$ 是标量。本文将其推广为 $x_t = M_t x_{t-1} + \eta_t$,其中 $M_t$ 是任意线性算子,$\eta_t \sim \mathcal{N}(0, \Sigma_{t|t-1})$。当选择下采样算子作为 $M_t$ 时,DDPM成为SSD的特例(对应恒等算子)。关键区别在于:(1)SSD的边际分布为 $q(x_t|x_0) = \mathcal{N}(M_{1:t}x_0, \sigma_t^2 I)$,其中 $M_{1:t} = M_t M_{t-1} \cdots M_1$ 是累积算子;(2)后验分布 $q(x_{t-1}|x_t, x_0)$ 的协方差 $\Sigma_{t \to t-1} = \sigma_{t-1}^2 I - \frac{\sigma_{t-1}^4}{\sigma_t^2} M_t^T M_t$ 是非各向同性的;(3)需要专门的采样技术(Lanczos算法)来处理非各向同性噪声。

方法步骤详情

Scale Space Diffusion的完整流程包括以下几个关键步骤:(1)**定义分辨率调度**:建立从扩散时间步 $t$ 到空间分辨率 $r(t)$ 的映射函数,使得分辨率随时间步单调递减。本文探索了多种调度策略,包括等间隔、凸衰减、Sigmoid形、Tanh形等,发现ConvexDecay 0.5在FID和训练时间之间取得最佳平衡。(2)**前向扩散**:使用累积算子 $M_{1:t}$ 和噪声标准差 $\sigma_t$ 生成噪声状态 $x_t = M_{1:t}x_0 + \sigma_t\epsilon$,其中 $M_{1:t}$ 包含多次2倍下采样操作。(3)**训练模型**:损失函数为 $\mathcal{L} = \mathbb{E}_{x_0,t,\epsilon}[\min(s^2(t), \gamma) \|x^{r(t-1)}_{0,\theta}(x_t, t) - \frac{M_{1:t-1}x_0}{a_{t-1}}\|^2_2]$,其中 $s^2(t)$ 是信噪比,$\gamma=5$ 是Min-SNR策略的截断参数。(4)**采样生成**:从最低分辨率 $r(T)$ 的随机高斯噪声开始,模型预测下一分辨率的干净图像 $\mu_{t-1} = a_{t-1} \cdot x^{r(t-1)}_{0,\theta}$,然后从后验分布采样 $x_{t-1} \sim \mathcal{N}(\mu_{t \to t-1}, \Sigma_{t \to t-1})$。(5)**非各向同性噪声采样**:对于分辨率变化的步骤,使用Lanczos算法计算 $\Sigma_{t \to t-1}^{1/2} \epsilon$ 来采样非各向同性噪声。

技术新颖性

本文的技术新颖性体现在多个层面:首先,在理论层面,建立了扩散模型在广义线性退化下的完整数学框架,包括Theorem 1(前向转移)、Theorem 2(后验分布)和Theorem 3(各向同性边际下的闭式解),这扩展了DDPM和Blurring Diffusion的理论基础。其次,在实现层面,提出了使用向量-雅可比积(vector-Jacobian product)来隐式计算线性算子的转置 $M_t^T v = \nabla_x \langle v, M_t x \rangle$,解决了图像缩放等隐式算子无法直接获得矩阵形式的难题。第三,在架构层面,Flexi-UNet通过动态激活UNet的不同层来处理多分辨率输入,高分辨率输入遍历完整UNet,低分辨率输入只经过深层,通过1×1卷积层调整通道维度。这种设计使得参数在不同分辨率间共享,同时支持分辨率保持和分辨率提升的去噪操作。第四,DDPM自然成为SSD在 $M_t = \sqrt{\bar{\alpha}_t} I$ 时的特例,这表明SSD是一个更一般的框架。

信息分析:扩散状态和尺度空间中的信息量
Figure 2: 信息分析:扩散状态和尺度空间中的信息量
Scale Space Diffusion的训练和推理流程
Figure 3: Scale Space Diffusion的训练和推理流程
分辨率调度策略及其对性能的影响
Figure 5: 分辨率调度策略及其对性能的影响
渐进式下采样退化下预测干净图像的动画
Figure 7: 渐进式下采样退化下预测干净图像的动画
渐进式下采样退化下噪声中间状态的动画
Figure 8: 渐进式下采样退化下噪声中间状态的动画

实验结果

本文在CelebA和ImageNet数据集上进行了全面的实验评估。在CelebA-64数据集上,SSD (2L)达到FID 2.14,训练时间62.63小时,而DDPM-ϵ基线FID为2.22,训练时间70.30小时。更重要的是,随着分辨率增加,效率优势更加明显:在CelebA-256上,SSD (6L)的FID为13.50,训练时间仅42.88小时,而DDPM基线FID为5.52,训练时间87.31小时——SSD的训练时间不到基线的一半。在计算量方面,SSD (6L)在256分辨率下仅需209.69 GFLOPs,而DDPM需要497.03 GFLOPs,计算量减少超过57%。在ImageNet-64这个更具挑战性的数据集上,SSD (4L)达到FID 17.89,与DDPM-ϵ的12.82相比仍有差距,但考虑到ImageNet包含1000个类别和130万张图像,这个结果表明SSD能够学习复杂的分布。时间扩展性分析显示,SSD的训练时间随分辨率增加而线性增长,而DDPM则呈二次增长,在256分辨率时SSD的优势最为显著。消融实验表明,Flexi-UNet在FID和推理速度上都略优于标准UNet架构,Lanczos采样的开销可以忽略不计(约0.03秒)。

DDPM、Blurring Diffusion和Scale Space Diffusion的公式对比
Table 1: DDPM、Blurring Diffusion和Scale Space Diffusion的公式对比
CelebA数据集上的主要结果
Table 2: CelebA数据集上的主要结果
ImageNet-64无条件图像生成结果
Table 3: ImageNet-64无条件图像生成结果
架构消融实验:Full UNet vs Flexi-UNet
Table 4: 架构消融实验:Full UNet vs Flexi-UNet
可视化生成样本和中间预测
Figure 4: 可视化生成样本和中间预测
时间扩展性:SSD vs DDPM的训练时间对比
Figure 6: 时间扩展性:SSD vs DDPM的训练时间对比
查看结构化数据
任务指标本文基线提升
CelebA-64 无条件图像生成 FID (50k samples, 1M iterations) SSD (2L): 2.14, SSD (3L): 3.61, SSD (4L): 4.28 DDPM-ϵ: 2.22, DDPM-x0: 2.98, Blurring Diffusion: 2.06 SSD (2L) 接近最优FID,训练时间减少11%(62.63 vs 70.30小时)
CelebA-128 无条件图像生成 FID (50k samples, 300k iterations) SSD (3L): 6.53, SSD (5L): 10.47 DDPM-ϵ: 4.16, DDPM-x0: 3.50 SSD (5L) 训练时间减少50%(25.41 vs 50.50小时)
CelebA-256 无条件图像生成 FID (50k samples, 300k iterations) SSD (3L): 7.79, SSD (4L): 10.52, SSD (6L): 13.50 DDPM-ϵ: 5.52, DDPM-x0: 5.47, Blurring Diffusion: 4.76 SSD (6L) 训练时间减少51%(42.88 vs 87.31小时),GFLOPs减少58%(209.69 vs 497.03)
ImageNet-64 无条件图像生成 FID (50k samples, 1M iterations) SSD (4L): 17.89 DDPM-ϵ: 12.82, DDPM-x0: 13.07 SSD在更具挑战性的数据集上仍能达到合理性能

局限与改进

尽管Scale Space Diffusion在效率上取得了显著提升,但存在几个重要局限性。首先,FID性能与DDPM基线相比存在一定差距,特别是在高分辨率设置下(CelebA-256的13.50 vs 5.52),这表明多尺度处理可能在某种程度上损失了生成质量。其次,ImageNet-64的结果(FID 17.89)相比DDPM基线(12.82)有明显差距,说明SSD在复杂多样化分布上的生成能力仍有待提升。第三,论文承认没有使用任何高级技术(如DDIM采样器、课程学习等)来优化性能,这可能限制了方法的最终表现。第四,实验仅在无条件图像生成任务上进行评估,未涉及条件生成、图像编辑等更广泛的应用场景。第五,论文使用1000步DDPM采样,虽然测试了250步和25步DDIM,但SSD的框架本身是为1000步设计的,减少采样步数可能导致性能下降。从个人观察来看,Flexi-UNet的设计虽然巧妙,但增加了架构复杂性,且不同分辨率下激活不同层可能导致参数利用不充分。

独立分析的弱点

基于对论文的深入分析,我识别出以下几个主要弱点:(1)**FID性能差距**:在高分辨率下,SSD的FID明显高于DDPM基线(例如CelebA-256上13.50 vs 5.52),这可能源于多尺度处理过程中的信息损失。改进方向:可以探索更精细的分辨率调度策略,或者在分辨率转换时引入残差连接来保留高频信息。(2)**ImageNet泛化能力不足**:SSD在ImageNet-64上FID为17.89,相比DDPM的12.82有显著差距,说明方法在复杂分布上的泛化能力有限。改进方向:可以结合类条件信息或使用更大的模型容量来处理ImageNet的多样性。(3)**架构设计过于复杂**:Flexi-UNet需要额外的1×1卷积层和动态层激活机制,增加了实现难度。改进方向:可以探索更简洁的多分辨率架构,如使用Neural Operator或自适应卷积。(4)**缺乏与其他多分辨率方法的直接比较**:论文没有与Cascaded Diffusion、Matryoshka Diffusion等方法进行公平比较,难以评估SSD的相对优势。改进方向:在相同计算预算下进行系统性比较。(5)**训练稳定性未讨论**:多分辨率训练可能导致梯度不稳定,但论文没有讨论训练过程中的稳定性问题。

未来方向

作者提出了几个有前景的研究方向:(1)**更先进的采样器**:使用DDIM、DPM-Solver等快速采样器替代DDPM采样器,可以同时提升推理速度和生成质量。(2)**课程学习**:借鉴Progressive GAN的思想,在训练初期使用低分辨率,逐步增加分辨率,可以改善训练优化。(3)**Transformer架构集成**:将Scale Space Diffusion与DiT、HDiT等Transformer架构结合,探索在像素空间和潜在空间的应用。基于SSD的成果,还可以延伸以下方向:(4)**潜在空间SSD**:将多分辨率降级作为潜在空间中的插值操作,实现更高效的Scale Space LDM。(5)**视频生成扩展**:将时间维度纳入尺度空间框架,实现时空联合的多尺度生成。(6)**条件生成**:将SSD扩展到文本到图像生成,利用低分辨率条件引导高分辨率生成。(7)**自适应分辨率调度**:根据输入图像的复杂度动态调整分辨率调度,而不是使用固定的调度函数。

复现评估

本文在复现性方面表现良好。代码和项目网站已公开,基于广泛使用的ADM代码库构建,降低了复现难度。所有实验细节都有详细记录:使用AdamW优化器,学习率1×10^{-4}(64×64和128×128)或5×10^{-5}(256×256),批量大小128或64,1000步线性噪声调度。数据集CelebA和ImageNet都是公开可用的标准数据集。然而,完全复现需要多GPU资源:实验在NVIDIA H100和RTX A4000上进行,256×256的实验需要2张H100。推理时间在单张NVIDIA GH200节点上测量。Lanczos算法和非各向同性噪声采样(Algorithm 3)的实现需要一定的数学背景,但论文提供了详细的伪代码。总体而言,对于有适当计算资源的研究者,该方法应该是可复现的。