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:
fuente
Respuestas:
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.
fuente
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:
Puede ejecutarlo en la ventana de ArcMap Python configurando
in_raster
yfc
así:donde
r_extent
es una clase de entidad poligonal existente. Luego simplemente copie el código y ejecútelo. Entiendo esto:fuente
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?
fuente
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.
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.
fuente
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.
fuente
Creo que, en el mundo de ArcGIS, puede lograrlo con la herramienta 3D Analyst \ Conversion \ From Raster \ Raster Domain .
fuente
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í.
fuente
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.
fuente