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?

— Jakub Sisak GeoGraphics
source

Réponses:

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.

— Matt Wilkie
source

très bonne méthode.

— Mapperz

ohhh je vais essayer ça, ça a l'air cool.

— Nathan W

J'espérais une solution ArcGIS mais j'essaierai cela aussi ...

— Jakub Sisak GeoGraphics

Je commence à penser que la seule façon de le faire doit être celle que vous avez suggérée dans un éditeur d'images comme GIMP, Photoshop, etc. (Étant donné que la rampe de couleur ArcGIS ne peut pas être définie du transparent au noir mais par défaut du blanc au noir)

— Jakub Sisak GeoGraphics

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 ).

— whuber
source

Merci Bill. L'ombrage détaillé de mon entrée est une taille de cellule 5, 5 mais l'ombrage généralisé est 20, 20. Le résultat de la calculatrice raster est donc également 20, 20. Il semble que cela améliore les ombres primaires mais les détails fins sont perdus. Que recommanderiez-vous pour conserver la taille de cellule d'origine? Rééchantillonner l'ombrage généralisé? Quelle technique de rééchantillonnage utiliser?

— Jakub Sisak GeoGraphics du

Je l'ai essayé, tout fonctionne mais le produit final de l'ombrage n'est pas très différent de l'original détaillé autre qu'un peu plus lisse ... (Plutôt que d'intensifier l'ombrage détaillé original)

— Jakub Sisak GeoGraphics

L'idée centrale mais inexprimée dans la réponse de Bill est que presque tout ce que vous pouvez faire dans Gimp / Photoshop / etc. en ce qui concerne le mélange des couches et des superpositions peut être fait avec un calcul raster - aussi longtemps que vous connaissez les mathématiques. Ce n'est pas facile, c'est pourquoi il n'y a pas beaucoup de solutions toutes faites. En partant de la théorie seule, je ne suis pas un chef des mathématiques, je pense que ce qui manque dans l'exemple ci-dessus est un troisième raster pour agir comme curseur d'opacité / transparence d'un pauvre.

— matt wilkie

@Jakub: Vous avez raison, vous devez conserver la taille de cellule la plus fine, sinon tout est perdu. Dans les produits ESRI, cela se fait généralement en définissant au préalable une variable "d'environnement" quelque part pour spécifier la taille, l'étendue et la projection de la cellule de sortie. Recherchez un menu "Options" Spatial Analyst ou un menu "Environnement" ArcToolBox dans ArcGIS.

— whuber

@Matt: J'ai peut-être raté quelque chose dans la description de Tom Patterson. (C'est difficile à suivre parce qu'il dit seulement quels boutons appuyer, pas pourquoi ni ce qu'ils font vraiment.) Pourriez-vous nous expliquer un peu ce que le "curseur d'opacité / transparence" fait pour nous?

— whuber