---
title: Signed Cumulative Distribution Transform
url: https://www.emergentmind.com/topics/signed-cumulative-distribution-transform-scdt
type: topic
---

# Signed Cumulative Distribution Transform

The Signed Cumulative Distribution Transform (SCDT) is an invertible, transport-based representation for one-dimensional signed signals that extends the cumulative distribution transform (CDT) from nonnegative unit-mass densities to arbitrary finite signed signals. Its defining construction is to decompose a signal into its positive and negative parts, apply CDT to the normalized components relative to a fixed positive reference, and retain the component masses. In this way, the SCDT separates sign, mass, and geometry, while preserving the quantile-based transport structure that makes translations, dilations, and more general monotone deformations tractable in transform space [2106.02146][2207.07989][2606.11432].

## 1. Definition and formal construction

For a nonnegative unit-mass signal or density \(s\), with cumulative distribution function \(F_s\), and a fixed strictly positive reference \(s_0\) with cumulative distribution function \(F_{s_0}\), the CDT is the monotone transport map
\[
\mathcal{C}(s)=F_s^{\dagger}\circ F_{s_0},
\]
where \(F^\dagger\) denotes the generalized inverse, or quantile function. In the notation used in several SCDT papers, the transformed signal is the transport map itself; for a uniform reference on \([0,1]\), it reduces to the quantile function \(F_s^{-1}\) [2207.07989][2110.05606].

For a signed signal \(s\in L^1(\mathbb{R})\), the SCDT begins with the Jordan decomposition
\[
s(x)=s^+(x)-s^-(x),\qquad s^+(x)=\max\{0,s(x)\},\qquad s^-(x)=\max\{0,-s(x)\}.
\]
Let
\[
m_+=\|s^+\|_1,\qquad m_-=\|s^-\|_1.
\]
When \(m_\pm>0\), define normalized nonnegative components \(p_\pm=s^\pm/m_\pm\). The SCDT then stores the CDT of each normalized component together with its mass:
\[
\widehat{s}=\big((s^+)^\star,\|s^+\|_1,(s^-)^\star,\|s^-\|_1\big),
\]
or equivalently,
\[
\mathcal{S}(s)=\bigl(m_+,\widehat p_+,m_-,\widehat p_-\bigr),
\]
depending on the notation adopted in a given paper [2207.07989][2606.11432].

A measure-theoretic formulation places the construction on finite signed measures on the extended real line. If \(\nu=\nu^+-\nu^-\) and \(\mu_0\) is a non-trivial, atomless, finite positive reference measure, then the SCDT is
\[
\mathcal{T}_{\mu_0}(\nu)=\big((\nu^+)^\star,\|\nu^+\|,(\nu^-)^\star,\|\nu^-\|\big),
\]
with each \((\nu^\pm)^\star\) obtained from the generalized inverse of the cumulation of \(\nu^\pm\) composed with the reference cumulation. This formulation makes the transform applicable beyond smooth densities and clarifies its status as a bijection between signed measures and pairs of monotone transport maps plus masses [2106.02146].

A recurring point in the literature is that the SCDT is not a single transport map for a signed signal. Rather, it is a paired representation: one transport map for the positive part, one for the negative part, and mass coordinates for both. This separation is the mechanism by which the transform extends optimal-transport-style geometry to signed data [2307.15339][2207.07989].

## 2. Inversion, metric structure, and optimal-transport interpretation

The SCDT is invertible. In the density setting, if \(T_+=(s^+)^\star\) and \(T_-=(s^-)^\star\), then the original signed signal is reconstructed from the inverse transport maps and the reference density \(s_0\) by
\[
s^+(x)=m_+\, s_0(T_+^{-1}(x))\,\frac{d}{dx}T_+^{-1}(x),\qquad
s^-(x)=m_-\, s_0(T_-^{-1}(x))\,\frac{d}{dx}T_-^{-1}(x),
\]
followed by
\[
s(x)=s^+(x)-s^-(x).
\]
Equivalent formulas appear in several variants, including the form
\[
s(x)=m_+\,\mathrm{CDT}^{-1}[\widehat p_+](x)-m_-\,\mathrm{CDT}^{-1}[\widehat p_-](x),
\]
which emphasizes inversion of the component CDTs before rescaling and recombination [2207.07989][2606.11432].

