Exact reconstruction formulas for a Radon transform over cones

Resumo Técnico: Fórmulas de Reconstrução Exata para uma Transformada de Radon sobre Cones

Enunciado do Problema
Este artigo aborda a inversão matemática de uma variante específica da transformada de Radon conhecida como transformada de Radon cônica, denotada por $R^{(p)}$. A transformada integra uma função suave $f$ suportada no semiespaço superior $\mathbb{R}^{d-1} \times (0, \infty)$ sobre superfícies cônicas unilaterais. Estes cones são caracterizados por vértices localizados no hiperplano $\mathbb{R}^{d-1} \times \{0\}$, um eixo de simetria ortogonal a este plano (alinhado com o eixo $y$), e um ângulo de abertura de meia-abertura $\theta \in (0, \pi/2)$. A integração inclui um peso radial $1/s^p$, onde $s$ representa a distância ao longo do gerador do cone.

A principal motivação para o estudo desta transformada reside na tomografia computadorizada (TC) que considera o espalhamento Compton. Em tais modalidades de imagem, a detecção de fótons espalhados permite rastrear a localização da fonte de volta à superfície de um cone, em vez de uma linha ou plano. Embora existam fórmulas de reconstrução de Fourier para casos específicos em três dimensões ($p \in \{0, 2\}$) e fórmulas de retroprojeção existam para duas dimensões, o caso geral para dimensões espaciais arbitrárias $d \geq 2$ e para um parâmetro de peso arbitrário $p \in \mathbb{R}$ não havia sido estudado anteriormente.

Metodologia
O autor, Markus Haltelier, emprega técnicas de análise harmônica, especificamente aproveitando a transformada de Fourier e a transformada de Hankel, para derivar as fórmulas de inversão. O núcleo da metodologia envolve estabelecer uma "identidade de fatia de Fourier" para a transformada de Radon cônica, análoga à bem conhecida teorema da fatia de Fourier para a transformada de Radon clássica.

A derivação procede em três etapas lógicas:

1. Redução ao Caso Não Ponderado: Um lema é estabelecido mostrando que a transformada $R^{(p)}$ para um peso $p$ arbitrário pode ser relacionada à transformada não ponderada $R^{(0)}$ via operadores multiplicativos simples envolvendo $\cos(\theta)$ e $y^{-p}$. Isso permite que os autores foquem sua derivação no caso $p=0$ e generalizem os resultados posteriormente.
2. Derivação da Identidade de Fatia de Fourier: Ao aplicar a transformada de Fourier em relação à variável espacial $u \in \mathbb{R}^{d-1}$ e utilizar a representação integral de funções de Bessel sobre a esfera $S^{d-2}$, o autor deriva uma identidade que liga a transformada de Fourier-Hankel da função original $f$ à transformada de Fourier dos dados da transformada de Radon cônica.
3. Inversão via Retroprojeção: Utilizando a identidade de fatia de Fourier derivada, o autor inverte a transformada aplicando as transformadas de Fourier e Hankel inversas. Este processo leva a fórmulas de reconstrução explícitas do tipo retroprojeção filtrada. A etapa de "filtragem" é identificada como a aplicação do potencial de Riesz $I^{(1-d)}$ em relação ao primeiro argumento.

Principais Contribuições e Resultados
O artigo fornece duas fórmulas de reconstrução explícitas para recuperar $f(x, y)$ a partir dos dados $R^{(p)}$:

1. Parametrização Angular (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$$
   Esta fórmula integra sobre todos os cones que passam pelo ponto de reconstrução $(x, y)$, parametrizados pelo ângulo de abertura $\theta$ e pela direção $n$ no círculo unitário.

2. Parametrização Espacial (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$$
   Esta forma integra sobre o plano de base $\mathbb{R}^{d-1}$, somando contribuições de cones cujos vértices estão em $u$.

Dimensões Específicas:

• Para $d=2$: O potencial de Riesz $I^{-1}$ corresponde à composição da derivada e da transformada de Hilbert ($\partial_u H_u$).
• Para $d=3$: O potencial de Riesz $I^{-2}$ corresponde ao Laplaciano negativo ($-\Delta_u$). Consequentemente, a fórmula de reconstrução 3D envolve a aplicação do Laplaciano aos dados seguida de uma retroprojeção cônica.

Significância e Alegações
O artigo afirma que estas são as primeiras fórmulas de reconstrução explícitas do tipo retroprojeção filtrada para a transformada de Radon cônica em dimensões arbitrárias $d \geq 2$ e para qualquer peso real $p$.

• Localidade: O autor nota uma distinção baseada na paridade da dimensão $d$. Para dimensões ímpares, a reconstrução é local, o que significa que o valor de $f$ em um ponto $(x, y)$ depende apenas de integrais sobre cones que passam por uma vizinhança arbitrariamente pequena de esse ponto. Para dimensões pares, a reconstrução é não-local, exigindo o conhecimento de integrais sobre todas as superfícies cônicas. Este comportamento espelha o da transformada de Radon clássica.
• Novidade: Embora fórmulas equivalentes para $d=2$ tenham sido derivadas anteriormente em outros trabalhos, as fórmulas para $d=3$ e dimensões superiores são apresentadas como novas contribuições.
• Praticidade: O artigo afirma que as fórmulas de reconstrução 3D são particularmente simples, consistindo na aplicação do Laplaciano seguida de retroprojeção cônica. O autor sugere que a implementação numérica deve ser direta, seguindo métodos estabelecidos para as transformadas de Radon clássica e esférica, embora estudos numéricos sejam deixados para pesquisas futuras.

O trabalho é estruturado como uma fundamentação teórica para a tomografia de emissão usando câmeras Compton, fornecendo as ferramentas matemáticas necessárias para inverter os dados gerados por tais dispositivos.