+/**************
+ * RadarConus *
+ **************/
+#define CONUS_NORTH 50.406626367301044
+#define CONUS_WEST -127.620375523875420
+#define CONUS_WIDTH 3400.0
+#define CONUS_HEIGHT 1600.0
+#define CONUS_DEG_PER_PX 0.017971305190311
+
+struct _RadarConus {
+ GisViewer *viewer;
+ GisHttp *http;
+ GtkWidget *config;
+ time_t time;
+ const gchar *message;
+ GStaticMutex loading;
+
+ gchar *path;
+ GisTile *tile[2];
+ gpointer *tile_ref[2];
+
+ guint time_id; // "time-changed" callback ID
+ guint refresh_id; // "refresh" callback ID
+};
+
+void _conus_update_loading(gchar *file, goffset cur,
+ goffset total, gpointer _conus)
+{
+ RadarConus *conus = _conus;
+ GtkWidget *progress_bar = gtk_bin_get_child(GTK_BIN(conus->config));
+ double percent = (double)cur/total;
+ gtk_progress_bar_set_fraction(GTK_PROGRESS_BAR(progress_bar), MIN(percent, 1.0));
+ gchar *msg = g_strdup_printf("Loading... %5.1f%% (%.2f/%.2f MB)",
+ percent*100, (double)cur/1000000, (double)total/1000000);
+ gtk_progress_bar_set_text(GTK_PROGRESS_BAR(progress_bar), msg);
+ g_free(msg);
+}
+
+/* Copy images to graphics memory */
+static void _conus_update_end_copy(GisTile *tile, guchar *pixels)
+{
+ if (!tile->data) {
+ tile->data = g_new0(guint, 1);
+ glGenTextures(1, tile->data);
+ }
+
+ guint *tex = tile->data;
+ glBindTexture(GL_TEXTURE_2D, *tex);
+
+ glPixelStorei(GL_UNPACK_ALIGNMENT, 1);
+ glPixelStorei(GL_PACK_ALIGNMENT, 1);
+ glTexImage2D(GL_TEXTURE_2D, 0, 4, 2048, 2048, 0,
+ GL_RGBA, GL_UNSIGNED_BYTE, NULL);
+ glTexSubImage2D(GL_TEXTURE_2D, 0, 1,1, CONUS_WIDTH/2,CONUS_HEIGHT,
+ GL_RGBA, GL_UNSIGNED_BYTE, pixels);
+ tile->coords.n = 1.0/(CONUS_WIDTH/2);
+ tile->coords.w = 1.0/ CONUS_HEIGHT;
+ tile->coords.s = tile->coords.n + CONUS_HEIGHT / 2048.0;
+ tile->coords.e = tile->coords.w + (CONUS_WIDTH/2) / 2048.0;
+ glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
+ glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
+ glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_CLAMP);
+ glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_CLAMP);
+ glFlush();
+}
+
+/* Split the pixbuf into east and west halves (with 2K sides)
+ * Also map the pixbuf's alpha values */
+static void _conus_update_end_split(guchar *pixels, guchar *west, guchar *east,
+ gint width, gint height, gint pxsize)
+{
+ g_debug("GisPluginRadar: _conus_update_end_split");
+ guchar *out[] = {west,east};
+ const guchar alphamap[][4] = {
+ {0x04, 0xe9, 0xe7, 0x30},
+ {0x01, 0x9f, 0xf4, 0x60},
+ {0x03, 0x00, 0xf4, 0x90},
+ };
+ for (int y = 0; y < height; y++)
+ for (int x = 0; x < width; x++) {
+ gint subx = x % (width/2);
+ gint idx = x / (width/2);
+ guchar *src = &pixels[(y*width+x)*pxsize];
+ guchar *dst = &out[idx][(y*(width/2)+subx)*4];
+ if (src[0] > 0xe0 &&
+ src[1] > 0xe0 &&
+ src[2] > 0xe0) {
+ dst[3] = 0x00;
+ } else {
+ dst[0] = src[0];
+ dst[1] = src[1];
+ dst[2] = src[2];
+ dst[3] = 0xff;
+ for (int j = 0; j < G_N_ELEMENTS(alphamap); j++)
+ if (src[0] == alphamap[j][0] &&
+ src[1] == alphamap[j][1] &&
+ src[2] == alphamap[j][2])
+ dst[3] = alphamap[j][3];
+ }
+ }
+}
+
+gboolean _conus_update_end(gpointer _conus)
+{
+ RadarConus *conus = _conus;
+ g_debug("GisPluginRadar: _conus_update_end");
+
+ /* Check error status */
+ if (conus->message) {
+ g_warning("GisPluginRadar: _conus_update_end - %s", conus->message);
+ _gtk_bin_set_child(GTK_BIN(conus->config), gtk_label_new(conus->message));
+ goto out;
+ }
+
+ /* Load and pixbuf */
+ GError *error = NULL;
+ GdkPixbuf *pixbuf = gdk_pixbuf_new_from_file(conus->path, &error);
+ guchar *pixels = gdk_pixbuf_get_pixels(pixbuf);
+ gint width = gdk_pixbuf_get_width(pixbuf);
+ gint height = gdk_pixbuf_get_height(pixbuf);
+ gint pxsize = gdk_pixbuf_get_has_alpha(pixbuf) ? 4 : 3;
+
+ /* Split pixels into east/west parts */
+ guchar *pixels_west = g_malloc(4*(width/2)*height);
+ guchar *pixels_east = g_malloc(4*(width/2)*height);
+ _conus_update_end_split(pixels, pixels_west, pixels_east,
+ width, height, pxsize);
+ g_object_unref(pixbuf);
+
+ /* Copy pixels to graphics memory */
+ _conus_update_end_copy(conus->tile[0], pixels_west);
+ _conus_update_end_copy(conus->tile[1], pixels_east);
+ g_free(pixels_west);
+ g_free(pixels_east);
+
+ /* Update GUI */
+ gchar *label = g_path_get_basename(conus->path);
+ _gtk_bin_set_child(GTK_BIN(conus->config), gtk_label_new(label));
+ gtk_widget_queue_draw(GTK_WIDGET(conus->viewer));
+ g_free(label);
+
+out:
+ g_free(conus->path);
+ g_static_mutex_unlock(&conus->loading);
+ return FALSE;
+}
+
+gpointer _conus_update_thread(gpointer _conus)
+{
+ RadarConus *conus = _conus;
+ conus->message = NULL;
+
+ /* Find nearest */
+ g_debug("GisPluginRadar: _conus_update_thread - nearest");
+ gboolean offline = gis_viewer_get_offline(conus->viewer);
+ gchar *conus_url = "http://radar.weather.gov/Conus/RadarImg/";
+ GList *files = gis_http_available(conus->http,
+ "^Conus_[^\"]*_N0Ronly.gif$", "",
+ NULL, (offline ? NULL : conus_url));
+ gchar *nearest = _find_nearest(conus->time, files, 6, "%Y%m%d_%H%M");
+ g_list_foreach(files, (GFunc)g_free, NULL);
+ g_list_free(files);
+ if (!nearest) {
+ conus->message = "No suitable files";
+ goto out;
+ }
+
+ /* Fetch the image */
+ g_debug("GisPluginRadar: _conus_update_thread - fetch");
+ gchar *uri = g_strconcat(conus_url, nearest, NULL);
+ conus->path = gis_http_fetch(conus->http, uri, nearest, GIS_ONCE,
+ _conus_update_loading, conus);
+ g_free(nearest);
+ g_free(uri);
+ if (!conus->path) {
+ conus->message = "Fetch failed";
+ goto out;
+ }
+
+out:
+ g_debug("GisPluginRadar: _conus_update_thread - done");
+ g_idle_add(_conus_update_end, conus);
+ return NULL;
+}
+
+void _conus_update(RadarConus *conus)
+{
+ if (!g_static_mutex_trylock(&conus->loading))
+ return;
+ conus->time = gis_viewer_get_time(conus->viewer);
+ g_debug("GisPluginRadar: _conus_update - %d",
+ (gint)conus->time);
+
+ /* Add a progress bar */
+ GtkWidget *progress = gtk_progress_bar_new();
+ gtk_progress_bar_set_text(GTK_PROGRESS_BAR(progress), "Loading...");
+ _gtk_bin_set_child(GTK_BIN(conus->config), progress);
+
+ g_thread_create(_conus_update_thread, conus, FALSE, NULL);
+}
+
+void radar_conus_toggle(RadarConus *conus)
+{
+ GIS_OBJECT(conus->tile[0])->hidden =
+ !GIS_OBJECT(conus->tile[0])->hidden;
+ GIS_OBJECT(conus->tile[1])->hidden =
+ !GIS_OBJECT(conus->tile[1])->hidden;
+ gtk_widget_queue_draw(GTK_WIDGET(conus->viewer));
+}
+
+RadarConus *radar_conus_new(GtkWidget *pconfig,
+ GisViewer *viewer, GisHttp *http)
+{
+ RadarConus *conus = g_new0(RadarConus, 1);
+ conus->viewer = g_object_ref(viewer);
+ conus->http = http;
+ conus->config = gtk_alignment_new(0, 0, 1, 1);
+ g_static_mutex_init(&conus->loading);
+
+ gdouble south = CONUS_NORTH - CONUS_DEG_PER_PX*CONUS_HEIGHT;
+ gdouble east = CONUS_WEST + CONUS_DEG_PER_PX*CONUS_WIDTH;
+ gdouble mid = CONUS_WEST + CONUS_DEG_PER_PX*CONUS_WIDTH/2;
+ conus->tile[0] = gis_tile_new(NULL, CONUS_NORTH, south, mid, CONUS_WEST);
+ conus->tile[1] = gis_tile_new(NULL, CONUS_NORTH, south, east, mid);
+ conus->tile[0]->zindex = 2;
+ conus->tile[1]->zindex = 1;
+ conus->tile_ref[0] = gis_viewer_add(viewer,
+ GIS_OBJECT(conus->tile[0]), GIS_LEVEL_WORLD, TRUE);
+ conus->tile_ref[1] = gis_viewer_add(viewer,
+ GIS_OBJECT(conus->tile[1]), GIS_LEVEL_WORLD, TRUE);
+
+ conus->time_id = g_signal_connect_swapped(viewer, "time-changed",
+ G_CALLBACK(_conus_update), conus);
+ conus->refresh_id = g_signal_connect_swapped(viewer, "refresh",
+ G_CALLBACK(_conus_update), conus);
+
+ gtk_notebook_append_page(GTK_NOTEBOOK(pconfig), conus->config,
+ _gtk_check_label_new("Conus", TRUE,
+ G_CALLBACK(radar_conus_toggle), conus));
+ _conus_update(conus);
+ return conus;
+}
+
+void radar_conus_free(RadarConus *conus)
+{
+ g_signal_handler_disconnect(conus->viewer, conus->time_id);
+ g_signal_handler_disconnect(conus->viewer, conus->refresh_id);
+
+ for (int i = 0; i < 2; i++) {
+ GisTile *tile = conus->tile[i];
+ gpointer ref = conus->tile_ref[i];
+ if (tile->data) {
+ glDeleteTextures(1, tile->data);
+ g_free(tile->data);
+ }
+ gis_viewer_remove(conus->viewer, ref);
+ }
+
+ g_object_unref(conus->viewer);
+ g_free(conus);
+}
+
+