Microlocal analysis of non-linear artifacts in cone beam CT

技術要約：コーンビームCTにおける非線形アーティファクトのマイクロローカル解析

問題提起
本論文は、コーンビームX線コンピュータ断層撮影（CT）において発生するビーム硬化アーティファクトのマイクロローカル解析を扱う。標準的な線形CTモデルとは異なり、CTにおける物理的なデータ収集は、X線源の多色性と材料のエネルギー依存的な減衰能に起因して、本質的に非線形である。順モデルはビール・ランベルトの法則に従い、測定データは、減衰係数 $\mu(x, E)$ の線積分を含む指数項を、光源のエネルギースペクトルで積分したものとなる。

著者らは、X線源の集合が $\mathbb{R}^3$ 内の1次元の滑らかな曲線 $\gamma$ に限定される特定のシナリオ（これはコーンビームCTの標準的なモデルである）に焦点を当てている。本研究の重要な点は、光源スペクトルが狭いエネルギービン内に限定されているという（先行文献でよく見られる簡略化された）制約的な仮定に依存していないことである。代わりに、光源スペクトル $P(E)$ が正の全質量を持つ非負の可積分関数であることのみを仮定している。再構成対象である減衰係数 $\mu$ は、表面 $S$（物体の境界を表す）を横切る際にジャンプ不連続性を持ち、それ以外の場所では滑らかであるとモデル化されている。主な問題は、この非線形な指数項が、線形X線変換には存在しない特異性をどのように生成し、それがフィルター補正逆投影法（FBP）を通じて再構成画像へとどのように伝播するかを特徴付けることである。

手法
著者らは、特異性の性質と強度を厳密に記述するために、共ノルマル分布（conormal distributions）の理論とマイクロローカル解析を用いている。手法は以下のステップを通じて進行する：

1. 物理モデリング： データ $g(L)$ は $-\log \int P(E) \exp(-\int_L \mu(x, E) dx) dE$ としてモデル化される。著者らは、この非線形演算子 $R\mu$ の、タンジェンシー（接線）を持つレイ（光線）に対応するシノグラム内の点付近における特異性を分析している。
2. マイクロローカル分解：
   • 単一接線レイ（Single Tangent Rays）： $S$ に対して単一点で接するレイについて、著者らは非線形データが局所的に共ノルマル分布として振る舞うことを示している。彼らは、この非線形性が平方根型の特異性を導入し、結果として $\Gamma$（接線レイの多様体）に対する分布 $\dot{I}^{-3/2}(\Gamma)$ を生成することを証明している。
   • 二重接線レイ（Double Tangent Rays）： 解析の核心は、2つの異なる点で同時に $S$ に接するレイに焦ブルースに焦点を当てている。この場合、非線形項は2つの共ノルマル分布（それぞれの接点に対応するもの）の積を生成する。対ペア・ラグランジュ分布（paired Lagrangian distributions）と分布の積（補題2.17）の理論を用いることで、著者らはこの積が、線形変換には存在しない「新しい」特異性をデータ内に生成することを証明している。これらの新しい特異性は、2つの接線レイの多様体の交差部分に位置する。
3. 再構成への伝播： 著者らは、これらの特異性に随伴演算子 $R^*$（FBP再構成 $R^*P$ の一部）を適用する。順演算子と随伴演算子の正準関係（canonical relations）の合成を分析することで、画像ドメインにおけるこれらのアーティファクトの位置と強度を決定する。
4. 幾何学的条件： マイクロローカル解析の妥当性を保証するために、特にボルカー条件（Bolker condition）に関して、特定の幾何学的条件が必要となる。著者らは仮定4.1および仮定5.1を導入しており、これらは本質的に、光源の軌跡 $\gamma$ への接平面が、物体の表面 $S$ の接平面と一致して退化を生じさせないことを要求している。彼らは、これらの条件が生成的に満たされると主張している。
5. 非滑らかな表面への拡張： 理論は、$S$ が「リッジ（稜線）」（例：立方体のエッジのように、2つの滑らかな表面パッチが横断的に出会う場所）を含むケースへと拡張される。著者らは「一般化された接線性（generalized tangency）」を定義し、解析が成立することを示すが、その際、接線の種類（滑らかな表面かリッジか）に応じて、結果として生じる分布の次数が変化する。

主要な貢献と結果

• アーティファクトの型の特徴付け： 本論文は、コーンビームCTにおけるビーム硬化アーティファクトが局所的に共ノルマル型であることを証明している。具体的には、非線形性に起因するアーティファクトは、光源の軌跡 $\gamma$ と交差する二重接線レイの和集合である2次元曲面 $L_a$ 上に存在する。
• アーティファクト強度の定量化： 共ノルマル分布のソボレフ次数を用いた強度の定量化は、重要な貢献である。
  • 望ましい特異性（$\mu$ の再構成されたジャンプ）は、線形変換に対応し、特定のソボレフ正則性を持つ。
  • 非線形アーティファクトは、望ましい特異性よりも1次分だけ滑らかである（ソボレフ尺度において）ことが示されている。例えば、望ましい特異性が $H^\alpha$（$\alpha < 1/2$）にある場合、非線形アーティファクトは $H^\beta$（$\beta < 3/2$）にある。これは、アーティファクトは存在するものの、理論的には本来の物体のエッジよりも弱いことを意味する。
• スペクトル仮定の緩和： マイクロローカル解析を行うために光源スペクトルが狭い区間に支持されている必要があるとする従来の著作（[25, 3]など）とは異なり、本論文は、一般的な可積分スペクトル $P(E)$ に対しても解析が成立することを実証している。
• リッジの特異性： 著者らは、非滑らかな表面（リッジ）を持つ表面へと理論を一般化し、非線形アーティファクトが依然として特定の曲面上に共ノルマル分布を形成することを示し、接線のタイプ（滑らかな表面かリッジか）に基づいて分布の次数を調整している。
• 検証： 理論的予測は、円錐形コーンビーム幾何学における金属物体（金および銀の球体）のシミュレーション再構成を通じて検証されている。シミュレーションは、非線形再構成において、予測された曲面 $L_a$（二重接線の和集合）に沿ったストリーク・アーティファクトの存在を確認しており、これらは線形再構成では見られないものである。また、シミュレーションは、これらのアーティファクトが本来の物体のエッジよりもかすかで、かつボルカー条件の違反によって生じるアーティファクト（制限されたレイ変換の幾何学自体から生じるもの）とは明確に異なることを示している。

意義
本論文は、実用的な設定（制限されたレイ変換と一般的なスペクトル）において、ヒューリスティックではない厳密なマイクロローカルな記述を提供した点に意義があると主張している。これらのアーティファクトが、真の物体のエッジとは異なる正則性クラスに属する共ノルマル分布であることを確立することで、本研究は、なぜこれらのアーティファクトが鋭いエッジではなく「かすかな」ストリークとして現れるのかを理解するための理論的基礎を提供している。この区別は、補正アルゴリズムの開発にとって極めて重要であり、アーティファクトが真の信号とは異なる正則性のクラスを占めていることを示唆している。非滑らかな境界（リッジ）への拡張は、理論の現実世界の物体への適用範囲をさらに広げている。著者らは、自らの理論が滑らかな領域や一般的なリッジには適用できる一方で、幾何学的仮定が崩れる特定の病理的なケースを構築可能であることも明記しているが、それらは実用的なアプリケーションにおいては起こりにくいものであるとしている。