depth_of_field.cpp

Go to the documentation of this file.

/*
 * Copyright (c) 1993-1997, Silicon Graphics, Inc.
 * ALL RIGHTS RESERVED 
 * Permission to use, copy, modify, and distribute this software for 
 * any purpose and without fee is hereby granted, provided that the above
 * copyright notice appear in all copies and that both the copyright notice
 * and this permission notice appear in supporting documentation, and that 
 * the name of Silicon Graphics, Inc. not be used in advertising
 * or publicity pertaining to distribution of the software without specific,
 * written prior permission. 
 *
 * THE MATERIAL EMBODIED ON THIS SOFTWARE IS PROVIDED TO YOU "AS-IS"
 * AND WITHOUT WARRANTY OF ANY KIND, EXPRESS, IMPLIED OR OTHERWISE,
 * INCLUDING WITHOUT LIMITATION, ANY WARRANTY OF MERCHANTABILITY OR
 * FITNESS FOR A PARTICULAR PURPOSE.  IN NO EVENT SHALL SILICON
 * GRAPHICS, INC.  BE LIABLE TO YOU OR ANYONE ELSE FOR ANY DIRECT,
 * SPECIAL, INCIDENTAL, INDIRECT OR CONSEQUENTIAL DAMAGES OF ANY
 * KIND, OR ANY DAMAGES WHATSOEVER, INCLUDING WITHOUT LIMITATION,
 * LOSS OF PROFIT, LOSS OF USE, SAVINGS OR REVENUE, OR THE CLAIMS OF
 * THIRD PARTIES, WHETHER OR NOT SILICON GRAPHICS, INC.  HAS BEEN
 * ADVISED OF THE POSSIBILITY OF SUCH LOSS, HOWEVER CAUSED AND ON
 * ANY THEORY OF LIABILITY, ARISING OUT OF OR IN CONNECTION WITH THE
 * POSSESSION, USE OR PERFORMANCE OF THIS SOFTWARE.
 * 
 * US Government Users Restricted Rights 
 * Use, duplication, or disclosure by the Government is subject to
 * restrictions set forth in FAR 52.227.19(c)(2) or subparagraph
 * (c)(1)(ii) of the Rights in Technical Data and Computer Software
 * clause at DFARS 252.227-7013 and/or in similar or successor
 * clauses in the FAR or the DOD or NASA FAR Supplement.
 * Unpublished-- rights reserved under the copyright laws of the
 * United States.  Contractor/manufacturer is Silicon Graphics,
 * Inc., 2011 N.  Shoreline Blvd., Mountain View, CA 94039-7311.
 *
 * OpenGL(R) is a registered trademark of Silicon Graphics, Inc.
 */

/*
 *  dof.c
 *  This program demonstrates use of the accumulation buffer to
 *  create an out-of-focus depth-of-field effect.  The teapots
 *  are drawn several times into the accumulation buffer.  The
 *  viewing volume is jittered, except at the focal point, where
 *  the viewing volume is at the same position, each time.  In
 *  this case, the gold teapot remains in focus.
 */
 
#include <GL/glut.h>
#include <GL/glu.h>
#include <stdlib.h>
#include <math.h>
#include "jitter.h"

#ifdef WIN32
#define near zNear
#define far zFar
#endif

#define PI_ 3.14159265358979323846


// File a little bit modified in order to give the user ability to 
// change the number of jitter point during runtime and focus distance.
// After some change scene first draw sharp and then jittered in accumulation
// buffer in order to get depth of field effect.

#include <stdio.h>


//focus
float z = 5.0;

bool changed = false;



//rotation variables
float AngleX=0, AngleY=0;
bool RotateUp = false, RotateDown = false, RotateLeft = false, RotateRight = false;

float CameraX = 0, CameraY = 0, CameraZ = 5.0f;

int LastTime = 0;
int TimerDelay;
float FPS = 0;

//helper function from Bastian
void glutString( char *pstr )
{
  while (*pstr!=char(0)) {
    glutBitmapCharacter( GLUT_BITMAP_8_BY_13, int( *pstr ));
    pstr++;
  }
}

