---
title: Matrix-Weighted X-Ray Transform
url: https://www.emergentmind.com/topics/matrix-weighted-x-ray-transform
type: topic
---

# Matrix-Weighted X-Ray Transform

The matrix-weighted X-ray transform generalizes the classical scalar geodesic ray transform by introducing a matrix-valued weight, thereby enabling the integration of vector- or matrix-valued functions along geodesics with multiplicative matrix attenuation or more general coupling. This framework underpins a range of inverse problems involving connections, Higgs fields, and non-Abelian scattering data, with significant implications for geometric tomography, quantum and polarization tomography, and microlocal analysis.

## 1. Definition and Basic Framework

Let $(M,g)$ be a compact, oriented Riemannian manifold of dimension $n\geq 2$ with smooth, strictly convex boundary. The unit tangent bundle is denoted $SM = \{(x,v)\in TM : |v|_g=1\}$. For $(x,v)\in SM$, let $\gamma_{x,v}(t)$ be the maximal unit-speed geodesic starting at $x$ in the direction $v$, with $\tau_+(x,v)$ its first exit time from $M$.

Given a continuous, typically smooth (or real-analytic in recent developments) map $W: SM \rightarrow \mathrm{GL}(N, \mathbb{C})$, called a *matrix weight*, and a function $f : M \to \mathbb{C}^N$, the matrix-weighted geodesic X-ray transform $I_Wf$ is defined as:
$$
I_W f(x,v) = \int_0^{\tau_+(x,v)} W(\gamma_{x,v}(t), \dot\gamma_{x,v}(t)) f(\gamma_{x,v}(t)) dt,\quad (x,v)\in\partial_+ SM
$$
where $\partial_+ SM = \{(x,v)\in SM : x\in \partial M, \langle v, \nu(x) \rangle_g \geq 0\}$ and $\nu(x)$ is the inward unit normal. In the non-Abelian (scattering) variant, one is given a matrix field $\Phi: M \to \mathfrak{gl}(N,\mathbb{C})$ and considers the fundamental matrix solution $U(t)$ of the ODE:
$$
\frac{d}{dt}U(t) + \Phi(\gamma(t)) U(t) = 0,\quad U(0) = I,
$$
with $C_\Phi(x,v) := U(\tau_+(x,v))$ forming the inverse problem’s data $I(\Phi) = C_\Phi$ [1905.00860][1605.07894][2506.24067].

## 2. Relation to Classical and Attenuated X-Ray Transforms

The matrix-weighted transform strictly generalizes the scalar attenuated X-ray transform. When $W$ is scalar, or $\Phi(x) = a(x) \cdot I$, the ODE becomes separable, yielding:
$$
C_\Phi(x,v) = \exp\Big(- \int_0^{\tau_+(x,v)} a(\gamma(t)) dt \Big) I.
$$
Thus, $\log C_\Phi(x,v)$ recovers the classical scalar attenuated X-ray transform of $a(x)$. In the general case, $W$ accounts for the effects of parallel transport in a vector bundle with connection, Higgs field, or more general matrix-valued attenuation or weight [1905.00860][2506.24067][1605.07894].

## 3. Injectivity and Uniqueness Results

### 3.1 Two-Dimensional and Real-Analytic Settings
On simple surfaces ($n=2$, no conjugate points, strictly convex boundary), the transform $I_W$ is injective for smooth matrix weights $W$ and for the scattering case $I: \Phi \mapsto C_\Phi$ with $\Phi$ valued in $\mathfrak{so}(n)$ or $\mathfrak{u}(n)$ [1905.00860][2006.02257]. For real-analytic metrics and weights, local injectivity near a strictly convex boundary point, as well as global injectivity on analytic surfaces, is established for any real-analytic $\mathrm{GL}(N,\mathbb{C})$-valued weight [2506.24067]. In the non-Abelian (connection/Higgs field) setting, the scattering data locally determines the Higgs field near strictly convex boundary, and global uniqueness holds by analytic continuation [2506.24067].

### 3.2 Higher Dimensions and Foliation
For $n\geq 3$, global injectivity holds for smooth matrix weights under the “foliation condition”: existence of a smooth strictly convex exhaustion function, whose sublevel sets form strictly convex hypersurfaces throughout $M$ [1605.07894][2008.07951][1901.03525]. This enables layer-stripping arguments that successively recover $f$ across layers from the boundary inward.

### 3.3 Piecewise Constant Functions and Minimal Regularity
For piecewise constant $f$ and continuous injective matrix weights $W$, injectivity holds in all dimensions $n\geq 2$ under nontrapping and strictly convex boundary (plus foliation for $n\geq 3$), without derivative bounds or assumption regarding conjugate points [1901.03525].

