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