An angular distortion matrix approach for joint wave-speed tomography and aberration correction in scattering media

技術要約：結合波速トモグラフィーおよび収差補正のための角度歪度行列アプローチ

問題提起
反射イメージング（超音波、光学、地震波）において、正確な反射率の再構成は、波速分布に支配される精密な波動伝搬モデルに依存している。しかし、複雑な散乱媒体においては、波速が事前に既知であることは稀である。想定された波速と実際の波速の不一致は位相収差を導入し、フォーカシングを劣化させ、点広がり関数（PSF）を広げ、散乱体の見かけ上の位置を歪ませる。これらの収差は、基礎となる波速の不均一性に関する情報をエンコードしているが、それらを媒体の反射率から分離することは困難であり、特にスペックルが支配的な媒体においては顕著である。既存の手法は根本的な結合問題に直面している。すなわち、収差を補正するには正確な波速モデルが必要であるが、このモデルを取得するには、収差によってバイアスを受けていない観測量が必要であるという点である。さらに、従来の歪度行列アプローチは、収差が入射経路または反射経路のいずれかに起因すると仮定しているが、実際には両方の寄与が本質的に混在しており、速度モデルを更新しない限り修正できない深度シフトを引き起こす。

手法
著者らは、波速分布を反復的に推定し、新しい角度歪度行列（$D_{\theta\theta}$）を用いて収差を補正する統一されたフレームワークを提案している。手法は以下のように進む：

1. 行列イメージングと周波数・波数表現：
   プロセスは、反射行列 $R_{uu}$（トランスデューサ間応答）を取得することから始まる。これは、フォーカスされた基底へと投影され、フォーカスされた反射行列 $R_{xx}$ を形成する。歪みを分離するために、著者らは周波数・波数表現の波の歪み $D_{kk}$ を導入する。この演算子は、任意の点 $x$ から散乱が発生したと仮定した際の、応答から期待される幾何学的位相成分を除去することによって構築される。これにより、ポストプロセッシングにおいて、あらゆる仮想的な実験に対する歪みの合成が可能となる。

2. 角度歪度行列の構築：
   $D_{kk}$ 演算子は、ガウスビーム基底上に投影され、角度歪度行列 $D_{\theta\theta}$ を形成する。この基底は、入射角（$\theta_{in}$）と反射角（$\theta_{out}$）によって定義されるステアード・ガウスビームを利用している。この投影は、角度メモリー効果を活用している。この効果により、波は無秩序な状態においても、特定の角度方向に沿って強い相関を示す。$D_{\theta\theta}$ の反対角成分（一定の中間角 $\theta_m = (\theta_{in} + \theta_{out})/2$）を解析することで、明確に定義された伝搬経路に関連する位相収差を分離できる。

3. 局所収差推定：
   視野（FOV）は、重なり合う「等方性パッチ（isoplanatic patches）」に分割される。各パッチ内では、角度歪度行列の要素から局所相関行列が計算される。波の歪みが小さな領域内で不変であるという仮定の下、この相関行列の第1固有ベクトルは、局所的な位相収差則（$u_m$）の推定値を提供する。位相エイリアシングを避けるため、手法は連続する半差角（$\theta_d$）間の相対的な位相差を利用する。

4. 波速トモグラフィーのための反復反転：
   推定された位相収差は、音速分布の再構成に使用される。問題は、測定された位相差と、伝搬モデルによって予測される位相差との間のミスマッチを最小化する最適化タスクとして定式化される。反転は以下の手順で反復的に実行される：

   • 初期の均質なモデル（例：$c_0 = 1540$ m/s）を仮定する。
   • 収差が推定され、補正される。
   • 線形化された逆問題（ヤコビアンには直線近似を用い、チホノフ正則化を使用）を用いて、スローネス摂動（$\Delta s$）を解く。
   • 波伝搬モデルが更新され、残留収差が無視できるほど小さくなるまでプロセスを繰り返す。

5. 多重散乱の処理：
   メモリ効果が劣化する高度に不均質な媒体（例：乳房組織）においては、反射行列に対して適応型コンフォーカルフィルタが適用される。これにより、単一散乱に関連する対角成分のエネルギーを保持しつつ、非干渉的なバックグラウンドノイズを抑制する。

主な結果
本アプローチは、3つのシナリオにおける超音波実験を用いて検証された：

• 組織模倣ファントム： 既知の層状構造（脂肪、筋肉、肝臓）と変化する音速（1420–1550 m/s）を持つファントムを使用した。反復法は、平均誤差が真値の1%以内となる音速マップの再構成に成功した。ストレル比（フォーカシング品質の尺度）は0.34から0.96へと向上し、反射率画像は鋭い界面と強化されたスペックルコントラストを示した。
• 生体内乳房イメージング： 浸潤癌患者に適用された。本手法は、周囲の脂肪および線維腺組織と区別される、局所的な高音速領域（~1574 m/s）を特定した。収差補正により、未補正の画像では解剖学的特徴の見かけ上の深度が偏っていたことが明らかになった。
• 生体内肝臓イメージング： 健康なボランティアに適用された。本手法は、腹壁（皮膚、脂肪、筋肉）の層状構造と肝実質を高軸方向分解能で解明した。回収された音速は健康な組織の文献値と一致しており、ストレル比は0.61から0.93へと上昇し、反射率画像を大幅に鮮明化させ、無エコーの静脈などの以前は不明瞭であった構造を明らかにした。

意義および主張
本論文は、角度歪度行列アプローチが、波速トモグラフィーと収差補正の結合問題に対する堅牢な解決策を提供すると主張している。その主な意義は以下の通りである：

• 統合： 歪度行列の概念（通常、収差補正に使用される）と、エコーモードにおける計算超音波トモグラフィー（CUTE）を、単一の反復的な行列イメージングフレームワーク内で統合している。
• モデル独立性： 角度歪度行列の定式化は、フォーカシングに使用される特定の波動伝搬モデルに依存しておらず、手法自体を漸進的に洗練させることが可能である。
• 高解像度： ガウスビーム基底への投影により、古典的な歪度行列アプローチと比較して、より小さな等方性パッチ（約1 mm）にわたって収差を推定することができ、トモグラフィーにおける高い空間分解能を実現している。
• 臨床的有用性： 本手法は、定量的な音速マップ（病理学のバイオマーカーとなり得る）を提供すると同時に、複雑で不均質な組織における構造画像の品質を向上させる能力を示している。
• 汎用性： 超音波で実証されているが、本フレームワークは、反射行列イメージングが実装可能なあらゆる波動様式（光学、地震波）に適用可能である。