-
Notifications
You must be signed in to change notification settings - Fork 46
/
debug_openstreetmap.m
33 lines (29 loc) · 1.4 KB
/
debug_openstreetmap.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
% example script for using the OpenStreetMap functions
% to illustrate usage of the OpenStreetMap functions for MATLAB
%
% See also PARSE_OPENSTREETMAP, PLOT_WAY, EXTRACT_CONNECTIVITY,
% GET_UNIQUE_NODE_XY.
%
% 2010.11.25 (c) Ioannis Filippidis, [email protected]
% download an OpenStreetMap XML Data file (extension .osm) from the
% OpenStreetMap website:
% http://www.openstreetmap.org/
% after zooming in the area of interest and using the "Export" option to
% save it as an OpenStreetMap XML Data file, selecting this from the
% "Format to Export" options. The OSM XML is specified in:
% http://wiki.openstreetmap.org/wiki/.osm
openstreetmap_filename = 'openstreetmap/genoa_small.osm';
%% convert XML -> MATLAB struct
% convert the OpenStreetMap XML Data file donwloaded as map.osm
% to a MATLAB structure containing part of the information describing the
% transportation network
[parsed_osm, osm_xml] = parse_openstreetmap(openstreetmap_filename);
%% plot
% plot the network, optionally a raster image can also be provided for the
% map under the vector graphics of the network
ax = gca;
plot_way(ax, parsed_osm)
%plot_way(ax, parsed_osm, map_img_filename) % if you also have a raster image
%% find connectivity
[connectivity_matrix, intersection_node_indices] = extract_connectivity(parsed_osm);
[uniquend] = get_unique_node_xy(parsed_osm, intersection_node_indices);