Revisiting XDoppler estimator for high spatiotemporal resolution volumetric axial velocity measurement using row-column arrays

1. Enunciazione del Problema

Una stima accurata della velocità volumetrica 3D è fondamentale per le applicazioni diagnostiche e terapeutiche degli ultrasuoni (ad es. elasticità dei tessuti, imaging funzionale e dinamica del flusso ematico). Tuttavia, le tradizionali matrici completamente popolate affrontano significativi ostacoli hardware e computazionali (costo, ingombro, potenza) che ne ostacolano l'adozione clinica.

Le matrici ad indirizzamento riga-colonna (RCA) offrono una soluzione riducendo la complessità hardware. Tuttavia, introducono sfide specifiche nell'imaging:

• Basso Contrasto e Risoluzione: Le matrici RCA soffrono di lobi laterali intensi e di Funzioni di Risposta al Punto (PSF) allungate a causa della trasmissione separata di righe e colonne.
• Limitazioni nella Stima della Velocità: Sebbene metodi come la composizione di Onde Piane Ortogonali (OPW) migliorino la qualità dell'immagine, spesso comportano una stima del moto inaccurata a causa della media spaziale causata dai lobi laterali.
• Aliasing: Gli autocorrelatori a spostamento di fase convenzionali (stimatore di Kasai) applicati ai dati RCA composti hanno una velocità di Nyquist limitata, portando all'aliasing nei flussi ad alta velocità (ad es. arterie carotidi).
• Lacuna nella Letteratura: I precedenti miglioramenti all'imaging RCA (ad es. XDoppler, RC-FMAS) si sono concentrati principalmente sul Power Doppler (contrasto) e sulla risoluzione, lasciando la stima della velocità largamente irrisolta.

2. Metodologia

Gli autori propongono un nuovo stimatore di velocità XDoppler che estende il concetto di correlazione incrociata del metodo Power Doppler XDoppler alla stima della velocità assiale.

A. Strategia di Acquisizione

• Hardware: Una sonda RCA da 128+128 elementi (frequenza centrale di 6 MHz) collegata a uno scanner Verasonics Vantage 256.
• Sequenza: Il sistema acquisisce due sottovolumi in sequenza:
  1. Sottovolume RC: Onde piane trasmesse dalle Righe, ricevute dalle Colonne.
  2. Sottovolume CR: Onde piane trasmesse dalle Colonne, ricevute dalle Righe.
• Ritardo Temporale: Esiste un ritardo temporale specifico ($\Delta T$) tra l'acquisizione dei sottovolumi RC e CR.

B. Derivazione Matematica

L'innovazione centrale risiede nella derivazione della velocità dalla fase della correlazione incrociata tra i sottovolumi ortogonali ($s_{RC}$ e $s_{CR}$).

• Il Problema del Bias: Una semplice correlazione incrociata $s_{RC}(i)s_{CR}^*(i+\Delta T)$ produce una stima di fase distorta perché i termini di fase statici ($\phi_{RC}$ e $\phi_{CR}$) specifici di ciascuna apertura non si annullano.
• La Soluzione: Gli autori propongono un stimatore combinato che somma due correlazioni incrociate complementari per isolare lo spostamento di fase Doppler ($\Delta \phi_{Doppler}$):
  $$v_z = -\frac{c}{4\pi f_0 T} \left[ \arg \left( \sum s_{RC}(i)s_{CR}^*(i+0.5) \right) + \arg \left( \sum s_{CR}(i+0.5)s_{RC}^*(i+1) \right) \right]$$
  Sommando gli argomenti di questi due termini, i bias di fase statici si annullano, lasciando solo lo spostamento Doppler.

C. Vantaggi Teorici

• Velocità di Nyquist: Dividendo il tempo di osservazione totale $T$ in due intervalli più brevi ($\Delta T_1 = \Delta T_2 = T/2$), il tempo di spostamento effettivo è dimezzato. Questo raddoppia l'intervallo di velocità non ambiguo (limite di Nyquist) rispetto allo stimatore di Kasai convenzionale applicato alla composizione OPW:
  $$v_{N,X} = 2 \times v_{N,comp}$$
• Risoluzione Spaziale: L'operazione di correlazione incrociata sopprime intrinsecamente i lobi laterali, portando a una PSF più nitida e a un migliore contrasto rispetto alla somma coerente OPW.

D. Validazione Sperimentale

