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

技術要約：均一な密度分布と逐次ビンニングのための改良型3D放射状フィロタキシス・トラジェクトリ

問題提起
3次元放射状トラジェクトリ、特にスパイラル・フィロタキシスは、その運動ロバスト性と自己ナビゲーション能力により、ダイナミックおよびモーション分解MRIにおいて広く用いられている。しかし、標準的なスパイラル・フィロタキシスは、主に2つの制限事項を抱えている：

1. 非一様なサンプリング密度： 標準的なトラジェクトリは、サンプリング密度（SDSD）の均一な球面分布を実現できない。赤道領域付近に読み出し方向が集中する傾向があり、これによりk空間サンプリングの異方性、実効的な視野および空間分解能の方向依存性の変化、そして潜在的なエイリアシングアーチファクトが生じる。
2. 逐次ビンニングによるアーチファクト： グローバルなk空間被覆は効率的であるが、多くのアプリケーション（例：ダイナミックイメージング、機能的MRI、モーション分解再構成）では、時間的に連続する読み出しのサブセットから画像を再構成すること（逐次ビンニング）が必要となる。標準的なフィロタキシスでは、連続するインターリーブは極角において系統的なクラスタリングを示す。これらの連続するサブセットをビンニングすると、k空間内に円形のギャップが生じ、それが再構成画像におけるリンギングアーチファクトや時間的な不整合として現れる。

手法
著者らは、3T MRIスキャナーを用い、2つのグラジエントエコーシーケンス（T1強調GREおよびLIBRE）を用いて10名の健康なボランティアに対してin vivoで評価された、2つの手法的な進展を提案している。

1. Uniform Phyllotaxis (UPhy):

   • 目的： 読み出し方向の均一な密度分布（均一なSDSD）を実現すること。
   • メカニズム： UPhyは、標準的なフィロタキシスの極角の定式化を再定義する。標準モデルでは $\theta_n = \frac{\pi}{2}\sqrt{\frac{n}{N}}$ を使用するが、UPhyは $\theta_n = \arccos(1 - \frac{n}{N})$ を利用する。
   • 効果： この修正により、元の方位角（ゴールデンアングル）の順序を保持しつつ、球面上への点の均一な分布を保証し、標準モデルで見られる赤道付近へのクラスタリングを排除する。

2. Flexible Phyllotaxis (FlexiPhy):

   • 目的： 任意の遡及的な逐次ビンニングに対するロバスト性を向上させること。
   • メカニズム： FlexiPhyはUPhyを基盤とし、インターリーブのランダムな置換（$\sigma$）を導入する。これにより、極角の順序と方位角のゴールデンアングル進行をデカップリング（分離）する。具体的には、$k$番目のインターリーブにおける$j$番目の点の極角は、連続的なインデックス$k$ではなく、ランダムに置換されたインデックス$\sigma(k)$によって決定される。
   • 効果： このランダム化により、極角の系統的なクラスタリングが打破される。その結果、連続する時間的なビン内のサンプリングパターンが良好に維持され、同心円状のk空間ギャップの形成と、それに伴うリンギングアーチファクトを防止する。

主な貢献

• 解析的および経験的な一様性： 本論文は、UPhyが解析的に証明された一様なSDSDを提供し、標準的なフィロタキシスの密度バイアスを補正することを実証している。
• ビンニングのためのデカップリング： FlexiPhyの導入は、ランダム化を通じて極角と方位角の順序をデカップリングする新しい戦略を提示しており、これは特に、時間的に連続するサブセット内でのサンプリング特性の劣化に対処するものである。
• In Vivo 検証： これらのトラジェクトリは、静止脳イメージングおよびダイナミック眼球イメージングのシナリオにおいて、2つの異なるシーケンス（標準GREおよび脂肪抑制LIBRE）を用いてテストされた。

結果

• 静止イメージング： 逐次的な時間ビン再構成の目視検査により、標準的なフィロタキシスとUPhyの両方が、目に見えるリンギングアーチファクトと時間的不整合を示すことが明らかになった。FlexiPhyはこれらのアーチファクトを大幅に軽減し、より時間的に安定した画像を作成した。
• ダイナミックイメージング： ダイナミック眼球イメージングにおいて、標準およびUPhyのトラジェクトリは、眼窩領域へと伝播する構造的なリンギングアーチファクトを発生させ、解剖学的詳細を不明瞭にした。FlexiPhyはこれらのアーチファクトを大幅に減少させ、眼球および周囲の構造の可視化を改善した。
• 定量的指標：
  • SSIMおよびL2誤差： GREおよびLIBREの両方のシーケンスにおいて、FlexiPhyは標準的なフィロタキシスおよびUPhyと比較して、有意に高い構造的類似性指数（SSIM）値と低い相対L2誤差を示した（ボンフェローニ補正済み $p < 0.05$）。
  • UPhy vs. 標準： UPhyは理論的な角度の一様性を改善したが、LIBREシーケンスにおける軽微なL2誤差の改善を除き、標準的なトラジェクトリと比較して逐次ビンニングにおけるリンギングアーチファクトを有意に減少させることはできなかった。これは、グローバルなSDSDの改善だけでは、逐次ビンニングのアーチファクトを解決するには不十分であることを裏付けている。

意義と主張
本論文は、FlexiPhyが、トラジェクトリに起因するリンギングアーチファクトと時間的不整合を低減することで、より信頼性の高い逐次ビンニング再構成を可能にすると結論づけている。そのランダム化された構成は特定のビンニング戦略に縛られず、ダイナミックおよびモーション分解MRIで使用される遡及的なビンニング手法と広く互換性がある。

著者らは、UPhyがグローバルな密度分布を補正する一方で、FlexiPhyにおける極角と方位角の順序のデカップリングこそが、時間的なサブセット内でのサンプリング品質を維持するために極めて重要であることを強調している。本研究は、FlexiPhyがスパイラル・フィロタキシスの実用的な利点（効率的なk空間被覆、滑らかなグラジエント遷移）を保持しつつ、時間分解能や遡及的なソートを必要とするアプリケーションに対して優れたロバスト性を提供することを示唆している。また、本研究の限界として、モーション推定の精度への影響を調査することや、これらのトラジェクトリとポール・トゥ・ポール（極間）取得技術との組み合わせを検討する必要があることを挙げている。