技術要約：薄い幾何形状における双極子拡散誤差

問題提起

コンピュータグラフィックスにおける表面下散乱（SSS）は、主に双極子近似（およびその派生である多極子、量子化拡散など）を用いてモデル化されている。これらのモデルは、平坦で半無限の平板に対して導出された放射反射プロファイル $R_d(r)$ を仮定している。厚みのある平面的な表面には効果的であるが、この仮定は、薄い、あるいは急激に湾曲した物体（例：葉の縁、宝石のエッジ）に適用される際、系統的な幾何学的誤差を生じさせる。このような幾何形状では、光が拡散長内で第2の境界に遭遇し、双極子モデルが予測するような再放出ではなく、物体から脱出してしまう。これにより、薄い特徴部が物理的にあるべき姿よりも明るく見える「オーバーグロー（過剰な輝き）」のアーティファクトが発生する。

多極子モデルのような既存の解決策は、平面的な薄いスラブには対処できるが、一般的な湾曲幾何には対応できない。さらに、一般的な幾何学を考慮した拡散エスティメータは、内部メッシュやグリッドを必要とするため計算コストが高く、符号付き距離関数（SDF）ベースのワークフローへの統合が困難である。本論文の目的は以下の通りである：

1. 双極子モデルの薄い幾何形状における誤差を、光学的厚さのスカラ関数として特性化すること。
2. 内部の離散化を必要とせず、SDF上で直接動作する、グリッドフリーで幾何学的に正確な拡散エスティメータを提供すること。
3. レンダリング、誤差予測、および逆問題における本エスティメータの有用性を実証すること。

手法

1. スラブ誤差の理論的分析

著者らは、同一の拡散定数と境界条件の下で、標準的な双極子（半無限スラブ）と多極子（有限スラブ）を比較することで、幾何学的誤差を分離している。

• 光学的厚さの法則: 著者らは、相対的な拡散アルベドの誤差が、光学的厚さ $\tau = T/\ell_d$ （ここで $T$ は物理的厚さ、$\ell_d$ は拡散長）に対して指数関数的に減衰することを導出した。
• 閉形式のレート: 主要な誤差項は $C e^{-2\tau}$ である。ここでのレート $2$は、スラブ内を光が「往復」することに対応しており、プレファクター$C$ は材料に依存する。
• 透過率: 薄いスラブによって失われる光（双極子が考慮できていない光）は、背面を透過する。この透過率は、単一通過のレート $e^{-\tau}$ で減衰する。
• 一般化: 空間的に変化する媒体に対しては、著者らは光学的距離（アグモン距離）の概念を用いてこれを拡張している。これは、遠方の境界への最小コスト経路 $\tau^*(x)$ を表す。これにより、スラブの法則を非均質な媒体へと一般化している。

2. グリッドフリー・ジオメトリ・エキザクト・エスティメータ

任意の形状上でメッシュを作成せずに拡散方程式を評価するために、著者らはWalk on Spheres (WoS) を用いて、以下の遮蔽ポアソン方程式を解いている：
$$\nabla^2 \phi - c \phi = -f$$

• SDF表現: 物体は符号付き距離関数のみによって定義される。WoSアルゴリズムは、ドメイン内に収まる最大の球体内の点をサンプリングし、距離場を使用して球の半径と境界への近接度を決定する。
• 半空間仮定の排除: 局所的な接平面半空間を仮定する双極子とは異なり、WoSウォークは、背面やタイトな曲率を含む真の幾何形状を尊重する。
• 単一散乱: 拡散モデルは多重散乱のみを扱うため、著者らはヘニエ・グリーンスタイン位相関数から導出された、個別の閉形式の単一散乱項を追加している。これは、鋭い近接場のピークを捉えるために、短いレイマーチングを介して計算される。
• 実装: ソルバーは、CPUリファレンスと共に、プレーン・フロートのMetal（GPU）で実装されており、均質および空間的に変化する媒体（ヌル散乱ウォークを介して）の両方をサポートしている。

3. 誤差予測とハイブリッドレンダリング

形状直径関数（局所的な厚さ）を $\ell_d$ で除した実用的な予測器が導出されている。

• 補正: この予測器はアルベドの損失を推定し、標準的な双極子に補正係数 $\kappa(\tau)$ を適用する。
• ハイブリッド・ディスパッチ: 予測器は、双極子の誤差が顕著なピクセル（例：逆光の薄い領域）を特定するために使用される。これらのピクセルは高価なグリッドフリーWoSエスティメータでレンダリングされ、それ以外のピクセルは補正された双極子が使用される。

4. 逆問題（拡散光トモグラフィ）

