The Textbook of RayTracing @TDU

[—ATOC—] 1 [—AUTO_SECTION_NUMBER—] 1

光の屈折

水やグラス，ビー玉など，光が透過する物質では別の物体との境界(媒質の境界)において光が直進せず屈折する(図1)．

完全鏡面反射と同様に，光が透過する物体を見ているとき，観察される色は物体そのものの色ではなく 物体を通過して届いた別の物体の色である(図2)．

図2. 透過する物体を通過して届いた別の物体の色を見ている

これをレイトレーシング法で再現するには，完全鏡面反射を実装したときと同様に， 物体表面の色を決める処理の中で，視線と物体の交点を仮想的な視点としてレイトレーシングを行えばよい(図3)．

図3. 屈折の再現

屈折方向および光の反射量・透過量には法則が存在する．次節ではそれらについて解説する．

光の屈折に関する物理法則

透過と完全鏡面反射の関係

光が透過する物体では，必ず完全鏡面反射も起こる．前章で扱った完全鏡面反射は物体が光を透過しない 透過現象の特別な場合である．

つまり，光が透過する物体の表面に光が入射すると，光の一部が物体内部に透過し，残りが反射する(図4)．

図4. 透過と反射

反射する場合は光の入射角の正反射方向に反射するのであるが，透過する場合には物質の屈折率に依存して屈折方向が決まる． この関係を表すのがスネルの法則である．

また，どれだけの割合の光が透過し，どれだけの割合の光が反射されるのかは物体の屈折率，入射角，屈折角に応じて決まる．この関係を表すのがフレネルの公式である．

スネルの法則

スネルの法則は，光が屈折する際の屈折方向に関する法則である．

絶対屈折率$$\eta_{1}$$の物質1から，絶対屈折率$$\eta_{2}$$の物質へ光が入射するとき， 入射角$$\theta_{1}$$と屈折角$$\theta_{2}$$には以下の関係が成立する(図5)．

図5. 入射角と屈折角

実際に屈折角を計算する場合には，入射ベクトル$$\vec{\bf i}$$，法線ベクトル$$\vec{\bf n}$$と物質1，2の 絶対屈折率$$\eta_{1},\eta_{2}$$から，屈折ベクトル$$\vec{\bf f}$$を求めることになる(図6)．

屈折ベクトル$$\vec{\bf f}$$は以下のように計算できる．

導出過程を以下に示す．

フレネルの公式

フレネルの公式は，反射と屈折の割合を表す式である．

物質の絶対屈折率が$$\eta_{1},\eta_{2}$$，入射角が$$\theta_{1}$$，屈折角が$$\theta_{2}$$である時， 光の完全鏡面反射率$$c_{r}$$と透過率$$c_{t}$$は以下のようになる．

$$\rho_{\parallel}$$，$$\rho_{\perp}$$は，それぞれp偏光反射率，s偏光反射率で，以下のように定義される．

(p偏光s偏光とは光が物質の境界面に進入する際の偏光成分であるが，実装の際は単に上記の式を使えばよいため，ここではこれ以上の解説はしない． 気になる者は調べてみると良い．)

反射光の計算式

以上の屈折に関する物理法則から反射光の計算式は以下のようになる．

屈折光を計算するための仮想的な視点位置と視線方向は以下のようになる．

• 仮想的な視点の位置
  • 物体と(大元の)視線との交点$$\vec{\bf p_{\cal i}}$$
• 仮想的な視点の視線方向
  • スネルの法則に基づいた屈折方向$$\vec{\bf f_{\cal e}}$$
  • $$\vec{\bf f_{\cal e}} = \frac{\eta_{1}}{\eta_{2}}\left(\vec{\bf d_{\cal e}}-\left(\sqrt{\left(\frac{\eta_{2}}{\eta_{1}}\right)^{2}-\left(1-cos^{2}\theta_{1}\right)}-cos\theta_{1}\right)\vec{\bf n}\right)$$
    • $$\vec{\bf d_{\cal e}}$$：（大元の）視線ベクトル（正規化する必要がある）
    • $$\vec{\bf n}$$：面の法線ベクトル

