![\"\"](\"https:\/\/4.bp.blogspot.com\/-E029tymfai4\/WgcyzsefLgI\/AAAAAAAAA20\/zzbv5ZpAtnQgL9OgN2w4s1Qw6oLHwNK5gCLcBGAs\/s1600\/IGN_25.png\")
<\/a>\nLe MNT dans son \u00e9tat cru<\/em><\/p>\n
Il faut, donc, lisser l\u2019affaire pour enlever les marches plates. Une approche (d\u00e9j\u00e0 utilis\u00e9e par certains) passe par la transformation du MNT en courbes de niveau et la r\u00e9g\u00e9n\u00e9ration du mod\u00e8le. Cependant, c\u2019est une approche lourde en calcul qui, de surcroit, introduit des potentielles erreurs dues aux d\u00e9calages des courbes vectorielles. Nous pouvons omettre l\u2019histoire des courbes de niveau vectorielles en nettoyant le MNT directement.<\/p>\n J\u2019ai fait un script d\u2019extraction des \u00ab marches \u00bb pour QGIS car j\u2019avais besoin de nettoyer le MNT pour mes analyses. Le script se trouve sur GitHub (lien ci-dessus<\/strong>).<\/p>\n Le script apparaitra sous la rubrique \u00ab Raster \u00bb de la boite d\u2019outils \u00ab Processing \u00bb. Il extrait tous les bordures des marches, c\u2019est-\u00e0-dire l\u2019empreinte des courbes de niveau en raster. Deux m\u00e9thodes sont int\u00e9gr\u00e9es pour le faire : standard, o\u00f9 les pixels des bordures se touchent aussi bien par les coins que par les c\u00f4t\u00e9s, et la \u00ab Ligne \u00e9paisse \u00bb qui produit des lignes plus \u00e9tanches o\u00f9 seulement les connections par les c\u00f4t\u00e9s sont admises. Le r\u00e9sultat \u00ab standard \u00bb est sur la figure ci-dessus.<\/p>\n La deuxi\u00e8me \u00e9tape passe par un outil d\u2019interpolation des valeurs dans les zones \u00e9limin\u00e9es. Le La derni\u00e8re touche, un lissage qui \u00e9liminera quelques effets sautillants r\u00e9siduels, peut \u00eatre faite par un bon nombre d\u2019outils. Le Enfin, l\u2019analyse de visibilit\u00e9 donne un r\u00e9sultat au moins visuellement satisfaisant ! Le lissage risque de d\u00e9grader quelque peu le MNT, mais l\u2019interpolation, si bien faite, ajoute des informations suppl\u00e9mentaires qui rendent l\u2019analyse plus pr\u00e9cise.<\/p>\n Script pour QGIS \u201cLissage BD ALTI\u201d<\/a> (versions 2 et 3).<\/p>\n<\/a>
\nL\u2019analyse de visibilit\u00e9 : le premier r\u00e9sultat<\/em><\/p>\n<\/a>
\nLes zones de transition dans le MNT (les isohypses)<\/em><\/p>\nr<\/span>.<\/span>surf<\/span>.<\/span>contour<\/span><\/code> de GRASS est fait pour ce type de calcul (l\u2019interpolation du MNT \u00e0 partir des courbes de niveau en format raster) : il donne un tr\u00e8s bon r\u00e9sultat.<\/p>\nSimple<\/span> filter<\/span><\/code> de SAGA a \u00e9t\u00e9 utilis\u00e9 ici pour obtenir une belle surface lisse.<\/p>\n<\/a>
\nLe resultat de r<\/span>.<\/span>surf<\/span>.<\/span>contour<\/span><\/code> de GRASS : tr\u00e8s bien mais pas tout \u00e0 fait lisse.<\/em><\/p>\n<\/a>
\nL\u2019analyse de visibilit\u00e9 sur un MNT liss\u00e9 : enfin un r\u00e9sultat acceptable !<\/em><\/p>\nSources<\/h2>\n