{"id":3308,"date":"2016-06-15T10:29:20","date_gmt":"2016-06-15T09:29:20","guid":{"rendered":"http:\/\/www.sigterritoires.fr\/?p=3308"},"modified":"2016-06-15T10:30:30","modified_gmt":"2016-06-15T09:30:30","slug":"developper-une-application-avec-pgrouting-sous-windows-4zones-de-chalandise-et-isochrones","status":"publish","type":"post","link":"https:\/\/www.sigterritoires.fr\/index.php\/developper-une-application-avec-pgrouting-sous-windows-4zones-de-chalandise-et-isochrones\/","title":{"rendered":"D\u00e9velopper une application avec pgrouting sous Windows (4):zones de chalandise et isochrones"},"content":{"rendered":"

pgRouting<\/a> comprend une fonction appel\u00e9e pgr_drivingDistance, qui, en conjonction avec trois autres fonctions, pgr_withPointsDD, pgr_alphaShape et pgr_pointsAsPolygon, permet de d\u00e9terminer quels chemins ont un certain co\u00fbt par rapport \u00e0 un point de d\u00e9part ou d’arriv\u00e9e.<\/p>\n

La zone qui entoure ces chemins est d\u00e9sign\u00e9e comme une zone de chalandise. pgr_drivingDistance, en conjonction avec pgr_pointsAsPolygon est souvent utilis\u00e9e pour d\u00e9terminer les exigences de couverture d’un service ou pour d\u00e9terminer les emplacements optimaux pour de nouvelles installations.<\/p>\n

Par exemple, si vous avez une caserne de pompiers, vous devez d\u00e9terminer toutes les parties de la ville qui peuvent \u00eatre atteintes dans les cinq minutes. Un probl\u00e8me inverse serait de trouver toutes les parties de la ville o\u00f9 les r\u00e9sidents peuvent conduire jusqu’\u00e0 un h\u00f4pital dans les cinq minutes.
\nLa fonction qui calcule les chemins pour les deux probl\u00e8mes est appel\u00e9e pgr_drivingDistance car elle est souvent utilis\u00e9e lorsque le co\u00fbt est mesur\u00e9 en temps de conduite. Cependant, les co\u00fbts ne sont pas forc\u00e9ment\u00a0 limit\u00e9s \u00e0 un temps ou une distance.
\nPour ce type d’application, conna\u00eetre tous les itin\u00e9raires vers ou \u00e0 partir d’un point n’est pas tr\u00e8s utile. C’est l\u00e0 que les fonctions\u00a0 pgr_alphaShape et pgr_pointsAsPolygon entrent en jeu. Le pgr_alphaShape prend en entr\u00e9e un ensemble de coordonn\u00e9es x, y et renvoie un ensemble de x, y qui correspondent aux sommets d’un polygone qui englobe tous les points en entr\u00e9e.
\npgr_pointsAsPolygon est similaire \u00e0 pgr_alphaShape, elle prend en entr\u00e9e un ensemble de coordonn\u00e9es x, y, mais au lieu de retourner une autre s\u00e9rie de points, renvoie un polygone\u00a0 qui entoure l’ensemble des points d’entr\u00e9e . En interne, il utilise la fonction pgr_alphaShape.<\/p>\n

Utilisation de pgr_drivingDistance<\/h2>\n
Penons le cas d’une caserne de pompiers et cherchons la zone couverte en 5 minutes. La commande sql correspondant \u00e0 ces crit\u00e8res se traduit par:<\/p>\n
SELECT * FROM pgr_drivingDistance(‘SELECT gid as id, source, target , cost_s as cost, reverse_cost_s as reverse_cost <\/strong><\/em>FROM public.ways’, 25736, 300, true)<\/strong><\/em><\/p>\n
Pour cette requ\u00eate nous utilisons dans la clause SELECT le co\u00fbt d’un tron\u00e7on en secondes (cost_s), et nous indiquons le reverse_cost_s car m\u00eame les camions de pompiers ne peuvent pas prendre les rues \u00e0 contresens. 25736<\/strong> <\/em>est le n\u0153ud correspondant \u00e0 la caserne de pompiers, 300<\/strong> <\/em>est le co\u00fbt recherch\u00e9 (5 minutes x 60 secondes), et la valeur true<\/strong><\/em> indique que le graphe est directionnel (prise en compte des sens interdits).<\/p>\n
$\"requ\u00eate$ <\/a><\/p>\n
Le r\u00e9sultat de la requ\u00eate retourne 6948 n\u0153uds qui sont dans un rayon de 5 minutes de la caserne de pompiers. Par contre le r\u00e9sultat de la requ\u00eate n’ayant pas de g\u00e9om\u00e9trie, nous ne pouvons pas les visualiser dans QGis.<\/p>\n
Pour pouvoir le faire il suffit de modifier la requ\u00eate comme suit:<\/p>\n
SELECT * FROM pgr_drivingDistance(‘SELECT gid as id, source, target , cost_s as cost, reverse_cost_s as reverse_cost <\/strong><\/em>FROM public.ways’, 25736, 300, true)<\/strong> as di JOIN ways_vertices_pgr on di.node=ways_vertices_pgr.id<\/strong><\/em><\/p>\n
Nous effectuons une jointure<\/a> avec la table des n\u0153uds (ways_vertices_pgr) et le r\u00e9sultat de la requ\u00eate contient alors les g\u00e9om\u00e9tries:<\/p>\n
$\"jointure$ <\/a>Nous pouvons alors visualiser le r\u00e9sultat dans QGis:<\/p>\n
$\"r\u00e9suktat$ <\/a><\/p>\n

Cr\u00e9ation de l’aire de chalandise avec pgr_pointsAsPolygon<\/h2>\n
Pour cr\u00e9er l’aire de chalandise, nous utilisons la requ\u00eate pr\u00e9c\u00e9dente comme sql de la fonction pgr_pointsAsPolygon, en cr\u00e9ant un polygone en sortie que nous pouvons charger dans QGis:
\nSELECT 1 As id, ST_SetSRID(pgr_pointsAsPolygon(<\/strong><\/em>
\n $$SELECT di.seq AS id, ST_X(v.the_geom) AS x, ST_Y(v.the_geom) As y<\/strong><\/em>
\n FROM pgr_drivingDistance(\u00a0\u00bbSELECT gid As id, source, target,<\/strong><\/em>
\n cost_s AS cost, reverse_cost_s AS reverse_cost<\/strong><\/em>
\n FROM public.ways\u00a0\u00bb,25736 , 300, true<\/strong><\/em>
\n ) AS di INNER JOIN public.ways_vertices_pgr AS v ON di.node = v.id$$<\/strong><\/em>
\n ), 4326) As geom;<\/strong><\/em><\/p>\n
$\"requ\u00eate$ <\/a>pgr_pointAsPolygon utilise les id,x et y des points g\u00e9n\u00e9r\u00e9s par pgr_drivingDistance. Pour les x et y on utilise les fonctions St_x et St_y sur la g\u00e9om\u00e9trie des points.<\/p>\n
Le r\u00e9sultat, une fois charg\u00e9 dans QGis est le suivant:<\/p>\n
$\"r\u00e9sultat$ <\/a>Construction d’isochrones<\/h2>\n
Les isochrones correspondent aux zones atteignables en une s\u00e9rie de pas de temps. On peut, dans cet exemple, vouloir conna\u00eetre quelles sont les zones qu’on peut atteindre en 1, 2, 3, 4 et 5 minutes.<\/p>\n
Pour cela nous allons proc\u00e9der en deux \u00e9tapes.<\/p>\n
La premi\u00e8re consiste \u00e0 cr\u00e9er une table avec nos noeuds r\u00e9sultant de pgr_drivingDistance pour le pas de temps le plus fort (5 minutes). Logiquement, dans cette table on trouvera aussi les n\u0153uds atteignables dans les laps de temps plus courts.<\/p>\n
Nous cr\u00e9ons cette table (dd_caserne) avec la commande sql suivante:<\/p>\n
CREATE TABLE dd_caserne AS<\/strong><\/em>
\n SELECT di.seq As id, di.node, di.edge,<\/strong><\/em>
\n di.agg_cost As temps_acces, <\/strong><\/em>v.the_geom As geom<\/strong><\/em>
\n FROM pgr_drivingDistance(‘SELECT gid As id, source, target,<\/strong><\/em>
\n cost_s AS cost, reverse_cost_s AS reverse_cost<\/strong><\/em>
\n FROM public.ways’,25736 , 300, true<\/strong><\/em>
\n ) AS di INNER JOIN public.ways_vertices_pgr AS v ON di.node = v.id;<\/strong><\/em><\/p>\n
La table r\u00e9sultante est la suivante:<\/p>\n
$\"table$ <\/a>A partir de cette table nous allons construire une s\u00e9rie de polygones avec pgr_pointsAsPolygon, pour les points situ\u00e9s \u00e0 60, 120,180,240 et 300 secondes de temps d’acc\u00e8s en utilisant generate_series<\/p>\n
SELECT i As temps_acces,ST_SetSRID(pgr_pointsAsPolygon(<\/strong><\/em>
\n‘SELECT id, ST_X(geom) AS x, ST_Y(geom) As y<\/strong><\/em>
\nFROM dd_caserne<\/strong><\/em>
\nWHERE temps_acces <= ‘ || i::text ) ,4326) As geom<\/strong><\/em>
\nFROM generate_series(60,300,60) As i<\/strong><\/em>
\nORDER BY i DESC<\/strong><\/em><\/p>\n
$\"cr\u00e9ation$ <\/a>Le r\u00e9sultat comprend les 5 polygones isochrones, que nous pouvons maintenant visualiser dans QGis:<\/p>\n
$\"isochrones$ <\/a><\/p>\n","protected":false},"excerpt":{"rendered":"
pgRouting comprend une fonction appel\u00e9e pgr_drivingDistance, qui, en conjonction avec trois autres fonctions, pgr_withPointsDD, pgr_alphaShape et pgr_pointsAsPolygon, permet de d\u00e9terminer quels chemins ont un certain co\u00fbt par rapport \u00e0 un point de d\u00e9part ou d’arriv\u00e9e. La…<\/p>\n","protected":false},"author":1,"featured_media":0,"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":[311,62],"tags":[333,326,332,334,58],"class_list":["post-3308","post","type-post","status-publish","format-standard","hentry","category-postgresql","category-qgis-2","tag-isochrones","tag-pgrouting","tag-pgr_drivingdistance","tag-pgr_pointsaspolygon","tag-qgis"],"aioseo_notices":[],"jetpack_featured_media_url":"","jetpack_shortlink":"https:\/\/wp.me\/p6XU0A-Rm","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/3308","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=3308"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/3308\/revisions"}],"wp:attachment":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media?parent=3308"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=3308"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=3308"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}