Exact reconstruction formulas for a Radon transform over cones

Resumen Técnico: Fórmulas de Reconstrucción Exacta para una Transformada de Radon sobre Conos

Planteamiento del Problema
Este artículo aborda la inversión matemática de una variante específica de la transformada de Radon conocida como la transformada de Radon cónica, denotada como $R^{(p)}$. La transformada integra una función suave $f$ soportada en el semiespacio superior $\mathbb{R}^{d-1} \times (0, \infty)$ sobre superficies cónicas de un solo lado. Estos conos se caracterizan por vértices ubicados en el hiperplano $\mathbb{R}^{d-1} \times \{0\}$, un eje de simetría ortogonal a este plano (alineado con el eje $y$), y un ángulo de apertura media $\theta \in (0, \pi/2)$. La integración incluye un peso radial $1/s^p$, donde $s$ representa la distancia a lo largo del generador del cono.

La motivación principal para estudiar esta transformada radica en la tomografía computarizada (TC) que tiene en cuenta la dispersión de Compton. En tales modalidades de imagen, la detección de fotones dispersados permite rastrear la ubicación de la fuente hacia la superficie de un cono en lugar de una línea o un plano. Aunque existen fórmulas de reconstrucción de Fourier para casos específicos en tres dimensiones ($p \in \{0, 2\}$) y fórmulas de retroproyección para dos dimensiones, el caso general para una dimensión espacial arbitraria $d \geq 2$ y un parámetro de peso arbitrario $p \in \mathbb{R}$ no había sido estudiado previamente.

Metodología
El autor, Markus Haltmeier, emplea técnicas de análisis armónico, específicamente aprovechando la transformada de Fourier y la transformada de Hankel, para derivar las fórmulas de inversión. El núcleo de la metodología consiste en establecer una "identidad de la rebanada de Fourier" para la transformada de Radon cónica, análoga al bien conocido teorema de la rebanada de Fourier para la transformada de Radon clásica.

La derivación procede en tres pasos lógicos:

1. Reducción al Caso Sin Peso: Se establece un lema que muestra que la transformada $R^{(p)}$ para un peso $p$ arbitrario puede relacionarse con la transformada sin peso $R^{(0)}$ mediante operadores multiplicativos simples que involucran $\cos(\theta)$ y $y^{-p}$. Esto permite a los autores centrarse en la derivación del caso $p=0$ y generalizar los resultados posteriormente.
2. Derivación de la Identidad de la Rebanada de Fourier: Aplicando la transformada de Fourier con respecto a la variable espacial $u \in \mathbb{R}^{d-1}$ y utilizando la representación integral de las funciones de Bessel sobre la esfera $S^{d-2}$, el autor deriva una identidad que vincula la transformada de Fourier-Hankel de la función original $f$ con la transformada de Fourier de los datos de la transformada de Radon cónica.
3. Inversión mediante Retroproyección: Utilizando la identidad de la rebra de Fourier derivada, el autor invierte la transformada aplicando las transformadas inversas de Fourier y de Hankel. Este proceso conduce a fórmulas de reconstrucción explícitas del tipo retroproyección filtrada. El paso de "filtrado" se identifica como la aplicación del potencial de Riesz $I^{(1-d)}$ con respecto al primer argumento.

Contribuciones Clave y Resultados
El artículo proporciona dos fórmulas de reconstrucción explícitas para recuperar $f(x, y)$ a partir de los datos $R^{(p)}f$:

1. Parametrización Angular (Teorema 1.1, Ec. 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$$
   Esta fórmula integra sobre todos los conos que pasan por el punto de reconstrucción $(x, y)$, parametrizados por el ángulo de apertura $\theta$ y la dirección $n$ en la esfera unidad.

2. Parametrización Espacial (Teorema 1.1, Ec. 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$$
   Esta forma integra sobre el plano base $\mathbb{R}^{d-1}$, sumando las contribuciones de los conos cuyos vértices se encuentran en $u$.

Dimensiones Específicas:

• Para $d=2$: El potencial de Riesz $I^{-1}$ corresponde a la composición de la derivada y la transformada de Hilbert ($\partial_u H_u$).
• Para $d=3$: El potencial de Riesz $I^{-2}$ corresponde al Laplaciano negativo ($-\Delta_u$). En consecuencia, la fórmula de reconstrucción 3D implica aplicar el Laplaciano a los datos seguido de una retroproyección cónica.

Significado y Reivindicaciones
El artículo afirma que estas son las primeras fórmulas de reconstrucción explícitas del tipo retroproyección filtrada para la transformada de Radon cónica en dimensiones arbitrarias $d \geq 2$ y para cualquier peso real $p$.

• Localidad: El autor señala una distinción basada en la paridad de la dimensión $d$. Para dimensiones impares, la reconstrucción es local, lo que significa que el valor de $f$ en un punto $(x, y)$ depende solo de las integrales sobre conos que pasan por un vecindad arbitrariamente pequeña de dicho punto. Para dimensiones pares, la reconstrucción es no local, requiriendo el conocimiento de las integrales sobre todas las superficies cónicas. Este comportamiento refleja el de la transformada de Radon clásica.
• Novedad: Si bien se han derivado fórmulas equivalentes para $d=2$ en otros trabajos, las fórmulas para $d=3$ y dimensiones superiores se presentan como nuevas contribuciones.
• Practicidad: El artículo afirma que las fórmulas de reconstrucción 3D son particularmente simples, consistiendo en la aplicación de un Laplaciano seguido de una retroproyección cónica. El autor sugiere que la implementación numérica debería ser sencilla, siguiendo los métodos establecidos para las transformadas de Radon clásica y esférica, aunque los estudios numéricos se dejan para investigaciones futuras.

El trabajo se enmarca en una base teórica para la tomografía de emisión mediante cámaras de Compton, proporcionando las herramientas matemáticas necesarias para invertir los datos generados por tales dispositivos.