ボリュームレンダリング入門(非一様媒質編) / Introduction to VPT (heterogeneous)

OSKのLT会で発表した際の資料。
非一様媒質でVolumtric Path Tracingを行う方法について説明しました。

Transcript

3. ! () " () 散乱係数[1/m] 吸収係数[1/m] 散乱しやすさ 吸収しやすさ 一様媒質と違う点 これらのパラメータが空間的に変化する (xに依存して変化する)
4. 非一様媒質であることによる弊害 , = ! ∫! " ## $!%& '% 透過率の解析的な評価が難しい。 距離サンプリングに必要な逆関数法も使えない。 ! ≔ " + #()
5. いくつかの解決策 ・デルタトラッキング ・レイマーチング ここではデルタトラッキングについて解説していく
6. デルタトラッキング & ∇ , = −! , + " # , + / $! % & , ' , ' ' まず、媒質中での放射輝度の変化を記述した放射伝達方程式(RTE)は 吸収・散乱による放射輝度の減少 発光による放射輝度の増加 散乱による放射輝度の増加 通常のボリュームレンダリング方程式はこの微分方程式を解くことで得られる。 そこで、この式の持つ意味を変えないで扱いやすい形に変形していくことを考える。 ※微分方程式の解き方はmemoRANDOMを参照してください。
7. デルタトラッキング そこで、衝突してきたレイの方向を変えず散乱させるような擬似的な粒子を導入する。 この衝突をnull-collisionと呼ぶ。
8. デルタトラッキング null-collisionによる放射輝度の変化は次のようになる −( , + ( / $! − ' , ' ' = 0 ! ()は" ()や# ()の役割と同様に、上の擬似的な散乱の起こりやすさを表す係数である。 減少 増加
9. デルタトラッキング −( , + ( / $! − ' , ' ' = 0 先ほど出てきた次の式 これを放射伝達方程式に加えても意味は変わらないことは明らかですね & ∇ , = −! , + " # , + / $! % & , ' , ' ' − ( , + ( ∫ $! − ' , ' '
10. デルタトラッキング 5 ! ≔ !() + (()とすれば、次のように書ける & ∇ , = −5 ! , + # $ , + + %! " & , ' , ' ' + ( + %! − ' , ' '
11. デルタトラッキング この微分方程式を解けば、null-collisionを考慮したボリュームレンダリング方程式を得ます。 , = 1 $ % 2 , 2 & # 2 & ' , + " 2 & 1 (! ) , * , * * + ! 2 & (, ) + 2 , (, ) 5 , ≔ : ∫ " # < =$ >:%? @%
12. デルタトラッキング ここで、 2 & ()の定義に戻れば、 ! ! ≔ ! () + " () ! ()の定め方は自由なので、 2 & が定数Cになるように! ()を定めれば、 ! , ≔ 0 ∫! " 2 3# 045 64 = 0 ∫! " 7 64= 078 透過率が解析的に求まっている！！(null-collisionを経路に組み込んだ上での透過率であることに注意) ちなみに定数Cは& ()の最大値に設定することが多い
13. デルタトラッキング また距離サンプリングに使う確率変数の分布関数は、 = 1 − ! (, ) = 1 − #$% これは逆関数法が適用でき、 = −log(1 − u) 非一様媒質を一様媒質のように扱うことができるようになった
14. デルタトラッキング 以上の内容から、nullな粒子が一様でない媒質を一様化させていることがわかる
15. デルタトラッキング モンテカルロ積分の推定値は次のようになる。 , = 2 +#, ≥ +#, , = 2 +#, 2 +#, , = , (" ≥ +#, ) , (" < +#, ) " ≔ 距離サンプリングで求めた値 , ≔ - . &" / 0# , 0$ , / 0# , ) &" 0$ , / 0# , 1% ,,3 4 - . &" / 0# , 0$ , / 0# , - . &" / 0# , 0$ , / 0# , 1% ,,3 41% ,,3 5&(3,3')- . &" / 0# , 0" , / 0# , )(3')) &" 0" , / 0# , 1 ,,3' 4 5&(3,3')- . &" / 0# , 0" , / 0# , 5&(3,3')- . &" / 0# , 0" , / 0# , 1 ,,3' 41 ,,3' 8(393)- . &" / 0# , 0" , / 0# , 8(393)) &" 0" , / 0# , 1 ,,3 4 - . &" / 0# , 0" , / 0# , - . &" / 0# , 0" , / 0# , 1 ,,3 41(,,3) ≔ 境界上の点 (吸収) (散乱) (null-collision)
16. 実際のレイトレースの流れ
17. レイがボリュームの境界である直方体にヒットする ※なんで直方体なのかは後で説明します！
18. 境界情報を取得する : $
19. 距離サンプリング 4 = 9 − log(1 − ) $ :
20. 距離サンプリング 4 ≥ : のときは、 そのまま通過 $ : "
21. $ : " 現在の点での、密度データ(%, ", ( )の取得
22. $ : " 散乱・吸収・null-collisionのどれが起きるかを確率的に決める =% < =$ + =& < =$ + =' < =$ = 1 なので、 確率 =% < =$ で、吸収・発光 確率 =& < =$ で、散乱 確率 =' < =$ で、null-collision
23. $ : " 吸収イベントが起こるとき 発光のある媒質(火、爆発など)なら、 return Throughput × ' (, ) 発光のない物質なら、 return Throughput × 0
24. $ : " 散乱イベントが起こるとき 散乱方向をサンプリングして、 そして、ループの最初に戻る 次に飛ばすレイを設定する。
25. $ : " null-collisionイベントが起こるとき レイの方向を変えず、 レイの始点だけを現在の点に変えて、 ループの最初に戻る。
26. 格子点に、密度データを格納する。ボリューム境界が直方体であることの理由はこれ。 (実際は3Dだが簡単のために2Dで書いている) データの表現方法
27. データの表現方法 の順番で線形補間していく 任意の点での密度データを取得するときは、補間を用いる。 上の図はトリリニア補間。
28. 以上のことを実装すると
33. NanoVDB
34. ・CUDAでVDBを扱うためのNVIDIA謹製のライブラリ ・外部ライブラリへの依存が少ない！ ・ヘッダーオンリー
35. ・vdbファイルを変換したnvdbファイルを使う。 使い方 ・まずはグリッドデータの読み込み
36. 使い方 ・GPU側へのグリッドデータのアロケーション
37. 使い方 デバイスメモリ上のグリッドデータは、 ホストメモリ上のグリッドデータは、
38. 使い方 ・NanoVDBはボリュームの境界であるAABBとレイの交差判定もやってくれる。
39. 使い方 ・密度のデータを取得する際のトリリニア補間はこんな感じ
40. このような感じで、ボリュームレンダリングに必要なものは大体揃ってます！！ みんなもNanoVDBでボリュームレンダリングに挑戦してみよう！！