Computer graphics III Rendering equation and its solution. Jaroslav Křivánek, MFF UK

Similar documents
Math 656 Midterm Examination March 27, 2015 Prof. Victor Matveev

:2;$-$(01*%<*=,-./-*=0;"%/;"-*

9.5 Complex variables

Review - Probabilistic Classification

Grand Canonical Ensemble

Heisenberg Model. Sayed Mohammad Mahdi Sadrnezhaad. Supervisor: Prof. Abdollah Langari

8-node quadrilateral element. Numerical integration

Jones vector & matrices

LECTURE 5 Guassian Wave Packet

Modern Physics. Unit 5: Schrödinger s Equation and the Hydrogen Atom Lecture 5.6: Energy Eigenvalues of Schrödinger s Equation for the Hydrogen Atom

Signals and Systems View Point

element k Using FEM to Solve Truss Problems

Reflection models. Rendering equation. Taxonomy 2. Taxonomy 1. Digital Image Synthesis Yung-Yu Chuang 11/16/2006

Chapter 6 : Gibbs Free Energy

Text: WMM, Chapter 5. Sections , ,

Topic 5: Discrete-Time Fourier Transform (DTFT)

September 27, Introduction to Ordinary Differential Equations. ME 501A Seminar in Engineering Analysis Page 1. Outline

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.

EE 204 Lecture 25 More Examples on Power Factor and the Reactive Power

ELEC 372 LECTURE NOTES, WEEK 11 Dr. Amir G. Aghdam Concordia University

Self-Adjointness and Its Relationship to Quantum Mechanics. Ronald I. Frank 2016

N J of oscillators in the three lowest quantum

Chapter 2 Linear Waveshaping: High-pass Circuits

The Hyperelastic material is examined in this section.

Types of Communication

Chapter 3, Solution 1C.

CHAPTER 7d. DIFFERENTIATION AND INTEGRATION

The graph of y = x (or y = ) consists of two branches, As x 0, y + ; as x 0, y +. x = 0 is the

5 Curl-free fields and electrostatic potential

Chapter 33 Gauss s Law

EE 119 Homework 6 Solution

Consider a system of 2 simultaneous first order linear equations

Even/Odd Mode Analysis of the Wilkinson Divider

Wp/Lmin. Wn/Lmin 2.5V

120~~60 o D 12~0 1500~30O, 15~30 150~30. ..,u 270,,,, ~"~"-4-~qno 240 2~o 300 v 240 ~70O 300

CHAPTER 33: PARTICLE PHYSICS

Physics 256: Lecture 2. Physics

The Fourier Transform

IGEE 401 Power Electronic Systems. Solution to Midterm Examination Fall 2004

Partial Derivatives: Suppose that z = f(x, y) is a function of two variables.

. This is made to keep the kinetic energy at outlet a minimum.

Microfacet models for reflection and refraction

Introduction to logistic regression

From Structural Analysis to FEM. Dhiman Basu

Study of Dynamic Aperture for PETRA III Ring K. Balewski, W. Brefeld, W. Decking, Y. Li DESY

Ερωτήσεις και ασκησεις Κεφ. 10 (για μόρια) ΠΑΡΑΔΟΣΗ 29/11/2016. (d)

Phy213: General Physics III 4/10/2008 Chapter 22 Worksheet 1. d = 0.1 m

A Note on Estimability in Linear Models

Lecture 26: Quadrature (90º) Hybrid.

Analyzing Frequencies

Abstract Interpretation: concrete and abstract semantics

Finite Element Models for Steady Flows of Viscous Incompressible Fluids

Static/Dynamic Deformation with Finite Element Method. Graphics & Media Lab Seoul National University

Chapter 6 Student Lecture Notes 6-1

Unit 3: Transistor at Low Frequencies

Math 656 March 10, 2011 Midterm Examination Solutions

Microfacet models for reflection and refraction

EAcos θ, where θ is the angle between the electric field and

Another Explanation of the Cosmological Redshift. April 6, 2010.

Section 11.6: Directional Derivatives and the Gradient Vector

167 T componnt oftforc on atom B can b drvd as: F B =, E =,K (, ) (.2) wr w av usd 2 = ( ) =2 (.3) T scond drvatv: 2 E = K (, ) = K (1, ) + 3 (.4).2.2

ELEG 413 Lecture #6. Mark Mirotznik, Ph.D. Professor The University of Delaware

Unfired pressure vessels- Part 3: Design

Lucas Imperfect Information Model

Introduction to Arithmetic Geometry Fall 2013 Lecture #20 11/14/2013

Math 34A. Final Review

Minimum Spanning Trees

Lecture 14. Relic neutrinos Temperature at neutrino decoupling and today Effective degeneracy factor Neutrino mass limits Saha equation

Random Process Part 1

Thomas Whitham Sixth Form

Division of Mechanics Lund University MULTIBODY DYNAMICS. Examination Name (write in block letters):.

TMMI37, vt2, Lecture 8; Introductory 2-dimensional elastostatics; cont.

Introduction to Condensed Matter Physics

Transient Conduction: Spatial Effects and the Role of Analytical Solutions

u x v x dx u x v x v x u x dx d u x v x u x v x dx u x v x dx Integration by Parts Formula

Physic 231 Lecture 33

Mathematics. Complex Number rectangular form. Quadratic equation. Quadratic equation. Complex number Functions: sinusoids. Differentiation Integration

10.5 Linear Viscoelasticity and the Laplace Transform

CIRCUIT ANALYSIS II Chapter 1 Sinusoidal Alternating Waveforms and Phasor Concept. Sinusoidal Alternating Waveforms and

MATHEMATICS FOR MANAGEMENT BBMP1103

6. Negative Feedback in Single- Transistor Circuits

10/7/14. Mixture Models. Comp 135 Introduction to Machine Learning and Data Mining. Maximum likelihood estimation. Mixture of Normals in 1D

Is current gain generally significant in FET amplifiers? Why or why not? Substitute each capacitor with a

The Language of SOCIAL MEDIA. Christine Dugan

Outlier-tolerant parameter estimation

SAMPLE LITANY OF THE SAINTS E/G. Dadd9/F. Aadd9. cy. Christ, have. Lord, have mer cy. Christ, have A/E. Dadd9. Aadd9/C Bm E. 1. Ma ry and. mer cy.

Exam 1. It is important that you clearly show your work and mark the final answer clearly, closed book, closed notes, no calculator.

ECE542, Fall 2004 Homework 7 Solutions

Thomas Whitham Sixth Form

CTN 2/23/16. EE 247B/ME 218: Introduction to MEMS Design Lecture 11m2: Mechanics of Materials. Copyright 2016 Regents of the University of California

Equivalent Voice-coil Models for Real-time Computation in Electromagnetic Actuation and Sensor Applications

Basic Polyhedral theory

Soft k-means Clustering. Comp 135 Machine Learning Computer Science Tufts University. Mixture Models. Mixture of Normals in 1D

Physics 107 HOMEWORK ASSIGNMENT #20

Background: We have discussed the PIB, HO, and the energy of the RR model. In this chapter, the H-atom, and atomic orbitals.

V. Electrostatics Lecture 27a: Diffuse charge at electrodes

Module B3. VLoad = = V S V LN

EEC 686/785 Modeling & Performance Evaluation of Computer Systems. Lecture 12

1973 AP Calculus AB: Section I

Function Spaces. a x 3. (Letting x = 1 =)) a(0) + b + c (1) = 0. Row reducing the matrix. b 1. e 4 3. e 9. >: (x = 1 =)) a(0) + b + c (1) = 0

