2025-11-19T04:04:12.889579

Sampling Density Compensation using Fast Fourier Deconvolution

Luo, Hu, Qi

Density Compensation Function (DCF) is widely used in non-Cartesian MRI reconstruction, either for direct Non-Uniform Fast Fourier Transform (NUFFT) reconstruction or for iterative undersampled reconstruction. Current state-of-the-art methods involve time-consuming tens of iterations, which is one of the main hurdles for widespread application of the highly efficient non-Cartesian MRI. In this paper, we propose an efficient, non-iterative method to calculate DCF for arbitrary non-Cartesian $k$-space trajectories using Fast Fourier Deconvolution. Simulation experiments demonstrate that the proposed method is able to yield DCF for 3D non-Cartesian reconstruction in around 20 seconds, achieving orders of magnitude speed improvement compared to the state-of-the-art method while achieving similar reconstruction quality.

academic

高速フーリエ逆畳み込みを用いたサンプリング密度補正

基本情報

論文ID: 2510.14873
タイトル: Sampling Density Compensation using Fast Fourier Deconvolution
著者: Rui Luo, Peng Hu, Haikun Qi (ShanghaiTech University)
分類: physics.med-ph
発表日: 2025年10月16日 (arXivプレプリント)
論文リンク: https://arxiv.org/abs/2510.14873

要約

密度補正関数(DCF)は、非直交MRI再構成において、直接的な非均一高速フーリエ変換(NUFFT)再構成および反復アンダーサンプリング再構成の両方で広く応用されている。現在の最先端手法は数十回の反復を要する耗時的なプロセスが必要であり、これが効率的な非直交MRI広範応用の主要な障壁の一つとなっている。本論文は、高速フーリエ逆畳み込みを用いて任意の非直交k空間軌跡に対するDCFを計算する効率的な非反復手法を提案する。シミュレーション実験により、本手法が約20秒以内に3D非直交再構成用のDCFを生成でき、最先端手法と比較して数桁の速度向上を実現しながら、同等の再構成品質を維持することが示された。

研究背景と動機

問題定義

非直交MRIサンプリングは低周波領域でのサンプリング密度が高周波領域よりはるかに高く、適切な重み付け調整なしに直接再構成すると画像ぼやけが生じる。密度補正関数(DCF)はサンプリング密度を平衡化するために使用され、非直交MRI再構成の重要な構成要素である。

問題の重要性

再構成品質: DCFは直接NUFFT再構成と反復再構成の両方に不可欠である
収束加速: DCFは問題の条件数を改善することで反復再構成と深層学習再構成の収束を加速できる
臨床応用: 効率的なDCF計算は非直交MRI臨床応用の重要なボトルネックである

既存手法の限界

ボロノイ図手法: 計算コストが高く数値的に不安定であり、特に3D軌跡に対して問題がある
反復手法: PipeとMenonが提案した古典的反復手法は数十回の反復を必要とし、各反復に数秒から数分を要する
最適化手法: より正確であるが、実行時間は反復手法より2桁増加する

研究動機

特に3D非直交サンプリングパターンを対象とした、高速で非反復的なDCF計算手法を開発し、計算効率のボトルネックを突破すること。

核心的貢献

高速フーリエ逆畳み込み(FFD)に基づく非反復DCF計算手法を提案
1～2桁の速度向上を実現、3D軌跡DCF計算時間を約10分から20秒以内に短縮
再構成品質を維持または若干改善しながら、より滑らかなDCFを生成
汎用ソリューションを提供、任意の非直交k空間軌跡に適用可能
オープンソース実装を提供、手法の再現性と広範な応用を促進

方法の詳細

タスク定義

非直交k空間サンプリングパターン $K = \{k_i\}_{i=1}^{N_k}$ が与えられたとき、密度補正関数 $D(k)$ を求め、加重サンプリングパターンの点拡がり関数(PSF)が視野内でインパルス関数に近似するようにする: $P(x) \approx \delta(x)$ for $\|x\| < L$ 。

核心的理論フレームワーク

1. 問題のモデル化

サンプリングプロセスは以下のように表現される: $S_1(k) = \text{III}(k) \cdot S_0(k)$

ここで $\text{III}(k) = \sum_{i=1}^{N_k} \delta(k - k_i)$ はインパルス列である。

加重後のk空間: $S_2(k) = D(k) \cdot \text{III}(k) \cdot S_0(k)$

加重サンプリングパターン(WSP): $E(k) = D(k) \cdot \text{III}(k)$

対応するPSF: $P(x) = \mathcal{F}^{-1}\{E(k)\}$

2. PSF分解戦略

初期PSF推定値 $\hat{P}(x)$ を以下のように分解する: $\hat{P}(x) = \hat{P}_{\text{in}}(x) + \hat{P}_{\text{out}}(x)$

ここで:

$\hat{P}_{\text{in}}(x) = \hat{P}(x) \cdot W(x)$
$\hat{P}_{\text{out}}(x) = \hat{P}(x) \cdot (1-W(x))$

$W(x)$ はウィンドウ関数であり、 $\|x\| \geq L$ のとき0である。

3. 最適加重サンプリングパターン

逆畳み込みにより最適WSPを得る: $E^*(k) = \hat{E}(k)/\hat{E}_{\text{in}}(k)$

これにより以下が実現される:

$P^*_{\text{in}}(x) = \mathcal{F}^{-1}\{1\} = \delta(x)$
$P^*_{\text{out}}(x) = \mathcal{F}^{-1}\{\hat{E}_{\text{out}}(k)/\hat{E}_{\text{in}}(k)\}$

技術的革新点

1. ウィンドウ関数の最適化

パラメータ化ウィンドウ関数 $W(x) = 1 - \|\bar{x}\|^p$ を採用し、ここで $\bar{x} = x/L$ 、最小-最大パラメータ探索により最適形状パラメータを決定する: $p^* = \arg\min_p \{\max_{i_{\text{test}}} \|P^*_{\text{out}}(x)\|/P^*_0\}$