[論文レビュー] Full-Wave Iterative Image Reconstruction in Photoacoustic Tomography with Acoustically Inhomogeneous Media

Chao Huang, Kun Wang+3 more|arXiv (Cornell University)|Mar 22, 2013

Photoacoustic and Ultrasonic Imaging参考文献 39被引用数 7

ひとこと要約

本稿では、音響不均一媒質における光声トモグラフィー（PACT）のための、完全波反復画像再構成法を提示する。離散的画像モデルは正確な光声波方程式に基づく。k-space擬スペクトル法を用いて整合的な前方および逆投影演算子を実装することで、ノイズ、不完全データ、音響不均一性に起因するアーチファクトを低減する強健な再構成が可能となり、シミュレーションおよび実験データの両方で検証された。

ABSTRACT

Existing approaches to image reconstruction in photoacoustic computed tomography (PACT) with acoustically heterogeneous media are limited to weakly varying media, are computationally burdensome, and/or cannot effectively mitigate the effects of measurement data incompleteness and noise. In this work, we develop and investigate a discrete imaging model for PACT that is based on the exact photoacoustic (PA) wave equation and facilitates the circumvention of these limitations. A key contribution of the work is the establishment of a procedure to implement a matched forward and backprojection operator pair associated with the discrete imaging model, which permits application of a wide-range of modern image reconstruction algorithms that can mitigate the effects of data incompleteness and noise. The forward and backprojection operators are based on the k-space pseudospectral method for computing numerical solutions to the PA wave equation in the time domain. The developed reconstruction methodology is investigated by use of both computer-simulated and experimental PACT measurement data.

研究の動機と目的

• 音響不均一媒質における既存のPACT再構成手法の限界を解消すること。具体的には、弱く変化する媒質の仮定や高い計算コストの問題を含む。
• 実際のPACT応用におけるデータ不完全性、ノイズ、モデル誤差に起因するアーチファクトを克服すること。
• 可変な音速、密度、減衰を有する不均一媒質における波動伝搬を正確に表現できる、正確な光声波方程式に基づく離散的画像モデルの構築。
• 複雑な媒質における反復再構成アルゴリズムの効率的実装を可能にする、整合的前方および逆投影演算子ペアの確立。
• プリマルデュアルハイブリッド勾配法に基づく現代の反復再構成技術を、現実的な音響不均一性を有する完全波PACTモデルに適用可能にする。

提案手法

• 可変音速、質量密度、周波数-べき乗則減衰を組み込んだ、正確な光声波方程式に基づく離散的画像モデルの構築。
• PA波方程式の時間領域数値解法として、k-space擬スペクトル法を用いて前方および逆投影演算子を実装。
• 反復再構成アルゴリズムの効率的適用を可能にする、整合的演算子ペア（前方、逆投影、およびその随伴演算）の開発。
• 不完全またはノイズ混在データ条件下でも逆問題を安定して解けるように、演算子をプリマルデュアルハイブリッド勾配（PDHG）アルゴリズムに統合。
• 主要な演算子を定義する：$\mathcal{P}_l$ は離散的発散類似演算を計算するためのもの、$\mathcal{P}_c$ は非負性制約のためのもの、$\mathcal{P}_l^T$ は随伴演算のためのもの、$\mathcal{P}_p$ は勾配更新におけるベクトル正規化のためのもの。
• アルゴリズムをコンピュータシミュレーションおよび実験的PACTデータに適用し、現実的な不均一媒質における性能を検証。

実験結果

リサーチクエスチョン

• RQ1 正確な光声波方程式に基づく完全波反復再構成フレームワークは、音響不均一性が存在する場合に、PACTにおけるアーチファクトを効果的に低減できるか？
• RQ2 提案された整合的前方および逆投影演算子ペアは、データ不完全性およびノイズが存在する状況で、再構成精度と収束性をどのように向上させるか？
• RQ3 k-space擬スペクトル法は、3次元不均一媒質におけるPACTのための効率的かつ正確な波動伝搬モデリングをどの程度可能にするか？
• RQ4 強力な音響歪み（例：頭蓋骨に起因する波面歪曲）を伴う実験的データにおいて、反復再構成法はどのように性能を発揮するか？
• RQ5 現実的な測定条件下で、従来の時間反転法やレイトレーシングに基づく手法と比較して、提案手法はアーチファクト低減および空間分解能の面で優位性を示せるか？

主な発見

• 提案された完全波反復再構成法は、頭蓋骨に起因する波面歪曲に起因するアーチファクトを顕著に低減することが、実験的頭部透過PACTデータから明らかになった。
• 不完全またはノイズ混在データの状況においても、従来の時間反転法やレイトレーシングに基づく手法と比較して、画像品質および空間分解能が向上した。
• k-space擬スペクトル法の使用により、強力な音響不均一性を有する3次元媒質でも、PA波方程式の正確で安定した数値解法が実現された。
• 整合的前方および逆投影演算子ペアのおかげで、反復アルゴリズムの効率的実装が可能となり、計算負荷は低減しながらも高い再構成忠実度を維持した。
• PDHGに基づくアルゴリズムは、スパースまたはノイズ混在測定データに対しても安定して収束し、正確な再構成を生成した。これは、シミュレーションおよび実験データの両方で検証された。
• 正確な波動方程式に基づくアプローチのおかげで、音速や減衰の変動といったモデル誤差や物理的パrameterの不確実性に対しても、本手法は頑健であることが示された。