src/osm.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 "viking.h"
  25 #include "coords.h"
  26 #include "vikcoord.h"
  27 #include "mapcoord.h"
  28 #include "vikmapslayer.h"
  29
  30 #include "osm.h"
  31
  32 static guint8 osm_zoom ( gdouble mpp );
  33
  34 static gboolean osm_coord_to_mapcoord ( const VikCoord *src, gdouble xzoom, gdouble yzoom, MapCoord *dest );
  35 static void osm_mapcoord_to_center_coord ( MapCoord *src, VikCoord *dest );
  36 static int osm_maplint_download ( MapCoord *src, const gchar *dest_fn );
  37 static int osm_mapnik_download ( MapCoord *src, const gchar *dest_fn );
  38 static int osm_osmarender_download ( MapCoord *src, const gchar *dest_fn );
  39 static int bluemarble_download ( MapCoord *src, const gchar *dest_fn );
  40
  41 static DownloadOptions osm_options = { NULL, 0 };
  42
  43 /* initialisation */
  44 void osm_init () {
  45   VikMapsLayer_MapType osmarender_type = { 12, 256, 256, VIK_VIEWPORT_DRAWMODE_MERCATOR, osm_coord_to_mapcoord, osm_mapcoord_to_center_coord, osm_osmarender_download };
  46   VikMapsLayer_MapType mapnik_type = { 13, 256, 256, VIK_VIEWPORT_DRAWMODE_MERCATOR, osm_coord_to_mapcoord, osm_mapcoord_to_center_coord, osm_mapnik_download };  VikMapsLayer_MapType maplint_type = { 14, 256, 256, VIK_VIEWPORT_DRAWMODE_MERCATOR, osm_coord_to_mapcoord, osm_mapcoord_to_center_coord, osm_maplint_download };
  47
  48   VikMapsLayer_MapType bluemarble_type = { 15, 256, 256, VIK_VIEWPORT_DRAWMODE_MERCATOR, osm_coord_to_mapcoord, osm_mapcoord_to_center_coord, bluemarble_download };
  49
  50   maps_layer_register_type("OpenStreetMap (Osmarender)", 12, &osmarender_type);
  51   maps_layer_register_type("OpenStreetMap (Mapnik)", 13, &mapnik_type);
  52   maps_layer_register_type("OpenStreetMap (Maplint)", 14, &maplint_type);
  53
  54   maps_layer_register_type("BlueMarble", 15, &bluemarble_type);
  55 }
  56
  57 /* 1 << (x) is like a 2**(x) */
  58 #define GZ(x) (1<<(x))
  59
  60 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),
  61                                       GZ(10), GZ(11), GZ(12), GZ(13), GZ(14), GZ(15), GZ(16), GZ(17) };
  62
  63 static const gint num_scales = (sizeof(scale_mpps) / sizeof(scale_mpps[0]));
  64
  65 #define ERROR_MARGIN 0.01
  66 guint8 osm_zoom ( gdouble mpp ) {
  67   gint i;
  68   for ( i = 0; i < num_scales; i++ ) {
  69     if ( ABS(scale_mpps[i] - mpp) < ERROR_MARGIN )
  70       return i;
  71   }
  72   return 255;
  73 }
  74
  75 static gboolean osm_coord_to_mapcoord ( const VikCoord *src, gdouble xzoom, gdouble yzoom, MapCoord *dest )
  76 {
  77   g_assert ( src->mode == VIK_COORD_LATLON );
  78
  79   if ( xzoom != yzoom )
  80     return FALSE;
  81
  82   dest->scale = osm_zoom ( xzoom );
  83   if ( dest->scale == 255 )
  84     return FALSE;
  85
  86   dest->x = (src->east_west + 180) / 360 * GZ(17) / xzoom;
  87   dest->y = (180 - MERCLAT(src->north_south)) / 360 * GZ(17) / xzoom;
  88   dest->z = 0;
  89   return TRUE;
  90 }
  91
  92 static void osm_mapcoord_to_center_coord ( MapCoord *src, VikCoord *dest )
  93 {
  94   gdouble socalled_mpp = GZ(src->scale);
  95   dest->mode = VIK_COORD_LATLON;
  96   dest->east_west = ((src->x+0.5) / GZ(17) * socalled_mpp * 360) - 180;
  97   dest->north_south = DEMERCLAT(180 - ((src->y+0.5) / GZ(17) * socalled_mpp * 360));
  98 }
  99
 100 /* Maplint tiles
 101  * Ex: http://dev.openstreetmap.org/~ojw/Tiles/maplint.php/10/517/375.png
 102  */
 103 static int osm_maplint_download ( MapCoord *src, const gchar *dest_fn )
 104 {
 105    int res = -1;
 106    gchar *uri = g_strdup_printf ( "/Tiles/maplint.php/%d/%d/%d.png", 17-src->scale, src->x, src->y );
 107    res = a_http_download_get_url ( "tah.openstreetmap.org", uri, dest_fn, &osm_options );
 108    g_free ( uri );
 109    return res;
 110 }
 111
 112 static int osm_mapnik_download ( MapCoord *src, const gchar *dest_fn )
 113 {
 114    int res = -1;
 115    gchar *uri = g_strdup_printf ( "/%d/%d/%d.png", 17-src->scale, src->x, src->y );
 116    res = a_http_download_get_url ( "tile.openstreetmap.org", uri, dest_fn, &osm_options );
 117    g_free ( uri );
 118    return res;
 119 }
 120
 121 static int osm_osmarender_download ( MapCoord *src, const gchar *dest_fn )
 122 {
 123    int res = -1;
 124    gchar *uri = g_strdup_printf ( "/Tiles/tile/%d/%d/%d.png", 17-src->scale, src->x, src->y );
 125    res = a_http_download_get_url ( "tah.openstreetmap.org", uri, dest_fn, &osm_options );
 126    g_free ( uri );
 127    return res;
 128 }
 129
 130 static int bluemarble_download ( MapCoord *src, const gchar *dest_fn )
 131 {
 132    int res = -1;
 133    gchar *uri = g_strdup_printf ( "/com.modestmaps.bluemarble/%d-r%d-c%d.jpg", 17-src->scale, src->y, src->x );
 134    res = a_http_download_get_url ( "s3.amazonaws.com", uri, dest_fn, &osm_options );
 135
 136    g_free ( uri );
 137    return res;
 138 }