---
title: 'Blind Denoising Diffusion Models: Theory & Applications'
url: https://www.emergentmind.com/topics/blind-denoising-diffusion-models-bddms
type: topic
---

# Blind Denoising Diffusion Models: Theory & Applications

Searching arXiv for recent papers on Blind Denoising Diffusion Models and closely related blind diffusion formulations.
Blind Denoising Diffusion Models (BDDMs) are diffusion-based models for settings in which a critical nuisance variable is unavailable at training time, sampling time, or both. In the narrow theoretical sense, a BDDM is a generative diffusion model whose denoiser is *blind* to the noise amplitude $\sigma$, so the network $f_\theta:\mathbb R^d\to\mathbb R^d$ never receives $\sigma$ as input yet is trained by marginalizing over $\sigma$ and can still follow an implicit reverse-time noise schedule [2602.09639]. In a broader applied sense found across recent arXiv literature, the term encompasses diffusion formulations for blind denoising and blind inverse problems in which the clean image, degradation parameters, or measurement operator are unknown, and the diffusion process is conditioned instead on auxiliary structure such as blind-spot predictions, degradation representations, or learned operator priors [2504.00264]. This dual usage makes BDDMs both a theoretical object in schedule-free generative modeling and a practical design pattern for ill-posed restoration problems.

## 1. Terminological scope and problem settings

Within the cited literature, *blindness* refers to several distinct unknowns. In the theoretical formulation of “Blind denoising diffusion models and the blessings of dimensionality” [2602.09639], the blind variable is the noise amplitude: the denoiser never sees $\sigma$ in either training or sampling. In medical self-supervised denoising, the blind variable is the unknown corruption process in observations of the form $x' = x + n$, where $n$ may come from an unknown, possibly non-i.i.d. distribution such as signal-dependent Gaussian, Poisson, or spatially correlated noise [2504.00264]. In blind inverse problems, both the target image and the measurement operator can be unknown, as in $y = H_\theta(x) + \varepsilon$ with unknown $\theta$ [2505.22923], or $y = Ax + z$ with unknown linear operator $A$ [2301.12686]. In blind face restoration, the unknown is the degradation type and degree affecting a low-quality face image [2411.10508]. In blind super-resolution for remote sensing images, the challenge is the unknown degradation kernel [2305.12170].

This suggests that BDDMs are better understood as a family of diffusion methods for *latent nuisance inference under uncertainty* rather than a single standardized architecture. A plausible implication is that comparisons across papers require care, because “blind” may mean unknown noise amplitude, unknown corruption law, unknown kernel, or unknown operator, and the conditioning strategy changes accordingly.

A compact view of the problem classes represented in the literature is given below.

| Setting | Unknown quantity | Representative paper |
|---|---|---|
| Schedule-free generative diffusion | Noise amplitude $\sigma$ | [2602.09639] |
| Self-supervised medical denoising | Clean reference and noise law | [2504.00264] |
| Blind deblurring / inverse problems | Measurement operator parameters | [2505.22923], [2301.12686] |
| Blind restoration | Degradation representation | [2411.10508] |
| Blind super-resolution | Blur kernel | [2305.12170] |
| Feature-space blind denoising | Noise in transferred quality-aware features | [2401.11949] |

## 2. Core mathematical formulation of blind denoising diffusion

The canonical theoretical BDDM replaces the standard denoiser $b_\theta:\mathbb R^d\times\mathbb R_+\to\mathbb R^d$, which is trained with explicit access to the true noise level $\sigma$, by a blind denoiser $f_\theta:\mathbb R^d\to\mathbb R^d$ trained by marginalizing over $\sigma$ drawn from a prior $\Theta$ on $[\sigma_T,\sigma_0]$ [2602.09639]. The population objective is
$$
f^* = \arg\min_f \mathbb E_{X\sim p_X,\;\sigma\sim\Theta,\;z\sim N(0,I)} \|X - f(X+\sigma z)\|^2.
$$
Proposition 3.1 in that work shows that the population minimizer admits the Bayesian form
$$
f^*(y) = y + \int \sigma^2 \nabla \log p_\sigma(y)\,d\mu(\sigma|y),
$$
so the blind residual
$$
s^*(y):=f^*(y)-y
$$
is a mixture of non-blind score functions [2602.09639]. This is the defining mathematical distinction: the score surrogate is not indexed by a provided noise amplitude, but by a posterior over amplitudes induced by the noisy sample itself.

