---
title: 'PICCS: Prior Image Constrained Compressed Sensing'
url: https://www.emergentmind.com/topics/prior-image-constrained-compressed-sensing-piccs
type: topic
---

# PICCS: Prior Image Constrained Compressed Sensing

Searching arXiv for PICCS and closely related prior-information compressed sensing work to ground the article in cited papers.
Prior Image Constrained Compressed Sensing (PICCS) is a prior-based compressed sensing framework for reconstructing an unknown target signal or image from undersampled measurements while exploiting a known, similar signal, typically a prior image. In the synthesis setting, it fits exactly into the general problem of compressed sensing with prior information, in which one reconstructs \(x\) from \(y = Ax\) while incorporating a side signal \(w\) through objectives such as \(\|x\|_{1} + \lambda \|x-w\|_{1}\); in CT and CBCT, the same idea is commonly instantiated in the analysis or total variation (TV) domain as a joint penalty on the image and on its difference from a prior image [1408.5250]. In repeated scanning procedures, especially under-sampled CBCT with only local changes relative to an earlier high-quality scan, PICCS uses that prior image to suppress undersampling artifacts and preserve genuine changes such as a new tumor or a surgical tool [2509.08574].

## 1. Definition and problem setting

PICCS addresses inverse problems of the form
\[
y = A x^{\star}
\]
in the noiseless case, or
\[
\mathbf{A}\boldsymbol{x} + \boldsymbol{e} = \boldsymbol{b}
\]
in the noisy CBCT formulation, where \(A\) is the sensing or system matrix, \(x^{\star}\) is the true target image, and \(w\) or \(x_p\) is a known similar image available beforehand [1408.5250]. The defining feature is that reconstruction is constrained not only by sparsity of the target but also by sparsity of the deviation from the prior.

In the formulation analyzed as compressed sensing with prior information, the target is sparse or analysis-sparse, and the prior is informative when it is similar to the target in support, sign pattern, or transform domain [1408.5250]. In CBCT repeated scanning, the relevant regime is one in which a sequence of under-sampled scans is acquired on the same object, the anatomy is largely unchanged, and only local and sparse changes occur. The assumed prior image is a good initial reconstruction from a previous over-sampled or high-dose scan and is well aligned to the current geometry [2509.08574].

PICCS is therefore not merely a denoising prior around a reference image. Its canonical form enforces sparsity both of the current image and of the change relative to the prior. In CT practice, this is typically implemented with TV penalties, so that genuine, salient differences are retained while spurious differences induced by undersampling are discouraged [2509.08574].

## 2. Convex formulations

In the synthesis setting, the central noiseless convex programs are the \(L_{1}\)-\(L_{1}\) and \(L_{1}\)-\(L_{2}\) formulations:
\[
\min_{x}\; \|x\|_{1} + \lambda \|x-w\|_{1}
\quad \text{s.t.} \quad Ax = y,
\]
and
\[
\min_{x}\; \|x\|_{1} + \lambda \|x-w\|_{2}
\quad \text{s.t.} \quad Ax = y,
\]
or, in the quadratic prior form that is analyzed,
\[
\min_{x}\; \|x\|_{1} + \frac{\beta}{2}\|x-w\|_{2}^{2}
\quad \text{s.t.} \quad Ax = y.
\]
The noisy constrained and Lagrangian versions replace the exact data-consistency constraint by either \(\|Ax-y\|_{2}\le \varepsilon\) or a quadratic data term \(\tfrac12\|Ax-y\|_{2}^{2}\) plus regularization [1408.5250].

PICCS in CT is the analysis-domain counterpart of the \(L_{1}\)-\(L_{1}\) template. In the notation used for analysis or TV regularization,
\[
\min_{x}\; \alpha \|\Psi x\|_{1} + \beta \|\Psi(x-w)\|_{1}
\quad \text{s.t.} \quad \|Ax-y\|_{2}\le \varepsilon,
\]
where \(\Psi\) is an analysis operator such as the discrete gradient for TV, \(A\) is the CT system matrix, and \(w\) is the prior image [1408.5250].

