]>
Commit | Line | Data |
---|---|---|
50a14534 EB |
1 | /* |
2 | * viking -- GPS Data and Topo Analyzer, Explorer, and Manager | |
3 | * | |
4 | * Copyright (C) 2003-2005, Evan Battaglia <gtoevan@gmx.net> | |
5 | * | |
6 | * Lat/Lon plotting functions calcxy* are from GPSDrive | |
7 | * GPSDrive Copyright (C) 2001-2004 Fritz Ganter <ganter@ganter.at> | |
8 | * | |
9 | * Multiple UTM zone patch by Kit Transue <notlostyet@didactek.com> | |
10 | * | |
11 | * This program is free software; you can redistribute it and/or modify | |
12 | * it under the terms of the GNU General Public License as published by | |
13 | * the Free Software Foundation; either version 2 of the License, or | |
14 | * (at your option) any later version. | |
15 | * | |
16 | * This program is distributed in the hope that it will be useful, | |
17 | * but WITHOUT ANY WARRANTY; without even the implied warranty of | |
18 | * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the | |
19 | * GNU General Public License for more details. | |
20 | * | |
21 | * You should have received a copy of the GNU General Public License | |
22 | * along with this program; if not, write to the Free Software | |
23 | * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA | |
24 | * | |
25 | */ | |
26 | ||
27 | #define DEFAULT_BACKGROUND_COLOR "#CCCCCC" | |
28 | ||
29 | #include <gtk/gtk.h> | |
30 | #include <math.h> | |
31 | ||
32 | #include "coords.h" | |
33 | #include "vikcoord.h" | |
34 | #include "vikviewport.h" | |
35 | ||
36 | #include "mapcoord.h" | |
37 | ||
38 | /* for ALTI_TO_MPP */ | |
39 | #include "globals.h" | |
40 | #include "googlemaps.h" | |
41 | #include "khmaps.h" | |
42 | ||
43 | static gdouble EASTING_OFFSET = 500000.0; | |
44 | ||
45 | static void viewport_class_init ( VikViewportClass *klass ); | |
46 | static void viewport_init ( VikViewport *vvp ); | |
47 | static void viewport_finalize ( GObject *gob ); | |
48 | static void viewport_utm_zone_check ( VikViewport *vvp ); | |
49 | ||
50 | static gboolean calcxy(double *x, double *y, double lg, double lt, double zero_long, double zero_lat, double pixelfact_x, double pixelfact_y, gint mapSizeX2, gint mapSizeY2 ); | |
51 | static gboolean calcxy_rev(double *lg, double *lt, gint x, gint y, double zero_long, double zero_lat, double pixelfact_x, double pixelfact_y, gint mapSizeX2, gint mapSizeY2 ); | |
52 | double calcR (double lat); | |
53 | ||
54 | static double Radius[181]; | |
55 | static void viewport_init_ra(); | |
56 | ||
57 | static GObjectClass *parent_class; | |
58 | ||
59 | static void viewport_google_rezoom ( VikViewport *vvp ); | |
60 | ||
61 | ||
62 | struct _VikViewport { | |
63 | GtkDrawingArea drawing_area; | |
64 | GdkPixmap *scr_buffer; | |
65 | gint width, height; | |
66 | VikCoord center; | |
67 | VikCoordMode coord_mode; | |
68 | gdouble xmpp, ympp; | |
69 | ||
70 | GdkPixbuf *alpha_pixbuf; | |
71 | guint8 alpha_pixbuf_width; | |
72 | guint8 alpha_pixbuf_height; | |
73 | ||
74 | gdouble utm_zone_width; | |
75 | gboolean one_utm_zone; | |
76 | ||
77 | GdkGC *background_gc; | |
78 | GdkColor background_color; | |
79 | ||
80 | /* subset of coord types. lat lon can be plotted in 2 ways, google or exp. */ | |
81 | VikViewportDrawMode drawmode; | |
82 | ||
83 | /* handy conversion factors which make google plotting extremely fast */ | |
84 | gdouble google_calcx_fact; | |
85 | gdouble google_calcy_fact; | |
86 | gdouble google_calcx_rev_fact; | |
87 | gdouble google_calcy_rev_fact; | |
88 | }; | |
89 | ||
90 | static gdouble | |
91 | viewport_utm_zone_width ( VikViewport *vvp ) | |
92 | { | |
93 | if ( vvp->coord_mode == VIK_COORD_UTM ) { | |
94 | struct LatLon ll; | |
95 | ||
96 | /* get latitude of screen bottom */ | |
97 | struct UTM utm = *((struct UTM *)(vik_viewport_get_center ( vvp ))); | |
98 | utm.northing -= vvp -> height * vvp -> ympp / 2; | |
99 | a_coords_utm_to_latlon ( &utm, &ll ); | |
100 | ||
101 | /* boundary */ | |
102 | ll.lon = (utm.zone - 1) * 6 - 180 ; | |
103 | a_coords_latlon_to_utm ( &ll, &utm); | |
104 | return fabs ( utm.easting - EASTING_OFFSET ) * 2; | |
105 | } else | |
106 | return 0.0; | |
107 | } | |
108 | ||
109 | ||
110 | GType vik_viewport_get_type (void) | |
111 | { | |
112 | static GType vvp_type = 0; | |
113 | ||
114 | if (!vvp_type) | |
115 | { | |
116 | static const GTypeInfo vvp_info = | |
117 | { | |
118 | sizeof (VikViewportClass), | |
119 | NULL, /* base_init */ | |
120 | NULL, /* base_finalize */ | |
121 | (GClassInitFunc) viewport_class_init, | |
122 | NULL, /* class_finalize */ | |
123 | NULL, /* class_data */ | |
124 | sizeof (VikViewport), | |
125 | 0, | |
126 | (GInstanceInitFunc) viewport_init, | |
127 | }; | |
128 | vvp_type = g_type_register_static ( GTK_TYPE_DRAWING_AREA, "VikViewport", &vvp_info, 0 ); | |
129 | } | |
130 | return vvp_type; | |
131 | } | |
132 | ||
133 | static void viewport_class_init ( VikViewportClass *klass ) | |
134 | { | |
135 | /* Destructor */ | |
136 | GObjectClass *object_class; | |
137 | ||
138 | object_class = G_OBJECT_CLASS (klass); | |
139 | ||
140 | object_class->finalize = viewport_finalize; | |
141 | ||
142 | parent_class = g_type_class_peek_parent (klass); | |
143 | } | |
144 | ||
145 | VikViewport *vik_viewport_new () | |
146 | { | |
147 | return VIK_VIEWPORT ( g_object_new ( VIK_VIEWPORT_TYPE, NULL ) ); | |
148 | } | |
149 | ||
150 | static void viewport_init ( VikViewport *vvp ) | |
151 | { | |
152 | viewport_init_ra(); | |
153 | ||
154 | /* TODO: not static */ | |
155 | vvp->xmpp = 4.0; | |
156 | vvp->ympp = 4.0; | |
157 | vvp->coord_mode = VIK_COORD_UTM; | |
158 | vvp->drawmode = VIK_VIEWPORT_DRAWMODE_UTM; | |
159 | vvp->center.north_south = 0; | |
160 | vvp->center.east_west = -166021; | |
161 | vvp->center.utm_zone = 31; | |
162 | vvp->center.utm_letter = 'N'; | |
163 | vvp->scr_buffer = NULL; | |
164 | vvp->alpha_pixbuf = NULL; | |
165 | vvp->alpha_pixbuf_width = vvp->alpha_pixbuf_height = 0; | |
166 | vvp->utm_zone_width = 0.0; | |
167 | vvp->background_gc = NULL; | |
168 | g_signal_connect (G_OBJECT(vvp), "configure_event", G_CALLBACK(vik_viewport_configure), NULL); | |
169 | } | |
170 | ||
171 | GdkColor *vik_viewport_get_background_gdkcolor ( VikViewport *vvp ) | |
172 | { | |
173 | GdkColor *rv = g_malloc ( sizeof ( GdkColor ) ); | |
174 | *rv = vvp->background_color; | |
175 | return rv; | |
176 | } | |
177 | ||
178 | /* returns pointer to internal static storage, changes next time function called, use quickly */ | |
179 | const gchar *vik_viewport_get_background_color ( VikViewport *vvp ) | |
180 | { | |
181 | static gchar color[8]; | |
182 | g_snprintf(color, sizeof(color), "#%.2x%.2x%.2x", (int)(vvp->background_color.red/256),(int)(vvp->background_color.green/256),(int)(vvp->background_color.blue/256)); | |
183 | return color; | |
184 | } | |
185 | ||
186 | void vik_viewport_set_background_color ( VikViewport *vvp, const gchar *colorname ) | |
187 | { | |
188 | g_assert ( vvp->background_gc ); | |
189 | gdk_color_parse ( colorname, &(vvp->background_color) ); | |
190 | gdk_gc_set_rgb_fg_color ( vvp->background_gc, &(vvp->background_color) ); | |
191 | } | |
192 | ||
193 | void vik_viewport_set_background_gdkcolor ( VikViewport *vvp, GdkColor *color ) | |
194 | { | |
195 | g_assert ( vvp->background_gc ); | |
196 | vvp->background_color = *color; | |
197 | gdk_gc_set_rgb_fg_color ( vvp->background_gc, color ); | |
198 | } | |
199 | ||
200 | ||
201 | GdkGC *vik_viewport_new_gc ( VikViewport *vvp, const gchar *colorname, gint thickness ) | |
202 | { | |
203 | GdkGC *rv; | |
204 | GdkColor color; | |
205 | ||
206 | rv = gdk_gc_new ( GTK_WIDGET(vvp)->window ); | |
207 | gdk_color_parse ( colorname, &color ); | |
208 | gdk_gc_set_rgb_fg_color ( rv, &color ); | |
209 | gdk_gc_set_line_attributes ( rv, thickness, GDK_LINE_SOLID, GDK_CAP_ROUND, GDK_JOIN_ROUND ); | |
210 | return rv; | |
211 | } | |
212 | ||
213 | GdkGC *vik_viewport_new_gc_from_color ( VikViewport *vvp, GdkColor *color, gint thickness ) | |
214 | { | |
215 | GdkGC *rv; | |
216 | ||
217 | rv = gdk_gc_new ( GTK_WIDGET(vvp)->window ); | |
218 | gdk_gc_set_rgb_fg_color ( rv, color ); | |
219 | gdk_gc_set_line_attributes ( rv, thickness, GDK_LINE_SOLID, GDK_CAP_ROUND, GDK_JOIN_ROUND ); | |
220 | return rv; | |
221 | } | |
222 | ||
223 | void vik_viewport_configure_manually ( VikViewport *vvp, gint width, guint height ) | |
224 | { | |
225 | vvp->width = width; | |
226 | vvp->height = height; | |
227 | if ( vvp->scr_buffer ) | |
228 | g_object_unref ( G_OBJECT ( vvp->scr_buffer ) ); | |
229 | vvp->scr_buffer = gdk_pixmap_new ( GTK_WIDGET(vvp)->window, vvp->width, vvp->height, -1 ); | |
230 | } | |
231 | ||
232 | ||
233 | GdkPixmap *vik_viewport_get_pixmap ( VikViewport *vvp ) | |
234 | { | |
235 | return vvp->scr_buffer; | |
236 | } | |
237 | ||
238 | gboolean vik_viewport_configure ( VikViewport *vvp ) | |
239 | { | |
240 | g_return_val_if_fail ( vvp != NULL, TRUE ); | |
241 | ||
242 | vvp->width = GTK_WIDGET(vvp)->allocation.width; | |
243 | vvp->height = GTK_WIDGET(vvp)->allocation.height; | |
244 | ||
245 | if ( vvp->scr_buffer ) | |
246 | g_object_unref ( G_OBJECT ( vvp->scr_buffer ) ); | |
247 | ||
248 | vvp->scr_buffer = gdk_pixmap_new ( GTK_WIDGET(vvp)->window, vvp->width, vvp->height, -1 ); | |
249 | ||
250 | /* this is down here so it can get a GC (necessary?) */ | |
251 | if ( ! vvp->background_gc ) | |
252 | { | |
253 | vvp->background_gc = vik_viewport_new_gc ( vvp, "", 1 ); | |
254 | vik_viewport_set_background_color ( vvp, DEFAULT_BACKGROUND_COLOR ); /* set to "backup" color in vvp->background_color */ | |
255 | } | |
256 | ||
257 | return FALSE; | |
258 | } | |
259 | ||
260 | static void viewport_finalize ( GObject *gob ) | |
261 | { | |
262 | VikViewport *vvp = VIK_VIEWPORT(gob); | |
263 | ||
264 | g_return_if_fail ( vvp != NULL ); | |
265 | ||
266 | if ( vvp->scr_buffer ) | |
267 | g_object_unref ( G_OBJECT ( vvp->scr_buffer ) ); | |
268 | ||
269 | if ( vvp->alpha_pixbuf ) | |
270 | g_object_unref ( G_OBJECT ( vvp->alpha_pixbuf ) ); | |
271 | ||
272 | if ( vvp->background_gc ) | |
273 | g_object_unref ( G_OBJECT ( vvp->background_gc ) ); | |
274 | ||
275 | G_OBJECT_CLASS(parent_class)->finalize(gob); | |
276 | } | |
277 | ||
278 | void vik_viewport_clear ( VikViewport *vvp ) | |
279 | { | |
280 | g_return_if_fail ( vvp != NULL ); | |
281 | gdk_draw_rectangle(GDK_DRAWABLE(vvp->scr_buffer), vvp->background_gc, TRUE, 0, 0, vvp->width, vvp->height); | |
282 | } | |
283 | ||
284 | void vik_viewport_sync ( VikViewport *vvp ) | |
285 | { | |
286 | g_return_if_fail ( vvp != NULL ); | |
287 | gdk_draw_drawable(GTK_WIDGET(vvp)->window, GTK_WIDGET(vvp)->style->bg_gc[0], GDK_DRAWABLE(vvp->scr_buffer), 0, 0, 0, 0, vvp->width, vvp->height); | |
288 | } | |
289 | ||
290 | void vik_viewport_pan_sync ( VikViewport *vvp, gint x_off, gint y_off ) | |
291 | { | |
292 | g_return_if_fail ( vvp != NULL ); | |
293 | gdk_draw_drawable(GTK_WIDGET(vvp)->window, GTK_WIDGET(vvp)->style->bg_gc[0], GDK_DRAWABLE(vvp->scr_buffer), 0, 0, x_off, y_off, vvp->width, vvp->height); | |
294 | } | |
295 | ||
296 | void vik_viewport_set_zoom ( VikViewport *vvp, gdouble xympp ) | |
297 | { | |
298 | g_return_if_fail ( vvp != NULL ); | |
299 | if ( xympp >= VIK_VIEWPORT_MIN_ZOOM && xympp <= VIK_VIEWPORT_MAX_ZOOM ) | |
300 | vvp->xmpp = vvp->ympp = xympp; | |
301 | ||
302 | if ( vvp->drawmode == VIK_VIEWPORT_DRAWMODE_UTM ) | |
303 | viewport_utm_zone_check(vvp); | |
304 | else if ( vvp->drawmode == VIK_VIEWPORT_DRAWMODE_GOOGLE ) | |
305 | viewport_google_rezoom ( vvp ); | |
306 | } | |
307 | ||
308 | /* or could do factor */ | |
309 | void vik_viewport_zoom_in ( VikViewport *vvp ) | |
310 | { | |
311 | g_return_if_fail ( vvp != NULL ); | |
312 | if ( vvp->xmpp >= (VIK_VIEWPORT_MIN_ZOOM*2) && vvp->ympp >= (VIK_VIEWPORT_MIN_ZOOM*2) ) | |
313 | { | |
314 | vvp->xmpp /= 2; | |
315 | vvp->ympp /= 2; | |
316 | ||
317 | if ( vvp->drawmode == VIK_VIEWPORT_DRAWMODE_GOOGLE ) | |
318 | viewport_google_rezoom ( vvp ); | |
319 | ||
320 | viewport_utm_zone_check(vvp); | |
321 | } | |
322 | } | |
323 | ||
324 | void vik_viewport_zoom_out ( VikViewport *vvp ) | |
325 | { | |
326 | g_return_if_fail ( vvp != NULL ); | |
327 | if ( vvp->xmpp <= (VIK_VIEWPORT_MAX_ZOOM/2) && vvp->ympp <= (VIK_VIEWPORT_MAX_ZOOM/2) ) | |
328 | { | |
329 | vvp->xmpp *= 2; | |
330 | vvp->ympp *= 2; | |
331 | ||
332 | if ( vvp->drawmode == VIK_VIEWPORT_DRAWMODE_GOOGLE ) | |
333 | viewport_google_rezoom ( vvp ); | |
334 | ||
335 | viewport_utm_zone_check(vvp); | |
336 | } | |
337 | } | |
338 | ||
339 | gdouble vik_viewport_get_zoom ( VikViewport *vvp ) | |
340 | { | |
341 | if ( vvp->xmpp == vvp->ympp ) | |
342 | return vvp->xmpp; | |
343 | return 0.0; | |
344 | } | |
345 | ||
346 | gdouble vik_viewport_get_xmpp ( VikViewport *vvp ) | |
347 | { | |
348 | return vvp->xmpp; | |
349 | } | |
350 | ||
351 | gdouble vik_viewport_get_ympp ( VikViewport *vvp ) | |
352 | { | |
353 | return vvp->ympp; | |
354 | } | |
355 | ||
356 | void vik_viewport_set_xmpp ( VikViewport *vvp, gdouble xmpp ) | |
357 | { | |
358 | if ( xmpp >= VIK_VIEWPORT_MIN_ZOOM && xmpp <= VIK_VIEWPORT_MAX_ZOOM ) { | |
359 | vvp->xmpp = xmpp; | |
360 | if ( vvp->drawmode == VIK_VIEWPORT_DRAWMODE_UTM ) | |
361 | viewport_utm_zone_check(vvp); | |
362 | if ( vvp->drawmode == VIK_VIEWPORT_DRAWMODE_GOOGLE ) | |
363 | viewport_google_rezoom ( vvp ); | |
364 | } | |
365 | } | |
366 | ||
367 | void vik_viewport_set_ympp ( VikViewport *vvp, gdouble ympp ) | |
368 | { | |
369 | if ( ympp >= VIK_VIEWPORT_MIN_ZOOM && ympp <= VIK_VIEWPORT_MAX_ZOOM ) { | |
370 | vvp->ympp = ympp; | |
371 | if ( vvp->drawmode == VIK_VIEWPORT_DRAWMODE_UTM ) | |
372 | viewport_utm_zone_check(vvp); | |
373 | if ( vvp->drawmode == VIK_VIEWPORT_DRAWMODE_GOOGLE ) | |
374 | viewport_google_rezoom ( vvp ); | |
375 | } | |
376 | } | |
377 | ||
378 | ||
379 | const VikCoord *vik_viewport_get_center ( VikViewport *vvp ) | |
380 | { | |
381 | g_return_val_if_fail ( vvp != NULL, NULL ); | |
382 | return &(vvp->center); | |
383 | } | |
384 | ||
385 | /* called every time we update coordinates/zoom */ | |
386 | static void viewport_utm_zone_check ( VikViewport *vvp ) | |
387 | { | |
388 | if ( vvp->coord_mode == VIK_COORD_UTM ) | |
389 | { | |
390 | struct UTM utm; | |
391 | struct LatLon ll; | |
392 | a_coords_utm_to_latlon ( (struct UTM *) &(vvp->center), &ll ); | |
393 | a_coords_latlon_to_utm ( &ll, &utm ); | |
394 | if ( utm.zone != vvp->center.utm_zone ) | |
395 | *((struct UTM *)(&vvp->center)) = utm; | |
396 | ||
397 | /* misc. stuff so we don't have to check later */ | |
398 | vvp->utm_zone_width = viewport_utm_zone_width ( vvp ); | |
399 | vvp->one_utm_zone = ( vik_viewport_rightmost_zone(vvp) == vik_viewport_leftmost_zone(vvp) ); | |
400 | } | |
401 | } | |
402 | ||
403 | void vik_viewport_set_center_latlon ( VikViewport *vvp, const struct LatLon *ll ) | |
404 | { | |
405 | vik_coord_load_from_latlon ( &(vvp->center), vvp->coord_mode, ll ); | |
406 | if ( vvp->coord_mode == VIK_COORD_UTM ) | |
407 | viewport_utm_zone_check ( vvp ); | |
408 | } | |
409 | ||
410 | void vik_viewport_set_center_utm ( VikViewport *vvp, const struct UTM *utm ) | |
411 | { | |
412 | vik_coord_load_from_utm ( &(vvp->center), vvp->coord_mode, utm ); | |
413 | if ( vvp->coord_mode == VIK_COORD_UTM ) | |
414 | viewport_utm_zone_check ( vvp ); | |
415 | } | |
416 | ||
417 | void vik_viewport_set_center_coord ( VikViewport *vvp, const VikCoord *coord ) | |
418 | { | |
419 | vvp->center = *coord; | |
420 | if ( vvp->coord_mode == VIK_COORD_UTM ) | |
421 | viewport_utm_zone_check ( vvp ); | |
422 | } | |
423 | ||
424 | void vik_viewport_corners_for_zonen ( VikViewport *vvp, int zone, VikCoord *ul, VikCoord *br ) | |
425 | { | |
426 | g_return_if_fail ( vvp->coord_mode == VIK_COORD_UTM ); | |
427 | ||
428 | /* get center, then just offset */ | |
429 | vik_viewport_center_for_zonen ( vvp, VIK_UTM(ul), zone ); | |
430 | ul->mode = VIK_COORD_UTM; | |
431 | *br = *ul; | |
432 | ||
433 | ul->north_south += (vvp->ympp * vvp->height / 2); | |
434 | ul->east_west -= (vvp->xmpp * vvp->width / 2); | |
435 | br->north_south -= (vvp->ympp * vvp->height / 2); | |
436 | br->east_west += (vvp->xmpp * vvp->width / 2); | |
437 | } | |
438 | ||
439 | void vik_viewport_center_for_zonen ( VikViewport *vvp, struct UTM *center, int zone) | |
440 | { | |
441 | if ( vvp->coord_mode == VIK_COORD_UTM ) { | |
442 | *center = *((struct UTM *)(vik_viewport_get_center ( vvp ))); | |
443 | center->easting -= ( zone - center->zone ) * vvp->utm_zone_width; | |
444 | center->zone = zone; | |
445 | } | |
446 | } | |
447 | ||
448 | gchar vik_viewport_leftmost_zone ( VikViewport *vvp ) | |
449 | { | |
450 | if ( vvp->coord_mode == VIK_COORD_UTM ) { | |
451 | VikCoord coord; | |
452 | g_assert ( vvp != NULL ); | |
453 | vik_viewport_screen_to_coord ( vvp, 0, 0, &coord ); | |
454 | return coord.utm_zone; | |
455 | } | |
456 | return '\0'; | |
457 | } | |
458 | ||
459 | gchar vik_viewport_rightmost_zone ( VikViewport *vvp ) | |
460 | { | |
461 | if ( vvp->coord_mode == VIK_COORD_UTM ) { | |
462 | VikCoord coord; | |
463 | g_assert ( vvp != NULL ); | |
464 | vik_viewport_screen_to_coord ( vvp, vvp->width, 0, &coord ); | |
465 | return coord.utm_zone; | |
466 | } | |
467 | return '\0'; | |
468 | } | |
469 | ||
470 | ||
471 | void vik_viewport_set_center_screen ( VikViewport *vvp, int x, int y ) | |
472 | { | |
473 | g_return_if_fail ( vvp != NULL ); | |
474 | if ( vvp->coord_mode == VIK_COORD_UTM ) { | |
475 | /* slightly optimized */ | |
476 | vvp->center.east_west += vvp->xmpp * (x - (vvp->width/2)); | |
477 | vvp->center.north_south += vvp->ympp * ((vvp->height/2) - y); | |
478 | viewport_utm_zone_check ( vvp ); | |
479 | } else { | |
480 | VikCoord tmp; | |
481 | vik_viewport_screen_to_coord ( vvp, x, y, &tmp ); | |
482 | vik_viewport_set_center_coord ( vvp, &tmp ); | |
483 | } | |
484 | } | |
485 | ||
486 | gint vik_viewport_get_width( VikViewport *vvp ) | |
487 | { | |
488 | g_return_val_if_fail ( vvp != NULL, 0 ); | |
489 | return vvp->width; | |
490 | } | |
491 | ||
492 | gint vik_viewport_get_height( VikViewport *vvp ) | |
493 | { | |
494 | g_return_val_if_fail ( vvp != NULL, 0 ); | |
495 | return vvp->height; | |
496 | } | |
497 | ||
498 | void vik_viewport_screen_to_coord ( VikViewport *vvp, int x, int y, VikCoord *coord ) | |
499 | { | |
500 | if ( vvp->coord_mode == VIK_COORD_UTM ) { | |
501 | int zone_delta; | |
502 | struct UTM *utm = (struct UTM *) coord; | |
503 | coord->mode = VIK_COORD_UTM; | |
504 | ||
505 | g_return_if_fail ( vvp != NULL ); | |
506 | ||
507 | utm->zone = vvp->center.utm_zone; | |
508 | utm->letter = vvp->center.utm_letter; | |
509 | utm->easting = ( ( x - ( vvp->width / 2) ) * vvp->xmpp ) + vvp->center.east_west; | |
510 | zone_delta = floor( (utm->easting - EASTING_OFFSET ) / vvp->utm_zone_width + 0.5 ); | |
511 | utm->zone += zone_delta; | |
512 | utm->easting -= zone_delta * vvp->utm_zone_width; | |
513 | utm->northing = ( ( ( vvp->height / 2) - y ) * vvp->ympp ) + vvp->center.north_south; | |
514 | } else if ( vvp->coord_mode == VIK_COORD_LATLON ) { | |
515 | coord->mode = VIK_COORD_LATLON; | |
516 | if ( vvp->drawmode == VIK_VIEWPORT_DRAWMODE_EXPEDIA ) | |
517 | calcxy_rev(&(coord->east_west), &(coord->north_south), x, y, vvp->center.east_west, vvp->center.north_south, vvp->xmpp * ALTI_TO_MPP, vvp->ympp * ALTI_TO_MPP, vvp->width/2, vvp->height/2); | |
518 | else if ( vvp->drawmode == VIK_VIEWPORT_DRAWMODE_GOOGLE ) { | |
519 | /* google */ | |
520 | coord->east_west = (x - (vvp->width/2)) * vvp->google_calcx_rev_fact + vvp->center.east_west; | |
521 | coord->north_south = ((vvp->height/2) - y) * vvp->google_calcy_rev_fact + vvp->center.north_south; | |
522 | } else if ( vvp->drawmode == VIK_VIEWPORT_DRAWMODE_KH ) { | |
523 | coord->east_west = vvp->center.east_west + (180.0 * vvp->xmpp / 65536 / 256 * (x - vvp->width/2)); | |
524 | coord->north_south = vvp->center.north_south + (180.0 * vvp->ympp / 65536 / 256 * (vvp->height/2 - y)); | |
525 | } else if ( vvp->drawmode == VIK_VIEWPORT_DRAWMODE_MERCATOR ) { | |
526 | /* FIXMERCATOR */ | |
527 | coord->east_west = vvp->center.east_west + (180.0 * vvp->xmpp / 65536 / 256 * (x - vvp->width/2)); | |
528 | coord->north_south = DEMERCLAT ( MERCLAT(vvp->center.north_south) + (180.0 * vvp->ympp / 65536 / 256 * (vvp->height/2 - y)) ); | |
529 | ||
530 | #if 0 | |
531 | --> THIS IS JUNK HERE. | |
532 | *y = vvp->height/2 + (65536.0 / 180 / vvp->ympp * (MERCLAT(center->lat) - MERCLAT(ll->lat)))*256.0; | |
533 | ||
534 | (*y - vvp->height/2) / 256 / 65536 * 180 * vvp->ympp = (MERCLAT(center->lat) - MERCLAT(ll->lat); | |
535 | DML((180.0 * vvp->ympp / 65536 / 256 * (vvp->height/2 - y)) + ML(cl)) = ll | |
536 | #endif | |
537 | } | |
538 | } | |
539 | } | |
540 | ||
541 | void vik_viewport_coord_to_screen ( VikViewport *vvp, const VikCoord *coord, int *x, int *y ) | |
542 | { | |
543 | static VikCoord tmp; | |
544 | g_return_if_fail ( vvp != NULL ); | |
545 | ||
546 | if ( coord->mode != vvp->coord_mode ) | |
547 | { | |
548 | g_warning ( "Have to convert in vik_viewport_coord_to_screen! This should never happen!"); | |
549 | vik_coord_copy_convert ( coord, vvp->coord_mode, &tmp ); | |
550 | coord = &tmp; | |
551 | } | |
552 | ||
553 | if ( vvp->coord_mode == VIK_COORD_UTM ) { | |
554 | struct UTM *center = (struct UTM *) &(vvp->center); | |
555 | struct UTM *utm = (struct UTM *) coord; | |
556 | if ( center->zone != utm->zone && vvp->one_utm_zone ) | |
557 | { | |
558 | *x = *y = VIK_VIEWPORT_UTM_WRONG_ZONE; | |
559 | return; | |
560 | } | |
561 | ||
562 | *x = ( (utm->easting - center->easting) / vvp->xmpp ) + (vvp->width / 2) - | |
563 | (center->zone - utm->zone ) * vvp->utm_zone_width / vvp->xmpp; | |
564 | *y = (vvp->height / 2) - ( (utm->northing - center->northing) / vvp->ympp ); | |
565 | } else if ( vvp->coord_mode == VIK_COORD_LATLON ) { | |
566 | struct LatLon *center = (struct LatLon *) &(vvp->center); | |
567 | struct LatLon *ll = (struct LatLon *) coord; | |
568 | double xx,yy; | |
569 | if ( vvp->drawmode == VIK_VIEWPORT_DRAWMODE_EXPEDIA ) { | |
570 | calcxy ( &xx, &yy, center->lon, center->lat, ll->lon, ll->lat, vvp->xmpp * ALTI_TO_MPP, vvp->ympp * ALTI_TO_MPP, vvp->width / 2, vvp->height / 2 ); | |
571 | *x = xx; *y = yy; | |
572 | } else if ( vvp->drawmode == VIK_VIEWPORT_DRAWMODE_GOOGLE ) { | |
573 | /* google */ | |
574 | *x = vvp->google_calcx_fact * (ll->lon - center->lon) + (vvp->width/2); | |
575 | *y = vvp->google_calcy_fact * (center->lat - ll->lat) + (vvp->height/2); | |
576 | } else if ( vvp->drawmode == VIK_VIEWPORT_DRAWMODE_KH ) { | |
577 | /* subtract, convert to KH coords; blow it up by 256 */ | |
578 | *x = vvp->width/2 + (65536.0 / 180 / vvp->xmpp * (ll->lon - center->lon))*256.0; | |
579 | *y = vvp->height/2 + (65536.0 / 180 / vvp->ympp * (center->lat - ll->lat))*256.0; | |
580 | } else if ( vvp->drawmode == VIK_VIEWPORT_DRAWMODE_MERCATOR ) { | |
581 | /* FIXMERCATOR: Optimize */ | |
582 | *x = vvp->width/2 + (65536.0 / 180 / vvp->xmpp * (ll->lon - center->lon))*256.0; | |
583 | *y = vvp->height/2 + (65536.0 / 180 / vvp->ympp * (MERCLAT(center->lat) - MERCLAT(ll->lat)))*256.0; | |
584 | } | |
585 | } | |
586 | } | |
587 | ||
588 | void vik_viewport_draw_line ( VikViewport *vvp, GdkGC *gc, gint x1, gint y1, gint x2, gint y2 ) | |
589 | { | |
590 | if ( ! ( ( x1 < 0 && x2 < 0 ) || ( y1 < 0 && y2 < 0 ) || | |
591 | ( x1 > vvp->width && x2 > vvp->width ) || ( y1 > vvp->height && y2 > vvp->height ) ) ) | |
592 | gdk_draw_line ( vvp->scr_buffer, gc, x1, y1, x2, y2); | |
593 | } | |
594 | ||
595 | void vik_viewport_draw_rectangle ( VikViewport *vvp, GdkGC *gc, gboolean filled, gint x1, gint y1, gint x2, gint y2 ) | |
596 | { | |
597 | if ( x1 > -10 && x1 < vvp->width + 10 && y1 > -10 && y1 < vvp->height + 10 ) | |
598 | gdk_draw_rectangle ( vvp->scr_buffer, gc, filled, x1, y1, x2, y2); | |
599 | } | |
600 | ||
601 | void vik_viewport_draw_string ( VikViewport *vvp, GdkFont *font, GdkGC *gc, gint x1, gint y1, const gchar *string ) | |
602 | { | |
603 | if ( x1 > -100 && x1 < vvp->width + 100 && y1 > -100 && y1 < vvp->height + 100 ) | |
604 | gdk_draw_string ( vvp->scr_buffer, font, gc, x1, y1, string ); | |
605 | } | |
606 | ||
607 | /* shouldn't use this -- slow -- change the alpha channel instead. */ | |
608 | void vik_viewport_draw_pixbuf_with_alpha ( VikViewport *vvp, GdkPixbuf *pixbuf, gint alpha, | |
609 | gint src_x, gint src_y, gint dest_x, gint dest_y, gint w, gint h ) | |
610 | { | |
611 | gint real_dest_x = MAX(dest_x,0); | |
612 | gint real_dest_y = MAX(dest_y,0); | |
613 | ||
614 | if ( alpha == 0 ) | |
615 | return; /* don't waste your time */ | |
616 | ||
617 | if ( w > vvp->alpha_pixbuf_width || h > vvp->alpha_pixbuf_height ) | |
618 | { | |
619 | if ( vvp->alpha_pixbuf ) | |
620 | g_object_unref ( G_OBJECT ( vvp->alpha_pixbuf ) ); | |
621 | vvp->alpha_pixbuf_width = MAX(w,vvp->alpha_pixbuf_width); | |
622 | vvp->alpha_pixbuf_height = MAX(h,vvp->alpha_pixbuf_height); | |
623 | vvp->alpha_pixbuf = gdk_pixbuf_new ( GDK_COLORSPACE_RGB, FALSE, 8, vvp->alpha_pixbuf_width, vvp->alpha_pixbuf_height ); | |
624 | } | |
625 | ||
626 | w = MIN(w,vvp->width - dest_x); | |
627 | h = MIN(h,vvp->height - dest_y); | |
628 | ||
629 | /* check that we are drawing within boundaries. */ | |
630 | src_x += (real_dest_x - dest_x); | |
631 | src_y += (real_dest_y - dest_y); | |
632 | w -= (real_dest_x - dest_x); | |
633 | h -= (real_dest_y - dest_y); | |
634 | ||
635 | gdk_pixbuf_get_from_drawable ( vvp->alpha_pixbuf, vvp->scr_buffer, NULL, | |
636 | real_dest_x, real_dest_y, 0, 0, w, h ); | |
637 | ||
638 | /* do a composite */ | |
639 | gdk_pixbuf_composite ( pixbuf, vvp->alpha_pixbuf, 0, 0, w, h, -src_x, -src_y, 1, 1, 0, alpha ); | |
640 | ||
641 | /* draw pixbuf_tmp */ | |
642 | vik_viewport_draw_pixbuf ( vvp, vvp->alpha_pixbuf, 0, 0, real_dest_x, real_dest_y, w, h ); | |
643 | } | |
644 | ||
645 | void vik_viewport_draw_pixbuf ( VikViewport *vvp, GdkPixbuf *pixbuf, gint src_x, gint src_y, | |
646 | gint dest_x, gint dest_y, gint w, gint h ) | |
647 | { | |
648 | gdk_draw_pixbuf ( vvp->scr_buffer, | |
649 | // GTK_WIDGET(vvp)->style->black_gc, | |
650 | NULL, | |
651 | pixbuf, | |
652 | src_x, src_y, dest_x, dest_y, w, h, | |
653 | GDK_RGB_DITHER_NONE, 0, 0 ); | |
654 | } | |
655 | ||
656 | void vik_viewport_draw_arc ( VikViewport *vvp, GdkGC *gc, gboolean filled, gint x, gint y, gint width, gint height, gint angle1, gint angle2 ) | |
657 | { | |
658 | gdk_draw_arc ( vvp->scr_buffer, gc, filled, x, y, width, height, angle1, angle2 ); | |
659 | } | |
660 | ||
661 | ||
662 | void vik_viewport_draw_polygon ( VikViewport *vvp, GdkGC *gc, gboolean filled, GdkPoint *points, gint npoints ) | |
663 | { | |
664 | gdk_draw_polygon ( vvp->scr_buffer, gc, filled, points, npoints ); | |
665 | } | |
666 | ||
667 | VikCoordMode vik_viewport_get_coord_mode ( const VikViewport *vvp ) | |
668 | { | |
669 | g_assert ( vvp ); | |
670 | return vvp->coord_mode; | |
671 | } | |
672 | ||
673 | static void viewport_set_coord_mode ( VikViewport *vvp, VikCoordMode mode ) | |
674 | { | |
675 | g_return_if_fail ( vvp != NULL ); | |
676 | vvp->coord_mode = mode; | |
677 | vik_coord_convert ( &(vvp->center), mode ); | |
678 | } | |
679 | ||
680 | /* Thanks GPSDrive */ | |
681 | static gboolean calcxy_rev(double *lg, double *lt, gint x, gint y, double zero_long, double zero_lat, double pixelfact_x, double pixelfact_y, gint mapSizeX2, gint mapSizeY2 ) | |
682 | { | |
683 | int px, py; | |
684 | gdouble dif, lat, lon; | |
685 | double Ra = Radius[90+(gint)zero_lat]; | |
686 | ||
687 | px = (mapSizeX2 - x) * pixelfact_x; | |
688 | py = (-mapSizeY2 + y) * pixelfact_y; | |
689 | ||
690 | lat = zero_lat - py / Ra; | |
691 | lat = zero_lat - py / Ra; | |
692 | lon = | |
693 | zero_long - | |
694 | px / (Ra * | |
695 | cos (lat * DEG2RAD)); | |
696 | ||
697 | dif = lat * (1 - (cos ((fabs (lon - zero_long)) * DEG2RAD))); | |
698 | lat = lat - dif / 1.5; | |
699 | lon = | |
700 | zero_long - | |
701 | px / (Ra * | |
702 | cos (lat * DEG2RAD)); | |
703 | ||
704 | *lt = lat; | |
705 | *lg = lon; | |
706 | return (TRUE); | |
707 | } | |
708 | ||
709 | /* Thanks GPSDrive */ | |
710 | static gboolean calcxy(double *x, double *y, double lg, double lt, double zero_long, double zero_lat, double pixelfact_x, double pixelfact_y, gint mapSizeX2, gint mapSizeY2 ) | |
711 | { | |
712 | double dif; | |
713 | double Ra; | |
714 | gint mapSizeX = 2 * mapSizeX2; | |
715 | gint mapSizeY = 2 * mapSizeY2; | |
716 | ||
717 | g_assert ( lt >= -90.0 && lt <= 90.0 ); | |
718 | // lg *= rad2deg; // FIXME, optimize equations | |
719 | // lt *= rad2deg; | |
720 | Ra = Radius[90+(gint)lt]; | |
721 | *x = Ra * | |
722 | cos (lt*DEG2RAD) * (lg - zero_long); | |
723 | *y = Ra * (lt - zero_lat); | |
724 | dif = Ra * RAD2DEG * (1 - (cos ((DEG2RAD * (lg - zero_long))))); | |
725 | *y = *y + dif / 1.85; | |
726 | *x = *x / pixelfact_x; | |
727 | *y = *y / pixelfact_y; | |
728 | *x = mapSizeX2 - *x; | |
729 | *y += mapSizeY2; | |
730 | if ((*x < 0)||(*x >= mapSizeX)||(*y < 0)||(*y >= mapSizeY)) | |
731 | return (FALSE); | |
732 | return (TRUE); | |
733 | } | |
734 | ||
735 | static void viewport_init_ra() | |
736 | { | |
737 | static gboolean done_before = FALSE; | |
738 | if ( !done_before ) | |
739 | { | |
740 | gint i; | |
741 | for ( i = -90; i <= 90; i++) | |
742 | Radius[i+90] = calcR ( (double)i ) * DEG2RAD; | |
743 | done_before = TRUE; | |
744 | } | |
745 | } | |
746 | ||
747 | double calcR (double lat) | |
748 | { | |
749 | double a = 6378.137, r, sc, x, y, z; | |
750 | double e2 = 0.081082 * 0.081082; | |
751 | /* | |
752 | * the radius of curvature of an ellipsoidal Earth in the plane of the | |
753 | * meridian is given by | |
754 | * | |
755 | * R' = a * (1 - e^2) / (1 - e^2 * (sin(lat))^2)^(3/2) | |
756 | * | |
757 | * | |
758 | * where a is the equatorial radius, b is the polar radius, and e is | |
759 | * the eccentricity of the ellipsoid = sqrt(1 - b^2/a^2) | |
760 | * | |
761 | * a = 6378 km (3963 mi) Equatorial radius (surface to center distance) | |
762 | * b = 6356.752 km (3950 mi) Polar radius (surface to center distance) e | |
763 | * = 0.081082 Eccentricity | |
764 | */ | |
765 | ||
766 | lat = lat * DEG2RAD; | |
767 | sc = sin (lat); | |
768 | x = a * (1.0 - e2); | |
769 | z = 1.0 - e2 * sc * sc; | |
770 | y = pow (z, 1.5); | |
771 | r = x / y; | |
772 | r = r * 1000.0; | |
773 | return r; | |
774 | } | |
775 | ||
776 | gboolean vik_viewport_is_one_zone ( VikViewport *vvp ) | |
777 | { | |
778 | return vvp->coord_mode == VIK_COORD_UTM && vvp->one_utm_zone; | |
779 | } | |
780 | ||
781 | void vik_viewport_draw_layout ( VikViewport *vvp, GdkGC *gc, gint x, gint y, PangoLayout *layout ) | |
782 | { | |
783 | if ( x > -100 && x < vvp->width + 100 && y > -100 && y < vvp->height + 100 ) | |
784 | gdk_draw_layout ( vvp->scr_buffer, gc, x, y, layout ); | |
785 | } | |
786 | ||
787 | void vik_gc_get_fg_color ( GdkGC *gc, GdkColor *dest ) | |
788 | { | |
789 | static GdkGCValues values; | |
790 | gdk_gc_get_values ( gc, &values ); | |
791 | gdk_colormap_query_color ( gdk_colormap_get_system(), values.foreground.pixel, dest ); | |
792 | } | |
793 | ||
794 | GdkFunction vik_gc_get_function ( GdkGC *gc ) | |
795 | { | |
796 | static GdkGCValues values; | |
797 | gdk_gc_get_values ( gc, &values ); | |
798 | return values.function; | |
799 | } | |
800 | ||
801 | void vik_viewport_set_drawmode ( VikViewport *vvp, VikViewportDrawMode drawmode ) | |
802 | { | |
803 | vvp->drawmode = drawmode; | |
804 | if ( drawmode == VIK_VIEWPORT_DRAWMODE_UTM ) | |
805 | viewport_set_coord_mode ( vvp, VIK_COORD_UTM ); | |
806 | else { | |
807 | viewport_set_coord_mode ( vvp, VIK_COORD_LATLON ); | |
808 | if ( drawmode == VIK_VIEWPORT_DRAWMODE_GOOGLE ) | |
809 | viewport_google_rezoom ( vvp ); | |
810 | } | |
811 | } | |
812 | ||
813 | VikViewportDrawMode vik_viewport_get_drawmode ( VikViewport *vvp ) | |
814 | { | |
815 | return vvp->drawmode; | |
816 | } | |
817 | ||
818 | static void viewport_google_rezoom ( VikViewport *vvp ) | |
819 | { | |
820 | vvp->google_calcx_fact = (GOOGLEMAPS_ZOOM_ONE_MPP * 65536.0 * 0.7716245833877 / vvp->xmpp); | |
821 | vvp->google_calcy_fact = (GOOGLEMAPS_ZOOM_ONE_MPP * 65536.0 / vvp->ympp); | |
822 | vvp->google_calcx_rev_fact = 1 / vvp->google_calcx_fact; | |
823 | vvp->google_calcy_rev_fact = 1 / vvp->google_calcy_fact; | |
824 | } |