{"id":11017,"date":"2023-07-06T09:00:00","date_gmt":"2023-07-06T07:00:00","guid":{"rendered":"https:\/\/www.sigterritoires.fr\/?p=11017"},"modified":"2023-07-05T15:48:31","modified_gmt":"2023-07-05T13:48:31","slug":"netwotk-analysis-with-qgis","status":"publish","type":"post","link":"https:\/\/www.sigterritoires.fr\/index.php\/en\/netwotk-analysis-with-qgis\/","title":{"rendered":"Netwotk Analysis with QGis"},"content":{"rendered":"\n

We’ll look at how to download and prepare OpenStreetMap data for a network analysis, then at all the QGis processes grouped together under the heading Network Analysis.<\/p>\n\n\n\n\n\n\n\n

We’ll look here at examples of road network analysis, but it’s clear that the same possibilities apply to any type of network: pipes, hydrographic, etc.<\/p>\n\n\n\n

One of the most common uses of network analysis is to find the shortest (or fastest) path. This requires a topologically flawless network, i.e. one in which all sections are well connected to each other, and in which the digitization of the section has been carried out in the direction of travel (from the starting point to the end point of the section).<\/p>\n\n\n\n

The hardest part of the job is bringing the network up to standard.<\/p>\n\n\n\n

As far as the road network is concerned, we have data in OpenStreetMap that meets all these conditions.<\/p>\n\n\n\n

For a given area, we have three sources of data, all available on the Geofabrik<\/a> website: shapefiles, compressed .osm files and pbf files.<\/p>\n\n\n\n

\"network<\/a><\/figure>\n\n\n\n

Preparing OSM data<\/h2>\n\n\n\n

Depending on the type of file you have chosen and your software configuration, the data preparation stage before you can use the network analysis tools can be more or less time-consuming and complicated.<\/p>\n\n\n\n

If you only have QGis<\/h3>\n\n\n\n

If you don’t have a PostgreSQL \/ Postgis<\/a> database, the format to download is the shapefile<\/a> format. You won’t get all the information contained in the other formats, but for road network (or hydrographic) analysis you’ll have everything you need.<\/p>\n\n\n

\n
\"layers<\/a><\/figure>\n<\/div>\n\n\n

To work with QGis network analysis processes, we’ll be using the gis_osm_roads_free_1<\/em><\/strong> layer.<\/p>\n\n\n\n

If you have QGis and Postgis<\/h3>\n\n\n\n

Depending on the tools you use to load an OSM file into Postgis, you may choose either OSM or pbf format.<\/p>\n\n\n\n

osm2pgrouting <\/em><\/strong>is a command-line utility for loading .osm format files into a PostgreSQL \/ Postgis database.<\/p>\n\n\n\n

If you have osm2pgrouting<\/strong><\/h3>\n\n\n\n

As its name suggests, it prepares OSM data by default for use with Postgis routing tools (pgrouting).<\/p>\n\n\n\n

An example command line is shown below:<\/p>\n\n\n\n

osm2pgrouting.exe \u2013clean=1 -f fichier_telecharg\u00e9.osm -d base_de_donn\u00e9es -p 5433 -h localhost --schema uruguay -U postgres -W ******  -c mapconfig_for_cars.xml\n<\/code><\/pre>\n\n\n\n
-d indicates which database on your server the data should be loaded into<\/p>\n\n\n\n
-p corresponds to the PostgreSQL server listening port<\/p>\n\n\n\n
-W is the user password<\/p>\n\n\n\n
-c the configuration file to be used for loading. osm2pgrouting offers three files, adapted for cars, pedestrians and bicycles.<\/p>\n\n\n\n
\"ligne<\/a><\/figure>\n\n\n\n
The loading result is as follows:<\/p>\n\n\n
\n
\"tables<\/a><\/figure>\n<\/div>\n\n\n
If you have osm2pgsql<\/h3>\n\n\n\n
osm2pgsql <\/em><\/strong>is a command-line utility for loading all OSM data from a .pbf file into a PostgreSQL \/ Postgis database.<\/p>\n\n\n\n
An example command line is<\/p>\n\n\n\n
C:\\osm2pgsql\\osm2pgsql-bin\\osm2pgsql.exe -c -d postgis_33_sample -U postgres -W -H localhost -P 5433 -S default.style uruguay-latest.osm.pbf<\/em><\/code><\/pre>\n\n\n\n
The parameters are similar to those used by osm2pgrouting<\/em><\/strong>, except for the -S parameter<\/em><\/strong> which defines the style file to be used. You’ll find this default style file in the osm2pgsql-bin<\/em><\/strong> directory, in the osm2pgsq<\/em><\/strong>l installation directory. It indicates which columns are to be created in the postgis table.<\/p>\n\n\n\n
For our network analysis work, it’s likely that at some point we’ll use the maximum authorized speed of the sections to calculate the fastest route. Unfortunately, this parameter (maxspeed) is not included in the list of attributes to be created. For it to be created, a line must be added to the default.style file:<\/p>\n\n\n\n
node,way maxspeed text linear<\/em><\/code><\/pre>\n\n\n\n
\"add<\/a><\/figure>\n\n\n\n
UOnce the command line has been executed:<\/p>\n\n\n\n
\"run<\/a><\/figure>\n\n\n\n
You’ll have four new tables in your Postgis database:<\/p>\n\n\n
\n
\"\npostgis<\/a><\/figure>\n<\/div>\n\n\n
Loading in QGis<\/h2>\n\n\n\n
For shp files<\/em><\/strong>, load the gis_osm_roads_free_1<\/em><\/strong> layer.<\/p>\n\n\n\n
\"load<\/a><\/figure>\n\n\n\n
For osm <\/em><\/strong>files processed with osm2pgrouting<\/em><\/strong>, load the ways <\/em><\/strong>table.<\/p>\n\n\n\n
\"load<\/a><\/figure>\n\n\n\n
For pbf files<\/strong> processed with osm2pgsql<\/em><\/strong>, load the planet_osm_line <\/em><\/strong>table. Contrary to popular belief, the planet_osm_roads<\/em><\/strong> table contains only the main roads. In the planet_osm_line<\/em><\/strong> table, you’ll have all linear entities, not just the road network. To have only the roads, you need to apply a filter on the highway <\/em><\/strong>attribute to retain only those rows where this attribute is not NULL. If you want to work on the hydrographic network, the IS NOT NULL filter must be applied to the waterway <\/em><\/strong>attribute, and if you want to work on the railway network, the railway <\/em><\/strong>attribute must be filtered.<\/p>\n\n\n\n
\"loading<\/a><\/figure>\n\n\n\n
If you’re working on the scale of a small town or village, you can use a specific symbology.<\/p>\n\n\n\n
\"network<\/a><\/figure>\n\n\n\n
You can download these symbols in xml format from this link<\/a>. To use them, open the style manager, click on import\/export and point to the downloaded oneway.xml<\/em><\/strong> file.<\/p>\n\n\n\n
\"symboles<\/a><\/figure>\n\n\n\n
To apply them, open the layer properties gis_osm_roadfs_free_1 -> symbology<\/em><\/strong>. Select Categorized, value-> oneway.<\/em><\/strong><\/p>\n\n\n\n
You’ll get three classes with values B, F and T.<\/p>\n\n\n\n
Apply the symbol<\/p>\n\n\n\n
    \n
  • oneway_both to value B,<\/li>\n\n\n\n
  • oneway_forward to value F,<\/li>\n\n\n\n
  • oneway_opposite to value T.<\/li>\n<\/ul>\n\n\n\n
    You also have network_speed <\/em><\/strong>tags, which display the maximum speed of the section and must use the layer’s maxspeed <\/em><\/strong>field.<\/p>\n","protected":false},"excerpt":{"rendered":"
    We’ll look at how to download and prepare OpenStreetMap data for a network analysis, then at all the QGis processes grouped together under the heading Network Analysis.<\/p>\n","protected":false},"author":1,"featured_media":10969,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"give_campaign_id":0,"_bbp_topic_count":0,"_bbp_reply_count":0,"_bbp_total_topic_count":0,"_bbp_total_reply_count":0,"_bbp_voice_count":0,"_bbp_anonymous_reply_count":0,"_bbp_topic_count_hidden":0,"_bbp_reply_count_hidden":0,"_bbp_forum_subforum_count":0,"sfsi_plus_gutenberg_text_before_share":"","sfsi_plus_gutenberg_show_text_before_share":"","sfsi_plus_gutenberg_icon_type":"","sfsi_plus_gutenberg_icon_alignemt":"","sfsi_plus_gutenburg_max_per_row":"","_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0,"jetpack_post_was_ever_published":false,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"footnotes":""},"categories":[1377],"tags":[2329,2327,2331,2333,2335,2337,2339,1349,1386,2341,2343],"class_list":["post-11017","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-qgis-en","tag-analysis-en","tag-download","tag-network-en","tag-ogr2pgsql-en","tag-osm-en","tag-osm2pgrouting-en","tag-pbf-en","tag-postgis-en","tag-qgis-en","tag-shp-en","tag-symbologie-en"],"aioseo_notices":[],"jetpack_featured_media_url":"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2023\/06\/symbologie_osm_qgis.png?fit=2643%2C1467&ssl=1","jetpack_shortlink":"https:\/\/wp.me\/p6XU0A-2RH","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/11017","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/comments?post=11017"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/11017\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media\/10969"}],"wp:attachment":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media?parent=11017"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=11017"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=11017"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}