The 2025 CBCT implementation uses the unconstrained quadratic form
\[
\min_{\boldsymbol{x}} \;\; \|\mathbf{A}\boldsymbol{x} - \boldsymbol{b}\|_2^2
\;+\; \alpha^2\,\mathrm{TV}(\boldsymbol{x})
\;+\; \lambda^2\,\mathrm{TV}(\boldsymbol{x}-\boldsymbol{x}_p),
\]
with 3D isotropic discrete TV
\[
\mathrm{TV}(\boldsymbol{x})
=
\left\|
\sqrt{(\mathbf{D}_x \boldsymbol{x})^2 + (\mathbf{D}_y \boldsymbol{x})^2 + (\mathbf{D}_z \boldsymbol{x})^2}
\right\|_1.
\]
That work also evaluates a PIPLE-style variant,
\[
\min_{\boldsymbol{x}} \;\; \|\mathbf{A}\boldsymbol{x} - \boldsymbol{b}\|_2^2
\;+\; \alpha^2\,\mathrm{TV}(\boldsymbol{x})
\;+\; \lambda^2\,\|\boldsymbol{x} - \boldsymbol{x}_p\|_2^2,
\]
which replaces the TV penalty on the difference by an \(\ell_2\) prior term [2509.08574].

Several implementation details delimit the scope of that CBCT formulation. No explicit non-negativity constraint is enforced, no spatial masks or adaptive weighting are used, and alignment of the prior image is assumed adequate. Parameters \(\alpha\) and \(\lambda\) are treated as given and tuned empirically to balance artifact suppression and change preservation [2509.08574].

## 3. Geometric recovery theory and measurement bounds

The theoretical analysis of PICCS-type reconstruction is geometric. For a convex objective \(f\), the descent cone at \(x^{\star}\) is
\[
\mathcal{D}(f,x^{\star}) := \mathrm{cone}\big(S_f(x^{\star})-x^{\star}\big),
\]
where \(S_f(x^{\star})=\{x: f(x)\le f(x^{\star})\}\). Exact recovery occurs when
\[
\mathcal{D}(f,x^{\star}) \cap \mathrm{null}(A) = \{0\}.
\]
For Gaussian sensing matrices, the probability of success is controlled by the Gaussian width
\[
w(C) := \mathbb{E}\big[\sup_{z\in C\cap \mathbb{S}^{n-1}} g^{\top}z\big]
= \mathbb{E}\big[\mathrm{dist}(g,C^{\circ})\big],
\]
with \(g\sim \mathcal{N}(0,I_n)\) [1408.5250].

The general high-probability recovery condition is
\[
m \ge w^{2}\big(\mathcal{D}(f,x^{\star})\big) + 1.
\]
In the noisy constrained case, if
\[
m \ge \frac{w^{2}\big(\mathcal{D}(f,x^{\star})\big)+3/2}{(1-\epsilon)^{2}},
\]
then every solution \(\hat{x}\) satisfies
\[
\|\hat{x}-x^{\star}\|_{2}\le \frac{2\sigma}{\epsilon}
\]
with high probability when \(\|Ax^{\star}-y\|_{2}\le \sigma\) [1408.5250].

For classical compressed sensing with plain \(\ell_1\) minimization, the bound is
\[
w^{2}\big(\mathcal{D}(\|x\|_{1},x^{\star})\big)
\le
2s\log\!\Big(\frac{n}{s}\Big) + \frac{7}{5}s,
\]
so
\[
m \ge 2s\log\!\Big(\frac{n}{s}\Big) + \frac{7}{5}s + 1
\]
suffices for an \(s\)-sparse signal [1408.5250].

