2 * viking -- GPS Data and Topo Analyzer, Explorer, and Manager
4 * Copyright (C) 2003-2005, Evan Battaglia <gtoevan@gmx.net>
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.
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.
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
30 #include "vikmapslayer.h"
35 VikMapsLayer_MapType map_type
= { 8, 256, 256, VIK_VIEWPORT_DRAWMODE_KH
, khmaps_coord_to_mapcoord
, khmaps_mapcoord_to_center_coord
, khmaps_download
};
37 maps_layer_register_type("Old KH Satellite Images", 8, &map_type
);
40 /* 1 << (x-1) is like a 2**(x-1) */
41 #define KH(x) ((1<<(x-1)))
43 static const gdouble scale_mpps
[] = { KH(1), KH(2), KH(3), KH(4), KH(5), KH(6), KH(7), KH(8), KH(9),
44 KH(10), KH(11), KH(12), KH(13), KH(14), KH(15) };
46 static const gint num_scales
= (sizeof(scale_mpps
) / sizeof(scale_mpps
[0])) - 1;
48 #define ERROR_MARGIN 0.01
49 guint8
khmaps_zoom ( gdouble mpp
) {
51 for ( i
= 0; i
<= num_scales
; i
++ ) {
52 if ( ABS(scale_mpps
[i
] - mpp
) < ERROR_MARGIN
)
58 gboolean
khmaps_coord_to_mapcoord ( const VikCoord
*src
, gdouble xzoom
, gdouble yzoom
, MapCoord
*dest
)
60 g_assert ( src
->mode
== VIK_COORD_LATLON
);
65 dest
->scale
= khmaps_zoom ( xzoom
);
66 if ( dest
->scale
== 255 )
70 dest
->x
= (gint
) ((1<<(16 - dest
->scale
)) + (src
->east_west
/ 180 * (1<<(16-dest
->scale
))));
71 dest
->y
= (gint
) ((1<<(16 - dest
->scale
)) - (src
->north_south
/ 180 * (1<<(16-dest
->scale
))));
76 void khmaps_mapcoord_to_center_coord ( MapCoord
*src
, VikCoord
*dest
)
78 dest
->mode
= VIK_COORD_LATLON
;
79 dest
->east_west
= (src
->x
+0.5 - (1<<(16 - src
->scale
))) * 180 / (1<<(16-src
->scale
));
80 dest
->north_south
= -((src
->y
+0.5 - (1<<(16 - src
->scale
))) * 180 / (1<<(16-src
->scale
)));
83 static char *kh_encode(guint32 x
, guint32 y
, guint8 scale
)
85 gchar
*buf
= g_malloc ( (20-scale
)*sizeof(gchar
) );
86 guint32 ya
= 1 << (17 - scale
);
89 if (y
< 0 || (ya
-1 < y
)) {
90 strcpy(buf
,"tqq"); /* BAD */
93 if (x
< 0 || ya
-1 < x
) {
100 for (j
= 1, i
= 16; i
>= scale
; i
--, j
++) {
124 void khmaps_download ( MapCoord
*src
, const gchar
*dest_fn
)
126 gchar
*tmp
= kh_encode(src
->x
, src
->y
, src
->scale
);
127 gchar
*uri
= g_strdup_printf ( "/kh?v=2&t=%s", tmp
);
128 g_print("%d %d %d = %s\n", src
->x
, src
->y
, src
->scale
, uri
);
130 a_http_download_get_url ( "kh.google.com", uri
, dest_fn
);
134 /* Popularity has its disadvantages ... */