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

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