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 */
24 #include <gdk/gdkkeysyms.h>
26 #include <gtk/gtkgl.h>
30 #include "gis-opengl.h"
33 #define FOV_DIST 2000.0
34 #define MPPX(dist) (4*dist/FOV_DIST)
41 static void set_visuals(GisOpenGL *self)
43 glMatrixMode(GL_MODELVIEW);
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);
55 float light_ambient[] = {0.7f, 0.7f, 0.7f, 1.0f};
56 float light_diffuse[] = {2.0f, 2.0f, 2.0f, 1.0f};
58 float light_ambient[] = {0.2f, 0.2f, 0.2f, 1.0f};
59 float light_diffuse[] = {5.0f, 5.0f, 5.0f, 1.0f};
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);
66 glEnable(GL_LIGHTING);
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);
80 glTranslatef(0, 0, -elev2rad(elev));
81 glRotatef(lat, 1, 0, 0);
82 glRotatef(-lon, 0, 1, 0);
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);
90 glDisable(GL_ALPHA_TEST);
92 glBlendFunc(GL_SRC_ALPHA, GL_ONE_MINUS_SRC_ALPHA);
97 glEnable(GL_CULL_FACE);
101 glDepthFunc(GL_LEQUAL);
102 glEnable(GL_DEPTH_TEST);
104 glEnable(GL_LINE_SMOOTH);
106 glPolygonMode(GL_FRONT_AND_BACK, GL_FILL);
107 //glShadeModel(GL_FLAT);
114 static void on_realize(GisOpenGL *self, gpointer _)
116 g_debug("GisOpenGL: on_realize");
118 roam_sphere_update_errors(self->sphere);
120 static gboolean on_configure(GisOpenGL *self, GdkEventConfigure *event, gpointer _)
122 g_debug("GisOpenGL: on_configure");
123 gis_opengl_begin(self);
125 double width = GTK_WIDGET(self)->allocation.width;
126 double height = GTK_WIDGET(self)->allocation.height;
127 glViewport(0, 0, width, height);
129 glMatrixMode(GL_PROJECTION);
131 double ang = atan(height/FOV_DIST);
132 gluPerspective(rad2deg(ang)*2, width/height, 1, 20*EARTH_R);
135 roam_sphere_update_errors(self->sphere);
138 gis_opengl_end(self);
142 static void on_expose_plugin(GisPlugin *plugin, gchar *name, GisOpenGL *self)
145 glMatrixMode(GL_PROJECTION); glPushMatrix();
146 glMatrixMode(GL_MODELVIEW); glPushMatrix();
147 gis_plugin_expose(plugin);
148 glMatrixMode(GL_PROJECTION); glPopMatrix();
149 glMatrixMode(GL_MODELVIEW); glPopMatrix();
151 static gboolean on_expose(GisOpenGL *self, GdkEventExpose *event, gpointer _)
153 g_debug("GisOpenGL: on_expose - begin");
154 gis_opengl_begin(self);
156 glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
159 gis_plugins_foreach(self->plugins, G_CALLBACK(on_expose_plugin), self);
161 if (self->wireframe) {
163 glPolygonMode(GL_FRONT_AND_BACK, GL_LINE);
164 roam_sphere_draw(self->sphere);
168 glColor4f(0.0, 0.0, 9.0, 0.6);
169 glDisable(GL_TEXTURE_2D);
170 glPolygonMode(GL_FRONT_AND_BACK, GL_LINE);
171 roam_sphere_draw(self->sphere);
173 //roam_sphere_draw_normals(self->sphere);
177 gis_opengl_end(self);
178 gis_opengl_flush(self);
179 g_debug("GisOpenGL: on_expose - end\n");
183 static gboolean on_button_press(GisOpenGL *self, GdkEventButton *event, gpointer _)
185 g_debug("GisOpenGL: on_button_press - Grabbing focus");
186 gtk_widget_grab_focus(GTK_WIDGET(self));
190 static gboolean on_key_press(GisOpenGL *self, GdkEventKey *event, gpointer _)
192 g_debug("GisOpenGL: on_key_press - key=%x, state=%x, plus=%x",
193 event->keyval, event->state, GDK_plus);
195 double lat, lon, elev, pan;
196 gis_view_get_location(self->view, &lat, &lon, &elev);
197 pan = MIN(elev/(EARTH_R/2), 30);
198 guint kv = event->keyval;
200 if (kv == GDK_Left || kv == GDK_h) gis_view_pan(self->view, 0, -pan, 0);
201 else if (kv == GDK_Down || kv == GDK_j) gis_view_pan(self->view, -pan, 0, 0);
202 else if (kv == GDK_Up || kv == GDK_k) gis_view_pan(self->view, pan, 0, 0);
203 else if (kv == GDK_Right || kv == GDK_l) gis_view_pan(self->view, 0, pan, 0);
204 else if (kv == GDK_minus || kv == GDK_o) gis_view_zoom(self->view, 10./9);
205 else if (kv == GDK_plus || kv == GDK_i) gis_view_zoom(self->view, 9./10);
206 else if (kv == GDK_H) gis_view_rotate(self->view, 0, 0, -2);
207 else if (kv == GDK_J) gis_view_rotate(self->view, 2, 0, 0);
208 else if (kv == GDK_K) gis_view_rotate(self->view, -2, 0, 0);
209 else if (kv == GDK_L) gis_view_rotate(self->view, 0, 0, 2);
212 else if (kv == GDK_w) {self->wireframe = !self->wireframe; gtk_widget_queue_draw(GTK_WIDGET(self));}
214 else if (kv == GDK_n) roam_sphere_split_one(self->sphere);
215 else if (kv == GDK_p) roam_sphere_merge_one(self->sphere);
216 else if (kv == GDK_r) roam_sphere_split_merge(self->sphere);
217 else if (kv == GDK_u) roam_sphere_update_errors(self->sphere);
219 gtk_widget_queue_draw(GTK_WIDGET(self));
227 static gboolean _update_errors_cb(gpointer sphere)
229 roam_sphere_update_errors(sphere);
232 static void on_view_changed(GisView *view,
233 gdouble _1, gdouble _2, gdouble _3, GisOpenGL *self)
235 g_debug("GisOpenGL: on_view_changed");
237 gis_opengl_begin(self);
240 g_idle_add_full(G_PRIORITY_HIGH_IDLE+30, _update_errors_cb, self->sphere, NULL);
241 //roam_sphere_update_errors(self->sphere);
243 gis_opengl_redraw(self);
244 gis_opengl_end(self);
248 static gboolean on_idle(GisOpenGL *self)
250 //g_debug("GisOpenGL: on_idle");
252 gis_opengl_begin(self);
253 if (roam_sphere_split_merge(self->sphere))
254 gis_opengl_redraw(self);
255 gis_opengl_end(self);
264 GisOpenGL *gis_opengl_new(GisWorld *world, GisView *view, GisPlugins *plugins)
266 g_debug("GisOpenGL: new");
267 GisOpenGL *self = g_object_new(GIS_TYPE_OPENGL, NULL);
270 self->plugins = plugins;
274 g_signal_connect(self->view, "location-changed", G_CALLBACK(on_view_changed), self);
275 g_signal_connect(self->view, "rotation-changed", G_CALLBACK(on_view_changed), self);
277 self->sphere = roam_sphere_new(self);
279 return g_object_ref(self);
282 void gis_opengl_center_position(GisOpenGL *self, gdouble lat, gdouble lon, gdouble elev)
284 glRotatef(lon, 0, 1, 0);
285 glRotatef(-lat, 1, 0, 0);
286 glTranslatef(0, 0, elev2rad(elev));
289 void gis_opengl_render_tile(GisOpenGL *self, GisTile *tile)
291 if (!tile || !tile->data)
293 GList *triangles = roam_sphere_get_intersect(self->sphere,
294 tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
296 g_warning("GisOpenGL: render_tiles - No triangles to draw: edges=%f,%f,%f,%f",
297 tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
298 //g_message("rendering %4d triangles for tile edges=%7.2f,%7.2f,%7.2f,%7.2f",
299 // g_list_length(triangles), tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
300 for (GList *cur = triangles; cur; cur = cur->next) {
301 RoamTriangle *tri = cur->data;
303 gdouble lat[3] = {tri->p.r->lat, tri->p.m->lat, tri->p.l->lat};
304 gdouble lon[3] = {tri->p.r->lon, tri->p.m->lon, tri->p.l->lon};
306 if (lon[0] < -90 || lon[1] < -90 || lon[2] < -90) {
307 if (lon[0] > 90) lon[0] -= 360;
308 if (lon[1] > 90) lon[1] -= 360;
309 if (lon[2] > 90) lon[2] -= 360;
312 gdouble n = tile->edge.n;
313 gdouble s = tile->edge.s;
314 gdouble e = tile->edge.e;
315 gdouble w = tile->edge.w;
317 gdouble londist = e - w;
318 gdouble latdist = n - s;
321 {(lon[0]-w)/londist, 1-(lat[0]-s)/latdist},
322 {(lon[1]-w)/londist, 1-(lat[1]-s)/latdist},
323 {(lon[2]-w)/londist, 1-(lat[2]-s)/latdist},
326 glEnable(GL_TEXTURE_2D);
327 glBindTexture(GL_TEXTURE_2D, *(guint*)tile->data);
328 glBegin(GL_TRIANGLES);
329 glNormal3dv(tri->p.r->norm); glTexCoord2dv(xy[0]); glVertex3dv((double*)tri->p.r);
330 glNormal3dv(tri->p.m->norm); glTexCoord2dv(xy[1]); glVertex3dv((double*)tri->p.m);
331 glNormal3dv(tri->p.l->norm); glTexCoord2dv(xy[2]); glVertex3dv((double*)tri->p.l);
334 g_list_free(triangles);
337 void gis_opengl_render_tiles(GisOpenGL *opengl, GisTile *tile)
339 /* Only render children if possible */
340 gboolean has_children = TRUE;
342 gis_tile_foreach(tile, child)
343 if (!child || !child->data)
344 has_children = FALSE;
346 /* Only render children */
347 gis_tile_foreach(tile, child)
348 gis_opengl_render_tiles(opengl, child);
350 /* No children, render this tile */
351 gis_opengl_render_tile(opengl, tile);
354 void gis_opengl_set_height_func(GisOpenGL *self, GisTile *tile,
355 RoamHeightFunc height_func, gpointer user_data, gboolean update)
359 /* TODO: get points? */
360 GList *triangles = roam_sphere_get_intersect(self->sphere,
361 tile->edge.n, tile->edge.s, tile->edge.e, tile->edge.w);
362 for (GList *cur = triangles; cur; cur = cur->next) {
363 RoamTriangle *tri = cur->data;
364 RoamPoint *points[] = {tri->p.l, tri->p.m, tri->p.r, tri->split};
365 for (int i = 0; i < G_N_ELEMENTS(points); i++) {
366 points[i]->height_func = height_func;
367 points[i]->height_data = user_data;
368 roam_point_update_height(points[i]);
371 g_list_free(triangles);
374 static void _gis_opengl_clear_height_func_rec(RoamTriangle *root)
378 RoamPoint *points[] = {root->p.l, root->p.m, root->p.r, root->split};
379 for (int i = 0; i < G_N_ELEMENTS(points); i++) {
380 points[i]->height_func = NULL;
381 points[i]->height_data = NULL;
382 roam_point_update_height(points[i]);
384 _gis_opengl_clear_height_func_rec(root->kids[0]);
385 _gis_opengl_clear_height_func_rec(root->kids[1]);
387 void gis_opengl_clear_height_func(GisOpenGL *self)
389 for (int i = 0; i < G_N_ELEMENTS(self->sphere->roots); i++)
390 _gis_opengl_clear_height_func_rec(self->sphere->roots[i]);
393 void gis_opengl_redraw(GisOpenGL *self)
395 g_debug("GisOpenGL: redraw");
396 gtk_widget_queue_draw(GTK_WIDGET(self));
398 void gis_opengl_begin(GisOpenGL *self)
400 g_assert(GIS_IS_OPENGL(self));
402 GdkGLContext *glcontext = gtk_widget_get_gl_context(GTK_WIDGET(self));
403 GdkGLDrawable *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
405 if (!gdk_gl_drawable_gl_begin(gldrawable, glcontext))
406 g_assert_not_reached();
408 void gis_opengl_end(GisOpenGL *self)
410 g_assert(GIS_IS_OPENGL(self));
411 GdkGLDrawable *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
412 gdk_gl_drawable_gl_end(gldrawable);
415 void gis_opengl_flush(GisOpenGL *self)
417 g_assert(GIS_IS_OPENGL(self));
418 GdkGLDrawable *gldrawable = gtk_widget_get_gl_drawable(GTK_WIDGET(self));
419 if (gdk_gl_drawable_is_double_buffered(gldrawable))
420 gdk_gl_drawable_swap_buffers(gldrawable);
423 gdk_gl_drawable_gl_end(gldrawable);
430 G_DEFINE_TYPE(GisOpenGL, gis_opengl, GTK_TYPE_DRAWING_AREA);
431 static void gis_opengl_init(GisOpenGL *self)
433 g_debug("GisOpenGL: init");
435 GdkGLConfig *glconfig = gdk_gl_config_new_by_mode(
436 GDK_GL_MODE_RGBA | GDK_GL_MODE_DEPTH |
437 GDK_GL_MODE_DOUBLE | GDK_GL_MODE_ALPHA);
439 g_error("Failed to create glconfig");
440 if (!gtk_widget_set_gl_capability(GTK_WIDGET(self),
441 glconfig, NULL, TRUE, GDK_GL_RGBA_TYPE))
442 g_error("GL lacks required capabilities");
443 g_object_unref(glconfig);
445 gtk_widget_set_size_request(GTK_WIDGET(self), 600, 550);
446 gtk_widget_set_events(GTK_WIDGET(self),
447 GDK_BUTTON_PRESS_MASK |
448 GDK_ENTER_NOTIFY_MASK |
450 g_object_set(self, "can-focus", TRUE, NULL);
453 self->sm_source[0] = g_timeout_add_full(G_PRIORITY_HIGH_IDLE+30, 33, (GSourceFunc)on_idle, self, NULL);
454 self->sm_source[1] = g_timeout_add_full(G_PRIORITY_HIGH_IDLE+10, 500, (GSourceFunc)on_idle, self, NULL);
457 g_signal_connect(self, "realize", G_CALLBACK(on_realize), NULL);
458 g_signal_connect(self, "configure-event", G_CALLBACK(on_configure), NULL);
459 g_signal_connect(self, "expose-event", G_CALLBACK(on_expose), NULL);
461 g_signal_connect(self, "button-press-event", G_CALLBACK(on_button_press), NULL);
462 g_signal_connect(self, "enter-notify-event", G_CALLBACK(on_button_press), NULL);
463 g_signal_connect(self, "key-press-event", G_CALLBACK(on_key_press), NULL);
465 static GObject *gis_opengl_constructor(GType gtype, guint n_properties,
466 GObjectConstructParam *properties)
468 g_debug("GisOpengl: constructor");
469 GObjectClass *parent_class = G_OBJECT_CLASS(gis_opengl_parent_class);
470 return parent_class->constructor(gtype, n_properties, properties);
472 static void gis_opengl_dispose(GObject *_self)
474 g_debug("GisOpenGL: dispose");
475 GisOpenGL *self = GIS_OPENGL(_self);
476 if (self->sm_source[0]) {
477 g_source_remove(self->sm_source[0]);
478 self->sm_source[0] = 0;
480 if (self->sm_source[1]) {
481 g_source_remove(self->sm_source[1]);
482 self->sm_source[1] = 0;
485 roam_sphere_free(self->sphere);
489 g_object_unref(self->world);
493 g_object_unref(self->view);
496 G_OBJECT_CLASS(gis_opengl_parent_class)->dispose(_self);
498 static void gis_opengl_finalize(GObject *_self)
500 g_debug("GisOpenGL: finalize");
501 GisOpenGL *self = GIS_OPENGL(_self);
502 G_OBJECT_CLASS(gis_opengl_parent_class)->finalize(_self);
504 static void gis_opengl_class_init(GisOpenGLClass *klass)
506 g_debug("GisOpenGL: class_init");
507 GObjectClass *gobject_class = G_OBJECT_CLASS(klass);
508 gobject_class->constructor = gis_opengl_constructor;
509 gobject_class->dispose = gis_opengl_dispose;
510 gobject_class->finalize = gis_opengl_finalize;