Estrarre elevazioni DEM raster in tabelle di attributi vettoriali in ArcGIS Desktop?

Ho un progetto ArcMap 10 con più gruppi di livelli. Un livello è un DEM raster, un altro gruppo di livelli è un gruppo di shapefile vettoriali (inclusi punti, linee e poligoni). Il client desidera che tutti i dati vettoriali contengano valori Z.

Il popolamento delle tabelle degli attributi per i dati vettoriali con i valori X e Y è abbastanza semplice ma come posso popolare i valori Z?

Esiste un modo in cui posso utilizzare Field Calculator per "afferrare" i valori di altezza dal livello DEM per ciascuna funzione vettoriale?

Dovrei aggiungere che non abbiamo una licenza per le estensioni Analista 3D o Analista spaziale.

In ArcGIS, l'ho mai fatto solo con i vettori punti: non sono sicuro di come lo faresti con i poligoni o le funzioni di linea. Purtroppo, però, hai bisogno dell'analista spaziale. C'è uno strumento in Strumenti Analista spaziale \ Estrazione \ Estrai valori multipli in punti ed Estrai valori in punti . Estrai più valori in punti è meglio perché ti consente di rinominare il nome del campo z in cui stai estraendo i tuoi valori: Estrai valori in punti non ti consente di farlo. Ho problemi a farlo in ArcMap (lo strumento si arresta in modo anomalo), ma funziona bene in ArcCatalog.

L' estensione (gratuita) di Hawths Analysis Tools per ArcMap aveva uno strumento Intersect Point ma lo sviluppo dell'estensione è stato interrotto un paio di anni fa. È stato sostituito da Geospatial Modeling Environment , che non ho mai usato.

SAGA-GIS (software standalone freeware) ha due funzioni in Forme: Griglia \ Aggiungi valori griglia a punti e forme - Griglia \ Aggiungi valori griglia a Forme che possono probabilmente fare quello che vuoi. Questa è la descrizione dello strumento Aggiungi valori griglia alle forme:

Join spaziale: recupera le informazioni dalle griglie selezionate nelle posizioni delle forme del livello forme selezionato e le aggiunge al livello forme risultante. Per i punti è simile al modulo "Aggiungi valori griglia ai punti". Per linee e poligoni, i valori medi verranno calcolati dalle celle della griglia interferenti. Per i poligoni il modulo "Grid Statistics for Polygons" offre opzioni più avanzate.

Secondo me, SAGA è la soluzione migliore.

Spero che sia di aiuto!

Non ho provato questo per i poligoni, ma funziona davvero bene per i punti. In Arcmap 10 se hai accesso a uno strumento di analisi 3D puoi cercare lo strumento Interpola forma Quindi seleziona il tuo file raster e il tuo livello vettoriale regola i valori z se vuoi elevazioni in piedi e il tuo DEM è in metri e lascia correre . Una volta completato, quindi cercare lo strumento AGGIUNGI XY e selezionare il layer appena creato e quindi calcolerà i valori X, Y, Z nella tabella degli attributi di tali layer.

forse potresti usare gdallocationinfo per recuperare tutti i valori Z passandogli i valori X, Y. una volta che hai i valori Z usa il calcolatore feild per aggiungere il valore Z ai punti che compongono ciascuna funzione. Ci sarebbero alcuni script coinvolti, ma in teoria dovrebbe funzionare.

A seconda delle unità o della risoluzione richiesta, è sempre possibile esportare il raster DEM in un file di forma poligonale, direttamente dalla casella degli strumenti, quindi eseguire un join spaziale: questo è più semplice per i punti. Se hai bisogno di valori Z per linee o poligoni, avrai ovviamente bisogno di segmenti definiti o suddivisioni di poligoni ad ogni cambiamento di elevazione come riportato dalla griglia DEM vettoriale. Ciò può essere fatto tramite strumenti di geoprocessing standard e join spaziali, ma è necessario prestare molta attenzione ai requisiti di risoluzione e alla quantità di elaborazione che consumerà un set di dati di grandi dimensioni.

Consiglio di ottenere una licenza per analista 3D e di utilizzare Aggiungi informazioni sulla superficie . In alternativa, utilizzare Interpolate Shape seguito da Add Z Information, che consente ad alcune funzionalità di restituire valori NA senza arrestare l'intero programma. Altrimenti, ci sono due metodi gratuiti comprovati che richiedono più legwork. È possibile implementarli all'interno di QGIS , in modo da non dover lasciare i vostri comodi ambienti GIS.

Lo strumento ESRI che vorresti utilizzare è "Aggiungi informazioni Z". Forse puoi ottenere una prova e vedere se lo strumento sarebbe di ulteriore utilità per te. Se ottenere l'estensione 3D Analyst non è un'opzione, l'unica cosa che mi viene in mente (oltre a una soluzione di programmazione) è estrarre i tuoi dati dall'ESRI e utilizzare una delle opzioni gratuite come SAGA-GIS / qGIS la loro "Aggiungi griglia Value to Shapes "o" Grid Statistics for Polygons "non è flessibile come ESRI ma fa la stessa cosa se stai cercando l'altezza media lungo il poligono / linea.

So come fare questo per punti. Linee e poligoni sono diversi perché attraversano più celle della griglia del DEM. Puoi fare punti lungo le linee o qualcosa del genere ed estrarre valori per quelli. In Arc, avrai bisogno dell'analista spaziale per ottenere i valori, ma puoi ottenerli tramite Qgis (gratuito), anche se non altrettanto utile. Dovrai essere abile nel lavorare con i file Excel e convertirli avanti e indietro dai file di forma. Qgis genererà un file di forma che ha solo le elevazioni dei punti estratti, quindi dovrai copiarlo e incollarlo nella tabella originale come file Excel (con i valori X, Y, in modo da poterlo ripristinare in un file di forma con "visualizza i dati x, y"). Dovrai tenere traccia delle tue proiezioni per assicurarti di scegliere quella giusta durante la riconversione in shapefile.