1 // oursun.hxx -- model earth's sun
3 // Written by Durk Talsma. Originally started October 1997, for distribution
4 // with the FlightGear project. Version 2 was written in August and
5 // September 1998. This code is based upon algorithms and data kindly
6 // provided by Mr. Paul Schlyter. (pausch@saaf.se).
8 // Separated out rendering pieces and converted to ssg by Curt Olson,
10 // This library is free software; you can redistribute it and/or
11 // modify it under the terms of the GNU Library General Public
12 // License as published by the Free Software Foundation; either
13 // version 2 of the License, or (at your option) any later version.
15 // This library is distributed in the hope that it will be useful,
16 // but WITHOUT ANY WARRANTY; without even the implied warranty of
17 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
18 // Library General Public License for more details.
20 // You should have received a copy of the GNU Library General Public
21 // License along with this library; if not, write to the
22 // Free Software Foundation, Inc., 59 Temple Place - Suite 330,
23 // Boston, MA 02111-1307, USA.
29 # include <simgear_config.h>
32 #include <simgear/compiler.h>
40 // define the following to enable a cheesy lens flare effect for the sun
41 // #define FG_TEST_CHEESY_LENS_FLARE
43 #ifdef FG_TEST_CHEESY_LENS_FLARE
44 # include <plib/ssgaLensFlare.h>
47 #include <simgear/screen/colors.hxx>
52 static double sun_exp2_punch_through;
54 // Set up sun rendering call backs
55 static int sgSunOrbPreDraw( ssgEntity *e ) {
56 /* cout << endl << "Sun orb pre draw" << endl << "----------------"
59 ssgLeaf *f = (ssgLeaf *)e;
60 if ( f -> hasState () ) f->getState()->apply() ;
62 glPushAttrib( GL_DEPTH_BUFFER_BIT | GL_FOG_BIT );
63 // cout << "push error = " << glGetError() << endl;
65 glDisable( GL_DEPTH_TEST );
71 static int sgSunOrbPostDraw( ssgEntity *e ) {
72 /* cout << endl << "Sun orb post draw" << endl << "----------------"
76 // cout << "pop error = " << glGetError() << endl;
81 static int sgSunHaloPreDraw( ssgEntity *e ) {
82 /* cout << endl << "Sun halo pre draw" << endl << "----------------"
85 ssgLeaf *f = (ssgLeaf *)e;
86 if ( f -> hasState () ) f->getState()->apply() ;
88 glPushAttrib( GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT | GL_FOG_BIT );
89 // cout << "push error = " << glGetError() << endl;
91 glDisable( GL_DEPTH_TEST );
92 // glDisable( GL_FOG );
93 glFogf (GL_FOG_DENSITY, sun_exp2_punch_through);
94 glBlendFunc ( GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA ) ;
99 static int sgSunHaloPostDraw( ssgEntity *e ) {
100 /* cout << endl << "Sun halo post draw" << endl << "----------------"
104 // cout << "pop error = " << glGetError() << endl;
111 SGSun::SGSun( void ) {
112 prev_sun_angle = -9999.0;
113 visibility = -9999.0;
118 SGSun::~SGSun( void ) {
123 // this might be nice to keep, just as an example of how to generate a
124 // texture on the fly ...
125 static GLuint makeHalo( GLubyte *sun_texbuf, int width ) {
132 // create a texture id
133 #ifdef GL_VERSION_1_1
134 glGenTextures(1, &texid);
135 glBindTexture(GL_TEXTURE_2D, texid);
136 #elif GL_EXT_texture_object
137 glGenTexturesEXT(1, &texid);
138 glBindTextureEXT(GL_TEXTURE_2D, texid);
143 glPixelStorei( GL_UNPACK_ALIGNMENT, 4 );
144 glTexParameteri( GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR );
145 glTexParameteri( GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_LINEAR );
146 glTexEnvi( GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_MODULATE ) ;
148 // create the actual texture contents
149 texSize = width * width;
152 SG_LOG( SG_EVENT, SG_ALERT,
153 "Could not allocate memroy for the sun texture");
159 radius = (double)(width / 2);
163 for ( i = 0; i < width; i++ ) {
164 for ( j = 0; j < width; j++ ) {
165 x = fabs((double)(i - (width / 2)));
166 y = fabs((double)(j - (width / 2)));
167 d = sqrt((x * x) + (y * y));
169 // t is 1.0 at center, 0.0 at edge
170 double t = 1.0 - (d / radius);
172 // inverse square looks nice
173 value = (int)((double) 0xff * (t*t));
186 /* glTexImage2D( GL_TEXTURE_2D,
191 GL_RGBA, GL_UNSIGNED_BYTE,
198 #define RGB 3 // 3 bytes of color info per pixel
199 #define RGBA 4 // 4 bytes of color+alpha info
200 void my_glWritePPMFile(const char *filename, GLubyte *buffer, int win_width, int win_height, int mode)
203 unsigned char *ibuffer;
205 int pixelSize = mode==GL_RGBA?4:3;
207 ibuffer = (unsigned char *) malloc(win_width*win_height*RGB);
209 fp = fopen(filename, "wb");
210 fprintf(fp, "P6\n# CREATOR: glReadPixel()\n%d %d\n%d\n",
211 win_width, win_height, UCHAR_MAX);
213 for (i = 0; i < win_height; i++) {
214 for (j = 0; j < win_width; j++) {
215 for (k = 0; k < RGB; k++) {
216 ibuffer[q++] = (unsigned char)
217 *(buffer + (pixelSize*((win_height-1-i)*win_width+j)+k));
222 // *(buffer + (pixelSize*((win_height-1-i)*win_width+j)+k));
224 fwrite(ibuffer, sizeof(unsigned char), RGB*win_width*win_height, fp);
228 printf("wrote file (%d x %d pixels, %d bytes)\n",
229 win_width, win_height, RGB*win_width*win_height);
234 // initialize the sun object and connect it into our scene graph root
235 ssgBranch * SGSun::build( SGPath path, double sun_size ) {
237 // set up the orb state
238 orb_state = new ssgSimpleState();
239 orb_state->setShadeModel( GL_SMOOTH );
240 orb_state->disable( GL_LIGHTING );
241 // orb_state->enable( GL_LIGHTING );
242 orb_state->disable( GL_CULL_FACE );
243 orb_state->disable( GL_TEXTURE_2D );
244 orb_state->enable( GL_COLOR_MATERIAL );
245 orb_state->setColourMaterial( GL_AMBIENT_AND_DIFFUSE );
246 orb_state->setMaterial( GL_EMISSION, 0, 0, 0, 1 );
247 orb_state->setMaterial( GL_SPECULAR, 0, 0, 0, 1 );
248 orb_state->disable( GL_BLEND );
249 orb_state->disable( GL_ALPHA_TEST );
251 cl = new ssgColourArray( 1 );
253 sgSetVec4( color, 1.0, 1.0, 1.0, 1.0 );
256 ssgBranch *orb = ssgMakeSphere( orb_state, cl, sun_size, 10, 10,
257 sgSunOrbPreDraw, sgSunOrbPostDraw );
259 // force a repaint of the sun colors with arbitrary defaults
263 // sun_texbuf = new GLubyte[64*64*3];
264 // sun_texid = makeHalo( sun_texbuf, 64 );
265 // my_glWritePPMFile("sunhalo.ppm", sun_texbuf, 64, 64, RGB);
267 // set up the halo state
268 path.append( "halo.rgba" );
269 halo_state = new ssgSimpleState();
270 halo_state->setTexture( (char *)path.c_str() );
271 halo_state->enable( GL_TEXTURE_2D );
272 halo_state->disable( GL_LIGHTING );
273 // halo_state->enable( GL_LIGHTING );
274 halo_state->setShadeModel( GL_SMOOTH );
275 halo_state->disable( GL_CULL_FACE );
276 halo_state->enable( GL_COLOR_MATERIAL );
277 halo_state->setColourMaterial( GL_AMBIENT_AND_DIFFUSE );
278 halo_state->setMaterial( GL_EMISSION, 0, 0, 0, 1 );
279 halo_state->setMaterial( GL_SPECULAR, 0, 0, 0, 1 );
280 halo_state->enable( GL_ALPHA_TEST );
281 halo_state->setAlphaClamp(0.01);
282 halo_state->enable ( GL_BLEND ) ;
284 // Build ssg structure
285 double size = sun_size * 10.0;
287 halo_vl = new ssgVertexArray;
288 sgSetVec3( v3, -size, 0.0, -size );
290 sgSetVec3( v3, size, 0.0, -size );
292 sgSetVec3( v3, -size, 0.0, size );
294 sgSetVec3( v3, size, 0.0, size );
298 halo_tl = new ssgTexCoordArray;
299 sgSetVec2( v2, 0.0f, 0.0f );
301 sgSetVec2( v2, 1.0, 0.0 );
303 sgSetVec2( v2, 0.0, 1.0 );
305 sgSetVec2( v2, 1.0, 1.0 );
309 new ssgVtxTable ( GL_TRIANGLE_STRIP, halo_vl, NULL, halo_tl, cl );
310 halo->setState( halo_state );
312 // build the ssg scene graph sub tree for the sky and connected
313 // into the provide scene graph branch
314 sun_transform = new ssgTransform;
316 halo->setCallback( SSG_CALLBACK_PREDRAW, sgSunHaloPreDraw );
317 halo->setCallback( SSG_CALLBACK_POSTDRAW, sgSunHaloPostDraw );
318 sun_transform->addKid( halo );
319 sun_transform->addKid( orb );
321 #ifdef FG_TEST_CHEESY_LENS_FLARE
323 sun_transform->addKid( new ssgaLensFlare );
326 return sun_transform;
330 // repaint the sun colors based on current value of sun_angle in
331 // degrees relative to verticle
332 // 0 degrees = high noon
333 // 90 degrees = sun rise/set
334 // 180 degrees = darkest midnight
335 bool SGSun::repaint( double sun_angle, double new_visibility ) {
336 if ( visibility != new_visibility ) {
337 visibility = new_visibility;
339 static double sqrt_m_log01 = sqrt( -log( 0.01 ) );
340 sun_exp2_punch_through = sqrt_m_log01 / (visibility * 15);
343 if (prev_sun_angle != sun_angle) {
344 prev_sun_angle = sun_angle;
346 double nv = (new_visibility > 5000.0) ? new_visibility : 5000.0;
347 double vis_factor = 10000.0 / (nv - 5000.0);
348 if ( vis_factor < 0.25 ) {
350 } else if ( vis_factor > 1.0) {
354 float sun_factor = 4 * (cos(sun_angle) + cos(sun_angle)/2) * vis_factor;
356 if (sun_factor > 1) sun_factor = 1.0;
357 if (sun_factor < -1) sun_factor = -1.0;
358 sun_factor = (sun_factor/2) + 0.5;
361 color[1] = sqrt(sun_factor);
362 color[0] = sqrt(color[1]);
363 color[2] = sun_factor * sun_factor;
364 color[2] *= color[2];
367 gamma_correct_rgb( color );
369 // cout << "color = " << color[0] << " " << color[1] << " "
370 // << color[2] << endl;
374 sgCopyVec4( ptr, color );
381 // reposition the sun at the specified right ascension and
382 // declination, offset by our current position (p) so that it appears
383 // fixed at a great distance from the viewer. Also add in an optional
384 // rotation (i.e. for the current time of day.)
385 bool SGSun::reposition( sgVec3 p, double angle,
386 double rightAscension, double declination,
389 sgMat4 T1, T2, GST, RA, DEC;
393 sgMakeTransMat4( T1, p );
395 sgSetVec3( axis, 0.0, 0.0, -1.0 );
396 sgMakeRotMat4( GST, angle, axis );
398 // xglRotatef( ((SGD_RADIANS_TO_DEGREES * rightAscension)- 90.0),
400 sgSetVec3( axis, 0.0, 0.0, 1.0 );
401 sgMakeRotMat4( RA, (rightAscension * SGD_RADIANS_TO_DEGREES) - 90.0, axis );
403 // xglRotatef((SGD_RADIANS_TO_DEGREES * declination), 1.0, 0.0, 0.0);
404 sgSetVec3( axis, 1.0, 0.0, 0.0 );
405 sgMakeRotMat4( DEC, declination * SGD_RADIANS_TO_DEGREES, axis );
407 // xglTranslatef(0,sun_dist);
408 sgSetVec3( v, 0.0, sun_dist, 0.0 );
409 sgMakeTransMat4( T2, v );
412 sgCopyMat4( TRANSFORM, T1 );
413 sgPreMultMat4( TRANSFORM, GST );
414 sgPreMultMat4( TRANSFORM, RA );
415 sgPreMultMat4( TRANSFORM, DEC );
416 sgPreMultMat4( TRANSFORM, T2 );
419 sgSetCoord( &skypos, TRANSFORM );
421 sun_transform->setTransform( &skypos );