反射光の計算式は完全鏡面反射の場合と同じである．

• $$L_{a}$$：環境光の反射光の輝度
• $$L^{i}_{d}$$：i番目の光源からの直接光の拡散反射光の輝度
• $$L^{i}_{s}$$：i番目の光源からの直接光の鏡面反射光の輝度
• $$L_{f}$$：完全鏡面反射光/屈折光の輝度

完全鏡面反射光/屈折光の輝度$$L_{f}$$を以下のように変更する．

• $$k_{f}$$：完全鏡面反射光/屈折光係数
• $$c_{r}$$：完全鏡面反射率
• $$c_{t}$$：透過率
• $$L_{re}$$：視線ベクトルの正反射方向$$\vec{\bf r_{\cal e}}$$方向で得られる輝度
• $$L_{fe}$$：$$\vec{\bf f_{\cal e}}$$方向で得られる輝度

光の屈折の実装

仮想的な視点の位置と方向

7-2.2. シャドウレイと同様に始点位置に注意する．
（大元の）視線ベクトルと物体の交点をそのままレイの始点として使用するとうまくいかない(説明省略)．

material構造体/scene構造体の変更

前節の屈折光の対応のため，material構造体にさらなる変更を加える．

typedef enum 
{ 
  MT_DEFAULT,       /* 通常の質感 */
  MT_PERFECT_REF,   /* 完全鏡面反射を使用する */
  MT_REFRACTION,    /* 完全鏡面反射/屈折を使用する */
} material_type;

typedef struct
{
  colorf_t ambient_ref;   /* 環境光反射率(RGB) */
  colorf_t diffuse_ref;   /* 拡散反射率(RGB) */
  colorf_t specular_ref;  /* 鏡面反射率(RGB) */
  float shininess;        /* 光沢度 */

  material_type type;     /* マテリアルタイプ */
  colorf_t reflect_ref;   /* 完全鏡面反射光/屈折光係数(RGB) */
  float refraction_index; /* 絶対屈折率 */
} material_t;

• type
  • 質感のタイプを指定する．
    • MT_DEFAULTなら通常の(完全鏡面反射を使用しない)質感
    • MT_PERFECT_REFなら完全鏡面反射を使用する質感
    • MT_REFRACTIONなら完全鏡面反射/屈折を使用する質感
• reflect_ref
  • 完全鏡面反射光/屈折光係数を指定する．
• refraction_index
  • その物体の絶対屈折率

また，scene構造体にも以下の変更を施す．

typedef struct
{
  shape_t *shapes;              /* 物体リストへのポインタ */
  size_t num_shapes_capacity;   /* 物体リストの最大格納数 */
  size_t num_shapes;            /* 物体リストに実際に格納されている数 */
  light_t *lights;              /* 光源リストへのポインタ */
  size_t num_lights_capacity;   /* 光源リストの最大格納数 */
  size_t num_lights;            /* 光源リストに実際に格納されている数 */
  colorf_t ambient_illuminance; /* 環境光の強さ(RGB) */
  float global_refraction_index;/* 大気の絶対屈折率 */
} scene_t;

• global_refraction_index
  • 大気の絶対屈折率

raytrace関数の変更

raytrace関数(recursive_raytrace関数)を再帰呼び出しすることで光の屈折を実装する部分は，完全鏡面反射の場合と同じである．

レイの経路に関する注意事項

光の屈折の実装では，今までの実装の中で初めてレイが物体面の裏側から進入する場合がある(図8)．

図8. レイが物体面の裏側から進入する

この場合，交点における法線ベクトルを反転する必要がある(図9)．

図9. 法線ベクトルを反転する

視線ベクトルの進入方向を判定するには，視線ベクトルの逆ベクトル$$\vec{\bf v}$$と (反転を施す前の元の)法線ベクトル$$\vec{\bf n}$$との内積$$(\vec{\bf v}\cdot\vec{\bf n})$$をとり， その値の正負をチェックすれば良い． 値が正ならばレイは物体の表面からの進入しており，負ならば物体の裏面から進入している．

