bopen / elevation

Python script to download global terrain digital elevation models, SRTM 30m DEM and SRTM 90m DEM.

Home Page:http://elevation.bopen.eu

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Mostly black resulting image

gingerbeardman opened this issue · comments

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

  • I was expecting a colour image
  • Any ideas what my problem is?

Setup

  • macOS 12.2.1 (21D62)
  • Python 3.9.10 (main, Jan 15 2022, 11:48:04)

I think that's correct for the sample command