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

Sintesi Tecnica: Validazione Sintetica Estesa della Ricostruzione del Campo di Velocità basata su Optical Flow per la Cine-RM di AAA

Problematica
La ricostruzione dei campi di velocità da una cine-RM standard di un aneurisma dell'aorta addominale (AAA) offre una via per la stratificazione del rischio emodinamico oltre le semplici misurazioni del diametro. Sebbene la RM di Phase-Contrast e la 4D flow siano gli standard di riferimento, esse richiedono acquisizioni specializzate. La cine-RM standard, acquisita routinamente, può teoricamente essere analizzata utilizzando l'optical flow (OF) senza costi aggiuntivi. Tuttavia, uno studio precedente dello stesso autore [9], che utilizzava l'optical flow sparso Pyramidal Lucas–Kanade (PLK) su un singolo paziente con AAA, era limitato dalla mancanza di un confronto con il flusso denso, dall'assenza di una validazione della velocità angolare e dalla dipendenza da un dataset di un singolo paziente. Questa estensione affronta tali lacune validando rigorosamente il PLK rispetto ai metodi di optical flow denso e stabilendo una verità fondamentale (ground truth) sintetica sia per le velocità lineari che per quelle angolari.

Metodologia
Lo studio impiega un approccio di validazione multifasettico che combina la riproduzione di dati reali, benchmark indipendenti e un nuovo fantoccio sintetico:

1. Fonti di Dati:

   • Dati Reali di Pazienti: Riproduzione dei risultati dello studio originale [9] utilizzando cine-RM in asse corto e lungo di un paziente con AAA (Rigshospitalet, 2020).
   • Benchmark del Disco Sintetico: Un benchmark indipendente che coinvolge un disco in traslazione a velocità irregolare per testare la sensibilità alla dimensione della finestra ($winSize$) e al rumore.
   • Fantoccio AAA Sintetico: Un nuovo fantoccio basato sulla letteratura, costruito con una verità fondamentale a forma chiusa sia per le velocità lineari che angolari. I parametri del fantoccio (strain della parete circonferenziale, spaziatura dei voxel, diametro del lume e velocità angolare) sono derivati da un nuovo dataset 4D-CTA ad accesso aperto di 10 pazienti con AAA (Jamshidian et al., 2025). Il fantoccio simula una dilatazione pulsatile e una rotazione rigida in senso antiorario, con l'aggiunta di rumore Rician per mimare le immagini di magnitudo MRI.

2. Algoritmi e Configurazioni:

   • Optical Flow Sparso (PLK): Pyramidal Lucas–Kanade con una piramide Gaussiana a 4 livelli. È stata eseguita una scansione di sensibilità per $winSize$ (da 4 a 30).
   • Optical Flow Denso (Farnebäck): Un metodo di flusso denso per pixel utilizzando espansioni polinomiali e una struttura a piramide.
   • Strategie: Sono state valutate entrambe le strategie di tracciamento "periodic-reset" (Euleriana, ri-ancoraggio dei punti alle posizioni iniziali) e "no-reset" (Lagrangiana).

3. Metriche:

   • Benchmark del Disco: Errore relativo medio (%).
   • Validazione del Fantoccio: Errore quadratico medio normalizzato (NRMSE) e correlazione di Pearson ($r$) rispetto alla verità fondamentale esatta per entrambe le velocità lineari e angolari.
   • Robustezza al Rumore: Prestazioni valutate attraverso diversi livelli di rumore ($\sigma \in [0, 0.12]$).

Contributi Chiave