//
enum { JITTER_2 = 0, JITTER_3, JITTER_4, JITTER_8, JITTER_15, JITTER_24, JITTER_66, JITTER_N};

int njitters[JITTER_N] = { 2, 3, 4, 8, 15, 24, 66 };

jitter_point* jitters[JITTER_N];

int current_jitter = JITTER_8;

GLuint teapotList;

/* accFrustum()
 * The first 6 arguments are identical to the glFrustum() call.
 *  
 * pixdx and pixdy are anti-alias jitter in pixels. 
 * Set both equal to 0.0 for no anti-alias jitter.
 * eyedx and eyedy are depth-of field jitter in pixels. 
 * Set both equal to 0.0 for no depth of field effects.
 *
 * focus is distance from eye to plane in focus. 
 * focus must be greater than, but not equal to 0.0.
 *
 * Note that accFrustum() calls glTranslatef().  You will 
 * probably want to insure that your ModelView matrix has been 
 * initialized to identity before calling accFrustum().
 */
void accFrustum(GLdouble left, GLdouble right, GLdouble bottom, 
   GLdouble top, GLdouble near, GLdouble far, GLdouble pixdx, 
   GLdouble pixdy, GLdouble eyedx, GLdouble eyedy, GLdouble focus)
{
   GLdouble xwsize, ywsize; 
   GLdouble dx, dy;
   GLint viewport[4];

   glGetIntegerv (GL_VIEWPORT, viewport);
        
   xwsize = right - left;
   ywsize = top - bottom;
        
   dx = -(pixdx*xwsize/(GLdouble) viewport[2] + eyedx*near/focus);
   dy = -(pixdy*ywsize/(GLdouble) viewport[3] + eyedy*near/focus);
        
   glMatrixMode(GL_PROJECTION);
   glLoadIdentity();
   glFrustum (left + dx, right + dx, bottom + dy, top + dy, near, far);
   glMatrixMode(GL_MODELVIEW);
   glLoadIdentity();
   glTranslatef (-eyedx, -eyedy, 0.0);
}

/*  accPerspective()
 * 
 *  The first 4 arguments are identical to the gluPerspective() call.
 *  pixdx and pixdy are anti-alias jitter in pixels. 
 *  Set both equal to 0.0 for no anti-alias jitter.
 *  eyedx and eyedy are depth-of field jitter in pixels. 
 *  Set both equal to 0.0 for no depth of field effects.
 *
 *  focus is distance from eye to plane in focus. 
 *  focus must be greater than, but not equal to 0.0.
 *
 *  Note that accPerspective() calls accFrustum().
 */
void accPerspective(GLdouble fovy, GLdouble aspect, 
   GLdouble near, GLdouble far, GLdouble pixdx, GLdouble pixdy, 
   GLdouble eyedx, GLdouble eyedy, GLdouble focus)
{
   GLdouble fov2,left,right,bottom,top;

   fov2 = ((fovy*PI_) / 180.0) / 2.0;

   top = near / (cos(fov2) / sin(fov2));
   bottom = -top;

   right = top * aspect;
   left = -right;

   accFrustum (left, right, bottom, top, near, far,
               pixdx, pixdy, eyedx, eyedy, focus);
}

