技術要約：VITO – 逆トポロジー最適化を用いた血管幾何学および血流推定

問題提起

心血管疾患（CVD）の評価は、患者固有の血管幾何学から導出される壁剪断応力（WSS）や冠血流予備量比（FFR）といった血行動態記述子に大きく依存している。標準的な臨床ワークフローでは、造影剤を用いたコンピュータ断層撮影（CE-CT）を用いてこれらの幾何学構造を再構成する。しかし、一般的なパイプライン（通常、フィルタ補正逆投影法（FBP）または反復再構成とその後のセグメンテーションを伴う）は、放射線量制約による疎な測定値や生理学的運動アーチファクトを特徴とする時間分解能の高い取得データに対して、課題を抱えている。これらの制限はアンダーサンプリングによるアーチファクトを引き起こし、枝の欠落や人工的な狭窄といった幾何学的な不正確さを招く。

さらに、既存の物理制約付きの血行動態逆問題手法の多くは、固定された事前セグメント化された幾何学を前提としている。誤った固定幾何学に対して流体物理を強制すると、妥当ではあるが代表性のない流場が得られ、重要な血行動態パラメータを捉えることができない。近年、幾何学と流れを共同で再構成しようとする研究も存在するが、それらは多くの場合、生（raw）のCT投影データ（サイノグラム）から直接幾何学を復元するのではなく、測定された速度画像（例：MRI）に依存している。したがって、事前の幾何学的仮定なしに、疎でノイズの多いCE-CTサイノグラムから血管トポロジーと流体パラメータを共同で復元できるフレームワークには、依然としてギャップが存在する。

手法

著者らは、血管幾何学の復元を密度ベースのトポロジー最適化（TO）問題として定式化した、物理制約付き再構成フレームワークであるVITOを提案する。幾何学を固定された入力として扱うのではなく、血管腔を未知のデザイン変数として扱う。

コアフレームワーク

本手法は、設計変数が支配的な物理法則に従う際に、取得された疎なサイノグラムと一致するシミュレーション投影データを生成するような流体・固体分布を見出すという逆問題を解く。

1. 設計表現:

   • 計算領域はメッシュ要素に離散化され、各要素には擬似密度変数 $\gamma_e \in [0, 1]$ が割り当てられる。
   • $\gamma = 0$ は流体（血管腔）、$\gamma = 1$ は固体（組織／壁）を表す。
   • 材料特性は、材料特性の有理近似（RAMP）スキームを用いて補間される。
     • ブリンクマン・ペナルティ ($\alpha$): 固体領域における大きなペナルティは流体速度を抑制し、明示的な境界追跡なしに滑りなし条件（no-slip conditions）を効果的に強制する。
     • 拡散係数 ($D$): 固体領域における造影剤輸送を制限し、人工的な拡散を防ぐ。

2. 支配物理（フォワードモデル）:

   • 定常流: 不可圧縮ナビエ・ストークス方程式をブリンクマン・ペナルティ項を用いて解き、流体速度場 $\mathbf{u}$ をモデル化する。
   • 非定常輸送: 時間変化する流入境界条件によって駆動される、定常速度場 $\mathbf{u}$ と造影剤濃度 $c$ の過渡的な移流拡散方程式をモデル化する。
   • 前方投影: 得られた時間分解濃度場は、ベール・ランベルト則に基づく微分可能な平行ビーム投影オペレータを使用して、測定空間へとマッピングされる。これにより、予測サイノグラム $\hat{g}$ が生成される。

