MATHEMATICAL TOOLS IN COMPUTER GRAPHICS WITH C# IMPLEMENTATIONS

Transcription

1 MATHEMATICAL TOOLS IN COMPUTER GRAPHICS WITH C# IMPLEMENTATIONS

2 This page intentionally left blank

3 World Scientific N E W J E R S E Y L O N D O N S I N G A P O R E B E I J I N G S H A N G H A I H O N G K O N G TA I P E I C H E N N A I

4 Published by World Scientific Publishing Co. Pte. Ltd. 5 Toh Tuck Link, Singapore USA office: 27 Warren Street, Suite , Hackensack, NJ UK office: 57 Shelton Street, Covent Garden, London WC2H 9HE British Library Cataloguing-in-Publication Data A catalogue record for this book is available from the British Library. MATHEMATICAL TOOLS IN COMPUTER GRAPHICS WITH C# IMPLEMENTATIONS Copyright 2008 by World Scientific Publishing Co. Pte. Ltd. All rights reserved. This book, or parts thereof, may not be reproduced in any form or by any means, electronic or mechanical, including photocopying, recording or any information storage and retrieval system now known or to be invented, without written permission from the Publisher. For photocopying of material in this volume, please pay a copying fee through the Copyright Clearance Center, Inc., 222 Rosewood Drive, Danvers, MA 01923, USA. In this case permission to photocopy is not required from the publisher. ISBN ISBN ISBN (pbk) ISBN (pbk) Printed in Singapore.

5 Preface Computer graphics is applied in many different fields including the entertainment industry, medical and geographic visualization and industrial design. Successful development of computer graphics programs and algorithms requires knowledge from a diverse set of fields to be combined and applied effectively. The most important part of any computer graphics algorithm is the underlying mathematical tools that are used to analyze and develop the algorithm. This book introduces computer graphics from a mathematical perspective. The basic computer graphics principles are introduced as well as more advanced topics. The mathematical techniques are developed in detail. Selected topics include: basic transforms, curves, surfaces and subdivision surfaces. New techniques such as wavelets, fractals, parameterization and fluid simulation are also included. Throughout a large portion of the text a new curve and surface algorithm are developed to illustrate the use of mathematics to develop computer graphics algorithms. C# implementations for many of the algorithms are provided throughout the book. The C# implementation allows the reader to obtain practical knowledge along with the theoretical knowledge. The provided C# implementations can be downloaded from our websites for use in other software projects. Book websites addresses of the authors alexandre.hardy@gmail.com steebwilli@gmail.com whsteeb@uj.ac.za v

6 This page intentionally left blank

7 Contents List of Figures Notation xiii xvii 1 Vectors, Matrices and Transforms Vector Spaces Points and Vectors Homogeneous Coordinates Representing Objects by Points Affine Transformations Introduction and Definitions Scaling Translation Rotation Concatenation of Transforms Projection Quaternions C# Implementation Lighting Shading Affine Transforms and Normal Vectors Local Lighting Models The Phong Lighting Model Emissive Properties Ambient Reflection Diffuse Reflection Specular Reflection Multiple Colored Light Sources Attenuation vii

8 viii Contents 2.4 Lights Spot Lights Transparent Objects Cook-Torrance Model Bidirectional Reflectivity Cook-Torrance Model Microfacet Distribution Term Geometric Surface Occlusion Term Fresnel Term Beer-Lambert Law C# Implementation Rasterization Pixels Drawing Lines Bresenham s Algorithm for Lines Drawing Circles Bresenham s Algorithm for Circles Filling Gouraud Shading Rasterization in C# Drawing Pixels Bresenham s Algorithms in C# Fractals Mandelbrot Set Julia Set Iterated Function Systems L-Systems and Fractals Kronecker Product and Fractals Definitions Kronecker Product Fractals Curves Introduction Affine Invariance Convex Hull Lagrange Interpolation C# Implementation Bézier Curves Affine Invariance Convex Hull Derivative at Edges

9 Contents ix Piecewise Continuous Bézier Curves Rendering Rational Bézier Curves Bézier Curves: Conic Sections C# Implementation Catmull-Rom Splines Bessel-Overhauser Splines Tension-Continuity-Bias Splines Uniform B-Splines Affine Invariance Convex Hull Cox-de Boor Formula C# Implementation Non-Uniform B-Splines Interpolating with B-Splines Periodic Interpolation Non-Uniform Rational B-Splines Trigonometric Interpolation METAPOST and Bézier Curves METAPOST Example Curvature and Torsion Harmonic Interpolation Interpolation Odd Case Even Case Examples Curvature Plots Numerical Stability Affine Invariance Convex Hull Property C# Implementation of Harmonic Interpolation Chebyshev Polynomials Odd Case Even Case Non-Uniform Harmonic Interpolation Wavelets Introduction One-Dimensional Wavelets Two-Dimensional Wavelets Curves C# Implementation

