Combinando raster ad alta risoluzione con sfumature generalizzate per produrre un miglior rilievo dell'ombra?

Aggiungo spesso i dati LiDAR a una mappa insieme a un rilievo ombreggiato. Il più delle volte devo generalizzare il DEM LiDAR a seconda della scala della mappa. Tuttavia, ci sono spesso troppi dettagli e le creste e le valli montane non sono ben pronunciate.

Di recente mi sono imbattuto in un articolo che descrive il processo di combinazione di 2 modelli a forma di collina, uno dettagliato e l'altro generalizzato. Il sollievo che ne risulta è molto bello. Vorrei produrre un tale effetto nelle mie mappe: supponendo di avere 2 raster di tipo Hillshade (dettagliati e generalizzati), qual è il modo migliore per combinarli per ottenere tale effetto?

— Jakub Sisak GeoGraphics
fonte

Risposte:

Sembra il lavoro di Tom Patterson su Resolution Bumping GTOPO 30 in Photoshop . La teoria è descritta abbastanza bene da essere adattabile ad altri software, sebbene il lavoro debba essere fatto elaborando le specifiche. L'idea di base è generalizzare (sfocare) un set di dati, molto, per enfatizzare la forma generale e nascondere dettagli specifici e quindi fondere insieme le immagini ad alta risoluzione e lo-res usando i canali alfa o le percentuali di opacità.

Mentre il processo descritto da Patterson è applicato ai modelli di elevazione stessi, l'ho adattato ai rilievi come segue.

Carica l'immagine in rilievo in GIMP:

strato duplicato
livello superiore == livello di dettaglio, lascia in pace, imposta la modalità su Overlay
bottom == il livello della forma, applica la sfocatura gaussiana con l'impostazione 20
salva una copia (relief_blurred.tif)

Riassegnare la proiezione:

gdalcopyproj shaded_relief.tif relief_blurred.tif

C'è molto spazio per l'adattamento nella scelta delle percentuali da utilizzare, se utilizzare i filtri Overlay o Screen o Darken ecc.

— Wilkie opaco
fonte

ottimo metodo.

— Mapperz

ohhh ho intenzione di provare questo, sembra bello.

— Nathan W,

speravo in una soluzione ArcGIS ma ci proverò anche io ...

— Jakub Sisak GeoGraphics,

Sto iniziando a pensare che l'unico modo per farlo sia quello che hai suggerito in un editor di immagini come GIMP, Photoshop, ecc. (Dato che la rampa di colori ArcGIS non può essere impostata da trasparente a nero, ma per impostazione predefinita va sempre dal bianco al nero)

— Jakub Sisak GeoGraphics,

Vedo che "ArcGIS" è un tag, Jakub. Usando Analista spaziale calcoleresti semplicemente una media ponderata delle due sfumature di collina. Ad esempio, la miscela 60-40 potrebbe essere generata con un calcolo come questo:

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

Se ne hai bisogno, la sfocatura gaussiana può essere eseguita eseguendo alcuni mezzi focali circolari su un raster. Il raggio effettivo di fare m mezzi focali, ciascuno del raggio r , è r * Sqrt ( m ).

— whuber
fonte

Grazie Bill. Il mio input dettagliato Hillshade è una dimensione della cella 5, 5 ma la generalizzata Hillshade è 20, 20. Il risultato del calcolatore raster è quindi anche 20, 20. Sembra che esalta le ombre primarie ma i dettagli fini vanno persi. Cosa consiglieresti di conservare la dimensione originale della cellula? Ricampionare la collina generalizzata? Quale tecnica di resapling usare?

— Jakub Sisak GeoGraphics,

L'ho provato, tutto funziona, ma il prodotto finale di Hillshade non è molto diverso dall'altro originale dettagliato che un po 'più liscio ... (Piuttosto che intensificare il dettaglio originale di Hillshade)

— Jakub Sisak GeoGraphics

l'idea centrale ma non espressa nella risposta di Bill è che praticamente tutto ciò che puoi fare in Gimp / Photoshop / ecc. per quanto riguarda la miscelazione di strati e sovrapposizioni, è possibile eseguire calcoli raster, purché si conosca la matematica. Non è facile, ma è per questo che non ci sono molte soluzioni pronte. Seguendo la teoria da sola, non sono una testa di matematica, penso che ciò che manca nell'esempio sopra sia un terzo quadro che funga da cursore per l'opacità / trasparenza di un uomo povero.

— matt wilkie,

@Jakub: hai ragione, devi conservare le dimensioni più sottili delle celle, altrimenti tutto è perduto. Nei prodotti ESRI ciò avviene di solito impostando preliminarmente una variabile "ambiente" da qualche parte per specificare la dimensione, l'estensione e la proiezione della cella di output. Cercare un menu "opzioni" di Analista spaziale o un menu "ambiente" di ArcToolBox in ArcGIS.

— whuber

@Matt: potrei essermi perso qualcosa nella descrizione di Tom Patterson. (È difficile da seguire perché dice solo quali pulsanti premere, non perché o cosa stanno realmente facendo.) Ti dispiacerebbe elaborare un po 'ciò che il "cursore opacità / trasparenza" sta facendo per noi?

— whuber