技術要約：固有パラメータ校正の摂動下における椎体DRRの投影ドメイン感度解析

問題提起

正確な幾何学的校正は、透視誘導下脊椎イメージング、デジタル再構成放射線画像（DRR）生成、および2D–3D椎体レジストレーションの前提条件である。従来、校正の品質は、再投影誤差や体積再構成精度といった再構成ドメインの指標を用いて評価されてきた。しかし、多くの透視誘導型レジストレーション・アルゴリズムは、主に投影ドメインで動作しており、取得されたX線写真と合成されたDRRとの間の整合性に依存しており、再構成の正確さのみに依存しているわけではない。

近年の研究では、校正の不確かさがDRR生成に使用されるイメージング・ジオメトリに影響を与えることが示唆されているが、固有パラメータ校正の摂動が投影ドメインの整合性に与える具体的な影響については、十分に理解されていない。これは、椎体の解剖学的構造が複雑な3D構造、重なり合う骨の特徴、および限られたビューを伴う脊椎イメージングにおいて特に重要である。イメージング・ジオメトリのわずかな変化は、椎体の輪郭、椎弓根の外観、および終板の整列を変化させ、潜在的に画像類似性指標やレジストレーションの収束を低下させる可能性がある。本研究は、基礎となる解剖学および取得ポーズが固定されている場合でも、固有パラメータ校正の摂動が椎体の投影整合性とダウンストリームの2D–3Dレジストレーション精度にどのように影響するかを調査することにより、再構成ベースの校正評価と投影ドメインのレジストレーション挙動との間のギャップに対処する。

手法

著者らは、固有パラメータ校正の摂動に対する椎体透視投影の感度を分析するための、制御された合成フレームワークを提示する。

1. 合成データセットおよびジオメトリ:

• モデル: 実験では、CT由来の椎体セグメンテーションマスクおよびメッシュモデル（VerSe、SpineWeb、TotalSegmentatorなどのデータセットから取得）を使用し、単一の椎体および多椎体腰椎セグメント（L1–L5）を含む。
• イメージング設定: 非アイソセントリックなCアーム・ジオメトリを、固定された光源・検出器間距離（SDD）および光源・アイソセンター間距離（SID）を用いてシミュレートする。2つの標準的なビュー（前後方向（AP）および側方（LAT））を生成する。椎体モデルは、解剖学的運動による影響を排除するために、ビュー間で固定されたポーズに保たれる。
• DRR生成: 線積分減衰蓄積に基づく、1024 × 1024 解像度のレイ駆動型投影レンダリングを使用してDRRを生成する。

2. 摂動モデル:

• 固有パラメータ: 外在パラメータ（$R, t$）は一定に保持される一方、制御された摂動は固有パラメータ行列（$K$）に対してのみ適用される。
• 特定の摂動:
  • 焦点スケール: 有効焦点距離（$f_x, f_y$）の乗法的スケーリングを行い、拡大率の不確かさをシミュレートする。
  • 貫通点の変位: 検出器平面上での主点（$c_x, c_y$）の2Dシフト。
• 比較: グランドトゥルース（真値）のパラメータを用いて生成されたDRRを、摂動を与えられたパラメータを用いて生成されたDRRと比較する。

3. 解析指標:

• 投影ドメインの整合性: 解剖学的ランドマークの変位、輪郭ベースの距離（Cham포ー、ハウスドルフ）、シルエットの重なり（Dice係数、IoU）、および生のDRR差分ヒートマップを用いて評価する。
• ダウンストリーム・レジストレーション: ランドマークベースの2D–3D剛体レジストレーション実験を行い、投影の不整合がポーズ推定誤差にどのように伝播するかを定量化する。指標には、並進誤差（mm）、回転誤差（度）、および残留再投影誤差が含まれる。

主な貢献

1. 合成フレームワーク: 椎体透視イメージングにおける固有パラメータ校正の摂動に対する投影ドメインの感度を分析するための、制御された手法の導入。これにより、校正の影響を解剖学的運動から分離する。
2. 体系的な定量化: ランドマーク変位、輪郭変形、シルエットの重なり、および画像類似性という補完的な指標を用いた、投影整合性の包括的な評価。
3. ビュー依存の感度特性: 側方（LAT）椎体投影が、前後方向（AP）投影よりも固有パラメータ校正の変動に対して大幅に高い感度を示すことの証明。
4. レジストレーションへの定量的リンク: 固有パラメータ校正の摂動、投影ドメインの不整合、およびダウンストリームの2D–3Dレジストレーション誤差の間の定量的関係の確立。
5. 補完的評価: 投影ドメインの解析が、従来の再構成ベースの校正評価とは異なる、かつそれを補完する情報を提供することを実証。

