I'm reprojecting and changing the format of a vrt to GeoTIFF and can't seem to control metadata options with gdalwarp.
Other tools like gdal_translate have the -mo option to specify metadata options but gdalwarp doesn't.
In this case, the original .dem (CDED format) have a metadata option that says
AREA_OR_POINT = POINT and the output GeoTIFF, using gdalwarp, defaults to
AREA_OR_POINT = AREA.
To be honest, I don't know the significance or importance of the
AREA_OR_POINT metadata option, but just like matching input to output as much as I can.
The command I'm using is:
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" \