]> Pileus Git - aweather/blobdiff - src/plugins/radar.c
Update for libgis
[aweather] / src / plugins / radar.c
index 334e0fb6b8f4598ccbaede664d220770e9bd45f2..4efa6cb57a126a408ed0102d078cb67283340286 100644 (file)
@@ -1,5 +1,5 @@
 /*
- * Copyright (C) 2009 Andy Spencer <spenceal@rose-hulman.edu>
+ * Copyright (C) 2009-2010 Andy Spencer <andy753421@gmail.com>
  *
  * This program is free software: you can redistribute it and/or modify
  * it under the terms of the GNU General Public License as published by
 #include <math.h>
 #include <rsl.h>
 
-#include <gis/gis.h>
+#include <gis.h>
 
 #include "radar.h"
 #include "marching.h"
+#include "../aweather-location.h"
 
 
 /**************************
@@ -55,7 +56,7 @@ static void _bscan_sweep(GisPluginRadar *self, Sweep *sweep, colormap_t *colorma
                        buf[buf_i+0] = colormap->data[val][0];
                        buf[buf_i+1] = colormap->data[val][1];
                        buf[buf_i+2] = colormap->data[val][2];
-                       buf[buf_i+3] = colormap->data[val][3]; // TESTING
+                       buf[buf_i+3] = colormap->data[val][3]*0.75; // TESTING
                        if (val == BADVAL     || val == RFVAL      || val == APFLAG ||
                            val == NOTFOUND_H || val == NOTFOUND_V || val == NOECHO) {
                                buf[buf_i+3] = 0x00; // transparent
@@ -313,10 +314,6 @@ static void _load_radar(GisPluginRadar *self, gchar *radar_file)
        _load_radar_gui(self, radar);
 }
 
-
-/*****************
- * ASync helpers *
- *****************/
 typedef struct {
        GisPluginRadar *self;
        gchar *radar_file;
@@ -396,23 +393,13 @@ static void _cache_done_cb(char *path, gboolean updated, gpointer _self)
        self->soup = NULL;
 }
 
-
-/*************
- * Callbacks *
- *************/
-static void _on_sweep_clicked(GtkRadioButton *button, gpointer _self)
+static void _set_radar(GisPluginRadar *self, const char *site, const char *time)
 {
-       GisPluginRadar *self = GIS_PLUGIN_RADAR(_self);
-       _load_colormap(self, g_object_get_data(G_OBJECT(button), "type" ));
-       _load_sweep   (self, g_object_get_data(G_OBJECT(button), "sweep"));
-}
+       if (!site || !time)
+               return;
+       g_debug("GisPluginRadar: set_radar - %s - %s", site, time);
 
