Pileus Git - grits/blob - src/gis-opengl.c

   1 /*
   2  * Copyright (C) 2009 Andy Spencer <spenceal@rose-hulman.edu>
   3  *
   4  * This program is free software: you can redistribute it and/or modify
   5  * it under the terms of the GNU General Public License as published by
   6  * the Free Software Foundation, either version 3 of the License, or
   7  * (at your option) any later version.
   8  *
   9  * This program is distributed in the hope that it will be useful,
  10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  12  * GNU General Public License for more details.
  13  *
  14  * You should have received a copy of the GNU General Public License
  15  * along with this program.  If not, see <http://www.gnu.org/licenses/>.
  16  */
  17
  18 /* Tessellation, "finding intersecting triangles" */
  19 /* http://research.microsoft.com/pubs/70307/tr-2006-81.pdf */
  20 /* http://www.opengl.org/wiki/Alpha_Blending */
  21
  22 #include <config.h>
  23 #include <math.h>
  24 #include <gdk/gdkkeysyms.h>
  25 #include <gtk/gtk.h>
  26 #include <gtk/gtkgl.h>
  27 #include <GL/gl.h>
  28 #include <GL/glu.h>
  29
  30 #include "gis-opengl.h"
  31 #include "roam.h"
  32 #include "wms.h"
  33
  34 #define FOV_DIST   2000.0
  35 #define MPPX(dist) (4*dist/FOV_DIST)
  36
  37 //#define ROAM_DEBUG
  38
  39 /*************
  40  * ROAM Code *
  41  *************/
  42 void roam_queue_draw(WmsCacheNode *node, gpointer _self)
  43 {
  44         gtk_widget_queue_draw(GTK_WIDGET(_self));
  45 }
  46
  47 static void set_camera(GisOpenGL *self)
  48 {
  49         glMatrixMode(GL_MODELVIEW);
  50         glLoadIdentity();
  51         double lat, lon, elev, rx, ry, rz;
  52         gis_view_get_location(self->view, &lat, &lon, &elev);
  53         gis_view_get_rotation(self->view, &rx, &ry, &rz);
  54         glRotatef(rx, 1, 0, 0);
  55         glRotatef(rz, 0, 0, 1);
  56         glTranslatef(0, 0, -elev2rad(elev));
  57         glRotatef(lat, 1, 0, 0);
  58         glRotatef(-lon, 0, 1, 0);
  59 }
  60
  61 static void set_visuals(GisOpenGL *self)
  62 {
  63         /* Lighting */
  64         glMatrixMode(GL_MODELVIEW);
  65         glLoadIdentity();
  66 #ifdef ROAM_DEBUG
  67         float light_ambient[]  = {0.7f, 0.7f, 0.7f, 1.0f};
  68         float light_diffuse[]  = {2.0f, 2.0f, 2.0f, 1.0f};
  69 #else
  70         float light_ambient[]  = {0.2f, 0.2f, 0.2f, 1.0f};
  71         float light_diffuse[]  = {5.0f, 5.0f, 5.0f, 1.0f};
  72 #endif
  73         float light_position[] = {-13*EARTH_R, 1*EARTH_R, 3*EARTH_R, 1.0f};
  74         glLightfv(GL_LIGHT0, GL_AMBIENT,  light_ambient);
  75         glLightfv(GL_LIGHT0, GL_DIFFUSE,  light_diffuse);
  76         glLightfv(GL_LIGHT0, GL_POSITION, light_position);
  77         glEnable(GL_LIGHT0);
  78         glEnable(GL_LIGHTING);
  79
  80         float material_ambient[]  = {0.2, 0.2, 0.2, 1.0};
  81         float material_diffuse[]  = {0.8, 0.8, 0.8, 1.0};
  82         float material_specular[] = {0.0, 0.0, 0.0, 1.0};
  83         float material_emission[] = {0.0, 0.0, 0.0, 1.0};
  84         glMaterialfv(GL_FRONT_AND_BACK, GL_AMBIENT,  material_ambient);
  85         glMaterialfv(GL_FRONT_AND_BACK, GL_DIFFUSE,  material_diffuse);
  86         glMaterialfv(GL_FRONT_AND_BACK, GL_SPECULAR, material_specular);
  87         glMaterialfv(GL_FRONT_AND_BACK, GL_EMISSION, material_emission);
  88         glDisable(GL_TEXTURE_2D);
  89         glDisable(GL_COLOR_MATERIAL);
  90
  91         /* Camera */
  92         set_camera(self);
  93
  94         /* Misc */
  95         gdouble lat, lon, elev;
  96         gis_view_get_location(self->view, &lat, &lon, &elev);
  97         gdouble rg   = MAX(0, 1-(elev/20000));
  98         gdouble blue = MAX(0, 1-(elev/50000));
  99         glClearColor(MIN(0.65,rg), MIN(0.65,rg), MIN(1,blue), 1.0f);
 100
 101         glDisable(GL_ALPHA_TEST);
 102
 103         glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
 104         glEnable(GL_BLEND);
 105
 106 #ifndef ROAM_DEBUG
 107         glCullFace(GL_BACK);
 108         glEnable(GL_CULL_FACE);
 109 #endif
 110
 111         glClearDepth(1.0);
 112         glDepthFunc(GL_LEQUAL);
 113         glEnable(GL_DEPTH_TEST);
 114
 115         glEnable(GL_LINE_SMOOTH);
 116
 117         glPolygonMode(GL_FRONT_AND_BACK, GL_FILL);
 118         //glShadeModel(GL_FLAT);
 119 }
 120
 121
 122 /*************
 123  * Callbacks *
 124  *************/
 125 static void on_realize(GisOpenGL *self, gpointer _)
 126 {
 127         set_visuals(self);
 128         roam_sphere_update_errors(self->sphere);
 129 }
 130 static gboolean on_configure(GisOpenGL *self, GdkEventConfigure *event, gpointer _)
 131 {
 132         g_debug("GisOpenGL: on_confiure");
 133         gis_opengl_begin(self);
 134
 135         double width  = GTK_WIDGET(self)->allocation.width;
 136         double height = GTK_WIDGET(self)->allocation.height;
 137         glViewport(0, 0, width, height);
 138
 139         glMatrixMode(GL_PROJECTION);
 140         glLoadIdentity();
 141         double ang = atan(height/FOV_DIST);
 142         gluPerspective(rad2deg(ang)*2, width/height, 1, 20*EARTH_R);
 143
 144 #ifndef ROAM_DEBUG
 145         roam_sphere_update_errors(self->sphere);
 146 #endif
 147
 148         gis_opengl_end(self);
 149         return FALSE;
 150 }
 151
 152 static void on_expose_plugin(GisPlugin *plugin, gchar *name, GisOpenGL *self)
 153 {
 154         set_visuals(self);
 155         glMatrixMode(GL_PROJECTION); glPushMatrix();
 156         glMatrixMode(GL_MODELVIEW);  glPushMatrix();
 157         gis_plugin_expose(plugin);
 158         glMatrixMode(GL_PROJECTION); glPopMatrix();
 159         glMatrixMode(GL_MODELVIEW);  glPopMatrix();
 160 }
 161 static gboolean on_expose(GisOpenGL *self, GdkEventExpose *event, gpointer _)
 162 {
 163         g_debug("GisOpenGL: on_expose - begin");
 164         gis_opengl_begin(self);
 165
 166         glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
 167
 168 #ifndef ROAM_DEBUG
 169         set_visuals(self);
 170         glEnable(GL_TEXTURE_2D);
 171         glPolygonMode(GL_FRONT_AND_BACK, GL_LINE);
 172         roam_sphere_draw(self->sphere);
 173 #else
 174         set_visuals(self);
 175         glColor4f(0.0, 0.0, 9.0, 0.6);
 176         glDisable(GL_TEXTURE_2D);
 177         glPolygonMode(GL_FRONT_AND_BACK, GL_LINE);
 178         roam_sphere_draw(self->sphere);
 179
 180         //roam_sphere_draw_normals(self->sphere);
 181 #endif
 182
 183         gis_plugins_foreach(self->plugins, G_CALLBACK(on_expose_plugin), self);
 184
 185         set_visuals(self);
 186         gis_opengl_end(self);
 187         gis_opengl_flush(self);
 188         g_debug("GisOpenGL: on_expose - end\n");
 189         return FALSE;
 190 }
 191
 192 static gboolean on_button_press(GisOpenGL *self, GdkEventButton *event, gpointer _)
 193 {
 194         g_debug("GisOpenGL: on_button_press - Grabbing focus");
 195         gtk_widget_grab_focus(GTK_WIDGET(self));
 196         return TRUE;
 197 }
 198
 199 static gboolean on_key_press(GisOpenGL *self, GdkEventKey *event, gpointer _)
 200 {
 201         g_debug("GisOpenGL: on_key_press - key=%x, state=%x, plus=%x",
 202                         event->keyval, event->state, GDK_plus);
 203
 204         double lat, lon, elev, pan;
 205         gis_view_get_location(self->view, &lat, &lon, &elev);
 206         pan = MIN(elev/(EARTH_R/2), 30);
 207         guint kv = event->keyval;
 208         if      (kv == GDK_Left  || kv == GDK_h) gis_view_pan(self->view,  0,  -pan, 0);
 209         else if (kv == GDK_Down  || kv == GDK_j) gis_view_pan(self->view, -pan, 0,   0);
 210         else if (kv == GDK_Up    || kv == GDK_k) gis_view_pan(self->view,  pan, 0,   0);
 211         else if (kv == GDK_Right || kv == GDK_l) gis_view_pan(self->view,  0,   pan, 0);
 212         else if (kv == GDK_minus || kv == GDK_o) gis_view_zoom(self->view, 10./9);
 213         else if (kv == GDK_plus  || kv == GDK_i) gis_view_zoom(self->view, 9./10);
 214         else if (kv == GDK_H) gis_view_rotate(self->view,  0,  0, -10);
 215         else if (kv == GDK_J) gis_view_rotate(self->view,  10, 0,  0);
 216         else if (kv == GDK_K) gis_view_rotate(self->view, -10, 0,  0);
 217         else if (kv == GDK_L) gis_view_rotate(self->view,  0,  0,  10);
 218
 219         /* Testing */
 220 #ifdef ROAM_DEBUG
 221         else if (kv == GDK_n) roam_sphere_split_one(self->sphere);
 222         else if (kv == GDK_p) roam_sphere_merge_one(self->sphere);
 223         else if (kv == GDK_r) roam_sphere_split_merge(self->sphere);
 224         else if (kv == GDK_u) roam_sphere_update_errors(self->sphere);
 225         gtk_widget_queue_draw(GTK_WIDGET(self));
 226 #endif
 227
 228         return TRUE;
 229 }
 230
 231 static void on_view_changed(GisView *view,
 232                 gdouble _1, gdouble _2, gdouble _3, GisOpenGL *self)
 233 {
 234         gis_opengl_begin(self);
 235         set_visuals(self);
 236 #ifndef ROAM_DEBUG
 237         roam_sphere_update_errors(self->sphere);
 238 #endif
 239         gis_opengl_redraw(self);
 240         gis_opengl_end(self);
 241 }
 242
 243 static gboolean on_idle(GisOpenGL *self)
 244 {
 245         gis_opengl_begin(self);
 246         if (roam_sphere_split_merge(self->sphere))
 247                 gis_opengl_redraw(self);
 248         gis_opengl_end(self);
 249         return TRUE;
 250 }
 251
 252
 253 /***********
 254  * Methods *
 255  ***********/
 256 GisOpenGL *gis_opengl_new(GisWorld *world, GisView *view, GisPlugins *plugins)
 257 {
 258         g_debug("GisOpenGL: new");
 259         GisOpenGL *self = g_object_new(GIS_TYPE_OPENGL, NULL);
 260         self->world   = world;
 261         self->view    = view;
 262         self->plugins = plugins;
 263         g_object_ref(world);
 264         g_object_ref(view);
 265
 266         g_signal_connect(self->view, "location-changed", G_CALLBACK(on_view_changed), self);
 267         g_signal_connect(self->view, "rotation-changed", G_CALLBACK(on_view_changed), self);
 268
 269         /* TODO: update point eights sometime later so we have heigh-res heights for them */
 270         self->sphere = roam_sphere_new(self);
 271
 272         return g_object_ref(self);
 273 }
 274
 275 void gis_opengl_center_position(GisOpenGL *self, gdouble lat, gdouble lon, gdouble elev)
 276 {
 277         set_camera(self);
 278         glRotatef(lon, 0, 1, 0);
 279         glRotatef(-lat, 1, 0, 0);
 280         glTranslatef(0, 0, elev2rad(elev));
 281 }
 282
 283 void gis_opengl_redraw(GisOpenGL *self)
 284 {
 285         g_debug("GisOpenGL: gl_redraw");
 286         gtk_widget_queue_draw(GTK_WIDGET(self));
 287 }
 288 void gis_opengl_begin(GisOpenGL *self)
 289 {
 290         g_assert(GIS_IS_OPENGL(self));
 291
 292         GdkGLContext   *glcontext  = gtk_widget_get_gl_context(GTK_WIDGET(self));
 293         GdkGLDrawable  *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
 294
 295         if (!gdk_gl_drawable_gl_begin(gldrawable, glcontext))
 296                 g_assert_not_reached();
 297 }
 298 void gis_opengl_end(GisOpenGL *self)
 299 {
 300         g_assert(GIS_IS_OPENGL(self));
 301         GdkGLDrawable *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
 302         gdk_gl_drawable_gl_end(gldrawable);
 303 }
 304 void gis_opengl_flush(GisOpenGL *self)
 305 {
 306         g_assert(GIS_IS_OPENGL(self));
 307         GdkGLDrawable *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
 308         if (gdk_gl_drawable_is_double_buffered(gldrawable))
 309                 gdk_gl_drawable_swap_buffers(gldrawable);
 310         else
 311                 glFlush();
 312         gdk_gl_drawable_gl_end(gldrawable);
 313 }
 314
 315
 316 /****************
 317  * GObject code *
 318  ****************/
 319 G_DEFINE_TYPE(GisOpenGL, gis_opengl, GTK_TYPE_DRAWING_AREA);
 320 static void gis_opengl_init(GisOpenGL *self)
 321 {
 322         g_debug("GisOpenGL: init");
 323         self->bmng = wms_info_new_for_bmng(NULL, NULL);
 324         self->srtm = wms_info_new_for_srtm(NULL, NULL);
 325
 326         /* OpenGL setup */
 327         GdkGLConfig *glconfig = gdk_gl_config_new_by_mode(
 328                         GDK_GL_MODE_RGBA   | GDK_GL_MODE_DEPTH |
 329                         GDK_GL_MODE_DOUBLE | GDK_GL_MODE_ALPHA);
 330         if (!glconfig)
 331                 g_error("Failed to create glconfig");
 332         if (!gtk_widget_set_gl_capability(GTK_WIDGET(self),
 333                                 glconfig, NULL, TRUE, GDK_GL_RGBA_TYPE))
 334                 g_error("GL lacks required capabilities");
 335         g_object_unref(glconfig);
 336
 337         gtk_widget_set_size_request(GTK_WIDGET(self), 600, 550);
 338         gtk_widget_set_events(GTK_WIDGET(self),
 339                         GDK_BUTTON_PRESS_MASK |
 340                         GDK_ENTER_NOTIFY_MASK |
 341                         GDK_KEY_PRESS_MASK);
 342         g_object_set(self, "can-focus", TRUE, NULL);
 343
 344 #ifndef ROAM_DEBUG
 345         self->sm_source = g_timeout_add(10, (GSourceFunc)on_idle, self);
 346 #endif
 347
 348         g_signal_connect(self, "realize",            G_CALLBACK(on_realize),      NULL);
 349         g_signal_connect(self, "configure-event",    G_CALLBACK(on_configure),    NULL);
 350         g_signal_connect(self, "expose-event",       G_CALLBACK(on_expose),       NULL);
 351
 352         g_signal_connect(self, "button-press-event", G_CALLBACK(on_button_press), NULL);
 353         g_signal_connect(self, "enter-notify-event", G_CALLBACK(on_button_press), NULL);
 354         g_signal_connect(self, "key-press-event",    G_CALLBACK(on_key_press),    NULL);
 355 }
 356 static GObject *gis_opengl_constructor(GType gtype, guint n_properties,
 357                 GObjectConstructParam *properties)
 358 {
 359         g_debug("GisOpengl: constructor");
 360         GObjectClass *parent_class = G_OBJECT_CLASS(gis_opengl_parent_class);
 361         return parent_class->constructor(gtype, n_properties, properties);
 362 }
 363 static void gis_opengl_dispose(GObject *_self)
 364 {
 365         g_debug("GisOpenGL: dispose");
 366         GisOpenGL *self = GIS_OPENGL(_self);
 367         if (self->sm_source) {
 368                 g_source_remove(self->sm_source);
 369                 self->sm_source = 0;
 370         }
 371         if (self->sphere) {
 372                 roam_sphere_free(self->sphere);
 373                 self->sphere = NULL;
 374         }
 375         if (self->world) {
 376                 g_object_unref(self->world);
 377                 self->world = NULL;
 378         }
 379         if (self->view) {
 380                 g_object_unref(self->view);
 381                 self->view = NULL;
 382         }
 383         G_OBJECT_CLASS(gis_opengl_parent_class)->dispose(_self);
 384 }
 385 static void gis_opengl_finalize(GObject *_self)
 386 {
 387         g_debug("GisOpenGL: finalize");
 388         GisOpenGL *self = GIS_OPENGL(_self);
 389         wms_info_free(self->bmng);
 390         wms_info_free(self->srtm);
 391         G_OBJECT_CLASS(gis_opengl_parent_class)->finalize(_self);
 392 }
 393 static void gis_opengl_class_init(GisOpenGLClass *klass)
 394 {
 395         g_debug("GisOpenGL: class_init");
 396         GObjectClass *gobject_class = G_OBJECT_CLASS(klass);
 397         gobject_class->constructor  = gis_opengl_constructor;
 398         gobject_class->dispose      = gis_opengl_dispose;
 399         gobject_class->finalize     = gis_opengl_finalize;
 400 }