[andy/viking.git] / src / osm.c

/*
 * viking -- GPS Data and Topo Analyzer, Explorer, and Manager
 *
 * Copyright (C) 2003-2005, Evan Battaglia <gtoevan@gmx.net>
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 *
 */

#include <gtk/gtk.h>
#include <math.h>
#include "viking.h"
#include "coords.h"
#include "vikcoord.h"
#include "mapcoord.h"
#include "vikmapslayer.h"

#include "osm.h"

static guint8 osm_zoom ( gdouble mpp );

static gboolean osm_coord_to_mapcoord ( const VikCoord *src, gdouble xzoom, gdouble yzoom, MapCoord *dest );
static void osm_mapcoord_to_center_coord ( MapCoord *src, VikCoord *dest );
static int osm_maplint_download ( MapCoord *src, const gchar *dest_fn );
static int osm_mapnik_download ( MapCoord *src, const gchar *dest_fn );
static int osm_osmarender_download ( MapCoord *src, const gchar *dest_fn );
static int osm_cycle_download ( MapCoord *src, const gchar *dest_fn );
static int bluemarble_download ( MapCoord *src, const gchar *dest_fn );
static int openaerialmap_download ( MapCoord *src, const gchar *dest_fn );

static DownloadOptions osm_options = { NULL, 0, a_check_map_file };

/* initialisation */
void osm_init () {
  VikMapsLayer_MapType osmarender_type = { 12, 256, 256, VIK_VIEWPORT_DRAWMODE_MERCATOR, osm_coord_to_mapcoord, osm_mapcoord_to_center_coord, osm_osmarender_download };
  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 };
  VikMapsLayer_MapType cycle_type = { 17, 256, 256, VIK_VIEWPORT_DRAWMODE_MERCATOR, osm_coord_to_mapcoord, osm_mapcoord_to_center_coord, osm_cycle_download };

  VikMapsLayer_MapType bluemarble_type = { 15, 256, 256, VIK_VIEWPORT_DRAWMODE_MERCATOR, osm_coord_to_mapcoord, osm_mapcoord_to_center_coord, bluemarble_download };

  VikMapsLayer_MapType openaerialmap_type = { 20, 256, 256, VIK_VIEWPORT_DRAWMODE_MERCATOR, osm_coord_to_mapcoord, osm_mapcoord_to_center_coord, openaerialmap_download };

  maps_layer_register_type("OpenStreetMap (Osmarender)", 12, &osmarender_type);
  maps_layer_register_type("OpenStreetMap (Mapnik)", 13, &mapnik_type);
  maps_layer_register_type("OpenStreetMap (Maplint)", 14, &maplint_type);
  maps_layer_register_type("OpenStreetMap (Cycle)", 17, &cycle_type);

  maps_layer_register_type("BlueMarble", 15, &bluemarble_type);
  maps_layer_register_type("OpenAerialMap", 20, &openaerialmap_type);
}

/* 1 << (x) is like a 2**(x) */
#define GZ(x) (1<<(x))

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),
                                      GZ(10), GZ(11), GZ(12), GZ(13), GZ(14), GZ(15), GZ(16), GZ(17) };

static const gint num_scales = (sizeof(scale_mpps) / sizeof(scale_mpps[0]));

#define ERROR_MARGIN 0.01
guint8 osm_zoom ( gdouble mpp ) {
  gint i;
  for ( i = 0; i < num_scales; i++ ) {
    if ( ABS(scale_mpps[i] - mpp) < ERROR_MARGIN )
      return i;
  }
  return 255;
}

static gboolean osm_coord_to_mapcoord ( const VikCoord *src, gdouble xzoom, gdouble yzoom, MapCoord *dest )
{
  g_assert ( src->mode == VIK_COORD_LATLON );

  if ( xzoom != yzoom )
    return FALSE;

  dest->scale = osm_zoom ( xzoom );
  if ( dest->scale == 255 )
    return FALSE;

  dest->x = (src->east_west + 180) / 360 * GZ(17) / xzoom;
  dest->y = (180 - MERCLAT(src->north_south)) / 360 * GZ(17) / xzoom;
  dest->z = 0;
  return TRUE;
}

static void osm_mapcoord_to_center_coord ( MapCoord *src, VikCoord *dest )
{
  gdouble socalled_mpp = GZ(src->scale);
  dest->mode = VIK_COORD_LATLON;
  dest->east_west = ((src->x+0.5) / GZ(17) * socalled_mpp * 360) - 180;
  dest->north_south = DEMERCLAT(180 - ((src->y+0.5) / GZ(17) * socalled_mpp * 360));
}

/* Maplint tiles
 * Ex: http://dev.openstreetmap.org/~ojw/Tiles/maplint.php/10/517/375.png
 */
static int osm_maplint_download ( MapCoord *src, const gchar *dest_fn )
{
   int res = -1;
   gchar *uri = g_strdup_printf ( "/Tiles/maplint.php/%d/%d/%d.png", 17-src->scale, src->x, src->y );
   res = a_http_download_get_url ( "tah.openstreetmap.org", uri, dest_fn, &osm_options );
   g_free ( uri );
   return res;
}

static int osm_mapnik_download ( MapCoord *src, const gchar *dest_fn )
{
   int res = -1;
   gchar *uri = g_strdup_printf ( "/%d/%d/%d.png", 17-src->scale, src->x, src->y );
   res = a_http_download_get_url ( "tile.openstreetmap.org", uri, dest_fn, &osm_options );
   g_free ( uri );
   return res;
}

static int osm_osmarender_download ( MapCoord *src, const gchar *dest_fn )
{
   int res = -1;
   gchar *uri = g_strdup_printf ( "/Tiles/tile/%d/%d/%d.png", 17-src->scale, src->x, src->y );
   res = a_http_download_get_url ( "tah.openstreetmap.org", uri, dest_fn, &osm_options );
   g_free ( uri );
   return res;
}

static int bluemarble_download ( MapCoord *src, const gchar *dest_fn )
{
   int res = -1;
   gchar *uri = g_strdup_printf ( "/com.modestmaps.bluemarble/%d-r%d-c%d.jpg", 17-src->scale, src->y, src->x );
   res = a_http_download_get_url ( "s3.amazonaws.com", uri, dest_fn, &osm_options );

   g_free ( uri );
   return res;

}
   
static int openaerialmap_download ( MapCoord *src, const gchar *dest_fn )
{
   int res = -1;
   gchar *uri = g_strdup_printf ( "/tiles/1.0.0/openaerialmap-900913/%d/%d/%d.jpg", 17-src->scale, src->x, src->y );
   res = a_http_download_get_url ( "tile.openaerialmap.org", uri, dest_fn, &osm_options );

   g_free ( uri );
   return res;
}

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