Pileus Git - aweather/blob - src/plugins/radar.c

   1 /*
   2  * Copyright (C) 2009-2012 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 #define _XOPEN_SOURCE
  19 #include <time.h>
  20 #include <config.h>
  21 #include <glib/gstdio.h>
  22 #include <gtk/gtk.h>
  23 #include <gio/gio.h>
  24 #include <math.h>
  25 #include <rsl.h>
  26
  27 #include <grits.h>
  28
  29 #include "radar.h"
  30 #include "level2.h"
  31 #include "../aweather-location.h"
  32
  33 #include "../compat.h"
  34
  35 static void aweather_bin_set_child(GtkBin *bin, GtkWidget *new)
  36 {
  37         GtkWidget *old = gtk_bin_get_child(bin);
  38         if (old)
  39                 gtk_widget_destroy(old);
  40         gtk_container_add(GTK_CONTAINER(bin), new);
  41         gtk_widget_show_all(new);
  42 }
  43
  44 static gchar *_find_nearest(time_t time, GList *files,
  45                 gsize offset)
  46 {
  47         g_debug("RadarSite: find_nearest ...");
  48         time_t  nearest_time = 0;
  49         char   *nearest_file = NULL;
  50
  51         struct tm tm = {};
  52         for (GList *cur = files; cur; cur = cur->next) {
  53                 gchar *file = cur->data;
  54                 sscanf(file+offset, "%4d%2d%2d_%2d%2d",
  55                                 &tm.tm_year, &tm.tm_mon, &tm.tm_mday,
  56                                 &tm.tm_hour, &tm.tm_min);
  57                 tm.tm_year -= 1900;
  58                 tm.tm_mon  -= 1;
  59                 if (ABS(time - mktime(&tm)) <
  60                     ABS(time - nearest_time)) {
  61                         nearest_file = file;
  62                         nearest_time = mktime(&tm);
  63                 }
  64         }
  65
  66         g_debug("RadarSite: find_nearest = %s", nearest_file);
  67         if (nearest_file)
  68                 return g_strdup(nearest_file);
  69         else
  70                 return NULL;
  71 }
  72
  73
  74 /**************
  75  * RadarSites *
  76  **************/
  77 typedef enum {
  78         STATUS_UNLOADED,
  79         STATUS_LOADING,
  80         STATUS_LOADED,
  81 } RadarSiteStatus;
  82 struct _RadarSite {
  83         /* Information */
  84         city_t         *city;
  85         GritsMarker    *marker;      // Map marker for grits
  86
  87         /* Stuff from the parents */
  88         GritsViewer    *viewer;
  89         GritsHttp      *http;
  90         GritsPrefs     *prefs;
  91         GtkWidget      *pconfig;
  92
  93         /* When loaded */
  94         gboolean        hidden;
  95         RadarSiteStatus status;      // Loading status for the site
  96         GtkWidget      *config;
  97         AWeatherLevel2 *level2;      // The Level2 structure for the current volume
  98
  99         /* Internal data */
 100         time_t          time;        // Current timestamp of the level2
 101         gchar          *message;     // Error message set while updating
 102         guint           time_id;     // "time-changed"     callback ID
 103         guint           refresh_id;  // "refresh"          callback ID
 104         guint           location_id; // "locaiton-changed" callback ID
 105         guint           idle_source; // _site_update_end idle source
 106 };
 107
 108 /* format: http://mesonet.agron.iastate.edu/data/nexrd2/raw/KABR/KABR_20090510_0323 */
 109 void _site_update_loading(gchar *file, goffset cur,
 110                 goffset total, gpointer _site)
 111 {
 112         RadarSite *site = _site;
 113         GtkWidget *progress_bar = gtk_bin_get_child(GTK_BIN(site->config));
 114         double percent = (double)cur/total;
 115         gtk_progress_bar_set_fraction(GTK_PROGRESS_BAR(progress_bar), MIN(percent, 1.0));
 116         gchar *msg = g_strdup_printf("Loading... %5.1f%% (%.2f/%.2f MB)",
 117                         percent*100, (double)cur/1000000, (double)total/1000000);
 118         gtk_progress_bar_set_text(GTK_PROGRESS_BAR(progress_bar), msg);
 119         g_free(msg);
 120 }
 121 gboolean _site_update_end(gpointer _site)
 122 {
 123         RadarSite *site = _site;
 124         if (site->message) {
 125                 g_warning("RadarSite: update_end - %s", site->message);
 126                 const char *fmt = "http://forecast.weather.gov/product.php?site=NWS&product=FTM&format=TXT&issuedby=%s";
 127                 char       *uri = g_strdup_printf(fmt, site->city->code+1);
 128                 GtkWidget  *box = gtk_box_new(GTK_ORIENTATION_VERTICAL, 0);
 129                 GtkWidget  *msg = gtk_label_new(site->message);
 130                 GtkWidget  *btn = gtk_link_button_new_with_label(uri, "View Radar Status");
 131                 gtk_box_set_homogeneous(GTK_BOX(box), TRUE);
 132                 gtk_box_pack_start(GTK_BOX(box), msg, TRUE, TRUE, 0);
 133                 gtk_box_pack_start(GTK_BOX(box), btn, TRUE, TRUE, 0);
 134                 aweather_bin_set_child(GTK_BIN(site->config), box);
 135                 g_free(uri);
 136         } else {
 137                 aweather_bin_set_child(GTK_BIN(site->config),
 138                                 aweather_level2_get_config(site->level2));
 139         }
 140         site->status = STATUS_LOADED;
 141         site->idle_source = 0;
 142         return FALSE;
 143 }
 144 gpointer _site_update_thread(gpointer _site)
 145 {
 146         RadarSite *site = _site;
 147         g_debug("RadarSite: update_thread - %s", site->city->code);
 148         site->message = NULL;
 149
 150         gboolean offline = grits_viewer_get_offline(site->viewer);
 151         gchar *nexrad_url = grits_prefs_get_string(site->prefs,
 152                         "aweather/nexrad_url", NULL);
 153
 154         /* Find nearest volume (temporally) */
 155         g_debug("RadarSite: update_thread - find nearest - %s", site->city->code);
 156         gchar *dir_list = g_strconcat(nexrad_url, "/", site->city->code,
 157                         "/", "dir.list", NULL);
 158         GList *files = grits_http_available(site->http,
 159                         "^\\w{4}_\\d{8}_\\d{4}$", site->city->code,
 160                         "\\d+ (.*)", (offline ? NULL : dir_list));
 161         g_free(dir_list);
 162         gchar *nearest = _find_nearest(site->time, files, 5);
 163         g_list_foreach(files, (GFunc)g_free, NULL);
 164         g_list_free(files);
 165         if (!nearest) {
 166                 site->message = "No suitable files found";
 167                 goto out;
 168         }
 169
 170         /* Fetch new volume */
 171         g_debug("RadarSite: update_thread - fetch");
 172         gchar *local = g_strconcat(site->city->code, "/", nearest, NULL);
 173         gchar *uri   = g_strconcat(nexrad_url, "/", local,   NULL);
 174         gchar *file  = grits_http_fetch(site->http, uri, local,
 175                         offline ? GRITS_LOCAL : GRITS_UPDATE,
 176                         _site_update_loading, site);
 177         g_free(nexrad_url);
 178         g_free(nearest);
 179         g_free(local);
 180         g_free(uri);
 181         if (!file) {
 182                 site->message = "Fetch failed";
 183                 goto out;
 184         }
 185
 186         /* Load and add new volume */
 187         g_debug("RadarSite: update_thread - load - %s", site->city->code);
 188         site->level2 = aweather_level2_new_from_file(
 189                         file, site->city->code, colormaps);
 190         g_free(file);
 191         if (!site->level2) {
 192                 site->message = "Load failed";
 193                 goto out;
 194         }
 195         grits_object_hide(GRITS_OBJECT(site->level2), site->hidden);
 196         grits_viewer_add(site->viewer, GRITS_OBJECT(site->level2),
 197                         GRITS_LEVEL_WORLD+3, TRUE);
 198
 199 out:
 200         if (!site->idle_source)
 201                 site->idle_source = g_idle_add(_site_update_end, site);
 202         return NULL;
 203 }
 204 void _site_update(RadarSite *site)
 205 {
 206         if (site->status == STATUS_LOADING)
 207                 return;
 208         site->status = STATUS_LOADING;
 209
 210         site->time = grits_viewer_get_time(site->viewer);
 211         g_debug("RadarSite: update %s - %d",
 212                         site->city->code, (gint)site->time);
 213
 214         /* Add a progress bar */
 215         GtkWidget *progress = gtk_progress_bar_new();
 216         gtk_progress_bar_set_text(GTK_PROGRESS_BAR(progress), "Loading...");
 217         aweather_bin_set_child(GTK_BIN(site->config), progress);
 218
 219         /* Remove old volume */
 220         g_debug("RadarSite: update - remove - %s", site->city->code);
 221         grits_object_destroy_pointer(&site->level2);
 222
 223         /* Fork loading right away so updating the
 224          * list of times doesn't take too long */
 225         g_thread_new("site-update-thread", _site_update_thread, site);
 226 }
 227
 228 /* RadarSite methods */
 229 void radar_site_unload(RadarSite *site)
 230 {
 231         if (site->status != STATUS_LOADED)
 232                 return; // Abort if it's still loading
 233
 234         g_debug("RadarSite: unload %s", site->city->code);
 235
 236         if (site->time_id)
 237                 g_signal_handler_disconnect(site->viewer, site->time_id);
 238         if (site->refresh_id)
 239                 g_signal_handler_disconnect(site->viewer, site->refresh_id);
 240         if (site->idle_source)
 241                 g_source_remove(site->idle_source);
 242         site->idle_source = 0;
 243
 244         /* Remove tab */
 245         if (site->config)
 246                 gtk_widget_destroy(site->config);
 247
 248         /* Remove radar */
 249         grits_object_destroy_pointer(&site->level2);
 250
 251         site->status = STATUS_UNLOADED;
 252 }
 253
 254 void radar_site_load(RadarSite *site)
 255 {
 256         g_debug("RadarSite: load %s", site->city->code);
 257
 258         /* Add tab page */
 259         site->config = gtk_alignment_new(0, 0, 1, 1);
 260         g_object_set_data(G_OBJECT(site->config), "site", site);
 261         gtk_notebook_append_page(GTK_NOTEBOOK(site->pconfig), site->config,
 262                         gtk_label_new(site->city->name));
 263         gtk_widget_show_all(site->config);
 264         if (gtk_notebook_get_current_page(GTK_NOTEBOOK(site->pconfig)) == 0)
 265                 gtk_notebook_set_current_page(GTK_NOTEBOOK(site->pconfig), -1);
 266
 267         /* Set up radar loading */
 268         site->time_id = g_signal_connect_swapped(site->viewer, "time-changed",
 269                         G_CALLBACK(_site_update), site);
 270         site->refresh_id = g_signal_connect_swapped(site->viewer, "refresh",
 271                         G_CALLBACK(_site_update), site);
 272         _site_update(site);
 273 }
 274
 275 void _site_on_location_changed(GritsViewer *viewer,
 276                 gdouble lat, gdouble lon, gdouble elev,
 277                 gpointer _site)
 278 {
 279         static gdouble min_dist = EARTH_R / 30;
 280         RadarSite *site = _site;
 281
 282         /* Calculate distance, could cache xyz values */
 283         gdouble eye_xyz[3], site_xyz[3];
 284         lle2xyz(lat, lon, elev, &eye_xyz[0], &eye_xyz[1], &eye_xyz[2]);
 285         lle2xyz(site->city->pos.lat, site->city->pos.lon, site->city->pos.elev,
 286                         &site_xyz[0], &site_xyz[1], &site_xyz[2]);
 287         gdouble dist = distd(site_xyz, eye_xyz);
 288
 289         /* Load or unload the site if necessasairy */
 290         if (dist <= min_dist && dist < elev*1.25 && site->status == STATUS_UNLOADED)
 291                 radar_site_load(site);
 292         else if (dist > 2*min_dist &&  site->status != STATUS_UNLOADED)
 293                 radar_site_unload(site);
 294 }
 295
 296 static gboolean on_marker_clicked(GritsObject *marker, GdkEvent *event, RadarSite *site)
 297 {
 298         GritsViewer *viewer = site->viewer;
 299         GritsPoint center = marker->center;
 300         grits_viewer_set_location(viewer, center.lat, center.lon, EARTH_R/35);
 301         grits_viewer_set_rotation(viewer, 0, 0, 0);
 302         /* Recursivly set notebook tabs */
 303         GtkWidget *widget, *parent;
 304         for (widget = site->config; widget; widget = parent) {
 305                 parent = gtk_widget_get_parent(widget);
 306                 if (GTK_IS_NOTEBOOK(parent)) {
 307                         gint i = gtk_notebook_page_num(GTK_NOTEBOOK(parent), widget);
 308                         gtk_notebook_set_current_page(GTK_NOTEBOOK(parent), i);
 309                 }
 310         }
 311         return TRUE;
 312 }
 313
 314 RadarSite *radar_site_new(city_t *city, GtkWidget *pconfig,
 315                 GritsViewer *viewer, GritsPrefs *prefs, GritsHttp *http)
 316 {
 317         RadarSite *site = g_new0(RadarSite, 1);
 318         site->viewer  = g_object_ref(viewer);
 319         site->prefs   = g_object_ref(prefs);
 320         //site->http    = http;
 321         site->http    = grits_http_new(G_DIR_SEPARATOR_S
 322                         "nexrad" G_DIR_SEPARATOR_S
 323                         "level2" G_DIR_SEPARATOR_S);
 324         site->city    = city;
 325         site->pconfig = pconfig;
 326         site->hidden  = TRUE;
 327
 328         /* Set initial location */
 329         gdouble lat, lon, elev;
 330         grits_viewer_get_location(viewer, &lat, &lon, &elev);
 331         _site_on_location_changed(viewer, lat, lon, elev, site);
 332
 333         /* Add marker */
 334         site->marker = grits_marker_new(site->city->name);
 335         GRITS_OBJECT(site->marker)->center = site->city->pos;
 336         GRITS_OBJECT(site->marker)->lod    = EARTH_R*0.75*site->city->lod;
 337         grits_viewer_add(site->viewer, GRITS_OBJECT(site->marker),
 338                         GRITS_LEVEL_HUD, FALSE);
 339         g_signal_connect(site->marker, "clicked",
 340                         G_CALLBACK(on_marker_clicked), site);
 341         grits_object_set_cursor(GRITS_OBJECT(site->marker), GDK_HAND2);
 342
 343         /* Connect signals */
 344         site->location_id  = g_signal_connect(viewer, "location-changed",
 345                         G_CALLBACK(_site_on_location_changed), site);
 346         return site;
 347 }
 348
 349 void radar_site_free(RadarSite *site)
 350 {
 351         radar_site_unload(site);
 352         grits_object_destroy_pointer(&site->marker);
 353         if (site->location_id)
 354                 g_signal_handler_disconnect(site->viewer, site->location_id);
 355         grits_http_free(site->http);
 356         g_object_unref(site->viewer);
 357         g_object_unref(site->prefs);
 358         g_free(site);
 359 }
 360
 361
 362 /**************
 363  * RadarConus *
 364  **************/
 365 #define CONUS_NORTH       50.406626367301044
 366 #define CONUS_WEST       -127.620375523875420
 367 #define CONUS_WIDTH       3400.0
 368 #define CONUS_HEIGHT      1600.0
 369 #define CONUS_DEG_PER_PX  0.017971305190311
 370
 371 struct _RadarConus {
 372         GritsViewer *viewer;
 373         GritsHttp   *http;
 374         GtkWidget   *config;
 375         time_t       time;
 376         const gchar *message;
 377         GMutex       loading;
 378
 379         gchar       *path;
 380         GritsTile   *tile[2];
 381
 382         guint        time_id;     // "time-changed"     callback ID
 383         guint        refresh_id;  // "refresh"          callback ID
 384         guint        idle_source; // _conus_update_end idle source
 385 };
 386
 387 void _conus_update_loading(gchar *file, goffset cur,
 388                 goffset total, gpointer _conus)
 389 {
 390         RadarConus *conus = _conus;
 391         GtkWidget *progress_bar = gtk_bin_get_child(GTK_BIN(conus->config));
 392         double percent = (double)cur/total;
 393         gtk_progress_bar_set_fraction(GTK_PROGRESS_BAR(progress_bar), MIN(percent, 1.0));
 394         gchar *msg = g_strdup_printf("Loading... %5.1f%% (%.2f/%.2f MB)",
 395                         percent*100, (double)cur/1000000, (double)total/1000000);
 396         gtk_progress_bar_set_text(GTK_PROGRESS_BAR(progress_bar), msg);
 397         g_free(msg);
 398 }
 399
 400 /* Copy images to graphics memory */
 401 static void _conus_update_end_copy(GritsTile *tile, guchar *pixels)
 402 {
 403         if (!tile->tex)
 404                 glGenTextures(1, &tile->tex);
 405
 406         gchar *clear = g_malloc0(2048*2048*4);
 407         glBindTexture(GL_TEXTURE_2D, tile->tex);
 408
 409         glPixelStorei(GL_UNPACK_ALIGNMENT, 1);
 410         glPixelStorei(GL_PACK_ALIGNMENT, 1);
 411         glTexImage2D(GL_TEXTURE_2D, 0, 4, 2048, 2048, 0,
 412                         GL_RGBA, GL_UNSIGNED_BYTE, clear);
 413         glTexSubImage2D(GL_TEXTURE_2D, 0, 1,1, CONUS_WIDTH/2,CONUS_HEIGHT,
 414                         GL_RGBA, GL_UNSIGNED_BYTE, pixels);
 415         tile->coords.n = 1.0/(CONUS_WIDTH/2);
 416         tile->coords.w = 1.0/ CONUS_HEIGHT;
 417         tile->coords.s = tile->coords.n +  CONUS_HEIGHT   / 2048.0;
 418         tile->coords.e = tile->coords.w + (CONUS_WIDTH/2) / 2048.0;
 419         glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_LINEAR);
 420         glTexParameterf(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_LINEAR);
 421         glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_CLAMP_TO_EDGE);
 422         glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_CLAMP_TO_EDGE);
 423         glFlush();
 424         g_free(clear);
 425 }
 426
 427 /* Split the pixbuf into east and west halves (with 2K sides)
 428  * Also map the pixbuf's alpha values */
 429 static void _conus_update_end_split(guchar *pixels, guchar *west, guchar *east,
 430                 gint width, gint height, gint pxsize)
 431 {
 432         g_debug("Conus: update_end_split");
 433         guchar *out[] = {west,east};
 434         const guchar alphamap[][4] = {
 435                 {0x04, 0xe9, 0xe7, 0x30},
 436                 {0x01, 0x9f, 0xf4, 0x60},
 437                 {0x03, 0x00, 0xf4, 0x90},
 438         };
 439         for (int y = 0; y < height; y++)
 440         for (int x = 0; x < width;  x++) {
 441                 gint subx = x % (width/2);
 442                 gint idx  = x / (width/2);
 443                 guchar *src = &pixels[(y*width+x)*pxsize];
 444                 guchar *dst = &out[idx][(y*(width/2)+subx)*4];
 445                 if (src[0] > 0xe0 &&
 446                     src[1] > 0xe0 &&
 447                     src[2] > 0xe0) {
 448                         dst[3] = 0x00;
 449                 } else {
 450                         dst[0] = src[0];
 451                         dst[1] = src[1];
 452                         dst[2] = src[2];
 453                         dst[3] = 0xff * 0.75;
 454                         for (int j = 0; j < G_N_ELEMENTS(alphamap); j++)
 455                                 if (src[0] == alphamap[j][0] &&
 456                                     src[1] == alphamap[j][1] &&
 457                                     src[2] == alphamap[j][2])
 458                                         dst[3] = alphamap[j][3];
 459                 }
 460         }
 461 }
 462
 463 gboolean _conus_update_end(gpointer _conus)
 464 {
 465         RadarConus *conus = _conus;
 466         g_debug("Conus: update_end");
 467
 468         /* Check error status */
 469         if (conus->message) {
 470                 g_warning("Conus: update_end - %s", conus->message);
 471                 aweather_bin_set_child(GTK_BIN(conus->config), gtk_label_new(conus->message));
 472                 goto out;
 473         }
 474
 475         /* Load and pixbuf */
 476         GError *error = NULL;
 477         GdkPixbuf *pixbuf = gdk_pixbuf_new_from_file(conus->path, &error);
 478         if (!pixbuf || error) {
 479                 g_warning("Conus: update_end - error loading pixbuf: %s", conus->path);
 480                 aweather_bin_set_child(GTK_BIN(conus->config), gtk_label_new("Error loading pixbuf"));
 481                 g_remove(conus->path);
 482                 goto out;
 483         }
 484
 485         /* Split pixels into east/west parts */
 486         guchar *pixels = gdk_pixbuf_get_pixels(pixbuf);
 487         gint    width  = gdk_pixbuf_get_width(pixbuf);
 488         gint    height = gdk_pixbuf_get_height(pixbuf);
 489         gint    pxsize = gdk_pixbuf_get_has_alpha(pixbuf) ? 4 : 3;
 490         guchar *pixels_west = g_malloc(4*(width/2)*height);
 491         guchar *pixels_east = g_malloc(4*(width/2)*height);
 492         _conus_update_end_split(pixels, pixels_west, pixels_east,
 493                         width, height, pxsize);
 494         g_object_unref(pixbuf);
 495
 496         /* Copy pixels to graphics memory */
 497         _conus_update_end_copy(conus->tile[0], pixels_west);
 498         _conus_update_end_copy(conus->tile[1], pixels_east);
 499         g_free(pixels_west);
 500         g_free(pixels_east);
 501
 502         /* Update GUI */
 503         gchar *label = g_path_get_basename(conus->path);
 504         aweather_bin_set_child(GTK_BIN(conus->config), gtk_label_new(label));
 505         grits_viewer_queue_draw(conus->viewer);
 506         g_free(label);
 507
 508 out:
 509         conus->idle_source = 0;
 510         g_free(conus->path);
 511         g_mutex_unlock(&conus->loading);
 512         return FALSE;
 513 }
 514
 515 gpointer _conus_update_thread(gpointer _conus)
 516 {
 517         RadarConus *conus = _conus;
 518         conus->message = NULL;
 519
 520         /* Find nearest */
 521         g_debug("Conus: update_thread - nearest");
 522         gboolean offline = grits_viewer_get_offline(conus->viewer);
 523         gchar *conus_url = "http://radar.weather.gov/Conus/RadarImg/";
 524         gchar *nearest;
 525         if (time(NULL) - conus->time < 60*60*5 && !offline) {
 526                 /* radar.weather.gov is full of lies.
 527                  * the index pages get cached and out of date */
 528                 /* gmtime is not thread safe, but it's not used very often so
 529                  * hopefully it'll be alright for now... :-( */
 530                 struct tm *tm = gmtime(&conus->time);
 531                 time_t onthe8 = conus->time - 60*((tm->tm_min+1)%10+1);
 532                 tm = gmtime(&onthe8);
 533                 nearest = g_strdup_printf("Conus_%04d%02d%02d_%02d%02d_N0Ronly.gif",
 534                                 tm->tm_year+1900, tm->tm_mon+1, tm->tm_mday,
 535                                 tm->tm_hour, tm->tm_min);
 536         } else {
 537                 GList *files = grits_http_available(conus->http,
 538                                 "^Conus_[^\"]*_N0Ronly.gif$", "", NULL, NULL);
 539                 nearest = _find_nearest(conus->time, files, 6);
 540                 g_list_foreach(files, (GFunc)g_free, NULL);
 541                 g_list_free(files);
 542                 if (!nearest) {
 543                         conus->message = "No suitable files";
 544                         goto out;
 545                 }
 546         }
 547
 548         /* Fetch the image */
 549         g_debug("Conus: update_thread - fetch");
 550         gchar *uri  = g_strconcat(conus_url, nearest, NULL);
 551         conus->path = grits_http_fetch(conus->http, uri, nearest,
 552                         offline ? GRITS_LOCAL : GRITS_ONCE,
 553                         _conus_update_loading, conus);
 554         g_free(nearest);
 555         g_free(uri);
 556         if (!conus->path) {
 557                 conus->message = "Fetch failed";
 558                 goto out;
 559         }
 560
 561 out:
 562         g_debug("Conus: update_thread - done");
 563         if (!conus->idle_source)
 564                 conus->idle_source = g_idle_add(_conus_update_end, conus);
 565         return NULL;
 566 }
 567
 568 void _conus_update(RadarConus *conus)
 569 {
 570         if (!g_mutex_trylock(&conus->loading))
 571                 return;
 572         conus->time = grits_viewer_get_time(conus->viewer);
 573         g_debug("Conus: update - %d",
 574                         (gint)conus->time);
 575
 576         /* Add a progress bar */
 577         GtkWidget *progress = gtk_progress_bar_new();
 578         gtk_progress_bar_set_text(GTK_PROGRESS_BAR(progress), "Loading...");
 579         aweather_bin_set_child(GTK_BIN(conus->config), progress);
 580
 581         g_thread_new("conus-update-thread", _conus_update_thread, conus);
 582 }
 583
 584 RadarConus *radar_conus_new(GtkWidget *pconfig,
 585                 GritsViewer *viewer, GritsHttp *http)
 586 {
 587         RadarConus *conus = g_new0(RadarConus, 1);
 588         conus->viewer  = g_object_ref(viewer);
 589         conus->http    = http;
 590         conus->config  = gtk_alignment_new(0, 0, 1, 1);
 591         g_mutex_init(&conus->loading);
 592
 593         gdouble south =  CONUS_NORTH - CONUS_DEG_PER_PX*CONUS_HEIGHT;
 594         gdouble east  =  CONUS_WEST  + CONUS_DEG_PER_PX*CONUS_WIDTH;
 595         gdouble mid   =  CONUS_WEST  + CONUS_DEG_PER_PX*CONUS_WIDTH/2;
 596         conus->tile[0] = grits_tile_new(NULL, CONUS_NORTH, south, mid, CONUS_WEST);
 597         conus->tile[1] = grits_tile_new(NULL, CONUS_NORTH, south, east, mid);
 598         conus->tile[0]->zindex = 2;
 599         conus->tile[1]->zindex = 1;
 600         grits_viewer_add(viewer, GRITS_OBJECT(conus->tile[0]), GRITS_LEVEL_WORLD+2, FALSE);
 601         grits_viewer_add(viewer, GRITS_OBJECT(conus->tile[1]), GRITS_LEVEL_WORLD+2, FALSE);
 602
 603         conus->time_id = g_signal_connect_swapped(viewer, "time-changed",
 604                         G_CALLBACK(_conus_update), conus);
 605         conus->refresh_id = g_signal_connect_swapped(viewer, "refresh",
 606                         G_CALLBACK(_conus_update), conus);
 607
 608         g_object_set_data(G_OBJECT(conus->config), "conus", conus);
 609         gtk_notebook_append_page(GTK_NOTEBOOK(pconfig), conus->config,
 610                         gtk_label_new("Conus"));
 611
 612         _conus_update(conus);
 613         return conus;
 614 }
 615
 616 void radar_conus_free(RadarConus *conus)
 617 {
 618         g_signal_handler_disconnect(conus->viewer, conus->time_id);
 619         g_signal_handler_disconnect(conus->viewer, conus->refresh_id);
 620         if (conus->idle_source)
 621                 g_source_remove(conus->idle_source);
 622
 623         for (int i = 0; i < 2; i++)
 624                 grits_object_destroy_pointer(&conus->tile[i]);
 625
 626         g_object_unref(conus->viewer);
 627         g_free(conus);
 628 }
 629
 630
 631 /********************
 632  * GritsPluginRadar *
 633  ********************/
 634 static void _draw_hud(GritsCallback *callback, GritsOpenGL *opengl, gpointer _self)
 635 {
 636         g_debug("GritsPluginRadar: _draw_hud");
 637         /* Setup OpenGL */
 638         glMatrixMode(GL_MODELVIEW ); glLoadIdentity();
 639         glMatrixMode(GL_PROJECTION); glLoadIdentity();
 640         glDisable(GL_TEXTURE_2D);
 641         glDisable(GL_ALPHA_TEST);
 642         glDisable(GL_CULL_FACE);
 643         glDisable(GL_LIGHTING);
 644         glEnable(GL_COLOR_MATERIAL);
 645
 646         GHashTableIter iter;
 647         gpointer name, _site;
 648         GritsPluginRadar *self = GRITS_PLUGIN_RADAR(_self);
 649         g_hash_table_iter_init(&iter, self->sites);
 650         while (g_hash_table_iter_next(&iter, &name, &_site)) {
 651                 /* Pick correct colormaps */
 652                 RadarSite *site = _site;
 653                 if (site->hidden || !site->level2)
 654                         continue;
 655                 AWeatherColormap *colormap = site->level2->sweep_colors;
 656
 657                 /* Print the color table */
 658                 glBegin(GL_QUADS);
 659                 int len = colormap->len;
 660                 for (int i = 0; i < len; i++) {
 661                         glColor4ubv(colormap->data[i]);
 662                         glVertex3f(-1.0, (float)((i  ) - len/2)/(len/2), 0.0); // bot left
 663                         glVertex3f(-1.0, (float)((i+1) - len/2)/(len/2), 0.0); // top left
 664                         glVertex3f(-0.9, (float)((i+1) - len/2)/(len/2), 0.0); // top right
 665                         glVertex3f(-0.9, (float)((i  ) - len/2)/(len/2), 0.0); // bot right
 666                 }
 667                 glEnd();
 668         }
 669 }
 670
 671 static void _load_colormap(gchar *filename, AWeatherColormap *cm)
 672 {
 673         g_debug("GritsPluginRadar: _load_colormap - %s", filename);
 674         FILE *file = fopen(filename, "r");
 675         if (!file)
 676                 g_error("GritsPluginRadar: open failed");
 677         guint8 color[4];
 678         GArray *array = g_array_sized_new(FALSE, TRUE, sizeof(color), 256);
 679         if (!fgets(cm->name, sizeof(cm->name), file)) goto out;
 680         if (!fscanf(file, "%f\n", &cm->scale))        goto out;
 681         if (!fscanf(file, "%f\n", &cm->shift))        goto out;
 682         int r, g, b, a;
 683         while (fscanf(file, "%d %d %d %d\n", &r, &g, &b, &a) == 4) {
 684                 color[0] = r;
 685                 color[1] = g;
 686                 color[2] = b;
 687                 color[3] = a;
 688                 g_array_append_val(array, color);
 689         }
 690         cm->len  = (gint )array->len;
 691         cm->data = (void*)array->data;
 692 out:
 693         g_array_free(array, FALSE);
 694         fclose(file);
 695 }
 696
 697 static void _update_hidden(GtkNotebook *notebook,
 698                 gpointer _, guint page_num, gpointer viewer)
 699 {
 700         g_debug("GritsPluginRadar: _update_hidden - 0..%d = %d",
 701                         gtk_notebook_get_n_pages(notebook), page_num);
 702
 703         for (gint i = 0; i < gtk_notebook_get_n_pages(notebook); i++) {
 704                 gboolean is_hidden = (i != page_num);
 705                 GtkWidget  *config = gtk_notebook_get_nth_page(notebook, i);
 706                 RadarConus *conus  = g_object_get_data(G_OBJECT(config), "conus");
 707                 RadarSite  *site   = g_object_get_data(G_OBJECT(config), "site");
 708
 709                 /* Conus */
 710                 if (conus) {
 711                         grits_object_hide(GRITS_OBJECT(conus->tile[0]), is_hidden);
 712                         grits_object_hide(GRITS_OBJECT(conus->tile[1]), is_hidden);
 713                 } else if (site) {
 714                         site->hidden = is_hidden;
 715                         if (site->level2)
 716                                 grits_object_hide(GRITS_OBJECT(site->level2), is_hidden);
 717                 } else {
 718                         g_warning("GritsPluginRadar: _update_hidden - no site or counus found");
 719                 }
 720         }
 721         grits_viewer_queue_draw(viewer);
 722 }
 723
 724 /* Methods */
 725 GritsPluginRadar *grits_plugin_radar_new(GritsViewer *viewer, GritsPrefs *prefs)
 726 {
 727         /* TODO: move to constructor if possible */
 728         g_debug("GritsPluginRadar: new");
 729         GritsPluginRadar *self = g_object_new(GRITS_TYPE_PLUGIN_RADAR, NULL);
 730         self->viewer = g_object_ref(viewer);
 731         self->prefs  = g_object_ref(prefs);
 732
 733         /* Setup page switching */
 734         self->tab_id = g_signal_connect(self->config, "switch-page",
 735                         G_CALLBACK(_update_hidden), viewer);
 736
 737         /* Load HUD */
 738         self->hud = grits_callback_new(_draw_hud, self);
 739         grits_viewer_add(viewer, GRITS_OBJECT(self->hud), GRITS_LEVEL_HUD, FALSE);
 740
 741         /* Load Conus */
 742         self->conus = radar_conus_new(self->config, self->viewer, self->conus_http);
 743
 744         /* Load radar sites */
 745         for (city_t *city = cities; city->type; city++) {
 746                 if (city->type != LOCATION_CITY)
 747                         continue;
 748                 RadarSite *site = radar_site_new(city, self->config,
 749                                 self->viewer, self->prefs, self->sites_http);
 750                 g_hash_table_insert(self->sites, city->code, site);
 751         }
 752
 753         return self;
 754 }
 755
 756 static GtkWidget *grits_plugin_radar_get_config(GritsPlugin *_self)
 757 {
 758         GritsPluginRadar *self = GRITS_PLUGIN_RADAR(_self);
 759         return self->config;
 760 }
 761
 762 /* GObject code */
 763 static void grits_plugin_radar_plugin_init(GritsPluginInterface *iface);
 764 G_DEFINE_TYPE_WITH_CODE(GritsPluginRadar, grits_plugin_radar, G_TYPE_OBJECT,
 765                 G_IMPLEMENT_INTERFACE(GRITS_TYPE_PLUGIN,
 766                         grits_plugin_radar_plugin_init));
 767 static void grits_plugin_radar_plugin_init(GritsPluginInterface *iface)
 768 {
 769         g_debug("GritsPluginRadar: plugin_init");
 770         /* Add methods to the interface */
 771         iface->get_config = grits_plugin_radar_get_config;
 772 }
 773 static void grits_plugin_radar_init(GritsPluginRadar *self)
 774 {
 775         g_debug("GritsPluginRadar: class_init");
 776         /* Set defaults */
 777         self->sites_http = grits_http_new(G_DIR_SEPARATOR_S
 778                         "nexrad" G_DIR_SEPARATOR_S
 779                         "level2" G_DIR_SEPARATOR_S);
 780         self->conus_http = grits_http_new(G_DIR_SEPARATOR_S
 781                         "nexrad" G_DIR_SEPARATOR_S
 782                         "conus"  G_DIR_SEPARATOR_S);
 783         self->sites      = g_hash_table_new_full(g_str_hash, g_str_equal,
 784                                 NULL, (GDestroyNotify)radar_site_free);
 785         self->config     = g_object_ref(gtk_notebook_new());
 786
 787         /* Load colormaps */
 788         for (int i = 0; colormaps[i].file; i++) {
 789                 gchar *file = g_build_filename(PKGDATADIR,
 790                                 "colors", colormaps[i].file, NULL);
 791                 _load_colormap(file, &colormaps[i]);
 792                 g_free(file);
 793         }
 794
 795         /* Need to position on the top because of Win32 bug */
 796         gtk_notebook_set_tab_pos(GTK_NOTEBOOK(self->config), GTK_POS_LEFT);
 797 }
 798 static void grits_plugin_radar_dispose(GObject *gobject)
 799 {
 800         g_debug("GritsPluginRadar: dispose");
 801         GritsPluginRadar *self = GRITS_PLUGIN_RADAR(gobject);
 802         if (self->viewer) {
 803                 GritsViewer *viewer = self->viewer;
 804                 self->viewer = NULL;
 805                 g_signal_handler_disconnect(self->config, self->tab_id);
 806                 grits_object_destroy_pointer(&self->hud);
 807                 radar_conus_free(self->conus);
 808                 g_hash_table_destroy(self->sites);
 809                 g_object_unref(self->config);
 810                 g_object_unref(self->prefs);
 811                 g_object_unref(viewer);
 812         }
 813         /* Drop references */
 814         G_OBJECT_CLASS(grits_plugin_radar_parent_class)->dispose(gobject);
 815 }
 816 static void grits_plugin_radar_finalize(GObject *gobject)
 817 {
 818         g_debug("GritsPluginRadar: finalize");
 819         GritsPluginRadar *self = GRITS_PLUGIN_RADAR(gobject);
 820         /* Free data */
 821         grits_http_free(self->conus_http);
 822         grits_http_free(self->sites_http);
 823         gtk_widget_destroy(self->config);
 824         G_OBJECT_CLASS(grits_plugin_radar_parent_class)->finalize(gobject);
 825
 826 }
 827 static void grits_plugin_radar_class_init(GritsPluginRadarClass *klass)
 828 {
 829         g_debug("GritsPluginRadar: class_init");
 830         GObjectClass *gobject_class = (GObjectClass*)klass;
 831         gobject_class->dispose  = grits_plugin_radar_dispose;
 832         gobject_class->finalize = grits_plugin_radar_finalize;
 833 }