3. 最適化ループ:

   • 目的関数: 予測されたサイノグラムと、取得された疎な／ノイズを含むサイノグラムデータ ($g_{data}$) との間の平均二乗誤差（MSE）。
   • 設計変数: 擬似密度場 $\gamma$（幾何学）および流入速度スケーリングパラメータ $a$（流動条件）。
   • 感度解析: 本フレームワークは、エンドツーエンドの微分可能なパイプラインを実現するために JAX エコシステムを利用する。
     • 内積関数定理（IFT）: 反復的なニュートン・ラフソン・ステップの展開に伴うメモリコストを回避するため、定常ナビエ・ストークスソルバの感度を計算するために使用される。
     • チェックポインティング: 時間ステップを跨いだ逆モード自動微分（AD）中のメモリ使用量を管理するために、非定常移流拡散ソルバに採用される。
   • オプティマイザ: 設計変数を更新するために移動漸近線法（MMA）が使用される。解をバイナリの流体・固体トポロジーへと導き、局所解を回避するために、継続戦略（ペナルティおよび投影の鋭さの段階的な増加）が用いられる。

主な貢献

• 結合逆再構成: 本論文は、事前セグメント化された幾何学を必要とせず、生のCTサイノグラムから血管トポロジーと流入流体パラメータを同時に復元する統一されたフレームワークを導入する。
• 物理駆動型正則化: 逆問題を生のナビエ・ストークス方程式および過渡的移流拡散方程式と結合させることで、単なる画像事前分布ではなく、血行動態の物理学を用いて再構成を正則化する。
• 微分可能なマルチフィジックス・パイプライン: 実装は、IFTやチェックポインティングを活用して、複雑に結合された流体・輸送・投影システムの複雑な感度を扱うために、自動微分（JAX）を利用したエンドツーエンドの微分可能なパイプラインを活用している。
• データの制限に対する堅牢性: 本フレームワークは、低線量または高速取得の臨床CTに典型的な、疎なビューおよびノイズの多い取得条件下での使用を明示的に想定して設計されている。

実験結果

本手法は、3つの血管シナリオを表す合成ファントムを用いて検証された：

1. 分岐動脈: 疎なサイノグラムから分岐トポロジーと流入流体スケーリングを正常に復元した。
2. 狭窄動脈: 局所的な50%の管腔狭窄（ステノーシス）と流体パラメータを正確に再構成した。
3. 狭窄頸動脈分岐: 分岐と非対称な狭窄を組み合わせた複雑な幾何学を復元した。

性能指標:

• 精度: 再構成された幾何学は、すべてのベンチマークにおいて、グラウンドトゥルースと比較して高いダイス係数（DSC > 0.99）および低い正規化平方根平均二乗誤差（nRMSE）を達成した。
• ノイズおよび疎性への堅牢性: 極めて疎なデータ（10投影）および高いノイズレベル（25%の相対ガウスノイズ）を含む実験において、提案手法はトポロジーの一貫性（DSC = 0.9830）を維持した。対照的に、閾値処理を伴うベースラインのFBPアプローチは、同様の条件下で境界の歪み、アーチファクト、および枝の連続性の喪失に苦しんだ。
• 流体推定: 本フレームワークは、流入流体スケーリングパラメータを真値に近い1（unity）付近に復元することに成功し、境界条件とともに幾何学を推論する能力を示した。

意義および主張

本論文は、再構成プロセスに基礎となる血行動態を統合することが、ロバストな幾何学的正則化を提供し、複雑な血管構造（分岐、狭窄）を疎な臨床投影データから直接復元することを可能にすると主張している。このアプローチは、幾何学的な誤差が血行動態解析に伝播するという現在のワークフローにおける決定的な限界に対処するものである。

著者らは、本研究を疎な臨床データからの信頼できる血管解析に向けた重要な一歩として位置づけている。彼らは、現在の研究が2次元定常流および剛体境界に限定されていることを明示的に述べている。今後の課題として、フレームワークを3次元の非定常流、流体・構造相互作用（FSI）を伴う弾性壁、流入・流出位置を未知数として扱う手法、不確実性の定量化、および投影オペレータを特定のスキャナー形状（例：ファンビームCT）に適応させることが挙げられている。本論文は、即時の臨床導入を主張するものではなく、そのようなアプリケーションのための数学的および計算的な基礎を確立するものである。