レイの進入方向によって，絶対屈折率$$\eta_{1},\eta_{2}$$の値の取り方も変わる． $$\eta_{1}$$は屈折前の媒質の絶対屈折率，$$\eta_{2}$$は屈折後の媒質の絶対屈折率である．

レイが物体の表面から進入した場合は，$$\eta_{1}$$＝大気の絶対屈折率，$$\eta_{2}$$＝物体の絶対屈折率となる．
レイが物体の裏面から進入した場合は，$$\eta_{1}$$＝物体の絶対屈折率，$$\eta_{2}$$＝大気の絶対屈折率となる．

(このことから分かるように，今回の実装では光線の経路は「大気→物体→大気」の順となることを前提としている． 物体から大気以外の物質に光線が出て行くことを想定していない．物体と視線との交点を見つけた場合に，それが大気との接点なのか他の物体との 接点なのかを判定するための処理はやや複雑になる．)

計算の効率化

屈折ベクトルおよび完全鏡面反射率の計算は，式を変形することで効率化できる．

$$\displaystyle \eta_{r}=\frac{\eta_{2}}{\eta_{1}}$$，$$\Omega=\eta_{r}cos\theta_{2}-cos\theta_{1}$$と置いて，

入射ベクトル$$\vec{\bf i}$$はここでは視線ベクトル$$\vec{\bf d_{\cal e}}$$なので，それに対する屈折ベクトル$$\vec{\bf f_{\cal e}}$$は以下となる．

また$$cos\theta{1}$$，$$cos\theta{2}$$は以下となる．

$$\vec{\bf v}$$は視線ベクトルの逆ベクトルである．すなわち，$$\vec{\bf v}=-\vec{\bf d_{\cal e}}$$となる．

課題 [—EXCLUDE—]

課題3-2. 光の屈折を実装する

光の屈折の処理を実装せよ．

課題3-2. 完全鏡面反射を実装するの ソースコードをディレクトリごとコピーしてから作業を行うこと(raytracing_tk_ch9という名前でコピー)．

シーンの設定は以下のようにせよ(図10)． 課題3-2. 完全鏡面反射を実装するとほぼ同じだが完全鏡面反射する球の位置と大きさが異なり，光を透過する球が追加されている．
scene_setting5.cとして新たにscene_setting関数を作成すること(Makefileの設定の追記も忘れないように注意する)．

図10. シーン設定

• 物体は全部で7つ
  1. 球1(完全鏡面反射)
     • 位置：$$(-0.4, -0.65, 3)$$，半径：$$0.35$$
     • 質感は完全鏡面反射のみ
  2. 球2(屈折)
     • 位置：$$(0.5, -0.65, 2)$$，半径：$$0.35$$
     • 絶対屈折率$$1.51$$
     • 質感は完全鏡面反射／屈折のみ
  3. 平面1(白い床)
     • 法線：$$(0,1,0)$$，位置：$$(0,-1,0)$$, 拡散反射のみ
  4. 平面2(白い天井)
     • 法線：$$(0,-1,0)$$，位置：$$(0,1,0)$$, 拡散反射のみ
  5. 平面3(右の緑の壁)
     • 法線：$$(-1,0,0)$$，位置：$$(1,0,0)$$, 拡散反射のみ
  6. 平面4(左の赤の壁)
     • 法線：$$(1,0,0)$$，位置：$$(-1,0,0)$$, 拡散反射のみ
  7. 平面5(白い奥壁)
     • 法線：$$(0,0,-1)$$，位置：$$(0,0,5)$$, 拡散反射のみ
• 光源は一つ
  1. 点光源
     • 位置：$$(0,0.9,2.5)$$，照度：$$1.0$$
• シーン
  • 大気の絶対屈折率は$$1.0$$

生成画像は以下のようになる．

図11. 生成画像

球2のチェックポイントとしては，

• 球の輪郭側にいくにつれて完全鏡面反射が強くなる
• 球全体としては後ろの風景の倒立像がみえる

という点に着目すること．

プログラムのビルド方法および実行方法

ヒント１

ヒント２