Exact reconstruction formulas for a Radon transform over cones

Sintesi Tecnica: Formule di Ricostruzione Esatta per una Trasformata di Radon su Coni

Enunciato del Problema
Questo articolo affronta l'inversione matematica di una specifica variante della trasformata di Radon nota come trasformata di Radon conica, indicata come $R^{(p)}$. La trasformata integra una funzione regolare $f$ supportata nel semispazio superiore $\mathbb{R}^{d-1} \times (0, \infty)$ su superfici coniche monodirezionali. Questi coni sono caratterizzati da vertici situati sull'iperpiano $\mathbb{R}^{d-1} \times \{0\}$, un asse di simmetria ortogonale a tale piano (allineato con l'asse $y$) e un angolo di apertura semilunare $\theta \in (0, \pi/2)$. L'integrazione include un peso radiale $1/s^p$, dove $s$ rappresenta la distanza lungo il generatore del cono.

La motivazione primaria per lo studio di questa trasformata risiede nella tomografia computerizzata (CT) che tiene conto dello scattering Compton. In tali modalità di imaging, la rilevazione dei fotoni diffusi permette di tracciare la posizione della sorgente all'indietro verso la superficie di un cono, piuttosto che lungo una linea o un piano. Sebbene esistano formule di ricostruzione tramite trasformata di Fourier per casi specifici in tre dimensioni ($p \in \{0, 2\}$) e formule di retroproiezione per due dimensioni, il caso generale per un'arbitraria dimensione spaziale $d \geq 2$ e un parametro di peso arbitrario $p \in \mathbb{R}$ non era stato precedentemente studiato.

Metodologia
L'autore, Markus Haltmeier, impiega tecniche di analisi armonica, sfruttando specificamente la trasformata di Fourier e la trasformata di Hankel per derivare le formule di inversione. Il nucleo della metodologia consiste nell'stabilire un "identità della fetta di Fourier" per la trasformata di Radon conica, analoga alla ben nota formula della fetta di Fourier per la classica trasformata di Radon.

La derivazione procede in tre passaggi logici:

1. Riduzione al Caso Non Pesato: Viene stabilito un lemma che mostra come la trasformata $R^{(p)}$ per un peso arbitrario $p$ possa essere correlata alla trasformata non pesata $R^{(0)}$ tramite operatori moltiplicativi semplici che coinvolgono $\cos(\theta)$ e $y^{-p}$. Ciò consente agli autori di concentrarsi sulla derivazione del caso $p=0$ e di generalizzare i risultati successivamente.
2. Derivazione dell'Identità della Fetta di Fourier: Applicando la trasformata di Fourier rispetto alla variabile spaziale $u \in \mathbb{R}^{d-1}$ e utilizzando la rappresentazione integrale delle funzioni di Bessel sulla sfera $S^{d-2}$, l'autore deriva un'identità che lega la trasformata di Fourier-Hankel della funzione originale $f$ alla trasformata di Fourier dei dati della trasformata di Radon conica.
3. Inversione tramite Retroproiezione: Utilizzando la derivata identità della fetta di Fourier, l'autore inverte la trasformata applicando le trasformate di Fourier e di Hankel inverse. Questo processo conduce a esplicite formule di ricostruzione del tipo retroproiezione filtrata. Il passaggio di "filtraggio" viene identificato con l'applicazione del potenziale di Riesz $I^{(1-d)}$ rispetto al primo argomento.

Contributi Chiave e Risultati
L'articolo fornisce due formule di ricostruzione esplicite per il recupero di $f(x, y)$ dai dati $R^{(p)}$:

1. Parametrizzazione Angolare (Teorema 1.1, Eq. 1.3):
   $$f(x, y) = \frac{y^p}{(2\pi)^{d-1}} \int_0^{\pi/2} \frac{1}{\cos(\theta)^{1+p}} \int_{S^{d-2}} \left( I^{(1-d)} R^{(p)} f \right)(x + y \tan(\theta) n, \theta) \, dn \, d\theta$$
   Questa formula integra su tutti i coni che passano per il punto di ricostruzione $(x, y)$, parametrizzati dall'angolo di apertura $\theta$ e dalla direzione $n$ sulla sfera unitaria.

2. Parametrizzazione Spaziale (Teorema 1.1, Eq. 1.4):
   $$f(x, y) = \frac{1}{(2\pi)^{d-1}} \int_{\mathbb{R}^{d-1}} \frac{(|u-x|^2 + y^2)^{\frac{p-1}{2}}}{|u-x|^{d-2}} \left( I^{(1-d)} R^{(p)} f \right)\left(u, \arctan\frac{|u-x|}{y}\right) \, du$$
   Questa forma integra sul piano di base $\mathbb{R}^{d-1}$, sommando i contributi dei coni i cui vertici giacciono in $u$.

Dimensioni Specifiche:

• Per $d=2$: Il potenziale di Riesz $I^{-1}$ corrisponde alla composizione della derivata e della trasformata di Hilbert ($\partial_u H_u$).
• Per $d=3$: Il potenziale di Riesz $I^{-2}$ corrisponde al Laplaciano negativo ($-\Delta_u$). Di conseguenza, la formula di ricostruzione 3D comporta l'applicazione del Laplaciano ai dati seguita da una retroproiezione conica.

Significato e Rivendicazioni
L'articolo sostiene che queste siano le prime formule di ricostruzione esplicite del tipo retroproiezione filtrata per la trasformata di Radon conica in arbitrarie dimensioni $d \geq 2$ e per qualsiasi peso reale $p$.

• Località: L'autore nota una distinzione basata sulla parità della dimensione $d$. Per le dimensioni dispari, la ricostruzione è locale, il che significa che il valore di $f$ in un punto $(x, y)$ dipende solo dagli integrali su coni che passano per un intorno arbitrariamente piccolo di quel punto. Per le dimensioni pari, la ricostruzione è non locale, richiedendo la conoscenza degli integrali su tutte le superfici coniche. Questo comportamento rispecchia quello della classica trasformata di Radon.
• Novità: Sebbene formule equivalenti per $d=2$ siano state precedentemente derivate in altri lavori, le formule per $d=3$ e dimensioni superiori sono presentate come nuovi contributi.
• Praticità: L'articolo afferma che le formule di ricostruzione 3D sono particolarmente semplici, consistendo nell'applicazione del Laplaciano seguita dalla retroproiezione conica. L'autore suggerisce che l'implementazione numerica dovrebbe essere agevole, seguendo i metodi stabiliti per le classiche trasformate di Radon e sferiche, sebbene gli studi numerici siano lasciati a ricerche future.

Il lavoro è inquadrato come una base teorica per la tomografia a emissione mediante camere Compton, fornendo gli strumenti matematici necessari per invertire i dati generati da tali dispositivi.