Microlocal analysis of a non-linear cone transform and applications to Compton camera imaging

Resumen Técnico: Análisis Microlocal de una Transformada de Cono No Lineal y Aplicaciones a la Imagen por Cámara Compton

Planteamiento del Problema
Este artículo aborda el problema inverso en la imagen por cámara Compton, específicamente la recuperación simultánea de la intensidad de la fuente $f: \mathbb{R}^n \to \mathbb{R}$ y del coeficiente de atenuación espacialmente variable $\mu: \mathbb{R}^n \to \mathbb{R}$. A diferencia de la literatura previa, que trata mayoritariamente las transformadas de cono como operadores lineales (asumiendo que $\mu = 0$), este trabajo analiza un modelo no lineal donde se tiene en cuenta la atenuación de los rayos. Los datos $h$ se modelan como la integral de $f$ sobre conos, ponderada por un término de atenuación exponencial $g = \exp(-G\mu)$, donde $G$ es una transformada de haz divergente. La no linealidad surge porque la atenuación depende de la integral de trayectoria de $\mu$, lo que hace que la transformada sea no lineal en $\mu$. Los autores pretenden caracterizar las singularidades introducidas por esta no linealidad, determinar la unicidad de $f$ y $\mu$, y proponer un método de reconstrucción que mitigue los artefactos causados por la interacción entre la fuente y el medio de atenuación.

Metodología
Los autores emplean el análisis microlocal, específicamente la teoría de distribuciones conormales y los Operadores Integrales de Fourier (FIO), para estudiar las singularidades de los datos. La metodología avanza a través de varias etapas clave:

1. Transformación a Integrales de Rayos Ponderadas:
   Utilizando argumentos geométricos similares al trabajo previo sobre transformadas de cono lineales, los autores convierten primero los datos de la integral de cono $R(f, \mu)$ en datos de integral de rayos ponderados, denotados como $\tilde{h} = D_k(f, \mu)$. Bajo condiciones geométricas específicas, se demuestra que estos datos se factorizan como:
   $$\tilde{h} = \exp(-X_{w_1}\mu) X_{w_2}f$$
   donde $X_w$ representa una transformada de rayos ponderada, y $w_1, w_2$ son pesos suaves. Esta formulación separa el problema en un producto de distribuciones.

2. Análisis de Singularidades Microlocales:
   Los autores analizan las singularidades del producto $\exp(-X_{w_1}\mu) X_{w_2}f$ utilizando la teoría de distribuciones conormales (específicamente las clases $I^{p,l}$ asociadas con variedades lagrangianas de intersección limpia).

   • Tangencia Simple: Cuando un rayo es tangente a los bordes del soporte de $f$ o $\mu$ en un único punto estrictamente convexo, los datos exhiben singularidades características de la transformada de rayos lineal.
   • Tangencia Doble: Cuando un rayo es tangente a los bordes en dos puntos distintos (por ejemplo, tangente tanto al borde de la fuente $\partial\Omega$ como al soporte de la atenuación $\Sigma$), el producto de las distribuciones genera nuevas singularidades. Estos artefactos aparecen en direcciones determinadas por la suma de los vectores normales en los dos puntos de tangencia, pero no están presentes en la transformada lineal.

3. Estrategia de Reconstrucción:
   Los autores proponen un algoritmo de reconstrucción basado en la regularidad de Sobolev de estas singularidades. Demuestran que las singularidades correspondientes al verdadero borde de la fuente $\partial\Omega$ son más fuertes (en términos de orden de Sobolev) que los artefactos generados por la interacción no lineal.

   • Recuperación del Envolvente Convexo: Primero se determina el envolvente convexo de la fuente, $\Omega_c$, identificando los ceros en los datos (rayos que no intersectan la fuente).
   • Supresión de Artefactos: Se aplica una función de corte suave $\psi$, con soporte en $\Omega_c$ y que se anula en el soporte de $\mu$, a la imagen reconstruida. Esto elimina las singularidades asociadas con $\mu$ y los artefactos no lineales, aislando las singularidades más fuertes de $\partial\Omega$.
   • Pruebas de Unicidad: Una vez recuperada $f$, los autores utilizan una versión generalizada del teorema de unicidad de Natterer para la transformada de rayos atenuada para demostrar que $\mu$ está determinado de forma única por los datos, siempre que se cumplan ciertas condiciones geométricas (por ejemplo, la existencia de rayos desde el borde de la fuente que no intersectan $\mu$).

Contribuciones Clave y Resultados

• Caracterización Microlocal de la No Linealidad: El artículo proporciona el primer análisis microlocal de una transformada de cono no lineal en este contexto. Cuantifica rigurosamente los artefactos introducidos por la atenuación, mostrando que se manifiestan como distribuciones conormales de un orden de Sobolev inferior ($H^\alpha$ para $\alpha < 3/2$) comparado con las singularidades reales de la fuente ($H^\alpha$ para $\alpha < 1/2$).
• Recuperación Única de la Fuente y la Atenuación: Bajo el supuesto de que los soportes de $f$ y $\mu$ son disjuntos y que la fuente está contenida dentro del envolvente convexo de la superficie del detector, los autores demuestran que tanto $f$ como $\mu$ están determinados de forma única por los datos.
• Distinción de Estabilidad: El análisis revela que, si bien $f$ puede recuperarse de forma estable (específicamente sus bordes), la recuperación de $\mu$ es altamente inestable. El problema para $\mu$ se identifica como un problema de tomografía exterior con datos limitados, lo que conduce a bordes difusos y características no resueltas en la reconstrucción, de forma similar a la tomografía de ángulo limitado.
• Validación mediante Simulación: Los autores presentan reconstrucciones simuladas en 2D ($n=2$) con ruido gaussiano añadido. Los resultados confirman que:
  • El envolvente convexo de la fuente puede identificarse con precisión a partir de los ceros de los datos.
  • El método de filtrado propuesto elimina con éxito los artefactos de rayas (streaking) causados por las tangencias dobles.
  • La intensidad de la fuente $f$ se recupera con alta precisión; el valor de la intensidad reconstruida tiene un error del 6% (1.06 frente a 1.0 real), mientras que la reconstrucción de la forma ($f_c$) logra un error de mínimos cuadrados relativo del 21%.
  • El coeficiente de atenuación $\mu$ solo se recupera parcialmente; solo se resuelven los bordes visibles para los datos, con una difusión significativa debido a la naturaleza mal planteada del problema exterior.

Significancia y Reivindicaciones
Los autores afirman que este trabajo hace avanzar el campo al ir más allá de las aproximaciones lineales en la imagen por cámara Compton. Al modelar explícitamente la atenuación, identifican y caracterizan artefactos específicos que de otro modo degradarían la calidad de la imagen. La principal significancia reside en el marco teórico que permite la separación de las singularidades de la fuente de los artefactos inducidos por la atenuación utilizando la regularidad de Sobolev y las restricciones geométricas.

El artículo señala modestamente que, aunque la unicidad teórica de $\mu$ está establecida, la recuperación práctica está limitada por el mal planteamiento del problema. Los autores no pretenden resolver el caso general de atenuación arbitraria, sino que proporcionan un método para escenarios donde la fuente y los materiales de atenuación están espacialmente separados (soportes disjuntos) y la geometría permite una cobertura de rayos suficiente. Reconocen que la recuperación de $\mu$ sigue siendo inestable y que se requiere trabajo futuro para abordar situaciones con menos datos (por ejemplo, vértices de cono en una curva en lugar de una superficie). El trabajo valida la utilidad del análisis microlocal para comprender y mitigar los efectos no lineales en las modalidades de imagen.