9 votos

¿Crear un geopackage a partir de varios archivos raster?

Tengo un montón de geotiffs que me gustaría combinar en un único geopackage. Utilizando QGIS, sólo he conseguido exportar una capa como geopackage. He intentado fusionar las capas y crear una capa virtual pero entonces QGIS desordena los colores. Los archivos ráster también están organizados en un conjunto de datos de Arcgis Mosaic (geodatabase de archivos) si eso sirve de ayuda.

Actualmente, hago lo siguiente (que probablemente no es la forma más rápida o eficiente):

  1. Generar un archivo .txt que contenga los nombres de archivo de todos los geotiffs (dir /b /s *.tif > list.txt)
  2. Utilice gdal_merge para fusionar todos los geotiffs de list.txt en un geotiff grande (gdal_merge -pct --optfile list.txt -o new_huge_geotiff.tif)
  3. Utilizar gdal_translate para convertir mi nuevo y enorme geotiff en un geopackage (gdal_translate -of GPKG new_huge_geotiff.tif new_tiny_geopackage.gpkg)
  4. utilizar gdaladdo para construir pirámides en mi nuevo geopackage, que supongo que también se podría hacer en el paso anterior al crear el geopackage. Es que soy muy novato. (gdaladdo -r cubic -oo tile_format=JPEG new_tiny_geopackage.gpkg 2 4 6 8 16 32 64)

El proceso de fusionar todos los geotiffs pequeños en uno enorme lleva mucho tiempo - me gustaría poder añadirlos todos a un geopackage a la vez. He intentado usar --optfile y la lista txt generada pero a gdal_translate no le gusta ese comando.

0 votos

¿Debería haber una tabla raster para cada imagen o todos los geotiffs irían a una tabla agregada?

0 votos

Perdona mi ignorancia, pero supongo que cualquiera de los dos funcionaría siempre que todas las imágenes se combinen en un solo geopackage. Tengo un software que lee geopackage pero no lee geotiff, esa es la razón por la que quiero combinarlos.

0 votos

El lío con los colores puede significar que tienes imágenes con paletas diferentes. ¿Son sus imágenes tal vez mapas topográficos o similares?

14voto

Joe Puntos 16

Las nuevas tablas ráster deben añadirse al GeoPackage existente como subconjuntos de datos, tal y como se documenta en la página del manual del controlador GDAL http://www.gdal.org/drv_geopackage_raster.html . Mira el quinto ejemplo

Adición de un nuevo subconjunto de datos a un GeoPackage existente, y elegir un nombre no predeterminado para la tabla ráster.

gdal_translate -of GPKG new.tif existing.gpkg -co APPEND_SUBDATASET=YES -co RASTER_TABLE=new_table

No hay una interfaz de usuario en QGIS para dar los parámetros adicionales, por lo que debe editar manualmente el comando gdal_translate generado automáticamente. Por ejemplo, si tiene un GeoPackage existente "test.gpkg" y quiere añadir la imagen "Img_Sample.tif" como una nueva tabla raster nombrada como "Img_Sample_2", edite el comando para que se vea así:

enter image description here

Si tienes muchas imágenes que convertir te sugiero que te olvides de QGIS y de la interfaz gráfica de usuario y utilices gdal_translate directamente desde la línea de comandos con un pequeño script.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X