The \(L_{1}\)-\(L_{1}\) theory introduces support and quality parameters that quantify how the prior aligns with the target:
\[
I := \{i: x_i^{\star}\neq 0\},
\qquad
J := \{j: x_j^{\star}\neq w_j\},
\qquad
s := |I|.
\]
It distinguishes sign-based good and bad components on \(I\):
\[
\text{good: } I_{+}J_{-}\cup I_{-}J_{+},
\qquad
\text{bad: } I_{+}J_{+}\cup I_{-}J_{-},
\]
with counts
\[
h := |I_{+}J_{-}| + |I_{-}J_{+}|,
\qquad
\overline{h} := |I_{+}J_{+}| + |I_{-}J_{-}|.
\]
It also uses
\[
q := |I\cup J|
\]
and the balance parameter
\[
\xi := |\{i: w_i\neq 0 = x_i^{\star}\}| - |\{i: w_i = x_i^{\star}\neq 0\}|.
\]
For \(\beta=1\), assuming \(\overline{h}>0\) and that there exists \(i\) with \(x_i^{\star}=w_i=0\), the simplified bound is
\[
m \ge 2\,\overline{h}\,\log\!\Big(\frac{n}{s+\tfrac{\xi}{2}}\Big)
+ \frac{7}{5}\Big(s+\frac{\xi}{2}\Big) + 1.
\]
The dominant logarithmic term scales with \(\overline{h}\) rather than \(s\). Since \(\overline{h}\le s\), good prior information can dramatically reduce the measurement requirement [1408.5250].

By contrast, the \(L_{1}\)-\(L_{2}\) bounds depend on magnitudes as well as supports and signs, through parameters such as \(K^{=}\), \(K^{\neq}\), \(v_{\beta}\), and \(\overline{w}\). Those bounds are explicitly described as looser and often track classical compressed sensing closely, which explains why \(L_{1}\)-\(L_{2}\) brings limited benefit relative to \(L_{1}\)-\(L_{1}\) unless the target is extremely sparse and the prior is very special [1408.5250].

## 4. Why PICCS helps in the TV or analysis setting

The geometric interpretation carries directly into analysis-sparse and TV-based PICCS. Adding the \(\|x-w\|_{1}\) or \(\|\Psi(x-w)\|_{1}\) term tightens the relevant sublevel set around directions that preserve similarity to the prior. When the prior has favorable sign alignment, the tangent or descent cone shrinks drastically, and only the bad components keep that cone wide [1408.5250].

For CT, this means that the benefit of PICCS is not determined solely by pixelwise similarity of the prior image. What matters is similarity in the sparsifying domain, typically the discrete gradient or another analysis operator. When \(\Psi x^{\star}\) and \(\Psi(x^{\star}-w)\) share support with favorable sign alignment, the number of projections required for accurate reconstruction can be substantially reduced relative to plain TV-regularized compressed sensing [1408.5250].

This point is central to the distinction between PICCS and quadratic prior methods. A TV or \(\ell_1\) penalty on the difference encodes sparsity of changes, whereas an \(\ell_2\) penalty around the prior acts more like averaging around a reference. The geometric analysis shows that the latter leaves the descent cone relatively wide, especially in high dimensions, so the measurement advantage is much smaller [1408.5250].

In repeated CBCT, the assumptions that make PICCS effective are explicit: a good prior image exists, the current scan is strongly under-sampled in angle, changes are local and sparse, and the prior is aligned to the current geometry. Under those conditions, PICCS suppresses undersampling artifacts and noise while preserving clinically relevant changes [2509.08574].

## 5. Parameterization, robustness, and algorithms

For \(L_{1}\)-\(L_{1}\) reconstruction, the theoretical bounds are minimized at \(\beta=1\), equivalently \(\lambda=1\) in the noiseless constrained form. This choice is described as universal, because it does not depend on unknown quantities in \(x^{\star}\) or \(w\), and phase transition curves for \(L_{1}\)-\(L_{1}\) typically achieve their best performance near \(\beta=1\) [1408.5250]. For \(L_{1}\)-\(L_{2}\), the bound-minimizing \(\beta\) depends on unknown magnitudes in \(x^{\star}\) and \(w\) and is therefore impractical in inverse problems; small \(\beta\) avoids the severe degradation seen for large \(\beta\) [1408.5250].