void init(void)
{
jitters[0] = j2;
jitters[1] = j3;
jitters[2] = j4;
jitters[3] = j8;
jitters[4] = j15;
jitters[5] = j24;
jitters[6] = j66;

GLenum performance = GL_FASTEST; //GL_NICEST
glHint(GL_POINT_SMOOTH_HINT, performance );
glHint(GL_LINE_SMOOTH_HINT, performance );
glHint(GL_POLYGON_SMOOTH_HINT, performance );
glHint(GL_FOG_HINT, performance );
glHint(GL_PERSPECTIVE_CORRECTION_HINT, performance );

/*
  glEnable(GL_DEPTH_TEST);
  glDepthFunc(GL_LESS);
  glShadeModel(GL_SMOOTH);
  glEnable(GL_LIGHTING);
  glDisable( GL_CULL_FACE );
  glEnable( GL_NORMALIZE );
  */
   GLfloat ambient[] = { 0.0, 0.0, 0.0, 1.0 };
   GLfloat diffuse[] = { 1.0, 1.0, 1.0, 1.0 };
   GLfloat specular[] = { 1.0, 1.0, 1.0, 1.0 };
   GLfloat position[] = { 0.0, 3.0, 3.0, 0.0 };
    
   GLfloat lmodel_ambient[] = { 0.2, 0.2, 0.2, 1.0 };
   GLfloat local_view[] = { 0.0 };

   glLightfv(GL_LIGHT0, GL_AMBIENT, ambient);
   glLightfv(GL_LIGHT0, GL_DIFFUSE, diffuse);
   glLightfv(GL_LIGHT0, GL_POSITION, position);
    
   glLightModelfv(GL_LIGHT_MODEL_AMBIENT, lmodel_ambient);
   glLightModelfv(GL_LIGHT_MODEL_LOCAL_VIEWER, local_view);

   glFrontFace (GL_CW);
   glEnable(GL_LIGHTING);
   glEnable(GL_LIGHT0);
   glEnable(GL_AUTO_NORMAL);
   glEnable(GL_NORMALIZE);
   glEnable(GL_DEPTH_TEST);

   glClearColor(0.0, 0.0, 0.0, 0.0);
   glClearAccum(0.0, 0.0, 0.0, 0.0); 
/*  make teapot display list */
   teapotList = glGenLists(1);
   glNewList (teapotList, GL_COMPILE);
   glutSolidTeapot (0.5);
   glEndList ();
}

void renderTeapot (GLfloat x, GLfloat y, GLfloat z, 
   GLfloat ambr, GLfloat ambg, GLfloat ambb, 
   GLfloat difr, GLfloat difg, GLfloat difb, 
   GLfloat specr, GLfloat specg, GLfloat specb, GLfloat shine)
{
   GLfloat mat[4];

   glPushMatrix();
   glTranslatef (x, y, z);
   mat[0] = ambr; mat[1] = ambg; mat[2] = ambb; mat[3] = 1.0;   
   glMaterialfv (GL_FRONT, GL_AMBIENT, mat);
   mat[0] = difr; mat[1] = difg; mat[2] = difb; 
   glMaterialfv (GL_FRONT, GL_DIFFUSE, mat);
   mat[0] = specr; mat[1] = specg; mat[2] = specb;
   glMaterialfv (GL_FRONT, GL_SPECULAR, mat);
   glMaterialf (GL_FRONT, GL_SHININESS, shine*128.0);
   glCallList(teapotList);
   glPopMatrix();
}

/*  display() draws 5 teapots into the accumulation buffer 
 *  several times; each time with a jittered perspective.
 *  The focal point is at z = 5.0, so the gold teapot will 
 *  stay in focus.  The amount of jitter is adjusted by the
 *  magnitude of the accPerspective() jitter; in this example, 0.33.
 *  In this example, the teapots are drawn 8 times.  See jitter.h
 */
void display(void)
{
   int jitter;
   GLint viewport[4];

//    glColor3f( 1.0f, 1.0f, 1.0f );

   glClearColor(0.0, 0.0, 0.0, 0.0);
   glClearAccum(0.0, 0.0, 0.0, 0.0); 

   glGetIntegerv (GL_VIEWPORT, viewport);
  // glClear(GL_ACCUM_BUFFER_BIT);
//glClear(GL_COLOR_BUFFER_BIT);
      glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT | GL_ACCUM_BUFFER_BIT);

