{"id":3961,"date":"2017-04-04T14:01:13","date_gmt":"2017-04-04T13:01:13","guid":{"rendered":"http:\/\/www.sigterritoires.fr\/?p=3961"},"modified":"2017-04-04T14:01:13","modified_gmt":"2017-04-04T13:01:13","slug":"lanalyse-spatiale-avec-sql-4-fonctions-de-traitement-spatial","status":"publish","type":"post","link":"https:\/\/www.sigterritoires.fr\/index.php\/lanalyse-spatiale-avec-sql-4-fonctions-de-traitement-spatial\/","title":{"rendered":"L’analyse spatiale avec SQL: 4-fonctions de traitement spatial"},"content":{"rendered":"

Dans l’article pr\u00e9c\u00e9dent nous avons vu les fonctions SQL de relation spatiale<\/a>. Voyons maintenant les fonctions de traitement spatial de SQL.Commen\u00e7ons par voir quelle est la diff\u00e9rence entre la notion de relation spatiale et de traitement spatial.<\/p>\n

Le plus simple est de prendre un exemple. Nous avons des g\u00e9om\u00e9tries qui partagent l’occupation de l’espace<\/p>\n

$\"\"$ <\/a>Une fonction de relation spatiale, telle que st_intersects<\/a><\/strong><\/em>, va r\u00e9pondre \u00e0 la question \u00ab\u00a0est-ce que ces g\u00e9om\u00e9tries s’intersectent?\u00a0\u00bb par VRAI ou FAUX.<\/p>\n

Une fonction de traitement spatial, telle que st_intersection<\/strong><\/em>, va extraire la partie commune des g\u00e9om\u00e9tries en renvoyant comme r\u00e9sultat une g\u00e9om\u00e9trie :<\/p>\n

$\"\"$ <\/a>Supposons que nous avons deux tables: rectangles et cercles, et que nous voulons extraire le r\u00e9sultat de cette image. Comment construit-on la requ\u00eate SQL appropri\u00e9e?<\/p>\n

En faisant la m\u00eame d\u00e9marche que d\u00e9crite dans le deuxi\u00e8me article de cette s\u00e9rie<\/a>:<\/p>\n

\n
quelles sont les donn\u00e9es que je veux obtenir?: je veux les g\u00e9om\u00e9tries r\u00e9sultantes de l’intersection des cercles et des rectangles->clause SELECT:<\/li>\n<\/ul>\n
SELECT st_intersection(cercles.geometry,rectangles.geometry) as geometry
\n<\/strong><\/em><\/p>\n
\n
o\u00f9 est-ce qu\u2019elles sont stock\u00e9es? : dans les tables cercles et rectangles->clause FROM:<\/li>\n<\/ul>\n
FROM cercles,rectangles<\/strong><\/em><\/p>\n
\n
quelles sont les conditions que je veux appliquer, spatiales ou pas? s\u00e9lectionner seulement les entit\u00e9s qui occupent un espace commun->clause WHERE:<\/li>\n<\/ul>\n
\u00a0WHERE st_intersects(cercles.geometry,rectangles.geometry)<\/strong><\/em><\/p>\n
Ce qui nous fait la requ\u00eate:<\/p>\n
SELECT st_intersection(cercles.geometry,rectangles.geometry) as geometry<\/strong><\/em>
\nFROM cercles,rectangles<\/strong><\/em>
\nWHERE st_intersects(cercles.geometry,rectangles.geometry)<\/strong><\/em><\/p>\n
Vous pouvez omettre la clause WHERE et vous obtiendrez le m\u00eame r\u00e9sultat, sauf que \u00e7a prendra beaucoup plus de temps. C’est logique de limiter l’op\u00e9ration d’intersection de st_intersection aux endroits o\u00f9 les g\u00e9om\u00e9tries s’intersectent.<\/p>\n
Les outils batch et les outils interactifs<\/h2>\n
Avant de voir les fonctions de traitement spatial dans le d\u00e9tail, faisons une petite digression. Dans les SIG bureautique, tel que ArcGis, vous avez l’habitude de travailler avec deux groupes d’outils diff\u00e9rents selon qu’il s’agisse de traiter la couche dans son ensemble ou bien des entit\u00e9s en particulier.<\/p>\n
Quand il s’ag\u00eet de traiter la couche dans son ensemble, vous utiliserez les outils de la Toolbox. Quand il s’ag\u00eet de traiter des entit\u00e9s en particulier vous utiliserez les outils interactifs d’ArcMap.<\/p>\n
Un avantage du SQL est que vous ne travaillez qu’avec un seul outil. Reprenons la requ\u00eate pr\u00e9c\u00e9dente pour l\u2019appliquer \u00e0 l’intersections entre une couche carto_risques(en beige)\u00a0 et zones_inondables (en bleu)<\/p>\n
$\"\"$ <\/a>Nous pouvons utiliser la requ\u00eate:<\/p>\n
SELECT st_intersection(carto_risque.geometry,zones_inondables.geometry) as geometry<\/strong><\/em>
\nFROM carto_risques,zones_inondables<\/strong><\/em>
\nWHERE st_intersects(carto_risque.geometry,zones_inondables.geometry)<\/strong><\/em><\/p>\n
Pour intersecter toutes les entit\u00e9s pr\u00e9sentes dans la cartographie<\/a> des risques avec les zones inondables:<\/p>\n
$\"\"$ <\/a>Ou seulement la zone avec un risque maximal avec les zones inondables<\/p>\n
SELECT st_intersection(carto_risque.geometry,zones_inondables.geometry) as geometry<\/strong><\/em>
\nFROM carto_risques,zones_inondables<\/strong><\/em>
\nWHERE st_intersects(carto_risque.geometry,zones_inondables.geometry)<\/strong><\/em> AND
\ncarto_risque.zone=’max’<\/p>\n
$\"\"$ <\/a>La seule diff\u00e9rence se situe dans la clause WHERE<\/strong> o\u00f9 nous ajoutons \u00ab\u00a0AND<\/strong><\/em>
\n carto_risque.zone=’max’<\/strong><\/em>\u00a0\u00bb<\/p>\n
Les fonctions de traitement spatial<\/h2>\nSi vous regardez la page d’aide d’ArcGis du jeu d’outils Superposition (https:\/\/pro.arcgis.com\/fr\/pro-app\/tool-reference\/analysis\/an-overview-of-the-overlay-toolset.htm<\/a>) vous trouverez la description d’une s\u00e9rie de fonctions de traitement spatial:<\/p>\n
$\"\"$ <\/a> $\"\"$ <\/a>Ils sont au nombre de sept\u00a0 : Effacer, Identit\u00e9, Intersection, Jointure<\/a> spatiale, Diff\u00e9rence sym\u00e9trique, Agr\u00e9gation et Mise \u00e0 jour.<\/p>\n
En SQL nous utiliserons seulement trois fonctions: st_difference, st_intersection et st_union pour construire les requ\u00eates \u00e9quivalentes.<\/p>\n
Voyons tout d’abord ces trois fonctions et leur \u00e9quivalent dans le Toolset d’ArcGis, puis nous verrons comment obtenir l’\u00e9quivalence des autres outilst du Toolset en combinant les trois fonctions de base.<\/p>\n
st_Difference = Effacer (Erase)<\/h2>\n
Elle est de la forme st_difference(g\u00e9om\u00e9trie A, g\u00e9om\u00e9trieB). Cette fonction retourne la partie de la g\u00e9om\u00e9trie A qui n’intersecte pas la g\u00e9om\u00e9trie B.<\/p>\n
$\"\"$ <\/a>Dans l’exemple de la cartographie des risques et les zones inondables, on peut souhaiter, par exemple, voir la partie des zones class\u00e9es comme risque maximal qui ne sont pas inondables :<\/p>\n
SELECT st_Difference(carto_risque.geometry,zones_inondables.geometry) as geometry
\nFROM zones_inondables,carto_risque
\nWHERE carto_risque.zone=’max’<\/p>\n
Le r\u00e9sultat de la commande appara\u00eet en jaune.<\/p>\n
$\"\"$ <\/a>st_intersection = Intersection (Intersection)<\/h2>\n
Elle est de la forme st_intersection(g\u00e9om\u00e9trie A, g\u00e9om\u00e9trie B). Cette fonction retourne la partie de la g\u00e9om\u00e9trie A qui intersecte\u00a0 la g\u00e9om\u00e9trie B.<\/p>\n
$\"\"$ <\/a>Nous avons d\u00e9j\u00e0 d\u00e9taill\u00e9 la requ\u00eate<\/p>\n
SELECT st_intersection(carto_risque.geometry,zones_inondables.geometry) as geometry<\/strong><\/em>
\nFROM carto_risques,zones_inondables<\/strong><\/em>
\nWHERE st_intersects(carto_risque.geometry,zones_inondables.geometry)<\/strong><\/em><\/p>\n
dans le paragraphe pr\u00e9c\u00e9dent sur les fonctions batch et interactives. Mais nous allons compl\u00e9ter l’usage des fonctions avec un autre \u00e9l\u00e9ment. Dans les exemples utilis\u00e9s jusqu’ici nous obtenons comme r\u00e9sultat de la requ\u00eate que des g\u00e9om\u00e9tries. Mais la plupart du temps nous avons besoin de r\u00e9cup\u00e9rer aussi des attributs des entit\u00e9s concern\u00e9es.<\/p>\n
Par exemple, dans les entit\u00e9s r\u00e9sultantes de la requ\u00eate nous voulons savoir de quelle version de \u00ab\u00a0zones_inondables\u00a0\u00bb il s’ag\u00eet et de quelle niveau de risque. Nous ajoutons ces champs dans la partie SELECT de la requ\u00eate<\/p>\n
SELECT<\/strong> zones_inondables.version, carto_risque.zone, st_intersection(carto_risque.geometry, zones_inondables.geometry) as geometry
\nFROM<\/strong> carto_risque,zones_inondables
\nWHERE<\/strong> st_intersects(carto_risque.geometry, zones_inondables.geometry)<\/p>\n
et le r\u00e9sultat sera:<\/p>\n
$\"\"$ <\/a>Les champs attributs sont ajout\u00e9s dans la table r\u00e9sultante. Chaque g\u00e9om\u00e9trie aura la valeur de l’attribut des g\u00e9om\u00e9tries qui s’intersectent \u00e0 cet endroit.<\/p>\n
st_union<\/h2>\n
Elle est de la forme st_union(g\u00e9om\u00e9trie A, g\u00e9om\u00e9trie B). Cette fonction prend chaque g\u00e9om\u00e9trie de A et l’unifie avec la ou les g\u00e9om\u00e9tries de B, en fusionnant les points qui les composent.<\/p>\n
$\"\"$ <\/a>Si vous avez x g\u00e9om\u00e9tries en A, vous aurez le m\u00eame nombre de g\u00e9om\u00e9tries en r\u00e9sultat. Ce n’est pas le nombre qui change mais bien la g\u00e9om\u00e9trie de chaque entit\u00e9 de A qui se voit adjoindre l’espace occup\u00e9 par la ou les g\u00e9om\u00e9tries de B.<\/p>\n
Cette fonction SQL pr\u00eate \u00e0 confusion, surtout en anglais, avec l’op\u00e9ration Union (en fran\u00e7ais Agr\u00e9gation) du Toolset d’ArcGis. Voyons tout de suite celle-ci pour \u00e9viter toute m\u00e9prise.<\/p>\n
Agr\u00e9gation (Union en anglais)<\/h2>\n
Cette fonction permet de r\u00e9cup\u00e9rer tout l’espace qui est occup\u00e9 par les g\u00e9om\u00e9tries A et par les g\u00e9om\u00e9tries B:<\/p>\n
\n
les zones dans lesquelles les g\u00e9om\u00e9trie s’intersectent<\/li>\n
les zones o\u00f9 l’on ne trouve que des g\u00e9om\u00e9tries A<\/li>\n
les zones o\u00f9 l’on ne trouve que des g\u00e9om\u00e9tries B<\/li>\n<\/ul>\n
Le sch\u00e9ma suivant montre, en entr\u00e9e une couche A avec deux rectangles, et une couche B avec un cercle.<\/p>\n
$\"\"$ <\/a><\/p>\n
Le r\u00e9sultat de l’union de ces deux couches comprend 5 g\u00e9om\u00e9tries:<\/p>\n
\n
les deux parties des rectangles qui ne s’intersectent pas avec B<\/li>\n
les deux parties des rectangles qui s’intersectent avec le cercle<\/li>\n
la partie du cercle qui ne s’intersecte pas avec les rectangles<\/li>\n<\/ul>\n
Nous avons vu que la fonction st_union garde le m\u00eame nombre de g\u00e9om\u00e9tries de A dans le r\u00e9sultat. Dans l’outil Union du Toolset, il y a toujours plus de g\u00e9om\u00e9tries dans le r\u00e9sultat que dans la couche A, le ,nombre d\u00e9pendant des zones d’intersections des deux couches.<\/p>\n
La suite…<\/h3>\n
Nous avons vu\u00a0 les trois fonctions de base de SQL pour\u00a0 le traitement spatial . Dans le prochain article nous verrons comment combiner ces fonctions de base pour obtenir les fonctions telles que Identit\u00e9, Agr\u00e9gation et Mettre \u00e0 jour.<\/p>\n","protected":false},"excerpt":{"rendered":"
Dans l’article pr\u00e9c\u00e9dent nous avons vu les fonctions SQL de relation spatiale. Voyons maintenant les fonctions de traitement spatial de SQL.Commen\u00e7ons par voir quelle est la diff\u00e9rence entre la notion de relation spatiale et de traitement…<\/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_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":"","jetpack_post_was_ever_published":false},"categories":[311,172],"tags":[377,15,312,58,207],"class_list":["post-3961","post","type-post","status-publish","format-standard","hentry","category-postgresql","category-tutoriel","tag-analyse-spatiale","tag-arcgis","tag-postgresql","tag-qgis","tag-requete-sql"],"aioseo_notices":[],"jetpack_featured_media_url":"","jetpack_shortlink":"https:\/\/wp.me\/p6XU0A-11T","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/3961","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=3961"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/3961\/revisions"}],"wp:attachment":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media?parent=3961"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=3961"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=3961"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}