2 * Copyright (C) 2009 Andy Spencer <spenceal@rose-hulman.edu>
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.
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.
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/>.
18 /* Tessellation, "finding intersecting triangles" */
19 /* http://research.microsoft.com/pubs/70307/tr-2006-81.pdf */
20 /* http://www.opengl.org/wiki/Alpha_Blending */
25 #include <gdk/gdkkeysyms.h>
27 #include <gtk/gtkgl.h>
31 #include "gis-opengl.h"
33 #include "gis-object.h"
36 #define FOV_DIST 2000.0
37 #define MPPX(dist) (4*dist/FOV_DIST)
44 static void _gis_opengl_begin(GisOpenGL *self)
46 g_assert(GIS_IS_OPENGL(self));
48 GdkGLContext *glcontext = gtk_widget_get_gl_context(GTK_WIDGET(self));
49 GdkGLDrawable *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
51 if (!gdk_gl_drawable_gl_begin(gldrawable, glcontext))
52 g_assert_not_reached();
55 static void _gis_opengl_end(GisOpenGL *self)
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);
66 static void _draw_marker(GisOpenGL *self, GisMarker *marker)
68 GisPoint *point = gis_object_center(GIS_OBJECT(marker));
70 gis_viewer_project(GIS_VIEWER(self),
71 point->lat, point->lon, point->elev,
74 g_debug("GisOpenGL: draw_marker - texture=%d", marker->tex);
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);
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);
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);
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);
102 static void _draw_objects(GisOpenGL *self)
104 g_debug("GisOpenGL: draw_objects");
106 for (GList *cur = self->objects; cur; cur = cur->next) {
107 glMatrixMode(GL_PROJECTION); glPushMatrix();
108 glMatrixMode(GL_MODELVIEW); glPushMatrix();
109 GisObject *object = cur->data;
110 switch (object->type) {
111 case GIS_TYPE_MARKER:
112 _draw_marker(self, GIS_MARKER(object));
117 glMatrixMode(GL_PROJECTION); glPopMatrix();
118 glMatrixMode(GL_MODELVIEW); glPopMatrix();
122 static void _load_object(GisOpenGL *self, GisObject *object)
124 g_debug("GisOpenGL: load_object");
125 switch (object->type) {
126 case GIS_TYPE_MARKER: {
127 GisMarker *marker = GIS_MARKER(object);
128 cairo_surface_t *surface = cairo_get_target(marker->cairo);
129 gdouble width = cairo_image_surface_get_width(surface);
130 gdouble height = cairo_image_surface_get_height(surface);
132 _gis_opengl_begin(self);
133 glEnable(GL_TEXTURE_2D);
134 glGenTextures(1, &marker->tex);
135 glBindTexture(GL_TEXTURE_2D, marker->tex);
137 glPixelStorei(GL_UNPACK_ALIGNMENT, 1);
138 glPixelStorei(GL_PACK_ALIGNMENT, 1);
139 glTexImage2D(GL_TEXTURE_2D, 0, 4, width, height, 0, GL_RGBA, GL_UNSIGNED_BYTE,
140 cairo_image_surface_get_data(surface));
141 glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
142 glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
143 g_debug("load_texture: %d", marker->tex);
144 _gis_opengl_end(self);
152 static void _free_object(GisOpenGL *self, GisObject *object)
154 g_debug("GisOpenGL: free_object");
155 switch (object->type) {
156 case GIS_TYPE_MARKER: {
157 GisMarker *marker = GIS_MARKER(object);
158 g_debug("delete_texture: %d", marker->tex);
159 glDeleteTextures(1, &marker->tex);
170 static void _set_visuals(GisOpenGL *self)
172 glMatrixMode(GL_MODELVIEW);
176 double lat, lon, elev, rx, ry, rz;
177 gis_viewer_get_location(GIS_VIEWER(self), &lat, &lon, &elev);
178 gis_viewer_get_rotation(GIS_VIEWER(self), &rx, &ry, &rz);
179 glRotatef(rx, 1, 0, 0);
180 glRotatef(rz, 0, 0, 1);
184 float light_ambient[] = {0.7f, 0.7f, 0.7f, 1.0f};
185 float light_diffuse[] = {2.0f, 2.0f, 2.0f, 1.0f};
187 float light_ambient[] = {0.2f, 0.2f, 0.2f, 1.0f};
188 float light_diffuse[] = {5.0f, 5.0f, 5.0f, 1.0f};
190 float light_position[] = {-13*EARTH_R, 1*EARTH_R, 3*EARTH_R, 1.0f};
191 glLightfv(GL_LIGHT0, GL_AMBIENT, light_ambient);
192 glLightfv(GL_LIGHT0, GL_DIFFUSE, light_diffuse);
193 glLightfv(GL_LIGHT0, GL_POSITION, light_position);
195 glEnable(GL_LIGHTING);
197 float material_ambient[] = {0.2, 0.2, 0.2, 1.0};
198 float material_diffuse[] = {0.8, 0.8, 0.8, 1.0};
199 float material_specular[] = {0.0, 0.0, 0.0, 1.0};
200 float material_emission[] = {0.0, 0.0, 0.0, 1.0};
201 glMaterialfv(GL_FRONT_AND_BACK, GL_AMBIENT, material_ambient);
202 glMaterialfv(GL_FRONT_AND_BACK, GL_DIFFUSE, material_diffuse);
203 glMaterialfv(GL_FRONT_AND_BACK, GL_SPECULAR, material_specular);
204 glMaterialfv(GL_FRONT_AND_BACK, GL_EMISSION, material_emission);
205 glDisable(GL_TEXTURE_2D);
206 glDisable(GL_COLOR_MATERIAL);
209 glTranslatef(0, 0, -elev2rad(elev));
210 glRotatef(lat, 1, 0, 0);
211 glRotatef(-lon, 0, 1, 0);
214 gdouble rg = MAX(0, 1-(elev/20000));
215 gdouble blue = MAX(0, 1-(elev/50000));
216 glClearColor(MIN(0.65,rg), MIN(0.65,rg), MIN(1,blue), 1.0f);
217 glColor4f(1, 1, 1, 1);
219 glDisable(GL_ALPHA_TEST);
221 glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
226 glEnable(GL_CULL_FACE);
230 glDepthFunc(GL_LEQUAL);
231 glEnable(GL_DEPTH_TEST);
233 glEnable(GL_LINE_SMOOTH);
235 glPolygonMode(GL_FRONT_AND_BACK, GL_FILL);
236 //glShadeModel(GL_FLAT);
238 roam_sphere_update_view(self->sphere);
245 static void on_realize(GisOpenGL *self, gpointer _)
247 g_debug("GisOpenGL: on_realize");
249 roam_sphere_update_errors(self->sphere);
251 static gboolean on_configure(GisOpenGL *self, GdkEventConfigure *event, gpointer _)
253 g_debug("GisOpenGL: on_configure");
254 _gis_opengl_begin(self);
256 double width = GTK_WIDGET(self)->allocation.width;
257 double height = GTK_WIDGET(self)->allocation.height;
259 /* Setup OpenGL Window */
260 glViewport(0, 0, width, height);
261 glMatrixMode(GL_PROJECTION);
263 double ang = atan(height/FOV_DIST);
264 gluPerspective(rad2deg(ang)*2, width/height, 1, 20*EARTH_R);
267 roam_sphere_update_errors(self->sphere);
270 _gis_opengl_end(self);
274 static void _on_expose_plugin(GisPlugin *plugin, gchar *name, GisOpenGL *self)
277 glMatrixMode(GL_PROJECTION); glPushMatrix();
278 glMatrixMode(GL_MODELVIEW); glPushMatrix();
279 gis_plugin_expose(plugin);
280 glMatrixMode(GL_PROJECTION); glPopMatrix();
281 glMatrixMode(GL_MODELVIEW); glPopMatrix();
283 static gboolean on_expose(GisOpenGL *self, GdkEventExpose *event, gpointer _)
285 g_debug("GisOpenGL: on_expose - begin");
286 _gis_opengl_begin(self);
288 glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
291 gis_plugins_foreach(GIS_VIEWER(self)->plugins,
292 G_CALLBACK(_on_expose_plugin), self);
294 if (self->wireframe) {
296 glPolygonMode(GL_FRONT_AND_BACK, GL_LINE);
297 roam_sphere_draw(self->sphere);
301 glColor4f(0.0, 0.0, 9.0, 0.6);
302 glDisable(GL_TEXTURE_2D);
303 glPolygonMode(GL_FRONT_AND_BACK, GL_LINE);
304 roam_sphere_draw(self->sphere);
305 //roam_sphere_draw_normals(self->sphere);
308 GdkGLDrawable *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
309 gdk_gl_drawable_swap_buffers(gldrawable);
311 _gis_opengl_end(self);
312 g_debug("GisOpenGL: on_expose - end\n");
316 static gboolean on_key_press(GisOpenGL *self, GdkEventKey *event, gpointer _)
318 g_debug("GisOpenGL: on_key_press - key=%x, state=%x, plus=%x",
319 event->keyval, event->state, GDK_plus);
321 guint kv = event->keyval;
325 self->wireframe = !self->wireframe;
326 gtk_widget_queue_draw(GTK_WIDGET(self));
329 else if (kv == GDK_n) roam_sphere_split_one(self->sphere);
330 else if (kv == GDK_p) roam_sphere_merge_one(self->sphere);
331 else if (kv == GDK_r) roam_sphere_split_merge(self->sphere);
332 else if (kv == GDK_u) roam_sphere_update_errors(self->sphere);
334 gtk_widget_queue_draw(GTK_WIDGET(self));
341 static gboolean _update_errors_cb(gpointer sphere)
343 roam_sphere_update_errors(sphere);
346 static void on_view_changed(GisOpenGL *self,
347 gdouble _1, gdouble _2, gdouble _3)
349 g_debug("GisOpenGL: on_view_changed");
351 _gis_opengl_begin(self);
354 g_idle_add_full(G_PRIORITY_HIGH_IDLE+30, _update_errors_cb, self->sphere, NULL);
355 //roam_sphere_update_errors(self->sphere);
357 _gis_opengl_end(self);
361 static gboolean on_idle(GisOpenGL *self)
363 //g_debug("GisOpenGL: on_idle");
365 _gis_opengl_begin(self);
366 if (roam_sphere_split_merge(self->sphere))
367 gtk_widget_queue_draw(GTK_WIDGET(self));
368 _gis_opengl_end(self);
374 /*********************
375 * GisViewer methods *
376 *********************/
377 GisViewer *gis_opengl_new(GisPlugins *plugins)
379 g_debug("GisOpenGL: new");
380 GisViewer *self = g_object_new(GIS_TYPE_OPENGL, NULL);
381 self->plugins = plugins;
385 static void gis_opengl_center_position(GisViewer *_self, gdouble lat, gdouble lon, gdouble elev)
387 GisOpenGL *self = GIS_OPENGL(_self);
388 glRotatef(lon, 0, 1, 0);
389 glRotatef(-lat, 1, 0, 0);
390 glTranslatef(0, 0, elev2rad(elev));
393 static void gis_opengl_project(GisViewer *_self,
394 gdouble lat, gdouble lon, gdouble elev,
395 gdouble *px, gdouble *py, gdouble *pz)
397 GisOpenGL *self = GIS_OPENGL(_self);
399 lle2xyz(lat, lon, elev, &x, &y, &z);
401 self->sphere->view->model,
402 self->sphere->view->proj,
403 self->sphere->view->view,
407 static void gis_opengl_render_tile(GisViewer *_self, GisTile *tile)
409 GisOpenGL *self = GIS_OPENGL(_self);
410 if (!tile || !tile->data)
412 GList *triangles = roam_sphere_get_intersect(self->sphere,
413 tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
415 g_warning("GisOpenGL: render_tiles - No triangles to draw: edges=%f,%f,%f,%f",
416 tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
417 //g_message("rendering %4d triangles for tile edges=%7.2f,%7.2f,%7.2f,%7.2f",
418 // g_list_length(triangles), tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
419 for (GList *cur = triangles; cur; cur = cur->next) {
420 RoamTriangle *tri = cur->data;
422 gdouble lat[3] = {tri->p.r->lat, tri->p.m->lat, tri->p.l->lat};
423 gdouble lon[3] = {tri->p.r->lon, tri->p.m->lon, tri->p.l->lon};
425 if (lon[0] < -90 || lon[1] < -90 || lon[2] < -90) {
426 if (lon[0] > 90) lon[0] -= 360;
427 if (lon[1] > 90) lon[1] -= 360;
428 if (lon[2] > 90) lon[2] -= 360;
431 gdouble n = tile->edge.n;
432 gdouble s = tile->edge.s;
433 gdouble e = tile->edge.e;
434 gdouble w = tile->edge.w;
436 gdouble londist = e - w;
437 gdouble latdist = n - s;
440 {(lon[0]-w)/londist, 1-(lat[0]-s)/latdist},
441 {(lon[1]-w)/londist, 1-(lat[1]-s)/latdist},
442 {(lon[2]-w)/londist, 1-(lat[2]-s)/latdist},
445 glEnable(GL_TEXTURE_2D);
446 glBindTexture(GL_TEXTURE_2D, *(guint*)tile->data);
447 glBegin(GL_TRIANGLES);
448 glNormal3dv(tri->p.r->norm); glTexCoord2dv(xy[0]); glVertex3dv((double*)tri->p.r);
449 glNormal3dv(tri->p.m->norm); glTexCoord2dv(xy[1]); glVertex3dv((double*)tri->p.m);
450 glNormal3dv(tri->p.l->norm); glTexCoord2dv(xy[2]); glVertex3dv((double*)tri->p.l);
453 g_list_free(triangles);
456 static void gis_opengl_render_tiles(GisViewer *_self, GisTile *tile)
458 GisOpenGL *self = GIS_OPENGL(_self);
459 /* Only render children if possible */
460 gboolean has_children = TRUE;
462 gis_tile_foreach(tile, child)
463 if (!child || !child->data)
464 has_children = FALSE;
466 /* Only render children */
467 gis_tile_foreach(tile, child)
468 gis_opengl_render_tiles(_self, child);
470 /* No children, render this tile */
471 gis_opengl_render_tile(_self, tile);
474 static void gis_opengl_set_height_func(GisViewer *_self, GisTile *tile,
475 RoamHeightFunc height_func, gpointer user_data, gboolean update)
477 GisOpenGL *self = GIS_OPENGL(_self);
480 /* TODO: get points? */
481 GList *triangles = roam_sphere_get_intersect(self->sphere,
482 tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
483 for (GList *cur = triangles; cur; cur = cur->next) {
484 RoamTriangle *tri = cur->data;
485 RoamPoint *points[] = {tri->p.l, tri->p.m, tri->p.r, tri->split};
486 for (int i = 0; i < G_N_ELEMENTS(points); i++) {
487 points[i]->height_func = height_func;
488 points[i]->height_data = user_data;
489 roam_point_update_height(points[i]);
492 g_list_free(triangles);
495 static void _gis_opengl_clear_height_func_rec(RoamTriangle *root)
499 RoamPoint *points[] = {root->p.l, root->p.m, root->p.r, root->split};
500 for (int i = 0; i < G_N_ELEMENTS(points); i++) {
501 points[i]->height_func = NULL;
502 points[i]->height_data = NULL;
503 roam_point_update_height(points[i]);
505 _gis_opengl_clear_height_func_rec(root->kids[0]);
506 _gis_opengl_clear_height_func_rec(root->kids[1]);
509 static void gis_opengl_clear_height_func(GisViewer *_self)
511 GisOpenGL *self = GIS_OPENGL(_self);
512 for (int i = 0; i < G_N_ELEMENTS(self->sphere->roots); i++)
513 _gis_opengl_clear_height_func_rec(self->sphere->roots[i]);
516 static void gis_opengl_begin(GisViewer *_self)
518 g_assert(GIS_IS_OPENGL(_self));
519 _gis_opengl_begin(GIS_OPENGL(_self));
522 static void gis_opengl_end(GisViewer *_self)
524 g_assert(GIS_IS_OPENGL(_self));
525 _gis_opengl_end(GIS_OPENGL(_self));
528 static void gis_opengl_add(GisViewer *_self, GisObject *object)
530 g_assert(GIS_IS_OPENGL(_self));
531 GisOpenGL *self = GIS_OPENGL(_self);
532 _load_object(self, object);
533 self->objects = g_list_prepend(self->objects, object);
536 static void gis_opengl_remove(GisViewer *_self, GisObject *object)
538 g_assert(GIS_IS_OPENGL(_self));
539 GisOpenGL *self = GIS_OPENGL(_self);
540 _free_object(self, object);
541 self->objects = g_list_remove(self->objects, object);
547 G_DEFINE_TYPE(GisOpenGL, gis_opengl, GIS_TYPE_VIEWER);
548 static void gis_opengl_init(GisOpenGL *self)
550 g_debug("GisOpenGL: init");
552 GdkGLConfig *glconfig = gdk_gl_config_new_by_mode(
553 GDK_GL_MODE_RGBA | GDK_GL_MODE_DEPTH |
554 GDK_GL_MODE_DOUBLE | GDK_GL_MODE_ALPHA);
556 g_error("Failed to create glconfig");
557 if (!gtk_widget_set_gl_capability(GTK_WIDGET(self),
558 glconfig, NULL, TRUE, GDK_GL_RGBA_TYPE))
559 g_error("GL lacks required capabilities");
560 g_object_unref(glconfig);
562 gtk_widget_set_size_request(GTK_WIDGET(self), 600, 550);
563 gtk_widget_set_events(GTK_WIDGET(self),
564 GDK_BUTTON_PRESS_MASK |
565 GDK_ENTER_NOTIFY_MASK |
567 g_object_set(self, "can-focus", TRUE, NULL);
569 self->sphere = roam_sphere_new(self);
572 self->sm_source[0] = g_timeout_add_full(G_PRIORITY_HIGH_IDLE+30, 33, (GSourceFunc)on_idle, self, NULL);
573 self->sm_source[1] = g_timeout_add_full(G_PRIORITY_HIGH_IDLE+10, 500, (GSourceFunc)on_idle, self, NULL);
576 g_signal_connect(self, "realize", G_CALLBACK(on_realize), NULL);
577 g_signal_connect(self, "configure-event", G_CALLBACK(on_configure), NULL);
578 g_signal_connect(self, "expose-event", G_CALLBACK(on_expose), NULL);
580 g_signal_connect(self, "key-press-event", G_CALLBACK(on_key_press), NULL);
582 g_signal_connect(self, "location-changed", G_CALLBACK(on_view_changed), NULL);
583 g_signal_connect(self, "rotation-changed", G_CALLBACK(on_view_changed), NULL);
585 static void gis_opengl_dispose(GObject *_self)
587 g_debug("GisOpenGL: dispose");
588 GisOpenGL *self = GIS_OPENGL(_self);
589 if (self->sm_source[0]) {
590 g_source_remove(self->sm_source[0]);
591 self->sm_source[0] = 0;
593 if (self->sm_source[1]) {
594 g_source_remove(self->sm_source[1]);
595 self->sm_source[1] = 0;
598 roam_sphere_free(self->sphere);
601 G_OBJECT_CLASS(gis_opengl_parent_class)->dispose(_self);
603 static void gis_opengl_finalize(GObject *_self)
605 g_debug("GisViewer: finalize");
606 GisOpenGL *self = GIS_OPENGL(_self);
607 G_OBJECT_CLASS(gis_opengl_parent_class)->finalize(_self);
609 static void gis_opengl_class_init(GisOpenGLClass *klass)
611 g_debug("GisOpenGL: class_init");
612 GObjectClass *gobject_class = G_OBJECT_CLASS(klass);
613 gobject_class->dispose = gis_opengl_dispose;
615 GisViewerClass *viewer_class = GIS_VIEWER_CLASS(klass);
616 viewer_class->center_position = gis_opengl_center_position;
617 viewer_class->project = gis_opengl_project;
618 viewer_class->clear_height_func = gis_opengl_clear_height_func;
619 viewer_class->set_height_func = gis_opengl_set_height_func;
620 viewer_class->render_tile = gis_opengl_render_tile;
621 viewer_class->render_tiles = gis_opengl_render_tiles;
622 viewer_class->begin = gis_opengl_begin;
623 viewer_class->end = gis_opengl_end;
624 viewer_class->add = gis_opengl_add;
625 viewer_class->remove = gis_opengl_remove;