1 // dome.cxx -- model sky with an upside down "bowl"
3 // Written by Curtis Olson, started December 1997.
4 // SSG-ified by Curtis Olson, February 2000.
6 // Copyright (C) 1997-2000 Curtis L. Olson - curt@flightgear.org
8 // This library is free software; you can redistribute it and/or
9 // modify it under the terms of the GNU Library General Public
10 // License as published by the Free Software Foundation; either
11 // version 2 of the License, or (at your option) any later version.
13 // This library is distributed in the hope that it will be useful,
14 // but WITHOUT ANY WARRANTY; without even the implied warranty of
15 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 // Library General Public License for more details.
18 // You should have received a copy of the GNU Library General Public
19 // License along with this library; if not, write to the
20 // Free Software Foundation, Inc., 59 Temple Place - Suite 330,
21 // Boston, MA 02111-1307, USA.
27 # include <simgear_config.h>
40 #include <simgear/debug/logstream.hxx>
46 # pragma global_optimizer off
50 // proportions of max dimensions fed to the build() routine
51 static const float center_elev = 1.0;
53 static const float upper_radius = 0.6;
54 static const float upper_elev = 0.15;
56 static const float middle_radius = 0.9;
57 static const float middle_elev = 0.08;
59 static const float lower_radius = 1.0;
60 static const float lower_elev = 0.0;
62 static const float bottom_radius = 0.8;
63 static const float bottom_elev = -0.1;
66 // Set up dome rendering callbacks
67 static int sgSkyDomePreDraw( ssgEntity *e ) {
68 /* cout << endl << "Dome Pre Draw" << endl << "----------------"
71 ssgLeaf *f = (ssgLeaf *)e;
72 if ( f -> hasState () ) f->getState()->apply() ;
74 glPushAttrib( GL_DEPTH_BUFFER_BIT | GL_FOG_BIT );
75 // cout << "push error = " << glGetError() << endl;
77 glDisable( GL_DEPTH_TEST );
83 static int sgSkyDomePostDraw( ssgEntity *e ) {
84 /* cout << endl << "Dome Post Draw" << endl << "----------------"
88 // cout << "pop error = " << glGetError() << endl;
95 SGSkyDome::SGSkyDome( void ) {
100 SGSkyDome::~SGSkyDome( void ) {
104 // initialize the sky object and connect it into our scene graph
105 ssgBranch * SGSkyDome::build( double hscale, double vscale ) {
112 dome_state = new ssgSimpleState();
113 dome_state->setShadeModel( GL_SMOOTH );
114 dome_state->disable( GL_LIGHTING );
115 dome_state->disable( GL_CULL_FACE );
116 dome_state->disable( GL_TEXTURE_2D );
117 dome_state->enable( GL_COLOR_MATERIAL );
118 dome_state->setColourMaterial( GL_AMBIENT_AND_DIFFUSE );
119 dome_state->setMaterial( GL_EMISSION, 0, 0, 0, 1 );
120 dome_state->setMaterial( GL_SPECULAR, 0, 0, 0, 1 );
121 dome_state->disable( GL_BLEND );
122 dome_state->disable( GL_ALPHA_TEST );
125 center_disk_vl = new ssgVertexArray( 14 );
126 center_disk_cl = new ssgColourArray( 14 );
128 upper_ring_vl = new ssgVertexArray( 26 );
129 upper_ring_cl = new ssgColourArray( 26 );
131 middle_ring_vl = new ssgVertexArray( 26 );
132 middle_ring_cl = new ssgColourArray( 26 );
134 lower_ring_vl = new ssgVertexArray( 26 );
135 lower_ring_cl = new ssgColourArray( 26 );
137 // initially seed to all blue
138 sgSetVec4( color, 0.0, 0.0, 1.0, 1.0 );
140 // generate the raw vertex data
141 sgVec3 center_vertex;
142 sgVec3 upper_vertex[12];
143 sgVec3 middle_vertex[12];
144 sgVec3 lower_vertex[12];
145 sgVec3 bottom_vertex[12];
147 sgSetVec3( center_vertex, 0.0, 0.0, center_elev * vscale );
149 for ( i = 0; i < 12; i++ ) {
150 theta = (i * 30.0) * SGD_DEGREES_TO_RADIANS;
152 sgSetVec3( upper_vertex[i],
153 cos(theta) * upper_radius * hscale,
154 sin(theta) * upper_radius * hscale,
155 upper_elev * vscale );
157 sgSetVec3( middle_vertex[i],
158 cos(theta) * middle_radius * hscale,
159 sin(theta) * middle_radius * hscale,
160 middle_elev * vscale );
162 sgSetVec3( lower_vertex[i],
163 cos(theta) * lower_radius * hscale,
164 sin(theta) * lower_radius * hscale,
165 lower_elev * vscale );
167 sgSetVec3( bottom_vertex[i],
168 cos(theta) * bottom_radius * hscale,
169 sin(theta) * bottom_radius * hscale,
170 bottom_elev * vscale );
173 // generate the center disk vertex/color arrays
174 center_disk_vl->add( center_vertex );
175 center_disk_cl->add( color );
176 for ( i = 11; i >= 0; i-- ) {
177 center_disk_vl->add( upper_vertex[i] );
178 center_disk_cl->add( color );
180 center_disk_vl->add( upper_vertex[11] );
181 center_disk_cl->add( color );
183 // generate the upper ring
184 for ( i = 0; i < 12; i++ ) {
185 upper_ring_vl->add( middle_vertex[i] );
186 upper_ring_cl->add( color );
188 upper_ring_vl->add( upper_vertex[i] );
189 upper_ring_cl->add( color );
191 upper_ring_vl->add( middle_vertex[0] );
192 upper_ring_cl->add( color );
194 upper_ring_vl->add( upper_vertex[0] );
195 upper_ring_cl->add( color );
197 // generate middle ring
198 for ( i = 0; i < 12; i++ ) {
199 middle_ring_vl->add( lower_vertex[i] );
200 middle_ring_cl->add( color );
202 middle_ring_vl->add( middle_vertex[i] );
203 middle_ring_cl->add( color );
205 middle_ring_vl->add( lower_vertex[0] );
206 middle_ring_cl->add( color );
208 middle_ring_vl->add( middle_vertex[0] );
209 middle_ring_cl->add( color );
211 // generate lower ring
212 for ( i = 0; i < 12; i++ ) {
213 lower_ring_vl->add( bottom_vertex[i] );
214 lower_ring_cl->add( color );
216 lower_ring_vl->add( lower_vertex[i] );
217 lower_ring_cl->add( color );
219 lower_ring_vl->add( bottom_vertex[0] );
220 lower_ring_cl->add( color );
222 lower_ring_vl->add( lower_vertex[0] );
223 lower_ring_cl->add( color );
225 // force a repaint of the sky colors with ugly defaults
227 sgSetVec4( fog_color, 1.0, 1.0, 1.0, 1.0 );
228 repaint( color, fog_color, 0.0, 5000.0 );
230 // build the ssg scene graph sub tree for the sky and connected
231 // into the provide scene graph branch
232 ssgVtxTable *center_disk, *upper_ring, *middle_ring, *lower_ring;
234 center_disk = new ssgVtxTable( GL_TRIANGLE_FAN,
235 center_disk_vl, NULL, NULL, center_disk_cl );
237 upper_ring = new ssgVtxTable( GL_TRIANGLE_STRIP,
238 upper_ring_vl, NULL, NULL, upper_ring_cl );
240 middle_ring = new ssgVtxTable( GL_TRIANGLE_STRIP,
241 middle_ring_vl, NULL, NULL, middle_ring_cl );
243 lower_ring = new ssgVtxTable( GL_TRIANGLE_STRIP,
244 lower_ring_vl, NULL, NULL, lower_ring_cl );
246 center_disk->setState( dome_state );
247 upper_ring->setState( dome_state );
248 middle_ring->setState( dome_state );
249 lower_ring->setState( dome_state );
251 dome_transform = new ssgTransform;
252 dome_transform->addKid( center_disk );
253 dome_transform->addKid( upper_ring );
254 dome_transform->addKid( middle_ring );
255 dome_transform->addKid( lower_ring );
257 // not entirely satisfying. We are depending here that the first
258 // thing we add to a parent is the first drawn
259 center_disk->setCallback( SSG_CALLBACK_PREDRAW, sgSkyDomePreDraw );
260 center_disk->setCallback( SSG_CALLBACK_POSTDRAW, sgSkyDomePostDraw );
262 upper_ring->setCallback( SSG_CALLBACK_PREDRAW, sgSkyDomePreDraw );
263 upper_ring->setCallback( SSG_CALLBACK_POSTDRAW, sgSkyDomePostDraw );
265 middle_ring->setCallback( SSG_CALLBACK_PREDRAW, sgSkyDomePreDraw );
266 middle_ring->setCallback( SSG_CALLBACK_POSTDRAW, sgSkyDomePostDraw );
268 lower_ring->setCallback( SSG_CALLBACK_PREDRAW, sgSkyDomePreDraw );
269 lower_ring->setCallback( SSG_CALLBACK_POSTDRAW, sgSkyDomePostDraw );
271 return dome_transform;
275 // repaint the sky colors based on current value of sun_angle, sky,
276 // and fog colors. This updates the color arrays for ssgVtxTable.
277 // sun angle in degrees relative to verticle
278 // 0 degrees = high noon
279 // 90 degrees = sun rise/set
280 // 180 degrees = darkest midnight
281 bool SGSkyDome::repaint( sgVec4 sky_color, sgVec4 fog_color, double sun_angle,
284 double diff, prev_sun_angle = 999.0;
285 sgVec3 outer_param, outer_amt, outer_diff;
286 sgVec3 middle_param, middle_amt, middle_diff;
289 if (prev_sun_angle == sun_angle)
292 prev_sun_angle = sun_angle;
294 // Check for sunrise/sunset condition
295 if (sun_angle > 80.0) // && (sun_angle < 100.0) )
298 sgSetVec3( outer_param,
299 (10.0 - fabs(90.0 - sun_angle)) / 20.0,
300 (10.0 - fabs(90.0 - sun_angle)) / 40.0,
301 -(10.0 - fabs(90.0 - sun_angle)) / 30.0 );
303 sgSetVec3( middle_param,
304 (10.0 - fabs(90.0 - sun_angle)) / 40.0,
305 (10.0 - fabs(90.0 - sun_angle)) / 80.0,
308 sgScaleVec3( outer_diff, outer_param, 1.0 / 6.0 );
310 sgScaleVec3( middle_diff, middle_param, 1.0 / 6.0 );
312 sgSetVec3( outer_param, 0.0, 0.0, 0.0 );
313 sgSetVec3( middle_param, 0.0, 0.0, 0.0 );
315 sgSetVec3( outer_diff, 0.0, 0.0, 0.0 );
316 sgSetVec3( middle_diff, 0.0, 0.0, 0.0 );
318 // printf(" outer_red_param = %.2f outer_red_diff = %.2f\n",
319 // outer_red_param, outer_red_diff);
321 // calculate transition colors between sky and fog
322 sgCopyVec3( outer_amt, outer_param );
323 sgCopyVec3( middle_amt, middle_param );
326 // First, recalulate the basic colors
330 sgVec4 upper_color[12];
331 sgVec4 middle_color[12];
332 sgVec4 lower_color[12];
333 sgVec4 bottom_color[12];
335 double vis_factor, cvf = vis;
340 vis_factor = (vis - 1000.0) / 2000.0;
341 if ( vis_factor < 0.0 ) {
343 } else if ( vis_factor > 1.0) {
347 for ( j = 0; j < 3; j++ ) {
348 diff = sky_color[j] - fog_color[j];
349 center_color[j] = sky_color[j]; // - diff * ( 1.0 - vis_factor );
351 center_color[3] = 1.0;
353 for ( i = 0; i < 6; i++ ) {
354 for ( j = 0; j < 3; j++ ) {
355 double saif = sun_angle/SG_PI;
356 diff = (sky_color[j] - fog_color[j]) * (0.8 + j * 0.2) * (0.8 + saif - ((6-i)/10));
358 // printf("sky = %.2f fog = %.2f diff = %.2f\n",
359 // l->sky_color[j], l->fog_color[j], diff);
361 upper_color[i][j] = sky_color[j] - diff *
362 ( 1.0 - vis_factor * (0.7 + 0.3 * cvf/45000) );
363 middle_color[i][j] = sky_color[j] - diff *
364 ( 1.0 - vis_factor * (0.1 + 0.85 * cvf/45000) )
366 lower_color[i][j] = fog_color[j] + outer_amt[j];
368 if ( upper_color[i][j] > 1.0 ) { upper_color[i][j] = 1.0; }
369 if ( upper_color[i][j] < 0.0 ) { upper_color[i][j] = 0.0; }
370 if ( middle_color[i][j] > 1.0 ) { middle_color[i][j] = 1.0; }
371 if ( middle_color[i][j] < 0.0 ) { middle_color[i][j] = 0.0; }
372 if ( lower_color[i][j] > 1.0 ) { lower_color[i][j] = 1.0; }
373 if ( lower_color[i][j] < 0.0 ) { lower_color[i][j] = 0.0; }
375 upper_color[i][3] = middle_color[i][3] = lower_color[i][3] = 1.0;
377 for ( j = 0; j < 3; j++ ) {
378 outer_amt[j] -= outer_diff[j];
379 middle_amt[j] -= middle_diff[j];
383 printf("upper_color[%d] = %.2f %.2f %.2f %.2f\n", i, upper_color[i][0],
384 upper_color[i][1], upper_color[i][2], upper_color[i][3]);
385 printf("middle_color[%d] = %.2f %.2f %.2f %.2f\n", i,
386 middle_color[i][0], middle_color[i][1], middle_color[i][2],
388 printf("lower_color[%d] = %.2f %.2f %.2f %.2f\n", i,
389 lower_color[i][0], lower_color[i][1], lower_color[i][2],
394 sgSetVec3( outer_amt, 0.0, 0.0, 0.0 );
395 sgSetVec3( middle_amt, 0.0, 0.0, 0.0 );
397 for ( i = 6; i < 12; i++ ) {
398 for ( j = 0; j < 3; j++ ) {
399 double saif = sun_angle/SG_PI;
400 diff = (sky_color[j] - fog_color[j]) * (0.8 + j * 0.2) * (0.8 + saif - ((-i+12)/10));
402 // printf("sky = %.2f fog = %.2f diff = %.2f\n",
403 // sky_color[j], fog_color[j], diff);
405 upper_color[i][j] = sky_color[j] - diff *
406 ( 1.0 - vis_factor * (0.7 + 0.3 * cvf/45000) );
407 middle_color[i][j] = sky_color[j] - diff *
408 ( 1.0 - vis_factor * (0.1 + 0.85 * cvf/45000) )
410 lower_color[i][j] = fog_color[j] + outer_amt[j];
412 if ( upper_color[i][j] > 1.0 ) { upper_color[i][j] = 1.0; }
413 if ( upper_color[i][j] < 0.0 ) { upper_color[i][j] = 0.0; }
414 if ( middle_color[i][j] > 1.0 ) { middle_color[i][j] = 1.0; }
415 if ( middle_color[i][j] < 0.0 ) { middle_color[i][j] = 0.0; }
416 if ( lower_color[i][j] > 1.0 ) { lower_color[i][j] = 1.0; }
417 if ( lower_color[i][j] < 0.0 ) { lower_color[i][j] = 0.0; }
419 upper_color[i][3] = middle_color[i][3] = lower_color[i][3] = 1.0;
421 for ( j = 0; j < 3; j++ ) {
422 outer_amt[j] += outer_diff[j];
423 middle_amt[j] += middle_diff[j];
427 printf("upper_color[%d] = %.2f %.2f %.2f %.2f\n", i, upper_color[i][0],
428 upper_color[i][1], upper_color[i][2], upper_color[i][3]);
429 printf("middle_color[%d] = %.2f %.2f %.2f %.2f\n", i,
430 middle_color[i][0], middle_color[i][1], middle_color[i][2],
432 printf("lower_color[%d] = %.2f %.2f %.2f %.2f\n", i,
433 lower_color[i][0], lower_color[i][1], lower_color[i][2],
438 for ( i = 0; i < 12; i++ ) {
439 sgCopyVec4( bottom_color[i], fog_color );
443 // Second, assign the basic colors to the object color arrays
449 // update the center disk color arrays
451 slot = center_disk_cl->get( counter++ );
453 // sgSetVec4( red, 1.0, 0.0, 0.0, 1.0 );
454 sgCopyVec4( slot, center_color );
455 for ( i = 11; i >= 0; i-- ) {
456 slot = center_disk_cl->get( counter++ );
457 sgCopyVec4( slot, upper_color[i] );
459 slot = center_disk_cl->get( counter++ );
460 sgCopyVec4( slot, upper_color[11] );
462 // generate the upper ring
464 for ( i = 0; i < 12; i++ ) {
465 slot = upper_ring_cl->get( counter++ );
466 sgCopyVec4( slot, middle_color[i] );
468 slot = upper_ring_cl->get( counter++ );
469 sgCopyVec4( slot, upper_color[i] );
471 slot = upper_ring_cl->get( counter++ );
472 sgCopyVec4( slot, middle_color[0] );
474 slot = upper_ring_cl->get( counter++ );
475 sgCopyVec4( slot, upper_color[0] );
477 // generate middle ring
479 for ( i = 0; i < 12; i++ ) {
480 slot = middle_ring_cl->get( counter++ );
481 sgCopyVec4( slot, lower_color[i] );
483 slot = middle_ring_cl->get( counter++ );
484 sgCopyVec4( slot, middle_color[i] );
486 slot = middle_ring_cl->get( counter++ );
487 sgCopyVec4( slot, lower_color[0] );
489 slot = middle_ring_cl->get( counter++ );
490 sgCopyVec4( slot, middle_color[0] );
492 // generate lower ring
494 for ( i = 0; i < 12; i++ ) {
495 slot = lower_ring_cl->get( counter++ );
496 sgCopyVec4( slot, bottom_color[i] );
498 slot = lower_ring_cl->get( counter++ );
499 sgCopyVec4( slot, lower_color[i] );
501 slot = lower_ring_cl->get( counter++ );
502 sgCopyVec4( slot, bottom_color[0] );
504 slot = lower_ring_cl->get( counter++ );
505 sgCopyVec4( slot, lower_color[0] );
511 // reposition the sky at the specified origin and orientation
512 // lon specifies a rotation about the Z axis
513 // lat specifies a rotation about the new Y axis
514 // spin specifies a rotation about the new Z axis (and orients the
515 // sunrise/set effects
516 bool SGSkyDome::reposition( sgVec3 p, double lon, double lat, double spin ) {
517 sgMat4 T, LON, LAT, SPIN;
520 // Translate to view position
521 // Point3D zero_elev = current_view.get_cur_zero_elev();
522 // xglTranslatef( zero_elev.x(), zero_elev.y(), zero_elev.z() );
523 sgMakeTransMat4( T, p );
525 // printf(" Translated to %.2f %.2f %.2f\n",
526 // zero_elev.x, zero_elev.y, zero_elev.z );
528 // Rotate to proper orientation
529 // printf(" lon = %.2f lat = %.2f\n",
530 // lon * SGD_RADIANS_TO_DEGREES,
531 // lat * SGD_RADIANS_TO_DEGREES);
532 // xglRotatef( lon * SGD_RADIANS_TO_DEGREES, 0.0, 0.0, 1.0 );
533 sgSetVec3( axis, 0.0, 0.0, 1.0 );
534 sgMakeRotMat4( LON, lon * SGD_RADIANS_TO_DEGREES, axis );
536 // xglRotatef( 90.0 - f->get_Latitude() * SGD_RADIANS_TO_DEGREES,
538 sgSetVec3( axis, 0.0, 1.0, 0.0 );
539 sgMakeRotMat4( LAT, 90.0 - lat * SGD_RADIANS_TO_DEGREES, axis );
541 // xglRotatef( l->sun_rotation * SGD_RADIANS_TO_DEGREES, 0.0, 0.0, 1.0 );
542 sgSetVec3( axis, 0.0, 0.0, 1.0 );
543 sgMakeRotMat4( SPIN, spin * SGD_RADIANS_TO_DEGREES, axis );
547 sgCopyMat4( TRANSFORM, T );
548 sgPreMultMat4( TRANSFORM, LON );
549 sgPreMultMat4( TRANSFORM, LAT );
550 sgPreMultMat4( TRANSFORM, SPIN );
553 sgSetCoord( &skypos, TRANSFORM );
555 dome_transform->setTransform( &skypos );