A Reciprocity-Law-Compliant Photoacoustic Forward-Adjoint Operator

Titel

Een aan de wederkerigheidswet (reciprocity law) gehoorzaamende fotoakoestische voorwaartse-geadjungeerde operator.

1. Probleemstelling

Fotoakoestische tomografie (PAT) is een hybride beeldvormingstechniek die de hoge contrasten van optische absorptie combineert met de hoge ruimtelijke resolutie van akoestische detectie. Het inverse probleem in PAT bestaat uit het reconstrueren van de initiële drukverdeling (die correleert met de optische absorptie) op basis van tijdsreeksen van drukmetingen aan de rand van het object.

Een cruciaal aspect bij het oplossen van dit inverse probleem, vooral binnen iteratieve minimaliseringsframeworks, is de nauwkeurige definitie van de voorwaartse operator (die de initiële druk naar de randmetingen mapt) en de bijbehorende geadjungeerde operator.

• Uitdaging: Bestaande methoden maken vaak aannames over de gladheid van het randoppervlak of behandelen de restrictie van het veld naar de rand als een enkelvoudige potentiaal. Dit kan leiden tot inconsistenties met de fundamentele fysica, specifiek de wederkerigheidswet (reciprocity law).
• Specifieke moeilijkheid: De fotoakoestische bron is instantane (een Dirac-delta in de tijd) en de restrictie van het drukveld naar de meetrand introduceert singulariteiten. Het is essentieel dat de voorwaartse en geadjungeerde operatoren een consistente inproduct-relatie (inner-product relation) voldoen om convergentie en stabiliteit in reconstructie-algoritmen te garanderen.

2. Methodologie

De auteur breidt een eerder ontwikkeld framework voor akoestische inverse problemen uit naar PAT. De kern van de methodologie bestaat uit drie stappen:

A. Definitie van de Voorwaartse Operator

De totale PAT voorwaartse operator wordt gedefinieerd als een samengestelde mapping ($A \circ Q$):

1. Bronoperator ($Q$): Deze mappert de initiële drukverdeling $p_0$ naar een tijds-geregulariseerde bron $s_Q$. Omdat de fysieke bron instantane is, wordt een geregulariseerde Dirac-delta distributie ($\delta_q(t)$) gebruikt in de tijd om singulariteiten te voorkomen.
2. Akoestische Voorwaartse Operator ($A$): Deze beschrijft de golfvoortplanting in de ruimte. Het cruciale component is de ontvangstoperator (reception operator) $R^D_{G,r}$.
   • Deze operator mappert het drukveld in de vrije ruimte naar de rand (ontvanger).
   • In plaats van een directe restrictie, wordt de druk op de rand geïnterpreteerd als een "uitgesmeerde singulariteit" in de normale afgeleide van het veld over de grens.
   • Er wordt een geschaalde variant van de ontvangstoperator gebruikt, die een geometrische schalingsfactor ($a_R$) bevat om afhankelijkheid van de regularisatiebandbreedte te elimineren.

B. Afleiding van de Geadjungeerde Operator

Om de geadjungeerde operator te vinden, wordt gebruikgemaakt van de inproduct-relatie:

• De geadjungeerde van de bronoperator ($Q^*$) wordt afgeleid via partiële integratie in de tijd. Dit resulteert in een operator die de tijdsafgeleide van het veld op $t=0$ neemt.
• De geadjungeerde van de akoestische operator ($A^*$) correspondeert met een tijd-omgekeerde dubbel-laag potentiaal integraal (time-reversed double-layer potential).
• De volledige geadjungeerde operator combineert deze twee componenten. Het resultaat is een operator die voldoet aan de wederkerigheidswet: de emissie-operator (van rand naar ruimte) en de ontvangst-operator (van ruimte naar rand) zijn elkaars linksinversen in een distributieve zin.

C. Discretisatie en Numerieke Implementatie

