パストレース：集約

パストレース：光伝達の統計的シミュレーションとしての物理的に正しいレンダリング

パストレースとは何ですか？

2次元の画像を、3次元シーン上で物理的に整合性のある方法で再構成することは、コンピュータグラフィックスにおける最も困難な問題の一つです。画像の各ピクセルに到達する光を正確に計算するには、物体の位置だけでなく、表面素材の光学特性、媒質の散乱挙動、シーン内の幾何学的関係をすべて考慮する必要があります。この文脈においては、光の伝達を統計的にモデル化する汎用のレンダリング技法であり、現代の物理ベースレンダリングの基盤をなすものの一つとして広く認識されています。

パストレーシングのアプローチは、シーンに送信されたフォトンやレイのランダムな経路を通じてエネルギー伝達をシミュレートします。これにより、シャドウ、アンビエント・オクルージョン（環境遮蔽）、カラー・ブリーディング（色の移り）、ディプス・フィールド効果、そして完璧な反射といった複雑な光学現象が、物理法則に従った形で表現されます。本記事では、パストレーシングの科学的基礎、数学的メカニズム、計算効率を提供する高速化構造、そして産業応用について検討します。

科学的基礎と歴史的発展

レンダリング方程式（Rendering Equation）

パストレースの理論的基盤は、1986年にジェームズ・カジアによって定義されたです。この方程式は、ある表面の一点から放出される総光量を、その点に入射する直接光源からの寄与と、その他の全ての表面からの反射による間接的（多重散乱）光との相互作用として表現します。

$$L_o(mathbf{x}, omega_o) = L_e(mathbf{x}, omega_o) + int_{Omega} f_r(mathbf{x}, omega_i, omega_o), L_i(mathbf{x}, omega_i), (omega_i cdot mathbf{n}), domega_i$$

この方程式はレンダリング方程式（Rendering Equation）と呼ばれ、1986年にケビン・トーマスが定式化した、コンピュータグラフィックスにおける光の伝達を記述する基礎的な方程式です。以下に各項を説明します。

## 各項の意味

| 記号 | 意味 |
|——|——|
| $L_o(mathbf{x}, omega_o)$ | 地点 $mathbf{x}$ から方向 $omega_o$ に向かう出射輝度（発光を除く） |
| $L_e(mathbf{x}, omega_o)$ | 地点 $mathbf{x}$ からの直接発光輝度 |
| $f_r(mathbf{x}, omega_i, omega_o)$ | BRDF（双反射分布関数） |
| $L_i(mathbf{x}, omega_i)$ | 地点 $mathbf{x}$ へ方向 $omega_i$ から入射する入射輝度 |
| $mathbf{n}$ | 地点 $mathbf{x}$ における法線ベクトル |
| $Omega$ | 半球（$mathbf{n}$ を軸とする）|
| $omega_i cdot mathbf{n}$ | コサイン項（ランベルトの余弦則）|

## 方程式の解釈

この方程式は、ある地点から観測される全出射輝度が、2つの成分の和であることを表しています。

### 1. 直接発光成分 $L_e$
物体自身が光を放出している場合の輝度です。光源（電球、ディスプレイ、炎など）で値が正になります。

### 2. 散乱成分（積分項）
他の場所から届いた光が物体表面で散乱され、観測方向へ向かう輝度です。

## 積分項の詳細

$$int_{Omega} f_r(mathbf{x}, omega_i, omega_o), L_i(mathbf{x}, omega_i), (omega_i cdot mathbf{n}), domega_i$$

この積分は、全半球方向からやってくる光を統合していることを意味します。

– $L_i(mathbf{x}, omega_i)$：各方向 $omega_i$ から入射する光の強さ
– $f_r(mathbf{x}, omega_i, omega_o)$：その入射光が、観測方向 $omega_o$ へどれだけ散乱されるか（材質の性質）
– $(omega_i cdot mathbf{n})$：光が斜めから当たるほど効率が下がることを表すコサイン係数
– $domega_i$：微小立体角要素

## 重要な性質：再帰的構造

レンダリング方程式の最も特徴的な点は、積分項の中に $L_i$ が含まれ、それが再びレンダリング方程式に従うことです。

$$L_i(mathbf{x}, omega_i) = L_o(mathbf{x}, omega_i)$$