フォワードソルバーが微分可能であるため、著者らは、拡散境界測定から隠れた吸収体（インクルージョン）の位置とコントラストを復元するという、グリッドフリーの逆問題を実証している。

主な貢献

1. スラブ法則と厚さ予測器:

   • 双極子のアルベド誤差が、平面的な場合において材料に依存しない指数レート $e^{-2\tau}$ に従うことを証明した。
   • 双極子のアルベド誤差を大幅に軽減する閉形式の補正 $\kappa(\tau)$ を提供した。
   • ブルートフォースのボリューム・パストレーサーを用いて、予測値（1.99 vs 2.0）を検証した。

2. 幾何学的に正確なグリッドフリー・エスティメータ:

   • 内部メッシュを必要とせず、SDF上の遮蔽ポアソン方程式のためのWoSソルバーを実装した。
   • 閉形式の解に対して0.75%の誤差を達成し、パストレーシングのリファレンスに対しても、双極子と比較して大幅に誤差を低減した。
   • 空間的に変化する媒体と複雑な境界条件を正常に扱うことができる。

3. 包括的な検証:

   • 背面光、前面光、サイドライトのレンズ、およびトーラスの4つのケースについて、パストレーシングのリファレンスと比較した。
   • グリッドフリー・エスティメータは、厚さに支配される背面光のケースでは改善が見られるものの、単一散乱や角度輸送が支配的な前面光や高度に湾曲したシナリオにおいて、普遍的に双極子を凌駕するわけではないことを明らかにした。

4. 微分可能な逆ソルバー:

   • 微分可能なフォワードマップを用い、物体の半径の5.4%以内の中心位置で、隠れたガウス分布のインクルージョンを高精度に復元できることを示した。

結果

• 精度: グリッドフリー・エスティメータは、閉形式のスラブ解と**0.75%**一致し、多極子解を正確に再現する。パストレーシングのリファレンスに対しては、背面光のレンズシナリオにおいて、平均生の誤差を72%から45%に低減させた。
• パフォーマンス: Apple M4 Proにおいて、GPUカーネルは毎秒 4.1 $\times$ 10$^8$ 個のボールウォークを維持する。
• ハイブリッド効率: （ハイブリッド・ディスパッチにより）フラグが立てられた背面光ピクセルに対してのみグリッドフリーソルバーを使用する手法は、純粋な双極子（0.180）や純粋なグリッドフリー（0.143）と比較して、最も低い平均トーンマップ誤差（0.132）を達成した。
• 特定された限界:
  • 曲率: 一般的な幾何形状において、誤差のプレファクターは厚さだけでなく、平均曲率（$\ell_d H$）に依存する。厚さだけでは誤差の分散の38%しか説明できず、曲率を加えることでこれが94%に上昇する。
  • モデル誤差: パストレーシングとの比較における残りの誤差は、拡散近似自体（放射輸送に対して）によるものであり、特に $\tau < 0.4$ の領域において、ソルバーの幾何学的誤差よりも支配的である。
  • コスト: 高コントラストの空間的に変化する媒体はモンテカルロノイズを増加させ、非厳密なSDF（緩い距離場）はウォークのステップ数を大幅に増加させる（$\alpha^{-2.1}$）。

意義と主張

本論文は、薄い特徴部における双極子近似に固有の特定の幾何学的誤差を定量化し、軽減することを主張している。これは、放射輸送に取って代わるものや、非常に薄い、あるいは光学的に薄い物体において単一散乱が支配的な状況を解決するためのものではない。

• 限定的な範囲: 著者らは、本手法は拡散領域内における幾何学的誤差を減少させるものであり、単一散達が支配的な非常に薄い、あるいは光学的に薄い物体において放射輸送に取って代わるものではないと明示している。
• 実用的な有用性: 導出された厚さ予測器は、計算リソースを割り当てるための有用なヒューリスティック（ハイブリッドレンダリング）として、またパストレーシングを完全に行うことなくプロダクションパイプラインでアルベドを補正するための手段として機能する。
• 基礎的な洞察: 本研究は、「双極子の失敗」がランダムなアーティファクトではなく、ランダムウォークが第2の境界に衝突する確率によって制御される、光学的距離の予測可能な関数であることを確立している。
• 逆問題の能力: グリッドフリー・ソルバーの微分可能性は、グリッドフリーの拡散光トモグラフィへの道を開いており、これはプロファイルベースのモデルでは利用不可能な能力である。

結論として、双極子は効率性の面で標準であり続けるが、その幾何学的誤差は系統的であり、光学的厚さを通じて補正可能である。また、グリッドフリーの拡散ソルバーは、特定の領域において幾何学的に正確な表面下散乱レンダリングを実現するための、メッシュフリーな代替手段となり得る。