RespGeomLib: A Reproducible Parametric Engine for Generating Analysis-Ready Human Airway Lumen Geometry

技術要約: RespGeomLib

問題提起
ヒトの気道幾何学構造を生成する現在の手法は、そのソースに応じて重大な限界に直面している。CT由来のモデルは、解剖学的な根拠はあるものの、遠位部のスキャン解像度に制約され、画像のノイズや部分容積効果の影響を受けやすく、形態計測分析や数値流体力学（CFD）に適した状態にするには、分岐部付近での大規模な手動クリーニングが必要となることが多い。一方、手続き型（プロシージャル）の代替手法は再現性を提供するが、管状のプリミティブを継ぎ合わせる手法に依存することが多い。このアプローチは、非平滑な接合部、鋭い継ぎ目、および不明瞭な開口境界を導入し、ダウンストリームの測定やシミュレーションの品質を低下させる。さらに、全領域に対してインプリシット（陰関数）表面モデリングを適用することは、近位から遠位の枝にかけての空間的な広がりとスケールの大きな変化により、計算コストが極めて高くなる場合が多い。

手法
著者らは、コンパクトで人間が読み取り可能なYAML仕様から、解析準備が整ったヒト気道ルーメン表面を生成する、再現可能なパラメトリックエンジンであるRespGeomLibを提示する。本フレームワークは、幾何学的な柔軟性と計算効率のバランスを取るために、ハイブリッド抽出戦略を採用している：

1. パラメトリック・プリミティブとアセンブリ: システムは、3つの幾何学的プリミティブ（Pipe（テーパー付き管）、Y2（2分岐接合部）、Y3（3分岐接合部））を利用する。これらは局所的な標準フレーム内で定義され、明示的な入口/出口ポートを持つ。アセンブリは、子プリミティブが仰角・方位角によって定義される剛体変換（回転および平行移動）を用いて親ポートに接続される、ポートベースのインターフェースを介して行われる。
2. 滑らかな接合部のブレンディング: 継ぎ目の問題を回避するため、接合部はインプリシットなスムース最小化（smooth-min）定式化を用いてモデル化される。幾何学は、以下のブレンディングされたスカラー場のゼロ・レベルセットとして定義される：
   $$\phi(x) = -\frac{1}{\kappa} \log \left( \sum_{i=1}^{N} e^{-\kappa \phi_i(x)} \right)$$
   ここで、$\phi_i(x)$ は第$i$番目の枝の中心線までの距離を表し、$\kappa$ はブレンディングの鋭さを制御する。これにより、内部壁のない、滑らかで連続的なルーメン接合部が保証される。
3. 局所インプリシット抽出: 全体ツリーのボクセル化に伴うコストを避けるため、エンジンは接合部の周囲のタイトなグリッド上でのみ、インプリシット等値面抽出（マーチングキューブ法による）を局所的に適用する。直線およびテーパー付きセグメントは解析的に生成される。局所的なグリッド間隔は、局所的な枝の半径に対して相対的に適応的に定義される。
4. 制御されたバリアントと成長: 本フレームワークは、形態計測に基づいたツリー生成（例：Weibel/ICRPモデル）および合成バリアントの作成をサポートする。局所的なルーメンの摂動（狭窄または拡張）は、パラメータ化された半径変調関数 $r(s) = r_0(s)(1 + \alpha w(s))$ を通じて適用され、基礎となるトポロジーとポートの定義を保持する。
5. 後処理: パイプラインは、ポート平面での表面クリッピングによる開口ポート境界の強制、それに続く決定論的なクリーニング（重複/退化の除去）および多様体性の検証を行い、CFDへの適合性を確保する。

主な貢献

• RespGeomLibエンジン: コンパクトな仕様から、明示的な入口/出口ポートを持つ三角形ルーメンメッシュを生成する、YAMLベースのパラメトリックシステム。
• 滑らかな接合部の構築: シームレスな分岐および三分岐を生成し、継ぎ目の継ぎ目や内部壁を排除する、スムース最小化インプリシット定式化。
• 効率的なハイブリッドメッシング: 全体的なボクセル化を回避する局所インプリシット/解析的パイプラインを実現し、グローバルなインプリシット抽出よりも大幅に高速かつメモリ効率が高いことを実証。
• 再利用可能なダウンストリームの有用性: 制御された合成バリアントの生成、形態計測に基づいたツリー合成、およびCFDシミュレーションへの正常なエクスポートのデモンストレーション。

結果

• 接合部の品質: Boolean/スティッチ（継ぎ合わせ）のベースラインと比較して、RespGeomLibは0個の非多様体エッジ（ベースラインは6個）を生成し、ほぼ完璧なポート平面性（最大誤差0.0010 vs 5.9796）を実現した。ベースラインは平面キャッピングテストに失敗したが、RespGeomLibは水密性の高い表面の生成に成功した。
• 計算のスケーリング: バランスの取れたY2ツリーの同深度比較において、局所抽出戦略は、深さ1〜3においてグローバルなインプリシット抽出よりも76.7倍から279.8倍高速であった。また、ピークメモリ使用量も1.11倍から1.31倍削減された。グローバルなベースラインは、浅い深度であっても計算コストが高くなっていた。
• CFDの適合性: RespGeomLibのY字接合部における定常非圧縮RANSシミュレーションは、安定したソルバー挙動を示した。このケースは、物理的に妥当な流れの再分配（46.5%/53.5%の非対称分割）、約 $1.0 \times 10^{-1}$ Paの圧力降下、および $10^{-5}$ オーダーの最終残差をもたらした。
• 合成バリアント: システムは、固定されたベースジオメトリから、一貫したトポロジーとポート定義を維持しながら、ラベル付けされた狭窄および拡張気道のコホートを正常に生成した。

意義と主張
本論文は、RespGeomLibを患者固有の再構成ツールとしてではなく、解析準備が整ったルーメン表面を生成するための再現可能な幾何学エンジンとして位置付けている。その主な意義は、解剖学的再構成と手続き型モデリングの間のギャップを埋える、コンパクトで編集可能かつシミュレーション準備が整った表現を提供することにある。

著者らは、ハイブリッドアプローチ（解析的セグメント + 局所インプリシットブレンディング）が、グローバルなインプリシットモデリングの計算コストに対する実行可能な解決策を提供すると同時に、スティッチされたプリミティブの幾何学的欠陥を克服することを主張している。これにより以下が可能となる：

1. 再現可能な形態計測: 標準的な形態計測ルールに基づく、一貫した気道ツリーの生成。
2. 制御されたベンチマーク: 特定のラベル付けされた幾何学的摂動（例：狭窄）を持つ合成バリアントの生成による、感度分析。
3. シミュレーションへの適合性: 即座にCFDワークフローで使用できる、水密で多様体な、明確な境界ポートを持つメッシュの生成。

本研究は、理想化されたルーメン幾何学に焦点を当てており、患者固有の非対称性、非円形断面、または壁厚についてはまだモデル化していないという限界を認めている。CFDの結果は、包括的な流体力学ベンチマークとしてではなく、ユーザビリティのデモンストレーションとして提示されている。