Exact reconstruction formulas for a Radon transform over cones

Technische Samenvatting: Exacte Reconstructieformules voor een Radon-transformatie over Kegels

Probleemstelling
Dit artikel behandelt de wiskundige inversie van een specifieke variant van de Radon-transformatie, bekend als de kegelvormige (conische) Radon-transformatie, aangeduid als $R^{(p)}$. De transformatie integreert een gladde functie $f$, ondersteund in de bovenste halve ruimte $\mathbb{R}^{d-1} \times (0, \infty)$, over eenzijdige kegelvormige oppervlakken. Deze kegels worden gekenmerkt door hun hoekpunten gelegen op het hypervlak $\mathbb{R}^{d-1} \times \{0\}$, een symmetrieas loodrecht op dit vlak (uitgelijnd met de $y$-as), en een halve openingshoek $\theta \in (0, \pi/2)$. De integratie omvat een radiale gewicht $1/s^p$, waarbij $s$ de afstand langs de generator van de kegel vertegenwoordigt.

De primaire motivatie voor het bestuderen van deze transformatie ligt in de computertomografie (CT) die rekening houdt met Comptonverstrooiing. In dergelijke beeldvormingsmodaliteiten maakt de detectie van verstrooide fotonen het mogelijk om de bronlocatie terug te traceren naar het oppervlak van een kegel, in plaats van naar een lijn of vlak. Hoewel er Fourier-reconstructieformules bestaan voor specifieke gevallen in drie dimensies ($p \in \{0, 2\}$) en back-projectieformules beschikbaar zijn voor twee dimensies, was het algemene geval voor een willekeurige ruimtelijke dimensie $d \geq 2$ en een willekeurige gewichtsparameter $p \in \mathbb{R}$ voorheen niet bestudeerd.

Methodologie
De auteur, Markus Haltmeier, maakt gebruik van harmonische analyse-technieken, specifiek door gebruik te maken van de Fourier-transformatie en de Hankel-transformatie, om de inversieformules af te leiden. De kern van de methodologie is het vaststellen van een "Fourier-snede-identiteit" voor de kegelvormige Radon-transformatie, analoog aan de bekende Fourier-snede-stelling voor de klassieke Radon-transformatie.

De afleiding verloopt in drie logische stappen:

1. Reductie tot het ongegewogen geval: Er wordt een lemma vastgesteld dat aantoont dat de transformatie $R^{(p)}$ voor een willekeurig gewicht $p$ kan worden gerelateerd aan de ongegewogen transformatie $R^{(0)}$ via eenvoudige multiplicatieve operatoren die $\cos(\theta)$ en $y^{-p}$ bevatten. Dit stelt de auteurs in staat om zich te concentreren op de afleiding voor het geval $p=0$ en de resultaten later te generaliseren.
2. Afleiding van de Fourier-snede-identiteit: Door de Fourier-transformatie toe te passen met betrekking tot de ruimtelijke variabele $u \in \mathbb{R}^{d-1}$ en gebruik te maken van de integraalrepresentatie van Bessel-functies over de sfeer $S^{d-2}$, leidt de auteur een identiteit af die de Fourier-Hankel-transformatie van de oorspronkelijke functie $f$ verbindt met de Fourier-transformatie van de gegevens van de kegelvormige Radon-transformatie.
3. Inversie via Back-Projectie: Door gebruik te maken van de afgeleide Fourier-snede-identiteit, keert de auteur de transformatie om door de inverse Fourier- en Hankel-transformaties toe te passen. Dit proces leidt tot expliciete reconstructieformules van het type gefilterde back-projectie. De "filtering"-stap wordt geïdentificeerd als de toepassing van de Riesz-potentiaal $I^{(1-d)}$ met betrekking tot het eerste argument.

Belangrijkste Bijdragen en Resultaten
Het artikel biedt twee primaire expliciete reconstructieformules voor het herstellen van $f(x, y)$ uit de gegevens $R^{(p)}$:

1. Angulaire Parametrisatie (Theorem 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$$
   Deze formule integreert over alle kegels die door het reconstructiepunt $(x, y)$ gaan, geparametriseerd door de openingshoek $\theta$ en de richting $n$ op de eenheidsbol.

2. Ruimtelijke Parametrisatie (Theorem 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$$
   Deze vorm integreert over het basisvlak $\mathbb{R}^{d-1}$ en somt de bijdragen op van kegels waarvan de hoekpunten liggen bij $u$.

Specifieke Dimensies:

• Voor $d=2$: De Riesz-potentiaal $I^{-1}$ komt overeen met de compositie van de afgeleide en de Hilbert-transformatie ($\partial_u H_u$).
• Voor $d=3$: De Riesz-potentiaal $I^{-2}$ komt overeen met de negatieve Laplaciaan ($-\Delta_u$). Bijgevolg omvat de 3D-reconstructieformule het toepassen van de Laplaciaan op de gegevens, gevolgd door een kegelvormige back-projectie.

Betekenis en Claims
Het artikel beweert dat dit de eerste expliciete reconstructieformules van het type gefilterde back-projectie zijn voor de kegelvormige Radon-transformatie in willekeurige dimensies $d \geq 2$ en voor elke reële gewicht $p$.

• Lokaliteit: De auteur merkt een onderscheid aan op basis van de pariteit van de dimensie $d$. Voor oneven dimensies is de reconstructie lokaal, wat betekent dat de waarde van $f$ in een punt $(x, y)$ alleen afhangt van integralen over kegels die door een willekeurig kleine omgeving van dat punt gaan. Voor even dimensies is de reconstructie niet-lokaal, waarbij kennis nodig is van integralen over alle kegelvormige oppervlakken. Dit gedrag weerspiegelt dat van de klassieke Radon-transformatie.
• Nieuwigheid: Hoewel equivalente formules voor $d=2$ eerder in andere werken zijn afgeleid, worden de formules voor $d=3$ en hogere dimensies als nieuwe bijdragen gepresenteerd.
• Praktische bruikbaarheid: De auteur stelt dat de 3D-reconstructieformules bijzonder eenvoudig zijn, bestaande uit een toepassing van de Laplaciaan gevolgd door kegelvormige back-projectie. De auteur suggereert dat de numerieke implementatie eenvoudig zou moeten zijn, volgens de gevestigde methoden voor klassieke en sferische Radon-transformaties, hoewel numerieke studies aan toekomstig onderzoek zijn overgelaten.

Het werk wordt gekaderd als een theoretische fundering voor emissietomografie met behulp van Compton-camera's, waarbij de noodzakelijke wiskundige instrumenten worden geleverd om de gegevens gegenereerd door dergelijke apparaten te inverteren.