¿Crear un shapefile que muestre huellas de rasters?

Tengo alrededor de 1,000 imágenes satelitales en formato tiff, y quiero crear un archivo de formas que sirva como índice para los rásteres. Esto es algo similar a un catálogo de ráster, pero no quiero crear un catálogo de ráster.

Algunos obstáculos que puedo prever es que la imagen está georreferenciada, por lo que no tiene forma rectangular (estoy hablando del área de datos).

Para aclarar, requiero que el polígono cubra solo los píxeles distintos de cero (o no nodatos) del ráster, y no todo el ráster rectangular. La mayoría de las respuestas hasta ahora, dan un polígono rectangular, que cubre los datos, así como los píxeles sin datos.

Mi imagen

Resultado dado por las herramientas que he examinado (como catlog de ráster, varios Arcscripts, script personalizado de Python en una de las respuestas):

Resultado que quiero:

Hay un complemento en QGIS llamado Límite de imagen . Es una gran herramienta Dentro de esta herramienta hay una opción para "Píxeles válidos" que omitirá los bordes negros de una imagen de satélite, por ejemplo.

El siguiente código tomará un ráster de entrada, obtendrá su extensión e insertará esa extensión en una clase de entidad de polígono:

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

Puede ejecutarlo en la ventana de ArcMap Python configurando in_rastery fcasí:

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

donde r_extentes una clase de entidad poligonal existente. Luego simplemente copie el código y ejecútelo. Entiendo esto:

Puede usar gdaltindex para esto: http://www.gdal.org/gdaltindex.html

Sin embargo, seguirá creando rectángulos (por ejemplo, 4 + 1 puntos) en el mismo sistema de referencia que las imágenes. Pero me pregunto si eso realmente es un problema: ¿qué tan grandes son sus imágenes?

Intenté la herramienta "Límite de imagen" en QGIS pero falla con algunos rásteres.

Otra opción que funcionó para mí es: 1. Convertir el ráster de entrada a un único número digital utilizando la utilidad de línea de comandos gdal_translate.

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

Esto creará un ráster con valores de 1 en todas partes donde inp.tif tenía valores distintos de cero.

2. Luego, en ArcMap, puede usar la herramienta "ráster a polígono" para obtener un polígono (shapefile).

El archivo de forma de salida tendrá límites no cuadriláteros, con múltiples vértices poligonales, siguiendo el límite real de la imagen.

Si está interesado en una solución con script, pruebe el comando gdal nearblack . También puede usar la bandera -blanco para recuperar datos delimitados por espacios en blanco. Puede poligonalizar la salida para generar una huella no rectangular de los datos ráster.

##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)

Creo que, en el mundo de ArcGIS, puede lograrlo con la herramienta 3D Analyst \ Conversion \ From Raster \ Raster Domain .

Existe una herramienta en ArcGIS 3D Analyst llamada Raster Domain que crea un polígono o una huella de polilínea de las porciones de datos de un dataset ráster, e ignora las celdas NoData. Funcionó muy bien para mí.

Puede crear un catálogo de ráster con una referencia espacial asignada, agregar el catálogo de ráster a ArcMap y luego "Exportar huella". Recuerde que existe la opción de crear un catálogo de ráster "NO MANEJADO" para que no realice un duplicado de cada dataset ráster.