In noisy problems, the constrained form \(\|Ax-y\|_{2}\le \varepsilon\) yields a stability guarantee under the Gaussian-width condition stated above. In practice, \(\beta=1\) or \(\lambda=1\) remains the recommended default for \(L_{1}\)-\(L_{1}\), while the Lagrangian weight \(\tau\) can be tuned by standard discrepancy or cross-validation rules [1408.5250]. In the CBCT implementation, practical guidance is more application-specific: increase \(\alpha\) with higher noise or undersampling to suppress streaking and noise; increase \(\lambda\) when prior fidelity is high and changes are sparse and local; decrease \(\lambda\) when misregistration or large anatomical changes are suspected; and keep the smoothing parameter \(\tau\) small but nonzero [2509.08574].

All of the formulations are convex. For the synthesis problems, the basic proximal operators are explicit:
\[
\mathrm{prox}_{\tau\|\cdot\|_{1}}(z)=\operatorname{soft}(z;\tau),
\]
and
\[
\mathrm{prox}_{\tau\lambda\|\cdot-w\|_{1}}(z)
=
w+\operatorname{soft}(z-w;\tau\lambda).
\]
For the quadratic prior form,
\[
\mathrm{prox}_{\tau(\|\cdot\|_{1}+\tfrac{\beta}{2}\|\cdot-w\|_{2}^{2})}(z)
=
\operatorname{soft}\!\Big(\mu;\frac{\tau}{1+\tau\beta}\Big),
\qquad
\mu=\frac{z+\tau\beta w}{1+\tau\beta}.
\]
Suitable solvers include proximal gradient or FISTA for Lagrangian forms, ADMM for constrained and analysis or TV formulations, and primal-dual methods such as Chambolle–Pock for TV or analysis regularization [1408.5250].

The 2025 real-time CBCT work implements PICCS through an iteratively reweighted norm majorization-minimization scheme. Non-smooth TV penalties are replaced by quadratically weighted \(\ell_2\) terms with weights updated across outer iterations, and the inner least-squares problem is solved by LSQR or, equivalently, CGLS; that work uses CGLS because of its short-recurrence, low-memory, and GPU-friendly nature [2509.08574]. At outer iteration \(k\), the weights are computed from the current iterate,
\[
\mathbf{W}_k^{(1)} = \widetilde{\mathbf{W}(\mathbf{D}\boldsymbol{x}_{k-1})},
\qquad
\mathbf{W}_k^{(2)} = \widetilde{\mathbf{W}(\mathbf{D}(\boldsymbol{x}_{k-1}-\boldsymbol{x}_p))},
\]
and the quadratic subproblem is
\[
\boldsymbol{x}_k
=
\arg\min_{\boldsymbol{x}}
\;\;
\|\mathbf{A}\boldsymbol{x}-\boldsymbol{b}\|_2^2
+
\alpha^2 \|\mathbf{W}_k^{(1)} \mathbf{D}\boldsymbol{x}\|_2^2
+
\lambda^2 \|\mathbf{W}_k^{(2)} \mathbf{D}(\boldsymbol{x}-\boldsymbol{x}_p)\|_2^2.
\]
This is cast as an augmented least-squares problem and solved without explicitly building the matrix, using forward projection, backprojection, finite-difference stencils, and diagonal weighting operations on the GPU through the TIGRE toolbox [2509.08574].

The computational profile is dominated by one forward projection and one backprojection per inner CGLS iteration. IRN converges to the minimizer of the smoothed objective, whereas CGLS can exhibit semi-convergence if iterated too long, so early stopping is used [2509.08574].

## 6. Empirical behavior, variants, and limitations

