とらりもんHOME  Index  Search  Changes  Login

とらりもん - SRTM Diff

  • Added parts are displayed like this.
  • Deleted parts are displayed like this.

* SRTMGL3: 90 m resolutionresolution (3 arc second)
* SRTMGL1: 30 m resoltutionresoltution (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()