Wavelet Methods For Computer Graphics

This thesis discusses how a wavelet basis can be used in the context of two computer graphics applications, realistic rendering and geometric modeling, to produce more e cient and exible algorithms. The goal of realistic rendering is to simulate the interre ection of light in some geometric environment to produce realistic images. Radiosity is a commonly used solution method for this problem. Recently Hanrahan et al. have introduced a hierarchical method that can solve radiosity problems in O(n) time instead of O(n2). This thesis explores how the hierarchical radiosity algorithm can be formally understood from the context of wavelet theory. When the radiosity problem is expressed with respect to a wavelet basis, the resulting linear system is sparse, with only O(n) signi cant terms. By casting the hierarchical method in this framework, a variety of wavelet basis functions can be used, resulting in more e cient radiosity methods. This thesis also discusses how wavelets can be used in the context of geometric modeling. Geometric modeling is the study of how geometric shapes can be represented and manipulated by a designer. This thesis explores the use of wavelets to represent parametric curves and surfaces within the context of geometric modeling interfaces. One intuitive modeling interface commonly used in geometric modeling allows the user to directly manipulate curves and surfaces. This manipulation de nes some set of constraints that the curve or surface must satisfy (such as interpolation and tangent constraints). Direct manipulation, however, usually leads to an underconstrained problem since there are, in general, many possible surfaces meeting some set of constraints. Therefore an optimization problem must be solved. iv This thesis discusses how geometric modeling optimization problems can be solved more e ciently by using a wavelet basis. Because the wavelet basis is hierarchical, iterative optimization methods converge rapidly. And because the wavelet coe cients indicate the degree of detail in the solution, they can be used to determine the number of basis functions needed to express the variational minimum, thus avoiding unnecessary computation. An implementation is discussed and experimental results are reported. v Acknowledgments None of the work discussed in this thesis would have been possible without the mostexcellent guidance of my advisor, Professor Michael Cohen. He has given me insight, guidance, and inspiration every step of the way. I also would like to thank the readers of this thesis, Professors Pat Hanrahan and Elisha Sacks for their important input to this thesis. The radiosity research described in this thesis was done in collaboration with Professor Hanrahan, and his comments at the fourth oor library meetings helped us focus on the important research issues. Professor Sacks' guidance during my early years at Princeton are sincerely appreciated. The radiosity research was done in close collaboration with Peter Schr oder. It was a delightful experience working with him. James Shaw implemented the user interface for the variational modeling tool. His e ort was a great help. I also enjoyed yapping with him and his wife Colleen. Chuck Rose and Dan Boneh rounded out the yapping crew, and made the department almost a fun place to be. Dan Boneh was also very patient with all my math questions. I would like to thank David Ohsie who always had good advice. My various o ce mates throughout the years have enriched my Princeton experience, but certainly none more than Ayellet Tal. It looks like our study meetings for general exams have nally paid o . I would like to thank my parents for their loving support. The research reported in this thesis was made possible in part through support provided by the NSF (Grant No. CCR92-96146) and DARPA (Grant No. DABT6392-C-0053.P00002). vi

[1]  R. Nagel,et al.  3-D Visual simulation , 1971 .

[2]  Steven A. Orszag,et al.  CBMS-NSF REGIONAL CONFERENCE SERIES IN APPLIED MATHEMATICS , 1978 .

[3]  E. Catmull,et al.  Recursively generated B-spline surfaces on arbitrary topological meshes , 1978 .

[4]  J. Meinguet Multivariate interpolation at arbitrary points made simple , 1979 .

[5]  Turner Whitted,et al.  An improved illumination model for shaded display , 1979, CACM.

[6]  Tom Lyche,et al.  Discrete B-splines and subdivision techniques in computer-aided geometric design and computer graphics , 1980 .

[7]  Philip E. Gill,et al.  Practical optimization , 1981 .

[8]  James F. Blinn,et al.  A Generalization of Algebraic Surface Drawing , 1982, TOGS.

[9]  Donald Meagher,et al.  Geometric modeling using octree encoding , 1982, Comput. Graph. Image Process..

[10]  Lance Williams,et al.  Pyramidal parametrics , 1983, SIGGRAPH.

[11]  Peter R. Atherton,et al.  A scan-line hidden surface removal procedure for constructive solid geometry , 1983, SIGGRAPH.

[12]  Donald P. Greenberg,et al.  Modeling the interaction of light between diffuse surfaces , 1998 .

[13]  Robert L. Cook,et al.  Distributed ray tracing , 1984, SIGGRAPH.

[14]  Donald P. Greenberg,et al.  The hemi-cube: a radiosity solution for complex environments , 1985, SIGGRAPH.

[15]  James T. Kajiya,et al.  The rendering equation , 1986, SIGGRAPH.

[16]  Demetri Terzopoulos,et al.  Regularization of Inverse Visual Problems Involving Discontinuities , 1986, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[17]  Demetri Terzopoulos,et al.  Image Analysis Using Multigrid Relaxation Methods , 1986, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[18]  Donald P. Greenberg,et al.  An Efficient Radiosity Approach for Realistic Image Synthesis , 1986, IEEE Computer Graphics and Applications.

[19]  Tim Van Hook Real-time shaded NC milling display , 1986, SIGGRAPH.

[20]  T. Boult Information-based complexity in nonlinear equations and computer vision (multivariate splines, optimal algorithms, zero-finding) , 1986 .

[21]  Harry Yserentant,et al.  On the multi-level splitting of finite element spaces , 1986 .

[22]  William L. Briggs,et al.  A multigrid tutorial , 1987 .

[23]  Andrew P. Witkin,et al.  Spacetime constraints , 1988, SIGGRAPH.

