Add z-index to tiles
[grits] / src / plugins / sat.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 /**
19  * SECTION:sat
20  * @short_description: Satellite plugin
21  *
22  * #GisPluginSat provides overlays using satellite imagery. This is mostly
23  * provided by NASA's Blue Marble Next Generation.
24  */
25
26 #include <time.h>
27 #include <string.h>
28 #include <glib/gstdio.h>
29 #include <GL/gl.h>
30
31 #include <gis.h>
32
33 #include "sat.h"
34
35 #define MAX_RESOLUTION 500
36 #define TILE_WIDTH     1024
37 #define TILE_HEIGHT    512
38
39 struct _LoadTileData {
40         GisPluginSat *sat;
41         GisTile      *tile;
42         GdkPixbuf    *pixbuf;
43 };
44 static gboolean _load_tile_cb(gpointer _data)
45 {
46         struct _LoadTileData *data = _data;
47         GisPluginSat *sat    = data->sat;
48         GisTile      *tile   = data->tile;
49         GdkPixbuf    *pixbuf = data->pixbuf;
50         g_free(data);
51
52         /* Create Texture */
53         g_debug("GisPluginSat: _load_tile_cb start");
54         guchar   *pixels = gdk_pixbuf_get_pixels(pixbuf);
55         gboolean  alpha  = gdk_pixbuf_get_has_alpha(pixbuf);
56         gint      width  = gdk_pixbuf_get_width(pixbuf);
57         gint      height = gdk_pixbuf_get_height(pixbuf);
58
59         /* Draw a border */
60         //gint border = 10;
61         //gint stride = gdk_pixbuf_get_rowstride(pixbuf);
62         //for (int i = 0; i < border; i++) {
63         //      memset(&pixels[(       i)*stride], 0xff, stride);
64         //      memset(&pixels[(height-i)*stride], 0xff, stride);
65         //}
66         //for (int i = 0; i < height; i++) {
67         //      memset(&pixels[(i*stride)], 0xff, border*4);
68         //      memset(&pixels[(i*stride)+((width-border)*4)], 0xff, border*4);
69         //}
70
71         guint *tex = g_new0(guint, 1);
72         glGenTextures(1, tex);
73         glBindTexture(GL_TEXTURE_2D, *tex);
74
75         glPixelStorei(GL_UNPACK_ALIGNMENT, 1);
76         glPixelStorei(GL_PACK_ALIGNMENT, 1);
77         glTexImage2D(GL_TEXTURE_2D, 0, 4, width, height, 0,
78                         (alpha ? GL_RGBA : GL_RGB), GL_UNSIGNED_BYTE, pixels);
79         glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_LINEAR);
80         glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR);
81         glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_CLAMP);
82         glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_CLAMP);
83         glFlush();
84
85         tile->data = tex;
86         gtk_widget_queue_draw(GTK_WIDGET(sat->viewer));
87         g_object_unref(pixbuf);
88         return FALSE;
89 }
90
91 static void _load_tile(GisTile *tile, gpointer _sat)
92 {
93         GisPluginSat *sat = _sat;
94         g_debug("GisPluginSat: _load_tile start %p", g_thread_self());
95         char *path = gis_wms_fetch(sat->wms, tile, GIS_ONCE, NULL, NULL);
96         struct _LoadTileData *data = g_new0(struct _LoadTileData, 1);
97         data->sat   = sat;
98         data->tile   = tile;
99         data->pixbuf = gdk_pixbuf_new_from_file(path, NULL);
100         if (data->pixbuf) {
101                 g_idle_add_full(G_PRIORITY_LOW, _load_tile_cb, data, NULL);
102         } else {
103                 g_warning("GisPluginSat: _load_tile - Error loading pixbuf %s", path);
104                 g_free(data);
105                 g_remove(path);
106         }
107         g_free(path);
108         g_debug("GisPluginSat: _load_tile end %p", g_thread_self());
109 }
110
111 static gboolean _free_tile_cb(gpointer data)
112 {
113         glDeleteTextures(1, data);
114         g_free(data);
115         return FALSE;
116 }
117 static void _free_tile(GisTile *tile, gpointer _sat)
118 {
119         GisPluginSat *sat = _sat;
120         g_debug("GisPluginSat: _free_tile: %p=%d", tile->data, *(guint*)tile->data);
121         g_idle_add_full(G_PRIORITY_LOW, _free_tile_cb, tile->data, NULL);
122 }
123
124 static gpointer _update_tiles(gpointer _sat)
125 {
126         g_debug("GisPluginSat: _update_tiles");
127         GisPluginSat *sat = _sat;
128         if (!g_mutex_trylock(sat->mutex))
129                 return NULL;
130         GisPoint eye;
131         gis_viewer_get_location(sat->viewer, &eye.lat, &eye.lon, &eye.elev);
132         gis_tile_update(sat->tiles, &eye,
133                         MAX_RESOLUTION, TILE_WIDTH, TILE_WIDTH,
134                         _load_tile, sat);
135         gis_tile_gc(sat->tiles, time(NULL)-10,
136                         _free_tile, sat);
137         g_mutex_unlock(sat->mutex);
138         return NULL;
139 }
140
141 /*************
142  * Callbacks *
143  *************/
144 static void _on_location_changed(GisViewer *viewer,
145                 gdouble lat, gdouble lon, gdouble elev, GisPluginSat *sat)
146 {
147         g_thread_create(_update_tiles, sat, FALSE, NULL);
148 }
149
150 static gpointer _threaded_init(GisPluginSat *sat)
151 {
152         _load_tile(sat->tiles, sat);
153         _update_tiles(sat);
154         return NULL;
155 }
156
157 /***********
158  * Methods *
159  ***********/
160 /**
161  * gis_plugin_sat_new:
162  * @viewer: the #GisViewer to use for drawing
163  *
164  * Create a new instance of the satellite plugin.
165  *
166  * Returns: the new #GisPluginSat
167  */
168 GisPluginSat *gis_plugin_sat_new(GisViewer *viewer)
169 {
170         g_debug("GisPluginSat: new");
171         GisPluginSat *sat = g_object_new(GIS_TYPE_PLUGIN_SAT, NULL);
172         sat->viewer = g_object_ref(viewer);
173
174         /* Load initial tiles */
175         g_thread_create((GThreadFunc)_threaded_init, sat, FALSE, NULL);
176
177         /* Connect signals */
178         sat->sigid = g_signal_connect(sat->viewer, "location-changed",
179                         G_CALLBACK(_on_location_changed), sat);
180
181         /* Add renderers */
182         gis_viewer_add(viewer, GIS_OBJECT(sat->tiles), GIS_LEVEL_WORLD, FALSE);
183
184         return sat;
185 }
186
187
188 /****************
189  * GObject code *
190  ****************/
191 /* Plugin init */
192 static void gis_plugin_sat_plugin_init(GisPluginInterface *iface);
193 G_DEFINE_TYPE_WITH_CODE(GisPluginSat, gis_plugin_sat, G_TYPE_OBJECT,
194                 G_IMPLEMENT_INTERFACE(GIS_TYPE_PLUGIN,
195                         gis_plugin_sat_plugin_init));
196 static void gis_plugin_sat_plugin_init(GisPluginInterface *iface)
197 {
198         g_debug("GisPluginSat: plugin_init");
199         /* Add methods to the interface */
200 }
201 /* Class/Object init */
202 static void gis_plugin_sat_init(GisPluginSat *sat)
203 {
204         g_debug("GisPluginSat: init");
205         /* Set defaults */
206         sat->mutex = g_mutex_new();
207         sat->tiles = gis_tile_new(NULL, NORTH, SOUTH, EAST, WEST);
208         sat->wms   = gis_wms_new(
209                 "http://www.nasa.network.com/wms", "bmng200406", "image/jpeg",
210                 "bmng/", "jpg", TILE_WIDTH, TILE_HEIGHT);
211 }
212 static void gis_plugin_sat_dispose(GObject *gobject)
213 {
214         g_debug("GisPluginSat: dispose");
215         GisPluginSat *sat = GIS_PLUGIN_SAT(gobject);
216         /* Drop references */
217         if (sat->viewer) {
218                 g_signal_handler_disconnect(sat->viewer, sat->sigid);
219                 g_object_unref(sat->viewer);
220                 sat->viewer = NULL;
221         }
222         G_OBJECT_CLASS(gis_plugin_sat_parent_class)->dispose(gobject);
223 }
224 static void gis_plugin_sat_finalize(GObject *gobject)
225 {
226         g_debug("GisPluginSat: finalize");
227         GisPluginSat *sat = GIS_PLUGIN_SAT(gobject);
228         /* Free data */
229         gis_tile_free(sat->tiles, _free_tile, sat);
230         gis_wms_free(sat->wms);
231         g_mutex_free(sat->mutex);
232         G_OBJECT_CLASS(gis_plugin_sat_parent_class)->finalize(gobject);
233
234 }
235 static void gis_plugin_sat_class_init(GisPluginSatClass *klass)
236 {
237         g_debug("GisPluginSat: class_init");
238         GObjectClass *gobject_class = (GObjectClass*)klass;
239         gobject_class->dispose  = gis_plugin_sat_dispose;
240         gobject_class->finalize = gis_plugin_sat_finalize;
241 }