• Ottimizzazione Dipendente dal Regime: Lo studio dimostra che la $winSize$ ottimale per il PLK non è universale ma dipende dal regime. Una finestra più piccola ($\approx 11\text{--}12$) è ottimale per grandi spostamenti discreti a singolo step (benchmark del disco), mentre una finestra più grande ($\approx 20$) è necessaria per il movimento pulsatile continuo (fantoccio AAA) per mediare il rumore.
• Verità Fondamentale Angolare Sintetica: L'introduzione di un fantoccio AAA sintetico con verità fondamentale esatta a forma chiusa consente la prima validazione quantitativa del recupero della velocità angolare in questo contesto.
• Confronto tra Sparso e Denso: Un confronto sistematico rivela che PLK e Farnebäck sono complementari piuttosto che intercambiabili. Il PLK è superiore per grandi spostamenti isolati e ambienti ad alto rumore, mentre il Farnebäck fornisce una copertura completa del frame ma fallisce su grandi spostamenti di singole caratteristiche contro sfondi statici ed è più sensibile al rumore nei regimi di piccolo movimento.
• Validazione della Coerenza Fisiologica: La velocità angolare di picco recuperata dal fantoccio ($\approx 0.20\text{--}0.23$ rad/s) si allinea con l'intervallo osservato nello studio originale su paziente reale ($+0.13$ a $+0.22$ rad/s) e con la letteratura esistente sulla 4D-flow.

Risultati

• Benchmark del Disco: Il PLK ha raggiunto un errore minimo $<0.5\%$ con $winSize \approx 11\text{--}12$, riproducendo i risultati dello studio originale. Il Farnebäck è fallito uniformemente ($\approx 83\text{--}85\%$ di errore) a causa delle violazioni del prior di regolarità locale quando un oggetto si sposta più del proprio diametro per frame.
• Velocità Lineare (Fantoccio): Entrambi i metodi hanno catturato correttamente il picco sistolico. Il PLK ha superato il Farnebäck con un NRMSE inferiore ($28.4\%$ vs $32.4\%$) e una correlazione più alta ($r=0.93$ vs $0.91$).
• Velocità Angolare (Fantoccio): Il PLK ha ottenuto un NRMSE del $15.4\%$ ($r=0.91$), superando significativamente il Farnebäck ($25.0\%$ NRMSE, $r=0.71$). Le velocità di picco recuperate corrispondevano all'intervallo fisiologico dello studio originale.
• Strategia di Reset: Sul fantoccio sintetico, le strategie "reset" e "no-reset" si sono comportate in modo quasi identico. Gli autori osservano che questa è una scoperta non confirmativa riguardo l'entità del drift, poiché la scala di movimento del fantoccio non ha riprodotto l'accumulo di drift visto nei dati reali. Tuttavia, la strategia di reset è comunque raccomandata per garantire la copertura dell'intero lume e prevenire la formazione di vuoti.
• Robustezza al Rumore: Il PLK ha dimostrato una robustezza superiore al rumore rispetto al Farnebäck. Ad alti livelli di rumore ($\sigma=0.12$), l'NRMSE della velocità lineare del PLK era del $94\%$ rispetto al $162\%$ del Farnebäck.

Significatività e Rivendicazioni
Il documento sostiene di fornire un percorso concreto e azionabile per la ricostruzione della velocità basata su optical flow nell'analisi dell'AAA. Validando i risultati dello studio originale rispetto a una verità fondamentale sintetica e confrontandoli con i metodi di flusso denso, gli autori stabiliscono che:

1. Il PLK è la scelta robusta per il tracciamento quantitativo nella cine-RM, particolarmente dove sono presenti grandi spostamenti o rumore, a patto che la dimensione della finestra sia calibrata sul regime di movimento specifico.
2. Il flusso denso (Farnebäck) è adatto per la visualizzazione qualitativa a campo completo, ma non è un sostituto diretto del tracciamento sparso in questa specifica applicazione.
3. La strategia di periodic-reset rimane valida per mantenere la copertura, anche se l'entità del drift non è stata pienamente riprodotta nel modello sintetico.
4. Il lavoro futuro è ancorato a un dataset specifico e ad accesso aperto di più pazienti (Jamshidian et al., 2025), trasformando il limite del "singolo paziente" dello studio originale in un passo successivo definito. Gli autori dichiarano esplicitamente che la validazione su più pazienti, l'integrazione con la 4D-flow e l'optical flow basato su apprendimento automatico rimangono aperti, ma sono ora radicati in risorse identificate.