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