Transcription:

Cmputr graphcs III Rndrng quatn and ts slutn Jarslav Křvánk MFF UK Jarslav.Krvank@mff.cun.cz

Glbal llumnatn GI Drct llumnatn 2

Rvw: Rflctn quatn Sum ntgral f cntrbutns vr th hmsphr: r d cs H f r d n r r Ttal utgng rad. Emttd radanc Rflctd rad.

Frm lcal rflctn t glbal lght transprt Rflctn quatn lcal rflctn Whr ds th ncmng radanc cm frm? Frm thr placs n th scn! r r - d cs H f r r Ra castng functn CG III NPGR010 - J. Křvánk 2015

Frm lcal rflctn t glbal lght transprt Plug fr nt th rflctn quatn Incmng radanc drps ut Outgng radanc at dscrbd n trms f at thr pnts n th scn d cs r H f r CG III NPGR010 - J. Křvánk 2015

Rndrng quatn Rmv th subscrpt frm th utgng radanc: Dscrptn f th stad stat nrg balanc n th scn Rndrng calculat fr all pnts vsbl thrugh pls such that t fulfls th rndrng quatn d cs r H f r CG III NPGR010 - J. Křvánk 2015

Rflctn quatn Dscrbs lcal lght rflctn at a sngl pnt Intgral that can b usd t calculat th utgng radanc f w knw th ncmng radanc Rndrng quatn Cndtn n th glbal dstrbutn f lght n scn Intgral quatn unknwn quantt n bth sds Rflctn quatn vs. Rndrng quatn d cs r H f r d cs H f r Smlar frm dffrnt manng

