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

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