---
title: Blind Deconvolution & Phase Retrieval
url: https://www.emergentmind.com/topics/blind-deconvolution-and-phase-retrieval-bdpr
type: topic
---

# Blind Deconvolution & Phase Retrieval

Blind Deconvolution and Phase Retrieval (BDPR) denotes a class of inverse problems in which both an unknown object and an unknown blur, kernel, probe, pupil, or illumination must be recovered from phaseless or intensity-only measurements. In the unified notation used for blind phase retrieval, the task is to find \((w,u)\) such that \( |\mathcal A(w,u)|^2=f \), which subsumes ptychography, Fourier ptychography, coherent diffractive imaging, FROG, and convolutional phase retrieval [2211.06619]. In this sense, BDPR combines the bilinear ambiguity of blind deconvolution with the phase loss of classical phase retrieval. Representative measurement models include phaseless Fourier measurements of circular convolution, \( \tilde y = |F(w\circledast x)| \) [1904.12680], low-resolution phaseless convolution measurements \( \boldsymbol y = |\boldsymbol F_{lo}(\boldsymbol x\circledast \boldsymbol h)|^2+\boldsymbol\eta \) [2111.13670], and incoherent turbulence imaging with coded aperture, \( g=f\otimes h_{\rm i} \) and \( h_{\rm i}=|\mathcal F^{-1}[P\times A]|^2 \) [2201.06267].

## 1. Core inverse-problem structure

BDPR problems are characterized by two coupled sources of nonconvexity. First, the forward map is bilinear or multiplicative in two unknowns. Second, the data are phaseless, so only magnitudes or intensities are observed. The survey formulation \( |\mathcal A(w,u)|^2=f \) makes this explicit and places a wide range of imaging models under a single algebraic template [2211.06619]. In convolutional settings, the model may be written as phaseless Fourier data of a circular convolution, \( y = |F(w\circledast x)| \), whereas in blind ptychography the phaseless data are \( \bm a_j = |\mathcal F(\omega\circ \mathcal S_j u)| \) for multiple scan positions \(j\) [1904.12680].

The role of additional physical structure varies by modality. In microscopy, origin-symmetric point spread functions imply that the blurred image and the original image share the same Fourier-transform phase [1503.04776]. In turbulence imaging with a coded aperture, the aberrated incoherent PSF is generated by a pupil field \(P\times A\) and the coded aperture acts as a known support mask in the pupil plane [2201.06267]. In low-resolution phaseless blind deconvolution, only low-frequency Fourier coefficients are retained, so the problem combines blind deconvolution, phase retrieval, and super-resolution in a single model [2111.13670].

| Setting | Representative measurement model | Structural ingredient |
|---|---|---|
| Convolutional BDPR | \( \tilde y=|F(w\circledast x)| \) | Known subspaces for both factors |
| Low-resolution BDPR | \( \boldsymbol y=|\boldsymbol F_{lo}(\boldsymbol x\circledast \boldsymbol h)|^2+\boldsymbol\eta \) | Partial Fourier sampling with \(m<n\) |
| Blind ptychography | \( \bm a_j=|\mathcal F(\omega\circ \mathcal S_j u)| \) | Multiple scan positions and overlap |
| Turbulence imaging | \( g=f\otimes |\mathcal F^{-1}[P\times A]|^2 \) | Coded aperture support on pupil plane |

This unification is not merely terminological. It indicates that many seemingly different imaging problems are governed by the same identifiability and optimization obstacles: bilinear factor ambiguity, lost phase, and strong sensitivity to side information, redundancy, or hardware design.

## 2. Ambiguities, identifiability, and structural priors

The most basic ambiguity in BDPR is bilinear scaling. In subspace-based convolutional formulations, \((h^\natural,m^\natural)\sim (\alpha h^\natural,\alpha^{-1}m^\natural)\), and the corresponding lifted matrices satisfy \((H^\natural,M^\natural)\sim(\alpha H^\natural,\alpha^{-1}M^\natural)\) [1904.12680]. Phase retrieval introduces additional global phase or sign ambiguities; for low-resolution phaseless blind deconvolution, the paper explicitly notes the inherent global phase ambiguity, while the tensorized formulation of simultaneous blind deconvolution and phase retrieval identifies the sign ambiguity \(x^\star\mapsto -x^\star\) [2111.13670].