Rndrng Equatn Kaja 1986 CG III NPGR010 - J. Křvánk 2015

Path tracng sktch

Rcursv unwndng f th RE Angular frm f th RE T calculat w nd t calculat r fr all drctns arund th pnt Fr th calculatn f ach r w nd t d th sam thng rcursvl tc. H r ' ' r r f r ' cs ' d' r CG III NPGR010 - J. Křvánk 2015

Path tracng v. zr rcursv frm gt ω: tracra ω rturn ω r ω // mttd radanc // rflctd radanc r ω: ω gnunfrmhmsphrrandmdr n rturn 2p * brdf ω ω * dtn ω * gt ω CG III NPGR010 - J. Křvánk 2015

Back t th thr: Angular and ara frm f th rndrng quatn

Angular vs. ara frm f th RE Angular frm ntgral vr th hmsphr n ncmng drctns Substtutn 2 cs d d r A d cs r H f r CG III NPGR010 - J. Křvánk 2015

Angular vs. ara frm f th RE Ara frm Intgral vr th scn surfac M r A V G f d 2 cs cs G vsblt 1 vsbl frm 0 thrws gmtr trm scn surfac CG III NPGR010 - J. Křvánk 2015

Angular frm Add radanc cntrbutns t a pnt frm all drctns Fr ach drctn fnd th narst surfac Implmntatn n stchastc path tracng: Fr a gvn gnrat randm drctns fr ach fnd th narst ntrsctn rturn th utgng radanc at that ntrsctn and multpl t wth th csn-wghtd BRDF. Avrag th rsult f ths calculatn vr all th gnratd drctns vr th hmsphr. CG III NPGR010 - J. Křvánk 2015

Ara frm Sum up cntrbutns t a pnt frm all thr pnts n th scn surfac Cntrbutn addd nl f th tw pnts ar mutuall vsbl Implmntatn n stchastc path tracng: Gnrat randml pnt n scn gmtr. Tst vsblt btwn and. If mutuall vsbl add th utgng radanc at mdulatd b th gmtr factr. Tpcal us: drct llumnatn calculatn fr ara lght surcs CG III NPGR010 - J. Křvánk 2015

Mst rndrng algrthms apprmat slutn f th RE cal llumnatn OpnG Onl pnt surcs ntgral bcms a sum Ds nt calculat qulbrum radanc s nt rall a slutn f th RE Fnt lmnt mthds radst [Gral 84] Dscrtz scn surfac fnt lmnts Dsrgard drctnalt f rflctns: vrthng s assumd t b dffus Cannt rprduc glss rflctns CG III NPGR010 - J. Křvánk 2015

