Suite

Fusionner correctement les rasters DEM et éliminer les effets de bord


J'utilise des rasters DEM de 10 m de la passerelle de données NRCS pour une grande partie de la Californie. Le problème est une mosaïque étrange qui apparaît après la fusion, mais n'est pas présente dans les données d'origine (avant la fusion). Les grilles créent des points bas, de sorte que l'accumulation de flux, etc. y est attirée.

J'ai essayé de mosaïquer vers un nouveau raster, de créer un catalogue de raster et de filtrer dans ArcGIS 10.2, de fusionner et de créer vrt dans Qgis2.4, mais le problème persiste.

Les grilles représentent les creux plats dans le raster de courbure sur les bords fusionnés des rasters DEM. Notez que ces lignes sont les bords des DEUX du MNT mosaïqué, créant 2 lignes.

Les lignes sont des dépressions plates, ce qui crée un problème en essayant de faire n'importe quelle sorte d'analyse de stabilité de pente ou d'analyse hydrologique… ce qui est tout l'intérêt de cet exercice.

Rasters DEM non mosaïques :

Raster de courbure généré à partir d'un raster en mosaïque à l'aide d'une interpolation bilinéaire :

EDIT : les données et la méthode de la mosaïque étaient les problèmes : les données elles-mêmes ont des effets de bord intégrés et j'utilisais le voisin le plus proche au lieu de l'interpolation bilinéaire


Il semble qu'il se passe deux choses ici. Premièrement, vos tuiles ne sont pas homogènes, dans la mesure où dans les zones de chevauchement aux bords des tuiles, les élévations ne sont pas identiques. Je peux le confirmer car j'ai numérisé plusieurs points le long de la zone de chevauchement et extrait les élévations dans les deux DEM raster et trouvé ceci :

Dans le cas des deux tuiles que vous avez téléchargées, il semble que les DEM soient très similaires près du bas, mais à mesure que vous vous déplacez vers le bord supérieur, les écarts d'élévation ont augmenté jusqu'à près d'un demi-mètre, ce qui est assez important. Vous ne pouvez pas vraiment vous attendre à ce qu'un DEM mosaïqué ne montre pas d'artefacts de bord si les tuiles composites ne sont pas exactement les mêmes dans la zone de chevauchement. C'est-à-dire, à moins que vous n'effectuiez un mosaïquage avec une sorte de mélange, et même dans ce cas, il est probable que certains effets de bord subsisteront étant donné les grands écarts. Lequel des types de mosaïque avez-vous choisi lorsque vous avez mosaïqué les carreaux ? ArcGIS propose un mode « fusionner ».

Le deuxième problème semble être lié à la manière dont ArcGIS gère réellement les données près des bords des tuiles composites pendant le mosaïquage. Je crains de ne pas pouvoir vous en dire plus sans avoir la possibilité d'examiner le code source d'Arc, ce qui n'arrivera jamais. Cependant, jetez un œil aux images suivantes, l'image de courbure du plan résultante dérivée de la mosaïque de vos tuiles dans un autre logiciel SIG open source :

Les bords ont été mis en évidence avec un vecteur superposé pour montrer les emplacements des bords (certains des points que j'ai utilisés précédemment s'y trouvent également). Maintenant, sans ces vecteurs de frontière de bord superposés :

Notez que les artefacts de bord que vous mettez en évidence dans votre question ne sont pas aussi évidents dans cette mosaïque, qui n'applique pas non plus les techniques de fusion. Ce sont les mêmes données source DEM, la seule différence est l'algorithme de mosaïque spécifique.

Donc, pour conclure, je dirais 1) que vos carreaux ne sont pas homogènes, donc sans utiliser un algorithme de mosaïquage plus sophistiqué impliquant un dégradé, vous devriez vous attendre à des effets de bord, et 2) dans une certaine mesure, l'étendue des artefacts que vous observez est le résultat de la façon dont l'algorithme de mosaïquage d'ArcGIS gère les bords des tuiles composites. J'espère que ça aide.


Nous vous suggérons d'essayer de créer un raster virtuel à partir des fichiers DEM TIFF : Raster -> Divers -> Build Virtual Raster (Catalog). Cela semble fonctionner avec vos exemples de données. Désolé, je devais juste mettre un peu d'ombrage sur le résultat. N.


Voir la vidéo: КОФТОЧКА,,ЛАВАНДА,, КРУГЛАЯ КОКЕТКА. SWEATER ROUND YOKEÖRGÜ KAZAKSUÉTER TEJIDOMAGLIONE A MAGLIA (Octobre 2021).