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

1. Probleemstelling

Accurate 3D volumetrische snelheidsschatting is cruciaal voor diagnostische en therapeutische echografie-toepassingen (bijv. weefselelasticiteit, functionele beeldvorming en bloedstroomdynamica). Traditionele volledig gepopuleerde matrixarrays ondervinden echter aanzienlijke hardware- en computationele barrières (kosten, omvang, vermogen) die hun klinische adoptie belemmeren.

Row-Column Addressed (RCA) arrays bieden een oplossing door de hardwarecomplexiteit te verminderen. Ze introduceren echter specifieke beeldvormingsuitdagingen:

• Slecht contrast en resolutie: RCA-arrays lijden aan sterke zijlobben en verlengde Point Spread Functions (PSF) als gevolg van de gescheiden transmissie van rijen en kolommen.
• Beperkingen in snelheidsschatting: Hoewel methoden zoals Orthogonal Plane Wave (OPW) compounding de beeldkwaliteit verbeteren, resulteren ze vaak in onnauwkeurige bewegingsschattingen door ruimtelijke middeling veroorzaakt door zijlobben.
• Aliasing: Conventionele faseverschuivingsautocorrelatoren (Kasai-schatter) toegepast op gecompounde RCA-gegevens hebben een beperkte Nyquist-snelheid, wat leidt tot aliasing bij hoge stroomsnelheden (bijv. halsslagaders).
• Gat in de literatuur: Eerdere verbeteringen van RCA-beeldvorming (bijv. XDoppler, RC-FMAS) richtten zich voornamelijk op Power Doppler (contrast) en resolutie, waardoor snelheidsschatting grotendeels onopgelost bleef.

2. Methodologie

De auteurs stellen een nieuwe XDoppler-snelheidsschatter voor die het cross-correlatieconcept van de XDoppler Power Doppler-methode uitbreidt naar axiale snelheidsschatting.

A. Acquisitiestrategie

• Hardware: Een 128+128 element RCA-probe (6 MHz centrale frequentie) aangesloten op een Verasonics Vantage 256 scanner.
• Sequentie: Het systeem acquireert twee sub-volumes achtereenvolgens:
  1. RC Sub-volume: Plangolven uitgezonden door rijen, ontvangen door kolommen.
  2. CR Sub-volume: Plangolven uitgezonden door kolommen, ontvangen door rijen.
• Temporale vertraging: Er is een specifieke tijdvertraging ($\Delta T$) tussen de acquisitie van de RC- en CR-sub-volumes.

B. Wiskundige Afleiding

De kerninnovatie ligt in het afleiden van de snelheid uit de fase van de cross-correlatie tussen de orthogonale sub-volumes ($s_{RC}$ en $s_{CR}$).

• Het Bias-probleem: Een eenvoudige cross-correlatie $s_{RC}(i)s_{CR}^*(i+\Delta T)$ levert een bevooroordeelde fase-schatting op omdat de statische fase-termen ($\phi_{RC}$ en $\phi_{CR}$) specifiek voor elke apertuur niet wegvallen.
• De Oplossing: De auteurs stellen een gecombineerde schatter voor die twee complementaire cross-correlaties optelt om de Doppler-faseverschuiving ($\Delta \phi_{Doppler}$) te isoleren:
  $$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]$$
  Door de argumenten van deze twee termen op te tellen, vallen de statische fase-biasen weg, waardoor alleen de Doppler-verschuiving overblijft.

C. Theoretische Voordelen

• Nyquist-snelheid: Door de totale observatietijd $T$ op te splitsen in twee kortere intervallen ($\Delta T_1 = \Delta T_2 = T/2$), wordt de effectieve verplaatsingstijd gehalveerd. Dit verdubbelt het ondubbelzinnige snelheidsbereik (Nyquist-limiet) in vergelijking met de conventionele Kasai-schatter toegepast op OPW-compounding:
  $$v_{N,X} = 2 \times v_{N,comp}$$
• Ruimtelijke resolutie: De cross-correlatie-operatie onderdrukt inherent zijlobben, wat leidt tot een scherpere PSF en beter contrast in vergelijking met coherente OPW-sommatie.

D. Experimentele Validatie