つまり、入射輝度は、その方向から見て物体が放出・散乱した出射輝度と等しくなります。これにより、

– 光が複数の物体間で多重反射する効果
– 間接照明（間接的に照らされた表面から放たれる光）

などが自然に表現できます。

## 視覚的な理解

$$underbrace{L_o}_{text{観測される輝度}} = underbrace{L_e}_{text{直接発光}} + underbrace{int_{Omega} underbrace{f_r}_{text{材質}} cdot underbrace{L_i}_{text{入射光}} cdot underbrace{(omega_i cdot mathbf{n})}_{text{角度}} , domega_i}_{text{半球全体からの散乱光}}$$

この方程式は、リアルなレンダリングを実現するための多くのアルゴリズム（レイトレース、モンテカルロ積分、輻射伝達など）の理論的基盤となっています。

この式で$L_o$は出力光、$L_e$は表面が自ら放出する光（発光）、$f_r$は（BRDF）として知られ、素材の光に対する反応を定義する関数です。積分項は、全方向からやってくる光（$L_i$）がどのように散乱されるかを計算します。パストレーシングは、レンダー式がシーン内の無限に続く反射の鎖を一つの恒等式に要約しているため、その数学的な核心と呼ばれています。

モンテカルロ積分とランダムサンプリング

レンダー方程式における積分は解析的に求解するのが複雑すぎるため、パストレーシングはを用いてその値を近似します。この方法は、ランダムにサンプリングされた方向に対して被積分関数の平均を取ることで動作し、積分の次元が大きければ大きいほど効率的になります。これはまさに、レンダー方程式に見られるような高次元問題に理想的です。サンプル数が増えるにつれ解はしますが、同時に画像に統計的なノイズが発生します。そのため、パストレーシングは実務的に数千から数万のサンプリングを必要とする技術であり、高品質な結果を得るには長い計算時間を要します。

メカニズム — 追跡ルートはどのように機能するのか？

基本ループ：トレースの生成と拡散

パストレースのアルゴリズムの基本ループは次のように動作します。画像の各ピクセルに対して、カメラからを生成します。この光線がシーンジオメトリと交差すると、その点でBRDFに従ってランダムに新しい方向が選択され、追跡用シーン内に「拡散」します。このプロセスは衝突のたびに繰り返され、光源に到達するか、特定の深さの境界に達するまで続きます。このモデルにより、多重反射が自然に計算されます。これは、従来のラスタライゼーション手法では手動でモデル化する必要があった効果を、ここでは自発的に導出することを意味します。

サンプリング戦略と統計的安定性

ノイズ低減のためには、さまざまなサンプリング手法が開発されてきました。層別サンプリング（stratified sampling）、インポータンスサンプリング、ブルーノイズマスク、テンポラル集積（temporal accumulation）などがその一例です。1997年にエリック・ヴィーチが発表した博士論文『Robust Monte Carlo Methods for 浅色 Transport』は、パストレーシングにおけるサンプリングの安定性に関する基礎的な参照文献として広く認識されています。特に「ラッソ・ルーレット（Russian roulette）」と呼ばれるランダムな打ち切りメカニズムは、遠くに位置したり寄与が小さいたりするパスが、計算を無駄に延長してしまうのを防ぎます。

高速化構造：交差テストの最適化

2次元の画像を生成するためには、シーン内の各光線が数千の物体と交差する必要があります。これにより直接計算は（O(n)）になります。パストレースは実用的に利用可能ですが、交差テストのコストはしばしばレンダリング時間の大部分を占めます。そのため、レイトレーサーの中核をなす要素であるは、不要な交差テストを削減するように設計されています。光線をシーン内で追跡するとき、その瞬間に出会うすべての物体と個別に試す必要がありますが、ほとんどのシーンでは光線はこれらの物体のごく一部に接近することさえないのです。

2つの基本的なアプローチがあります。は、3次元空間を軸に整列した箱や領域に分割し、どのオブジェクトがどの領域と重なるかを記録します。交差が見つかったとき、光が通る領域に属するオブジェクトのみが検査されます。このアプローチの例がという高速化構造です。この構造は空間を適応的に分割し、各領域に属するオブジェクトの数に応じて深くなっていきます。