Beyond these trivial ambiguities, one-dimensional linear convolution admits a richer algebraic ambiguity structure. In the \(z\)-domain, convolution becomes polynomial multiplication, and ambiguities are in one-to-one correspondence with factorization ambiguities generated by swapping zeros between the two factors [1701.04890]. Exact recovery from correlation data is therefore tied to a coprimeness condition: the two \(z\)-transforms must have no common factors. The same paper shows that adding both autocorrelations and cross-correlations turns blind deconvolution into a deterministic convex recovery problem, whereas omitting one autocorrelation destroys uniqueness [1701.04890].

Several BDPR approaches attack ambiguity through engineered priors or known context rather than through generic regularization alone. In microscopy, the convex set
\[
C_\phi=\{x[n_1,n_2]\mid \measuredangle X(w_1,w_2)=\measuredangle X_o(w_1,w_2)\}
\]
encodes the fact that symmetric blur preserves Fourier phase [1503.04776]. In turbulence imaging, the coded aperture \(A\) reduces the unknown pupil support to the transmitting coordinates \( \{u:A(u)=1\} \), which reduces the number of unknowns and improves conditioning [2201.06267]. In Fourier phase retrieval with background information, known background pixels are embedded into the reconstruction domain; for \(d\)-dimensional signals, the paper states that background lengths satisfying
\[
\prod_{i=1}^d(n_i+k_i)\ge 2\prod_{i=1}^d n_i+\prod_{i=1}^d(2n_i-1)
\]
are sufficient to ensure uniqueness with probability \(1\), assuming i.i.d. Gaussian background entries [2308.08328].

A common misconception is that BDPR is identifiable once enough measurements are taken. The literature summarized here does not support that claim. Instead, identifiability is repeatedly tied to concrete side conditions: random known subspaces, no common zeros, sufficient scan overlap, known background information, known probe magnitudes, or hardware support constraints [1808.05802].

## 3. Convex formulations and exact-recovery theory

One major line of BDPR research converts the original bilinear phaseless inverse problem into a convex program by lifting the unknown factors to positive semidefinite rank-one matrices. For phaseless Fourier measurements of a circular convolution, \(w=Bh\) and \(x=Cm\) are assumed to lie in known subspaces, and the lifted matrices \(H=hh^*\) and \(M=mm^*\) yield bilinear measurements in nonnegative scalar functionals. The nonconvex equalities are then relaxed to first-quadrant hyperbolic inequalities, and the rank constraints are replaced by trace minimization:
\[
\min_{H,M}\ \operatorname{Tr}(H)+\operatorname{Tr}(M)
\]
subject to
\[
\frac{1}{m}\langle b_\ell b_\ell^*,H\rangle\langle c_\ell c_\ell^*,M\rangle \ge y_\ell^2,\quad
H\succeq 0,\quad M\succeq 0.
\]
Under random Gaussian subspaces of dimensions \(k\) and \(n\), exact recovery up to scaling holds with high probability when \( m \gg (k+n)\log^2 m \), and the method does not require initialization [1904.12680].

