Combiner des rasters ombrage haute résolution et généralisés pour produire un relief ombré amélioré?

J'ajoute souvent des données LiDAR à une carte avec un relief ombré. La plupart du temps, je dois généraliser le DEM LiDAR en fonction de l'échelle de la carte. Pourtant, il y a souvent trop de détails et les crêtes et les vallées des montagnes ne sont pas bien prononcées.

Je suis tombé récemment sur un article qui décrit le processus de combinaison de 2 modèles d'ombrage, l'un détaillé et l'autre généralisé. Le relief qui en résulte est très agréable. Je voudrais produire un tel effet dans mes cartes: en supposant que j'ai 2 rasters ombrage (détaillés et généralisés), quelle est la meilleure façon de les combiner pour obtenir un tel effet?

Cela ressemble au travail de Tom Patterson sur Resolution Bumping GTOPO 30 dans Photoshop . La théorie est suffisamment bien décrite pour être adaptable à d'autres logiciels, bien que des travaux doivent être effectués pour arriver aux détails. L'idée de base est de généraliser (brouiller) un ensemble de données, beaucoup, pour souligner la forme générale et masquer des détails spécifiques, puis mélanger les images haute résolution et basse résolution ensemble à l'aide de canaux alpha ou de pourcentages d'opacité.

Alors que le processus décrit par Patterson est appliqué aux modèles d'élévation eux-mêmes, je l'ai adapté aux reliefs comme suit.

Chargez l'image en relief dans GIMP:

• calque en double
• calque supérieur == le calque de détail, laissez-le tranquille, réglez Mode sur Overlay
• en bas == le calque de forme, appliquer un flou gaussien avec le réglage 20
• enregistrer une copie (relief_blurred.tif)

Réaffecter la projection:

• gdalcopyproj shaded_relief.tif relief_blurred.tif

Il y a beaucoup de place pour l'adaptation dans le choix des pourcentages à utiliser, que ce soit pour utiliser des filtres Overlay ou Screen ou Darken, etc.

Je vois "ArcGIS" est une balise, Jakub. À l'aide de Spatial Analyst, vous calculez simplement une moyenne pondérée des deux ombrages. Par exemple, le mélange 60-40 pourrait être généré avec un calcul comme celui-ci:

(60*[Detailed hillshade] + 40*[Generalized hillshade]) / 100

Si vous en avez besoin, le flou gaussien peut être exécuté en exécutant quelques moyens focaux circulaires sur un raster. Le rayon effectif de faire m moyennes focales, chacun de rayon r , est r * Sqrt ( m ).