src/google.c

   1 /*
   2  * viking -- GPS Data and Topo Analyzer, Explorer, and Manager
   3  *
   4  * Copyright (C) 2003-2005, Evan Battaglia <gtoevan@gmx.net>
   5  *
   6  * This program is free software; you can redistribute it and/or modify
   7  * it under the terms of the GNU General Public License as published by
   8  * the Free Software Foundation; either version 2 of the License, or
   9  * (at your option) any later version.
  10  *
  11  * This program is distributed in the hope that it will be useful,
  12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  14  * GNU General Public License for more details.
  15  *
  16  * You should have received a copy of the GNU General Public License
  17  * along with this program; if not, write to the Free Software
  18  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
  19  *
  20  */
  21
  22 #include <gtk/gtk.h>
  23 #include <math.h>
  24 #include <string.h>
  25 #include "coords.h"
  26 #include "vikcoord.h"
  27 #include "mapcoord.h"
  28 #include "download.h"
  29 #include "globals.h"
  30 #include "google.h"
  31 #include "vikmapslayer.h"
  32
  33 #define GOOGLE_VERSION "w2.46"
  34 #define GOOGLE_TRANS_VERSION "w2t.47"
  35 #define GOOGLE_KH_VERSION "17"
  36
  37 static int google_download ( MapCoord *src, const gchar *dest_fn );
  38 static int google_trans_download ( MapCoord *src, const gchar *dest_fn );
  39 static int google_kh_download ( MapCoord *src, const gchar *dest_fn );
  40 static void google_mapcoord_to_center_coord ( MapCoord *src, VikCoord *dest );
  41 static gboolean google_coord_to_mapcoord ( const VikCoord *src, gdouble xzoom, gdouble yzoom, MapCoord *dest );
  42
  43 void google_init () {
  44   VikMapsLayer_MapType google_1 = { 7, 256, 256, VIK_VIEWPORT_DRAWMODE_MERCATOR, google_coord_to_mapcoord, google_mapcoord_to_center_coord, google_download };
  45   VikMapsLayer_MapType google_2 = { 10, 256, 256, VIK_VIEWPORT_DRAWMODE_MERCATOR, google_coord_to_mapcoord, google_mapcoord_to_center_coord, google_trans_download };
  46   VikMapsLayer_MapType google_3 = { 11, 256, 256, VIK_VIEWPORT_DRAWMODE_MERCATOR, google_coord_to_mapcoord, google_mapcoord_to_center_coord, google_kh_download };
  47
  48   maps_layer_register_type("Google Maps", 7, &google_1);
  49   maps_layer_register_type("Transparent Google Maps", 10, &google_2);
  50   maps_layer_register_type("Google Satellite Images", 11, &google_3);
  51 }
  52
  53 /* 1 << (x) is like a 2**(x) */
  54 #define GZ(x) ((1<<x))
  55
  56 static const gdouble scale_mpps[] = { GZ(0), GZ(1), GZ(2), GZ(3), GZ(4), GZ(5), GZ(6), GZ(7), GZ(8), GZ(9),
  57                                            GZ(10), GZ(11), GZ(12), GZ(13), GZ(14), GZ(15), GZ(16), GZ(17) };
  58
  59 static const gint num_scales = (sizeof(scale_mpps) / sizeof(scale_mpps[0]));
  60
  61 #define ERROR_MARGIN 0.01
  62 guint8 google_zoom ( gdouble mpp ) {
  63   gint i;
  64   for ( i = 0; i < num_scales; i++ ) {
  65     if ( ABS(scale_mpps[i] - mpp) < ERROR_MARGIN )
  66       return i;
  67   }
  68   return 255;
  69 }
  70
  71 gboolean google_coord_to_mapcoord ( const VikCoord *src, gdouble xzoom, gdouble yzoom, MapCoord *dest )
  72 {
  73   g_assert ( src->mode == VIK_COORD_LATLON );
  74
  75   if ( xzoom != yzoom )
  76     return FALSE;
  77
  78   dest->scale = google_zoom ( xzoom );
  79   if ( dest->scale == 255 )
  80     return FALSE;
  81
  82   dest->x = (src->east_west + 180) / 360 * GZ(17) / xzoom;
  83   dest->y = (180 - MERCLAT(src->north_south)) / 360 * GZ(17) / xzoom;
  84   dest->z = 0;
  85
  86   return TRUE;
  87 }
  88
  89 void google_mapcoord_to_center_coord ( MapCoord *src, VikCoord *dest )
  90 {
  91   gdouble socalled_mpp = GZ(src->scale);
  92   dest->mode = VIK_COORD_LATLON;
  93   dest->east_west = ((src->x+0.5) / GZ(17) * socalled_mpp * 360) - 180;
  94   dest->north_south = DEMERCLAT(180 - ((src->y+0.5) / GZ(17) * socalled_mpp * 360));
  95 }
  96
  97 static int real_google_download ( MapCoord *src, const gchar *dest_fn, const char *verstr )
  98 {
  99    int res;
 100    gchar *uri = g_strdup_printf ( "/mt?v=%s&x=%d&y=%d&zoom=%d", verstr, src->x, src->y, src->scale );
 101    res = a_http_download_get_url ( "mt.google.com", uri, dest_fn );
 102    g_free ( uri );
 103    return res;
 104 }
 105
 106 static int google_download ( MapCoord *src, const gchar *dest_fn )
 107 {
 108    return(real_google_download ( src, dest_fn, GOOGLE_VERSION ));
 109 }
 110
 111 static int google_trans_download ( MapCoord *src, const gchar *dest_fn )
 112 {
 113    return(real_google_download ( src, dest_fn, GOOGLE_TRANS_VERSION ));
 114 }
 115
 116 static char *kh_encode(guint32 x, guint32 y, guint8 scale)
 117 {
 118   gchar *buf = g_malloc ( (20-scale)*sizeof(gchar) );
 119   guint32 ya = 1 << (17 - scale);
 120   gint8 i, j;
 121
 122   if (y < 0 || (ya-1 < y)) {
 123     strcpy(buf,"tqq"); /* BAD */
 124     return buf;
 125   }
 126   if (x < 0 || ya-1 < x) {
 127     x %= ya;
 128     if (x < 0)
 129       x += ya;
 130   }
 131
 132   buf[0] = 't';
 133   for (j = 1, i = 16; i >= scale; i--, j++) {
 134     ya /= 2;
 135     if (y < ya) {
 136       if (x < ya)
 137         buf[j]='q';
 138       else {
 139         buf[j]='r';
 140         x -= ya;
 141       }
 142     } else {
 143       if (x < ya) {
 144         buf[j] = 't';
 145         y -= ya;
 146       } else {
 147         buf[j] = 's';
 148         x -= ya;
 149         y -= ya;
 150       }
 151     }
 152   }
 153   buf[j] = '\0';
 154   return buf;
 155 }
 156
 157 static int google_kh_download ( MapCoord *src, const gchar *dest_fn )
 158 {
 159    int res;
 160    gchar *khenc = kh_encode( src->x, src->y, src->scale );
 161    gchar *uri = g_strdup_printf ( "/kh?v=%s&t=%s", GOOGLE_KH_VERSION, khenc );
 162    g_free ( khenc );
 163    res = a_http_download_get_url ( "kh.google.com", uri, dest_fn );
 164    g_free ( uri );
 165    return(res);
 166 }