• De operatoren worden geïmplementeerd op een gestaggerd rooster (staggered grid) in ruimte en tijd.
• Ruimtelijke discretisatie: Een k-ruimte pseudospectrale methode wordt gebruikt voor de golfvergelijking. Perfect Matched Layers (PML) worden toegepast om reflecties aan de randen van het rekenrooster te voorkomen.
• Regularisatie: De Dirac-delta distributies ($\delta_b$) voor de ruimtelijke restrictie worden geregulariseerd met een parameter $b$ (vaak gelijk aan de roosterafstand $\Delta x$). Waarden onder een bepaalde drempel worden genegeerd om de rekentijd te beperken.
• Minimalisatie: Het inverse probleem wordt opgelost via een projectie-gradientafdalingsalgoritme (projected gradient descent) met een niet-negativiteitsbeperking ($p_0 \geq 0$).

3. Belangrijkste Bijdragen

1. Wederkerigheidswet-Compliant Framework: De paper introduceert een voorwaartse-geadjungeerde operatorpaar dat strikt voldoet aan de fysica van wederkerigheid. De ontvangstoperator fungeert als linksinvers van de emissie-operator, wat fundamenteel is voor de stabiliteit van iteratieve reconstructies.
2. Regularisatie van Singulariteiten: Er wordt een systematische aanpak geboden voor het regulariseren van zowel de instantane bron (tijds-domein) als de restrictie naar de rand (ruimtelijk domein) met behulp van geregulariseerde Dirac-delta distributies.
3. Geschaalde Operator: De introductie van een schalingsfactor ($1/a_R$) in de operator definities elimineert de afhankelijkheid van de bandbreedte-parameter van de regularisatie, wat leidt tot een unieke en robuuste formulering.
4. Numerieke Validatie: De consistentie van de afgeleide operatoren wordt wiskundig bewezen en numeriek geverifieerd via een inproduct-test (inner-product test) op zowel een rooster-gealigneerde ("on-grid") als een niet-rooster-gealigneerde ("off-grid") configuratie.

4. Resultaten

De auteurs hebben de methode getest in twee scenario's:

• Validatie van de Inproduct-Relatie:

  • On-grid scenario: Ontvangers liggen exact op de roosterpunten. De fout in de inproduct-relatie was extreem laag ($9.32 \times 10^{-4}\%$).
  • Off-grid scenario: Ontvangers liggen op een cirkel die niet samenvalt met de roosterpunten (realistischer scenario). Ook hier was de consistentie zeer hoog, met een fout van $2.07 \times 10^{-5}\%$. Dit bewijst dat de methode robuust is voor niet-gealigneerde meetopstellingen.

• Beeldreconstructie:

  • Een simulatie werd uitgevoerd met een "vessel phantom" (vaatjes-phantom) als initiële drukbron.
  • Om "inverse crime" (het gebruik van exact dezelfde discretisatie voor simulatie en reconstructie) te voorkomen, werden verschillende roosterafstanden en tijdstappen gebruikt voor de simulatie van de data en de reconstructie.
  • Conclusie: Na 25 iteraties convergerde het algoritme monotoon. De gereconstrueerde afbeelding was zowel kwalitatief (visuele overeenkomst) als kwantitatief (lage relatieve fout) nauwkeurig in vergelijking met de grondwahrheid (ground truth).

5. Betekenis en Conclusie

Deze studie biedt een wiskundig rigoureuze en fysiek correcte basis voor iteratieve reconstructie in fotoakoestische tomografie.

• Theoretische impact: Het lost het probleem op van inconsistenties tussen voorwaartse en geadjungeerde operatoren die vaak voorkomen in discrete implementaties, door strikt vast te houden aan de wederkerigheidswet.
• Praktische impact: De methode maakt het mogelijk om hoge-kwaliteit beelden te reconstrueren zelfs bij complexe meetopstellingen (off-grid) en met realistische ruis.
• Toekomstperspectief: De auteurs zien de vertaling van dit framework naar klinische toepassingen als de volgende stap, wat kan leiden tot betere diagnostische beeldvorming in de geneeskunde.

Kortom, het artikel presenteert een geavanceerd wiskundig model dat de brug slaat tussen de continue fysica van geluidsgolven en de discrete numerieke implementatie, waarbij de fundamentele wetten van de natuurkunde (wederkerigheid) centraal staan voor een stabiele en accurate beeldreconstructie.