Extended Synthetic Validation of Optical Flow-Based Velocity Field Reconstruction for AAA Cine-MRI: Angular Velocity, Noise Robustness, and Dense Flow Comparison

Technische Samenvatting: Uitgebreide Synthetische Validatie van Optical Flow-gebaseerde Snelheidsveldreconstructie voor AAA Cine-MRI

Probleemstelling
De reconstructie van snelheidsvelden uit standaard abdominale aorta-aneurysma (AAA) cine-MRI biedt een pad naar hemodynamische risicostratificatie die verder gaat dan eenvoudige diameters metingen. Hoewel Phase-Contrast en 4D flow MRI referentiestandaarden zijn, vereisen deze gespecialiseerde acquisitie. Standaard cine-MRI, dat routinematig wordt verkregen, kan theoretisch worden geanalyseerd met behulp van optical flow (OF) zonder extra kosten. Echter, een eerdere studie door dezelfde auteur [9], die gebruikmaakte van Pyramidal Lucas–Kanade (PLK) sparse optical flow op een enkele AAA-patiënt, werd beperkt door een gebrek aan vergelijking met dense-flow, de afwezigheid van validatie van hoeksnelheid en de afhankelijkheid van een dataset van één patiënt. Deze uitbreiding adresseert deze hiaten door PLK rigoureus te valideren tegenover dense optical flow-methoden en een synthetische grondwaarheid (ground truth) vast te stellen voor zowel lineaire als hoeksnelheden.

Methodologie
De studie hanteert een veelzijdige validatiebenadering die de reproductie van real-world data, onafhankelijke benchmarks en een nieuwe synthetische fantoom combineert:

1. Databronnen:

   • Real-world Patiëntgegevens: Reproductie van de resultaten uit de oorspronkelijke studie [9] met behulp van short-axis en long-axis cine-MRI van één AAA-patiënt (Rigshospitalet, 2020).
   • Synthetische Schijf Benchmark: Een onafhankelijke benchmark bestaande uit een onregelmatig bewegende, translerende schijf om de gevoeligheid voor venstergrootte ($winSize$) en ruis te testen.
   • Synthetische AAA Fantoom: Een nieuwe, op de literatuur gebaseerde fantoom geconstrueerd met closed-form grondwaarheid voor zowel lineaire als hoeksnelheden. De fantoomparameters (circumferentiële wandrek, voxel-spacing, lumen diameter en hoeksnelheid) zijn afgeleid van een echte, open-access 4D-CTA dataset van 10 AAA-patiënten (Jamshidian et al., 2025). De fantoom simuleert pulsatiele dilatatie en rigide anti-klokse rotatie, waarbij Rician-ruis is toegevoegd om MRI-magnitudebeelden na te bootsen.

2. Algoritmen & Configuraties:

   • Sparse Optical Flow (PLK): Pyramidal Lucas–Kanade met een 4-niveau Gaussische piramide. Een sensitiviteitsmeting van $winSize$ (variërend van 4 tot 30) werd uitgevoerd.
   • Dense Optical Flow (Farnebäck): Een dense per-pixel flow methode met behulp van polynomiale expansies en een piramide-structuur.
   • Strategieën: Zowel "periodic-reset" (Euleriaans, het opnieuw verankeren van punten naar initiële posities) als "no-reset" (Lagrangiaanse) tracking-strategieën werden geëvalueerd.

3. Metrieken:

   • Schijf Benchmark: Gemiddelde relatieve fout (%).
   • Fantoom Validatie: Genormaliseerde Root Mean Square Error (NRMSE) en Pearson correlatie ($r$) tegen de exacte grondwaarheid voor zowel lineaire als hoeksnelheden.
   • Ruisrobuustheid: Prestaties geëvalueerd over verschillende ruisniveaus ($\sigma \in [0, 0.12]$).

Belangrijkste Bijdragen

