[論文レビュー] Weighted Diffeomorphic Density Matching with Applications to Thoracic Image Registration

Caleb Rottman, Martin Bauer+2 more|arXiv (Cornell University)|Aug 11, 2015

Medical Imaging Techniques and Applications参考文献 7被引用数 3

ひとこと要約

本稿では、微分同相の密度作用と、ヤコビアン行列式に空間的に変化するペナルティを組み合わせることで、質量保存と空間的に変化する組織の圧縮性を強制することにより、計算的に効率的な解剖的胸郭CT画像登録のための微分同相的密度マッチングフレームワークを提案する。この手法は、標準LDDMMに比べて質量を保存し、体積変化を圧縮可能な組織に限定することで、10倍の高速化（1分間に400反復対45反復）を達成する。

ABSTRACT

In this article we study the problem of thoracic image registration, in particular the estimation of complex anatomical deformations associated with the breathing cycle. Using the intimate link between the Riemannian geometry of the space of diffeomorphisms and the space of densities, we develop an image registration framework that incorporates both the fundamental law of conservation of mass as well as spatially varying tissue compressibility properties. By exploiting the geometrical structure, the resulting algorithm is computationally efficient, yet widely general.

研究の動機と目的

• 標準LDDMMが質量を保存しないことと、組織固有の圧縮性をモデル化しないという限界に対処する。
• 呼吸周期中に肺が高圧縮性であり、軟部組織・骨がほぼ非圧縮性であるという、胸郭解剖の物理的現実を再現する。
• 質量保存と解剖的組織特性を尊重する計算的に効率的な登録フレームワークを開発する。
• 微分同相のリーマン幾何学と密度を統合し、安定で幾何学的に一貫した変換を可能にする。
• 動的胸郭画像における変形推定の向上により、放射線オンコロジーにおける正確な4次元治療計画を実現する。

提案手法

• 微分同相が体積形式（すなわち密度）に作用する引き戻し作用を用いて、画像登録を密度マッチング問題として定式化し、質量保存を保証する。
• 密度空間におけるフィッシャー＝ライアの度量を用いてリーマン構造を定義し、微分同相の幾何学と密度変換を結びつける。
• 微分同相のヤコビアン行列式に空間的に変化するペナルティ関数 $ f(x) $ を適用し、$ f(x) $ は初期CT強度をもとにロジスティック関数によって導出され、組織の圧縮性を反映する。
• Sobolev $ H^1 $ 度量を用いて微分同相空間における測地線路を構築し、マッチング適合性と滑らかさのバランスを取るエネルギー汎関数を最小化する。
• PyCAソフトウェアパッケージを用いてGPU上でアルゴリズムを実装し、測地線方程式の統合を回避することで顕著な高速化を実現する。
• 推定された微分同相による呼気画像のプッシュフォワードを用いて、吸入画像を合成し、真値と直接比較可能にする。

実験結果

リサーチクエスチョン

• RQ1 大規模な肺容量変化が生じる状況下でも、質量を保存する微分同相的登録フレームワークを設計可能か？
• RQ2 肺の高圧縮性と骨の非圧縮性といった、空間的に変化する組織の圧縮性を、画像登録において効果的にモデル化できるか？
• RQ3 標準的な $ L^2 $ 画像作用の代わりに微分同相の密度作用を用いることで、より正確で物理的に妥当な変形場が得られるか？
• RQ4 ヤコビアン行列式に空間的に変化するペナルティを適用することで、登録精度を向上させつつ計算効率を維持できるか？
• RQ5 動的胸郭画像データにおいて、本手法は標準LDDMMに比べて、速度、正確性、解剖学的妥当性の観点でどの程度優れているか？

主な発見

• 提案手法は標準LDDMMに比べ10倍の高速化を達成し、1分間に400反復（LDDMMは45反復）の速度を実現した。
• 空間的に変化するペナルティを用いた密度作用により、体積変化が主に肺および体外に限定され、非圧縮性組織での不自然な圧縮・膨張が回避された。
• 合成された吸入画像と真の吸入画像との絶対差がLDDMMに比べ顕著に減少し、肺領域においては最小限の誤差が生じた。
• 空間的に変化するペナルティを導入しない場合、LDDMMに比べ改善は見られるが、肺以外の領域で物理的に不自然な体積変化が生じ、アダプティブペナルティの必要性が明確になった。
• ヤコビアン行列式マップから、本手法が高密度構造（例：骨）で非圧縮性を維持するとともに、低密度領域（例：肺）では圧縮性を許容していることが確認された。
• エネルギープロットとロジスティックペナルティ関数から、安定した収束が示され、初期CT強度に基づいて圧縮可能な組織と非圧縮可能な組織を効果的に区別していることが確認された。