• In Vitro: Fantocci di flusso (diametri di 4 mm e 8 mm) con fluido mimante il sangue sono stati utilizzati per testare accuratezza, tempo di convergenza e robustezza all'aliasing rispetto a un modello di flusso di Poiseuille.
• In Vivo: Le acquisizioni sono state eseguite sull'arteria carotide e sulla vena giugulare di un volontario sano per valutare le prestazioni in flussi pulsati e scenari di aliasing.
• Metriche: Errore Quadratico Medio (RMSE), bias della portata, deviazione standard e tempo di convergenza ($t_{90\%}$).

3. Contributi Chiave

1. Nuovo Stimatore di Velocità: La prima derivazione di uno stimatore di velocità assiale basato sulla fase specificamente per matrici RCA che utilizza il framework di correlazione incrociata XDoppler.
2. Correzione del Bias: Una dimostrazione matematica e un'implementazione di un metodo per annullare i bias di fase statici intrinseci alle correlazioni incrociate di aperture ortogonali.
3. Limite di Nyquist Raddoppiato: Dimostrazione che il metodo raggiunge teoricamente e sperimentalmente una velocità di Nyquist doppia rispetto allo stimatore di Kasai standard sui dati OPW.
4. Contrasto e Risoluzione Migliorati: Il metodo sfrutta la soppressione dei lobi laterali di XDoppler per fornire profili di velocità più accurati, in particolare vicino alle pareti vasali dove i metodi convenzionali falliscono a causa della media spaziale.

4. Risultati

Prestazioni In Vitro

• Accuratezza: Lo stimatore XDoppler ha corrisposto strettamente al profilo di velocità parabolico teorico di Poiseuille. Al contrario, lo stimatore convenzionale Kasai/OPW ha sottostimato le velocità centrali di circa il 10% e sovrastimato le velocità vicino alle pareti.
• Stima della Portata: XDoppler ha mostrato un bias medio di -0.17 mL/min e una deviazione standard di 3.43 mL/min, superando significativamente lo stimatore di Kasai (Bias: -9.33 mL/min; Deviazione Standard: 6.19 mL/min).
• Convergenza: Lo stimatore si è stabilizzato entro 0.05s e 0.2s, raggiungendo il 90% del valore finale rapidamente.
• Resistenza all'Aliasing: Nei test in cui le velocità di flusso superavano il limite di Nyquist dello stimatore di Kasai (7.7 cm/s) ma rimanevano al di sotto del limite XDoppler (15.4 cm/s), lo stimatore di Kasai ha mostrato un grave aliasing, mentre XDoppler ha fornito misurazioni corrette.

Prestazioni In Vivo (Arteria Carotide)

• Flusso Dinamico: Il metodo ha tracciato con successo la dinamica del flusso pulsato su più cicli cardiaci (circa 44 bpm).
• Riduzione dell'Aliasing: Nella vena giugulare, lo stimatore convenzionale ha mostrato aliasing durante la diastole a causa della geometria vascolare, mentre XDoppler no.
• Entità della Velocità: Durante la sistole, lo stimatore XDoppler ha misurato velocità di picco più elevate (1.54 cm/s) rispetto a OPW (1.32 cm/s), attribuito alla ridotta aliasing e alla migliore conservazione dei segnali ad alta velocità.

5. Significato e Conclusione

Questo studio colma con successo il divario tra imaging 3D di alta qualità e quantificazione accurata della velocità per le matrici RCA.

• Impatto Clinico: Il metodo consente un imaging Doppler a colori 3D più affidabile in ambienti di flusso complessi (ad es. stenosi, rigurgito valvolare) dove sono comuni alte velocità e aliasing.
• Avanzamento Tecnico: Dimostra che la correlazione incrociata di aperture ortogonali non è solo uno strumento per il contrasto Power Doppler, ma un meccanismo robusto per la stima della velocità, offrendo un miglioramento di 2x nella velocità di Nyquist senza sacrificare la risoluzione spaziale.
• Prospettive Future: Gli autori suggeriscono di integrare questo stimatore con l'Imaging del Flusso Vettoriale (VFI) per risolvere i vettori di velocità 3D completi e applicare tecniche di de-aliasing per applicazioni con velocità ancora più elevate (ad es. imaging cardiaco).

In sintesi, lo stimatore di velocità XDoppler rappresenta un passo significativo in avanti nel rendere l'imaging ecografico 4D (volumetrico + tempo) clinicamente fattibile per l'analisi emodinamica quantitativa.