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 <string.h>
  25 #include <gdk/gdkkeysyms.h>
  26 #include <gtk/gtk.h>
  27 #include <gtk/gtkgl.h>
  28 #include <GL/gl.h>
  29 #include <GL/glu.h>
  30
  31 #include "gis-opengl.h"
  32 #include "gis-util.h"
  33 #include "gis-object.h"
  34 #include "roam.h"
  35
  36 #define FOV_DIST   2000.0
  37 #define MPPX(dist) (4*dist/FOV_DIST)
  38
  39 // #define ROAM_DEBUG
  40
  41 /***********
  42  * Helpers *
  43  ***********/
  44 static void _gis_opengl_begin(GisOpenGL *self)
  45 {
  46         g_assert(GIS_IS_OPENGL(self));
  47
  48         GdkGLContext   *glcontext  = gtk_widget_get_gl_context(GTK_WIDGET(self));
  49         GdkGLDrawable  *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
  50
  51         if (!gdk_gl_drawable_gl_begin(gldrawable, glcontext))
  52                 g_assert_not_reached();
  53 }
  54
  55 static void _gis_opengl_end(GisOpenGL *self)
  56 {
  57         g_assert(GIS_IS_OPENGL(self));
  58         GdkGLDrawable *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
  59         gdk_gl_drawable_gl_end(gldrawable);
  60 }
  61
  62
  63 /********************
  64  * Object handleing *
  65  ********************/
  66 static void _draw_marker(GisOpenGL *self, GisMarker *marker)
  67 {
  68         GisPoint *point = gis_object_center(GIS_OBJECT(marker));
  69         gdouble px, py, pz;
  70         gis_viewer_project(GIS_VIEWER(self),
  71                         point->lat, point->lon, point->elev,
  72                         &px, &py, &pz);
  73
  74         g_debug("GisOpenGL: draw_marker - texture=%d", marker->tex);
  75
  76         cairo_surface_t *surface = cairo_get_target(marker->cairo);
  77         gdouble width  = cairo_image_surface_get_width(surface);
  78         gdouble height = cairo_image_surface_get_height(surface);
  79
  80         glMatrixMode(GL_PROJECTION); glLoadIdentity();
  81         glMatrixMode(GL_MODELVIEW);  glLoadIdentity();
  82         glOrtho(0, GTK_WIDGET(self)->allocation.width,
  83                 0, GTK_WIDGET(self)->allocation.height, -1, 1);
  84         glTranslated(px - marker->xoff,
  85                      py - marker->yoff, 0);
  86
  87         glDisable(GL_LIGHTING);
  88         glDisable(GL_COLOR_MATERIAL);
  89         glDisable(GL_DEPTH_TEST);
  90         glEnable(GL_TEXTURE_2D);
  91         glBindTexture(GL_TEXTURE_2D, marker->tex);
  92         g_debug("bind_texture: %d", marker->tex);
  93         glBegin(GL_QUADS);
  94         glTexCoord2f(1, 1); glVertex3f(width, 0     , 0);
  95         glTexCoord2f(1, 0); glVertex3f(width, height, 0);
  96         glTexCoord2f(0, 0); glVertex3f(0    , height, 0);
  97         glTexCoord2f(0, 1); glVertex3f(0    , 0     , 0);
  98         glEnd();
  99 }
 100
 101 static void _draw_callback(GisOpenGL *self, GisCallback *callback)
 102 {
 103         callback->callback(callback, callback->user_data);
 104 }
 105
 106 static void _draw_object(GisOpenGL *self, GisObject *object)
 107 {
 108         //g_debug("GisOpenGL: draw_object");
 109         /* Skip out of range objects */
 110         if (object->lod > 0) {
 111                 gdouble eye[3], obj[3];
 112                 gis_viewer_get_location(GIS_VIEWER(self), &eye[0], &eye[1], &eye[2]);
 113                 lle2xyz(eye[0], eye[1], eye[2], &eye[0], &eye[1], &eye[2]);
 114                 lle2xyz(object->center.lat, object->center.lon, object->center.elev,
 115                         &obj[0], &obj[1], &obj[2]);
 116                 gdouble dist = distd(obj, eye);
 117                 if (object->lod < dist)
 118                         return;
 119         }
 120
 121         /* Draw */
 122         glMatrixMode(GL_PROJECTION); glPushMatrix();
 123         glMatrixMode(GL_MODELVIEW);  glPushMatrix();
 124         glPushAttrib(GL_ALL_ATTRIB_BITS);
 125         switch (object->type) {
 126         case GIS_TYPE_MARKER:
 127                 _draw_marker(self, GIS_MARKER(object));
 128                 break;
 129         case GIS_TYPE_CALLBACK:
 130                 _draw_callback(self, GIS_CALLBACK(object));
 131                 break;
 132         default:
 133                 break;
 134         }
 135         glPopAttrib();
 136         glMatrixMode(GL_PROJECTION); glPopMatrix();
 137         glMatrixMode(GL_MODELVIEW);  glPopMatrix();
 138 }
 139
 140 static void _load_object(GisOpenGL *self, GisObject *object)
 141 {
 142         g_debug("GisOpenGL: load_object");
 143         switch (object->type) {
 144         case GIS_TYPE_MARKER: {
 145                 GisMarker *marker = GIS_MARKER(object);
 146                 cairo_surface_t *surface = cairo_get_target(marker->cairo);
 147                 gdouble width  = cairo_image_surface_get_width(surface);
 148                 gdouble height = cairo_image_surface_get_height(surface);
 149
 150                 _gis_opengl_begin(self);
 151                 glEnable(GL_TEXTURE_2D);
 152                 glGenTextures(1, &marker->tex);
 153                 glBindTexture(GL_TEXTURE_2D, marker->tex);
 154
 155                 glPixelStorei(GL_UNPACK_ALIGNMENT, 1);
 156                 glPixelStorei(GL_PACK_ALIGNMENT, 1);
 157                 glTexImage2D(GL_TEXTURE_2D, 0, 4, width, height, 0, GL_RGBA, GL_UNSIGNED_BYTE,
 158                                 cairo_image_surface_get_data(surface));
 159                 glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
 160                 glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
 161                 g_debug("load_texture: %d", marker->tex);
 162                 _gis_opengl_end(self);
 163                 break;
 164         }
 165         default:
 166                 break;
 167         }
 168 }
 169
 170 static void _unload_object(GisOpenGL *self, GisObject *object)
 171 {
 172         g_debug("GisOpenGL: unload_object");
 173         switch (object->type) {
 174         case GIS_TYPE_MARKER: {
 175                 GisMarker *marker = GIS_MARKER(object);
 176                 g_debug("delete_texture: %d", marker->tex);
 177                 glDeleteTextures(1, &marker->tex);
 178                 break;
 179         }
 180         default:
 181                 break;
 182         }
 183 }
 184
 185 /*************
 186  * ROAM Code *
 187  *************/
 188 static void _set_visuals(GisOpenGL *self)
 189 {
 190         glMatrixMode(GL_MODELVIEW);
 191         glLoadIdentity();
 192
 193         /* Camera 1 */
 194         double lat, lon, elev, rx, ry, rz;
 195         gis_viewer_get_location(GIS_VIEWER(self), &lat, &lon, &elev);
 196         gis_viewer_get_rotation(GIS_VIEWER(self), &rx, &ry, &rz);
 197         glRotatef(rx, 1, 0, 0);
 198         glRotatef(rz, 0, 0, 1);
 199
 200         /* Lighting */
 201 #ifdef ROAM_DEBUG
 202         float light_ambient[]  = {0.7f, 0.7f, 0.7f, 1.0f};
 203         float light_diffuse[]  = {2.0f, 2.0f, 2.0f, 1.0f};
 204 #else
 205         float light_ambient[]  = {0.2f, 0.2f, 0.2f, 1.0f};
 206         float light_diffuse[]  = {5.0f, 5.0f, 5.0f, 1.0f};
 207 #endif
 208         float light_position[] = {-13*EARTH_R, 1*EARTH_R, 3*EARTH_R, 1.0f};
 209         glLightfv(GL_LIGHT0, GL_AMBIENT,  light_ambient);
 210         glLightfv(GL_LIGHT0, GL_DIFFUSE,  light_diffuse);
 211         glLightfv(GL_LIGHT0, GL_POSITION, light_position);
 212         glEnable(GL_LIGHT0);
 213         glEnable(GL_LIGHTING);
 214
 215         float material_ambient[]  = {0.2, 0.2, 0.2, 1.0};
 216         float material_diffuse[]  = {0.8, 0.8, 0.8, 1.0};
 217         float material_specular[] = {0.0, 0.0, 0.0, 1.0};
 218         float material_emission[] = {0.0, 0.0, 0.0, 1.0};
 219         glMaterialfv(GL_FRONT_AND_BACK, GL_AMBIENT,  material_ambient);
 220         glMaterialfv(GL_FRONT_AND_BACK, GL_DIFFUSE,  material_diffuse);
 221         glMaterialfv(GL_FRONT_AND_BACK, GL_SPECULAR, material_specular);
 222         glMaterialfv(GL_FRONT_AND_BACK, GL_EMISSION, material_emission);
 223         glDisable(GL_TEXTURE_2D);
 224         glDisable(GL_COLOR_MATERIAL);
 225
 226         /* Camera 2 */
 227         glTranslatef(0, 0, -elev2rad(elev));
 228         glRotatef(lat, 1, 0, 0);
 229         glRotatef(-lon, 0, 1, 0);
 230
 231         /* Misc */
 232         gdouble rg   = MAX(0, 1-(elev/20000));
 233         gdouble blue = MAX(0, 1-(elev/50000));
 234         glClearColor(MIN(0.65,rg), MIN(0.65,rg), MIN(1,blue), 1.0f);
 235         glColor4f(1, 1, 1, 1);
 236
 237         glDisable(GL_ALPHA_TEST);
 238
 239         glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
 240         glEnable(GL_BLEND);
 241
 242 #ifndef ROAM_DEBUG
 243         glCullFace(GL_BACK);
 244         glEnable(GL_CULL_FACE);
 245 #endif
 246
 247         glClearDepth(1.0);
 248         glDepthFunc(GL_LEQUAL);
 249         glEnable(GL_DEPTH_TEST);
 250
 251         glEnable(GL_LINE_SMOOTH);
 252
 253         glPolygonMode(GL_FRONT_AND_BACK, GL_FILL);
 254         //glShadeModel(GL_FLAT);
 255
 256         roam_sphere_update_view(self->sphere);
 257 }
 258
 259
 260 /*************
 261  * Callbacks *
 262  *************/
 263 /* The unsorted/sroted GLists are blank head nodes,
 264  * This way us we can remove objects from the level just by fixing up links
 265  * I.e. we don't need to do a lookup to remove an object if we have its GList */
 266 struct RenderLevel {
 267         GList unsorted;
 268         GList sorted;
 269 };
 270
 271 static void on_realize(GisOpenGL *self, gpointer _)
 272 {
 273         g_debug("GisOpenGL: on_realize");
 274         _set_visuals(self);
 275         roam_sphere_update_errors(self->sphere);
 276 }
 277 static gboolean on_configure(GisOpenGL *self, GdkEventConfigure *event, gpointer _)
 278 {
 279         g_debug("GisOpenGL: on_configure");
 280         _gis_opengl_begin(self);
 281
 282         double width  = GTK_WIDGET(self)->allocation.width;
 283         double height = GTK_WIDGET(self)->allocation.height;
 284
 285         /* Setup OpenGL Window */
 286         glViewport(0, 0, width, height);
 287         glMatrixMode(GL_PROJECTION);
 288         glLoadIdentity();
 289         double ang = atan(height/FOV_DIST);
 290         gluPerspective(rad2deg(ang)*2, width/height, 1, 20*EARTH_R);
 291
 292 #ifndef ROAM_DEBUG
 293         roam_sphere_update_errors(self->sphere);
 294 #endif
 295
 296         _gis_opengl_end(self);
 297         return FALSE;
 298 }
 299
 300 static gboolean _draw_level(gpointer key, gpointer value, gpointer user_data)
 301 {
 302         g_debug("GisOpenGL: _draw_level - level=%-4d", (int)key);
 303         GisOpenGL *self = user_data;
 304         struct RenderLevel *level = value;
 305         int nsorted = 0, nunsorted = 0;
 306         GList *cur = NULL;
 307
 308         /* Draw opaque objects without sorting */
 309         glDepthMask(TRUE);
 310         glClear(GL_DEPTH_BUFFER_BIT);
 311         for (cur = level->unsorted.next; cur; cur = cur->next, nunsorted++)
 312                 _draw_object(self, GIS_OBJECT(cur->data));
 313
 314         /* Freeze depth buffer and draw transparent objects sorted */
 315         /* TODO: sorting */
 316         glDepthMask(FALSE);
 317         for (cur = level->sorted.next; cur; cur = cur->next, nsorted++)
 318                 _draw_object(self, GIS_OBJECT(cur->data));
 319
 320         /* TODO: Prune empty levels */
 321
 322         g_debug("GisOpenGL: _draw_level - drew %d,%d objects",
 323                         nunsorted, nsorted);
 324         return FALSE;
 325 }
 326
 327 static gboolean on_expose(GisOpenGL *self, GdkEventExpose *event, gpointer _)
 328 {
 329         g_debug("GisOpenGL: on_expose - begin");
 330         _gis_opengl_begin(self);
 331
 332         glClear(GL_COLOR_BUFFER_BIT);
 333
 334 #ifndef ROAM_DEBUG
 335         g_tree_foreach(self->objects, _draw_level, self);
 336         if (self->wireframe) {
 337                 _set_visuals(self);
 338                 glPolygonMode(GL_FRONT_AND_BACK, GL_LINE);
 339                 roam_sphere_draw(self->sphere);
 340         }
 341 #else
 342         _set_visuals(self);
 343         glColor4f(0.0, 0.0, 9.0, 0.6);
 344         glDisable(GL_TEXTURE_2D);
 345         glPolygonMode(GL_FRONT_AND_BACK, GL_LINE);
 346         roam_sphere_draw(self->sphere);
 347         //roam_sphere_draw_normals(self->sphere);
 348 #endif
 349
 350         GdkGLDrawable *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
 351         gdk_gl_drawable_swap_buffers(gldrawable);
 352
 353         _gis_opengl_end(self);
 354         g_debug("GisOpenGL: on_expose - end\n");
 355         return FALSE;
 356 }
 357
 358 static gboolean on_key_press(GisOpenGL *self, GdkEventKey *event, gpointer _)
 359 {
 360         g_debug("GisOpenGL: on_key_press - key=%x, state=%x, plus=%x",
 361                         event->keyval, event->state, GDK_plus);
 362
 363         guint kv = event->keyval;
 364         gdk_threads_leave();
 365         /* Testing */
 366         if (kv == GDK_w) {
 367                 self->wireframe = !self->wireframe;
 368                 gtk_widget_queue_draw(GTK_WIDGET(self));
 369         }
 370 #ifdef ROAM_DEBUG
 371         else if (kv == GDK_n) roam_sphere_split_one(self->sphere);
 372         else if (kv == GDK_p) roam_sphere_merge_one(self->sphere);
 373         else if (kv == GDK_r) roam_sphere_split_merge(self->sphere);
 374         else if (kv == GDK_u) roam_sphere_update_errors(self->sphere);
 375         gdk_threads_enter();
 376         gtk_widget_queue_draw(GTK_WIDGET(self));
 377 #else
 378         gdk_threads_enter();
 379 #endif
 380         return FALSE;
 381 }
 382
 383 static gboolean _update_errors_cb(gpointer sphere)
 384 {
 385         roam_sphere_update_errors(sphere);
 386         return FALSE;
 387 }
 388 static void on_view_changed(GisOpenGL *self,
 389                 gdouble _1, gdouble _2, gdouble _3)
 390 {
 391         g_debug("GisOpenGL: on_view_changed");
 392         gdk_threads_enter();
 393         _gis_opengl_begin(self);
 394         _set_visuals(self);
 395 #ifndef ROAM_DEBUG
 396         g_idle_add_full(G_PRIORITY_HIGH_IDLE+30, _update_errors_cb, self->sphere, NULL);
 397         //roam_sphere_update_errors(self->sphere);
 398 #endif
 399         _gis_opengl_end(self);
 400         gdk_threads_leave();
 401 }
 402
 403 static gboolean on_idle(GisOpenGL *self)
 404 {
 405         //g_debug("GisOpenGL: on_idle");
 406         gdk_threads_enter();
 407         _gis_opengl_begin(self);
 408         if (roam_sphere_split_merge(self->sphere))
 409                 gtk_widget_queue_draw(GTK_WIDGET(self));
 410         _gis_opengl_end(self);
 411         gdk_threads_leave();
 412         return TRUE;
 413 }
 414
 415
 416 /*********************
 417  * GisViewer methods *
 418  *********************/
 419 GisViewer *gis_opengl_new(GisPlugins *plugins)
 420 {
 421         g_debug("GisOpenGL: new");
 422         GisViewer *self = g_object_new(GIS_TYPE_OPENGL, NULL);
 423         self->plugins = plugins;
 424         return self;
 425 }
 426
 427 static void gis_opengl_center_position(GisViewer *_self, gdouble lat, gdouble lon, gdouble elev)
 428 {
 429         GisOpenGL *self = GIS_OPENGL(_self);
 430         glRotatef(lon, 0, 1, 0);
 431         glRotatef(-lat, 1, 0, 0);
 432         glTranslatef(0, 0, elev2rad(elev));
 433 }
 434
 435 static void gis_opengl_project(GisViewer *_self,
 436                 gdouble lat, gdouble lon, gdouble elev,
 437                 gdouble *px, gdouble *py, gdouble *pz)
 438 {
 439         GisOpenGL *self = GIS_OPENGL(_self);
 440         gdouble x, y, z;
 441         lle2xyz(lat, lon, elev, &x, &y, &z);
 442         gluProject(x, y, z,
 443                 self->sphere->view->model,
 444                 self->sphere->view->proj,
 445                 self->sphere->view->view,
 446                 px, py, pz);
 447 }
 448
 449 static void gis_opengl_render_tile(GisViewer *_self, GisTile *tile)
 450 {
 451         GisOpenGL *self = GIS_OPENGL(_self);
 452         if (!tile || !tile->data)
 453                 return;
 454         GList *triangles = roam_sphere_get_intersect(self->sphere,
 455                         tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
 456         if (!triangles)
 457                 g_warning("GisOpenGL: render_tiles - No triangles to draw: edges=%f,%f,%f,%f",
 458                         tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
 459         //g_message("rendering %4d triangles for tile edges=%7.2f,%7.2f,%7.2f,%7.2f",
 460         //              g_list_length(triangles), tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
 461         for (GList *cur = triangles; cur; cur = cur->next) {
 462                 RoamTriangle *tri = cur->data;
 463
 464                 gdouble lat[3] = {tri->p.r->lat, tri->p.m->lat, tri->p.l->lat};
 465                 gdouble lon[3] = {tri->p.r->lon, tri->p.m->lon, tri->p.l->lon};
 466
 467                 if (lon[0] < -90 || lon[1] < -90 || lon[2] < -90) {
 468                         if (lon[0] > 90) lon[0] -= 360;
 469                         if (lon[1] > 90) lon[1] -= 360;
 470                         if (lon[2] > 90) lon[2] -= 360;
 471                 }
 472
 473                 gdouble n = tile->edge.n;
 474                 gdouble s = tile->edge.s;
 475                 gdouble e = tile->edge.e;
 476                 gdouble w = tile->edge.w;
 477
 478                 gdouble londist = e - w;
 479                 gdouble latdist = n - s;
 480
 481                 gdouble xy[][3] = {
 482                         {(lon[0]-w)/londist, 1-(lat[0]-s)/latdist},
 483                         {(lon[1]-w)/londist, 1-(lat[1]-s)/latdist},
 484                         {(lon[2]-w)/londist, 1-(lat[2]-s)/latdist},
 485                 };
 486
 487                 glEnable(GL_TEXTURE_2D);
 488                 glBindTexture(GL_TEXTURE_2D, *(guint*)tile->data);
 489                 glBegin(GL_TRIANGLES);
 490                 glNormal3dv(tri->p.r->norm); glTexCoord2dv(xy[0]); glVertex3dv((double*)tri->p.r);
 491                 glNormal3dv(tri->p.m->norm); glTexCoord2dv(xy[1]); glVertex3dv((double*)tri->p.m);
 492                 glNormal3dv(tri->p.l->norm); glTexCoord2dv(xy[2]); glVertex3dv((double*)tri->p.l);
 493                 glEnd();
 494         }
 495         g_list_free(triangles);
 496 }
 497
 498 static void gis_opengl_render_tiles(GisViewer *_self, GisTile *tile)
 499 {
 500         GisOpenGL *self = GIS_OPENGL(_self);
 501         /* Only render children if possible */
 502         gboolean has_children = TRUE;
 503         GisTile *child;
 504         gis_tile_foreach(tile, child)
 505                 if (!child || !child->data)
 506                         has_children = FALSE;
 507         if (has_children)
 508                 /* Only render children */
 509                 gis_tile_foreach(tile, child)
 510                         gis_opengl_render_tiles(_self, child);
 511         else
 512                 /* No children, render this tile */
 513                 gis_opengl_render_tile(_self, tile);
 514 }
 515
 516 static void gis_opengl_set_height_func(GisViewer *_self, GisTile *tile,
 517                 RoamHeightFunc height_func, gpointer user_data, gboolean update)
 518 {
 519         GisOpenGL *self = GIS_OPENGL(_self);
 520         if (!tile)
 521                 return;
 522         /* TODO: get points? */
 523         GList *triangles = roam_sphere_get_intersect(self->sphere,
 524                         tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
 525         for (GList *cur = triangles; cur; cur = cur->next) {
 526                 RoamTriangle *tri = cur->data;
 527                 RoamPoint *points[] = {tri->p.l, tri->p.m, tri->p.r, tri->split};
 528                 for (int i = 0; i < G_N_ELEMENTS(points); i++) {
 529                         points[i]->height_func = height_func;
 530                         points[i]->height_data = user_data;
 531                         roam_point_update_height(points[i]);
 532                 }
 533         }
 534         g_list_free(triangles);
 535 }
 536
 537 static void _gis_opengl_clear_height_func_rec(RoamTriangle *root)
 538 {
 539         if (!root)
 540                 return;
 541         RoamPoint *points[] = {root->p.l, root->p.m, root->p.r, root->split};
 542         for (int i = 0; i < G_N_ELEMENTS(points); i++) {
 543                 points[i]->height_func = NULL;
 544                 points[i]->height_data = NULL;
 545                 roam_point_update_height(points[i]);
 546         }
 547         _gis_opengl_clear_height_func_rec(root->kids[0]);
 548         _gis_opengl_clear_height_func_rec(root->kids[1]);
 549 }
 550
 551 static void gis_opengl_clear_height_func(GisViewer *_self)
 552 {
 553         GisOpenGL *self = GIS_OPENGL(_self);
 554         for (int i = 0; i < G_N_ELEMENTS(self->sphere->roots); i++)
 555                 _gis_opengl_clear_height_func_rec(self->sphere->roots[i]);
 556 }
 557
 558 static void gis_opengl_begin(GisViewer *_self)
 559 {
 560         g_assert(GIS_IS_OPENGL(_self));
 561         _gis_opengl_begin(GIS_OPENGL(_self));
 562 }
 563
 564 static void gis_opengl_end(GisViewer *_self)
 565 {
 566         g_assert(GIS_IS_OPENGL(_self));
 567         _gis_opengl_end(GIS_OPENGL(_self));
 568 }
 569
 570 static gpointer gis_opengl_add(GisViewer *_self, GisObject *object,
 571                 gint key, gboolean sort)
 572 {
 573         g_assert(GIS_IS_OPENGL(_self));
 574         GisOpenGL *self = GIS_OPENGL(_self);
 575         _load_object(self, object);
 576         struct RenderLevel *level = g_tree_lookup(self->objects, (gpointer)key);
 577         if (!level) {
 578                 level = g_new0(struct RenderLevel, 1);
 579                 g_tree_insert(self->objects, (gpointer)key, level);
 580         }
 581         GList *list = sort ? &level->sorted : &level->unsorted;
 582         list->next = g_list_prepend(list->next, object);
 583         return list->next;
 584 }
 585
 586 static void gis_opengl_remove(GisViewer *_self, gpointer _link)
 587 {
 588         g_assert(GIS_IS_OPENGL(_self));
 589         GList *link = _link;
 590         GisOpenGL *self = GIS_OPENGL(_self);
 591         _unload_object(self, link->data);
 592         /* Just unlink and free it (blowup link to avoid warnings) */
 593         link = g_list_delete_link(NULL, link);
 594 }
 595
 596 /****************
 597  * GObject code *
 598  ****************/
 599 static int _objects_cmp(gconstpointer _a, gconstpointer _b)
 600 {
 601         gint a = (int)_a, b = (int)_b;
 602         return a < b ? -1 :
 603                a > b ?  1 : 0;
 604 }
 605
 606 G_DEFINE_TYPE(GisOpenGL, gis_opengl, GIS_TYPE_VIEWER);
 607 static void gis_opengl_init(GisOpenGL *self)
 608 {
 609         g_debug("GisOpenGL: init");
 610         /* OpenGL setup */
 611         GdkGLConfig *glconfig = gdk_gl_config_new_by_mode(
 612                         GDK_GL_MODE_RGBA   | GDK_GL_MODE_DEPTH |
 613                         GDK_GL_MODE_DOUBLE | GDK_GL_MODE_ALPHA);
 614         if (!glconfig)
 615                 g_error("Failed to create glconfig");
 616         if (!gtk_widget_set_gl_capability(GTK_WIDGET(self),
 617                                 glconfig, NULL, TRUE, GDK_GL_RGBA_TYPE))
 618                 g_error("GL lacks required capabilities");
 619         g_object_unref(glconfig);
 620
 621         gtk_widget_set_size_request(GTK_WIDGET(self), 600, 550);
 622         gtk_widget_set_events(GTK_WIDGET(self),
 623                         GDK_BUTTON_PRESS_MASK |
 624                         GDK_ENTER_NOTIFY_MASK |
 625                         GDK_KEY_PRESS_MASK);
 626         g_object_set(self, "can-focus", TRUE, NULL);
 627
 628         self->objects = g_tree_new(_objects_cmp);
 629         self->sphere = roam_sphere_new(self);
 630
 631 #ifndef ROAM_DEBUG
 632         self->sm_source[0] = g_timeout_add_full(G_PRIORITY_HIGH_IDLE+30, 33,  (GSourceFunc)on_idle, self, NULL);
 633         self->sm_source[1] = g_timeout_add_full(G_PRIORITY_HIGH_IDLE+10, 500, (GSourceFunc)on_idle, self, NULL);
 634 #endif
 635
 636         g_signal_connect(self, "realize",          G_CALLBACK(on_realize),      NULL);
 637         g_signal_connect(self, "configure-event",  G_CALLBACK(on_configure),    NULL);
 638         g_signal_connect(self, "expose-event",     G_CALLBACK(on_expose),       NULL);
 639
 640         g_signal_connect(self, "key-press-event",  G_CALLBACK(on_key_press),    NULL);
 641
 642         g_signal_connect(self, "location-changed", G_CALLBACK(on_view_changed), NULL);
 643         g_signal_connect(self, "rotation-changed", G_CALLBACK(on_view_changed), NULL);
 644 }
 645 static void gis_opengl_dispose(GObject *_self)
 646 {
 647         g_debug("GisOpenGL: dispose");
 648         GisOpenGL *self = GIS_OPENGL(_self);
 649         if (self->sm_source[0]) {
 650                 g_source_remove(self->sm_source[0]);
 651                 self->sm_source[0] = 0;
 652         }
 653         if (self->sm_source[1]) {
 654                 g_source_remove(self->sm_source[1]);
 655                 self->sm_source[1] = 0;
 656         }
 657         if (self->sphere) {
 658                 roam_sphere_free(self->sphere);
 659                 self->sphere = NULL;
 660         }
 661         /* TODO: Cleanup/free objects tree */
 662         G_OBJECT_CLASS(gis_opengl_parent_class)->dispose(_self);
 663 }
 664 static void gis_opengl_finalize(GObject *_self)
 665 {
 666         g_debug("GisViewer: finalize");
 667         GisOpenGL *self = GIS_OPENGL(_self);
 668         G_OBJECT_CLASS(gis_opengl_parent_class)->finalize(_self);
 669 }
 670 static void gis_opengl_class_init(GisOpenGLClass *klass)
 671 {
 672         g_debug("GisOpenGL: class_init");
 673         GObjectClass *gobject_class = G_OBJECT_CLASS(klass);
 674         gobject_class->dispose     = gis_opengl_dispose;
 675
 676         GisViewerClass *viewer_class = GIS_VIEWER_CLASS(klass);
 677         viewer_class->center_position   = gis_opengl_center_position;
 678         viewer_class->project           = gis_opengl_project;
 679         viewer_class->clear_height_func = gis_opengl_clear_height_func;
 680         viewer_class->set_height_func   = gis_opengl_set_height_func;
 681         viewer_class->render_tile       = gis_opengl_render_tile;
 682         viewer_class->render_tiles      = gis_opengl_render_tiles;
 683         viewer_class->begin             = gis_opengl_begin;
 684         viewer_class->end               = gis_opengl_end;
 685         viewer_class->add               = gis_opengl_add;
 686         viewer_class->remove            = gis_opengl_remove;
 687 }