Improved 3D Radial Phyllotaxis Trajectories for Uniform Density Distribution of Readout Directions and Sequential Binning

Technische Samenvatting: Verbeterde 3D Radiale Fyllotaxis Trajecten voor Uniforme Dichtheidsdistributie en Sequentiële Binning

Probleemstelling
Driedimensionale radiale trajecten, specifiek spiraalvormige fyllotaxis, worden breed gebruikt in dynamische en bewegings-opgeloste MRI vanwege hun robuustheid tegen beweging en zelf-navigatiecapaciteiten. De standaard spiraalvormige fyllotaxis-trajecten lijden echter onder twee primaire beperkingen:

1. Niet-uniforme bemonsteringsdichtheid: Het standaard traject bereikt geen uniforme sferische verdeling van de bemonsteringsdichtheid (SDSD). Het heeft de neiging om uitleesrichtingen te concentreren nabij de equatoriale regio, wat leidt tot anisotrope k-ruimte bemonstering, richtingafhankelijke variaties in het effectieve gezichtsveld en de ruimtelijke resolutie, en potentiële aliasing-artefacten.
2. Artefacten door sequentiële binning: Hoewel de globale k-ruimte dekking efficiënt is, vereisen veel toepassingen (bijv. dynamische beeldvorming, functionele MRI, bewegings-opgeloste reconstructie) het reconstrueren van afbeeldingen uit tijdelijk opeenvolgende subsets van uitlezingen (sequentiële binning). In standaard fyllotaxis vertonen opeenvolgende interleaves een systematische clustering in polaire hoeken. Wanneer deze opeenvolgende subsets worden gebinned, creëren ze circulaire gaten in de k-ruimte, wat zich manifesteert als ringartefacten en temporele inconsistenties in de gereconstrueerde beelden.

Methodologie
De auteurs stellen twee methodologische verbeteringen voor om deze beperkingen aan te pakken, geëvalueerd in vivo op 10 gezonde vrijwilligers met behulp van een 3T MRI-scanner met twee gradiënt-echo-sequenties (T1-gewogen GRE en LIBRE).

1. Uniforme Fyllotaxis (UPhy):

   • Doel: Om een uniforme dichtheidsverdeling van uitleesrichtingen te bereiken (uniforme SDSD).
   • Mechanisme: UPhy herdefinieert de polaire hoekformulering van de standaard fyllotaxis. Terwijl het standaard model $\theta_n = \frac{\pi}{2}\sqrt{\frac{n}{N}}$ gebruikt, maakt UPhy gebruik van $\theta_n = \arccos(1 - \frac{n}{N})$.
   • Effect: Deze modificatie zorgt voor een uniforme verdeling van punten over de sfeer, waardoor de equatoriale clustering die in het standaard model wordt waargenomen, wordt geëlimineerd, terwijl de oorspronkelijke azimutale (gouden hoek) ordening behouden blijft.

2. Flexibele Fyllotaxis (FlexiPhy):

   • Doel: Om de robuustheid tegen willekeurige retrospectieve sequentiële binning te verbeteren.
   • Mechanisme: FlexiPhief bouwt voort op UPhy maar introduceert een gerandomiseerde permutatie ($\sigma$) van de interleaves. Dit ontkoppelt de ordening van de polaire hoek van de azimutale gouden hoek-progressie. Specifiek wordt de polaire hoek voor het $j$-de punt in de $k$-de interleaf bepaald door een gerandomiseerde gepermuteerde index $\sigma(k)$, in plaats van de sequentiële index $k$.
   • Effect: Deze randomisatie verbreekt de systematische clustering van polaire hoeken binnen opeenvolgende temporele bins. Bijgevolg behouden opeenvolgende subsets van uitlezingen een goed verdeeld bemonsteringspatroon, wat de vorming van concentrische k-ruimte gaten en de bijbehorende ringartefacten voorkomt.

Belangrijkste Bijdragen

• Analytische en Empirische Uniformiteit: De paper demonstreert dat UPhy een analytisch bewezen uniforme SDSD biedt, waarmee de dichtheidsbias van de standaard fyllotaxis wordt gecorrigeerd.
• Ontkoppeling voor Binning: De introductie van FlexiPhy biedt een nieuwe strategie om de polaire en azimutale ordening te ontkoppelen via randomisatie, specifiek gericht op het verbeteren van de bemonsteringskwaliteit binnen temporele subsets.
• In Vivo Validatie: De trajecten werden getest met twee verschillende sequenties (standaard GRE en vetonderdrukkende LIBRE) tijdens statische hersenbeeldvorming en dynamische oogbeeldvorming scenario's.

Resultaten

• Statische Beeldvorming: Visuele inspectie van sequentiële temporele-bin reconstructies toonde aan dat zowel de standaard fyllotaxis als UPhy zichtbare ringartefacten en temporele inconsistenties vertoonden. FlexiPhy verminderde deze artefacten aanzienlijk, wat resulteerde in meer temporeel stabiele beelden.
• Dynamische Beeldvorming: Bij dynamische oogbeeldvorming produceerden de standaard en UPhy trajecten gestructureerde ringartefacten die doorwerkten in de orbitale regio en anatomische details maskeerden. FlexiPhy verminderde deze artefacten aanzienlijk, wat de visualisatie van het oog en de omliggende structuren verbeterde.
• Kwantitatieve Metrieken:
  • SSIM en L2-fout: FlexiPhy vertoonde significant hogere Structural Similarity Index Measure (SSIM) waarden en lagere relatieve L2-fouten vergeleken met zowel de standaard fyllotaxis als UPhy voor zowel GRE als LIBRE sequenties (Bonferroni-gecorrigeerde $p < 0,05$).
  • UPhy versus Standaard: Hoewel UPhy de theoretische hoekuniformiteit verbeterde, verminderde het de ringartefacten bij sequentiële binning niet significant vergeleken met de standaard traject, met uitzondering van een kleine verbetering in de L2-fout voor de LIBRE-sequentie. Dit bevestigt dat verbetering van de globale SDSD alleen onvoldoende is om problemen met sequentiële binning op te lossen.

Betekenis en Claims
De paper concludeert dat FlexiPhy meer betrouwbare sequentiële binning reconstructies mogelijk maakt door traject-geïnduceerde ringartefacten en temporele inconsistenties te verminderen. De gerandomiseerde constructie is niet gebonden aan een specifieke binningstrategie, waardoor het breed compatibel is met retrospectieve binning-benaderingen die worden gebruikt in dynamische en bewegings-opgeloste MRI.

De auteurs benadrukken dat hoewel UPhy de globale dichtheidsverdeling corrigeert, het de ontkoppeling van de polaire en azimutale ordening in FlexiPhy is die cruciaal is voor het behoud van de bemonsteringskwaliteit binnen temporele subsets. Het werk suggereert dat FlexiPhy de praktische voordelen van spiraalvormige fyllotaxis (efficiënte k-ruimte dekking, vloeiende gradiëntovergangen) behoudt terwijl het een superieure robuustheid biedt voor toepassingen die een temporele resolutie of retrospectieve sortering vereisen. De studie erkent beperkingen en merkt op dat toekomstig werk de impact op downstream taken zoals bewegingsschattingsnauwkeurigheid moet onderzoeken en de combinatie van deze trajecten met pole-to-pole acquisitietechnieken moet verkennen.