Surface Extraction for Industrial CT Based On Adaptive Particle System

技術要約：適応型粒子システムに基づく産業用CTの表面抽出

問題提起
産業用コンピュータ断層撮影（ICT）は、複雑なワークピースの内部および外部の幾レトリカル情報を取得するための重要な非破壊検査技術である。しかし、ICTデータから高精度かつ高品質な表面メッシュを抽出することは依然として困難である。従来の静的な抽出アルゴリズム（マーチングキューブ法やデュアルマーチングキューブ法など）は、しばしば位相的な曖昧さ、メッシュ品質の低下（細長い三角形など）、および不均一なサンプリング分布に悩まされる。さらに、これらの手法は、閾値選択や、ノイズやビーム硬化といったCTデータ特有のアーティファクトに対して非常に敏感である。粒子システムに基づく動的な手法は、より優れた分布を実現する可能性があるが、広範な粒子間力の計算による計算効率の低さ、サンプリングレート増加時におけるシステム平衡の維持の困難さ、および精度とデータボリュームのバランスを取るための適応能力の不足といった課題に直面している。

手法
著者らは、マルチスケール・ボクセルピラミッド・フレームワークと統合された、粒子システムに基づく適応型表面サンプリングアルゴリズムを提案している。本手法は以下の段階を経て進行する。

1. マルチスケール適応型精緻化:

   • エッジ保存型ガイダンスフィルタリングとダウンサンプリングを用いて、生のCTデータからボクセルピラミッドを構築する。
   • 離散グリッドにおける方向バイアスを軽減するため、各スケールにおいて3D Scharrオペレータを用いて勾配および曲率場を計算する。
   • コース・トゥ・ファイン（粗から密へ）の適応型精緻化戦略を採用する。高曲率領域を統計的閾値（平均 + $k \cdot \sigma$）を用いて特定し、形態学的膨張によって拡張することで「候補領域」を形成する。これらの領域はより細かいスケールへと精緻化され、これにより、幾何学的に複雑な領域ではサンプリング密度が高まり、平坦な領域では疎な状態に保たれる。

2. デカップルされた力を備えた拡張粒子システム:

   • 初期化: 粒子は、局所解を避けるためのランダムな摂動を伴うコーススケールの応答ユニット内に初期化される。
   • 勾配加重投影（引力）: 単純な等値面への吸引ではなく、本アルゴリズムは、グレースケール情報と二次勾配（ラプラシアン）のゼロ交差を組み合わせた、デカップルされた適応型加重投影を使用する。
     • グレースケール項は、表面近傍への迅速な収束を駆動する。
     • 二次勾配項は、位置を精密なエッジ位置へと精緻化する（サブピクセル精度）。
     • 重みは動的に調整される：初期反復ではグレースケールが支配的であり、後半では二次勾配が微調整を担う。
   • 局所短距離斥力: 粒子の分布を最適化するため、斥力は局所的な短距離相互作用として扱われる。粒子は、スケール依存のカットオフ半径によって定義される26近傍内でのみ相互作用する。
   • 接線成分の抽出: 斥力を表面の接平面上に投影する。分布の均一性を最適化するために、接線成分のみが適用される。これにより、斥力が粒子を等値面から押し出し、投影精度を損なうことを防ぐ。

3. 増分粒子挿入:

   • 従来の粒子分割メカニズムによって引き起こされる不安定性に、論文ではボロノイ図に基づく増分挿入戦略を導入している。
   • システムが平衡状態に近づくと、高曲率領域の接平面上に2D局所ボロノイ図を構築し、最大の未カバー空隙を特定する。
   • 新しい粒子は、これらの空隙の重心に挿入され、3D空間へとマッピングされる。このメカニズムにより、システムのエネルギ変動を最小限に抑え、密度増加時における平衡を維持する。

主な貢献

• マルチスケール曲率適応型サンプリング: ボクセルピラミッドと曲率ベースの精緻化の統合により、アルゴリズムはサンプリング密度を自動的に調整でき、複雑な領域の詳細を保持しつつ、平坦な領域でのデータボリュームを削減できる。
• デカップルされた力メカニズム: 時間領域において引力（投影）と斥力（分布）を分離し、グレースケールと二次勾配の加重投影を利用することで、分布の均一性を犠牲にすることなく、精密なエッジ位置特定を実現する。
• 効率的な局所相互作用: 斥力を26近傍内の局所的な短距離力として扱うことで、グローバルな相互作用モデルと比較して計算効率を大幅に向上させている。
• 安定した増分挿入: ボロノイ誘導型の挿入メカニズムにより、従来の粒子システムにおける失敗点である、システムの機械的平衡を乱すことなく適応的なサンプリングレートの増加が可能となる。

実験結果
本手法は、2つのアルミニウムワークピースの産業用CTデータを用いて検証され、その結果はマーチングキューブ法（MC）、デュアルマーチングキューブ法（DMC）、および古典的な粒子システムと比較された。高精度な三次元測定機（CMM）が真値として使用された。

• サンプリング効率: 提案手法は、再構成精度を維持しながら、MCおよびDMCと比較してサンプリングポイントの数を約70%削減し、ストレージおよび処理の負担を大幅に軽減した。
• 分布品質: サンプリング分布の統計的分析（六角形タイリングの類似性により評価）により、提案手法はMCおよびDMCよりも有意に高い均一性と、より少ないクラスタリングアーティファクトを実現したことが示された。
• 寸法精度: アルゴリズムは、二次勾配補正により、ほとんどの形状において、より優れた寸法精度（円筒径測定）を示した。
• 幾何学的精度: 特定の形状に対して、一部の従来手法と比較して幾何学的精度（円筒度や平面度）においてわずかなトレードオフが見られたものの、結果は計測用途として許容範囲内に留まった。マルチスケールアプローチは、ノイズ感受性の緩和に寄与した。
• 計算コスト: 反復的な最適化を行うため、実行時間は静的な手法（MC/DMC）よりも長いが、局所的な力計算と適応型サンプリングにより、従来の粒子システム手法よりも大幅に高速である。

意義および主張
本論文は、提案手法が産業用計測CTアプリケーションにおいて優れた工学的実用価値を提供することを主張している。適応型サンプリングによって、より少ないポイント数で実現することで、サンプリング精度とデータボリュームのトレードオフを成功裏にバランスさせている。本手法は、ノイズやアーティファクトに対する堅牢性を大幅に向上させ、従来の静的および動的な表面抽出アルゴリズムと比較して、より均一な粒子分布を保証する。

著者らは以下の限界も認めている：出力は現在ポイントクラウドであり、トポロジー的な接続性を欠いているため、メッシュ生成のための後処理が必要である。また、マルチスケール・ダウンサンプリングへの依存は、ピラミッドのレベル数が過剰な場合に情報の損失を招く可能性があり、マルチスケール場のデータを保持するためメモリ要件が高い。今後の課題として、初期化時におけるトポロジー構造の保持を組み込むこと、および有効な領域のみを保持することによるメモリ使用量の最適化が挙げられている。