一方、もう一つのアプローチではが、シーン内のオブジェクトを近くの小さな要素のグループに分割します。部屋モデルは、まず4つの壁、天井、椅子として分割されます。レイが部屋の全体を囲む境界体積を見逃した場合、その中の全てのオブジェクトが一発で除去（culling）されます。そうでない場合は、各部分に対して順次テストが行われます。椅子の境界体積に命中した場合、今度は脚、座面、背もたれといった要素が個別に試されます。このアプローチの現代的な対応物が構造であり、パストレース応用で広く使用されています。エリア分割もオブジェクトベースの分割も、交差コストを削減する上で成功しており、どちらかが他方より優先されるべきという根本的な理由はありません。

セクター別・応用における次元：ゲームフィルムとリアルタイム・ライティング

シネマ・ビジュアルエフェクトにおけるパストレース

映画産業において、パストレーシングは長年にわたりフォトリアルなアニメーション映画の標準的な技術として定着してきました。レンダリング時間に許容性がこの分野では、光の伝達、拡散素材における色の変化、大気散乱などの効果が物理的に正確にモデル化されており、パストレーシングはシネマティックレンダリングの品質を決定する主要な要因の一つとしてその地位を維持してきました。

リアルタイムレイトレーシング：ゲーム産業に変革をもたらす技術

ゲーム産業では、リアルタイムパフォーマンスへの要件から従来ラスタライゼーションが採用されてきましたが、近年はハードウェアレベルでのレイトレーシングユニット（RTコア）の開発に伴い、がますます普及しています。現代のゲームエンジンは、影・反射・環境散乱をリアルタイムで計算することができ、これによりプレイヤーにはシネマティックな品質ながらインタラクションに対応する速度での体験が提供されています。この移行には、ハードウェアアーキテクチャと、サンプリング・ノイズ低減ソフトウェア（テンポラル・アンチエイリアシングやデノイジング技術など）の双方の協調的な発展が必要でした。なぜなら、パストレーシング由来の統計的ノイズは、リアルタイム応用において許容可能な水準以下に低減しなければならないからです。

学術的な結果

パストレースは、光の伝播を物理法則の枠組み内でモデル化する際の最も一般的で、数学的に洗練された解法である。カジアのレンダリング方程式からヴィーチの統計的安定性に関する研究に至る理論的基盤；モンテカルロ積分による高次元問題解決能力；高速化構造物によって実用的な解となり、ハードウェアの発展によりリアルタイム応用へと持ち込まれたこと。これらの要因が、パストレースをコンピュータグラフィックスの最重要研究分野の一つへと押し上げた。ゲームから映画まで幅広い分野でこの技術が採用されていることは、物理ベースのレンダリングが学術的な話題であるだけでなく、産業標準として定着した基盤技術であることを示している。

学術的評価と参考文献

パストレースは、文献では「経路積分（path integral）」としても知られる光伝達方程式の統計解に基づいて構築されています。このアプローチの強みは、光学現象を何らかの前仮定なしに自発的に導出できる点にありますが、弱みは高いサンプリング数に依存する計算コストと、それに伴うノイズです。このため、最新の研究は主にサンプリング効率の向上、ノイズ低減、ハードウェア 친화的な高速化構造の設計に焦点を当てています。

参考文献:
– Kajiya, J. T. (1986). . SIGGRAPH ’86 Proceedings。— レンダー方程式の定義とパストレーシングの理論的基盤の確立。
– Pharr, M., & Humphreys, G. (2010). （第3版）。Morgan Kaufmann。— パストレーシングのメカニズム、高速化構造、実装詳細に関する包括的な参考文献。
– Veach, E. (1997). 。Stanford University doctoral thesis。— サンプリングの安定性、インポータンスサンプリング、ルーシアンルーレットなどの統計的手法の基本的な参考文献。
– Jensen, H. W., & Majoring, P. (2000). 。A K Peters。— パストレーシングの実践的応用とリアルタイムの照明モデリングに関する包括的な研究。

ソースノート：本研究は、https://www.pbr-book.org/4ed/Primitives_and_Intersection_Acceleration/Aggregates のアーカイブされたドキュメントに基づき、学術的な技術レビュー形式で編纂・更新されました。