Suite

Obtenir d'étranges artefacts de «grille» lors de la création d'ombrages à l'aide de gdaldem

Obtenir d'étranges artefacts de «grille» lors de la création d'ombrages à l'aide de gdaldem


J'essaie de créer des ombrages à partir de données d'altitude tchèques (Base fondamentale de données géographiques de la République tchèque (ZABAGED®) - altimétrie - grille 10x10 m) - un fichier de démonstration est disponible ici : http://geoportal.cuzk.cz/ UKAZKOVA_DATA/GRID10x10.zip

Je combine certains de leurs fichiers txt à l'aide de scripts bash, puis je crée un GeoTiff à l'aide de gdal_grid. Le GeoTiff résultant ressemble à ceci lorsqu'il est importé dans QGIS :

Comme prochaine étape, j'aimerais créer des ombrages à l'aide de Raster->Analyse->DEM et le résultat ressemble à ceci :

Je me suis assuré d'utiliser l'option bilinéaire lors de la déformation du raster et j'ai déjà essayé pratiquement tous les algorithmes disponibles de gdal_grid.

Je ne sais pas si cela est pertinent, mais c'est à quoi ressemble le TIFF d'ombre de colline lorsqu'il est ouvert dans OS X Preview :

Quelle est la source de ces artefacts et comment les éviter ?


Utilisez QGIS Vector to Raster pour convertir les points du shapefile en raster, je vais essayer d'expliquer pourquoi :

L'utilisation de l'utilitaire GDAL_Grid a interpolé de manière incorrecte, c'est de là que vient le pas à pas, vous ne le voyez tout simplement pas dans un moteur de rendu noir à blanc. Voici comment je vois les exemples de données interpolées à l'aide de GDAL_Grid dans Esri : Notez les bandes horizontales.

Utilisation des paramètres :

gdal_grid -ot float32 -of GTIFF -zfield Z -l grid10x10 -outsize 517 422

L'interpolation fonctionne mieux (517 par 422 a été calculé à partir de l'étendue divisée par 10) produisant l'ombrage :

Remarque : le cerclage est meilleur mais reste visible.

Le baguage est introduit par le programme GDAL_Grid ! En utilisant QGIS Raster::Conversion::Rasterize (Vector to Raster), remplissez la taille de la cellule afin qu'elle ne soit pas interpolée, puis l'ombrage de Raster::Terrain Analysis::Hillshade:

Comme l'image n'est pas interpolée, il y a des points qui ne sont pas remplis et peuvent nécessiter une moyenne focale. Après avoir réglé ça, ça devrait aller !


Je vous suggère d'utiliser les outils de création de raster SAGA dans la boîte à outils de traitement. Vous devrez peut-être les activer dans vos paramètres. Ces outils ne montrent pas des artefacts aussi horribles.