Allineare due griglie raster equi-risoluzione non coincidenti in ArcGIS Desktop?

Ho due raster: un modello digitale di elevazione (DEM) e una banda Landsat. La risoluzione della griglia di entrambi è di 30 m, ma le griglie non sono coincidenti (vedi immagine: il dark raster in cima è il DEM; il light raster sotto è la banda di Landsat).

Voglio eseguire qualche trasformazione della banda Landsat in modo che la sua griglia si allinei perfettamente con la griglia DEM. Ho provato a ricampionare la banda Landsat a 30 m in ArcGIS usando DEM come Snap Raster nelle impostazioni dell'ambiente, ma questo non ha alcun effetto. Ho anche provato a tagliare la banda di Landsat nella misura di un poligono vettoriale (a cui anche il DEM è troncato), impostando il DEM come Snap Raster, ma senza risultati.

Ho il sospetto che ciò accada perché entrambi i raster sono già alla stessa risoluzione (solo su griglie diverse), ma ho davvero bisogno che le loro griglie siano allineate per poter fare ulteriori analisi.

Posso farlo in ArcGIS Desktop?

Dopo ulteriori ricerche, ho scoperto che l'impostazione di Snap Raster non funzionante era effettivamente causata da un bug in ArcGIS 10. La soluzione sembra essere quella di disattivare il Geoprocessing in background ( Geoprocessing \ Geoprocessing Options ... quindi deselezionare Abilita nella sezione Elaborazione in background ), ed eseguire l'operazione di ritaglio, con Snap Raster specificato, in ArcCatalog.

C'è una discussione sui forum ArcGIS qui: http://forums.arcgis.com/threads/50808-Another-Snap-to-Raster-problem

Apparentemente il bug è stato corretto in ArcGIS SP 4. Ho installato il service pack ma non ho ancora testato.

È interessante (o meglio, fastidiosamente), se eseguo la clip sul mio poligono dell'area di interesse con lo Snap Raster specificato (Snap Raster è DEM) come per la soluzione nel post del forum, sopra, lo Snap Raster funziona ma il Landsat la band ha una colonna in più e una riga in più rispetto al raster DEM, che è stato anche ritagliato in ArcGIS nel poligono dell'area di interesse. Questo è fastidioso perché per ulteriori analisi ho bisogno che i due raster abbiano esattamente lo stesso numero di righe e colonne. Una semplice soluzione consiste nell'eseguire un calcolo con la Con()funzione che crea un nuovo raster con lo stesso numero di righe e colonne del DEM, ma che ha i valori della banda Landsat.

In realtà ho lottato con questo problema per molto tempo e ho trovato quello che considero un rimedio molto efficace e abbastanza facile. Scegli uno dei tuoi raster, quello che ritieni sia nel posto "corretto", ingrandisci fino a quando ArcGIS lo consentirà e crea un nuovo file di forma punti direttamente sopra l'intersezione di quattro celle raster. L'unico trucco qui è quello di essere sicuri di aver ingrandito il più possibile e di creare il nuovo punto direttamente sopra la "croce" formata da quattro celle raster. Quindi registra le coordinate del punto appena creato. Infine, riproiettare il secondo raster (ma lasciare invariate le dimensioni della cella di proiezione e di output) e utilizzare le coordinate del punto appena creato come "Punto di riferimento". Questo allineerà perfettamente i tuoi due raster. Nessun errore,

Potresti provare a usare lo strumento di mascheramento per sovrapporre il DEM a Landsat.

Potrebbe essere necessario creare una casella dell'area di studio per catturare pixel che si sovrappongono l'uno sull'altro.

Misurare lo spazio tra un raster e l'altro, quindi spostare di conseguenza un raster.

È possibile utilizzare lo strumento spostamento dalla casella degli strumenti di gestione dei dati.

Ho avuto lo stesso problema. Ho provato la soluzione Louis Longchamps , che ha funzionato meravigliosamente su un piccolo campione di prova, ma il mio set di dati era troppo grande per usare praticamente questo metodo. Dopo un po 'di tempo e molti modi per provare a farlo, alla fine, per un capriccio, ho usato "Estrai per rettangolo" dopo un ricampionamento, che ha funzionato proprio come avrebbe dovuto funzionare la funzione snap raster nello strumento clip. Ho creato e provato ad usare un modellista, ma mi ha dato un errore, quindi ho dovuto fare questo un passo alla volta, ma alla fine ho trovato tutto il mio raster abbinato.

Non mi sentivo a mio agio con l'idea di spostare il raster perché, secondo me, cambiate le informazioni. Per evitarlo, ho usato una sequenza di conversioni.

Ecco i passaggi:

1. Per il raster che si desidera utilizzare come modello, andare su Properties/Sourcee scrivere Cellsizesu un pezzo di carta con tutti i decimali disponibili.
2. Converti questo raster in punti usando Conversion Tools/Raster to Point.
3. Estrai i valori da tutti gli altri raster che ti interessano Spatial Analyst Tools/Extraction/Extract Multi Values to Points, specificando la funzione Punto appena creata (nel passaggio 2) come funzione del punto di input e tutti i raster che desideri campionare nei raster di input.
4. Crea i raster da questi punti uno per uno usando Conversion Tools/To Raster/Point to Raster, specificando quale livello stai creando nel campo Valore e usando un nome file corrispondente E il Cellsizevalore che hai annotato sul pezzo di carta.

Dovrebbero allinearsi perfettamente e puoi trasferirli su R per esempio e avere matrice e vettori delle stesse identiche dimensioni.

Eravamo così preoccupati per la nostra immagine DEM che non era allineata all'immagine raster classificata supervisionata. Abbiamo anche usato il metodo di "Matthew Long" .

abbiamo aggiunto le istantanee di immagini non allineate e allineate.

Dopo il metodo di Matthew Long, la nostra immagine DEM era totalmente allineata con la nostra immagine raster classificata supervisionata.

Questa è l'immagine dello strumento che abbiamo usato per allineare la nostra immagine DEM, in ARCGIS 10.2.1 in Strumenti di gestione dei dati:

Questa è l'immagine prima dell'allineamento, l'immagine nera è DEM e l'immagine colorata è immagine raster classificata supervisionata.

Questa è l'immagine dopo l'allineamento, l'immagine nera è DEM e l'immagine colorata è un'immagine raster classificata supervisionata.

Risolvo il mio problema con il metodo di "Matthew Long". Tuttavia, per renderlo ancora più semplice e preciso al 100%. Senza alcun ingrandimento, basta usare la coordinata dell'estensione in alto a sinistra del raster allineato bene.

Per fare ciò: tasto destro, proprietà, pannello Sorgente in Proprietà estensione

Inoltre, utilizzando la coordinata dell'estensione in alto a sinistra, è possibile creare una soluzione Python senza interventi manuali.

Ho avuto un problema simile con due file da fonti diverse. Innanzitutto, un file raster per il livello di riferimento della deforestazione a livello nazionale (Perù). In secondo luogo, un file consolidato sulla piattaforma GEE è stato scaricato e riproiettato nel sistema del primo elemento. In conclusione, la sovrapposizione di entrambi mi ha mostrato che c'era uno spostamento tra le celle (30 * 30 celle) non superiore a 5 metri. Al fine di allineare e adattare il secondo file al primo, mi è venuto in mente di stabilire una * base estesa del primo file in grado di coprire entrambi i file * e mantenere l'allineamento dei pixel del primo file con il suo sistema di proiezione e altre funzionalità . Pertanto, al momento di riproiettare il mio secondo file o qualsiasi sottoprodotto con un offset simile,