In the measure-theoretic setting, inversion is expressed as a pushforward:
\[
(f,r,g,s)\mapsto r\,f_{\#}\!\left(\frac{\mu_0}{\|\mu_0\|}\right)-s\,g_{\#}\!\left(\frac{\mu_0}{\|\mu_0\|}\right).
\]
This formulation is important because it does not rely on classical derivatives and therefore extends to finite signed measures and generalized inverse maps [2106.02146].

The geometric content of the SCDT is inherited from one-dimensional monotone optimal transport. For positive densities, the CDT is the unique monotone rearrangement mapping the reference to the target. For signed signals, the same transport geometry is applied channel-wise to the positive and negative parts. This leads to a Wasserstein-type metric
\[
\begin{aligned}
D_S^2(r,s)&=
d_{W^2}^2\!\left(\frac{r^+}{\|r^+\|_1},\frac{s^+}{\|s^+\|_1}\right)
+\left|\|r^+\|_1-\|s^+\|_1\right|^2 \\
&\quad+
d_{W^2}^2\!\left(\frac{r^-}{\|r^-\|_1},\frac{s^-}{\|s^-\|_1}\right)
+\left|\|r^-\|_1-\|s^-\|_1\right|^2,
\end{aligned}
\]
which is exactly the squared norm difference of the SCDT representations:
\[
D_S^2(s,r)=\|\widehat{s}-\widehat{r}\|_{(L^2(s_0(x)\,dx)\times\mathbb{R})^2}^2.
\]
Thus, transport distances between signed signals reduce to Euclidean distances between transform coordinates, with separate contributions from transport geometry and mass mismatch [2207.07989][2106.02146].

This Hilbert-space embedding explains why subspace methods, least-squares estimators, and orthogonal projections recur throughout the SCDT literature. The transform replaces a nonlinear comparison problem in signal space with an \(L^2\)-type problem on monotone maps and mass coordinates [2307.15339][2510.00148].

## 3. Transformation laws and linearization properties

The central structural property of the SCDT is its behavior under mass-preserving monotone warps. If \(g:\mathbb{R}\to\mathbb{R}\) is a strictly increasing differentiable bijection and
\[
s_g(t)=g'(t)\,s(g(t)),
\]
then the SCDT satisfies
\[
\widehat{s}_g=
\left(g^{-1}\circ (s^+)^\star,\|s^+\|_1,\;g^{-1}\circ (s^-)^\star,\|s^-\|_1\right).
\]
The masses are preserved, while the transport maps are composed with \(g^{-1}\). This is the precise sense in which domain deformations become simple operations in transform space [2207.07989][2106.02146].

Several important special cases follow immediately. For an affine warp \(g(t)=\omega t-\tau\) with \(\omega>0\),
\[
\widehat{s}_g=
\left(\frac{(s^+)^\star+\tau}{\omega},\|s^+\|_1,\frac{(s^-)^\star+\tau}{\omega},\|s^-\|_1\right).
\]
Hence translations act additively on each transport component, and dilations act multiplicatively. Amplitude scaling behaves differently: multiplying the signal by a positive scalar leaves the CDT shape coordinates unchanged and scales only the masses [2207.07989][2308.12259].

For the positive-density CDT, this linearization is exact and especially simple for rigid shifts. If \(w_s(x)=w(x-s)\), then
\[
\widehat{w_s}(\alpha)=\widehat{w}(\alpha)+s,
\]
so a translation family becomes an affine line parallel to the constant mode in \(L_r^2\). The recent analysis of additive perturbations in CDT coordinates shows that, under a local nondegeneracy condition, additive noise in physical space induces a nonlocal perturbation in transform space through the primitive of the noise weighted by the reciprocal density. In particular, perturbations are amplified in low-density regions. For signed signals, that paper considers an SCDT analogue and uses numerical feature matching for shift estimation rather than a closed-form constant-mode projection [2606.11432].

A further consequence is convexification of deformation classes. If a signal class is generated from a template by increasing warps \(g_j\), then in transform space the class has the form \(g_j^{-1}\circ\widehat{\phi}\). The literature states that the transformed class is convex for every template if and only if the inverse deformation family is convex. In finite-data settings, this often produces “thin” or approximately low-dimensional sets, which motivates nearest-subspace and local-subspace models for classification and regression [2110.05606][2205.00348].

The transform therefore does not merely encode signals differently; it reorganizes deformation variability into additive, multiplicative, or compositional structure. This suggests why linear estimation procedures in SCDT space can succeed even when the native-domain problem is nonlinear or nonconvex.

## 4. Discrete, numerical, and algorithmic formulations

Although the classical CDT is stated for continuous densities, recent work has developed a fully discrete CDT and discrete SCDT for atomic measures on the real line. With a fixed atomic reference
\[
\sigma=\sum_{j=1}^m q_j\delta_{y_j}
\]
and a target probability measure
\[
\mu=\sum_{i=1}^n p_i\delta_{x_i},
\]
the discrete CDT is defined by
\[
T_\mu=F_\mu^{-1}\circ F_\sigma,\qquad T_\mu(y_j)=F_\mu^{-1}(Q_j),
\]
where \(Q_j\) are the reference cumulative masses. The inverse reconstruction is the pushforward
\[
(T_\mu)_\#\sigma=\sum_{j=1}^m q_j\delta_{T_\mu(y_j)}.
\]
In the signed case, if \(f=\sum_i a_i\delta_{x_i}\), with positive and negative masses \(m^\pm\) and normalized channels \(\mu^\pm\), then the discrete SCDT is
\[
S_\sigma(f)=(m^+,T^+,m^-,T^-),\qquad T^\pm=F_{\mu^\pm}^{-1}\circ F_\sigma,
\]
and reconstruction is
\[
\widetilde f=m^+(T^+)_\#\sigma-m^-(T^-)_\#\sigma.
\]
This framework exposes a genuine finite-resolution obstruction: deterministic atomic transport cannot split mass, so exact reconstruction holds only under a cumulative-mass compatibility criterion, while weak convergence is recovered under reference refinement [2606.12131].

The discrete theory also introduces thresholded stabilization near zero crossings. Because a small perturbation can flip a discrete coefficient from the positive channel to the negative channel, a dead-zone rule around zero can suppress unstable channel switching. This is a specifically signed phenomenon and highlights that zero crossings are a numerical singular locus for the split-channel representation [2606.12131].

For sampled continuous signals, practical SCDT computation is based on positive/negative splitting, normalization, cumulative sums, generalized inverse evaluation, and interpolation. The numerical workload is typically linear in the number of samples for CDF construction and quantile inversion. In the 2026 shift-recovery study, CDT/SCDT computation on a grid is described as \(O(n)\) per signal for CDF/quantile construction, with known-template CDT shift estimation also \(O(n)\), unknown-template de-shift-and-average \(O(Nn)\), and SCDT numerical matching \(O(|\mathcal{G}|\,n)\) per observation per iteration [2606.11432].

These computational descriptions show a persistent pattern across the literature: the main numerical ingredients are monotone cumulative functions and their inverses. No transport plan matrix is needed in the one-dimensional setting, which is one reason the SCDT remains attractive for large-scale or repeatedly solved problems [2606.12131][2510.00148].

## 5. Applications in signal analysis, estimation, and learning

One major application is alignment and shift recovery. For positive densities, CDT coordinates permit exact linear recovery of translations through projection onto the constant mode when the template is known, and de-shift-and-average procedures when multiple shifted observations are available. In the signed case, the 2026 study uses the SCDT as a transport feature map, estimates shifts numerically by feature matching over a grid, and recovers unknown templates by alternating alignment and averaging. Its numerical experiments report effective recovery for both density-valued and signed signals, with stronger difficulty for discontinuous square-wave signals under stronger noise [2606.11432].

A second major application is classification. The nearest-subspace search method in SCDT space represents each class by the span of transformed training samples and classifies by projection residual. On synthetic data generated by polynomial warps of prototype signals, the reported result is that the SCDT-NS method achieves near-perfect accuracy with only 16 training samples per class, while none of the compared CNNs attained perfect accuracy even with 256 training samples per class. On an ECG heartbeat classification problem with three classes, the reported results are: DeepConvNet \(47.57\%\), F1 \(0.4065\); ShallowConvNet \(33.68\%\), F1 \(0.2618\); CompactConvNet \(29.59\%\), F1 \(0.2466\); SCDT-NS \(61.50\%\), F1 \(0.5979\) [2110.05606].

An end-to-end extension replaces global class subspaces by nearest local subspaces enriched with analytic deformation directions. On ten 1D time-series datasets, the reported method achieves the top accuracy in \(4/10\) datasets, matches the average arithmetic ranking of \(1\)NN-DTW at \(2.2\), attains mean per class error \(0.038\), and is described as more data efficient than several deep-learning baselines. The same study also reports that on a synthetic out-of-distribution setup the proposed method achieved near-perfect accuracy with approximately 16 training samples per class [2205.00348].

A third application is parameter estimation. For composition models \(s_g(x)=g'(x)s(g(x))\), the SCDT converts Wasserstein-type distance minimization into linear least squares in transform space. For polynomial \(g_p\), the composition \(g_p\circ \widehat r\) becomes linear in the unknown coefficients, so the estimation problem admits a global minimizer via linear least squares. The reported examples include time delay and dispersion estimation, where native-domain objectives are described as nonconvex and multimodal while the SCDT-domain objectives become convex with a unique global minimum [2207.07989].

This least-squares perspective is developed further for dynamical systems and PDE parameter identification. In that setting, one measures a time signal at a fixed sensor, models it as a warped template induced by the governing PDE, and uses SCDT nearest-local-subspace regression for coarse parameter recovery. Reported results include nonlinearity detection accuracy of \(98.0\%\) for SCDT-NLS versus \(95.5\%\) for FT-SVM and \(91.9\%\) for 1D-VGG, dispersion detection accuracy of \(99.0\%\), and near-perfect damage-level identification on the UNESP-CONCEPT structural health monitoring dataset, including \(100.0\%\) for one sensor and \(99.67\%\) for another [2308.12259].

The framework has also been lifted to images through the Radon Signed Cumulative Distribution Transform (RSCDT), which applies the one-dimensional SCDT to Radon projections. Reported results include \(100\%\) accuracy on a simulated signed-image task for RSCDT-NS versus \(49\%\) for an unsigned RCDT-NS baseline that used absolute values, as well as strong performance on geometric shape and sign-language datasets [2307.15339].

More recent work applies SCDT to hyperspectral anomaly detection. In that setting, each spectrum is modeled as a signed or preprocessed one-dimensional signal, mapped to SCDT coordinates, and background variability is learned by PCA in transform space. The reported AUC values are \(0.8477\) on AVIRIS-I, \(0.6452\) on AVIRIS-II, \(0.8091\) on Urban, \(0.8525\) on Pavia, and \(0.6447\) on Forest for the FPR range \(\le 10^{-2}\), with superiority on four datasets at low false-positive rates and the highest full-range AUC on all five [2510.00148].

## 6. Assumptions, limitations, and current scope

The SCDT is inherently a one-dimensional construction. Its exact optimal-transport interpretation depends on monotone rearrangement on the line, and most theoretical statements assume a strictly positive reference density or an atomless reference measure. For practical density formulas, the positive and negative components are usually assumed absolutely continuous on their supports, while generalized inverses are used to handle non-strict monotonicity [2106.02146][2207.07989].

The transform is best matched to variability generated by monotone, mass-preserving deformations. Non-monotone warps, severe multipath interference, or signal classes not well described by transport of a template can degrade convexity and subspace structure. One study explicitly notes that gearbox vibration signals in the raw time domain do not satisfy the generative assumptions well, whereas performance improves after transformation to the Fourier domain, where finite-support event structure is more compatible with the model [2205.00348].

Noise sensitivity is structured rather than benign. In the CDT perturbation analysis, additive noise in physical space becomes nonlocal in transform space through integration, and perturbations are amplified in low-density regions by the factor \(1/u(\widehat u(\alpha))\). Even white noise becomes correlated after transformation. In density recovery under strong noise, positivity clipping and renormalization may be needed to preserve admissibility, which departs from a purely additive model [2606.11432].

The signed extension introduces its own limitations. Positive and negative parts are transported independently, so interactions between signs are not represented by a single coupled transport map. Near zero crossings, small perturbations can trigger channel switching, motivating thresholded stabilization in the discrete setting [2606.12131]. In alignment problems, the positive-density CDT admits a closed-form constant-mode projection for translation estimation, whereas the SCDT workflow described in the 2026 shift-recovery study uses numerical feature matching rather than a closed-form projection [2606.11432].

These constraints delimit the present scope of the SCDT rather than diminish its utility. Within its natural regime—one-dimensional signed signals whose dominant variability is transport-like—the transform provides a rare combination of bijectivity, explicit inversion, metric structure, deformation linearization, and computational simplicity. That combination explains its continuing use across classification, parameter estimation, inverse problems, signed-image analysis, and anomaly detection [2106.02146][2510.00148].

Source: https://www.emergentmind.com/topics/signed-cumulative-distribution-transform-scdt