10 x Contents 6 Surfaces Parametric Surfaces Tensor Product Surfaces Bézier Surfaces Tensor Product Bézier Surfaces Triangular Bézier Surfaces Rational Bézier Surfaces Bézier Surface Interpolation B-Spline Tensor Product Surfaces B-Spline Surface Interpolation Subdivision Surfaces Loop Subdivision Modified Butterfly Subdivision Subdivision Interpolating 3 Subdivision Catmull-Clark Subdivision Doo-Sabin Subdivision Comparison Interpolation with Subdivision Surfaces Curvature of Surfaces Harmonic Surfaces Tensor Product Surface Harmonic Subdivision Local Harmonic Subdivision Local Harmonic Interpolation for Curves Parametric Distance Subdivision Rules Irregular Vertices Boundaries Geometry Images and Parameterization Cutting a Mesh into a Disk Parameterization Rasterization of the Geometry Image Interpolation of Geometry Images Geometry Images Approximation Rendering Approximating Basis Functions Combined Results Curvature C# Implementation

11 Contents xi 7 Raytracing Raytracing Process Representation of a Ray Reflection Refraction Intersections Sphere Infinite Plane Triangles Effect of Transforms C# Implementation of a Raytracer Implicit Surfaces Sphere Tracing Distance Functions C# Implementation CSG Objects C# Implementation Parametric Surfaces Interval Arithmetic Interval Root Finding Bisection Interval Root Finding Newton-Raphson Ray Tracing Harmonic Surfaces Lighting Models Supersampling Regular Supersampling Stochastic Supersampling Adaptive Supersampling Ambient Occlusion Ray Marching Photon Mapping Transport Notation Path Tracing Creating the Photon Map Photon Tracing Photon Map Data Structure Radiance Estimate C# Implementation Radiosity Light Transport Notation Radiosity Matrix Solving for Radiosity Values

12 xii Contents Solving: Jacobi Method Solving: Gauss-Seidel Iteration Solving: Shooting Method Form Factors Numerical Solution Raytracing Method Hemicube Method Rendering C# Implementation Animation Traditional Animation Techniques Keyframing Motion Capture Physics Models Animation of Position Arc length parameterization Orientation Articulated Structures (Kinematics) Forward Kinematics Vertex Blending Inverse Kinematics Mass Spring Systems Particle Systems Free Form Deformations Fluids Navier-Stokes Equations Advection Diffusion Projection Boundary Conditions C# Implementation Free Surface C# Implementation of Free Surfaces Bibliography 459 Index 471

13 List of Figures 1.1 Rotation in the two-dimensional Cartesian plane Projection onto a plane The view frustum and canonical view volume Diffuse reflection Specular reflection and the reflection vector Phong lighting model Microfacets V shaped grooves Shadowing and masking Blinn s geometric lemma Phong compared to Cook-Torrance lighting model Mandelbrot Set Julia Set IFS Fern Tree generated with an L-System First four steps in the construction of the Hilbert curve First five steps in the construction of the Koch snowflake Sierpinski carpet after seven iterations First six iterations of the Sierpinski triangle Kronecker product fractal: Gray scale fractal First four iterations of checkerboard Basis functions for uniform Lagrange interpolation of degree Uniform Lagrange interpolation of 4 points Basis functions for Bézier curves of degree A Bézier curve and its convex hull Piecewise continuous Bézier curves de Casteljau s method Subdivision of figure 4.4 into two Bézier curves Section of a circle, drawn with a rational Bézier curve Basis functions of B-splines of degree Blending functions of B-splines of degree xiii

14 xiv List of figures 4.11 Bézier curve of degree B-spline curve of degree Interpolation using B-splines Periodic interpolation using B-splines Interpolation using METAPOST Periodic interpolation using METAPOST Catenary and its curvature Spiral curve σ(t) function for n = σ k (t) basis functions for n = σ(t) function for n = σ k (t) basis functions for n = σ(t) function for n = σ(t) function for n = Circle drawn by METAPOST and harmonic interpolation Circle drawn by B-spline and harmonic interpolation Seven sided polygon: METAPOST, harmonic interpolation Seven sided polygon: B-spline, harmonic interpolation Five sided polygon: METAPOST, harmonic interpolation Five sided polygon: B-spline, harmonic interpolation Four point polygon: Trigonometric interpolation D curves: Comparison Curvature plot of figure Curvature plot of figure Curvature plot of figure Curvature plot of an eight sided figure Stability of harmonic interpolation Non-uniform harmonic interpolation of seven points Seven sided polygon: Non-uniform harmonic interpolation Five sided polygon: Non-uniform harmonic interpolation Eight sided polygon: Non-uniform harmonic interpolation The box functions for V The Haar wavelets for W Example parametric surfaces Example Bézier surface patches Surface constructed from piecewise continuous Bézier patches A Bézier triangle control mesh B-spline surface fitting: Curve network Stages in subdivision using Chaikin s scheme Stages in interpolatory subdivision using the 4-point scheme Notation for Loop subdivision Triangle subdivision for Loop subdivision