The empirical results separate three issues: the effect of prior information itself, the difference between \(L_{1}\)-\(L_{1}\) and \(L_{1}\)-\(L_{2}\)-type priors, and the effect of algorithmic implementation. In a representative synthetic example with \(n=1000\) and \(s=70\), with prior \(w\) having about \(50\%\) relative \(\ell_2\) error and mixed support overlap, the \(L_{1}\)-\(L_{1}\) bound predicted successful recovery with \(m\approx 136\) measurements versus \(\approx 472\) for classical compressed sensing, and the empirical phase transition confirmed a drastic reduction. \(L_{1}\)-\(L_{2}\) tracked classical compressed sensing closely and provided negligible savings [1408.5250].

In the 2025 synthetic head phantom experiment with 20 projections and a prior equal to the head without tumor, FDK was qualitatively poor due to extreme undersampling, SIRT was smooth but blurry, and CGLS retained artifacts. For the prior-based methods, IRN-PICCS with 25 iterations achieved PSNR \(22.40\), SSIM \(0.726\), HaarPSI \(0.530\), and time \(19\) s; IRN-PICCS with 100 iterations and 4 outer cycles achieved PSNR \(24.18\), SSIM \(0.807\), HaarPSI \(0.571\), and time \(1{:}16\); ASD-POCS-PICCS with 20 iterations achieved PSNR \(23.88\), SSIM \(0.801\), HaarPSI \(0.657\), and time \(2{:}46\) [2509.08574]. In that synthetic setting, IRN-PIPLE produced the strongest quantitative results and preserved tumor and texture very well, whereas IRN-PICCS identified changes but was more aggressive on texture.

For the thorax phantom with real data, a metal needle, and 180, 50, and 20 projections, IRN-PICCS with 20 inner iterations yielded PSNR \(41.70\), SSIM \(0.918\), HaarPSI \(0.697\), and time \(1{:}46\) at 180 projections; PSNR \(43.76\), SSIM \(0.945\), HaarPSI \(0.777\), and time \(1{:}43\) at 50 projections; and PSNR \(43.73\), SSIM \(0.948\), HaarPSI \(0.816\), and time \(1{:}41\) at 20 projections [2509.08574]. At all undersampling levels, the proposed IRN-PIPLE and IRN-PICCS materially outperformed FDK, CGLS, and TV-only methods in PSNR and SSIM and reduced metal artifacts and noise much better. Runtime speedups over ASD-POCS-PICCS were substantial: about \(60\times\) at 180 projections, about \(4.4\times\) at 50 projections, and about \(1.6\times\) at 20 projections [2509.08574].

These results clarify a common confusion. The geometric theory strongly favors \(L_{1}\)-\(L_{1}\) over \(L_{1}\)-\(L_{2}\) as a general compressed sensing design principle, because only the former replaces sparsity \(s\) by the smaller effective bad-component count \(\overline{h}\) in the dominant measurement term [1408.5250]. However, the CBCT study found that a PIPLE-style \(\ell_2\) prior could outperform PICCS in some heavily undersampled practical settings, particularly in clean synthetic data, where it had a more averaging effect and could be more robust when the prior was strong but changes were small [2509.08574]. This suggests that theoretical recovery geometry and application-level image-quality tradeoffs are related but not identical.

The strengths and limitations are correspondingly regime-dependent. PICCS excels when changes are small and localized, the prior is trustworthy, and undersampling artifacts must be suppressed without erasing true changes [2509.08574]. Its benefit diminishes when the prior is poor, when there are many bad components in the transform domain, or when misregistration introduces structured discrepancies [1408.5250]. Large global changes or overly large \(\lambda\) can bias the reconstruction toward the prior, and very strong metal artifacts remain challenging [2509.08574]. Best practices therefore include registration of the prior to the current geometry, moderate initialization of regularization weights, inspection of difference images \(x-x_p\), and comparison with TV-only baselines to verify that the prior term is not suppressing genuine changes [2509.08574].

Source: https://www.emergentmind.com/topics/prior-image-constrained-compressed-sensing-piccs