A deterministic convex alternative arises when additional autocorrelation information is available. In the one-dimensional aperiodic convolution setting, blind deconvolution with the autocorrelations \(x_1*x_1^-\) and \(x_2*x_2^-\), together with the cross-correlations, can be written as a feasibility SDP over the lifted rank-one matrix \(X=xx^*\). The main theorem states that if the \(z\)-transforms of \(x_1\) and \(x_2\) have no common factors and \(x_i\in\mathbb C_{0,0}^{L_i}\), then the feasible convex program has the unique solution \(X^\#=xx^*\), hence exact recovery up to global phase [1701.04890].

Convexity also appears at the level of modular constraints rather than whole-program lifting. For microscopy deblurring, the closed convex set of images with prescribed Fourier phase is combined with the epigraph set of total variation,
\[
\text{C}_{ESTV}=\{\underline{\mathbf x}=[x^T\ z]^T\mid \mathrm{TV}(\mathbf x)\le z\},
\]
so that blind deconvolution alternates between blur updates, phase-consistency projection, TV-epigraph projection, and positivity or support constraints [1503.04776]. The paper argues that these projections are non-expansive and can be inserted into many blind deconvolution algorithms.

Convexification is therefore not a single method but a family of strategies. It may act globally, through lifted semidefinite recovery with exact-recovery theorems, or locally, through closed convex feasibility sets embedded inside alternating schemes. This suggests that the boundary between “convex BDPR” and “nonconvex BDPR” is often determined by which part of the original bilinear phase-loss structure is being relaxed.

## 4. Nonconvex algorithms, initialization, and optimization geometry

A second major line of work solves BDPR directly in factorized variables. In low-resolution phaseless blind deconvolution, BliPhaSu minimizes a non-convex least-squares objective over low-dimensional coefficients \( (\boldsymbol g,\boldsymbol z) \) after imposing known subspaces \( \boldsymbol h=\boldsymbol B\boldsymbol g \) and \( \boldsymbol x=\boldsymbol C\boldsymbol z \) [2111.13670]. The algorithm consists of a spectral initialization based on the leading eigenvectors of two moment matrices,
\[
\boldsymbol H_g=\frac{1}{m}\sum_{\ell=1}^m \boldsymbol y[\ell]\hat{\boldsymbol b}_\ell \hat{\boldsymbol b}_\ell^H,\qquad
\boldsymbol H_z=\frac{1}{m}\sum_{\ell=1}^m \boldsymbol y[\ell]\hat{\boldsymbol c}_\ell \hat{\boldsymbol c}_\ell^H,
\]
followed by stochastic Wirtinger-gradient iterations. Under \(m\ge \beta(k+s)\), bounded noise, and proper step sizes, the paper proves linear convergence in expectation to the true pair up to a noise floor [2111.13670].

For single-shot turbulence imaging, the reconstruction is based on the ptychographic iterative engine. The method alternates among updating the object spectrum \(\widehat F\), the incoherent transfer function \(\widehat H_{\rm i}\), and the pupil function \(\widehat P\) on the coded-aperture support. The amplitude constraint
\[
\widehat h_{\rm c}= \sqrt{\mathcal F^{-1}[\widehat H_{\rm i}]}\,
\exp\!\big(i\,\arg(\widehat h_{\rm c}')\big)
\]
and the phase-only pupil constraint
\[
\widehat P=\exp\!\big(i\,\arg(\widehat P')\big)
\]
enforce the non-absorbing turbulence model during iteration [2201.06267].

In blind ptychographic phase retrieval, generalized ADMM introduces split variables \(z=\mathcal A(\omega,u)\) and optimizes a Poisson-noise-motivated nonlinear least-squares or maximum-likelihood objective with amplitude constraints on both \(\omega\) and \(u\). The augmented-Lagrangian subproblems for \(\omega\) and \(u\) admit fast element-wise updates because of the bilinear structure and orthogonality of the Fourier transform, while the \(z\)-update is a proximal step for penalized AGM or pIPM data terms [1808.05802]. The same paper proves convergence of the generalized ADMM to a stationary point under sufficient-overlap and bounded-preconditioner conditions.

Several theoretical works analyze why such nonconvex methods can nevertheless behave regularly. For phase retrieval and blind deconvolution, vanilla gradient descent can remain inside a region of incoherence and contraction without explicit regularization, a phenomenon called implicit regularization [1711.10467]. For the tensorized surrogate of simultaneous blind deconvolution and phase retrieval, the population-risk landscape on the unit sphere has global minimizers at \(x=\pm x^\star,\ h=h^\star\), and Riemannian gradient descent converges linearly under local initialization conditions; analogous local linear convergence and noise-floor guarantees are proved under a tensor restricted isometry property for the surrogate sensing model [2509.10834].

These results do not imply that generic nonconvex BDPR is globally benign. They show, rather, that benign geometry can emerge under very specific models: random subspaces, sufficient overlap, proper initialization, or TRIP-like sensing conditions. The survey literature treats convergence for many practical projection methods as still incomplete [2211.06619].

## 5. Experimental modalities and domain-specific formulations

BDPR appears in substantially different physical regimes, and the literature repeatedly adapts the inverse problem to the governing optics. In microscopy, the key assumption is that many point spread functions are approximately origin-symmetric, so the Fourier phase of the object is preserved and can be used as a convex constraint inside a blind deconvolution loop [1503.04776]. In turbulence imaging under spatially incoherent light, the object and the turbulence-induced pupil aberration are estimated jointly from a single captured image, with a binary random coded aperture on the pupil plane acting as a support prior [2201.06267].

Ptychographic BDPR treats the sample and the probe jointly as unknowns. The survey formulation \( f_j^{P}=|\mathcal F(w\circ \mathcal S_j u)|^2 \) unifies conventional ptychography with related models such as Fourier ptychography and FROG [2211.06619]. One blind-pycthography line formulates the task directly as constrained nonlinear optimization under Poisson noise and solves it with convergent ADMM [1808.05802]. Another line rewrites far-field blind ptychography as a blind deconvolution problem over autocorrelations,
\[
(\mathbf Y^T\mathbf F_d^T)_k
=
d\cdot (x\circ S_k\bar x)*_d(\tilde m\circ S_{-k}\bar{\tilde m})+\text{noise},
\]
then combines blind deconvolution, angular synchronization, and regularized Wirtinger descent [2305.04060]. This algebraic recasting makes the ptychography-to-BDPR connection explicit.

Phase-retrieval imaging can also be blind in a parametric sense. In the regularization view of Paganin-type phase retrieval, the recovered quantity \(p\) is obtained from the convolutional filter
\[
\widehat{\,_\ell}(\mathbf q)=\frac{1}{1+4\pi^2\ell\|\mathbf q\|^2},
\]
where \(\ell=d\delta/\mu\) depends on physical sample parameters and propagation geometry [1702.05092]. The paper treats the unknown \(\ell\) as the blur parameter of a Tikhonov-regularized deconvolution problem and determines it by a one-dimensional optimization based on L-curve curvature. This suggests a broader BDPR interpretation in which the “blind” variable is not necessarily a free-form kernel but may be a physically meaningful parameter controlling a convolution model.

Learned priors have introduced a different kind of structure. Deep PBD assumes that the sharp image and blur kernel lie in the ranges of two pretrained generators and minimizes a phaseless measurement-misfit in latent space with alternating gradient descent and latent \(\ell_2\) regularization [2002.12578]. The paper reports performance on MNIST, CelebA, and Shoes, under both oversampled Fourier and subsampled Fourier ptychography forward models. This replaces analytical priors such as support, TV, or symmetry with class-specific generator manifolds.

## 6. Limitations, misconceptions, and open directions

A persistent misconception is that BDPR is simply “blind deblurring plus phase retrieval.” The literature indicates a more heterogeneous picture. Some problems are blind because both factors are unknown vectors in known random subspaces; some are blind because an unknown probe modulates an unknown object; some are blind because a physical blur parameter such as \(\ell\) is unknown; and some are made identifiable only by additional autocorrelations, background information, or coded apertures [1702.05092]. The commonality lies in joint recovery under phaseless data, not in a single universal forward model.

Another misconception is that convex theory has resolved BDPR in general. The strongest exact-recovery results rely on restrictive side conditions: random Gaussian subspaces with \( m\gg (k+n)\log^2 m \), deterministic correlation designs with no common factors, or sufficient known background information [1904.12680]. The survey literature explicitly notes that uniqueness and convergence theory remain subtle for realistic blind ptychographic and Fourier ptychographic models, especially under noise, probe drift, partial coherence, or scan pathologies [2211.06619].

Hardware assistance improves conditioning but introduces its own tuning problems. In coded-aperture single-shot blind deconvolution, successful experimental reconstruction occurs for transmittance ratios between \(30\%\) and \(60\%\); lower transmittance leads to photon-starved noisy measurements, whereas higher transmittance weakens the support constraint [2201.06267]. In background-assisted phase retrieval, decreased background length improves practicality, but the theory and algorithms still depend on sufficient known context and on the interaction between measurement noise and background perturbation [2308.08328]. In blind ptychography, extra phaseless measurements of the probe can suppress periodic artifacts, but obtaining them requires additional calibration [1808.05802].

The main open directions identified in the survey include stronger convergence theory for blind problems, sharper identifiability results beyond idealized random or subspace assumptions, more noise-robust formulations, scalable algorithms for large 2D and 3D imaging, and better handling of experimental artifacts such as probe drift, background scattering, and partial coherence [2211.06619]. The recent tensor-landscape analysis suggests one possible path: exposing structured low-rank tensor geometry and proving local linear convergence for Riemannian algorithms under tractable surrogate models [2509.10834]. A plausible implication is that future BDPR theory will continue to alternate between physically faithful formulations and analytically simpler surrogates, with hardware design, side information, and optimization geometry treated as inseparable parts of the same inverse problem.

Source: https://www.emergentmind.com/topics/blind-deconvolution-and-phase-retrieval-bdpr