if (changed) 
        {
        changed = false;
      glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
      accPerspective (45.0, 
         (GLdouble) viewport[2]/(GLdouble) viewport[3], 
         1.0, 15.0, 0.0, 0.0,
         0, 0, z);
      renderTeapot (-1.1, -0.5, -4.5, 0.1745, 0.01175, 
                    0.01175, 0.61424, 0.04136, 0.04136, 
                    0.727811, 0.626959, 0.626959, 0.6);
      renderTeapot (-0.5, -0.5, -5.0, 0.24725, 0.1995, 
                    0.0745, 0.75164, 0.60648, 0.22648, 
                    0.628281, 0.555802, 0.366065, 0.4);
      renderTeapot (0.2, -0.5, -5.5, 0.19225, 0.19225, 
                    0.19225, 0.50754, 0.50754, 0.50754, 
                    0.508273, 0.508273, 0.508273, 0.4);
      renderTeapot (1.0, -0.5, -6.0, 0.0215, 0.1745, 0.0215, 
                    0.07568, 0.61424, 0.07568, 0.633, 
                    0.727811, 0.633, 0.6);
      renderTeapot (1.8, -0.5, -6.5, 0.0, 0.1, 0.06, 0.0, 
                    0.50980392, 0.50980392, 0.50196078, 
                    0.50196078, 0.50196078, .25);
       glAccum (GL_ACCUM, 1.0);
       glAccum (GL_RETURN, 1.0);
       
       
       //and now making depth of field
        }
        else
        {//if there are no changes in camera position,
        // then accumulate to the buffer

   jitter_point *tmp;
   for (jitter = 0; jitter < njitters[current_jitter]; jitter++) {
      glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
      tmp = (jitters[current_jitter] + jitter);
      accPerspective (45.0, 
         (GLdouble) viewport[2]/(GLdouble) viewport[3], 
         1.0, 15.0, 0.0, 0.0,
         .33 * tmp->x, .33 * tmp -> y, z);
//        glMatrixMode(GL_PROJECTION);
        //  glTranslated(CameraX, CameraY, CameraZ);
//      glRotated( AngleY, 0, 1, 0);
          glMatrixMode(GL_MODELVIEW);
//        glTranslated(-CameraX, -CameraY, -CameraZ);
                glPushMatrix();
/*      ruby, gold, silver, emerald, and cyan teapots   */
      renderTeapot (-1.1, -0.5, -4.5, 0.1745, 0.01175, 
                    0.01175, 0.61424, 0.04136, 0.04136, 
                    0.727811, 0.626959, 0.626959, 0.6);
      renderTeapot (-0.5, -0.5, -5.0, 0.24725, 0.1995, 
                    0.0745, 0.75164, 0.60648, 0.22648, 
                    0.628281, 0.555802, 0.366065, 0.4);
      renderTeapot (0.2, -0.5, -5.5, 0.19225, 0.19225, 
                    0.19225, 0.50754, 0.50754, 0.50754, 
                    0.508273, 0.508273, 0.508273, 0.4);
      renderTeapot (1.0, -0.5, -6.0, 0.0215, 0.1745, 0.0215, 
                    0.07568, 0.61424, 0.07568, 0.633, 
                    0.727811, 0.633, 0.6);
      renderTeapot (1.8, -0.5, -6.5, 0.0, 0.1, 0.06, 0.0, 
                    0.50980392, 0.50980392, 0.50196078, 
                    0.50196078, 0.50196078, .25);
                glPopMatrix();
      glAccum (GL_ACCUM, (1.0f/njitters[current_jitter]) ); //0.125
   }//end for
   glAccum (GL_RETURN, 1.0);
   }//end if 
   
    
/*       // Display framerate
glPushMatrix();
    char buf[100];
    sprintf( buf, "READY %0.2f fps,  X=%f, Y=%f, Z=%f, alpha=%f, beta=%f ", FPS, CameraX, CameraY, CameraZ, 
        AngleX, AngleY );
    glMatrixMode( GL_PROJECTION );
    glLoadIdentity();
    glMatrixMode( GL_MODELVIEW );
    glLoadIdentity();
    glDisable( GL_DEPTH_TEST );
    glDisable( GL_LIGHTING );
    glColor3f( 1.0f, 1.0f, 1.0f );
    glRasterPos2f( -1,-1 );
    glutString( buf );
    glEnable( GL_DEPTH_TEST );
    glEnable(GL_LIGHTING);
glPopMatrix();
  */ 
   glutSwapBuffers();
}

