Tenere conto dei passaggi a ponte in un DEM di terra per consentire il flusso d'acqua sotto il cavalcavia

Ho un DEM di terra nuda ricavato da dati lidar (A).

Funziona bene e ha superato tutti i test di precisione. È utilizzato nei modelli di mareggiate. Il problema che stiamo riscontrando è sui sottopassi. Devo consentire all'acqua di fluire attraverso i sottopassi quando una strada ne attraversa un'altra. Il mio pensiero è di:

1) Ottieni tutte le strade contrassegnate come ponti da OSM.

2) Seleziona tutte le strade trasversali che attraversano sotto di esse (si intersecano).

3) Quindi, sulla nuova selezione di nuovi segmenti stradali (2), convertili in raster.

4) Estrarre tutti i valori DEM (A) in questi segmenti di strada.

5) Seleziona il minimo di questi valori nel nuovo segmento di strada (2).

6) Bruciare questo valore minimo nel DEM Devo tagliare un canale attraverso il sottopasso.

Aperto a tutti i suggerimenti, miglioramenti, idee su come raggiungere questo obiettivo, comprese soluzioni completamente diverse.

Altre opzioni sono solo ottenere il punto di intersezione e ampliarlo o forse identificare i lavandini e tornare indietro da questo.

Martin ha ragione nel dire che mentre il tuo flusso di lavoro andrà bene per un caso utente specifico, non tiene conto di molti dei problemi che gli argini stradali creano per la modellazione del percorso di flusso utilizzando dati LiDAR a risoluzione fine, come i problemi con il flusso discontinuo nei fossati lungo la strada e gli effetti di canali sotterranei non mappati minori (che possono alterare notevolmente i percorsi di flusso). Raccomando due approcci per affrontare l'apparente danneggiamento dei percorsi di flusso risultante dagli argini stradali (ovvero l'incapacità di rappresentare adeguatamente canali sotterranei e ponti in DEM raster). Ho implementato entrambi gli approcci negli strumenti di analisi geospaziale GIS Whitebox gratuiti e open source. Innanzitutto, se disponi di un buon set di dati su strade e stream, puoi utilizzare Burn Streams at Roadsstrumento per imporre il flusso da un lato dell'argine all'altro attraverso la masterizzazione del flusso limitato. (Si noti che in genere non sono un fan della masterizzazione di stream full-out per imporre percorsi di flusso quando si utilizzano i dati LiDAR perché le reti di stream mappate sono raramente accurate come i flussi derivati ​​da LiDAR, cioè il flusso mappato a volte cade dal flusso digitale derivato da DEM .) Questa è solo una soluzione adeguata quando i tuoi livelli di strade e stream sono di alta qualità ed è spesso il caso che i dati LiDAR rappresentino flussi migliori delle reti della linea blu derivate dai dati delle mappe.

Il secondo approccio è più inclusivo ed è quello di applicare lo strumento Breach Depressions al DEM. L'implementazione della violazione della depressione in Whitebox utilizza un approccio di percorso meno costoso per determinare il percorso di violazione. Questo metodo è particolarmente adatto per rinforzare il flusso lungo i fossati lungo la strada e per identificare le posizioni in cui canali sotterranei o ponti tagliano gli argini stradali (vedi illustrazione sotto). Ciò può fornire una soluzione molto più adatta anche in assenza di dati accessori (ad esempio una rete stradale o di rete di alta qualità). Risolverà anche il problema comune che gli argini stradali causano quando si utilizza il riempimento della depressione per modellare i percorsi di flusso superficiale, ovvero le grandi dighe di manufatti che spesso compaiono all'interno delle valli fluviali dietro gli argini.

Ed ecco un esempio di quel dannato artefatto che può verificarsi dietro gli argini stradali quando si usa il riempimento della depressione:

AGGIORNARE

Ho appena rilasciato l'ultima versione del software, che include lo strumento Burn Streams At Roads descritto in questo post.

Sembra un flusso di lavoro valido per questa attività. Tuttavia, a seconda del livello di precisione con cui lavori, ti perderai cose come le tubature sotto le strade che consentono a fossati e piccole insenature di fluire attraverso la strada. Ciò richiederebbe probabilmente qualche lavoro sul campo per identificare però.

Non menzionate quale software utilizzate, ma per ArcMap point (6) assicuratevi di usare Calcola le statistiche delle celle piuttosto che Mosaic raster , poiché non interpola inutilmente.

Buona fortuna e per favore condividi le tue esperienze qui dopo averle provate!

Ok ecco l'immagine

Blu sono i ponti di OSM. Il DEM è in scala di grigi con edifici (voglio edifici)

Ho selezionato le strade che si intersecano che attraversano i ponti.

Memorizza le strade selezionate in base alla larghezza approssimativa di una strada.

Ho eseguito "Minimo zonale" sul buffer.

Ora hai il raster minimo zonale (rosso) con i valori bassi.

Ho appena usato CON per fondere i due raster (originale e minimo zonale), quindi ho il DEM originale con i nuovi valori più bassi nell'area rossa.

Sembra funzionare alla grande.

Ora da fare in batch.

Il lotto ha funzionato e il processo sembra modificare il nostro modello di alluvione nel modo previsto. Le nostre uscite del nuovo modello ora si allineano molto meglio con i segni d'acqua misurati dalla FEMA nelle aree dei sottopassi.

Ho pensato che l'apertura fosse totale per i ponti, che penso sia un presupposto sicuro.