I work on ETOPO1, a 466MB file, and want to move to SRTM (~20GB Earth). My workflow involve several time-consuming processing steps so I strongly want to work on a minimum raster. It usually want to take 10.000 to 100.000px width GIS rasters to scale them down to ~1000px width, like my final map.
How to correctly resize a raster GIS file to a given width via terminal ? (gdal prefered)
It should keep negative and positive z-values with precision.
Note: I already tried the following…
- Crop: to crop my area of insterest using
gdal_translate -projwin $(WEST) $(NORTH) $(EAST) $(SOUTH) ETOPO1_Ice_g_geotiff.tif crop.tmp.tif
, it still give me files 50 times too big for a processing 50 times slower. - Early resizing: early resizing with Imagemagick + reinjecting the GIS metadata for following processings destroys the negative values (altitudes under sea level become 0m).