15 List of figures xv 6.10 Masks for Loop subdivision Cube: Loop subdivision Star: Loop subdivision Limit surfaces (Loop subdivision) Masks for modified butterfly subdivision Cube: Modified butterfly subdivision Star: Modified butterfly subdivision Limit surfaces (modified butterfly subdivision) Masks for 3 subdivision Triangle subdivision for 3 subdivision Cube: 3 subdivision Star: 3 subdivision Limit surfaces ( 3 subdivision) Masks for interpolating 3 subdivision Subdivision for the Catmull-Clark scheme Cube (triangulated): Catmull-Clark subdivision Cube (quadrilateral): Catmull-Clark subdivision Star (triangulated): Catmull-Clark subdivision Star (quadrilateral): Catmull-Clark subdivision Limit surfaces (Catmull-Clark subdivision) Cube (triangulated): Doo-Sabin subdivision Cube (quadrilaterals): Doo-Sabin subdivision Star (triangulated): Doo-Sabin subdivision Star (quadrilaterals): Doo-Sabin subdivision Limit surfaces (Doo-Sabin subdivision) Subdivision of the mannequin head Tensor product harmonic surfaces Stages in harmonic subdivision Stages in subdivision to produce a torus Stages in subdivision to produce a deformed torus Basis functions for the 2D and 3D harmonic interpolation Harmonic interpolation of nearest 4 points Harmonic subdivision with local support Stencil for regular vertices Stencil for irregular vertices Subdivision of a simple chess piece Subdivision of a simple diamond model with a valence 12 vertex Subdivision for refining an icosahedron Subdivision for refining the mannequin head Parameterization of Mannequin Shape-preserving parameterization: Computation of weights 225

16 xvi List of figures 6.51 Mean value coordinates: Computation of weights Application of theorems 6.1 and 6.3 (with m = 1) Scanline-rendered harmonic surface (torus) Approximation of harmonic tensor product surface Bézier patches used to approximate figure Final rendering of mannequin head Bézier patches Final rendering of Spock formed from Bézier patches Control points of various parametric surfaces Comparison of plane and harmonic surface Comparison of partial sphere and harmonic surface Comparison of torus and harmonic surface Comparison of sphere and harmonic surface Comparison of Möbius strip and harmonic surface Comparison of surface of revolution and harmonic surface Comparison of shell and harmonic surface Comparison of another shell and harmonic surface Comparison of partial Klein bottle and harmonic surface Comparison of overhand knot and harmonic surface Comparison of Solomon s Seal knot and harmonic surface A scene and viewpoint Raytracing by casting rays from the viewer into the scene View of the scene from the camera Computation of the reflected vector r Computation of the refracted vector t Sphere tracing: Hit Sphere tracing: Miss Distance to a cone A superquadric Ray traced harmonic surface (torus) Supersampling Adaptive supersampling Ambient occlusion Global illumination with photon mapping Caustics The Hemicube method Determining radiosity values Examples of radiosity solutions An articulated structure A fluid cell in the MAC-grid

17 Notation N the set of natural numbers Z the set of integer numbers Q the set of rational numbers R the set of real numbers C the set of complex numbers (of the form a + bi, a, b R) R n the Cartesian product space, i.e. R n = { (x 0, x 1,..., x n 1 ) x j R } z complex number z complex conjugate of z C i i = 1 I n n n unit matrix M T transpose of the matrix M M, det(m) determinant of the matrix M x = (x 0 x 1... x n 1 ) T column vector x R n x T row vector x R n 0 zero column vector n normal vector t tangent vector x Euclidean norm of x, x = x T x x, y the inner product (scalar product) of x and y x y dot product, an inner product defined by x T y = n 1 k=0 x ky k, where x, y R n x y the vector product (cross product) of x R 3 and y R 3, xvii

18 xviii Notation ang(a, b, c) the angle between the vectors a b and c b [0, 1] unit interval δ jk Kronecker delta, 1 if j = k, 0 otherwise I n n n identity matrix t curve parameter u, v surface parameters C k continuity of a curve or surface G k geometric continuity of a curve or surface x floor of x R, the largest y such that y Z and y x x ceiling of x R, the smallest y such that y Z and y x Kronecker product