3.1 Viewing and Projection

Transcription

1 Fall 2017 CSCI 420: Computer Graphics 3.1 Viewing and Projection Hao Li 1

2 Recall: Affine Transformations Given a point [xyz] > form homogeneous coordinates [xyz1] > The transformed point is [x 0 y 0 z 0 ] > 2

3 Transformation Matrices in OpenGL Transformation matrices in OpenGL are vectors of 16 values (column-major matrices) In glloadmatrixf(glfloat *m); m > =[m 1,m 2,...,m 16 ] > represents Some books transpose all matrices! 3

4 Shear Transformations x-shear scales x proportional to y Leaves y and z values fixed 4

5 Specification via Shear Angle cot( ) =(x 0 x)/y x 0 = x + y cot( ) y 0 = y z 0 = z y [x, y] x 0 x [x 0,y 0 ] y = shear angle x 5

6 Specification via Ratios For example, shear in both x and z direction Leave y fixed Slope for x-shear, for z-shear Solve Yields 6

7 Composing Transformations Let p = Aq, and q = Bs Then p =(AB)s s q p B A AB matrix multiplication 7

8 Composing Transformations Every affine transformation is a composition of rotations, scalings, and translations So, how do we compose these to form an x-shear? Exercise! 8

9 Outline Shear Transformation Camera Positioning Simple Parallel Projections Simple Perspective Projections 9

10 Transform Camera = Transform Scene Camera position is identified with a frame Either move and rotate the objects Or move and rotate the camera Initially, camera at origin, pointing in negative z-direction 10

11 The Look-At Function Convenient way to position camera glulookat(ex, ey, ez, fx, fy, fz, ux, uy, uz); e = eye point f = focus point u e u f u = up vector e f view plane 11

12 OpenGL code void display() { glclear (GL_COLOR_BUFFER_BIT GL_DEPTH_BUFFER_BIT); glmatrixmode (GL_MODELVIEW); glloadidentity(); glulookat (ex, ey, ez, fx, fy, fz, ux, uy, uz); gltranslatef(x, y, z);... renderbunny(); } glutswapbuffers(); 12

13 Implementing the Look-At Function 1. Transform world frame to camera frame - Compose a rotation with translation - 2. Invert to obtain viewing transformation - W = TR W R V = W 1 =(TR) 1 = R 1 T 1 - Derive, then, then R T R 1 T 1 T V 13

14 World Frame to Camera Frame I Camera points in negative z direction n =(f e)/kf ek is unit normal to view plane R [0 0 1] > [n x n y n z ] > Therefore, maps to u e n f view plane 14

15 World Frame to Camera Frame II R [0 1 0] > maps to projection of u onto view plane This projection v equals: = u > n/knk = u > n v 0 = u n v = v 0 /kv 0 k v 0 u e n f view plane 15

16 World Frame to Camera Frame III Set w to be orthogonal to n and v, w = n v, is right-handed [wv n] > v w e n f view plane 16

17 Summary of Rotation glulookat(ex, ey, ez, fx, fy, fz, ux, uy, uz); n =(f e)/kf ek, v =(u (u > n)n)/ku (u > n)nk, w = n v. Rotation must map: - [1 0 0] to - [0 1 0] to w v - [0 0 1] to n 17

18 World Frame to Camera Frame IV Translation of origin to e > =[e x e y e z 1] > 18

19 Camera Frame to Rendering Frame, V = W 1 =(TR) 1 = R 1 T 1 R R 1 = R > is rotation, so T T 1 is translation, so negates displacement 19

20 Putting it Together Calculate V = R 1 T 1 This is different from book [Angel, Ch ] There, u, v, n are right-handed (here: u, v, n ) 20

21 Other Viewing Functions Roll (about z), pitch (about x), yaw (about y) Assignment 2 poses a related problem 21

22 Outline Shear Transformation Camera Positioning Simple Parallel Projections Simple Perspective Projections 22

23 Projection Matrices Recall geometric pipeline Projection takes 3D to 2D Projections are not invertible Projections are described by a 4x4 matrix Homogenous coordinates crucial Parallel and perspective projections 23

24 Parallel Projection Project 3D object to 2D via parallel lines The lines are not necessarily orthogonal to projection plane source:wikipedia 24

25 Parallel Projection Problem: objects far away do not appear smaller Can lead to impossible objects : Penrose stairs source:wikipedia 25

26 Orthographic Projection A special kind of parallel projection: projectors perpendicular to projection plane Simple, but not realistic Used in blueprints (multiview projections) 26

27 Orthographic Projection Matrix Project onto z =0 x p = x y p = y z p =0,, In homogenous coordinates 27

28 Perspective Perspective characterized by foreshortening More distant objects appear smaller Parallel lines appear to converge Rudimentary perspective in cave drawings: Lascaux, France source: Wikipedia 28

29 Discovery of Perspective Foundation in geometry (Euclid) Mural from Pompeii, Italy 29

30 Middle Ages Art in the service of religion Perspective abandoned or forgotten Ottonian manuscript, ca

31 Renaissance Rediscovery, systematic study of perspective Filippo Brunelleschi Florence,

32 Projection (Viewing) in OpenGL Remember: camera is pointing in the negative z direction 32

33 Orthographic Viewing in OpenGL glortho(xmin, xmax, ymin, ymax, near, far) z min = near,z max = far 33

34 Perspective Viewing in OpenGL Two interfaces: glfrustum and gluperspective glfrustum(xmin, xmax, ymin, ymax, near, far); z min = near,z max = far 34

35 Field of View Interface gluperspective(fovy, aspectratio, near, far); near and far as before aspectratio = w/h Fovy specifies field of view as height ( y) angle 35

36 OpenGL code void reshape(int x, int y) { glviewport(0, 0, x, y); glmatrixmode(gl_projection); glloadidentity(); gluperspective(60.0, 1.0 * x / y, 0.01, 10.0); } glmatrixmode(gl_modelview); 36

37 Perspective Viewing Mathematically d = focal length y/z = y p /d so y p = y/(z/d) =yd/z Note that y p is non-linear in the depth z! 37

38 Exploiting the 4 th Dimension Perspective projection is not affine: has no solution for M Idea: exploit homogeneous coordinates for arbitrary w 6= 0 38

39 Perspective Projection Matrix Use multiple of point Solve with 39

40 Projection Algorithm [xyz] > Input: 3D point to project Form [xyz1] > Multiply with ; obtaining M Perform perspective division:,, X/W Y/W Z/W [xyz1] > [XY ZW] > Output: [X/W, Y/W, Z/W ] > (last coordinate will be d ) 40

41 Perspective Division [X Y ZW] > [X/W, Y/W, Z/W, 1] > Normalize to Perform perspective division after projection Projection in OpenGL is more complex (includes clipping) 41

42 Thanks! 42