On Skylight and Aerial Perspective. A.J. Preetham ATI Research

Transcription

1 On Skylight and Aerial Perspective A.J. Preetham ATI Research

2 Introduction

3 Outline Atmosphere Skylight Simulation models Analytic models Aerial Perspective Scattering using graphics hardware. Interactive and Real time techniques

4 Sun Solar constant: 1.37 kw/m 2 Radius: 680x10 6 m (109 times earth s radius) Solid angle: 7.2x10-5 sr -1 Color: White

5 Atmosphere A very thin layer and extends up to 600 km. 10% at 18 km 1% at 30 km 600km 6400km Earth Atmosphere

6 Atmosphere Layers Troposphere (12 km) Stratosphere (50 km) Mesosphere (85 km) Thermosphere (600 km)

7 US Standard Atmosphere Courtesy of Windows to the Universe

8 Atmosphere Sunlight Transmitted Scattered Absorbed Multiple Scattered Reflected Earth Atmosphere

9 Atmosphere Composition Mostly N 2 and O 2 Mixed gases, H 2 0, dust Absorption UV O 2, O 3 IR H 2 0, CH 4, C0 2, N 2 0 Oxygen 21% Others 1% Nitrogen 78%

10 Atmosphere Optical Mass. Mass of medium in path of unit cross-section. Optical Mass = Optical Length Optical Mass / Density at earth s surface. Optical Length = ρ(x) dx 1 ρ ρ 0 ( xdx )

11 Optical Length l 1 = ρ 0 ρ( xdx ) r(x) x r 0 r 0 Optical length

