bopen / elevation

Python script to download global terrain digital elevation models, SRTM 30m DEM and SRTM 90m DEM.
http://elevation.bopen.eu
Apache License 2.0
290 stars 74 forks source link

Mostly black resulting image #52

Closed gingerbeardman closed 2 years ago

gingerbeardman commented 2 years ago

Steps

$ eio selfcheck
Your system is ready.

$ eio clip -o Rome-30m-DEM.tif --bounds 12.35 41.8 12.65 42
curl -s -o spool/N41/N41E012.hgt.gz.temp https://s3.amazonaws.com/elevation-tiles-prod/skadi/N41/N41E012.hgt.gz && mv spool/N41/N41E012.hgt.gz.temp spool/N41/N41E012.hgt.gz
gunzip spool/N41/N41E012.hgt.gz 2>/dev/null || touch spool/N41/N41E012.hgt
gdal_translate -q -co TILED=YES -co COMPRESS=DEFLATE -co ZLEVEL=9 -co PREDICTOR=2 spool/N41/N41E012.hgt cache/N41/N41E012.tif 2>/dev/null || touch cache/N41/N41E012.tif
rm spool/N41/N41E012.hgt
gdalbuildvrt -q -overwrite SRTM1.vrt cache/N41/N41E012.tif
cp SRTM1.vrt SRTM1.26be21147c87403fbf2dffe59bafacfb.vrt
gdal_translate -q -co TILED=YES -co COMPRESS=DEFLATE -co ZLEVEL=9 -co PREDICTOR=2 -projwin 12.35 42.0 12.65 41.8 SRTM1.26be21147c87403fbf2dffe59bafacfb.vrt /Users/matt/Downloads/Rome-30m-DEM.tif
rm -f SRTM1.26be21147c87403fbf2dffe59bafacfb.vrt

$ open Rome-30m-DEM.tif

Output

Rome-30m-DEM

Expected

Setup

gingerbeardman commented 2 years ago

I think that's correct for the sample command