ADAPTIVE MATRIX COMPLETION FOR FAST VISIBILITY COMPUTATIONS WITH - - PowerPoint PPT Presentation

adaptive matrix completion
SMART_READER_LITE
LIVE PREVIEW

ADAPTIVE MATRIX COMPLETION FOR FAST VISIBILITY COMPUTATIONS WITH - - PowerPoint PPT Presentation

ADAPTIVE MATRIX COMPLETION FOR FAST VISIBILITY COMPUTATIONS WITH MANY LIGHTS RENDERING Sunrise Wang and Nicolas Holzschuch 1 Many-Lights Problem Direct lighting problem Large numbers of point lights Instant Radiosity (Keller


slide-1
SLIDE 1

ADAPTIVE MATRIX COMPLETION FOR FAST VISIBILITY COMPUTATIONS WITH MANY LIGHTS RENDERING

Sunrise Wang and Nicolas Holzschuch

1

slide-2
SLIDE 2

Many-Lights Problem

■ Direct lighting problem ■ Large numbers of point lights – Instant Radiosity (Keller '97)

■ VPLs approximate indirect lighting

■ Can range from thousands to millions of lights – Expensive – Accurate visibility needs ray-casting

2

slide-3
SLIDE 3

Contributions

■ Fast Many-Lights framework – Clustering – Visibility Approximation ■ Adaptive matrix completion to approximate visibility ■ Improvements to Adaptive Matrix Completion for boolean visibility

3

slide-4
SLIDE 4

Scalable Many-Lights Methods

■ Hierarchical clustering (Walter et al. '05, Walter et al. '06, Walter et al. '12, Bus et al. '15.) – Hierarchically clusters the VPLs – Traverses tree to extract clusters ■ Matrix-based methods (Hasan et al. '07, Ou & Pellacini '11) – Formulates problem as matrix – Clusters VPLs with information from sampling matrix

4

slide-5
SLIDE 5

Many-Lights Matrix

■ Lights as columns, receiving points as rows ■ Final colour of receiver = sum of values in row

Lights

R e c e i v e r s

5

slide-6
SLIDE 6

Many-Lights Matrix

■ Approximately low-rank – Most information lies within few basis columns ■ Some existing work – Generalized Nystrom method for inverse rendering (Wang et al. '09) – Low-rank and sparse matrix recovery on Many-Lights Matrix (Huo et al. '15) – Matrix completion on participating media Many-Lights Matrix (Huo et al. '16)

6

slide-7
SLIDE 7

Low-rank Matrix Completion

■ Complete matrix from partially observed coefficients ■ Many or infinite solutions, recover lowest rank

Matrix Completion

7

slide-8
SLIDE 8

Low-rank Matrix Completion

■ Matrix Completion Algorithms – Nuclear-norm minimization (Candes & Recht '09, Cai & Candes '10) – Alternating Least Squares (Haldar & Hernando '09, Wen et al. '12, Tanner & Wei '16) – Low-rank approximation (Goreinovet al. '97, Krishnamurthy& Singh '14) ■ Factors that influence required initial observations – Rank (number of linearly independent dimensions) – Coherence (localization of features)

8

slide-9
SLIDE 9

Naive Matrix Completion of the Many- Lights Matrix

■ Limited to matrices with lower rank – Diffuse scenes with relatively simple occlusion – Participating media ■ Difficult to determine rank and coherence for arbitrary scenes ■ Materials and geometry influence rank and coherence – Glossier scenes & complex visibility = higher rank and coherence

9

slide-10
SLIDE 10

Completing visibility

■ Can separate the many-lights matrix as M = S • V – S is shading matrix without visibility – V is a boolean visibility matrix –

  • is element-wise multiplication

■ Advantages – Shading easy to hardware accelerate – Fewer factors impacting rank and coherence of V – Allows for improvements to matrix completion algorithm

10

slide-11
SLIDE 11

Our framework

Generate VPLs Cluster VPLs Generate Receivers Slice Receivers Shade without visibility For each h slice Combine Slices Generate final slice image Complete visibility

11

Our contri ribution

  • n
slide-12
SLIDE 12

Adaptive Matrix Completion (Krishnamurty A. & Singh A. '14)[7]

■ Iteratively completes sub-sampled columns – Using basis – If can't, full sample & expand basis ■ Advantages – Relaxes row space coherence constrants – Operates on smaller basis -> faster

12

slide-13
SLIDE 13

Adaptive Matrix Completion for Boolean visibility

■ Blue = low rank, red = high rank, green = high coherence ■ Low-rank areas have near identical columns

13

slide-14
SLIDE 14

Adaptive Matrix Completion for Boolean visibility

■ Can approximate columnusing a single basis vector – Replace pseudo-inverse with matching – O(mn) rather than O(mn2) ■ Compared against pseudo-inverse of Q over ℝ and Gauss-Jordan elimination

  • f matrix over GF(2)

14

slide-15
SLIDE 15

Adaptive Matrix Completion for Boolean visibility

15

slide-16
SLIDE 16

Adaptive Matrix Completion for Boolean visibility

16

slide-17
SLIDE 17

Importance sampling

■ More samples where visibility changes ■ Use information from previous columns ■ Maintain discrete distribution d – Update d when column fully sampled

17

slide-18
SLIDE 18

Results

■ Over 12 scenes of varying complexity – VPLs in more diffuse scenes, VSLs in glossier scenes ■ Ground truth obtained by brute force ■ Compared to Low-rank and sparse separation (Huo et al. '15), IlluminationCut(Bus et

  • al. '15), Lightslice (Ou & Pellacini '11), and our framework without matrix completion

  • nly LightSlice and our method without matrix completion in VSL scenes

18

slide-19
SLIDE 19

Results – Bathroom (185s, 23% samples)

19

slide-20
SLIDE 20

Results– San Miguel (125s, 16.64% samples)

20

slide-21
SLIDE 21

Results – Sponza (140s, 11.81%)

21

slide-22
SLIDE 22

Results – San Miguel

  • L1 equal time error images of Low rank and sparse recovery (Huo et al. '15) (first) versus our method

(second)

22

1

slide-23
SLIDE 23

Visibility matrix rank

■ Rank is higher in penumbra – Receivers visibile varies slightly per light ■ Compared to Many-Lights matrix – Rank higher in diffuse scenes – Lower in glossier scenes

23

slide-24
SLIDE 24

Drawbacks

■ Inefficient column sampling

24

slide-25
SLIDE 25

Drawbacks

■ Clustering method does not consider visibility

25

slide-26
SLIDE 26

Conclusions

■ Method for speeding up Many-Lights renders – clustering & visibility approximation ■ Approximates visibility with Adaptive Matrix Completion ■ Improve Adaptive Matrix Completion for visibility ■ Over 3 times faster for equal quality

26

slide-27
SLIDE 27

Future Work

■ Address to drawbacks ■ Extend algorithm to area lights ■ Integrating over other domains?

– Time, Spectrum, Lens etc.

■ Matrix completion in other contexts

– Path-tracing, denoising etc.

27

slide-28
SLIDE 28

Quest stions? s?

28

slide-29
SLIDE 29

Adapting sample-rate

■ Slices have different ranks ■ Inefficient to use same sample-rate ■ Use verification samples to adapt

29