実験結果

1. 定性的およびグローバルな整合性:

• ビュー依存性: 同等の摂動量（例：100–500ピクセルの焦点距離シフト）の下で、LAT投影はAP投影と比較して、生のDRR差分ヒートマップにおいて著しく強く、空間的に広がった不一致を示す。
• シルエットの劣化: 摂動の大きさが大きくなるにつれ、LAT投影はシルエットの整合性の顕著な低下を示す。焦点長摂動が100から500ピクセルに増加するにつれ、LATのDice係数は約0.995から約0.973に低下し、IoUは約0.991から約0.948に低下する。対照的に、AP投影は高い安定性を維持する（Dice > 0.997、IoU ~0.995）。

2. ランドマーク変位:

• シフトの大きさ: LAT投影は、実質的に大きな解剖学的ランドマーク変位を示す。例えば、特定の摂動下において、LATビューでの下終板および棘突起の先端は約15–18ピクセルのシフトを示したが、APのランドマークは2ピクセル未満であった。
• 幾何学的要因: LATビューにおけるこの高い感度は、非アイソセントリックなセットアップにより椎体が検出器に近いこと、その結果として投影スケールが大きくなることに起因する。したがって、固有パラメータの不正確さが、見かけ上のランドマーク変位を増幅させる。

3. 2D–3D レジストレーションへの影響:

• 誤差の伝播: 校正の摂動によって誘発された投影の不整合は、ダウンストリームのレジストレーション誤差へと伝播する。
• 回転の感度: 回転の整列精度は特に敏感である。直接レジストレーション実験では、焦点スケールの摂動が増加するにつれて、回転誤差は約0.5°から約3.8°へと増加した。Kabschベースのレジストレーションも、摂動に伴い回転誤差が増加した（0.05°から0.37°）。
• 並進の安定性: 並進誤差は比較的安定していたが、摂動の大きさに伴い単調に増加する傾向を示した。

意義および主張

本論文は、再投影誤差のような再構成ベースの指標のみに頼ることは、透視誘導下での椎体イメージングにおける校正品質を評価するには不十分であると主張している。本研究は以下を実証している：

• 小さな摂動の影響: 相対的に小さな固有パラメータ校正の変動であっても、椎体の投影ジオメトリ、輪郭形態、およびランドマークの位置特定に測定可能な変化をもたらす。
• ビュー依存性の重要性: イメージングシステムの感度は一様ではない。側方ビューは、奥行きの圧縮と解剖学的重なりにより、本質的に校正エラーに対して脆弱である。
• ダウンストリームへの影響: 投影ドメインの不整合は、特に回転の整列に影響を与え、ダウンストリームの2D–3Dレジストレーション精度を直接低下させる。

著者らは、本研究を、単なる体積再構成の忠実度ではなく、投影ドメインのタスク（DRR生成およびレジストレーション）の特定の要件を考慮した、より包括的な校正評価へのステップとして位置づけている。提案されたフレームワークは、幾何学的感度を特徴付けるための制御された環境を提供しており、これは将来の校正分析およびDRR生成の評価に有用であると著者らは示唆している。

限界

著者らは以下の限界を明示的に認めている：

• 合成的な性質: 実験では、散乱、検出器ノイズ、ビーム硬化、または不均一な検出器応答といった物理的効果をモデル化していない、CT由来のメッシュモデルを使用している。
• 剛体解剖学: 本研究は剛体としての椎体解剖を想定しており、軟部組織の相互作用、変形、または患者の動きは除外されている。
• 摂動の範囲: 解析は、解析的に定義された固有パラメータの摂動に限定されており、実際のシステムでは、固有、外在、および時間的な変動の複雑な組み合わせが含まれる可能性がある。
• レジストレーションの範囲: レジストレーションの評価は、制御された合成設定と定義されたランドマークに限定されており、完全な強度ベースのレジストレーション・パイプラインや、実際の臨床的な画像の劣化は考慮されていない。

今後の課題として、実世界の透視データセットを用いた知見の検証、現実的な検出器物理の組み込み、および変形可能な解剖学的シミュレーションの評価が提案されている。