• Regime-Afhankelijke Optimalisatie: De studie toont aan dat de optimale $winSize$ voor PLK niet universeel is maar regime-afhankelijk. Een kleiner venster ($\approx 11\text{--}12$) is optimaal voor grote, discrete enkelvoudige stapsgewijze verplaatsingen (schijf benchmark), terwijl een groter venster ($\approx 20$) nodig is voor kleine, continue pulsatiele beweging (AAA fantoom) om ruis te middelen.
• Synthetische Grondwaarheid voor Hoeksnelheid: De introductie van een synthetische AAA fantoom met exacte closed-form grondwaarheid maakt de eerste kwantitatieve validatie van de reconstructie van hoeksnelheid in deze context mogelijk.
• Vergelijking Sparse vs. Dense: Een systematische vergelijking laat zien dat PLK en Farnebäck complementair zijn in plaats van vervangbaar. PLK is superieur voor grote geïsoleerde verplaatsingen en omgevingen met hoge ruis, terwijl Farnebäck volledige frame-dekking biedt maar faalt bij grote single-feature verplaatsingen tegen een statische achtergrond en gevoeliger is voor ruis in kleine bewegingsregimes.
• Validatie van Fysiologische Consistentie: De gereconstrueerde piek hoeksnelheid uit de fantoom ($\approx 0.20\text{--}0.23$ rad/s) komt overeen met het bereik dat in de oorspronkelijke real-patient studie werd waargenomen ($+0.13$ tot $+0.22$ rad/s) en de bestaande 4D-flow literatuur.

Resultaten

• Schijf Benchmark: PLK bereikte een minimale fout van $<0.5\%$ bij $winSize \approx 11\text{--}12$, waarmee de bevindingen van de oorspronkelijke studie werden gereproduceerd. Farnebäck faalde uniform ($\approx 83\text{--}85\%$ fout) door schendingen van de lokale gladheidsprior wanneer een object meer dan zijn diameter per frame verplaatst.
• Lineaire Snelheid (Fantoom): Beide methoden legden de systolische piek correct vast. PLK presteerde beter dan Farnebäck met een lagere NRMSE ($28.4\%$ vs. $32.4\%$) en een hogere correlatie ($r=0.93$ vs. $0.91$).
• Hoeksnelheid (Fantoom): PLK behaalde een NRMSE van $15.4\%$ ($r=0.91$), wat significant beter is dan Farnebäck ($25.0\%$ NRMSE, $r=0.71$). De gereconstrueerde piekvelociteiten kwamen overeen met het fysiologische bereik van de oorspronkelijke studie.
• Reset Strategie: Op de synthetische fantoom presteerden de "reset" en "no-reset" strategieën vrijwel identiek. De auteurs merken op dat dit een niet-bevestigende bevinding is met betrekking tot de grootte van de drift, aangezien de bewegingsschaal van de fantoom de geaccumuleerde drift die in echte data wordt gezien, niet reproduceerde. De reset-strategie wordt echter nog steeds aanbevolen om volledige lumen-dekking te garanderen en het ontstaan van gaten (voids) te voorkomen.
• Ruisrobuustheid: PLK vertoonde een superieure robuustheid tegen ruis vergeleken met Farnebäck. Bij hoge ruisniveaus ($\sigma=0.12$) was de lineaire snelheid NRMSE van PLK $94\%$ versus $162\%$ voor Farnebäck.

Significantie en Claims
Het artikel claimt een concreet, actiegericht pad te bieden voor optical flow-gebaseerde snelheidsreconstructie in AAA-analyse. Door de bevindingen van de oorspronkelijke studie te valideren tegen een synthetische grondwaarheid en ze te vergelijken met dense flow-methoden, stellen de auteurs vast dat:

1. PLK de robuuste keuze is voor kwantitatieve tracking in cine-MRI, met name waar grote verplaatsingen of ruis aanwezig zijn, mits de venstergrootte wordt afgestemd op het specifieke bewegingsregime.
2. Dense flow (Farnebäck) geschikt is voor kwalitatieve volledige veldvisualisatie, maar geen directe vervanging is voor sparse tracking in deze specifieke toepassing.
3. De periodic-reset strategie geldig blijft voor het behoud van dekking, zelfs als de exacte magnitude van de drift niet volledig werd gereproduceerd in het synthetische model.
4. Toekomstig werk is verankerd aan een specifieke, open-access multi-patiënt dataset (Jamshidian et al., 2025), waardoor de beperking van de "single-patient" uit de oorspronkelijke studie wordt omgezet in een gedefinieerde volgende stap. De auteurs stellen expliciet dat multi-patiënt validatie, 4D flow integratie en learning-based optical flow openstaan, maar nu gegrond zijn in geïdentificeerde bronnen.