Modelling the Directionality of Light Scattered in Translucent - PowerPoint PPT Presentation

Modelling the Directionality of Light Scattered in Translucent Materials Jeppe Revall Frisvad Joint work with Toshiya Hachisuka, Aarhus University Thomas Kjeldsen, The Alexandra Institute March 2014

Materials (scattering and absorption of light) ◮ Optical properties (index of refraction, n ( λ ) = n ′ ( λ ) + i n ′′ ( λ ) ). ◮ Reflectance distribution functions, S ( x i , � ω i ; x o , � ω o ). BSSRDF n 1 x i x o n 2 BRDF n 1 n 2

Subsurface scattering ◮ Behind the rendering equation [Nicodemus et al. 1977] : n 1 d L r ( x o , � ω o ) ω i ) = S ( x i , � ω i ; x o , � ω o ) . x i x o dΦ i ( x i , � n 2 ◮ An element of reflected radiance d L r is proportional to an element of incident flux dΦ i . ◮ S (the BSSRDF) is the factor of proportionality. d 2 Φ ◮ Using the definition of radiance L = cos θ d A d ω , we have � � L r ( x o , � ω o ) = S ( x i , � ω i ; x o , � ω o ) L i ( x i , � ω i ) cos θ d ω i d A . A 2 π References - Nicodemus, F. E., Richmond, J. C., Hsia, J. J., Ginsberg, I. W., and Limperis, T. Geometrical considerations and nomenclature for reflectance. Tech. rep., National Bureau of Standards (US), 1977.

BRDF BSSRDF [Jensen et al. 2001]

[Donner and Jensen 2006]

Splitting up the BSSRDF ◮ Bidirectional Scattering-Surface Reflectance Distribution Function: S = S ( x i , � ω i ; x o , � ω o ) . ◮ Away from sources and boundaries, we can use diffusion. ◮ Splitting up the BSSRDF S = T 12 ( S (0) + S (1) + S d ) T 21 . where ◮ T 12 and T 21 are Fresnel transmittance terms (using � ω i , � ω o ). ◮ S (0) is the direct transmission part (using Dirac δ -functions). ◮ S (1) is the single scattering part (using all arguments). ◮ S d is the diffusive part (multiple scattering, using | x o − x i | ). ◮ We distribute the single scattering to the other terms using the delta-Eddington approximation: S = T 12 ( S δ E + S d ) T 21 , and generalize the model such that S d = S d ( x i , � ω i ; x o ).

Diffusion theory ◮ Think of multiple scattering as a diffusion process. ◮ In diffusion theory, we use quantities that describe the light field in an element of volume of the scattering medium. ◮ Total flux, or fluence, is defined by x d z y � d y φ ( x ) = L ( x , � ω ) d ω . 4 π d x x z ◮ We find an expression for φ by solving the diffusion equation ( D ∇ 2 − σ a ) φ ( x ) = − q ( x ) + 3 D ∇· Q ( x ) , where σ a and D are absorption and diffusion coefficients, while q and Q are zeroth and first order source terms.

Deriving a BSSRDF ◮ Assume that emerging light is diffuse due to a large number of scattering events: S d ( x i , � ω i ; x o , � ω o ) = S d ( x i , � ω i ; x o ). ◮ Integrating emerging diffuse radiance over outgoing directions, we find S d = C φ ( η ) φ − C E ( η ) D � n o ·∇ φ , Φ 4 π C φ (1 /η ) where ◮ Φ is the flux entering the medium at x i . ◮ � n o is the surface normal at the point of emergence x o . ◮ C φ and C E depend on the relative index of refraction η and are polynomial fits of different hemispherical integrals of the Fresnel transmittance. ◮ This connects the BSSRDF and the diffusion theory. ◮ To get an analytical model, we use a special case solution for the diffusion equation (an expression for φ ). ◮ Then, “all” we need to do is to find ∇ φ (do the math) and deal with boundary conditions (build a plausible model).

Point source diffusion or ray source diffusion our model standard dipole ◮ Ray source diffusion ◮ Point source diffusion [Menon et al. 2005a; 2005b] [Bothe 1941; 1942] e − σ tr r e − σ tr r Φ Φ φ ( r , θ ) = φ ( r ) = , 4 π D r 4 π D r 1 + 3 D 1+ σ tr r � � cos θ , r where r = | x o − x i | and � σ tr = σ a / D is the where θ is the angle effective transport between the refracted ray coefficient. and x o − x i .

Our BSSRDF when disregarding the boundary d ◮ Using x = x o − x i , r = | x | , cos θ = x · � ω 12 / r , we take the gradient of φ ( r , θ ) (the expression for ray source diffusion) and insert to find � � � 4 π 2 e − σ tr r r 2 1 1 S ′ d ( x , � ω 12 , r ) = C φ ( η ) D + 3(1 + σ tr r ) x · � ω 12 4 C φ (1 /η ) r 3 � � � � (1 + σ tr r ) + 3 D 3(1+ σ tr r )+( σ tr r ) 2 � − C E ( η ) 3 D (1 + σ tr r ) � ω 12 · � n o − x · � ω 12 x · � , n o r 2 which would be the BSSRDF if we neglect the boundary.

Dipole configuration (method of mirror images) d ◮ We place the “real” ray source at the boundary and reflect it in an extrapolated boundary to place the “virtual” ray source. ◮ Distance to the extrapolated boundary [Davison 1958]: � d e = 2 . 131 D / 1 − 3 D σ a . ◮ In case of a refractive boundary ( η 1 � = η 2 ), the distance is A = 1 − C E ( η ) with . Ad e 2 C φ ( η )

Modified tangent plane d ◮ The dipole assumes a semi-infinite medium. ◮ We assume that the boundary contains the vector x o − x i and that it is perpendicular to the plane spanned by � n i and x o − x i . ◮ The normal of the assumed boundary plane is then i = x o − x i | x o − x i | × � n i × ( x o − x i ) � n ∗ or � n ∗ i = � n i × ( x o − x i ) | , n i if x o = x i . | � and the virtual source is given by n ∗ n ∗ n ∗ x v = x i +2 Ad e � i , d v = | x v − x i | , � ω v = � ω 12 − 2( � ω 12 · � i ) � i .

Distance to the real source (handling the singularity) z -axis z r real source n dr d n 2 z = 0 r n 1 ω i ω o incident light dv observer -zv virtual source our model standard dipole r 2 + z 2 � d r = r . d r = r ? ◮ Emergent radiance is an integral over z of a Hankel transform of a Green function which is Fourier transformed in x and y . ◮ Approximate analytic evaluation is possible if r is corrected to R 2 = r 2 + ( z ′ + d e ) 2 . ◮ The resulting model for z ′ = 0 corresponds to the standard dipole where z ′ = z r and d e is replaced by the virtual source.

Distance to the real source (handling the singularity) � r 2 + d 2 e � n o d e β x o x i r θ − � n o � ω 12 θ 0 ◮ Since we neither have normal incidence nor x o in the tangent plane, we modify the distance correction: R 2 = r 2 + z ′ 2 + d 2 e − 2 z ′ d e cos β . ◮ It is possible to reformulate the integral over z to an integral along the refracted ray. ◮ We can approximate this integral by choosing an offset D ∗ along the refracted ray. Then z ′ = D ∗ | cos θ 0 | .

Our BSSRDF when considering boundary conditions ◮ Our final distance to the real source becomes � r 2 + D µ 0 ( D µ 0 − 2 d e cos β ) for µ 0 > 0 (frontlit) d 2 r = r 2 + 1 / (3 σ t ) 2 otherwise (backlit) , with µ 0 = cos θ 0 = − � n o · � ω 12 and � r 2 − ( x · ω 12 ) 2 r cos β = − sin θ = − . r 2 + d 2 r 2 + d 2 � e e ◮ The diffusive part of our BSSRDF is then S d ( x i , � ω i ; x o ) = S ′ d ( x o − x i , � ω 12 , d r ) − S ′ d ( x o − x v , � ω v , d v ) , while the full BSSRDF is as before: S = T 12 ( S δ E + S d ) T 21 .

Previous Models ◮ Previous models are based on the point source solution of the diffusion equation and have the problems listed below. 1. Ignore incoming light direction: ◮ Standard dipole [Jensen et al. 2001]. ◮ Multipole [Donner and Jensen 2005]. ◮ Quantized diffusion [d’Eon and Irving 2011]. 2. Require precomputation: ◮ Precomputed BSSRDF [Donner et al. 2009, Yan et al. 2012]. 3. Rely on numerical integration: ◮ Photon diffusion [Donner and Jensen 2007, Habel et al. 2013]. ◮ Using ray source diffusion, we can get rid of those problems.

Results (Grapefruit Bunnies) dipole ours reference quantized

Results (marble Bunnies) dipole ours reference quantized

Results (Simple Scene) ◮ Path traced single scattering RMSE dipole was added to the existing 0.3 btpole qntzd 0.2 ours models but not to ours. 0.1 0 ◮ Faded bars show quality apple marble potato skin1 choc milk soy milk w. grape milk (r.) SSIM measurements when single 1 0.9 scattering is not added. 0.8 ◮ The four leftmost materials 0.7 apple marble potato skin1 choc milk soy milk w. grape milk (r.) scatter light isotropically.

Results (2D plots, 30 ◦ Oblique Incidence) quantized ours ◮ Our model is significantly different ◮ when the angle of incidence changes ◮ when the direction toward the point of emergence changes.

Results (2D plots, 45 ◦ Oblique Incidence) quantized ours ◮ Our model is significantly different ◮ when the angle of incidence changes ◮ when the direction toward the point of emergence changes.

Results (2D plots, 60 ◦ Oblique Incidence) quantized ours ◮ Our model is significantly different ◮ when the angle of incidence changes ◮ when the direction toward the point of emergence changes.

Results (Diffuse Reflectance Curves) 10 1 dipole dipole btpole btpole qntzd 10 0 ours qntzd ptrace ours 10 -1 ptrace R (x) 10 -2 d 10 -3 10 -4 10 -5 -16 -12 -8 -4 0 4 8 12 16 ◮ Our model comes closer than the existing analytical models to measured and simulated diffuse reflectance curves.

Results (Image Based Lighting) quantized ours

The 3Shape Buddha! (scanned with a TRIOS Scanner) matte milk-coloured mini milk