### 3.4 Generic and Twisted Matrix Weights
Generic injectivity for the X-ray transform twisted by a unitary connection on a Hermitian vector bundle (i.e., parallel transport matrix weights) has been proven in dimension $n\geq 3$ on Anosov or strictly convex no-conjugate-point manifolds [2107.05119]. The microlocal-perturbative method used extends to matrix weights of finite Fourier degree in velocity.

## 4. Stability, Reconstruction, and Statistical Inversion

Quantitative stability results are established both in two and higher dimensions. In the two-dimensional non-Abelian case, there exist constants depending on $M$ and $C^1$-norms of $\Phi, \Psi$ such that:
$$
\|\Phi - \Psi\|_{L^2(M)} \leq c(\Phi, \Psi)\| C_\Phi C_\Psi^{-1} - I \|_{H^1(\partial_+ SM)}
$$
with an explicit form for $c(\Phi,\Psi)$. This enables statistical inversion: using a Bayesian framework with Gaussian process priors for $\Phi$, noisy measurements of $C_\Phi$ are inverted by (infinite-dimensional) pCN-MCMC, with contraction rates for both the posterior and the mean estimator. For smooth $\Phi$, the $L^2$-error in recovery decays with the optimal $N^{-1/2}$ rate in the number $N$ of measurements [1905.00860]. 

In higher dimensions, H\"older-type stability estimates relate the $L^2$-norm of the error in $\Phi$ on compact sets to a H\"older power of the $L^2$-difference of the data $C_{\Phi_1} - C_{\Phi_2}$ on suitable geodesic classes, with explicit control of exponents and constants in terms of $C^k$-norms of the weights [2008.07951].

## 5. Microlocal and Analytic Techniques

Recent advances rely on analytic microlocal analysis and the double fibration framework: the matrix-weighted transform is realized as a Fourier integral operator associated to the twisted conormal bundle of the incidence variety in $G\times X$. Under the Bolker condition (no conjugate points), analytic or smooth singularities of $f$ are propagated to those of $I_W f$. In the real-analytic setting, propagation of analytic wavefront sets via the matrix-weighted transform allows for unique analytic continuation, and hence global injectivity, even from partial data near a strictly convex boundary point [2506.24067].

## 6. Applications, Obstructions, and Extensions

Typical applications include quantum state tomography, where a matrix-valued Hamiltonian is probed via its propagator, and polarization tomography, involving tensor weights and recovery of perturbations. In the setting involving both connection and Higgs field, the kernel of the matrix-weighted X-ray transform consists of gauge transformations, i.e., $h \mapsto h - (\nabla + A + \Phi)p$ for suitable $p$ vanishing at the boundary. Uniqueness is then modulo this natural gauge [1605.07894].

The methodology accommodates Banach-space valued $f$ and operator-valued, merely continuous weights, provided injectivity in the fiber [1901.03525]. Inverse problems with minimal regularity, e.g., piecewise-constant targets and continuous weights, are covered. However, for smooth weights and functions, non-injectivity can occur on simple 2-manifolds (cf. Boman counterexample), indicating the sharpness of regularity conditions.

## 7. Overview Table: Notions and Key Results

| Setting                    | Hypotheses            | Result Type                  |
|----------------------------|-----------------------|------------------------------|
| 2D, real-analytic          | Strictly convex, analytic boundary | Local/global injectivity for analytic weights [2506.24067] |
| 2D, smooth/simple surface  | Strictly convex, no conjugate points | Injectivity for $\mathfrak{u}(n)$ / $\mathfrak{so}(n)$ weights [1905.00860][2006.02257] |
| $n \geq 3$, smooth         | Foliation, strictly convex boundary | Global injectivity for smooth weights [1605.07894][1901.03525] |
| Any $n$, piecewise constant $f$ | Strictly convex boundary, foliation ($n\ge3$) | Injectivity for continuous weights [1901.03525] |
| Non-Abelian scattering     | Above + gauge constraint              | Unique recovery up to gauge    |
| Generic matrix weights     | Anosov or strictly convex, no conjugate points | Generic injectivity [2107.05119] |

The matrix-weighted X-ray transform thus provides the analytic core for a broad spectrum of inverse boundary and tomography problems, with the interplay of geometric, analytic, microlocal, and statistical techniques yielding the current landscape of injectivity, stability, and practical inversion methodologies.

Source: https://www.emergentmind.com/topics/matrix-weighted-x-ray-transform