¿Cómo configurar las opciones de metadatos para el GeoTIFF de salida cuando se utiliza gdalwarp?

Question

Estoy reproyectando y cambiando el formato de un vrt a GeoTIFF y no puedo controlar las opciones de metadatos con gdalwarp.

Otras herramientas como gdal_translate tienen la opción -mo para especificar opciones de metadatos pero gdalwarp no.

En este caso, los originales .dem (formato CDED) tienen una opción de metadatos que dice AREA_OR_POINT = POINT y el GeoTIFF de salida, utilizando gdalwarp, por defecto es AREA_OR_POINT = AREA .

Para ser honesto, no sé el significado o la importancia de la AREA_OR_POINT opción de metadatos, pero sólo me gusta hacer coincidir la entrada con la salida tanto como pueda.

El comando que estoy usando es:
gdalwarp -of GTiff \ -ot Int16 \ -s_srs "EPSG:4269" -t_srs "EPSG:3005" \ -tr 15 15 -r bilinear \ -wm 4096 -multi \ -co COMPRESS=DEFLATE -co PREDICTOR=2 -co ZLEVEL=9 \ -co NUM_THREADS=ALL_CPUS \ -srcnodata "-32767" \ mosaic_dem_cded_4269.vrt dem_3005.tiff

Answer 1

Puedes deformar a VRT y luego traducir a tu GeoTIFF de salida. Esto también evitará cuestiones que pueden surgir al utilizar gdalwarp para dar salida a los archivos comprimidos.

gdalwarp -of VRT \
          -ot Int16 \
          -s_srs "EPSG:4269" -t_srs "EPSG:3005" \
          -tr 15 15 -r bilinear \
          -wm 4096 -multi \
          -co NUM_THREADS=ALL_CPUS \
          -srcnodata "-32767" \
          mosaic_dem_cded_4269.vrt dem_3005.vrt

gdal_translate -mo "your metadata" dem_3005.vrt dem_3005.tif

También puede canalizar la salida de gdalwarp a gdal_translate utilizando la opción de creación STREAMABLE_OUTPUT=YES y los nombres de archivo virtuales /vsistdout/ y /vsistdin/ (GDAL 2x):

 gdalwarp -of GTiff \
          -ot Int16 \
          -co STREAMABLE_OUTPUT=YES \
          -s_srs "EPSG:4269" -t_srs "EPSG:3005" \
          -tr 15 15 -r bilinear \
          -wm 4096 -multi \
          -co NUM_THREADS=ALL_CPUS \
          -srcnodata "-32767" \
          mosaic_dem_cded_4269.vrt /vsistdout/ | gdal_translate /vsistdin/ dem_3005.tif

Notas: debe utilizar una barra diagonal final para /vsistdout/ y /vsistdin/ (es decir, no /vsistdout y /vsistdin) y he no idea de si el streaming es compatible con el warping multihilo...