Stai creando lo shapefile che mostra le impronte dei raster?

Ho circa 1.000 immagini satellitari in formato tiff e voglio creare un file di forma che fungerà da indice per i raster. Questo è qualcosa di simile a un catalogo raster, ma non voglio creare un catalogo raster.

Alcuni ostacoli che posso prevedere, è che l'immagine è georefrangiata, quindi non ha una forma rettangolare (sto parlando dell'area dati).

Per chiarire, ho bisogno che il poligono copra solo i pixel diversi da zero (o non nodati) del raster e non l'intero raster rettangolare. La maggior parte delle risposte finora fornisce un poligono rettangolare, che copre i dati, nonché i pixel non dati.

La mia immagine

Risultato fornito dagli strumenti che ho esaminato (come catlog raster, vari arcscript, script Python personalizzati forniti in una delle risposte):

Risultato che voglio:

C'è un plugin in QGIS chiamato Image Boundary . E 'un ottimo strumento. All'interno di questo strumento c'è un'opzione per "Pixel validi" che, ad esempio, ometterà i bordi neri di un'immagine satellitare.

Il seguente codice prenderà un input raster, ne prenderà l'estensione e la inserirà in una poligonale poligonale:

import arcpy

r = arcpy.Raster(in_raster)
point = arcpy.Point()
array = arcpy.Array()

corners = ["lowerLeft", "lowerRight", "upperRight", "upperLeft"]

cursor = arcpy.InsertCursor(fc)
feat = cursor.newRow()
for corner in corners:    
    point.X = getattr(r.extent, "%s" % corner).X
    point.Y = getattr(r.extent, "%s" % corner).Y
    array.add(point)
array.add(array.getObject(0))
print len(array)
polygon = arcpy.Polygon(array)
feat.shape = polygon
cursor.insertRow(feat)
array.removeAll()
del feat
del cursor

Puoi eseguirlo nella finestra di ArcMap Python impostando in_rastere in questo fcmodo:

>>> fc = 'r_extent'
>>> in_raster = 'CaliCoast10mNED_HavCurvedPCS'

dove si r_extenttrova una poligonale esistente poligonale. Quindi copia il codice ed eseguilo. Capisco questo:

Puoi usare gdaltindex per questo: http://www.gdal.org/gdaltindex.html

Creerà comunque rettangoli (ad es. 4 + 1 punti) nello stesso sistema di riferimento delle immagini. Ma mi chiedo se questo sia davvero un problema: quanto sono grandi le tue immagini?

Ho provato lo strumento "Image Boundary" in QGIS ma non funziona con alcuni raster.

Un'altra opzione che ha funzionato per me è: - 1. Convertire l'input raster in un singolo numero digitale utilizzando l'utilità della riga di comando gdal_translate.

gdal_translate -co compress = lzw -b 1 -ot byte -scale 1 1 inp.tif out.tif

Questo creerà raster con valori di 1 ovunque inp.tif avesse valori diversi da zero.

2. Quindi in ArcMap puoi usare lo strumento "raster to poligono" per ottenere il poligono (shapefile).

Il file di forma di output avrà confini non quadrilatero, con più vertici poligonali, seguendo il confine effettivo dell'immagine.

Se sei interessato a una soluzione con script, prova il comando gdal nearblack . È possibile utilizzare il flag -white per recuperare anche i dati delimitati da spazi bianchi. È possibile poligonizzare l'output per generare un'impronta non rettangolare dei dati raster.

##run nearblack
shellcmds = "nearblack -white -o " + outpath + " " + inpath
os.system(shellcmds)

##run polygonize
shellcmds = "gdal_polygonize.py " + outpath + " -f" + """ "ESRI Shapefile"  """ + outpath2 + " nearblack"
os.system(shellcmds)

Credo che nel mondo ArcGIS sia possibile ottenerlo con lo strumento Analista 3D \ Conversione \ Da dominio raster \ raster .

Nell'analista 3D ArcGIS è presente uno strumento chiamato Raster Domain che crea un footprint poligonale o polilinea delle porzioni di dati di un set di dati raster e ignora le celle NoData. Ha funzionato alla grande per me.

È possibile creare un catalogo raster con un riferimento spaziale assegnato, aggiungere il catalogo raster ad ArcMap e quindi "Esporta footprint". Ricorda che esiste la possibilità di creare un catalogo raster "UNMANAGED" in modo da non creare un duplicato di ciascun set di dati raster.