[Gdal-dev] Re: GDAL SRTM hgt code
Even Rouault
even.rouault at mines-paris.org
Thu Nov 15 14:56:20 EST 2007
Hi,
A specific driver for SRTM hgt files has been added a few weeks ago to GDAL
trunk (future 1.5.0). You may test it with nightly snapshot or latest
FWTools-2.0.0
Le Thursday 15 November 2007 20:14:34 Rick Niles, vous avez écrit :
> Since SRTM hgt-files is the best public available terrain data around
> (that I know of..), it seems like it's worth it to have a specific
> front-end to handle it directly via GDALOpen(). I can see why some
> would argue that's just a special case of a raw 2-byte file, but I don't
> think it would be ugly to have a specific filter to handle this
> extremely popular data. There's 3-second and 1-second SRTM data and I'm
> not sure you can auto-detect which is which.
>
> What do you think?
>
> Rick Niles.
>
> WolfgangZ wrote:
> > list67 at netscape.net schrieb:
> >> Has anyone added code to GDAL so that a SRTM (1 or 3 arcsecond) hgt
> >> file can be passed to GDALOpen(...) for use?
> >>
> >> I did some searching and noticed that there are some format
> >> conversion workarounds for using the SRTM files with GDAL, but I
> >> didn't find any info about wether or not a more direct open
> >> capability had been implemented?
> >>
> >> Thank you.
> >>
> >> Mike
> >
> > hi Mike,
> >
> > you can have a look at the grass r.in.srtm script. It basically unzips
> > and renames the file to .bil and furthermore you need to create a
> > textfile with the same file name and the fiel ending .hdr with the
> > following content:
> > For 3arcsec data:
> > BYTEORDER M
> > LAYOUT BIL
> > NROWS 1201
> > NCOLS 1201
> > NBANDS 1
> > NBITS 16
> > BANDROWBYTES 2402
> > TOTALROWBYTES 2402
> > BANDGAPBYTES 0
> > NODATA -32768
> > ULXMAP $ULXMAP
> > ULYMAP $ULYMAP
> > XDIM 0.000833333333333
> > YDIM 0.000833333333333
> >
> > of for 1arcsec data:
> > BYTEORDER M
> > LAYOUT BIL
> > NROWS 3601
> > NCOLS 3601
> > NBANDS 1
> > NBITS 16
> > BANDROWBYTES 7202
> > TOTALROWBYTES 7202
> > BANDGAPBYTES 0
> > NODATA -32768
> > ULXMAP $ULXMAP
> > ULYMAP $ULYMAP
> > XDIM 0.000277777777777778
> > YDIM 0.000277777777777778
> >
> > You need to replace $ULXMAP and $ULYMAP with the upper left
> > coordinates of that tile. As s_srs you can use EPSG:4326.
> >
> > Greetings
> > Wolfgang
> >
> >
> > _______________________________________________
> > Gdal-dev mailing list
> > Gdal-dev at lists.maptools.org
> > http://lists.maptools.org/mailman/listinfo/gdal-dev
>
> _______________________________________________
> Gdal-dev mailing list
> Gdal-dev at lists.maptools.org
> http://lists.maptools.org/mailman/listinfo/gdal-dev
More information about the gdal-dev
mailing list