Mst rndrng algrthms apprmat slutn f th RE Ra tracng [Whttd 80] Drct llumnatn n dffus and glss surfacs du t pnt surcs Indrct llumnatn nl n dal mrrr rflctn / rfractns Cannt calculat ndrct llumnatn n dffus and glss scns sft shadws tc. Dstrbutd ra tracng [Ck 84] Estmat th lcal rflctn usng th MC mthd Can calculat sft shadws glss rflctns camra dfcus blur tc. CG III NPGR010 - J. Křvánk 2015

Mst rndrng algrthms apprmat slutn f th RE Path tracng [Kaja 86] Tru slutn f th RE va th Mnt Carl mthd Tracng f randm paths randm walks frm th camra Can calculat ndrct llumnatn f hghr rdr CG III NPGR010 - J. Křvánk 2015

Frm th rndrng quatn t fnt lmnt radst

Frm th rndrng quatn t radst Start frm th ara frm f th RE: Th Radst mthd assumptns Onl dffus surfacs BRDF cnstant n and Radst.. radant tanc s spatall cnstant flat vr th ndvdual lmnts M r A V G f d CG III NPGR010 - J. Křvánk 2015

Frm th rndrng quatn t radst Dffus surfacs nl Th BRDF s cnstant n and Outgng radanc s ndpndnt f and t s qual t radst B dvdd b p M A V G d p M A V G B B B d p ' G CG III NPGR010 - J. Křvánk 2015

Frm th rndrng quatn t radst Spatall cnstant flat radst B f th cntrbutng surfac lmnts N B B B j G' da j j 1 Aj Radst f th j-th lmnt Gmtr factr btwn surfac lmnt j and pnt CG III NPGR010 - J. Křvánk 2015

Frm th rndrng quatn t radst Spatall cnstant flat radst f th rcvng surfac lmnt : B 1 A Avrag radst vr th lmnt A B da N 1 B B j G' da j da j 1 A A Aj da A j da j A j F j frm factr CG III NPGR010 - J. Křvánk 2015

Classc radst quatn Sstm f lnar quatns B B N j 1 B j F j Frm factrs 1 F j G' da j da A A A j Cnclusn: th radst mthd s nthng but a wa t slv th RE undr a spcfc st f assumptns CG III NPGR010 - J. Křvánk 2015

Radst mthd Classcal radst 1. Frm fact calculatn Mnt Carl hmcub 2. Slv th lnar sstm Gathrng Shtng Stchastc radst Avds plct calculatn f frm factrs Mtda Mnt Carl Radst s nt practcal nt usd Scn subdvsn -> snstv t th qualt f th gmtr mdl but n ralt mdls ar alwas brkn Hgh mmr cnsumptn cmpl mplmntatn CG III NPGR010 - J. Křvánk 2015

Th pratr frm f th RE

RE s a Frdhm ntgral quatn f th 2 nd knd Gnral frm th Frdhlm ntgral quatn f th 2 nd knd f g k f d unknwn functn Rndrng quatn: knwn functns quatn krnl H r f r cs d CG III NPGR010 - J. Křvánk 2015

nar pratrs nar pratrs act n functns as matrcs act n vctrs h f Th pratr s lnar f th actng s a lnar pratn af bg a f b g Eampls f lnar pratrs FIX NOTATION shul K f k f d f D f CG III NPGR010 - J. Křvánk 2015

Transprt pratr d cs H f r T Rndrng quatn T CG III NPGR010 - J. Křvánk 2015

Slutn f th RE n th pratr frm Rndrng quatn T Frmal slutn I T 1 I T unusabl n practc th nvrs cannt b plctl calculatd CG III NPGR010 - J. Křvánk 2015

Epansn f th rndrng quatn Rcursv substtutn T T T T T 2 n-fld rpttn lds th Numann srs n T 0 T n 1 CG III NPGR010 - J. Křvánk 2015

