[論文レビュー] Inverse anisotropic conductivity from power densities in dimension $n\ge 3$

François Monard, Guillaume Bal|arXiv (Cornell University)|Aug 29, 2012

Numerical methods in inverse problems参考文献 17被引用数 10

ひとこと要約

本稿では、次元 $ n \geq 3 $ における内部の電力密度測定値から、完全に異方的な伝導度テンソルを構成的・明示的に再構成する手法を提示する。著者らは、線形独立な勾配を持つ $ n $ 重の解と電力密度の偏光を活用し、スカラー倍率 $ \tau $ と異方的構造 $ \tilde{\gamma} $ の両方を回復する。スカラー倍率 $ \tau $ についてはスケールの損失がなく、$ \tilde{\gamma} $ については1階の正則性損失を伴う局所的再構成が、関数の汎関数のランク最大化条件のもとで達成される。

ABSTRACT

We investigate the problem of reconstructing a fully anisotropic conductivity tensor $γ$ from internal functionals of the form $ abla u\cdotγ abla u$ where $u$ solves $ abla\cdot(γ abla u) = 0$ over a given bounded domain $X$ with prescribed Dirichlet boundary condition. This work motivated by hybrid medical imaging methods covers the case $n\ge 3$, following the previously published case $n=2$ \cite{Monard2011}. Under knowledge of enough such functionals, and writing $γ= τ ilde γ$ ($\det ildeγ= 1$) with $τ$ a positive scalar function, we show that all of $γ$ can be explicitely and locally reconstructed, with no loss of scales for $τ$ and loss of one derivative for the anisotropic structure $ ildeγ$. The reconstruction algorithms presented require rank maximality conditions that must be satisfied by the functionals or their corresponding solutions, and we discuss different possible ways of ensuring these conditions for $\C^{1,α}$-smooth tensors ($0

研究の動機と目的

• 内部の電力密度測定値 $ \nabla u \cdot \gamma \nabla u $ から完全に異方的な伝導度テンソル $ \gamma $ を再構成する逆問題を、ハイブリッド医療画像法の動機づけのもとで取り扱う。
• 2次元における異方的伝導度の再構成フレームワークを高次元へ拡張し、明示的かつ安定した回復を保証する。
• 正則性損失の観点から再構成の安定性を特徴づける：スカラー係数 $ \tau $ については損失なし、正規化された異方的構造 $ \tilde{\gamma} $ については1階の損失。
• 再構成アルゴリズムの可逆性を保証するための、解またはその勾配に関するランク最大化条件を同定・分析する。
• 関数の $ C^{1,\alpha} $-滑らかさ（$ 0 < \alpha < 1 $）の設定において、適切な境界条件を用いることで、ランク最大化仮定を満たす十分条件を提示する。

提案手法

• 境界領域 $ X \subset \mathbb{R}^n $ において $ \nabla \cdot (\gamma \nabla u) = 0 $ を満たす $ u $ を解とするとき、$ \gamma = \tau \tilde{\gamma} $（ただし $ \det \tilde{\gamma} = 1 $）の再構成問題を、電力密度汎関数 $ \mathcal{H}_\gamma[u](x) = \nabla u \cdot \gamma \nabla u $ の知識に基づいて定式化する。
• 相互電力密度 $ \nabla u \cdot \gamma \nabla v $ を用いた偏光を導入し、解の勾配を含む方程式系の構築を可能にする。
• 対称行列の空間における一般化された外積構成を用いて、解の勾配の $ n $ 重組み合わせから異方的構造 $ \tilde{\gamma} $ を抽出する。
• 特に、行列による左乗法の下での一般化外積の変換則を適用し、基底変換における安定性の挙動を導出する。
• 電力密度データと解の勾配のランク最大化条件に基づいて導出された方程式系を解くことにより、明示的な再構成公式を導出する。
• 再構成がデータに依存する様子を分析することで安定性推定を確立し、$ \det \gamma $ についてはリプシッツ安定性、$ \tilde{\gamma} $ についてはホルダー型安定性を示す。

実験結果

リサーチクエスチョン

• RQ1 十分なデータが与えられた場合、次元 $ n \geq 3 $ における内部電力密度測定値から完全に異方的な伝導度テンソル $ \gamma $ を明示的に再構成できるか？
• RQ2 異方的構造 $ \tilde{\gamma} $ の再構成における正則性損失は何か？スカラー係数 $ \tau $ の安定性と比較するとどうなるか？
• RQ3 解 $ u $, $ v $ などの条件下で、逆問題のためのランク最大化条件を事前に保証できるか？
• RQ4 データの摂動に対して再構成アルゴリズムはどのように振る舞い、どのような安定性推定が得られるか？
• RQ5 $ C^{1,\alpha} $-滑らかさの設定において、適切なディリクレ境界条件を選び、ランク最大化条件を満たすことができるか？

主な発見

• 著者らは、2次元の結果を拡張し、次元 $ n \geq 3 $ における内部電力密度測定値から、完全に異方的伝導度テンソル $ \gamma $ を構成的・明示的に再構成するアルゴリズムを確立した。
• 分解 $ \gamma = \tau \tilde{\gamma} $ におけるスカラー係数 $ \tau $ は正則性の損失なし（リプシッツ安定性）で回復され、異方的構造 $ \tilde{\gamma} $ は1階の正則性損失（ホルダー型安定性）を受ける。
• 再構成は、各点で $ \mathbb{R}^n $ の基底をなす勾配を持つ $ n $ 個の伝導度方程式の解が存在することに依存しており、これをランク最大化条件と呼ぶ。
• 適切なディリクレ境界条件を選び、事前にランク最大化条件を満たす十分条件を提示した。
• 対称行列空間における一般化された外積構成を用いて、異方的構造 $ \bar{\gamma} $ を抽出し、線形写像の下での変換則を用いて安定性解析を可能にした。
• 安定性推定は、電力密度汎関数が $ W^{1,\infty}(X) $ に属すると仮定して導出されており、$ 0 < \alpha < 1 $ における $ C^{1,\alpha} $-滑らかな伝導度テンソルに対しても結果が成り立つ。