12 Optical length (cnt d)( Transmittance (T) T T = = e e 1 ρ K( x) dx 0 K 0 ρ ( xdx ) Scattering coefficient: Density: ρ Optical length: l Optical depth: τ K T = e Kl 0 T = e τ

13 Relative optical length Optical length varies with direction. Relative optical length l( θ ) m = l (0) Approximation by Kasten, m = (cosθ *( θ))

14 Relative optical lengths(cnt d) 8.4 km, 1.25 km s 11.9 km, 1.8 km 16.9 km, 2.54 km θ 45 o 60 o 75 o 32.6 km, 4.9 km 306 km, 45.6 km 90 o Earth Atmosphere Optical lengths for air, aerosols

15 Skylight

16 Skylight (cnt d)( Sunlight Atmosphere Earth

17 Skylight & Atmospheric Scattering Measurement (Meteorology) Simulation (Computer Graphics) Data Data Fit Model Fit Model Rendering Luminance CIE Clear sky CIE Overcast sky CIE Intermediate sky Perez all weather sky Spectral Radiance Preetham Spectral Radiance Blinn Dobashi Kaneda Klassen Nishita Max

18 Skylight (cnt d)( E s A D ω s B R dx Atmosphere L( ω) ω x C Earth Kl dl= Ee 0 ABKP ωω e 0 BCdx s 0 (, s ) Kl

19 Skylight (cnt d)( E s A D ω s B R dx Atmosphere L( ω) ω x C Earth Kl dl= Ee 0 ABKP ωω e 0 BCdx s 0 (, s ) Kl

20 Skylight (cnt d)( E s A D ω s B R dx Atmosphere L( ω) ω x C Earth Kl dl= Ee 0 ABKP ωω e 0 BCdx s 0 (, s ) Kl

21 Skylight (cnt d)( E s A D ω s B R dx Atmosphere L( ω) ω x C Earth Kl dl= Ee 0 ABKP ωω e 0 BCdx s 0 (, s ) Kl

22 Skylight (cnt d)( E s A D ω s B R dx Atmosphere L( ω) ω x C Earth Kl dl= Ee 0 ABKP ωω e 0 BCdx s 0 (, s ) Kl

23 Skylight (cnt d)( E s A D ω s B R dx Atmosphere L( ω) ω x C Earth Kl dl= Ee 0 ABKP ωω e 0 BCdx s 0 (, s ) Kl

24 Skylight (cnt d)( E s A D ω s B R dx Atmosphere L( ω) ω x C Earth D Kl 1 = s 0 C Kl L ( ω) Ee 0 ABKP( ωω, ) e 0 BC dx s

25 Skylight (cnt d)( E s A D ω s B R dx Atmosphere L( ω) ω x C Earth D Kl 1 = s 0 C Kl L ( ω) Ee 0 ABKP( ωω, ) e 0 BC dx s

26 Skylight (cnt d)( Atmosphere E s Earth L( ω) C A ω s ω L( ω,x) ω R B x dx D D 4π L ( ) L( ', xkp ) (, ') d ' e 0 ABdx ω ω ωω ω = 2 0 C 0 Kl

27 Skylight (cnt d)( Atmosphere E s Earth L( ω) C A ω s ω L( ω,x) ω R B x dx D D 4π L ( ) L( ', xkp ) (, ') d ' e 0 ABdx ω ω ωω ω = 2 0 C 0 Kl

28 Skylight (cnt d)( Atmosphere E s Earth L( ω) C A ω s ω L( ω,x) ω R B x dx D D 4π L ( ) L( ', xkp ) (, ') d ' e 0 ABdx ω ω ωω ω = 2 0 C 0 Kl

29 Skylight (cnt d)( Atmosphere E s Earth L( ω) C A ω s ω L( ω,x) ω R B x dx D D 4π L ( ) L( ', xkp ) (, ') d ' e 0 ABdx ω ω ωω ω = 2 0 C 0 Kl

30 Skylight (cnt d)( Atmosphere E s Earth L( ω) C A ω s ω L( ω,x) ω R B x dx D L( ω) = L ( ω) + L ( ω)

31 Atmospheric Scattering simulation Nishita et al, Display of The Earth Taking into account Atmospheric Scattering, Siggraph 1993

32 Atmospheric Scattering simulation Dobashi et al, Fast, accurate rendering method of Outdoor scenes using An All-weather skylight model, Journal IEEE of Japan, 1999

33 Skylight & Atmospheric Scattering Measurement (Meteorology) Simulation (Computer Graphics) Data Data Fit Model Fit Model Rendering Luminance CIE Clear sky CIE Overcast sky CIE Intermediate sky Perez all weather sky Spectral Radiance Preetham Spectral Radiance Blinn Dobashi Kaneda Klassen Nishita Max

34 Skylight Meteorological measurements IDMP (International Daylight Measurement Programme), 1991 was setup by CIE. Measurements included solar radiation and distribution of light over sky. Analytic models fit to real world measurement data.

35 Luminance Luminous flux (lumen) Energy of light emitted per sec in all directions Luminous intensity (cd) Luminous flux in a direction per unit solid angle Luminance (cd/m 2 ) Luminous intensity emitted per unit surface area Range: 1.6X10 9 cd/m 2 (sun at noon) to 8000cd/m 2 (clear sky) to 10-3 cd/m 2 (night sky)

36 Sky dome variables θ s θ v W γ N S φ s φ E

37 Overcast sky Luminance model: Moon & Spencer (1942), CIE (1955)

38 Overcast sky (cnt d)( CIE Luminance model: Y oc = Y z 1+ 2cosθ 3

39 Clear sky Luminance model: Pokrowski, Kittler (1967), CIE (1973).

40 Clear Sky (cnt d)( CIE luminance model Y c = Y z 3γ /cosγ e cos γ 1 e ( )( ) e cos 3θ s θs 1 e

41 All weather sky Luminance model: Perez et al Based on 5 different parameters Darkening or brightening of horizon Luminance gradient near horizon. Relative intensity of circum solar region Width of circum solar region Relative backscattered light

42 All weather sky (cnt d)( Perez luminance model B/cosθ Dγ 2 F( θγ, ) = (1 + Ae )(1+ Ce + Ecos γ) F( θγ, ) YP = YZ F (0, θ ) s

43 Skylight ASRC-CIE luminance model Linear combination of CIE clear sky Turbid clear sky (Gusev) Intermediate sky (Nakamura) CIE overcast sky The weights are computed using sky clearness and sky brightness factors.

44 Clear Sky Spectral radiance model: Preetham et al, Based on simulation data. Simulation based on USSA. Up-to second order scattering. Provides Luminance Y and chromaticities x,y similar to Perez model.

45 Clear sky (cnt d)( Spectral radiance model B/cosθ Dγ 2 F( θγ, ) = (1 + Ae )(1+ Ce + Ecos γ) F( θγ, ) Y = YZ F (0, θ ) s F( θγ, ) x= xz F (0, θ ) T : Turbidity Yz, xz, yz : functions of T, θs ABCDE,,,, : polynomials in s T F( θγ, ) y = yz F (0, θ ) s

46 Clear sky (cnt d)( Chromaticity to spectral curve x y M1 = x y M x y = x y S S MS M S D( λ) = 0( λ) + 1 1( λ) + 2 2( λ)

47 Clear Sky (cnt d)( Spectral curve & Luminance to spectral radiance Y = Km L( λ) V( λ) dλ V ( λ): Spectral sensitivity curve K m : Luminance efficacy(680 lm/w)

48 Clear Sky (cnt d)( Preetham et al, A Practical Analytical model for Daylight, Siggraph 1999.

49 Aerial Perspective Courtesy: Irene Alora & Pandromeda

50 Aerial Perspective (cnt d) Sunlight Inscattering Extinction

51 Aerial Perspective Simple model Earth is flat. Rays travel through constant density. Single scattering only.

52 Aerial Perspective (cnt d) L = fl0 + Lin C L 0 x L( ω' ) ω' B L in inscatter s ω L A

53 Aerial Perspective (cnt d) Kl f = e AC = e Ks C L 0 x L( ω' ) ω' B L in inscatter s ω L A

54 Aerial Perspective (cnt d) s Kl L e AB L( ω') KP( ωω, ') dω' dx in = 0 4π C L 0 x L( ω' ) ω' B L in inscatter s ω L A

55 Sunlight at earth s surface Transmittance Rayleigh scattering Aerosol scattering Water vapor absorption Ozone absorption Mixed gases absorption

56 Sunlight at earth s surface (cnt d)( Spectral attenuation Luminance v angle

57 Using graphics hardware

58 Interactive Rendering Dobashi et al (2000, 2002) Model atmosphere scattering by rendering virtual planes and accumulate terms using alpha blending Visibility at any virtual plane is calculated using standard shadow map technique.

59 Interactive Rendering (cnt d) D x E s x w s ω DL in Virtual planes Kl x in s s L = Vis( xe ) ( xe ) KP( ωω, ) x

60 Interactive Rendering (cnt d) x k + x L ( ) ( ) ( ) x in xk = Vis xes xe KP( ωω, s) dx x k Kl D x L ( x ) f ( x ) f ( x ) in k h k l k E s x ws x k + x fh( xk) = Vis( xe ) s( xdx ) x k ω DL in x k + x Kl f ( ) x l xk = e KP( ωω, s) dx x k fl planes f h planes

61 Interactive Rendering (cnt d) Dobashi et al, Interactive Rendering of Atmospheric Effects Using Graphics Hardware, Hardware Workshop 2002.

62 Real time rendering Programmable graphics hardware User specifies vertex and/or fragment program. Programmability exposed through DX8/9, GL ARB vertex and pixel shaders. Many offline shading models ported to graphics processor

63 Real time rendering (cnt d) Hoffman et al, 2002 Aerial Perspective - Evaluate f and L in in vertex shader Multiply f with L 0 and add L in in pixel shader L= fl0 + Lin C L 0 x L( ω' ) s ω' B L in inscatter ω L A

64 Real time rendering (cnt d) Transmittance term Kl f = e AC = e Ks In the presence of 2 types of particles ( K K ) s f e r + r+ m = m

65 Real time rendering (cnt d) Inscattering term. s K( s x) in = s ωωs 0 L e EKP(, ) dx Ks L = EP( ωω, )(1 e ) in s s In the presence of 2 types of particles KP Lin = E ( r m) s e + K + K r r ( ωω, s ) + KmPm ( ωω, s ) (1 ( Kr+ Km) s ) r m

66 Real time rendering (cnt d) Position L 0 K r, K E s ω s m Vertex Shader L f in Pixel Shader L= fl + L 0 in

67 Real time rendering (cnt d) Hoffman et al, Rendering Outdoor Light Scattering in Real Time, GDC 2002.

68 Real time rendering (cnt d) DEMO

69 Summary Atmosphere Skylight Simulation models Analytic models Aerial Perspective Scattering using graphics hardware. Interactive and Real time techniques

70 THANK YOU!

71