The forward process in that paper is a variance-exploding SDE,
$$
dX_t = \sqrt{2a_t}\,dB_t,\qquad X_0\sim p_X,
$$
with $\mathrm{Law}(X_t)=p_{\sigma_t}=p_X * \mathcal N(0,\sigma_t^2 I)$ and $\sigma_t^2=2\int_0^t a_s\,ds$ [2602.09639]. The blind reverse process replaces the true $\sigma$-conditioned score by the blind estimator $s_\theta(x)=f_\theta(x)-x$ and uses
$$
dY_t = s_\theta(Y_t)\,dt + \sqrt{2a_t}\,dB_t.
$$
No $\sigma$ input is required at sampling time [2602.09639].

This differs materially from conventional DDPM and EDM parameterizations in the applied papers, which usually retain timestep conditioning $t$ and sometimes additional blind-side information. For example, DiffDenoise adopts the standard DDPM forward process
$$
q(X_t|X_{t-1}) = \mathrm{Normal}(X_t; \sqrt{\alpha_t}X_{t-1}, \beta_t I),
$$
with reverse kernel
$$
p_\theta(X_{t-1}|X_t,x_c)=\mathrm{Normal}(X_{t-1};\mu_\theta(X_t,x_c,t),\tilde\beta_t I),
$$
where $x_c=B(x')$ is the blind-spot conditioning image [2504.00264]. Blind-Spot Guided Diffusion uses two reverse kernels, one unconditional branch and one blind-spot branch, and fuses their noise estimates in a classifier-free-style guidance [2509.16091]. DR-BFR works in latent space and conditions on both low-quality content features and a degradation representation learned by a Degradation Representation Module [2411.10508].

## 3. Implicit schedules, dimensionality, and theoretical guarantees

The central theoretical claim of the strict BDDM formulation is that blind denoisers automatically track an implicit noise schedule along the reverse process [2602.09639]. Under the heuristic ansatz that $\mathrm{Law}(Y_t)\approx p_{\sigma_t}$ and replacing the posterior on $\sigma$ by $\delta_{\sigma_t}$, the ideal SDE becomes
$$
dX_t = \sigma_t^2\nabla\log p_{\sigma_t}(X_t)\,dt + \sqrt{2a_t}\,dB_t.
$$
A Fokker–Planck calculation shows that consistency $\mathrm{Law}(X_t)=p_{\sigma_t}$ holds exactly if
$$
\tfrac12\,d(\sigma_t^2)/dt = -\sigma_t^2 + a_t,
$$
equivalently
$$
\sigma_t^2 = \sigma_0^2 e^{-2t} + 2\int_0^t a_s e^{-2(t-s)}\,ds
$$
[2602.09639]. Lemma 3.3 and 3.4 in that paper show that $\sigma_t$ decreases to $0$ if $a_t\to0$.

The main rigorous error bound depends on bounded support, low intrinsic dimension $k\ll d$, and a score-estimation error assumption. Theorem 3.6 shows
$$
\mathrm{KL}(p_{\sigma_T}\|\mathrm{Law}(Y_T))
\lesssim R^2/\sigma_0^2 + \epsilon_{\mathrm{score}}^2 + ((R^2+k)k^2/d)\int_0^T a_t^{-1}dt
$$
[2602.09639]. The last term comes from Bayesian noise-level estimation error, and the paper’s posterior concentration argument yields
$$
\mathbb E[|1/\sigma_t^2-1/\sigma^2|^2] \lesssim (1/d + k^2/d^2)
$$
[2602.09639]. The phrase “blessings of dimensionality” refers precisely to this regime: when the ambient dimension is sufficiently large relative to intrinsic dimension, blind estimation of the residual noise level becomes accurate.

For the special schedule $a_t=a\,\sigma_t^2$ with $0<a<1$, the paper derives $\sigma_t=\sigma_0e^{-(1-a)t}$ and gives a discretization result for the exponential–Euler integrator. Theorem 3.11 states
$$
\mathrm{Disc}(h)\lesssim_{\log} R^2k^3 h^2 A_T\cdot\mathbf 1_{a\neq 1/2} + k^2 h,
$$
with $A_T\approx 1/(2a(1-a)\sigma_T^2)$, and notes that if $a=1/2$ the $h^2$ term vanishes [2602.09639]. Corollary 3.14 then gives a sampling guarantee in the perceptual metric $D_X$: with
$$
\sigma_T\sim \epsilon/\sqrt{k},\qquad \sigma_0\sim R/\epsilon,\qquad h\sim \epsilon^2/k^2,\qquad N=T/h\sim k^2/\epsilon^2,
$$
one obtains $D_X(\mathrm{Law}(Y_N),p_X)\lesssim \epsilon$ provided $d\gg (R^2+k)k^3\epsilon^{-4}$ up to logs [2602.09639].

A common misconception is that omitting the explicit noise amplitude necessarily destroys control of the reverse process. The theoretical analysis argues the opposite under low-intrinsic-dimensional structure: the blind denoiser can infer the residual noise variance from the noisy image itself and thereby recover an implicit schedule [2602.09639].

## 4. Conditioning mechanisms in practical BDDMs

Applied BDDMs rarely remain fully unconditional. Instead, they compensate for blind uncertainty by injecting auxiliary signals that are informative about structure, degradation, or operator state.

DiffDenoise conditions a DDPM-style denoiser on the output of a pretrained Blind-Spot Network, $x_c=B(x')$, obtained under $J$-invariance [2504.00264]. The conditioning image is concatenated channel-wise with $X_t$ and injected at multiple resolutions via simple concatenation in the encoder and decoder of a U-Net–style backbone; the paper notes that cross-attention or FiLM layers are also possible, but concatenation sufficed in its experiments [2504.00264]. The training objective is the standard simple denoising loss
$$
L_{\mathrm{simple}}(\theta)=\mathbb E_{x',x_c,t,\epsilon}\|\epsilon-\epsilon_\theta(X_t,x_c,t)\|^2,
$$
with $X_t=\sqrt{\bar\alpha_t}x' + \sqrt{1-\bar\alpha_t}\epsilon$ [2504.00264].

Blind-Spot Guided Diffusion likewise combines blind-spot structure with a conventional diffusion branch. It learns two Gaussian reverse kernels, an unconditional branch $p_\theta^u(x_{t-1}|x_t)$ and a blind-spot branch $p_\theta^b(x_{t-1}|x_t,x_0^c)$, and fuses their noise predictions as
$$
\epsilon(x_t)=w\cdot \epsilon_{fb}(x_t,x_0^c|t) + (1-w)\cdot \epsilon_{fu}(x_t),\qquad w\in[0,1]
$$
[2509.16091]. The blind-spot branch is based on the PUCA blind-spot network modified to be time-dependent, while the non-blind branch is a standard U-shaped architecture as in Ho et al. [2509.16091].

DR-BFR uses a different conditioning logic. Its Degradation Representation Module extracts a degradation prompt $f_{DR}$ from low-quality faces via unsupervised contrastive learning, LQ reconstruction, and a distribution-matching loss, with total loss
$$
L_{DR}=L_{\mathrm{recon}}+\lambda_1L_{\mathrm{contrast}}+\lambda_2L_{\mathrm{dist}}
$$
[2411.10508]. The Latent Diffusion Restoration Module then conditions on both latent low-quality content features $f_{LQ}$ and the degradation representation $f_{DR}$. The former is injected by simple channel-wise concatenation at the U-Net input, while the latter is injected via cross-attention layers, with time-adaptive prompt strength $\phi_t(f_{DR})=\alpha_t\cdot f_{DR}+\beta_t$ where $(\alpha_t,\beta_t)=\mathrm{MLP}(t)$ [2411.10508].

Blind super-resolution in DDSR splits conditioning across two DDPMs: a kernel-estimation DDPM conditioned on an LR encoding $u=f_e(I_{LR})$, and a reconstruction DDPM conditioned on both $u$ and the predicted kernel $v$ [2305.12170]. Blind inverse-problem samplers use yet another strategy. Blind-PnPDM alternates between image and operator updates, each formulated as a Gaussian denoising task solved by a separate diffusion prior, one for clean images and one for operator parameters [2505.22923]. GibbsDDRM keeps a pretrained unconditional diffusion model for the signal prior and couples it with explicit online updates of the operator through a partially collapsed Gibbs sampler and Langevin dynamics [2301.12686].

Taken together, these papers suggest that practical BDDMs often replace direct access to the unknown nuisance with a *surrogate observable*: a blind-spot estimate, a latent degradation code, a current operator sample, or a pseudo-clean conditioning image.

## 5. Sampling, stabilization, and blind posterior inference

Sampling behavior is a major differentiator within the BDDM literature. The theoretical schedule-free BDDM trains a blind denoiser with Algorithm 1 and samples with Algorithm 2 by initializing $Y_0\sim N(0,\sigma_0^2I)$, iterating
$$
s_k=f_\theta(Y_k)-Y_k,\qquad
Y_{k+1}=Y_k + h\,s_k + \xi_k,
$$
with $\xi_k\sim N(0,2\int_{kh}^{(k+1)h}a_t\,dt\cdot I)$, and stopping when $\|s_k\|/\sqrt d \le \sigma_T$ [2602.09639]. The stopping rule itself reflects the model’s ability to estimate residual noise from the sample.

DiffDenoise uses DDIM sampling for speed, but introduces Symmetric-Noise Stabilization (SRDS). Two DDIM trajectories are run from initial noises $\epsilon$ and $-\epsilon$:
$$
\hat x^+ = \mathrm{DDIM}(\epsilon; x_c),\qquad
\hat x^- = \mathrm{DDIM}(-\epsilon; x_c),
$$
and the final output is averaged,
$$
\hat x=(\hat x^+ + \hat x^-)/2
$$
[2504.00264]. The stated rationale is that score-model errors are roughly odd functions in $\epsilon$, so averaging cancels leading-order bias terms [2504.00264].

Blind-Spot Guided Diffusion uses Random & Complementary Replacement during sampling. At each step it forms
$$
x_0^c|t = (x_t-\sqrt{1-\alpha_t}\,\epsilon(x_t))/\sqrt{\alpha_t},
$$
replaces pixels of $x_0^c|t$ by $x_0$ with probability $p$ in Base Replacement, repeats sampling for $R$ rounds, and averages the resulting estimates [2509.16091]. The paper states that random and complementary pixel replacement reintroduces noisy-pixel information in controlled fashion, preventing over-smoothing and preserving spatial continuity [2509.16091].

Blind inverse-problem samplers use alternating posterior updates rather than a single reverse chain. Blind-PnPDM employs a split-Gibbs sampler: an image-update alternates a Gaussian likelihood step with a prior denoising step implemented by an EDM image denoiser $D_x$, and an operator-update alternates an analogous likelihood step with a prior denoising step implemented by a second diffusion denoiser $D_\theta$ [2505.22923]. GibbsDDRM similarly samples from $p(x_{0:T},A|y)$ by interleaving modified DDRM updates for the latent signal states with Langevin updates for the unknown operator, using SVD spectral coordinates for efficient Gaussian conditionals [2301.12686].

A plausible implication is that “sampling” in BDDMs spans three distinct computational idioms: schedule-free reverse SDE integration [2602.09639], guided diffusion with stabilization heuristics [2504.00264; 2509.16091], and MCMC-style alternating posterior simulation for blind inverse problems [2505.22923; 2301.12686].

## 6. Empirical domains and reported performance

The empirical record of BDDMs is diverse because the blind variable differs by domain.

In self-supervised medical image denoising, DiffDenoise was tested on FastMRI knee MRI with Gaussian, Poisson, Gamma pixel noise, SIIM chest X-rays, spatially correlated variants, and real low-field brain MRI (M4Raw) [2504.00264]. The paper reports that on independent Gaussian pixel noise, DiffDenoise achieves $+4\,\mathrm{dB}$ PSNR gain over the corrupted input and sits only $1\,\mathrm{dB}$ below a fully supervised NAF-Net while using no clean training images, and that SSIM matches or exceeds supervised baselines [2504.00264]. Against Neighbor2Neighbor, PUCA, LG-BPN, and Noise2Self, it is reported as $+1$–$2\,\mathrm{dB}$ PSNR better and as preserving vessel and cartilage edges that other methods blur [2504.00264]. Under spatially correlated noise, competing BSN-based approaches degrade by $1$–$2\,\mathrm{dB}$ as noise correlation grows, whereas DiffDenoise stays within $0.2\,\mathrm{dB}$ [2504.00264]. On real M4Raw brain MRI, it outperforms PUCA by $0.3$–$0.5\,\mathrm{dB}$ PSNR and surpasses NAF-Net in SSIM [2504.00264].

In real-world image denoising, Blind-Spot Guided Diffusion was evaluated on SIDD and DND. On SIDD, the paper reports 37.98 dB / 0.944, compared with approximately 37.69 dB / 0.937 for the best prior BSN, SelfFormer; on DND it reports 38.99 dB / 0.943, compared with approximately 38.92 dB / 0.943 for the best prior method [2509.16091]. Ablation results indicate guidance weight $w$ performs best around $0.7$–$0.8$ for PSNR, replacement probability $p$ is optimal at approximately $0.25$, sampling steps are best at $8$–$16$, and $8$ rounds $\times$ $8$ steps are recommended [2509.16091].

In blind face restoration, DR-BFR reports on CelebA-Test that DR-BFR attains FID $=13.69$ and NIQE $=5.01$, both best, with LPIPS $=0.2499$ and PSNR $\approx24.55$, SSIM $\approx0.643$ [2411.10508]. The paper states that the next best SOTA methods, CodeFormer, RestoreFormer, and DiffBIR, obtain FID approximately $16$–$19$ and NIQE approximately $5.3$–$6.1$, and that on real-world sets DR-BFR consistently attains top-2 NIQE and FID [2411.10508]. Its ablation indicates that removing degradation representation worsens FID and NIQE, with DR-None at FID $\approx20.1$, NIQE $\approx5.89$, and the full DR-ALL at FID $\approx13.7$, NIQE $\approx5.01$ [2411.10508].

In blind deblurring, Blind-PnPDM was tested on 100 FFHQ images under Gaussian kernels with random sigma in $[1,5]$ px and motion kernels with random length in $[5,25]$ px [2505.22923]. For Gaussian blur it reports 27.13 dB PSNR, 0.802 SSIM, and 0.180 LPIPS; for motion blur it reports 27.42 dB PSNR, 0.795 SSIM, and 0.176 LPIPS [2505.22923]. These exceed the listed baselines Pan-DCP, DeblurGAN-v2, BlindDPS, and GibbsDDRM on those test conditions [2505.22923]. GibbsDDRM itself, on FFHQ $256\times256$ blind image deblurring with random motion kernels and additive Gaussian noise $\sigma_y=0.02$, reports FID 38.71, LPIPS 0.115, and PSNR 25.80 [2301.12686]. On vocal dereverberation it reports FAD 4.21, SI-SDR improvement $+0.59$, and SRMR 8.40 [2301.12686].

In ultrasound denoising, the DDPM-based unsupervised method of “Deep Ultrasound Denoising Using Diffusion Probabilistic Models” reports at $T=20$ a PSNR of approximately 23.6 dB versus 22.0 dB for NLM and 22.4 dB for BM3D, with GCNR approximately 92.6% versus 91.6% and 90.4%, while preserving speckles [2306.07440]. In blind BIQA, PFD-IQA reports PLCC 0.937 on KADID and 0.922 on LIVE-C, with overall average SRCC/PLCC of 0.931/0.935 versus previous SOTA 0.905/0.887 [2401.11949].

## 7. Conceptual implications, misconceptions, and open directions

One recurring misconception is to equate BDDMs exclusively with blind *image denoising*. The cited literature shows a broader scope: BDDMs appear in generative sampling without explicit noise schedules [2602.09639], self-supervised denoising [2504.00264; 2509.16091; 2306.07440], blind restoration [2411.10508], blind super-resolution [2305.12170], blind inverse problems [2505.22923; 2301.12686], and even feature-space refinement for BIQA [2401.11949]. This suggests that the essential ingredient is not denoising in the narrow pixel-space sense, but diffusion-based inference when the corruption descriptor is latent or partially observed.

A second misconception is that blind diffusion must sacrifice fidelity because it lacks explicit conditioning on the true degradation amplitude or operator. The theoretical work argues that schedule-free blind denoisers can outperform non-blind counterparts by correcting mismatch between the true residual noise of the image and the noise assumed by a prescribed schedule [2602.09639]. The empirical papers make analogous domain-specific claims: blind-spot guidance is introduced because pure BSNs sacrifice local detail and pure diffusion branches tend to blur [2509.16091]; degradation prompts in DR-BFR are introduced because diffusion models otherwise lack awareness of specific degradation and may produce unnatural details and inaccurate textures [2411.10508]; DiffDenoise uses blind-spot conditioning and symmetric-noise averaging to preserve high-frequency structures that self-supervised methods tend to over-smooth [2504.00264].

The limitations are equally consistent. The theoretical guarantees in [2602.09639] require low intrinsic dimension and large ambient dimension, and the paper notes that quantifying intrinsic dimension for natural images remains an empirical challenge. Blind-PnPDM highlights higher per-sample cost because each outer iteration requires two diffusion runs and depends on well-trained diffusion models [2505.22923]. GibbsDDRM identifies SVD cost, sampling cost, and non-convexity of the joint posterior as practical difficulties [2301.12686]. DR-BFR notes that unseen corruptions such as watermarks or real compression artifacts may not be fully captured because its Degradation Representation Module is trained on synthetic degradations [2411.10508].

The trajectory of the field points toward richer blind conditioning and broader task coverage. The papers explicitly mention integrating blind denoisers into latent or text-conditioned models [2602.09639], extending blind diffusion to deblurring, inpainting, and other inverse problems with unknown or spatially varying $\sigma$ [2602.09639], retraining degradation representation modules for deblurring, denoising, and super-resolution [2411.10508], and applying blind plug-and-play diffusion to blind compressive sensing or unified joint annealing schemes [2505.22923]. Taken together, these works depict BDDMs not as a fixed model class but as a rapidly diversifying methodology for diffusion-based inference under hidden corruption structure.

Source: https://www.emergentmind.com/topics/blind-denoising-diffusion-models-bddms