Convergence of Ray- and Pixel-Driven Discretization Frameworks in the Strong Operator Topology

1. Problématique

La tomographie par ordinateur (CT) repose sur la transformée de Radon pour modéliser le processus de mesure (intégrales de lignes) et sa transposée, la rétroprojection, pour la reconstruction. En pratique, ces opérateurs infinis doivent être discrétisés pour être traités numériquement.

Le problème central abordé dans cet article est l'absence de justification théorique rigoureuse concernant les erreurs d'approximation des schémas de discrétisation les plus couramment utilisés en pratique :

• Approche « Ray-Driven » (pilotée par les rayons) pour la transformée de Radon (avant-projection).
• Approche « Pixel-Driven » (pilotée par les pixels) pour la rétroprojection.

Cette combinaison hybride (souvent notée rd-pd*) est anecdotiquement considérée comme offrant les meilleures performances d'approximation, contrairement aux paires « appariées » (rd-rd* ou pd-pd*). Cependant, les paires non appariées posent un problème théorique : elles ne sont pas strictement adjointes l'une de l'autre, ce qui peut théoriquement nuire à la convergence des solveurs itératifs. De plus, la plupart des analyses existantes supposent des résolutions déséquilibrées (ex: résolution spatiale $\delta_x \ll$ résolution détecteur $\delta_s$), alors que la pratique courante utilise des résolutions équilibrées ($\delta_x \approx \delta_s$).

L'objectif est de prouver rigoureusement la convergence de ces méthodes hybrides sous des résolutions équilibrées.

2. Méthodologie

L'auteur propose une nouvelle interprétation mathématique des schémas de discrétisation pour permettre une analyse structurée :

• Interprétation Convolutionnelle et Opérateurs de Rang Fin : Les discrétisations sont reformulées comme des opérateurs de rang fini agissant sur des espaces $L^2$. Cela permet de définir des fonctions de poids (weight functions) explicites :
  • $\omega^{rd}_{\delta_x}$ : Fonction de poids pour la méthode ray-driven (basée sur la longueur d'intersection rayon-pixel).
  • $\omega^{pd}_{\delta_s}$ : Fonction de poids pour la méthode pixel-driven (basée sur l'interpolation linéaire / anterpolation).
• Analyse de Convergence : L'étude se concentre sur la convergence dans la topologie de l'opérateur fort (Strong Operator Topology - SOT). Cela signifie que pour toute fonction $f$ (ou $g$) donnée, la suite des opérateurs discrétisés converge vers l'opérateur continu lorsque les paramètres de discrétisation ($\delta_x, \delta_\phi, \delta_s$) tendent vers zéro.
• Hypothèses de Résolution : L'analyse distingue deux cas :
  • Résolutions équilibrées : $\delta_x \approx \delta_s$.
  • Résolutions déséquilibrées : $\delta_s \ll \delta_x$ ou $\delta_x \ll \delta_s$.
• Preuves Théoriques : L'auteur utilise des inégalités de Taylor, des estimations de sommes de Riemann et le principe de l'unicité bornée (Uniform Boundedness Principle) pour démontrer la convergence. Des lemmes techniques établissent les propriétés des sommes de poids (exactitude pour la méthode ray-driven sur les centres de pixels, interpolation pour la méthode pixel-driven).

3. Contributions Clés

1. Justification Théorique de l'Approche Hybride (rd-pd*) : Le résultat principal (Théorème 3.1) prouve que l'approche combinée ray-driven (avant) et pixel-driven (arrière) converge dans la topologie de l'opérateur fort sous des résolutions équilibrées ($\delta_x \approx \delta_s$).
2. Analyse des Cas Déséquilibrés :
   • La rétroprojection ray-driven converge uniquement si $\delta_s \ll \delta_x$ (résolution détecteur beaucoup plus fine que l'espace).
   • La transformée de Radon pixel-driven converge uniquement si $\delta_x \ll \delta_s$.
3. Démonstration de la Non-Convergence des Paires Appariées en Équilibre : Les résultats implicites montrent que les paires appariées (rd-rd* ou pd-pd*) ne convergent pas nécessairement dans la topologie forte sous des résolutions équilibrées, ce qui remet en question leur utilisation théorique dans ce contexte, malgré leur popularité pour des raisons de cohérence algébrique.
4. Extension aux Angles Limités et Épars : Les résultats de convergence sont étendus aux cas de tomographie à angles limités et aux projections sur des ensembles d'angles discrets.

4. Résultats

• Théorème de Convergence (SOT) :
  • Pour toute fonction $f \in L^2(\Omega)$, $\|R f - R^{rd}_{\delta_n} f\| \to 0$ si $\delta_s / \delta_x \leq c$ (cas équilibré).
  • Pour toute fonction $g \in L^2(S)$, $\|R^* g - R^{pd*}_{\delta_n} g\| \to 0$ si $\delta_x / \delta_s \leq c$ (cas équilibré).
  • Cela valide l'utilisation de la combinaison ray-driven / pixel-driven pour des résolutions équilibrées.
• Comportement des Erreurs :
  • La méthode pixel-driven pour la transformée de Radon présente des erreurs importantes (oscillations) pour certains angles spécifiques (multiples de $\pi/4$) lorsque les résolutions sont équilibrées, mais ces erreurs diminuent si $\delta_x \ll \delta_s$.
  • La méthode ray-driven pour la rétroprojection ne converge pas vers la valeur exacte pour chaque angle individuel sous résolutions équilibrées, mais converge en moyenne (sur l'ensemble des angles).
• Expériences Numériques :
  • Utilisation du fantôme FORBILD (tête) pour valider les théories.
  • Confirmation que l'erreur pixel-driven stagne sous résolutions équilibrées pour certains angles, tandis que l'erreur ray-driven diminue avec la résolution.
  • Pour la rétroprojection, l'erreur ray-driven reste constante si le rapport $\delta_s/\delta_x$ est fixe (même si la résolution globale augmente), mais diminue lorsque $\delta_s/\delta_x \to 0$.

5. Signification et Impact

Cet article comble un fossé critique entre la pratique numérique et la théorie mathématique en tomographie :

• Validation Pratique : Il fournit une base théorique solide à l'usage empirique répandu de la combinaison ray-driven / pixel-driven dans les logiciels de reconstruction CT modernes, confirmant qu'elle est mathématiquement justifiée même avec des résolutions équilibrées.
• Choix de Stratégie : Il suggère que l'utilisation d'opérateurs non appariés mais convergents (rd-pd*) est préférable à l'utilisation d'opérateurs appariés mais potentiellement non convergents (rd-rd* ou pd-pd*) dans le cadre des résolutions équilibrées.
• Limites et Perspectives : L'article met en lumière que la convergence en norme d'opérateur (plus forte que la convergence SOT) n'est pas garantie pour les résolutions équilibrées, ce qui pourrait expliquer certaines instabilités dans les solveurs itératifs. Il ouvre la voie à des recherches futures sur les transformées en éventail (fanbeam) et conique (conbeam) et sur l'analyse de la convergence en norme d'opérateur dans des régimes de résolution déséquilibrée.

En résumé, ce travail démontre que l'approche hybride, bien que mathématiquement "non appariée", est la seule à garantir la convergence de l'approximation de la transformée de Radon et de sa rétroprojection dans les conditions de résolution standard de l'industrie médicale.