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)
43 static void _gis_opengl_draw_triangle(GisOpenGL *self, GisTriangle *tri)
47 static void _gis_opengl_draw_quad(GisOpenGL *self, GisQuad *quad)
51 static void _gis_opengl_draw_callback(GisOpenGL *self, GisCallback *cb)
55 static void _gis_opengl_draw_marker(GisOpenGL *self, GisMarker *marker)
57 GisProjection *proj = (GisProjection*)self->sphere->view;
58 GisPoint *point = gis_object_center(marker);
59 gis_point_project(point, proj);
61 double width = GTK_WIDGET(self)->allocation.width;
62 double height = GTK_WIDGET(self)->allocation.height;
64 cairo_set_source_rgba(self->canvas, 1, 1, 1, 1);
65 cairo_arc(self->canvas, point->px, height-point->py, 4, 0, 2*G_PI);
66 cairo_fill(self->canvas);
67 cairo_move_to(self->canvas, point->px+4, height-point->py-8);
68 cairo_set_font_size(self->canvas, 10);
69 cairo_show_text(self->canvas, marker->label);
72 static void gis_opengl_draw_object(GisOpenGL *self, GisObject *object)
74 g_debug("GisOpenGL: draw_object - Drawing object of type %d", object->type);
75 switch (object->type) {
76 case GIS_TYPE_TRIANGLE: _gis_opengl_draw_triangle(self, GIS_TRIANGLE(object)); break;
77 case GIS_TYPE_QUAD : _gis_opengl_draw_quad (self, GIS_QUAD (object)); break;
78 case GIS_TYPE_CALLBACK: _gis_opengl_draw_callback(self, GIS_CALLBACK(object)); break;
79 case GIS_TYPE_MARKER : _gis_opengl_draw_marker (self, GIS_MARKER (object)); break;
80 default: g_warning("GisOpenGL: draw_object - invalid type %d", object->type);
83 static void gis_opengl_draw_objects(GisOpenGL *self)
85 g_debug("GisOpenGL: draw_objects");
87 double width = GTK_WIDGET(self)->allocation.width;
88 double height = GTK_WIDGET(self)->allocation.height;
89 cairo_surface_t *surface = cairo_get_target(self->canvas);
90 int stride = cairo_image_surface_get_stride(surface);
91 guchar *data = cairo_image_surface_get_data(surface);
92 memset(data, 0, height*stride);
95 for (GList *cur = GIS_VIEWER(self)->objects; cur; cur = cur->next)
96 gis_opengl_draw_object(self, cur->data);
98 /* Copy canvas to opengl */
99 glMatrixMode(GL_PROJECTION); glPushMatrix(); glLoadIdentity();
100 glMatrixMode(GL_MODELVIEW); glPushMatrix(); glLoadIdentity();
101 glDisable(GL_CULL_FACE);
102 glDisable(GL_DEPTH_TEST);
103 glDisable(GL_LIGHTING);
104 glDisable(GL_COLOR_MATERIAL);
105 glEnable(GL_TEXTURE_2D);
106 glBindTexture(GL_TEXTURE_2D, self->canvas_tex);
107 //glTexSubImage2D(GL_TEXTURE_2D, 0, 0,0, width,height,
108 // GL_RGBA, GL_UNSIGNED_BYTE, data);
110 glTexCoord2d(0, 0); glVertex3f(-1, 1, 1);
111 glTexCoord2d(1, 0); glVertex3f( 1, 1, 1);
112 glTexCoord2d(1, 1); glVertex3f( 1, -1, 1);
113 glTexCoord2d(0, 1); glVertex3f(-1, -1, 1);
115 glMatrixMode(GL_PROJECTION); glPopMatrix();
116 glMatrixMode(GL_MODELVIEW); glPopMatrix();
123 static void _gis_opengl_begin(GisOpenGL *self)
125 g_assert(GIS_IS_OPENGL(self));
127 GdkGLContext *glcontext = gtk_widget_get_gl_context(GTK_WIDGET(self));
128 GdkGLDrawable *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
130 if (!gdk_gl_drawable_gl_begin(gldrawable, glcontext))
131 g_assert_not_reached();
134 static void _gis_opengl_end(GisOpenGL *self)
136 g_assert(GIS_IS_OPENGL(self));
137 GdkGLDrawable *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
138 gdk_gl_drawable_gl_end(gldrawable);
144 static void _set_visuals(GisOpenGL *self)
146 glMatrixMode(GL_MODELVIEW);
150 double lat, lon, elev, rx, ry, rz;
151 gis_viewer_get_location(GIS_VIEWER(self), &lat, &lon, &elev);
152 gis_viewer_get_rotation(GIS_VIEWER(self), &rx, &ry, &rz);
153 glRotatef(rx, 1, 0, 0);
154 glRotatef(rz, 0, 0, 1);
158 float light_ambient[] = {0.7f, 0.7f, 0.7f, 1.0f};
159 float light_diffuse[] = {2.0f, 2.0f, 2.0f, 1.0f};
161 float light_ambient[] = {0.2f, 0.2f, 0.2f, 1.0f};
162 float light_diffuse[] = {5.0f, 5.0f, 5.0f, 1.0f};
164 float light_position[] = {-13*EARTH_R, 1*EARTH_R, 3*EARTH_R, 1.0f};
165 glLightfv(GL_LIGHT0, GL_AMBIENT, light_ambient);
166 glLightfv(GL_LIGHT0, GL_DIFFUSE, light_diffuse);
167 glLightfv(GL_LIGHT0, GL_POSITION, light_position);
169 glEnable(GL_LIGHTING);
171 float material_ambient[] = {0.2, 0.2, 0.2, 1.0};
172 float material_diffuse[] = {0.8, 0.8, 0.8, 1.0};
173 float material_specular[] = {0.0, 0.0, 0.0, 1.0};
174 float material_emission[] = {0.0, 0.0, 0.0, 1.0};
175 glMaterialfv(GL_FRONT_AND_BACK, GL_AMBIENT, material_ambient);
176 glMaterialfv(GL_FRONT_AND_BACK, GL_DIFFUSE, material_diffuse);
177 glMaterialfv(GL_FRONT_AND_BACK, GL_SPECULAR, material_specular);
178 glMaterialfv(GL_FRONT_AND_BACK, GL_EMISSION, material_emission);
179 glDisable(GL_TEXTURE_2D);
180 glDisable(GL_COLOR_MATERIAL);
183 glTranslatef(0, 0, -elev2rad(elev));
184 glRotatef(lat, 1, 0, 0);
185 glRotatef(-lon, 0, 1, 0);
188 gdouble rg = MAX(0, 1-(elev/20000));
189 gdouble blue = MAX(0, 1-(elev/50000));
190 glClearColor(MIN(0.65,rg), MIN(0.65,rg), MIN(1,blue), 1.0f);
191 glColor4f(1, 1, 1, 1);
193 glDisable(GL_ALPHA_TEST);
195 glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
200 glEnable(GL_CULL_FACE);
204 glDepthFunc(GL_LEQUAL);
205 glEnable(GL_DEPTH_TEST);
207 glEnable(GL_LINE_SMOOTH);
209 glPolygonMode(GL_FRONT_AND_BACK, GL_FILL);
210 //glShadeModel(GL_FLAT);
212 roam_sphere_update_view(self->sphere);
219 static void on_realize(GisOpenGL *self, gpointer _)
221 g_debug("GisOpenGL: on_realize");
223 roam_sphere_update_errors(self->sphere);
225 static gboolean on_configure(GisOpenGL *self, GdkEventConfigure *event, gpointer _)
227 g_debug("GisOpenGL: on_configure");
228 _gis_opengl_begin(self);
230 double width = GTK_WIDGET(self)->allocation.width;
231 double height = GTK_WIDGET(self)->allocation.height;
233 /* Setup OpenGL Window */
234 glViewport(0, 0, width, height);
235 glMatrixMode(GL_PROJECTION);
237 double ang = atan(height/FOV_DIST);
238 gluPerspective(rad2deg(ang)*2, width/height, 1, 20*EARTH_R);
240 /* Recreate the canvas */
242 cairo_surface_destroy(cairo_get_target(self->canvas));
243 cairo_destroy(self->canvas);
244 glDeleteTextures(1, &self->canvas_tex);
246 g_message("creating %fx%f canvas", width, height);
247 cairo_surface_t *surface = cairo_image_surface_create(CAIRO_FORMAT_ARGB32, width, height);
248 self->canvas = cairo_create(surface);
249 glEnable(GL_TEXTURE_2D);
250 glPixelStorei(GL_UNPACK_ALIGNMENT, 1);
251 glPixelStorei(GL_PACK_ALIGNMENT, 1);
252 glGenTextures(1, &self->canvas_tex);
253 glBindTexture(GL_TEXTURE_2D, self->canvas_tex);
254 glTexImage2D(GL_TEXTURE_2D, 0, GL_RGBA, width, height, 0, GL_RGBA, GL_UNSIGNED_BYTE,
255 cairo_image_surface_get_data(cairo_get_target(self->canvas)));
256 glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
257 glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
261 roam_sphere_update_errors(self->sphere);
264 _gis_opengl_end(self);
268 static void _on_expose_plugin(GisPlugin *plugin, gchar *name, GisOpenGL *self)
271 glMatrixMode(GL_PROJECTION); glPushMatrix();
272 glMatrixMode(GL_MODELVIEW); glPushMatrix();
273 gis_plugin_expose(plugin);
274 glMatrixMode(GL_PROJECTION); glPopMatrix();
275 glMatrixMode(GL_MODELVIEW); glPopMatrix();
277 static gboolean on_expose(GisOpenGL *self, GdkEventExpose *event, gpointer _)
279 g_debug("GisOpenGL: on_expose - begin");
280 _gis_opengl_begin(self);
282 glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
285 gis_plugins_foreach(GIS_VIEWER(self)->plugins,
286 G_CALLBACK(_on_expose_plugin), self);
287 gis_opengl_draw_objects(self);
288 if (self->wireframe) {
290 glPolygonMode(GL_FRONT_AND_BACK, GL_LINE);
291 roam_sphere_draw(self->sphere);
295 glColor4f(0.0, 0.0, 9.0, 0.6);
296 glDisable(GL_TEXTURE_2D);
297 glPolygonMode(GL_FRONT_AND_BACK, GL_LINE);
298 roam_sphere_draw(self->sphere);
299 //roam_sphere_draw_normals(self->sphere);
302 GdkGLDrawable *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
303 gdk_gl_drawable_swap_buffers(gldrawable);
305 _gis_opengl_end(self);
306 g_debug("GisOpenGL: on_expose - end\n");
310 static gboolean on_key_press(GisOpenGL *self, GdkEventKey *event, gpointer _)
312 g_debug("GisOpenGL: on_key_press - key=%x, state=%x, plus=%x",
313 event->keyval, event->state, GDK_plus);
315 guint kv = event->keyval;
319 self->wireframe = !self->wireframe;
320 gtk_widget_queue_draw(GTK_WIDGET(self));
323 else if (kv == GDK_n) roam_sphere_split_one(self->sphere);
324 else if (kv == GDK_p) roam_sphere_merge_one(self->sphere);
325 else if (kv == GDK_r) roam_sphere_split_merge(self->sphere);
326 else if (kv == GDK_u) roam_sphere_update_errors(self->sphere);
328 gtk_widget_queue_draw(GTK_WIDGET(self));
335 static gboolean _update_errors_cb(gpointer sphere)
337 roam_sphere_update_errors(sphere);
340 static void on_view_changed(GisOpenGL *self,
341 gdouble _1, gdouble _2, gdouble _3)
343 g_debug("GisOpenGL: on_view_changed");
345 _gis_opengl_begin(self);
348 g_idle_add_full(G_PRIORITY_HIGH_IDLE+30, _update_errors_cb, self->sphere, NULL);
349 //roam_sphere_update_errors(self->sphere);
351 _gis_opengl_end(self);
355 static gboolean on_idle(GisOpenGL *self)
357 //g_debug("GisOpenGL: on_idle");
359 _gis_opengl_begin(self);
360 if (roam_sphere_split_merge(self->sphere))
361 gtk_widget_queue_draw(GTK_WIDGET(self));
362 _gis_opengl_end(self);
368 /*********************
369 * GisViewer methods *
370 *********************/
371 GisViewer *gis_opengl_new(GisPlugins *plugins)
373 g_debug("GisOpenGL: new");
374 GisViewer *self = g_object_new(GIS_TYPE_OPENGL, NULL);
375 self->plugins = plugins;
379 static void gis_opengl_center_position(GisViewer *_self, gdouble lat, gdouble lon, gdouble elev)
381 GisOpenGL *self = GIS_OPENGL(_self);
382 glRotatef(lon, 0, 1, 0);
383 glRotatef(-lat, 1, 0, 0);
384 glTranslatef(0, 0, elev2rad(elev));
387 static void gis_opengl_project(GisViewer *_self,
388 gdouble lat, gdouble lon, gdouble elev,
389 gdouble *px, gdouble *py, gdouble *pz)
391 GisOpenGL *self = GIS_OPENGL(_self);
393 lle2xyz(lat, lon, elev, &x, &y, &z);
395 self->sphere->view->model,
396 self->sphere->view->proj,
397 self->sphere->view->view,
401 static void gis_opengl_render_tile(GisViewer *_self, GisTile *tile)
403 GisOpenGL *self = GIS_OPENGL(_self);
404 if (!tile || !tile->data)
406 GList *triangles = roam_sphere_get_intersect(self->sphere,
407 tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
409 g_warning("GisOpenGL: render_tiles - No triangles to draw: edges=%f,%f,%f,%f",
410 tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
411 //g_message("rendering %4d triangles for tile edges=%7.2f,%7.2f,%7.2f,%7.2f",
412 // g_list_length(triangles), tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
413 for (GList *cur = triangles; cur; cur = cur->next) {
414 RoamTriangle *tri = cur->data;
416 gdouble lat[3] = {tri->p.r->lat, tri->p.m->lat, tri->p.l->lat};
417 gdouble lon[3] = {tri->p.r->lon, tri->p.m->lon, tri->p.l->lon};
419 if (lon[0] < -90 || lon[1] < -90 || lon[2] < -90) {
420 if (lon[0] > 90) lon[0] -= 360;
421 if (lon[1] > 90) lon[1] -= 360;
422 if (lon[2] > 90) lon[2] -= 360;
425 gdouble n = tile->edge.n;
426 gdouble s = tile->edge.s;
427 gdouble e = tile->edge.e;
428 gdouble w = tile->edge.w;
430 gdouble londist = e - w;
431 gdouble latdist = n - s;
434 {(lon[0]-w)/londist, 1-(lat[0]-s)/latdist},
435 {(lon[1]-w)/londist, 1-(lat[1]-s)/latdist},
436 {(lon[2]-w)/londist, 1-(lat[2]-s)/latdist},
439 glEnable(GL_TEXTURE_2D);
440 glBindTexture(GL_TEXTURE_2D, *(guint*)tile->data);
441 glBegin(GL_TRIANGLES);
442 glNormal3dv(tri->p.r->norm); glTexCoord2dv(xy[0]); glVertex3dv((double*)tri->p.r);
443 glNormal3dv(tri->p.m->norm); glTexCoord2dv(xy[1]); glVertex3dv((double*)tri->p.m);
444 glNormal3dv(tri->p.l->norm); glTexCoord2dv(xy[2]); glVertex3dv((double*)tri->p.l);
447 g_list_free(triangles);
450 static void gis_opengl_render_tiles(GisViewer *_self, GisTile *tile)
452 GisOpenGL *self = GIS_OPENGL(_self);
453 /* Only render children if possible */
454 gboolean has_children = TRUE;
456 gis_tile_foreach(tile, child)
457 if (!child || !child->data)
458 has_children = FALSE;
460 /* Only render children */
461 gis_tile_foreach(tile, child)
462 gis_opengl_render_tiles(_self, child);
464 /* No children, render this tile */
465 gis_opengl_render_tile(_self, tile);
468 static void gis_opengl_set_height_func(GisViewer *_self, GisTile *tile,
469 RoamHeightFunc height_func, gpointer user_data, gboolean update)
471 GisOpenGL *self = GIS_OPENGL(_self);
474 /* TODO: get points? */
475 GList *triangles = roam_sphere_get_intersect(self->sphere,
476 tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
477 for (GList *cur = triangles; cur; cur = cur->next) {
478 RoamTriangle *tri = cur->data;
479 RoamPoint *points[] = {tri->p.l, tri->p.m, tri->p.r, tri->split};
480 for (int i = 0; i < G_N_ELEMENTS(points); i++) {
481 points[i]->height_func = height_func;
482 points[i]->height_data = user_data;
483 roam_point_update_height(points[i]);
486 g_list_free(triangles);
489 static void _gis_opengl_clear_height_func_rec(RoamTriangle *root)
493 RoamPoint *points[] = {root->p.l, root->p.m, root->p.r, root->split};
494 for (int i = 0; i < G_N_ELEMENTS(points); i++) {
495 points[i]->height_func = NULL;
496 points[i]->height_data = NULL;
497 roam_point_update_height(points[i]);
499 _gis_opengl_clear_height_func_rec(root->kids[0]);
500 _gis_opengl_clear_height_func_rec(root->kids[1]);
503 static void gis_opengl_clear_height_func(GisViewer *_self)
505 GisOpenGL *self = GIS_OPENGL(_self);
506 for (int i = 0; i < G_N_ELEMENTS(self->sphere->roots); i++)
507 _gis_opengl_clear_height_func_rec(self->sphere->roots[i]);
510 static void gis_opengl_begin(GisViewer *_self)
512 g_assert(GIS_IS_OPENGL(_self));
513 _gis_opengl_begin(GIS_OPENGL(_self));
516 static void gis_opengl_end(GisViewer *_self)
518 g_assert(GIS_IS_OPENGL(_self));
519 _gis_opengl_end(GIS_OPENGL(_self));
525 G_DEFINE_TYPE(GisOpenGL, gis_opengl, GIS_TYPE_VIEWER);
526 static void gis_opengl_init(GisOpenGL *self)
528 g_debug("GisOpenGL: init");
530 GdkGLConfig *glconfig = gdk_gl_config_new_by_mode(
531 GDK_GL_MODE_RGBA | GDK_GL_MODE_DEPTH |
532 GDK_GL_MODE_DOUBLE | GDK_GL_MODE_ALPHA);
534 g_error("Failed to create glconfig");
535 if (!gtk_widget_set_gl_capability(GTK_WIDGET(self),
536 glconfig, NULL, TRUE, GDK_GL_RGBA_TYPE))
537 g_error("GL lacks required capabilities");
538 g_object_unref(glconfig);
540 gtk_widget_set_size_request(GTK_WIDGET(self), 600, 550);
541 gtk_widget_set_events(GTK_WIDGET(self),
542 GDK_BUTTON_PRESS_MASK |
543 GDK_ENTER_NOTIFY_MASK |
545 g_object_set(self, "can-focus", TRUE, NULL);
547 self->sphere = roam_sphere_new(self);
550 self->sm_source[0] = g_timeout_add_full(G_PRIORITY_HIGH_IDLE+30, 33, (GSourceFunc)on_idle, self, NULL);
551 self->sm_source[1] = g_timeout_add_full(G_PRIORITY_HIGH_IDLE+10, 500, (GSourceFunc)on_idle, self, NULL);
554 g_signal_connect(self, "realize", G_CALLBACK(on_realize), NULL);
555 g_signal_connect(self, "configure-event", G_CALLBACK(on_configure), NULL);
556 g_signal_connect(self, "expose-event", G_CALLBACK(on_expose), NULL);
558 g_signal_connect(self, "key-press-event", G_CALLBACK(on_key_press), NULL);
560 g_signal_connect(self, "location-changed", G_CALLBACK(on_view_changed), NULL);
561 g_signal_connect(self, "rotation-changed", G_CALLBACK(on_view_changed), NULL);
563 static void gis_opengl_dispose(GObject *_self)
565 g_debug("GisOpenGL: dispose");
566 GisOpenGL *self = GIS_OPENGL(_self);
567 if (self->sm_source[0]) {
568 g_source_remove(self->sm_source[0]);
569 self->sm_source[0] = 0;
571 if (self->sm_source[1]) {
572 g_source_remove(self->sm_source[1]);
573 self->sm_source[1] = 0;
576 roam_sphere_free(self->sphere);
579 G_OBJECT_CLASS(gis_opengl_parent_class)->dispose(_self);
581 static void gis_opengl_finalize(GObject *_self)
583 g_debug("GisViewer: finalize");
584 GisOpenGL *self = GIS_OPENGL(_self);
586 cairo_surface_destroy(cairo_get_target(self->canvas));
587 cairo_destroy(self->canvas);
590 G_OBJECT_CLASS(gis_opengl_parent_class)->finalize(_self);
592 static void gis_opengl_class_init(GisOpenGLClass *klass)
594 g_debug("GisOpenGL: class_init");
595 GObjectClass *gobject_class = G_OBJECT_CLASS(klass);
596 gobject_class->dispose = gis_opengl_dispose;
598 GisViewerClass *viewer_class = GIS_VIEWER_CLASS(klass);
599 viewer_class->center_position = gis_opengl_center_position;
600 viewer_class->project = gis_opengl_project;
601 viewer_class->clear_height_func = gis_opengl_clear_height_func;
602 viewer_class->set_height_func = gis_opengl_set_height_func;
603 viewer_class->render_tile = gis_opengl_render_tile;
604 viewer_class->render_tiles = gis_opengl_render_tiles;
605 viewer_class->begin = gis_opengl_begin;
606 viewer_class->end = gis_opengl_end;