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 <stdlib.h>
  23 #include <string.h>
  24 #include <gtk/gtk.h>
  25 #include <math.h>
  26 #include <string.h>
  27 #include "coords.h"
  28 #include "vikcoord.h"
  29 #include "mapcoord.h"
  30 #include "download.h"
  31 #include "curl_download.h"
  32 #include "globals.h"
  33 #include "google.h"
  34 #include "vikmapslayer.h"
  35
  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 static DownloadOptions google_options = { "http://maps.google.com/", 0 };
  44
  45 void google_init () {
  46   VikMapsLayer_MapType google_1 = { 7, 256, 256, VIK_VIEWPORT_DRAWMODE_MERCATOR, google_coord_to_mapcoord, google_mapcoord_to_center_coord, google_download };
  47   VikMapsLayer_MapType google_2 = { 10, 256, 256, VIK_VIEWPORT_DRAWMODE_MERCATOR, google_coord_to_mapcoord, google_mapcoord_to_center_coord, google_trans_download };
  48   VikMapsLayer_MapType google_3 = { 11, 256, 256, VIK_VIEWPORT_DRAWMODE_MERCATOR, google_coord_to_mapcoord, google_mapcoord_to_center_coord, google_kh_download };
  49
  50   maps_layer_register_type("Google Maps", 7, &google_1);
  51   maps_layer_register_type("Transparent Google Maps", 10, &google_2);
  52   maps_layer_register_type("Google Satellite Images", 11, &google_3);
  53 }
  54
  55 /* 1 << (x) is like a 2**(x) */
  56 #define GZ(x) ((1<<x))
  57
  58 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),
  59                                            GZ(10), GZ(11), GZ(12), GZ(13), GZ(14), GZ(15), GZ(16), GZ(17) };
  60
  61 static const gint num_scales = (sizeof(scale_mpps) / sizeof(scale_mpps[0]));
  62
  63 #define ERROR_MARGIN 0.01
  64 guint8 google_zoom ( gdouble mpp ) {
  65   gint i;
  66   for ( i = 0; i < num_scales; i++ ) {
  67     if ( ABS(scale_mpps[i] - mpp) < ERROR_MARGIN )
  68       return i;
  69   }
  70   return 255;
  71 }
  72
  73 typedef enum {
  74         TYPE_GOOGLE_MAPS = 0,
  75         TYPE_GOOGLE_TRANS,
  76         TYPE_GOOGLE_SAT,
  77
  78         TYPE_GOOGLE_NUM
  79 } GoogleType;
  80
  81 static gchar *parse_version_number(gchar *text)
  82 {
  83   int i;
  84   gchar *vers;
  85   gchar *s = text;
  86
  87   for (i = 0; (s[i] != '\\') && (i < 8); i++)
  88     ;
  89   if (s[i] != '\\') {
  90     return NULL;
  91   }
  92
  93   return vers = g_strndup(s, i);
  94 }
  95
  96 static const gchar *google_version_number(MapCoord *mapcoord, GoogleType google_type)
  97 {
  98   static gboolean first = TRUE;
  99   static char *vers[] = { "w2.60", "w2t.60", "20" };
 100   FILE *tmp_file;
 101   int tmp_fd;
 102   gchar *tmpname;
 103   gchar *uri;
 104   VikCoord coord;
 105   gchar *text, *pat, *beg;
 106   GMappedFile *mf;
 107   gsize len;
 108   gchar *gvers, *tvers, *kvers, *tmpvers;
 109   static DownloadOptions dl_options = { "http://maps.google.com/", 0 };
 110   static const char *gvers_pat = "http://mt0.google.com/mt?n\\x3d404\\x26v\\x3d";
 111   static const char *kvers_pat = "http://kh0.google.com/kh?n\\x3d404\\x26v\\x3d";
 112
 113   g_assert(google_type < TYPE_GOOGLE_NUM);
 114
 115   if (!first)
 116     return (vers[google_type]);
 117
 118
 119   first = FALSE;
 120   gvers = tvers = kvers = NULL;
 121   if ((tmp_fd = g_file_open_tmp ("vikgvers.XXXXXX", &tmpname, NULL)) == -1) {
 122     g_critical("couldn't open temp file %s\n", tmpname);
 123     exit(1);
 124   }
 125
 126   google_mapcoord_to_center_coord(mapcoord, &coord);
 127   uri = g_strdup_printf("http://maps.google.com/maps?f=q&hl=en&q=%f,%f", coord.north_south, coord.east_west);
 128   tmp_file = fdopen(tmp_fd, "r+");
 129
 130   if (curl_download_uri(uri, tmp_file, &dl_options)) {  /* error */
 131     g_warning("Failed downloading %s\n", tmpname);
 132   } else {
 133     if ((mf = g_mapped_file_new(tmpname, FALSE, NULL)) == NULL) {
 134       g_critical("couldn't map temp file\n");
 135       exit(1);
 136     }
 137     len = g_mapped_file_get_length(mf);
 138     text = g_mapped_file_get_contents(mf);
 139
 140     if ((beg = g_strstr_len(text, len, "GLoadApi")) == NULL) {
 141       g_warning("Failed fetching Google numbers (\"GLoadApi\" not found)\n");
 142       goto failed;
 143     }
 144
 145     pat = beg;
 146     while (!gvers || !tvers) {
 147       if ((pat = g_strstr_len(pat, &text[len] - pat, gvers_pat)) != NULL) {
 148         pat += strlen(gvers_pat);
 149         if ((tmpvers = parse_version_number(pat)) != NULL) {
 150           if (strstr(tmpvers, "t."))
 151             tvers = tmpvers;
 152           else
 153             gvers = tmpvers;
 154         }
 155       }
 156       else
 157         break;
 158     }
 159
 160     if ((pat = g_strstr_len(beg, &text[len] - beg, kvers_pat)) != NULL)
 161         kvers = parse_version_number(pat + strlen(kvers_pat));
 162
 163     if (gvers && tvers && kvers) {
 164       vers[TYPE_GOOGLE_MAPS] = gvers;
 165       vers[TYPE_GOOGLE_TRANS] = tvers;
 166       vers[TYPE_GOOGLE_SAT] = kvers;
 167     }
 168     else
 169       g_warning("Failed getting google version numbers");
 170
 171     if (gvers)
 172       fprintf(stderr, "DEBUG gvers=%s\n", gvers);
 173     if (tvers)
 174       fprintf(stderr, "DEBUG tvers=%s\n", tvers);
 175     if (kvers)
 176       fprintf(stderr, "DEBUG kvers=%s\n", kvers);
 177
 178 failed:
 179     g_mapped_file_free(mf);
 180   }
 181
 182   fclose(tmp_file);
 183   g_free(tmpname);
 184   return (vers[google_type]);
 185 }
 186
 187 gboolean google_coord_to_mapcoord ( const VikCoord *src, gdouble xzoom, gdouble yzoom, MapCoord *dest )
 188 {
 189   g_assert ( src->mode == VIK_COORD_LATLON );
 190
 191   if ( xzoom != yzoom )
 192     return FALSE;
 193
 194   dest->scale = google_zoom ( xzoom );
 195   if ( dest->scale == 255 )
 196     return FALSE;
 197
 198   dest->x = (src->east_west + 180) / 360 * GZ(17) / xzoom;
 199   dest->y = (180 - MERCLAT(src->north_south)) / 360 * GZ(17) / xzoom;
 200   dest->z = 0;
 201
 202   return TRUE;
 203 }
 204
 205 void google_mapcoord_to_center_coord ( MapCoord *src, VikCoord *dest )
 206 {
 207   gdouble socalled_mpp = GZ(src->scale);
 208   dest->mode = VIK_COORD_LATLON;
 209   dest->east_west = ((src->x+0.5) / GZ(17) * socalled_mpp * 360) - 180;
 210   dest->north_south = DEMERCLAT(180 - ((src->y+0.5) / GZ(17) * socalled_mpp * 360));
 211 }
 212
 213 static int real_google_download ( MapCoord *src, const gchar *dest_fn, const char *verstr )
 214 {
 215    int res;
 216    gchar *uri = g_strdup_printf ( "/mt?n=404&v=%s&x=%d&y=%d&zoom=%d", verstr, src->x, src->y, src->scale );
 217    res = a_http_download_get_url ( "mt.google.com", uri, dest_fn, &google_options );
 218    g_free ( uri );
 219    return res;
 220 }
 221
 222 static int google_download ( MapCoord *src, const gchar *dest_fn )
 223 {
 224    const gchar *vers_str = google_version_number(src, TYPE_GOOGLE_MAPS);
 225    return(real_google_download ( src, dest_fn, vers_str ));
 226 }
 227
 228 static int google_trans_download ( MapCoord *src, const gchar *dest_fn )
 229 {
 230    const gchar *vers_str = google_version_number(src, TYPE_GOOGLE_TRANS);
 231    return(real_google_download ( src, dest_fn, vers_str ));
 232 }
 233
 234 static char *kh_encode(guint32 x, guint32 y, guint8 scale)
 235 {
 236   gchar *buf = g_malloc ( (20-scale)*sizeof(gchar) );
 237   guint32 ya = 1 << (17 - scale);
 238   gint8 i, j;
 239
 240   if (y < 0 || (ya-1 < y)) {
 241     strcpy(buf,"tqq"); /* BAD */
 242     return buf;
 243   }
 244   if (x < 0 || ya-1 < x) {
 245     x %= ya;
 246     if (x < 0)
 247       x += ya;
 248   }
 249
 250   buf[0] = 't';
 251   for (j = 1, i = 16; i >= scale; i--, j++) {
 252     ya /= 2;
 253     if (y < ya) {
 254       if (x < ya)
 255         buf[j]='q';
 256       else {
 257         buf[j]='r';
 258         x -= ya;
 259       }
 260     } else {
 261       if (x < ya) {
 262         buf[j] = 't';
 263         y -= ya;
 264       } else {
 265         buf[j] = 's';
 266         x -= ya;
 267         y -= ya;
 268       }
 269     }
 270   }
 271   buf[j] = '\0';
 272   return buf;
 273 }
 274
 275 static int google_kh_download ( MapCoord *src, const gchar *dest_fn )
 276 {
 277    int res;
 278    gchar *khenc = kh_encode( src->x, src->y, src->scale );
 279    const gchar *vers_str = google_version_number(src, TYPE_GOOGLE_SAT);
 280    gchar *uri = g_strdup_printf ( "/kh?n=404&v=%s&t=%s", vers_str, khenc );
 281    g_free ( khenc );
 282    res = a_http_download_get_url ( "kh.google.com", uri, dest_fn, &google_options );
 283    g_free ( uri );
 284    return(res);
 285 }