开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > GIS编程 > 查看源码
smoke.c
资源名称：glut23-src.zip　[点击查看]
上传用户：xk288cn
上传日期：2007-05-28
资源大小：4876k
文件大小：10k
源码类别：

GIS编程

开发平台：

Visual C++

smoke.c：源码内容
  1. #include <math.h>
  2. #include <stdio.h>
  3. #include <stdlib.h>
  4. #include <GL/glut.h>
  5. #include "texture.h"
  7. /* Some <math.h> files do not define M_PI... */
  8. #ifndef M_PI
  9. #define M_PI 3.14159265358979323846
  10. #endif
  12. #ifndef __sgi
  13. /* Most math.h's do not define float versions of the trig functions. */
  14. #define sinf sin
  15. #define cosf cos
  16. #define atan2f atan2
  17. #endif
  19. static int texture = 1;
  20. static float rot = 0;
  21. static float opacity = 1.0;
  22. static float intensity = 1.0;
  23. static float size = .001, delta = 0;
  24. static float scale = 1.;
  25. static float transx, transy, rotx, roty;
  26. static int ox = -1, oy = -1;
  27. static int mot = 0;
  28. #define PAN 1
  29. #define ROT 2
  31. void
  32. pan(int x, int y) {
  33.     transx +=  (x-ox)/500.;
  34.     transy -= (y-oy)/500.;
  35.     ox = x; oy = y;
  36.     glutPostRedisplay();
  37. }
  39. void
  40. rotate(int x, int y) {
  41.     rotx += x-ox;
  42.     if (rotx > 360.) rotx -= 360.;
  43.     else if (rotx < -360.) rotx += 360.;
  44.     roty += y-oy;
  45.     if (roty > 360.) roty -= 360.;
  46.     else if (roty < -360.) roty += 360.;
  47.     ox = x; oy = y;
  48.     glutPostRedisplay();
  49. }
  51. void
  52. motion(int x, int y) {
  53.     if (mot == PAN) pan(x, y);
  54.     else if (mot == ROT) rotate(x,y);
  55. }
  57. void
  58. mouse(int button, int state, int x, int y) {
  59.     if(state == GLUT_DOWN) {
  60. switch(button) {
  61. case GLUT_LEFT_BUTTON:
  62.     mot = PAN;
  63.     motion(ox = x, oy = y);
  64.     break;
  65. case GLUT_RIGHT_BUTTON:
  66.     mot = ROT;
  67.     motion(ox = x, oy = y);
  68.     break;
  69. case GLUT_MIDDLE_BUTTON:
  70.     break;
  71. }
  72.     } else if (state == GLUT_UP) {
  73. mot = 0;
  74.     }
  75. }
  77. void afunc(void) {
  78.     static int state;
  79.     if (state ^= 1) {
  80. glAlphaFunc(GL_GREATER, .01);
  81. glEnable(GL_ALPHA_TEST);
  82.     } else {
  83. glDisable(GL_ALPHA_TEST);
  84.     }
  85. }
  87. void bfunc(void) {
  88.     static int state;
  89.     if (state ^= 1) {
  90. glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
  91. glEnable(GL_BLEND);
  92.     } else {
  93. glDisable(GL_BLEND);
  94.     }
  95. }
  97. void tfunc(void) {
  98.     texture ^= 1;
  99. }
  101. void up(void) { scale += .1; }
  102. void down(void) { scale -= .1; }
  103. void left(void) { intensity -= .05f; if (intensity < 0.f) intensity = 0.0f; }
  104. void right(void) { intensity += .05f; if (intensity > 1.f) intensity = 1.0f; }
  106. void help(void) {
  107.     printf("Usage: smoke [image]n");
  108.     printf("'h'            - helpn");
  109.     printf("'a'            - toggle alpha testn");
  110.     printf("'b'            - toggle blendn");
  111.     printf("'t'            - toggle texturingn");
  112.     printf("'UP'           - scale upn");
  113.     printf("'DOWN'         - scale downn");
  114.     printf("'LEFT'    - darkenn");
  115.     printf("'RIGHT'    - brightenn");
  116.     printf("left mouse     - pann");
  117.     printf("right mouse    - rotaten");
  118. }
  120. void init(char *filename) {
  121.     static unsigned *image;
  122.     static int width, height, components;
  123.     if (filename) {
  124. image = read_texture(filename, &width, &height, &components);
  125. if (image == NULL) {
  126.     fprintf(stderr, "Error: Can't load image file "%s".n",
  127.     filename);
  128.     exit(EXIT_FAILURE);
  129. } else {
  130.     printf("%d x %d image loadedn", width, height);
  131. }
  132. #if 0
  133. if (components == 1) {
  134.     GLubyte *p = (GLubyte *)image;
  135.     int i;
  136.     for (i = 0; i < width*height; i++) {
  137. p[i*4+3] = p[i*4+0];
  138.     }
  139.     components = 2;
  140. }
  141. #endif
  142. if (components != 2 && components != 4) {
  143.     printf("must be an RGBA or LA imagen");
  144.     exit(EXIT_FAILURE);
  145. }
  146.     } else {
  147. int i, j;
  148. unsigned char *img;
  149. components = 4; width = height = 512;
  150. image = (unsigned *) malloc(width*height*sizeof(unsigned));
  151. img = (unsigned char *)image;
  152. for (j = 0; j < height; j++)
  153.     for (i = 0; i < width; i++) {
  154. int w2 = width/2, h2 = height/2;
  155. if (i & 32)
  156.     img[4*(i+j*width)+0] = 0xff;
  157. else
  158.     img[4*(i+j*width)+1] = 0xff;
  159. if (j&32)
  160.     img[4*(i+j*width)+2] = 0xff;
  161. if ((i-w2)*(i-w2) + (j-h2)*(j-h2) > 64*64 &&
  162.     (i-w2)*(i-w2) + (j-h2)*(j-h2) < 300*300) img[4*(i+j*width)+3] = 0xff;
  163.     }
  165.     }
  166.     glPixelStorei(GL_UNPACK_ALIGNMENT, 1);
  167.     glTexEnvf(GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_MODULATE);
  168.     glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR);
  169.     glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_LINEAR);
  170.     glTexImage2D(GL_TEXTURE_2D, 0, components, width,
  171.                  height, 0, GL_RGBA, GL_UNSIGNED_BYTE,
  172.                  image);
  173.     glEnable(GL_TEXTURE_2D);
  174.     glMatrixMode(GL_PROJECTION);
  175.     glLoadIdentity();
  176.     gluPerspective(50.,1.,.1,20.);
  177.     glMatrixMode(GL_MODELVIEW);
  178.     glLoadIdentity();
  179.     glTranslatef(0.,0.,-5.5);
  180.     glClearColor(.25f, .25f, .75f, .25f);
  182.     glAlphaFunc(GL_GREATER, 0.016);
  183.     glEnable(GL_ALPHA_TEST);
  185.     glEnable(GL_BLEND);
  186.     glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
  187.     glEnable(GL_DEPTH_TEST);
  188.     glEnable(GL_LIGHT0);
  189.     glEnable(GL_NORMALIZE);
  190. }
  192. void
  193. animate(void) {
  194.     if (delta > 8) {
  195. delta = 0.f;
  196. size = 0.f;
  197.         opacity = 1.f;
  198. rot = 0.f;
  199.     }
  200.     size += .02f;
  201.     delta += .03f;
  202.     rot += .9f;
  203.     opacity -= .005f;
  205.     glutPostRedisplay();
  206. }
  208. void
  209. cube(void) {
  210.     glBegin(GL_QUADS);
  211. glNormal3f(0.f, 0.f, -1.f);
  212. glTexCoord2f(0.0, 0.0); glVertex3f(-1.0, -1.0, -1.0);
  213. glTexCoord2f(1.0, 0.0); glVertex3f( 1.0, -1.0, -1.0);
  214. glTexCoord2f(1.0, 1.0); glVertex3f( 1.0,  1.0, -1.0);
  215. glTexCoord2f(0.0, 1.0); glVertex3f(-1.0,  1.0, -1.0);
  217. glNormal3f(0.f, 0.f, 1.f);
  218. glTexCoord2f(0.0, 0.0); glVertex3f(-1.0, -1.0,  1.0);
  219. glTexCoord2f(1.0, 0.0); glVertex3f( 1.0, -1.0,  1.0);
  220. glTexCoord2f(1.0, 1.0); glVertex3f( 1.0,  1.0,  1.0);
  221. glTexCoord2f(0.0, 1.0); glVertex3f(-1.0,  1.0,  1.0);
  223. glNormal3f(0.f, 1.f, 0.f);
  224. glTexCoord2f(0.0, 0.0); glVertex3f(-1.0,  1.0, -1.0);
  225. glTexCoord2f(1.0, 0.0); glVertex3f( 1.0,  1.0, -1.0);
  226. glTexCoord2f(1.0, 1.0); glVertex3f( 1.0,  1.0,  1.0);
  227. glTexCoord2f(0.0, 1.0); glVertex3f(-1.0,  1.0,  1.0);
  229. glNormal3f(0.f, -1.f, 0.f);
  230. glTexCoord2f(0.0, 0.0); glVertex3f(-1.0, -1.0, -1.0);
  231. glTexCoord2f(1.0, 0.0); glVertex3f( 1.0, -1.0, -1.0);
  232. glTexCoord2f(1.0, 1.0); glVertex3f( 1.0, -1.0,  1.0);
  233. glTexCoord2f(0.0, 1.0); glVertex3f(-1.0, -1.0,  1.0);
  235. glNormal3f( 1.f, 0.f, 0.f);
  236. glTexCoord2f(0.0, 0.0); glVertex3f( 1.0, -1.0, -1.0);
  237. glTexCoord2f(1.0, 0.0); glVertex3f( 1.0,  1.0, -1.0);
  238. glTexCoord2f(1.0, 1.0); glVertex3f( 1.0,  1.0,  1.0);
  239. glTexCoord2f(0.0, 1.0); glVertex3f( 1.0, -1.0,  1.0);
  241. glNormal3f(-1.f, 0.f, 0.f);
  242. glTexCoord2f(0.0, 0.0); glVertex3f(-1.0, -1.0, -1.0);
  243. glTexCoord2f(1.0, 0.0); glVertex3f(-1.0,  1.0, -1.0);
  244. glTexCoord2f(1.0, 1.0); glVertex3f(-1.0,  1.0,  1.0);
  245. glTexCoord2f(0.0, 1.0); glVertex3f(-1.0, -1.0,  1.0);
  246.     glEnd();
  247. }
  248. static void calcMatrix(void);
  250. void display(void) {
  251.     glClear(GL_COLOR_BUFFER_BIT|GL_DEPTH_BUFFER_BIT);
  253.     glLoadIdentity();
  254. #define RAD(x) (((x)*M_PI)/180.)
  255.     gluLookAt(-sinf(RAD(rotx))*5.5,transy,cosf(RAD(rotx))*5.5, 0.,0.,0., 0.,1.,0.);
  257.     /* floor */
  258.     glColor4f(0.f,.2f,0.f,1.f);
  259.     glBegin(GL_POLYGON);
  260. glVertex3f(-4.0, -1.0, -4.0);
  261. glVertex3f( 4.0, -1.0, -4.0);
  262. glVertex3f( 4.0, -1.0,  4.0);
  263. glVertex3f(-4.0, -1.0,  4.0);
  264.     glEnd();
  266.     glEnable(GL_COLOR_MATERIAL);
  267.     glEnable(GL_LIGHTING);
  268.     glColor3f(.3f,.3f,.3f);
  269.     glPushMatrix();
  270.     glTranslatef(-1.f, -1.+.2f, -1.5f);
  271.     glScalef(.2f,.2f, .2f);
  272.     cube();
  273.     glDisable(GL_LIGHTING);
  274.     glPopMatrix();
  276.     glPushMatrix();
  277.     glTranslatef(delta/2.f-1.f, delta-1.f, -1.5f);
  278.     calcMatrix();
  279.     glScalef(size,size,1.);
  280.     if (texture) glEnable(GL_TEXTURE_2D);
  281.     glColor4f(intensity, intensity, intensity, opacity);
  282.     glRotatef(rot, 0., 0., 1.);
  283.     glDepthMask(0);
  284.     glBegin(GL_POLYGON);
  285. glTexCoord2f(0.0, 0.0); glVertex2f(-1.0, -1.0);
  286. glTexCoord2f(1.0, 0.0); glVertex2f(1.0, -1.0);
  287. glTexCoord2f(1.0, 1.0); glVertex2f(1.0, 1.0);
  288. glTexCoord2f(0.0, 1.0); glVertex2f(-1.0, 1.0);
  289.     glEnd();
  290.     glDepthMask(1);
  291.     glPopMatrix();
  292.     glDisable(GL_TEXTURE_2D);
  294.     glutSwapBuffers();
  295. }
  297. void reshape(int w, int h) {
  298.     glViewport(0, 0, w, h);
  299. }
  301. /* ARGSUSED1 */
  302. void
  303. key(unsigned char key, int x, int y) {
  304.     switch(key) {
  305.     case 'a': afunc(); break;
  306.     case 'b': bfunc(); break;
  307.     case 'h': help(); break;
  308.     case 't': tfunc(); break;
  309.     case '33': exit(EXIT_SUCCESS); break;
  310.     default: break;
  311.     }
  312.     glutPostRedisplay();
  313. }
  315. /* ARGSUSED1 */
  316. void
  317. special(int key, int x, int y) {
  318.     switch(key) {
  319.     case GLUT_KEY_UP: up(); break;
  320.     case GLUT_KEY_DOWN: down(); break;
  321.     case GLUT_KEY_LEFT: left(); break;
  322.     case GLUT_KEY_RIGHT:right(); break;
  323.     }
  324. }
  326. int main(int argc, char** argv) {
  327.     glutInitWindowSize(512, 512);
  328.     glutInit(&argc, argv);
  329.     glutInitDisplayMode(GLUT_RGBA|GLUT_DOUBLE|GLUT_DEPTH);
  330.     (void)glutCreateWindow("smoke");
  331.     init(argc == 1 ? "../data/smoke.bw" : argv[1]);
  332.     glutDisplayFunc(display);
  333.     glutKeyboardFunc(key);
  334.     glutSpecialFunc(special);
  335.     glutReshapeFunc(reshape);
  336.     glutMouseFunc(mouse);
  337.     glutMotionFunc(motion);
  338.     glutIdleFunc(animate);
  339.     glutMainLoop();
  340.     return 0;
  341. }
  343. void
  344. printmat(float *m) {
  345.     int i;
  346.     for(i = 0; i < 4; i++) {
  347. printf("%f %f %f %fn", m[4*i+0], m[4*i+1], m[4*i+2], m[4*i+3]);
  348.     }
  349. }
  351. void
  352. buildRot(float theta, float x, float y, float z, float m[16]) {
  353.     float d = x*x + y*y + z*z;
  354.     float ct = cosf(RAD(theta)), st = sinf(RAD(theta));
  356.     /* normalize */
  357.     if (d > 0) {
  358. d = 1/d;
  359. x *= d;
  360. y *= d;
  361. z *= d;
  362.     }
  364.     m[ 0] = 1; m[ 1] = 0; m[ 2] = 0; m[ 3] = 0;
  365.     m[ 4] = 0; m[ 5] = 1; m[ 6] = 0; m[ 7] = 0;
  366.     m[ 8] = 0; m[ 9] = 0; m[10] = 1; m[11] = 0;
  367.     m[12] = 0; m[13] = 0; m[14] = 0; m[15] = 1;
  369.     /* R = uu' + cos(theta)*(I-uu') + sin(theta)*S
  370.      *
  371.      * S =  0  -z   y    u' = (x, y, z)
  372.      *     z   0  -x
  373.      *    -y   x   0
  374.      */
  376.      m[0] = x*x + ct*(1-x*x) + st*0;
  377.      m[4] = x*y + ct*(0-x*y) + st*-z;
  378.      m[8] = x*z + ct*(0-x*z) + st*y;
  380.      m[1] = y*x + ct*(0-y*x) + st*z;
  381.      m[5] = y*y + ct*(1-y*y) + st*0;
  382.      m[9] = y*z + ct*(0-y*z) + st*-x;
  384.      m[2] = z*x + ct*(0-z*x) + st*-y;
  385.      m[6] = z*y + ct*(0-z*y) + st*x;
  386.      m[10]= z*z + ct*(1-z*z) + st*0;
  387. }
  389. static void
  390. calcMatrix(void) {
  391.     float mat[16];
  393.     glGetFloatv(GL_MODELVIEW_MATRIX, mat);
  395.     buildRot(-180*atan2f(mat[8], mat[10])/M_PI, 0, 1, 0, mat);
  396.     glMultMatrixf(mat);
  397. }