Fast Eikonal Phase Retrieval for High-Throughput Beamlines

Technische Samenvatting: Snelle Eikonal Faseherstel voor High-Throughput Beamlines

Probleemstelling
Propagatie-gebaseerde fasecontrast micro-tomografie (PPC-µCT) bij vierde-generatie synchrotronbronnen maakt hooggevoelige beeldvorming van zwak brekende kenmerken mogelijk door gebruik te maken van lange afstand tussen monster en detector. Echter, dit regime introduceert aanzienlijke uitdagingen voor beeldreconstructie. Standaard faseherstelmethoden met één afstand (bijv. Paganin-filters) vertrouwen op gelineraliseerde transport-van-intensiteit aannames. In heterogene specimens die zowel zwakke zachte weefsels als sterke interfaces bevatten (bijv. bot–lucht of bot–vloeistof), genereren steile fasegradiënten sterke intensiteitsmodulaties die deze lineaire aannames schenden. Dit leidt tot karakteristieke niet-lineaire artefacten, zoals langdurige strepen (streaks), die de beeldkwaliteit verslechteren en downstream analyse compromitteren.

Hoewel volledige Fresnel wave-optics formuleringen deze effecten kunnen modelleren, vereisen ze een expliciete pixel-gewijze representatie van het complexe golfveld. Dit legt strikte bemonsteringsbeperkingen op; wanneer fasevariaties snel zijn (zoals veroorzaakt door scherpe interfaces), wordt standaard Fresnel-propagatie computationeel prohibitief vanwege aliasing en enorme geheugenvereisten. Daarentegen leed het eerder geïntroduceerde Eikonal Phase Retrieval (EPR) algoritme aan hoge computationele kosten (tientallen uren voor grote datasets), wat incompatibel was met de hoge throughput-acquisitieratens van moderne beamlines.

Methodologie
De auteurs introduceren een "Fast Eikonal Phase Retrieval" (Fast EPR) formulering die het oorspronkelijke EPR-algoritme met meer dan twee ordes van grootte versnelt terwijl de gecontroleerde nauwkeurigheid behouden blijft. De methodologie is gebaseerd op een tweede-orde asymptotische expansie van de Fresnel-propagator in de propagatiedistantie $L$, aangevuld met de leidende Wentzel–Kramers–Brillouin (WKB) wave-optics correctie.

De aanpak maakt gebruik van een verenigd framework met twee complementaire solvers om verschillende fysieke regimes te behandelen:

1. Lokale Solver (Sub-pixel verschuivingen): Voor regimes waar de eikonal-geïnduceerde transversale straalverschuiving $\Delta(r)$ kleiner is dan één pixel ($\eta < 1$), gebruikt de methode een lokale $O(L^2)$ sluiting. Dit omvat een tweede-orde differentiële expansie van de intensiteitstransport, inclusend de kwadratische eikonal term (WKB0) en de leidende diffractie-druk term (WKB1). Dit lokale model wordt gepreconditioneerd door FFT-diagonale, energie-afhankelijke inverse operatoren (Paganin-type filters) om convergentie te versnellen.
2. Niet-lokale Solver (Multi-pixel verschuivingen): Wanneer verschuivingen groter zijn dan één pixel ($\eta > 1$), faalt de lokale Taylor-expansie. De methode schakelt over naar een niet-lokale formulering waarbij intensiteit wordt gepropageerd via een expliciete eikonal straalmapping. Dit wordt geïmplementeerd als een massabehoudende bilineaire herverdeling op het detectorrooster. De detectorresiduen worden teruggetransfereerd naar het objectrooster met de bijbehorende adjoint (transponeer) operator, geïmplementeerd als bilineaire interpolatie, voordat de FFT-diagonale preconditioner wordt toegepast.

