{"id":9708,"date":"2021-10-26T10:14:47","date_gmt":"2021-10-26T08:14:47","guid":{"rendered":"https:\/\/www.sigterritoires.fr\/?p=9708"},"modified":"2021-10-26T10:14:48","modified_gmt":"2021-10-26T08:14:48","slug":"classification-superviseelissage-des-frontieres-de-classe-avec-arcmap","status":"publish","type":"post","link":"https:\/\/www.sigterritoires.fr\/index.php\/classification-superviseelissage-des-frontieres-de-classe-avec-arcmap\/","title":{"rendered":"Classification supervis\u00e9e:lissage des fronti\u00e8res de classe avec ArcMap"},"content":{"rendered":"\n

Dans cette s\u00e9rie d’articles<\/a> nous avons vu un tutoriel qui aborde la classification d\u2019images a\u00e9riennes de drone en milieu littoral. Il utilise QGis dans sa version 3.20 et le plugin SCP dans sa version 7.3.<\/p>\n\n\n\n

Le but de ce tutoriel est de trouver un moyen de multiplier les trois bandes RGB des images a\u00e9riennes en utilisant des calculs de texture de sol et de l\u2019analyse en composantes principales . Dans le dernier article nous avons vu comment traiter les sorties d’une classification supervis\u00e9e pour int\u00e9grer les r\u00e9sultats dans un SIG<\/a>. <\/p>\n\n\n\n

Ici nous verrons comment lisser les fronti\u00e8res de classe du r\u00e9sultat raster de la classification.<\/p>\n\n\n\n

Cette op\u00e9ration n\u2019est pas disponible dans QGis. Elle peut se faire avec ArcMap avec l\u2019outil \u00ab Nettoyage de limites \u00bb. Avec QGis il est n\u00e9cessaire de vectoriser le raster class\u00e9 et de lisser les polygones r\u00e9sultants avec les outils de traitement vecteur.<\/p>\n\n\n\n\n\n\n\n

Lissage des fronti\u00e8res de classe<\/a><\/h1>\n\n\n\n
L’outil Nettoyage de limites<\/strong> lisse les bords irr\u00e9guliers de limites de classe et regroupe les classes. Cette op\u00e9ration augmente la coh\u00e9rence spatiale dans les classes. Les r\u00e9gions adjacentes qui appartiennent \u00e0 la m\u00eame classe peuvent \u00eatre connect\u00e9es.<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Outre la d\u00e9finition du raster en entr\u00e9e et du raster en sortie, vous avez deux options. La premi\u00e8re est la technique de tri.<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Ce menu sp\u00e9cifie le type de tri \u00e0 utiliser dans le processus de lissage.<\/p>\n\n\n\n
Cela d\u00e9termine la priorit\u00e9 par laquelle les cellules peuvent se d\u00e9velopper dans leurs voisins.<\/p>\n\n\n\n
NO_SORT\u2014 Ne fait pas de tri par taille. Les zones avec des valeurs plus grandes ont une priorit\u00e9 plus \u00e9lev\u00e9e pour se d\u00e9velopper en zones avec des valeurs plus petites. C\u2019est la valeur par d\u00e9faut.<\/p>\n\n\n\n
DESCEND \u2013 Trie les zones par ordre d\u00e9croissant de taille. Les zones ayant des zones totales plus grandes ont une priorit\u00e9 plus \u00e9lev\u00e9e pour s\u2019\u00e9tendre en zones ayant des zones totales plus petites.<\/p>\n\n\n\n
ASCEND \u2013 Trie les zones par ordre croissant de taille. Les zones dont la superficie totale est plus petite ont une priorit\u00e9 plus \u00e9lev\u00e9e pour s\u2019\u00e9tendre \u00e0 des zones dont la superficie totale est plus grande.<\/p>\n\n\n\n
Le lissage des bords se fait au d\u00e9triment soit du patch, soit de la zone qui l\u2019entoure. L\u2019option ASCEND pr\u00e9serve les patchs en rognant sur la zone environnante. L\u2019option DESCEND rogne plut\u00f4t les petits patchs.<\/p>\n\n\n\n
Sur la zone pr\u00e9c\u00e9dente, voici les r\u00e9sultats :<\/p>\n\n\n\n
$\"\"$ <\/a>
Avec l\u2019option ASCEND<\/figcaption><\/figure>\n\n\n\n
$\"\"$ <\/a>
Avec l\u2019option DESCEND<\/figcaption><\/figure>\n\n\n\n
La case \u00e0 cocher \u00ab Expansion et r\u00e9tr\u00e9cissement deux fois (facultatif) \u00bb sp\u00e9cifie le nombre de directions dans lesquelles le processus de lissage aura lieu.<\/p>\n\n\n\n
Coch\u00e9 : Effectue l\u2019expansion et le r\u00e9tr\u00e9cissement selon le type de tri sp\u00e9cifi\u00e9, puis effectue un r\u00e9tr\u00e9cissement et une expansion suppl\u00e9mentaires avec la priorit\u00e9 invers\u00e9e.<\/p>\n\n\n\n
C\u2019est la valeur par d\u00e9faut.<\/p>\n\n\n\n
Non coch\u00e9 \u2013 Effectue une expansion et une contraction une fois, selon le type de tri.<\/p>\n\n\n\n
G\u00e9n\u00e9ralisation de la sortie en supprimant de petites r\u00e9gions isol\u00e9es<\/a><\/h1>\n\n\n\n
Ce processus supprime de petites r\u00e9gions isol\u00e9es d’une image class\u00e9e. Les r\u00e9gions de taille sup\u00e9rieure \u00e0 un certain nombre de pixels restent sur l’image. Les outils d\u2019ArcMap Groupes par r\u00e9gions, SetNull et Grignotage permettent d’accomplir cette t\u00e2che.<\/p>\n\n\n\n
Cette op\u00e9ration est plus compliqu\u00e9e que ce qu\u2019elle para\u00eet au d\u00e9part. En effet, on l\u2019associe \u00e0 la suppression de polygones avec une taille inf\u00e9rieure \u00e0 une limite donn\u00e9e. Mais en r\u00e9alit\u00e9 elle consiste \u00e0 d\u00e9terminer des zones de pixels ayant un nombre inf\u00e9rieur de pixels \u00e0 une valeur donn\u00e9e et de leur affecter la valeur des pixels qui entourent cette zone.<\/p>\n\n\n\n
C\u2019est pourquoi il nous faut trois op\u00e9rations pour r\u00e9aliser ce traitement :<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Grouper par r\u00e9gion : va nous permettre de donner un identifiant commun aux pixels appartenant \u00e0 chaque zone de notre image<\/li>
La table attributaire du raster en sortie comporte une colonne \u00ab Total \u00bb qui indique le nombre de pixels de chaque groupe. C\u2019est \u00e0 partir de ces informations que vous d\u00e9terminerez le seuil que vous allez appliquer, c\u2019est-\u00e0-dire le nombre minimum de pixels que devront avoir les zones pour \u00eatre conserv\u00e9es au final.<\/li>
L\u2019outil SetNull va nous permettre de cr\u00e9er un masque avec une valeur de 1 pour les pixels des zones \u00e0 conserver.<\/li>
L\u2019outil \u00ab Grignotage \u00bb utilise ce masque pour remplir les pixels \u00e0 \u00e9liminer avec la valeur des pixels environnants. Les pixels correspondants \u00e0 la valeur 1 du masque ne sont pas touch\u00e9s, les autres pixels sont modifi\u00e9s pour leur affecter la valeur du pixel masqu\u00e9 le plus proche.<\/li><\/ol>\n\n\n\n
Exemple<\/a><\/h2>\n\n\n\n
Nous partons de l\u2019image liss\u00e9e :<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
La premi\u00e8re \u00e9tape est le calcul des r\u00e9gions de notre image class\u00e9e :<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Ce qui donne :<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Chaque \u00ab paquet \u00bb de pixels a maintenant une seule valeur qui l\u2019identifie.<\/p>\n\n\n\n
Nous ouvrons la table attributaire du raster r\u00e9sultant et nous pouvons trier en ordre croissant la colonne Total, puis s\u00e9lectionner des enregistrements pour visualiser les paquets ayant moins d\u2019une certaine quantit\u00e9 de pixels. Dans l\u2019image suivante nous avons s\u00e9lectionn\u00e9 les paquets ayant moins de 50 pixels :<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Supposons, pour la suite, que nous choisissons d\u2019\u00e9liminer les zones ayant moins de 15 pixels.<\/p>\n\n\n\n
Nous utilisons l\u2019outil SetNull pour cr\u00e9er un masque :<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
On ex\u00e9cute cet outil sur notre raster de r\u00e9gions (pas sur l\u2019image class\u00e9e).<\/p>\n\n\n\n
L\u2019expression contient la condition souhait\u00e9e : nombre de pixels du groupe inf\u00e9rieur \u00e0 15.<\/p>\n\n\n\n
La valeur pour le r\u00e9sultat de l\u2019\u00e9valuation de l\u2019expression comme fausse (les paquets ayants plus de 15 pixels) est mise \u00e0 1.<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
Les zones correspondantes aux paquets inf\u00e9rieurs \u00e0 15 pixels ont la valeur NoData.<\/p>\n\n\n\n
Nous allons utiliser maintenant l\u2019outil Grignotage pour affecter \u00e0 ces zones la valeur du pixel le plus proche. Le masque d\u00e9fini sur quelles parties de l\u2019image s\u2019effectue cette op\u00e9ration. Seules les parties \u00ab trou\u00e9es \u00bb du masque permettent l\u2019op\u00e9ration.<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
On ex\u00e9cute cet outil sur l\u2019image class\u00e9e (pas sur l\u2019image des r\u00e9gions).<\/p>\n\n\n\n
Le r\u00e9sultat ne contient plus les zones de moins de 15 pixels :<\/p>\n\n\n\n
$\"\"$ <\/a><\/figure>\n\n\n\n
L’image est pr\u00eate pour \u00eatre vectoris\u00e9e en polygones pour son int\u00e9gration dans un SIG.<\/p>\n","protected":false},"excerpt":{"rendered":"
Dans cette s\u00e9rie d’articles nous avons vu un tutoriel qui aborde la classification d\u2019images a\u00e9riennes de drone en milieu littoral. Il utilise QGis dans sa version 3.20 et le plugin SCP dans sa version 7.3. Le…<\/p>\n","protected":false},"author":1,"featured_media":9723,"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":[145,209,146,147,152,26,62,172],"tags":[65,1480,1690,1693,1691,1689,1692,1486],"class_list":["post-9708","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-arcgis","category-arcgisqgis","category-arcmap","category-images","category-outils_divers","category-outils-et-methodes","category-qgis-2","category-tutoriel","tag-arcmap","tag-classification","tag-generalisation","tag-grignotage","tag-groupe-par-regions","tag-lissage-classes","tag-setnull","tag-supervisee"],"aioseo_notices":[],"jetpack_featured_media_url":"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2021\/10\/i1-3.png?fit=1037%2C492&ssl=1","jetpack_shortlink":"https:\/\/wp.me\/p6XU0A-2wA","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/9708","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=9708"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/9708\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media\/9723"}],"wp:attachment":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media?parent=9708"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=9708"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=9708"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}