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()
Keyword(s):
References:[トレーニングコース] [とらりもんHOME]