2 * Hic Est Leo -- Global Map Operations
4 * (c) 2014 Martin Mares <mj@ucw.cz>
10 #include <ucw/mempool.h>
11 #include <ucw/simple-lists.h>
23 double map_min_x, map_min_y;
24 double map_max_x, map_max_y;
25 double page_width, page_height;
26 uns map_clip, map_rotate, map_draw_border;
29 char *map_style_sheet;
33 static struct cf_section map_style_cf = {
34 #define P(x) PTR_TO(struct data_source_style, x)
35 CF_TYPE(struct data_source_style),
37 CF_STRING("Name", P(name)),
43 static const char * const map_formats[] = {
49 static struct cf_section map_source_cf = {
50 #define P(x) PTR_TO(struct data_source, x)
51 CF_TYPE(struct data_source),
53 CF_STRING("File", P(file)),
54 CF_LOOKUP("Format", P(format), map_formats),
55 CF_LIST("StyleSheet", P(styles), &map_style_cf),
61 static struct cf_section map_cf = {
63 CF_DOUBLE("MinX", &map_min_x),
64 CF_DOUBLE("MinY", &map_min_y),
65 CF_DOUBLE("MaxX", &map_max_x),
66 CF_DOUBLE("MaxY", &map_max_y),
67 CF_DOUBLE("PageWidth", &page_width),
68 CF_DOUBLE("PageHeight", &page_height),
69 CF_UNS("Clip", &map_clip),
70 CF_UNS("Rotate", &map_rotate),
71 CF_UNS("DrawBorder", &map_draw_border),
72 CF_LIST("Source", &map_sources, &map_source_cf),
73 CF_STRING("Projection", &map_projection),
74 CF_STRING("SVGOutput", &map_svg_output),
79 static void CONSTRUCTOR map_preinit(void)
81 cf_declare_section("Map", &map_cf, 0);
86 double page_offset_x, page_offset_y;
87 double page_map_width, page_map_height;
89 void map_set_scale(void)
91 double x_range = map_max_x - map_min_x;
92 double y_range = map_max_y - map_min_y;
93 double x_scale = page_width / x_range;
94 double y_scale = page_height / y_range;
95 map_scale = MIN(x_scale, y_scale);
96 page_map_width = x_range * map_scale;
97 page_map_height = y_range * map_scale;
98 page_offset_x = (page_width - page_map_width) / 2;
99 page_offset_y = (page_height - page_map_height) / 2;
101 msg(L_INFO, "Setting scale %.3g (orig window [%.6g,%.6g], page window [%.6g,%.6g]+[%.6g,%.6g] on [%.6g,%.6g])",
104 page_map_width, page_map_height,
105 page_offset_x, page_offset_y,
106 page_width, page_height);
108 double pmin_x = INFINITY, pmax_x = -INFINITY;
109 double pmin_y = INFINITY, pmax_y = -INFINITY;
110 double rmin_x = INFINITY, rmax_x = -INFINITY;
111 double rmin_y = INFINITY, rmax_y = -INFINITY;
112 CLIST_FOR_EACH(struct data_source *, ds, map_sources)
114 CLIST_FOR_EACH(struct osm_node *, n, ds->osm->obj_list[OSM_TYPE_NODE])
116 pmin_x = MIN(pmin_x, n->x);
117 pmax_x = MAX(pmax_x, n->x);
118 pmin_y = MIN(pmin_y, n->y);
119 pmax_y = MAX(pmax_y, n->y);
120 n->x = (n->x - map_min_x) * map_scale + page_offset_x;
121 n->y = page_height - page_offset_y - (n->y - map_min_y) * map_scale;
122 rmin_x = MIN(rmin_x, n->x);
123 rmax_x = MAX(rmax_x, n->x);
124 rmin_y = MIN(rmin_y, n->y);
125 rmax_y = MAX(rmax_y, n->y);
128 msg(L_INFO, "Bounds before scaling: [%.10g,%.10g] x [%.10g,%.10g]", pmin_x, pmax_x, pmin_y, pmax_y);
129 msg(L_INFO, "Bounds after scaling: [%.10g,%.10g] x [%.10g,%.10g]", rmin_x, rmax_x, rmin_y, rmax_y);
131 if (debug_dump_after_scaling)
133 puts("=== Map after scaling ===");
134 CLIST_FOR_EACH(struct data_source *, ds, map_sources)
142 bool map_object_visible_p(struct osm_object *o)
150 struct osm_node *n = (struct osm_node *) o;
151 return (n->x >= page_offset_x - margin && n->x <= page_offset_x + page_map_width + margin &&
152 n->y >= page_offset_y - margin && n->y <= page_offset_y + page_map_height + margin);
156 struct osm_way *w = (struct osm_way *) o;
158 OSM_FOR_EACH_BEGIN(struct osm_object *, n, w->nodes)
160 ok |= map_object_visible_p(n);
165 case OSM_TYPE_RELATION:
167 struct osm_relation *r = (struct osm_relation *) o;
169 OSM_FOR_EACH_BEGIN(struct osm_object *, n, r->members)
171 ok |= map_object_visible_p(n);
176 case OSM_TYPE_MULTIPOLYGON:
177 return map_object_visible_p(&((struct osm_multipolygon *) o)->rel->o);
183 void map_load_styles(void)
186 clist_init(&style_cache);
188 CLIST_FOR_EACH(struct data_source *, ds, map_sources)
190 CLIST_FOR_EACH(struct data_source_style *, ss, ds->styles)
192 CLIST_FOR_EACH(struct simp_node *, n, style_cache)
194 struct data_source_style *x = n->p;
195 if (!strcmp(x->name, ss->name))
203 msg(L_DEBUG, "Loading style sheet %s", ss->name);
204 ss->css = css_load(ss->name);
208 printf("=== Stylesheet %s ===", ss->name);
212 simp_append(cf_get_pool(), &style_cache)->p = ss;
218 static void map_load_source(struct data_source *ds)
220 ds->osm = osm_init();
227 case DATA_SOURCE_OSMXML:
228 msg(L_INFO, "Parsing %s as OSM XML", ds->file);
230 die("OSM XML data sources must have a file name");
231 osm_xml_parse(ds->file);
235 case DATA_SOURCE_FIXED:
236 msg(L_INFO, "Adding fixed objects");
240 die("Invalid data source format");
243 if (debug_dump_source)
245 puts("=== Source data ===");
249 osm_make_multipolygons();
253 msg(L_INFO, "Projecting");
254 osm_project(map_projection);
255 if (debug_dump_after_proj)
257 puts("=== Map after projection ===");
263 void map_load_sources(void)
265 CLIST_FOR_EACH(struct data_source *, ds, map_sources)
269 void map_apply_styles(struct svg *svg)
271 struct style_results r;
274 CLIST_FOR_EACH(struct data_source *, ds, map_sources)
276 msg(L_INFO, "Applying stylesheet on %s", ds->file);
277 for (uns i = OSM_TYPE_NODE; i <= OSM_TYPE_MULTIPOLYGON; i++)
278 CLIST_FOR_EACH(struct osm_object *, o, ds->osm->obj_list[i])
280 if (debug_dump_styling)
282 puts("===============================");
285 if (!map_object_visible_p(o))
287 if (debug_dump_styling)
288 printf("--> invisible\n");
292 CLIST_FOR_EACH(struct data_source_style *, ss, ds->styles)
293 css_apply(ss->css, &r);
294 if (debug_dump_styling)
296 sym_from_style(o, &r, svg);
301 // FIXME: Ought to destroy the style_results