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
  33 #define FOV_DIST   2000.0
  34 #define MPPX(dist) (4*dist/FOV_DIST)
  35
  36 // #define ROAM_DEBUG
  37
  38 /*************
  39  * ROAM Code *
  40  *************/
  41 static void set_visuals(GisOpenGL *self)
  42 {
  43         glMatrixMode(GL_MODELVIEW);
  44         glLoadIdentity();
  45
  46         /* Camera 1 */
  47         double lat, lon, elev, rx, ry, rz;
  48         gis_view_get_location(self->view, &lat, &lon, &elev);
  49         gis_view_get_rotation(self->view, &rx, &ry, &rz);
  50         glRotatef(rx, 1, 0, 0);
  51         glRotatef(rz, 0, 0, 1);
  52
  53         /* Lighting */
  54 #ifdef ROAM_DEBUG
  55         float light_ambient[]  = {0.7f, 0.7f, 0.7f, 1.0f};
  56         float light_diffuse[]  = {2.0f, 2.0f, 2.0f, 1.0f};
  57 #else
  58         float light_ambient[]  = {0.2f, 0.2f, 0.2f, 1.0f};
  59         float light_diffuse[]  = {5.0f, 5.0f, 5.0f, 1.0f};
  60 #endif
  61         float light_position[] = {-13*EARTH_R, 1*EARTH_R, 3*EARTH_R, 1.0f};
  62         glLightfv(GL_LIGHT0, GL_AMBIENT,  light_ambient);
  63         glLightfv(GL_LIGHT0, GL_DIFFUSE,  light_diffuse);
  64         glLightfv(GL_LIGHT0, GL_POSITION, light_position);
  65         glEnable(GL_LIGHT0);
  66         glEnable(GL_LIGHTING);
  67
  68         float material_ambient[]  = {0.2, 0.2, 0.2, 1.0};
  69         float material_diffuse[]  = {0.8, 0.8, 0.8, 1.0};
  70         float material_specular[] = {0.0, 0.0, 0.0, 1.0};
  71         float material_emission[] = {0.0, 0.0, 0.0, 1.0};
  72         glMaterialfv(GL_FRONT_AND_BACK, GL_AMBIENT,  material_ambient);
  73         glMaterialfv(GL_FRONT_AND_BACK, GL_DIFFUSE,  material_diffuse);
  74         glMaterialfv(GL_FRONT_AND_BACK, GL_SPECULAR, material_specular);
  75         glMaterialfv(GL_FRONT_AND_BACK, GL_EMISSION, material_emission);
  76         glDisable(GL_TEXTURE_2D);
  77         glDisable(GL_COLOR_MATERIAL);
  78
  79         /* Camera 2 */
  80         glTranslatef(0, 0, -elev2rad(elev));
  81         glRotatef(lat, 1, 0, 0);
  82         glRotatef(-lon, 0, 1, 0);
  83
  84         /* Misc */
  85         gdouble rg   = MAX(0, 1-(elev/20000));
  86         gdouble blue = MAX(0, 1-(elev/50000));
  87         glClearColor(MIN(0.65,rg), MIN(0.65,rg), MIN(1,blue), 1.0f);
  88         glColor4f(1, 1, 1, 1);
  89
  90         glDisable(GL_ALPHA_TEST);
  91
  92         glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
  93         glEnable(GL_BLEND);
  94
  95 #ifndef ROAM_DEBUG
  96         glCullFace(GL_BACK);
  97         glEnable(GL_CULL_FACE);
  98 #endif
  99
 100         glClearDepth(1.0);
 101         glDepthFunc(GL_LEQUAL);
 102         glEnable(GL_DEPTH_TEST);
 103
 104         glEnable(GL_LINE_SMOOTH);
 105
 106         glPolygonMode(GL_FRONT_AND_BACK, GL_FILL);
 107         //glShadeModel(GL_FLAT);
 108
 109         roam_sphere_update_view(self->sphere);
 110 }
 111
 112
 113 /*************
 114  * Callbacks *
 115  *************/
 116 static void on_realize(GisOpenGL *self, gpointer _)
 117 {
 118         g_debug("GisOpenGL: on_realize");
 119         set_visuals(self);
 120         roam_sphere_update_errors(self->sphere);
 121 }
 122 static gboolean on_configure(GisOpenGL *self, GdkEventConfigure *event, gpointer _)
 123 {
 124         g_debug("GisOpenGL: on_configure");
 125         gis_opengl_begin(self);
 126
 127         double width  = GTK_WIDGET(self)->allocation.width;
 128         double height = GTK_WIDGET(self)->allocation.height;
 129         glViewport(0, 0, width, height);
 130
 131         glMatrixMode(GL_PROJECTION);
 132         glLoadIdentity();
 133         double ang = atan(height/FOV_DIST);
 134         gluPerspective(rad2deg(ang)*2, width/height, 1, 20*EARTH_R);
 135
 136 #ifndef ROAM_DEBUG
 137         roam_sphere_update_errors(self->sphere);
 138 #endif
 139
 140         gis_opengl_end(self);
 141         return FALSE;
 142 }
 143
 144 static void on_expose_plugin(GisPlugin *plugin, gchar *name, GisOpenGL *self)
 145 {
 146         set_visuals(self);
 147         glMatrixMode(GL_PROJECTION); glPushMatrix();
 148         glMatrixMode(GL_MODELVIEW);  glPushMatrix();
 149         gis_plugin_expose(plugin);
 150         glMatrixMode(GL_PROJECTION); glPopMatrix();
 151         glMatrixMode(GL_MODELVIEW);  glPopMatrix();
 152 }
 153 static gboolean on_expose(GisOpenGL *self, GdkEventExpose *event, gpointer _)
 154 {
 155         g_debug("GisOpenGL: on_expose - begin");
 156         gis_opengl_begin(self);
 157
 158         glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
 159
 160 #ifndef ROAM_DEBUG
 161         gis_plugins_foreach(self->plugins, G_CALLBACK(on_expose_plugin), self);
 162
 163         if (self->wireframe) {
 164                 set_visuals(self);
 165                 glPolygonMode(GL_FRONT_AND_BACK, GL_LINE);
 166                 roam_sphere_draw(self->sphere);
 167         }
 168 #else
 169         set_visuals(self);
 170         glColor4f(0.0, 0.0, 9.0, 0.6);
 171         glDisable(GL_TEXTURE_2D);
 172         glPolygonMode(GL_FRONT_AND_BACK, GL_LINE);
 173         roam_sphere_draw(self->sphere);
 174
 175         //roam_sphere_draw_normals(self->sphere);
 176 #endif
 177
 178         set_visuals(self);
 179         gis_opengl_end(self);
 180         gis_opengl_flush(self);
 181         g_debug("GisOpenGL: on_expose - end\n");
 182         return FALSE;
 183 }
 184
 185 static gboolean on_button_press(GisOpenGL *self, GdkEventButton *event, gpointer _)
 186 {
 187         g_debug("GisOpenGL: on_button_press - Grabbing focus");
 188         gtk_widget_grab_focus(GTK_WIDGET(self));
 189         return TRUE;
 190 }
 191
 192 static gboolean on_key_press(GisOpenGL *self, GdkEventKey *event, gpointer _)
 193 {
 194         g_debug("GisOpenGL: on_key_press - key=%x, state=%x, plus=%x",
 195                         event->keyval, event->state, GDK_plus);
 196
 197         double lat, lon, elev, pan;
 198         gis_view_get_location(self->view, &lat, &lon, &elev);
 199         pan = MIN(elev/(EARTH_R/2), 30);
 200         guint kv = event->keyval;
 201         gdk_threads_leave();
 202         if      (kv == GDK_Left  || kv == GDK_h) gis_view_pan(self->view,  0,  -pan, 0);
 203         else if (kv == GDK_Down  || kv == GDK_j) gis_view_pan(self->view, -pan, 0,   0);
 204         else if (kv == GDK_Up    || kv == GDK_k) gis_view_pan(self->view,  pan, 0,   0);
 205         else if (kv == GDK_Right || kv == GDK_l) gis_view_pan(self->view,  0,   pan, 0);
 206         else if (kv == GDK_minus || kv == GDK_o) gis_view_zoom(self->view, 10./9);
 207         else if (kv == GDK_plus  || kv == GDK_i) gis_view_zoom(self->view, 9./10);
 208         else if (kv == GDK_H) gis_view_rotate(self->view,  0, 0, -2);
 209         else if (kv == GDK_J) gis_view_rotate(self->view,  2, 0,  0);
 210         else if (kv == GDK_K) gis_view_rotate(self->view, -2, 0,  0);
 211         else if (kv == GDK_L) gis_view_rotate(self->view,  0, 0,  2);
 212
 213         /* Testing */
 214         else if (kv == GDK_w) {self->wireframe = !self->wireframe; gtk_widget_queue_draw(GTK_WIDGET(self));}
 215 #ifdef ROAM_DEBUG
 216         else if (kv == GDK_n) roam_sphere_split_one(self->sphere);
 217         else if (kv == GDK_p) roam_sphere_merge_one(self->sphere);
 218         else if (kv == GDK_r) roam_sphere_split_merge(self->sphere);
 219         else if (kv == GDK_u) roam_sphere_update_errors(self->sphere);
 220         gdk_threads_enter();
 221         gtk_widget_queue_draw(GTK_WIDGET(self));
 222 #else
 223         gdk_threads_enter();
 224 #endif
 225
 226         return TRUE;
 227 }
 228
 229 static gboolean _update_errors_cb(gpointer sphere)
 230 {
 231         roam_sphere_update_errors(sphere);
 232         return FALSE;
 233 }
 234 static void on_view_changed(GisView *view,
 235                 gdouble _1, gdouble _2, gdouble _3, GisOpenGL *self)
 236 {
 237         g_debug("GisOpenGL: on_view_changed");
 238         gdk_threads_enter();
 239         gis_opengl_begin(self);
 240         set_visuals(self);
 241 #ifndef ROAM_DEBUG
 242         g_idle_add_full(G_PRIORITY_HIGH_IDLE+30, _update_errors_cb, self->sphere, NULL);
 243         //roam_sphere_update_errors(self->sphere);
 244 #endif
 245         gis_opengl_redraw(self);
 246         gis_opengl_end(self);
 247         gdk_threads_leave();
 248 }
 249
 250 static gboolean on_idle(GisOpenGL *self)
 251 {
 252         //g_debug("GisOpenGL: on_idle");
 253         gdk_threads_enter();
 254         gis_opengl_begin(self);
 255         if (roam_sphere_split_merge(self->sphere))
 256                 gis_opengl_redraw(self);
 257         gis_opengl_end(self);
 258         gdk_threads_leave();
 259         return TRUE;
 260 }
 261
 262
 263 /***********
 264  * Methods *
 265  ***********/
 266 GisOpenGL *gis_opengl_new(GisWorld *world, GisView *view, GisPlugins *plugins)
 267 {
 268         g_debug("GisOpenGL: new");
 269         GisOpenGL *self = g_object_new(GIS_TYPE_OPENGL, NULL);
 270         self->world   = world;
 271         self->view    = view;
 272         self->plugins = plugins;
 273         g_object_ref(world);
 274         g_object_ref(view);
 275
 276         g_signal_connect(self->view, "location-changed", G_CALLBACK(on_view_changed), self);
 277         g_signal_connect(self->view, "rotation-changed", G_CALLBACK(on_view_changed), self);
 278
 279         self->sphere = roam_sphere_new(self);
 280
 281         return g_object_ref(self);
 282 }
 283
 284 void gis_opengl_center_position(GisOpenGL *self, gdouble lat, gdouble lon, gdouble elev)
 285 {
 286         glRotatef(lon, 0, 1, 0);
 287         glRotatef(-lat, 1, 0, 0);
 288         glTranslatef(0, 0, elev2rad(elev));
 289 }
 290
 291 void gis_opengl_project(GisOpenGL *self,
 292                 gdouble lat, gdouble lon, gdouble elev,
 293                 gdouble *px, gdouble *py, gdouble *pz)
 294 {
 295         gdouble x, y, z;
 296         lle2xyz(lat, lon, elev, &x, &y, &z);
 297         gluProject(x, y, z,
 298                 self->sphere->view->model,
 299                 self->sphere->view->proj,
 300                 self->sphere->view->view,
 301                 px, py, pz);
 302 }
 303
 304 void gis_opengl_render_tile(GisOpenGL *self, GisTile *tile)
 305 {
 306         if (!tile || !tile->data)
 307                 return;
 308         GList *triangles = roam_sphere_get_intersect(self->sphere,
 309                         tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
 310         if (!triangles)
 311                 g_warning("GisOpenGL: render_tiles - No triangles to draw: edges=%f,%f,%f,%f",
 312                         tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
 313         //g_message("rendering %4d triangles for tile edges=%7.2f,%7.2f,%7.2f,%7.2f",
 314         //              g_list_length(triangles), tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
 315         for (GList *cur = triangles; cur; cur = cur->next) {
 316                 RoamTriangle *tri = cur->data;
 317
 318                 gdouble lat[3] = {tri->p.r->lat, tri->p.m->lat, tri->p.l->lat};
 319                 gdouble lon[3] = {tri->p.r->lon, tri->p.m->lon, tri->p.l->lon};
 320
 321                 if (lon[0] < -90 || lon[1] < -90 || lon[2] < -90) {
 322                         if (lon[0] > 90) lon[0] -= 360;
 323                         if (lon[1] > 90) lon[1] -= 360;
 324                         if (lon[2] > 90) lon[2] -= 360;
 325                 }
 326
 327                 gdouble n = tile->edge.n;
 328                 gdouble s = tile->edge.s;
 329                 gdouble e = tile->edge.e;
 330                 gdouble w = tile->edge.w;
 331
 332                 gdouble londist = e - w;
 333                 gdouble latdist = n - s;
 334
 335                 gdouble xy[][3] = {
 336                         {(lon[0]-w)/londist, 1-(lat[0]-s)/latdist},
 337                         {(lon[1]-w)/londist, 1-(lat[1]-s)/latdist},
 338                         {(lon[2]-w)/londist, 1-(lat[2]-s)/latdist},
 339                 };
 340
 341                 glEnable(GL_TEXTURE_2D);
 342                 glBindTexture(GL_TEXTURE_2D, *(guint*)tile->data);
 343                 glBegin(GL_TRIANGLES);
 344                 glNormal3dv(tri->p.r->norm); glTexCoord2dv(xy[0]); glVertex3dv((double*)tri->p.r);
 345                 glNormal3dv(tri->p.m->norm); glTexCoord2dv(xy[1]); glVertex3dv((double*)tri->p.m);
 346                 glNormal3dv(tri->p.l->norm); glTexCoord2dv(xy[2]); glVertex3dv((double*)tri->p.l);
 347                 glEnd();
 348         }
 349         g_list_free(triangles);
 350 }
 351
 352 void gis_opengl_render_tiles(GisOpenGL *opengl, GisTile *tile)
 353 {
 354         /* Only render children if possible */
 355         gboolean has_children = TRUE;
 356         GisTile *child;
 357         gis_tile_foreach(tile, child)
 358                 if (!child || !child->data)
 359                         has_children = FALSE;
 360         if (has_children)
 361                 /* Only render children */
 362                 gis_tile_foreach(tile, child)
 363                         gis_opengl_render_tiles(opengl, child);
 364         else
 365                 /* No children, render this tile */
 366                 gis_opengl_render_tile(opengl, tile);
 367 }
 368
 369 void gis_opengl_set_height_func(GisOpenGL *self, GisTile *tile,
 370                 RoamHeightFunc height_func, gpointer user_data, gboolean update)
 371 {
 372         if (!tile)
 373                 return;
 374         /* TODO: get points? */
 375         GList *triangles = roam_sphere_get_intersect(self->sphere,
 376                         tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
 377         for (GList *cur = triangles; cur; cur = cur->next) {
 378                 RoamTriangle *tri = cur->data;
 379                 RoamPoint *points[] = {tri->p.l, tri->p.m, tri->p.r, tri->split};
 380                 for (int i = 0; i < G_N_ELEMENTS(points); i++) {
 381                         points[i]->height_func = height_func;
 382                         points[i]->height_data = user_data;
 383                         roam_point_update_height(points[i]);
 384                 }
 385         }
 386         g_list_free(triangles);
 387 }
 388
 389 static void _gis_opengl_clear_height_func_rec(RoamTriangle *root)
 390 {
 391         if (!root)
 392                 return;
 393         RoamPoint *points[] = {root->p.l, root->p.m, root->p.r, root->split};
 394         for (int i = 0; i < G_N_ELEMENTS(points); i++) {
 395                 points[i]->height_func = NULL;
 396                 points[i]->height_data = NULL;
 397                 roam_point_update_height(points[i]);
 398         }
 399         _gis_opengl_clear_height_func_rec(root->kids[0]);
 400         _gis_opengl_clear_height_func_rec(root->kids[1]);
 401 }
 402 void gis_opengl_clear_height_func(GisOpenGL *self)
 403 {
 404         for (int i = 0; i < G_N_ELEMENTS(self->sphere->roots); i++)
 405                 _gis_opengl_clear_height_func_rec(self->sphere->roots[i]);
 406 }
 407
 408 void gis_opengl_redraw(GisOpenGL *self)
 409 {
 410         g_debug("GisOpenGL: redraw");
 411         gtk_widget_queue_draw(GTK_WIDGET(self));
 412 }
 413 void gis_opengl_begin(GisOpenGL *self)
 414 {
 415         g_assert(GIS_IS_OPENGL(self));
 416
 417         GdkGLContext   *glcontext  = gtk_widget_get_gl_context(GTK_WIDGET(self));
 418         GdkGLDrawable  *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
 419
 420         if (!gdk_gl_drawable_gl_begin(gldrawable, glcontext))
 421                 g_assert_not_reached();
 422 }
 423 void gis_opengl_end(GisOpenGL *self)
 424 {
 425         g_assert(GIS_IS_OPENGL(self));
 426         GdkGLDrawable *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
 427         gdk_gl_drawable_gl_end(gldrawable);
 428 }
 429 void gis_opengl_flush(GisOpenGL *self)
 430 {
 431         g_assert(GIS_IS_OPENGL(self));
 432         GdkGLDrawable *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
 433         if (gdk_gl_drawable_is_double_buffered(gldrawable))
 434                 gdk_gl_drawable_swap_buffers(gldrawable);
 435         else
 436                 glFlush();
 437         gdk_gl_drawable_gl_end(gldrawable);
 438 }
 439
 440
 441 /****************
 442  * GObject code *
 443  ****************/
 444 G_DEFINE_TYPE(GisOpenGL, gis_opengl, GTK_TYPE_DRAWING_AREA);
 445 static void gis_opengl_init(GisOpenGL *self)
 446 {
 447         g_debug("GisOpenGL: init");
 448         /* OpenGL setup */
 449         GdkGLConfig *glconfig = gdk_gl_config_new_by_mode(
 450                         GDK_GL_MODE_RGBA   | GDK_GL_MODE_DEPTH |
 451                         GDK_GL_MODE_DOUBLE | GDK_GL_MODE_ALPHA);
 452         if (!glconfig)
 453                 g_error("Failed to create glconfig");
 454         if (!gtk_widget_set_gl_capability(GTK_WIDGET(self),
 455                                 glconfig, NULL, TRUE, GDK_GL_RGBA_TYPE))
 456                 g_error("GL lacks required capabilities");
 457         g_object_unref(glconfig);
 458
 459         gtk_widget_set_size_request(GTK_WIDGET(self), 600, 550);
 460         gtk_widget_set_events(GTK_WIDGET(self),
 461                         GDK_BUTTON_PRESS_MASK |
 462                         GDK_ENTER_NOTIFY_MASK |
 463                         GDK_KEY_PRESS_MASK);
 464         g_object_set(self, "can-focus", TRUE, NULL);
 465
 466 #ifndef ROAM_DEBUG
 467         self->sm_source[0] = g_timeout_add_full(G_PRIORITY_HIGH_IDLE+30, 33,  (GSourceFunc)on_idle, self, NULL);
 468         self->sm_source[1] = g_timeout_add_full(G_PRIORITY_HIGH_IDLE+10, 500, (GSourceFunc)on_idle, self, NULL);
 469 #endif
 470
 471         g_signal_connect(self, "realize",            G_CALLBACK(on_realize),      NULL);
 472         g_signal_connect(self, "configure-event",    G_CALLBACK(on_configure),    NULL);
 473         g_signal_connect(self, "expose-event",       G_CALLBACK(on_expose),       NULL);
 474
 475         g_signal_connect(self, "button-press-event", G_CALLBACK(on_button_press), NULL);
 476         g_signal_connect(self, "enter-notify-event", G_CALLBACK(on_button_press), NULL);
 477         g_signal_connect(self, "key-press-event",    G_CALLBACK(on_key_press),    NULL);
 478 }
 479 static GObject *gis_opengl_constructor(GType gtype, guint n_properties,
 480                 GObjectConstructParam *properties)
 481 {
 482         g_debug("GisOpengl: constructor");
 483         GObjectClass *parent_class = G_OBJECT_CLASS(gis_opengl_parent_class);
 484         return parent_class->constructor(gtype, n_properties, properties);
 485 }
 486 static void gis_opengl_dispose(GObject *_self)
 487 {
 488         g_debug("GisOpenGL: dispose");
 489         GisOpenGL *self = GIS_OPENGL(_self);
 490         if (self->sm_source[0]) {
 491                 g_source_remove(self->sm_source[0]);
 492                 self->sm_source[0] = 0;
 493         }
 494         if (self->sm_source[1]) {
 495                 g_source_remove(self->sm_source[1]);
 496                 self->sm_source[1] = 0;
 497         }
 498         if (self->sphere) {
 499                 roam_sphere_free(self->sphere);
 500                 self->sphere = NULL;
 501         }
 502         if (self->world) {
 503                 g_object_unref(self->world);
 504                 self->world = NULL;
 505         }
 506         if (self->view) {
 507                 g_object_unref(self->view);
 508                 self->view = NULL;
 509         }
 510         G_OBJECT_CLASS(gis_opengl_parent_class)->dispose(_self);
 511 }
 512 static void gis_opengl_finalize(GObject *_self)
 513 {
 514         g_debug("GisOpenGL: finalize");
 515         GisOpenGL *self = GIS_OPENGL(_self);
 516         G_OBJECT_CLASS(gis_opengl_parent_class)->finalize(_self);
 517 }
 518 static void gis_opengl_class_init(GisOpenGLClass *klass)
 519 {
 520         g_debug("GisOpenGL: class_init");
 521         GObjectClass *gobject_class = G_OBJECT_CLASS(klass);
 522         gobject_class->constructor  = gis_opengl_constructor;
 523         gobject_class->dispose      = gis_opengl_dispose;
 524         gobject_class->finalize     = gis_opengl_finalize;
 525 }