[論文レビュー] Harmonic analysis of spherical sampling in diffusion MRI

Alessandro Daducci, Jason D. McEwen+3 more|arXiv (Cornell University)|Jun 1, 2011

Advanced Neuroimaging Techniques and Applications被引用数 10

ひとこと要約

本稿では、球面調和解析を用いた理論的枠組みを提示し、拡散MRIにおける球面サンプリングを分析する。特に、q空間データにおける aliasing アーチファクトに及ぼす角度分解能の影響に焦点を当てる。球面調和（SH）変換の厳密な解析と、方向分布関数（ODF）の計算におけるその役割を通じて、サンプリング不足が再構成誤差を引き起こすメカニズムを定量化し、マルチシェル測定の最適化に役立てる基盤を提供する。これにより、正確性を損なわず測定時間を短縮できる。

ABSTRACT

In the last decade diffusion MRI has become a powerful tool to non-invasively study white-matter integrity in the brain. Recently many research groups have focused their attention on multi-shell spherical acquisitions with the aim of effectively mapping the diffusion signal with a lower number of q-space samples, hence enabling a crucial reduction of acquisition time. One of the quantities commonly studied in this context is the so-called orientation distribution function (ODF). In this setting, the spherical harmonic (SH) transform has gained a great deal of popularity thanks to its ability to perform convolution operations efficiently and accurately, such as the Funk-Radon transform notably required for ODF computation from q-space data. However, if the q-space signal is described with an unsuitable angular resolution at any b-value probed, aliasing (or interpolation) artifacts are unavoidably created. So far this aspect has been tackled empirically and, to our knowledge, no study has addressed this problem in a quantitative approach. The aim of the present work is to study more theoretically the efficiency of multi-shell spherical sampling in diffusion MRI, in order to gain understanding in HYDI-like approaches, possibly paving the way to further optimization strategies.

研究の動機と目的

• マルチシェル拡散MRIにおける角度サンプリングの制限に関する定量的理解の不足に対処すること。
• 球面サンプリングが不十分な場合にq空間データにどのような aliasing アーチファクトが生じるかを調査すること。
• HYDIに類似した拡散MRI手法における球面サンプリング戦略の最適化のための理論的基盤を提供すること。
• 球面サンプリングの効率を調和解析および信号再構成の忠実度の観点から評価するフレームワークを確立すること。

提案手法

• 著者らは、球面上での拡散信号をモデル化するために球面調和（SH）解析を用い、効率的な畳み込み演算の実現にSH変換を活用する。
• ODF計算のための主要な演算であるファンク＝ラドン変換を、離散的球面サンプリングの文脈で分析する。
• 調和解析を用いて、信号再構成に及ぼすアンダーサンプリングの影響を定量化し、特に不十分な角度分解能に起因する aliasing に注目する。
• スペクトル漏れや再構成誤差を回避するための最小サンプリング要件を特定する理論的境界を導出する。
• フレームワークをマルチシェル測定に適用し、異なるサンプリング方式がODF推定の正確性に与える影響を評価する。
• 球面調和関数およびサンプリング定理の理論に裏付けられた分析により、実用的最適化のための数学的基盤を提供する。

実験結果

リサーチクエスチョン

• RQ1 マルチシェル拡散MRI測定において、q空間データの aliasing を回避するための最小角度サンプリングレートは何か？
• RQ2 球面サンプリング方式の選択が、方向分布関数（ODF）再構成の正確性にどのように影響するか？
• RQ3 不十分なサンプリングを用いた場合、aliasing アーチファクトは拡散信号推定の品質をどの程度劣化させるか？
• RQ4 拡散MRIにおける球面サンプリングの効率を予測・最適化できる理論的フレームワークを確立できるか？
• RQ5 拡散信号のスペクトル特性は、球面上での必要サンプリング密度にどのように影響するか？

主な発見

• 角度サンプリングが不十分な場合、q空間データがアンダーサンプリングされていると、再構成された拡散信号に避けがたい aliasing アーチファクトが生じる。
• 球面調和変換は、ファンク＝ラドン変換を介してODFを効率的かつ正確に計算するための方法を提供するが、これはサンプリング分解能が十分である場合に限る。
• 理論的分析により、元の拡散信号のバンド幅を捉えるために十分なサンプリングレートが必要であることが明らかになった。これによりスペクトル漏れを防げる。
• 本研究は、サンプリング密度と再構成忠実度の間の定量的関係を確立し、aliasing の程度を予測可能にする。
• 結果として得られた理論的基盤により、測定時間を短縮しつつ信号の正確性を維持できるマルチシェル拡散MRIプロトコルの最適化が可能になる。
• 本フレームワークにより、サンプリング方式の体系的評価が可能となり、拡散MRIにおけるより効率的な測定シーケンスの設計を支援する。