Polychromatische Ondersteuning
Het framework ondersteunt polychromatische data via een compacte spectrale discretisatie. Het spectrum wordt gemodelleerd als een incoherente som van $N_s$ spectrale componenten. Alle per-lijn grootheden (fase/absorptie parameters en verschuivingen) worden energie-afhankelijk behandeld, waardoor energie-afhankelijke transport en inversie mogelijk is terwijl FFT-efficiëntie behouden blijft.

Belangrijkste Bijdragen

• Algoritmische Versnelling: De nieuwe formulering vermindert de computationele kosten van faseherstel met ongeveer 580 keer (2,8 ordes van grootte) vergeleken met de oorspronkelijke EPR-implementatie wanneer genormaliseerd naar dezelfde GPU-bronnen.
• Hybride Solver Strategie: De introductie van een dual-solver aanpak (lokale $O(L^2)$ voor sub-pixel verschuivingen en niet-lokale straalmapping voor multi-pixel verschuivingen) zorgt voor robuustheid over een breed scala aan experimentele condities, van zwak brekende zachte weefsels tot sterke interfaces.
• Tweede-Orde Fysica: Door de volledige $O(L^2)$ inhoud van de WKB-expansie te behouden (inclusief de diffractie-druk term), onderdrukt de methode niet-lineaire streep-artefacten die persistent zijn in eerste-orde lineaire modellen.
• Polychromatische Efficiëntie: De methode maakt nauwkeurige reconstructie van polychromatische data mogelijk zonder de computationele efficiëntie op te offeren die vereist is voor high-throughput workflows.

Resultaten
De auteurs hebben de methode gevalideerd met behulp van een sheep-head HiP-CT dataset en een bamboo sample stress test:

• Artefact Reductie: Vergelijkingen op de sheep-head dataset toonden aan dat het tweede-orde ($L^2$) model de resterende niet-lineaire propagatie-artefacten aanzienlijk vermindert vergeleken met eerste-orde ($L^1$) modellen, wat zorgt voor een schonere weergave van zacht-weefsel contrast nabij sterke gradiënten.
• Convergentie: Het algoritme convergeert snel; het merendeel van de verbetering wordt bereikt in de eerste iteratie zodra de tweede-orde fysica geactiveerd zijn.
• Robuustheid: Op de bamboo sample, waar straalverschuivingen de detector pixelgrootte overschreden, faalde de lokale solver (divergerend na de eerste iteratie), terwijl de niet-lokale solver stabiel en accuraat bleef.
• Prestaties: Voor een dataset van 8.000 radiografieën ($256 \times 3104$ pixels), voltooide de Fast EPR methode de faseherstel in ongeveer 1,1 tot 1,9 minuten op een enkele GPU (voor monochromatische data) en 1,47 tot 1,9 minuten voor polychromatische data (5 spectrale lijnen). Dit staat in contrast met de oorspronkelijke EPR benchmark, die ongeveer 870 minuten vereiste voor een vergelijkbaar volume.

Significantie
Het artikel stelt dat deze verenigde aanpak nauwkeurige en computationeel efficiënte faseherstel mogelijk maakt voor propagatiecondities die relevant zijn voor high-throughput PPC-µCT experimenten. Door de reconstructietijden te reduceren tot compatibiliteit met beamline acquisitieratens (één monster elke paar uur), faciliteert de methode het routinematig gebruik van niet-lineaire faseherstel in high-throughput workflows. Bovendien maakt de beschikbaarheid van een snelle tweede-orde near-field model met expliciete wave-optics correctie onderzoekers in staat om dichter bij de near-field grens te opereren en regimes met sterkere fasegradiënten te verkennen (bijv. lagere energie of hogere resolutie) waar standaard lineaire benaderingen falen, zonder de prohibitieve kosten van volledige Fresnel-propagatie te dragen. De auteurs merken op dat de algoritmische complexiteit van de nieuwe methode verborgen kan worden achter tomografische back-projectie in concurrente hardware-opstellingen, wat verdere integratie in praktische reconstructie pipelines mogelijk maakt.