-static void _on_time_changed(GisViewer *viewer, const char *time, gpointer _self)
-{
-       GisPluginRadar *self = GIS_PLUGIN_RADAR(_self);
-       g_debug("GisPluginRadar: on_time_changed - setting time=%s", time);
        // format: http://mesonet.agron.iastate.edu/data/nexrd2/raw/KABR/KABR_20090510_0323
-       char *site = gis_viewer_get_site(viewer);
        char *path = g_strdup_printf("nexrd2/raw/%s/%s_%s", site, site, time);
 
        /* Set up progress bar */
@@ -450,9 +437,62 @@ static void _on_time_changed(GisViewer *viewer, const char *time, gpointer _self
        g_free(path);
 }
 
+
+/*************
+ * Callbacks *
+ *************/
+static void _on_sweep_clicked(GtkRadioButton *button, gpointer _self)
+{
+       GisPluginRadar *self = GIS_PLUGIN_RADAR(_self);
+       _load_colormap(self, g_object_get_data(G_OBJECT(button), "type" ));
+       _load_sweep   (self, g_object_get_data(G_OBJECT(button), "sweep"));
+}
+
+static void _on_time_changed(GisViewer *viewer, const char *time, gpointer _self)
+{
+       g_debug("GisPluginRadar: _on_time_changed");
+       GisPluginRadar *self = GIS_PLUGIN_RADAR(_self);
+       g_free(self->cur_time);
+       self->cur_time = g_strdup(time);
+       _set_radar(self, self->cur_site, self->cur_time);
+}
+
+static void _on_location_changed(GisViewer *viewer,
+               gdouble lat, gdouble lon, gdouble elev,
+               gpointer _self)
+{
+       g_debug("GisPluginRadar: _on_location_changed");
+       GisPluginRadar *self = GIS_PLUGIN_RADAR(_self);
+       gdouble min_dist = EARTH_R / 5;
+       city_t *city, *min_city = NULL;
+       for (city = cities; city->type; city++) {
+               if (city->type != LOCATION_CITY)
+                       continue;
+               gdouble city_loc[3] = {};
+               gdouble eye_loc[3]  = {lat, lon, elev};
+               lle2xyz(city->lat, city->lon, city->elev,
+                               &city_loc[0], &city_loc[1], &city_loc[2]);
+               lle2xyz(lat, lon, elev,
+                               &eye_loc[0], &eye_loc[1], &eye_loc[2]);
+               gdouble dist = distd(city_loc, eye_loc);
+               if (dist < min_dist) {
+                       min_dist = dist;
+                       min_city = city;
+               }
+       }
+       static city_t *last_city = NULL;
+       if (min_city && min_city != last_city) {
+               self->cur_site = min_city->code;
+               _set_radar(self, self->cur_site, self->cur_time);
+       }
+       last_city = min_city;
+}
+
 static gpointer _draw_radar(GisCallback *callback, gpointer _self)
 {
        GisPluginRadar *self = GIS_PLUGIN_RADAR(_self);
+
+       /* Draw wsr88d */
        if (self->cur_sweep == NULL)
                return NULL;
        g_debug("GisPluginRadar: _draw_radar");
@@ -497,14 +537,16 @@ static gpointer _draw_radar(GisCallback *callback, gpointer _self)
 
        glDisable(GL_ALPHA_TEST);
        glDisable(GL_CULL_FACE);
-       glDisable(GL_DEPTH_TEST);
        glDisable(GL_LIGHTING);
        glEnable(GL_TEXTURE_2D);
+       glEnable(GL_POLYGON_OFFSET_FILL);
+       glPolygonOffset(1.0, -2.0);
        glColor4f(1,1,1,1);
 
        /* Draw the rays */
        glBindTexture(GL_TEXTURE_2D, self->cur_sweep_tex);
-       glBegin(GL_QUAD_STRIP);
+       g_message("Tex = %d", self->cur_sweep_tex);
+       glBegin(GL_TRIANGLE_STRIP);
        for (int ri = 0; ri <= sweep->h.nrays; ri++) {
                Ray  *ray = NULL;
                double angle = 0;
@@ -590,14 +632,26 @@ GisPluginRadar *gis_plugin_radar_new(GisViewer *viewer, GisPrefs *prefs)
        self->prefs  = prefs;
        self->time_changed_id = g_signal_connect(viewer, "time-changed",
                        G_CALLBACK(_on_time_changed), self);
-       /* Add renderers */
-       GisCallback *callback;
+       self->location_changed_id = g_signal_connect(viewer, "location-changed",
+                       G_CALLBACK(_on_location_changed), self);
+
+       for (city_t *city = cities; city->type; city++) {
+               if (city->type != LOCATION_CITY)
+                       continue;
+               g_debug("Adding marker for %s %s", city->code, city->label);
+               GisMarker *marker = gis_marker_new(city->label);
+               gis_point_set_lle(gis_object_center(GIS_OBJECT(marker)),
+                               city->lat, city->lon, city->elev);
+               GIS_OBJECT(marker)->lod = EARTH_R/2;
+               gis_viewer_add(self->viewer, GIS_OBJECT(marker), GIS_LEVEL_OVERLAY, FALSE);
+       }
 
-       callback = gis_callback_new(_draw_radar, self);
-       gis_viewer_add(viewer, GIS_OBJECT(callback), GIS_LEVEL_WORLD, TRUE);
+       /* Add renderers */
+       GisCallback *radar_cb = gis_callback_new(_draw_radar, self);
+       GisCallback *hud_cb   = gis_callback_new(_draw_hud, self);
 
-       callback = gis_callback_new(_draw_hud, self);
-       gis_viewer_add(viewer, GIS_OBJECT(callback), GIS_LEVEL_HUD, FALSE);
+       gis_viewer_add(viewer, GIS_OBJECT(radar_cb),    GIS_LEVEL_WORLD, TRUE);
+       gis_viewer_add(viewer, GIS_OBJECT(hud_cb),      GIS_LEVEL_HUD,   FALSE);
 
        return self;
 }