[24]  Fujio Yamaguchi,et al.  Curves and Surfaces in Computer Aided Geometric Design , 1988, Springer Berlin Heidelberg.

[25]  David R. Forsey,et al.  Hierarchical B-spline refinement , 1988, SIGGRAPH.

[26]  Donald P. Greenberg,et al.  A progressive refinement approach to fast radiosity image generation , 1988, SIGGRAPH.

[27]  Demetri Terzopoulos,et al.  The Computation of Visible-Surface Representations , 1988, IEEE Trans. Pattern Anal. Mach. Intell..

[28]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[29]  Stéphane Mallat,et al.  A Theory for Multiresolution Signal Decomposition: The Wavelet Representation , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

[30]  Claude Puech,et al.  A general two-pass method integrating specular and diffuse reflection , 1989, SIGGRAPH '89.

[31]  John C. Beatty,et al.  A technique for the direct manipulation of spline curves , 1989 .

[32]  Daniel R. Baum,et al.  Improving radiosity solutions through the use of analytically determined form-factors , 1989, SIGGRAPH.

[33]  Richard Szeliski,et al.  Fast Surface Interpolation Using Hierarchical Basis Functions , 1990, IEEE Trans. Pattern Anal. Mach. Intell..

[34]  R. Coifman,et al.  Fast wavelet transforms and numerical algorithms I , 1991 .

[35]  B. K. Alpert Construction of Simple Multiscale Bases for Fast Matrix Operations , 1991 .

[36]  J. Roulier,et al.  Designing faired parametric surfaces , 1991, Comput. Aided Des..

[37]  Pat Hanrahan,et al.  A rapid hierarchical radiosity algorithm , 1991, SIGGRAPH.

[38]  P. Heckbert Simulating Global Illumination Using Adaptive Meshing , 1991 .

[39]  Barry Fowler,et al.  Geometric manipulation of tensor product surfaces , 1992, I3D '92.

[40]  GermanyNumerische Mathematik,et al.  Multilevel Preconditioning , 1992 .

[41]  Ingrid Daubechies,et al.  Ten Lectures on Wavelets , 1992 .

[42]  C. Chui,et al.  Wavelets on a Bounded Interval , 1992 .

[43]  David Salesin Dani Lischinski Tony DeRose Reconstructing Illumination Functions with Selected Discontinuities , 1992 .

[44]  Carlo H. Séquin,et al.  Functional optimization for fair surface design , 1992, SIGGRAPH.

[45]  Andrew P. Witkin,et al.  Variational surface modeling , 1992, SIGGRAPH.

[46]  D. Forsey,et al.  Multi-Resolution Surface Approximation for Animation , 1992 .

[47]  I. Daubechies,et al.  Biorthogonal bases of compactly supported wavelets , 1992 .

[48]  Pat Hanrahan,et al.  Wavelet radiosity , 1993, SIGGRAPH.

[49]  Dani Lischinski,et al.  Combining hierarchical radiosity and discontinuity meshing , 1993, SIGGRAPH.

[50]  Pat Hanrahan,et al.  On the form factor between two polygons , 1993, SIGGRAPH.

[51]  Ronald R. Coifman,et al.  Wavelet-Like Bases for the Fast Solution of Second-Kind Integral Equations , 1993, SIAM J. Sci. Comput..

[52]  B. Alpert A class of bases in L 2 for the sparse representations of integral operators , 1993 .

[53]  Tony DeRose,et al.  Efficient, fair interpolation using Catmull-Clark surfaces , 1993, SIGGRAPH.

[54]  Ming-Haw Yaou,et al.  Fast surface interpolation using multi-resolution wavelet transform , 1993, Proceedings of IEEE Systems Man and Cybernetics Conference - SMC.

[55]  Michael F. Cohen,et al.  Radiosity and realistic image synthesis , 1993 .

[56]  Harold R. Zatz Galerkin radiosity: a higher order solution method for global illumination , 1993, SIGGRAPH.

[57]  Christophe Schlick,et al.  A Customizable Reflectance Model for Everyday Rendering , 1993 .

[58]  Sam Qian,et al.  Wavelets and the Numerical Solution of Partial Differential Equations , 1993 .

[59]  S. Jaffard,et al.  Orthonormal wavelets, analysis of operators, and applications to numerical analysis , 1993 .

[60]  Zicheng Liu,et al.  Hierarchical spacetime control , 1994, SIGGRAPH.

[61]  Pat Hanrahan,et al.  Textures and radiosity: controlling emission and reflection with texture maps , 1994, SIGGRAPH.

[62]  James Arvo,et al.  A clustering algorithm for radiosity in complex environments , 1994, SIGGRAPH.

[63]  Pat Hanrahan,et al.  Wavelet Projections for Radiosity , 1994, Comput. Graph. Forum.

[64]  Philipp Slusallek,et al.  Radiosity and relaxation methods , 1994, IEEE Computer Graphics and Applications.

[65]  David Salesin,et al.  Multiresolution curves , 1994, SIGGRAPH.

[66]  Wen-Thong Chang,et al.  Fast Surface Interpolation using Multiresolution Wavelet Transform , 1994, IEEE Trans. Pattern Anal. Mach. Intell..

[67]  E. Quak,et al.  Decomposition and Reconstruction Algorithms for Spline Wavelets on a Bounded Interval , 1994 .

[68]  Robert R. Lewis,et al.  Making Shaders More Physically Plausible , 1994, Comput. Graph. Forum.

[69]  François X. Sillion Clustering and Volume Scattering for Hierarchical Radiosity Calculations , 1995 .

[70]  Thomas Ertl,et al.  Computer Graphics - Principles and Practice, 3rd Edition , 2014 .