[論文レビュー] On an inverse boundary problem arising in brain imaging

Youjun Deng, Hongyu Liu+1 more|arXiv (Cornell University)|Feb 1, 2017

Numerical methods in inverse problems参考文献 19被引用数 5

ひとこと要約

本稿は、時間領域Maxwell系に基づく非線形逆問題枠組みを用いて、境界での電場および磁場測定値から、脳に類似した媒質内の電流源および区分的定数な電磁パラメータ（誘電率、透磁率、導電率）の両方の固有の回復を確立する。著者らは、一般的な条件下で同定可能性を証明し、特に電流源が方向に依存しないか、調和的である場合に、EEG/MEG逆画像化における完全な媒質特性化を進める。

ABSTRACT

We consider the inverse problem of recovering both an unknown electric current and the surrounding electromagnetic parameters of a medium from boundary measurements. This inverse problem arises in brain imaging. We show that under generic conditions one can recover both the source and the electromagnetic parameters if these are piecewise constant and the current source is invariant in a fixed direction or a harmonic function.

研究の動機と目的

• 脳イメージングにおける境界測定値から未知の電流源と周囲の電磁パラメータを同時に回復する逆問題に対処すること。
• 電流源と媒質パラメータの両方を含む逆問題の同定可能性および一意性を確立する十分条件を提示すること。
• EEG/MEGにおける既存の逆源問題を拡張し、電流源 J とともに誘電率 ε、透磁率 μ、導電率 σ の媒質パラメータの回復を組み込むこと。
• 電流源が固定方向に依存しないか、調和的である場合の解析を、区分的定数な電磁パラメータのもとで行うこと。
• 媒質を未知とし、同時に回復可能であるとみなすことで、より良い脳イメージング技術の理論的基盤を提供すること。

提案手法

• 有界領域 B ⊂ ℝ³ における因果的電流源および区分的定数な電磁パラメータを伴う時間領域Maxwell系を用いて逆問題を定式化する。
• 境界 ∂B × ℝ₊ における電場および磁場の測定値を未知パラメータへマップする境界測定作用素 Π を定義する。
• 層間ポテンシャル技術および境界積分方程式を用いてシステムを解析し、特に背景グリーン関数を用いた単層および二重層ポテンシャルを用いる。
• 跳躍関係およびトレース作用素を適用して境界上での積分方程式を導出し、未知数 g, h, l, m が ∂B および ∂Σ₁ 上に現れる方程式系に帰着させる。
• 得られた境界積分方程式系の解の一意性と、その構造を用いて、すべてのパラメータが一致しない限り矛盾が生じることを示す。
• 導出された方程式系の可解性および一意性を保証するため、二層媒質構造に適正性条件を課す。

実験結果

リサーチクエスチョン

• RQ1 境界での電場および磁場測定値から、電流源および電磁パラメータ（ε, μ, σ）を両方とも一意に回復できるか？
• RQ2 電流源と媒質パラメータを同時に回復する逆問題が同定可能となる条件は何か？
• RQ3 電流源が任意ではなく、固定方向に依存しないか、調和的である場合、一意性の結果は保持されるか？
• RQ4 媒質パラメータの区分的定数構造が、逆問題の同定可能性にどのように影響するか？
• RQ5 同様の一意性保証が得られるように、理論枠組みを N 層の区分的定数媒質に拡張できるか？

主な発見

• 区分的定数パラメータおよび方向に依存しないか調和的な電流源を仮定した一般的な条件下で、電流源 𝒥 および電磁パラメータ（ε, μ, σ）の両方の逆問題は一意に同定可能である。
• 証明は境界積分方程式技術に依拠し、得られた方程式系の一意的可解性に基づき、すべてのパラメータが一致しない限り矛盾が生じることを示す。
• 適正性条件 (3.45) を満たす二層媒質構造では、境界測定値から電流源および媒質パラメータが一意に回復可能である。
• この手法は N 層の区分的定数媒質に拡張可能であるが、適正性条件は N > 2 の場合に複雑化し、明示的に示されていない。
• 条件 (3.45) は一般的であり、球対称幾何における検証が可能であることが示され、理論的結果の実用的妥当性を支持する。
• 解析により、境界測定作用素 Π が (𝒥, ε, μ, σ) の完全な構成を一意に決定することが確認され、媒質特性が未知である場合のEEG/MEGイメージングにおける長年の課題が解決された。