Quelle est la source des rayures horizontales et verticales dans les DEM USGS?

10

Lors du traitement des données DEM de 30 m et 10 m téléchargées via la visionneuse de carte nationale du jeu de données national d'élévation, nous avons remarqué des rayures horizontales et verticales non seulement dans les résultats produits, mais dans les nuances analytiques des collines des DEM bruts. Quelqu'un connaît-il la source? Si ce n'est pas la source, peut-être comment supprimer ces artefacts? Ces artefacts deviennent très prononcés lors de l'utilisation des DEM pour calculer les indices topographiques. Ces artefacts restent même après le remplissage de la dépression.

Vous trouverez ci-dessous des images montrant le rayage sur 30 m et 10 m de données de bassins versants en Pennsylvanie et au Colorado, et un calcul d'index topographique terminé montrant les artefacts d'un bassin versant à Syracuse, NY.

Colorado - HUC8 - 10190004 - 10m

Colorado - HUC8 - 10190004 - 10m

Colorado - HUC8 - 10190004 - 30m

Colorado - HUC8 - 10190004 - 30m

Pennsylvanie - HUC8 - 02040103 - 10m

Pennsylvanie - HUC8 - 02040103 - 10m

Pennsylvanie - HUC8 - 02040103 - 30m

Pennsylvanie - HUC8 - 02040103 - 30m

Calcul TI terminé pour le bassin versant du ruisseau Onondaga à Syracuse, NY

entrez la description de l'image ici

traggatmot
la source

Réponses:

15

Essayez de répondre à ma propre question:

La cause de l'entrelacement dans les exemples que j'ai fournis est entièrement due à mon flux de travail, pas à un problème hérité de la façon dont les données ont été initialement assemblées ou mosaïquées ensemble. Les DEM auxquels je faisais face étaient tous générés à partir de techniques plus récentes, comme en témoigne cette carte:

entrez la description de l'image ici

Les deux méthodes qui couvrent les domaines avec lesquels je travaillais sont le LIDAR et d'autres capteurs actifs ou l'interpolation linéaire complexe. Les anciennes techniques référencées par @Dan Patterson sont les techniques de profilage manuel et de photomappeur Gestalt. En effet, l'USGS fait référence à cela dans le lien NED @Dan Patterson partage:

Les DEM de source plus ancienne produits par des méthodes qui sont désormais obsolètes ont été filtrés pendant le processus d'assemblage NED pour minimiser les artefacts que l'on trouve couramment dans les données produites par ces méthodes. L'élimination des artefacts améliore considérablement la qualité des informations sur la pente, le relief ombré et le drainage synthétique qui peuvent être dérivées des données d'altitude. Le processus de filtrage de la suppression des artefacts n'élimine pas tous les artefacts. Dans les zones où le seul DEM disponible est produit par des méthodes plus anciennes, le "striping" peut encore se produire. Le traitement du NED comprend également des étapes pour ajuster les valeurs là où les DEM adjacents ne correspondent pas bien, et pour remplir les zones de ruban de données manquantes entre les DEM. Ces étapes de traitement garantissent que le NED n'a pas de zones vides et des discontinuités artificielles minimales.

Alors, qu'est-ce qui a causé mes problèmes de rayures?

Alors que, pour calculer correctement les valeurs TI dans SAGA GIS, nous avons besoin que les unités cellulaires soient en mètres, pas la mesure de degré du système de coordonnées géographiques d'origine, et donc la première étape de notre flux de travail a consisté à utiliser ArcMAP (je déteste le jeu d'outils de projection de SAGA) pour projetez le DEM dans la projection UTM correcte. Dans cette étape, il existe différentes options pour rééchantillonner le DEM. Dans tous les DEM et les sorties résultantes qui avaient une répartition, nous avons incorrectement laissé la technique de rééchantillonnage par défaut comme notre choix- l'algorithme de rééchantillonnage par défaut est le plus proche voisin, qui ne devrait jamais être utilisé avec un ensemble de données continu comme les données d'évolution présentes dans un DEM. Lorsque les DEM ont été projetés à l'aide du rééchantillonnage d'interpolation bi-linéaire, aucun artefact horizontal ou vertical n'a été observé dans le DEM ou dans aucun des produits résultants.

ESRI le savait:

Les DEM sont sensibles aux artefacts. De nombreux DEM ont déjà des artefacts introduits lors de la création; les ombrages de ces MNT amplifieront les anomalies et les rendront visibles. Si le DEM ne contient aucun artefact avant d'être rendu en ombrage, le problème peut être dû à l'utilisation d'une méthode de rééchantillonnage incorrecte lors de la projection des données DEM. Un DEM est des données raster continues. La méthode de rééchantillonnage bilinéaire doit être utilisée dans les projections raster ou toutes les transformations raster. Lorsque vous projetez des données raster à l'aide de l'outil Project Raster GP, n'utilisez pas la méthode de rééchantillonnage par défaut. Choisissez plutôt le rééchantillonnage bilinéaire ou la méthode de rééchantillonnage par convolution cubique.

Source: http://support.esri.com/en/knowledgebase/techarticles/detail/29127

Et l'USGS le sait, déclarant dans la FAQ:

Q: Quelles méthodes de rééchantillonnage sont les meilleures pour préserver la précision des données NED et les caractéristiques du terrain?

R: La convolution cubique et l' interpolation bilinéaire sont les méthodes préférées de rééchantillonnage des données d'élévation numériques, et donneront une apparence plus lisse. Le voisin le plus proche a tendance à laisser des artefacts tels que des marches d'escalier et des rayures périodiques dans les données qui peuvent ne pas être visibles lors de la visualisation des données d'élévation mais peuvent affecter les dérivés, tels que le relief ombré ou les rasters de pente. *

Source: http://ned.usgs.gov/faq.html#RESAMPLE

Donc, mon acceptation stupide des paramètres par défaut dans ArcMap (et mon ignorance des résultats) a provoqué cela. Une erreur très évidente probablement.

Vivre et apprendre.

traggatmot
la source
1

Il y a quelques articles sur la source possible de l'entrelacement qui sont trop longs pour être copiés ici et publiés, par exemple

https://geonet.esri.com/message/248734?sr=search&searchId=8194652f-cac8-4737-93a2-c5dccdeb29ff&searchIndex=5#248734

http://ned.usgs.gov/about.html

http://www.ctmap.com/assets/pdfprojects/destripe.pdf

Certains problèmes sont liés à la nature des données elles-mêmes.


la source
Merci pour ces liens utiles. Ils m'ont aidé à comprendre comment j'observerais les données de segmentation ou de profilage manuel (pour utiliser la terminologie ici: ctmap.com/assets/pdfprojects/destripe.pdf ) ou les données traitées automatiquement à l'aide de GPM (Gestalt Photomapper). Mais ils ont également conduit à des sources indiquant que les DEM que j'ai téléchargés sont dérivés de techniques de traitement plus récentes. Par exemple, le jeu de données du Colorado a été généré à l'aide de "l'interpolation linéaire complexe et le jeu de données PA à l'aide des données LIDAR. Ces techniques, je pense, devraient considérablement réduire l'entrelacement.
N'est
Veuillez essayer de résumer les informations dans les liens.
Aaron