Epansn f th rndrng quatn If T s a cntractn tj. T < 1 whch hlds fr th RE thn 1 lm T n n 0 Slutn f th rndrng quatn s thn gvn b 0 T CG III NPGR010 - J. Křvánk 2015

A dffrnt drvatn f th Numann srs Frmal slutn f th rndrng quatn Prpstn Prf... 2 1 T T I T I I T T T T T I T T I T I T I T I......... 3 2 2 2 1 1 T I CG III NPGR010 - J. Křvánk 2015

Rndrng quatn T T Slutn: Numann srs T 2 3 T T... CG III NPGR010 - J. Křvánk 2015

Prgrssv apprmatn T T T T T T 2 3 T T T... T CG III NPGR010 - J. Křvánk 2015

Prgrssv apprmatn Each applcatn f T crrspnds t n stp f rflctn & lght prpagatn T 2 3 T T... mssn n-bunc ndrct llumnatn drct llumnatn OpnG shadng CG III NPGR010 - J. Křvánk 2015 tw-bunc ndrct llumnatn

Cntractvt f T Hlds fr all phscall crrct mdls Fllws frm th cnsrvatn f nrg It mans that rpttv applcatn f th pratr lwr th rmanng lght nrg maks sns snc rflctn/rfractn cannt crat nrg Scns wth wht r hghl spcular surfacs rflctvt cls t 1 t achv cnvrgnc w nd t smulat mr buncs f lght CG III NPGR010 - J. Křvánk 2015

Alrght s what hav w achvd? Rndrng quatn Slutn thrugh th Numann srs T 0 T W hav rplacd an ntgral quatn b a sum f smpl ntgrals Grat w knw hw t calculat ntgrals numrcall th Mnt Carl mthd whch mans that w knw hw t slv th RE and that mans that w can rndr mags a! Rcursv applcatn t T crrspnds t th rcursv ra tracng frm th camra CG III NPGR010 - J. Křvánk 2015

What act ntgral ar w valuatng thn? CG III NPGR010 - J. Křvánk 2015 M r M r M r A A A A A V G f V G f A V G f ž z z z ž z z z z d d...d d d d

Paths vs. rcursn: Sam thng dpnds n hw w lk at t Paths n a hgh-dmnsnal path spac T 2 3 T T... Rcursv slutn f a srs f nstd hmsphrcal ntgrals: T T T... CG III NPGR010 - J. Křvánk 2015

Rcursv ntrprtatn Angular frm f th RE T calculat I nd t calculat r fr all drctns arund th pnt. Fr th calculatn f ach r I nd t d th sam thng rcursvl tc. H r ' ' r r W v sn ths alrad rght? But unlk at th bgnnng f th lctur b nw w knw ths actuall slvs th RE. f r ' cs ' d' r CG III NPGR010 - J. Křvánk 2015

Path tracng v. 0.1 rcursv frm gt ω: narstintrsct ω rturn gt ω gtr ω // mttd radanc // rflctd radanc gtr ω nc : [ω gn pdf gn ] gnrnddrbrdfisω nc nrmal rturn 1/pdf gn * gt ω gn * brdf ω nc ω gn * dtnrmal ω gn CG III NPGR010 - J. Křvánk

Path tracng v. 2012 Arnld Rndrr 2012 Clumba Pcturs Industrs Inc. All Rghts Rsrvd. CG III NPGR010 - J. Křvánk 2015

Path tracng [Kaja86] Onl n scndar ra at ach ntrsctn Randm slctn f ntractn dffus rflctn rfractn tc. Drct llumnatn: tw stratgs Hp that th gnratd scndar ra hts th lght surc r Eplctl pck a pnt n th lght surc Trac hundrds f paths thrugh ach pl and avrag th rsult Advantag vr dstrbutd ra tracng: nw branchng f th ra tr mans n plsn f th numbr f ras wth rcursn dpth CG III NPGR010 - J. Křvánk 2015