+/**************
+ * 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 {
+ GritsViewer *viewer;
+ GritsHttp *http;
+ GtkWidget *config;
+ time_t time;
+ const gchar *message;
+ GMutex loading;
+
+ gchar *path;
+ GritsTile *tile[2];
+
+ guint time_id; // "time-changed" callback ID
+ guint refresh_id; // "refresh" callback ID
+ guint idle_source; // _conus_update_end idle source
+};
+
+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(GritsTile *tile, guchar *pixels)
+{
+ if (!tile->tex)
+ glGenTextures(1, &tile->tex);
+
+ gchar *clear = g_malloc0(2048*2048*4);
+ glBindTexture(GL_TEXTURE_2D, tile->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, clear);
+ 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_LINEAR);
+ glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR);
+ glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_CLAMP_TO_EDGE);
+ glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_CLAMP_TO_EDGE);
+ glFlush();
+ g_free(clear);
+}
+
+/* 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("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 * 0.75;
+ 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("Conus: update_end");
+
+ /* Check error status */
+ if (conus->message) {
+ g_warning("Conus: update_end - %s", conus->message);
+ aweather_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);
+ if (!pixbuf || error) {
+ g_warning("Conus: update_end - error loading pixbuf: %s", conus->path);
+ aweather_bin_set_child(GTK_BIN(conus->config), gtk_label_new("Error loading pixbuf"));
+ g_remove(conus->path);
+ goto out;
+ }
+
+ /* Split pixels into east/west parts */
+ 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;
+ 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);
+ aweather_bin_set_child(GTK_BIN(conus->config), gtk_label_new(label));
+ grits_viewer_queue_draw(conus->viewer);
+ g_free(label);
+
+out:
+ conus->idle_source = 0;
+ g_free(conus->path);
+ g_mutex_unlock(&conus->loading);
+ return FALSE;
+}
+
+gpointer _conus_update_thread(gpointer _conus)
+{
+ RadarConus *conus = _conus;
+ conus->message = NULL;
+
+ /* Find nearest */
+ g_debug("Conus: update_thread - nearest");
+ gboolean offline = grits_viewer_get_offline(conus->viewer);
+ gchar *conus_url = "http://radar.weather.gov/Conus/RadarImg/";
+ gchar *nearest;
+ if (time(NULL) - conus->time < 60*60*5 && !offline) {
+ /* radar.weather.gov is full of lies.
+ * the index pages get cached and out of date */
+ /* gmtime is not thread safe, but it's not used very often so
+ * hopefully it'll be alright for now... :-( */
+ struct tm *tm = gmtime(&conus->time);
+ time_t onthe8 = conus->time - 60*((tm->tm_min+1)%10+1);
+ tm = gmtime(&onthe8);
+ nearest = g_strdup_printf("Conus_%04d%02d%02d_%02d%02d_N0Ronly.gif",
+ tm->tm_year+1900, tm->tm_mon+1, tm->tm_mday,
+ tm->tm_hour, tm->tm_min);
+ } else {
+ GList *files = grits_http_available(conus->http,
+ "^Conus_[^\"]*_N0Ronly.gif$", "", NULL, NULL);
+ nearest = _find_nearest(conus->time, files, 6);
+ 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("Conus: update_thread - fetch");
+ gchar *uri = g_strconcat(conus_url, nearest, NULL);
+ conus->path = grits_http_fetch(conus->http, uri, nearest,
+ offline ? GRITS_LOCAL : GRITS_ONCE,
+ _conus_update_loading, conus);
+ g_free(nearest);
+ g_free(uri);
+ if (!conus->path) {
+ conus->message = "Fetch failed";
+ goto out;
+ }
+
+out:
+ g_debug("Conus: update_thread - done");
+ if (!conus->idle_source)
+ conus->idle_source = g_idle_add(_conus_update_end, conus);
+ return NULL;
+}
+
+void _conus_update(RadarConus *conus)