Pourquoi l'outil «Contour» renvoie-t-il des contours avec un motif rectangulaire étrange?

8

Je génère des contours de 2 pieds à partir d'un MNT lissé d'une résolution horizontale de 3,2808 pieds (1 m). Le résultat montre ces carrés impairs qui font environ 80 pieds de côté. Le lissage a été fait en utilisant 8 voisins ( EDIT: je pensais que j'utilisais la majorité focale mais j'ai vérifié mes notes et j'ai effectivement utilisé une moyenne focale . J'avais aussi tort sur les 8 voisins. C'était un quartier 9x9 .). Le DEM n'affiche pas ces carrés (voir la deuxième capture d'écran ci-dessous).

J'ai dessiné l'un des carrés avec l'outil de dessin pour le démontrer, mais je pense qu'ils sont assez évidents.

Q) Quelqu'un peut-il fournir une explication à ce sujet ou une solution potentielle?

Ce résultat provient de l'outil "Contour" dans "Outils Spatial Analyst -> Surface -> Contour" et j'utilise 10.3.1

Contours montrant un motif rectangulaire

Deuxième capture d'écran montrant le DEM lissé (même zone générale mais zoomée un peu) entrez la description de l'image ici

jbchurchill
la source
Il le fait, en particulier dans les zones plates, il a quelque chose à voir avec la logique de cellule oui / non et les blocs de lecture. J'ai trouvé qu'un bon lissage les rend moins évidents (rééchantillonnage et / ou moyenne focale). Vous pouvez essayer d'utiliser GDAL_Contour ou si vous disposez d'un terrain, d'un jeu de données LAS ou d'un travail TIN à partir de ceux-ci - mais de mémoire, ces méthodes sont identiques ou pires.
Michael Stimson
J'adorerais apporter mon soutien, mais la seule action que je puisse suggérer est de renvoyer DEM au fournisseur. Quelque chose de terriblement mal avec le recouvrement / la manipulation des chevauchements (?) Qu'ils ont appliqués
FelixIP

Réponses:

2

Par défaut, toute interpolation (projection / rééchantillonnage) sera le plus proche voisin et c'est le problème la plupart du temps (destiné aux données catégorielles NON continues). Je dirais que la majorité focale est un mauvais choix car j'imagine que les cellules ont toutes des valeurs flottantes uniques. Mieux vaut la moyenne. Vous pouvez également modifier l'interpolation en bilinéaire et utiliser simplement le rééchantillonnage (Gestion des données> Raster> Traitement raster> Rééchantillonner) par opposition à la moyenne focale.

Mike
la source
Je suis assez certain d'avoir utilisé bilinéaire pour la projection. Je vais devoir revoir mes notes. De plus, j'ai peut-être oublié le filtre que j'ai utilisé (la logique de l'utilisation de la moyenne focale est évidente pour moi et c'est peut-être ce que j'ai réellement utilisé. Merci pour les idées.
jbchurchill
Vous pouvez essayer d'aller dans C: \ Program Files (x86) \ ArcGIS \ Desktop10.1 \ Utilities \ AdvancedArcMapSettings.exe (ou où que ce soit pour vous) et simplement changer la valeur par défaut en bilinéaire dans l'onglet Raster jusqu'à ce que cela fonctionne . Juste au cas où quelque chose se passe en arrière-plan. Je me souviens vaguement d'un dossier de dems fourni comme valeur entière dans un système géographique que je devais convertir en float avant de projeter ou bien il voulait sortir un entier ... ce qui rendait le pas plus important ... pourrait être lié au type de pixel peut-être aussi
Mike