とらりもんHOME  Index  Search  Changes  Login

SRTM

  • SRTMGL3: 90 m resolution (3 arc second)
  • SRTMGL1: 30 m resoltution (1 arc second)
for lat in `seq 12 21`; do for lon in 097 098 099 100 101 102 103 104 105 106; do
wget --user=**** --password=***** https://e4ftl01.cr.usgs.gov/MEASURES/SRTMGL1.003/2000.02.11/N${lat}E${lon}.SRTMGL1.num.zip
wget --user=**** --password=***** https://e4ftl01.cr.usgs.gov/MEASURES/SRTMGL1.003/2000.02.11/N${lat}E${lon}.SRTMGL1.hgt.zip
wget --user=**** --password=***** https://e4ftl01.cr.usgs.gov/MEASURES/SRTMGL1.003/2000.02.11/N${lat}E${lon}.SRTMGL1.hgt.zip.xml
done; done

1 arcsec (30 m) products:

http://e4ftl01.cr.usgs.gov/SRTM/SRTMGL1.003/2000.02.11/

wget http://e4ftl01.cr.usgs.gov/MEASURES/SRTMGL1.003/2000.02.11/S01E100.SRTMGL1.hgt.zip

How to read into python

$ unzip N36E140.hgt.zip
$ dd conv=swab if=N36E140.hgt of=N36E140sw.hgt    # byte swap for endian conversion
$ ipython3
import numpy as np
fn="N36E140sw.hgt"
c=open(fn, 'rb').read()
d=struct.unpack('H'*(len(c)//2), c)
d=np.resize(np.array(d), (1201,1201))
ind=np.where(d>40000)
d=d/10.0
d[ind]=np.nan
plt.imshow(d)
plt.show()
Last modified:2019/03/05 16:32:07
Keyword(s):
References:[とらりもんHOME] [トレーニングコース]