• In Vitro: Stroomfantomen (4 mm en 8 mm diameter) met bloed-imitatievloeistof werden gebruikt om nauwkeurigheid, convergentietijd en aliasing-robuustheid te testen tegen een Poiseuille-stroommodel.
• In Vivo: Acquisities werden uitgevoerd op de halsslagader en de vena jugularis van een gezonde vrijwilliger om prestaties te beoordelen bij pulserende stroming en aliasing-scenario's.
• Metingen: Root Mean Squared Error (RMSE), stroomsnelheidsbias, standaardafwijking en convergentietijd ($t_{90\%}$).

3. Belangrijkste Bijdragen

1. Nieuwe Snelheidsschatter: De eerste afleiding van een op fase gebaseerde axiale snelheidsschatter specifiek voor RCA-arrays met behulp van het XDoppler cross-correlatiekader.
2. Bias-correctie: Een wiskundig bewijs en implementatie van een methode om statische fase-biasen inherent aan orthogonale apertuur cross-correlaties te annuleren.
3. Verdubbelde Nyquist-limiet: Aantonen dat de methode theoretisch en experimenteel een Nyquist-snelheid bereikt die twee keer zo hoog is als de standaard Kasai-schatter op OPW-gegevens.
4. Versterkt contrast en resolutie: De methode maakt gebruik van de zijlob-onderdrukking van XDoppler om nauwkeurigere snelheidsprofielen te bieden, met name nabij vaatwanden waar conventionele methoden falen door ruimtelijke middeling.

4. Resultaten

In Vitro Prestaties

• Nauwkeurigheid: De XDoppler-schatter kwam nauwkeurig overeen met het theoretische Poiseuille-parabolische snelheidsprofiel. Daarentegen onderschatte de conventionele Kasai/OPW-schatter centrale snelheden met ongeveer 10% en overschatte snelheden nabij de wanden.
• Stroomsnelheidsschatting: XDoppler toonde een gemiddelde bias van -0,17 mL/min en een standaardafwijking van 3,43 mL/min, wat aanzienlijk beter presteerde dan de Kasai-schatter (Bias: -9,33 mL/min; Std Dev: 6,19 mL/min).
• Convergentie: De schatter stabiliseerde binnen 0,05s tot 0,2s en bereikte snel 90% van de uiteindelijke waarde.
• Aliasing-resistentie: In tests waarbij stroomsnelheden de Nyquist-limiet van de Kasai-schatter (7,7 cm/s) overschreden maar onder de XDoppler-limiet (15,4 cm/s) bleven, vertoonde de Kasai-schatter ernstige aliasing, terwijl XDoppler correcte metingen leverde.

In Vivo Prestaties (Halsslagader)

• Dynamische stroming: De methode slaagde erin de dynamiek van pulserende stroming over meerdere hartcycli (ongeveer 44 bpm) succesvol te volgen.
• Aliasing-reductie: In de vena jugularis vertoonde de conventionele schatter aliasing tijdens diastole vanwege de vaatgeometrie, terwijl XDoppler dit niet deed.
• Snelheidsmagnitude: Tijdens systole mat de XDoppler-schatter hogere pieksnelheden (1,54 cm/s) in vergelijking met OPW (1,32 cm/s), toegeschreven aan de verminderde aliasing en betere behoud van hoog-snelheidssignalen.

5. Betekenis en Conclusie

Deze studie overbrugt succesvol de kloof tussen hoogwaardige 3D-beeldvorming en accurate snelheidskwantificering voor RCA-arrays.

• Klinische impact: De methode maakt betrouwbaardere 3D kleur-Doppler-beeldvorming mogelijk in complexe stromingsomgevingen (bijv. stenose, klepinsufficiëntie) waar hoge snelheden en aliasing veel voorkomen.
• Technische vooruitgang: Het bewijst dat de cross-correlatie van orthogonale aperturen niet alleen een hulpmiddel is voor Power Doppler-contrast, maar een robuust mechanisme voor snelheidsschatting, wat een 2x verbetering in Nyquist-snelheid biedt zonder in te leveren op ruimtelijke resolutie.
• Toekomstperspectieven: De auteurs suggereren het integreren van deze schatter met Vector Flow Imaging (VFI) om volledige 3D-snelheidsvectoren op te lossen en het toepassen van de-aliasing-technieken voor toepassingen met nog hogere snelheden (bijv. cardiale beeldvorming).

Kortom, de XDoppler-snelheidsschatter vertegenwoordigt een belangrijke stap voorwaarts om 4D (volumetrisch + tijd) echografie-beeldvorming klinisch haalbaar te maken voor kwantitatieve hemodynamische analyse.