Generare contorni da DEM usando GDAL?

16

Abbiamo appena ricevuto una vasta serie di DEM al lavoro e vorrei generare contorni da essi. I DEM hanno una risoluzione di 1m e una dimensione di 1kmx1km.

Uscita da gdalinfo:

Driver: AAIGrid/Arc/Info ASCII Grid
Files: 380000_6888000_1k_1m_DEM_ESRI.asc
Size is 1000, 1000
Coordinate System is `'
Origin = (380000.000000000000000,6888000.000000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Corner Coordinates:
Upper Left  (  380000.000, 6888000.000)
Lower Left  (  380000.000, 6887000.000)
Upper Right (  381000.000, 6888000.000)
Lower Right (  381000.000, 6887000.000)
Center      (  380500.000, 6887500.000)
Band 1 Block=1000x1 Type=Float32, ColorInterp=Undefined
  NoData Value=-9999

So di poter usare gdal_contour per generare i contorni (il mio post sul blog sull'argomento ) ma mi chiedo quali siano le migliori pratiche per generare contorni.

Ci sono delle regole da seguire per ottenere il massimo dai contorni ma non inventare o perdere troppe informazioni?

Supponiamo che io voglia generare tre serie di contorni:

250 millimetri
1m
5m

C'è qualcosa che dovrei fare al DEM prima di ogni set?

Post smoothing delle linee è un buon modo per andare o smoothing raster è un'opzione migliore?

— Nathan W
fonte

1

Per cosa userete i contorni? Quali sono le statistiche sulla posizione e l'accuratezza del valore del DEM?

— whuber

I DEM sono stati generati da LiDAR, quindi direi che sono abbastanza buoni. Beh, vorrei un set per la visualizzazione pubblica per 1m e 5m e 250mm con in house per proposte di design.

— Nathan W,

Per le proposte di design e il lavoro interno probabilmente stai meglio usando lo stesso DEM: rispetto a una serie di polilinee di contorno, si presta meglio a quasi tutti i calcoli. In effetti, per la visualizzazione pubblica puoi spesso fare bene sovrapponendo le funzioni chiave (flussi, creste, picchi, ecc.) Su un DEM con colline. Hai davvero bisogno di generare contorni?

— whuber

Sono personalmente d'accordo con @whuber, ma ho ottenuto il pushback sull'uso di un DEM a causa delle differenze culturali tra i diversi utenti finali. Gli ingegneri (almeno nel mio settore di mercato) sono reticenti a usare qualsiasi cosa tranne i contorni. Inoltre, i contorni possono essere sovrapposti alle immagini aeree più facilmente di un DEM.

— DPierce

14

Regole cartografiche per rappresentare il rilievo durante la presentazione dei contorni famoso libro di Imhof sulla rappresentanza sollievo, capitolo C . Alcune di queste regole sono riportate in questa pagina di Wikipedia . Il consiglio principale quando si semplifica i contorni per preservare le caratteristiche principali del terreno.

La levigatura dei contorni in modo indipendente non impedisce loro di sovrapporsi: si consiglia di levigare prima il DTM. Un livellamento gaussiano tradizionale con un parametro adatto a seconda della risoluzione target consente di cancellare i piccoli dettagli. Uno svantaggio è che il livellamento DTM riempie le valli e le depressioni e appiattisce le creste e le cime. L'uso dell'algoritmo filtro douglass-peucker come in questo documento può essere una soluzione. Esistono anche numerosi metodi basati sull'uso di algoritmi di livellamento del contorno vincolati dalla rete di drenaggio o da uno scheletro. Infine, per evitare che i contorni si sovrappongano in parti sciatte, è possibile cancellarli localmente o anche deformarli .

— julien
fonte

1

Voglio commentare il secondo commento di @ Whuber. L'analisi quantitativa è sempre migliore direttamente da un DEM e l'analisi visiva è spesso migliore se eseguita da una sfumatura piuttosto che dai contorni.

Per rispondere direttamente alla domanda: In ArcGIS vorrei utilizzare Focal Statistics o Aggregate [Spatial Analyst Toolbox] per smussare le linee di contorno risultanti. Poiché i contorni sono una funzione di analisi visiva, la quantità di levigatura varia in base alle proprie esigenze. Quindi dovrai sperimentare e vedere cosa funziona meglio per il tuo progetto.

Il livellamento delle linee dopo averle generate funziona ma è un po 'goffo rispetto alla modifica prima del raster. Una pulizia di generazione post-contorno che faccio spesso è selezionare le linee di una certa lunghezza (ad es. <10 'di lunghezza) ed eliminarle. Questo elimina i dati di piccoli frammenti "rumorosi" di contorni chiusi (cioè piccoli anelli) che difficilmente riflettono il carattere reale della superficie modellata dai contorni.

Altre cose che potresti guardare [utenti ArcGIS] sono:

Realizzare un mosaico raster e / o statistiche focali Raster Mosaic e utilizzare un modello per generare contorni per l'intero set di dati.

Creazione di contorni 3D per AutoCAD.

Fai attenzione alle aree con dati DEM insoliti (ad es. Ampie distese di terreno pianeggiante che necessitano di intervalli di contorno diversi per rappresentare accuratamente la superficie, aree di fitta vegetazione che danno valori di Terra Nuda falsi, aree di rilievo verticale - scogliere, ecc.).

— user23715
fonte

0

C'è un modo semplice usando gdal_contour. Dopo aver impostato tutte le opzioni nella finestra di dialogo è quindi possibile modificare la riga di comando e invece "intervallo -i" è possibile utilizzare livelli fissi "-fl livelli". Come l'immagine mostra muggito. Puoi controllare altre opzioni qui http://www.gdal.org/gdal_contour.html

— christian.gobel
fonte