GIS编程

开发平台：
Visual C++

tess_demo.c：源码内容
							
/* tess_demo.c */
/* A demo of the GLU polygon tesselation functions written by Bogdan Sikorski. */
/* Bug fixes by Mark Kilgard 11/8/96 */
#include <GL/glut.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
/* Win32 calling conventions. */
#ifndef CALLBACK
#define CALLBACK
#endif
#define MAX_POINTS 200
#define MAX_CONTOURS 50
#define HALF_WIDTH 5     /* half of the width of a grid box */
int menu;
typedef enum {
  QUIT, TESSELATE, CLEAR
} menu_entries;
typedef enum {
  DEFINE, TESSELATED
} mode_type;
struct {
  GLint p[MAX_POINTS][2];
  GLuint point_cnt;
} contours[MAX_CONTOURS];
int contour_cnt;
GLsizei width, height;
mode_type mode;
#ifdef GLU_VERSION_1_2
typedef struct listmem {
   GLint *data;
   struct listmem *next;
} ListMem;
ListMem root;
#endif
void CALLBACK
my_error(GLenum err)
{
  int len, i;
  const char *str;
  glColor3f(0.9, 0.9, 0.9);
  glRasterPos2i(5, 5);
  str = (char *) gluErrorString(err);
  len = (int) strlen(str);
  for (i = 0; i < len; i++)
    glutBitmapCharacter(GLUT_BITMAP_9_BY_15, str[i]);
}
#ifdef GLU_VERSION_1_2
/* ARGSUSED */
void CALLBACK
myCombine(GLdouble coords[3], void *d[4], GLfloat w[4], void **dataOut)
{
   ListMem *ptr = &root;
   while (ptr->next) ptr = ptr->next;
   ptr->next = malloc(sizeof(ListMem));
   ptr = ptr->next;
   ptr->next = 0;
   ptr->data = malloc(3 * sizeof(GLint));
   ptr->data[0] = (GLint) coords[0];
   ptr->data[1] = (GLint) coords[1];
   ptr->data[2] = (GLint) coords[2];
   *dataOut = ptr->data;
}
#endif
void 
set_screen_wh(GLsizei w, GLsizei h)
{
  width = w;
  height = h;
}
void 
tesse(void)
{
  static GLUtriangulatorObj *tobj = NULL;
  GLdouble data[3];
  int i, j, point_cnt;
  if (tobj == NULL) {
    tobj = gluNewTess();
  }
  if (tobj != NULL) {
    glClear(GL_COLOR_BUFFER_BIT);
    glColor3f(0.7, 0.5, 0.0);
    gluTessCallback(tobj, GLU_BEGIN, (void (CALLBACK*)())glBegin);
    gluTessCallback(tobj, GLU_END, (void (CALLBACK*)())glEnd);
    gluTessCallback(tobj, GLU_ERROR, (void (CALLBACK*)()) my_error);
    gluTessCallback(tobj, GLU_VERTEX, (void (CALLBACK*)()) glVertex2iv);
#ifdef GLU_VERSION_1_2
    root.data = 0;
    root.next = 0;
    gluTessCallback(tobj, GLU_TESS_COMBINE, (void (CALLBACK*)()) myCombine);
#endif
    gluBeginPolygon(tobj);
    for (j = 0; j <= contour_cnt; j++) {
      point_cnt = contours[j].point_cnt;
      gluNextContour(tobj, GLU_UNKNOWN);
      for (i = 0; i < point_cnt; i++) {
        data[0] = (GLdouble) (contours[j].p[i][0]);
        data[1] = (GLdouble) (contours[j].p[i][1]);
        data[2] = 0.0;
        gluTessVertex(tobj, data, contours[j].p[i]);
      }
    }
    gluEndPolygon(tobj);
    mode = TESSELATED;
  }
#ifdef GLU_VERSION_1_2
  {
    ListMem *fptr = root.next;
    ListMem *tmp;
    while (fptr) {
      tmp = fptr->next;
      free(fptr->data);
      free(fptr);
      fptr = tmp;
    }
  }
#endif
}
void 
left_down(int x1, int y1)
{
  GLint P[2];
  GLuint point_cnt;
  /* translate GLUT into GL coordinates */
  P[0] = x1;
  P[1] = height - y1;
  point_cnt = contours[contour_cnt].point_cnt;
  contours[contour_cnt].p[point_cnt][0] = P[0];
  contours[contour_cnt].p[point_cnt][1] = P[1];
  glBegin(GL_LINES);
  if (point_cnt) {
    glVertex2iv(contours[contour_cnt].p[point_cnt - 1]);
    glVertex2iv(P);
  } else {
    glVertex2iv(P);
    glVertex2iv(P);
  }
  glEnd();
  glFlush();
  ++(contours[contour_cnt].point_cnt);
  if (contours[contour_cnt].point_cnt >= MAX_POINTS) {
    printf("Too many points specified.n");
    exit(1);
  }
}
/* ARGSUSED */
void 
middle_down(int x1, int y1)
{
  GLuint point_cnt;
  point_cnt = contours[contour_cnt].point_cnt;
  if (point_cnt > 2) {
    glBegin(GL_LINES);
    glVertex2iv(contours[contour_cnt].p[0]);
    glVertex2iv(contours[contour_cnt].p[point_cnt - 1]);
    contours[contour_cnt].p[point_cnt][0] = -1;
    glEnd();
    glFlush();
    contour_cnt++;
    if (contour_cnt >= MAX_CONTOURS) {
      printf("Too many contours specified.n");
      exit(1);
    }
    contours[contour_cnt].point_cnt = 0;
  }
}
void 
mouse_clicked(int button, int state, int x, int y)
{
  x += HALF_WIDTH;
  y += HALF_WIDTH;
  x -= x % (2*HALF_WIDTH);
  y -= y % (2*HALF_WIDTH);
  switch (button) {
  case GLUT_LEFT_BUTTON:
    if (state == GLUT_DOWN)
      left_down(x, y);
    break;
  case GLUT_MIDDLE_BUTTON:
    if (state == GLUT_DOWN)
      middle_down(x, y);
    break;
  }
}
void 
display(void)
{
  GLint i, j;
  int point_cnt;
  glClear(GL_COLOR_BUFFER_BIT);
  switch (mode) {
  case DEFINE:
    /* draw grid */
    glColor3f(0.6, 0.5, 0.5);
    glBegin(GL_LINES);
    for (i = 0; i < width; i += (2*HALF_WIDTH)) {
        glVertex2i(i, height);
        glVertex2i(i, 0);
    }
    for (j = 0; j < height; j += (2*HALF_WIDTH)) {
        glVertex2i(0, j);
        glVertex2i(width, j);
    }
    glColor3f(1.0, 1.0, 0.0);
    for (i = 0; i <= contour_cnt; i++) {
      point_cnt = contours[i].point_cnt;
      glBegin(GL_LINES);
      switch (point_cnt) {
      case 0:
        break;
      case 1:
        glVertex2iv(contours[i].p[0]);
        glVertex2iv(contours[i].p[0]);
        break;
      case 2:
        glVertex2iv(contours[i].p[0]);
        glVertex2iv(contours[i].p[1]);
        break;
      default:
        --point_cnt;
        for (j = 0; j < point_cnt; j++) {
          glVertex2iv(contours[i].p[j]);
          glVertex2iv(contours[i].p[j + 1]);
        }
        if (contours[i].p[j + 1][0] == -1) {
          glVertex2iv(contours[i].p[0]);
          glVertex2iv(contours[i].p[j]);
        }
        break;
      }
      glEnd();
    }
    glFlush();
    break;
  case TESSELATED:
    /* draw lines */
    tesse();
    break;
  }
  glColor3f(1.0, 1.0, 0.0);
}
void 
clear(void)
{
  contour_cnt = 0;
  contours[0].point_cnt = 0;
  glutMouseFunc(mouse_clicked);
  mode = DEFINE;
  display();
}
void 
quit(void)
{
  exit(0);
}
void 
menu_selected(int entry)
{
  switch (entry) {
  case CLEAR:
    clear();
    break;
  case TESSELATE:
    tesse();
    break;
  case QUIT:
    quit();
    break;
  }
}
/* ARGSUSED1 */
void 
key_pressed(unsigned char key, int x, int y)
{
  switch (key) {
  case 't':
  case 'T':
    tesse();
    glFlush();
    break;
  case 'q':
  case 'Q':
    quit();
    break;
  case 'c':
  case 'C':
    clear();
    break;
  }
}
void 
myinit(void)
{
  /* clear background to gray */
  glClearColor(0.4, 0.4, 0.4, 0.0);
  glShadeModel(GL_FLAT);
  menu = glutCreateMenu(menu_selected);
  glutAddMenuEntry("clear", CLEAR);
  glutAddMenuEntry("tesselate", TESSELATE);
  glutAddMenuEntry("quit", QUIT);
  glutAttachMenu(GLUT_RIGHT_BUTTON);
  glutMouseFunc(mouse_clicked);
  glutKeyboardFunc(key_pressed);
  contour_cnt = 0;
  glPolygonMode(GL_FRONT, GL_FILL);
  mode = DEFINE;
}
static void 
reshape(GLsizei w, GLsizei h)
{
  glViewport(0, 0, w, h);
  glMatrixMode(GL_PROJECTION);
  glLoadIdentity();
  glOrtho(0.0, (GLdouble) w, 0.0, (GLdouble) h, -1.0, 1.0);
  glMatrixMode(GL_MODELVIEW);
  glLoadIdentity();
  set_screen_wh(w, h);
}
static void 
usage(void)
{
  printf("Use left mouse button to place vertices.n");
  printf("Press middle mouse button when done.n");
  printf("Select tesselate from the pop-up menu.n");
}
/*  Main Loop
 *  Open window with initial window size, title bar, 
 *  RGBA display mode, and handle input events.
 */
int 
main(int argc, char **argv)
{
  usage();
  glutInitWindowSize(400, 400);
  glutInit(&argc, argv);
  glutInitDisplayMode(GLUT_SINGLE | GLUT_RGB);
  glutCreateWindow(argv[0]);
  myinit();
  glutDisplayFunc(display);
  glutReshapeFunc(reshape);
  glutMainLoop();
  return 0;             /* ANSI C requires main to return int. */
}

						