void print_jitter_points (int n)
{
jitter_point * tmp;
   for (int jitter = 0; jitter < njitters[n]; jitter++) 
      {
      tmp = (jitters[n] + jitter);
      printf(" coordinates of %d-jitter point: x = %f, y = %f \n", jitter + 1, tmp->x, tmp->y );
      }

}

void reshape(int w, int h)
{
changed = true;
   glViewport(0, 0, (GLsizei) w, (GLsizei) h);
   glutPostRedisplay();
}

void keyboard(unsigned char key, int x, int y)
{
   switch (key) {
      case 27:
      case 'q':
         exit(0);
         break;
          case ' ':
                 glutPostRedisplay();
                 break;
          case 'a':
                RotateLeft = true;
                RotateRight = false;
                break;
          case 'd':
                RotateRight = true;
                RotateLeft = false;
                break;
          case 'n': 
                //increase focus
                z+= 0.1;
                changed = true;
                break;
          case 'm': z-=0.1;
                changed = true;
                break;
         case '1': 
                changed = true; 
                current_jitter = JITTER_2; 
                printf(" From now will use %d jitter points. \n", 2);
                print_jitter_points(JITTER_2);
                break;
         case '2': 
                changed = true; 
                current_jitter = JITTER_3; 
                printf(" From now will use %d jitter points. \n", 3);
                print_jitter_points(JITTER_3);
                break;
         case '3': 
                changed = true; 
                current_jitter = JITTER_4; 
                printf(" From now will use %d jitter points. \n", 4);
                print_jitter_points(JITTER_4);
                break;
         case '4': 
                changed = true; 
                current_jitter = JITTER_8; 
                printf(" From now will use %d jitter points. \n", 8);
                print_jitter_points(JITTER_8);
                break;
         case '5': 
                changed = true; 
                current_jitter = JITTER_15; 
                printf(" From now will use %d jitter points. \n", 15);
                print_jitter_points(JITTER_15);
                break;
         case '6': 
                changed = true; 
                current_jitter = JITTER_24; 
                printf(" From now will use %d jitter points. \n", 24);
                print_jitter_points(JITTER_24);
                break;
         case '7': 
                changed = true; 
                current_jitter = JITTER_66; 
                printf(" From now will use %d jitter points. \n", 66);
                print_jitter_points(JITTER_66);
                break;
   }
}

static void timerCallback(int)
{

int nc = glutGet(GLUT_ELAPSED_TIME);
float fSec = float(nc-LastTime) / 1000.0f;
  // Compute average frame rate from the last twenty frames
  FPS = (19.0f*FPS + 1.0f / fSec) / 20.0f;
  LastTime = nc;
  
  float fCamArc = fSec * 180.0f;
  if (RotateLeft) AngleY -= fCamArc/4;
  if (RotateRight) AngleY += fCamArc/4;
  
  glutPostRedisplay();
  glutTimerFunc(TimerDelay, timerCallback, 0);

}



/*  Main Loop
 *  Be certain you request an accumulation buffer.
 */

void keyboardUp(unsigned char key, int, int)
{
switch (key)
        {
        case 'a':
                RotateLeft = false;
                break;
        case 'd': RotateRight = false; break;
        }

}
 
int main(int argc, char** argv)
{

printf(" Keys: \n");
printf("\tn/m  -- increase/decrease focus distance \n");
printf("\t1..7 -- change number of jitter points   \n");
printf("\tq    -- to Quit.\n");

   glutInit(&argc, argv);
   glutInitDisplayMode (GLUT_DOUBLE | GLUT_RGB
                        | GLUT_ACCUM | GLUT_DEPTH); 
   glutInitWindowSize (400, 400);
   glutInitWindowPosition (100, 100);
   glutCreateWindow (argv[0]);
   init();
   glutReshapeFunc(reshape);
   glutDisplayFunc(display);
   glutKeyboardFunc(keyboard);
   glutKeyboardUpFunc(keyboardUp);
   LastTime = glutGet (GLUT_ELAPSED_TIME );
   glutTimerFunc ( TimerDelay, timerCallback, 0);
   glutMainLoop();
   return 0;
}

---