[GRASS-dev] fseeko() in libiostream

Markus Neteler neteler at osgeo.org
Tue Oct 7 08:34:56 EDT 2008


On Tue, Oct 7, 2008 at 11:53 AM, Paul Kelly
<paul-grass at stjohnspoint.co.uk> wrote:
> Are there any really huge files anybody can test it with?

For easily create a huge DEM file, do this:

SRTM V4 CGIAR Europe:

-------------------- snip fetch_europe.sh -------------
#!/bin/sh

# fetch Europe block
# http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp

for e in `seq 01 07 | awk '{printf "%02d\n", $1}' ` ; do
        for s in `seq 34 45` ; do
                #
http://hypersphere.telascience.org/elevation/cgiar_srtm_v4/tiff/zip/srtm_33_08.ZIP
                wget -c
http://hypersphere.telascience.org/elevation/cgiar_srtm_v4/tiff/zip/srtm_${s}_${e}.ZIP
        done
done
-------------------- snap -------------

# extract:
for i in *.ZIP ; do unzip -o $i ; NAME=`basename $i .ZIP` ; mv
$NAME.TIF $NAME.tif; done

# merge into single LatLong SRTM:
gdalwarp -s_srs epsg:4326 srtm_40_03.tif srtm_40_04.tif srtm_35_03.tif
srtm_40_05.tif srtm_35_04.tif srtm_40_06.tif srtm_35_05.tif
srtm_35_06.tif srtm_41_03.tif srtm_41_04.tif srtm_36_03.tif
srtm_41_05.tif srtm_36_04.tif srtm_41_06.tif srtm_36_05.tif
srtm_36_06.tif srtm_37_03.tif srtm_37_04.tif srtm_37_05.tif
srtm_37_06.tif srtm_38_03.tif srtm_38_04.tif srtm_38_05.tif
srtm_38_06.tif srtm_39_03.tif srtm_39_04.tif srtm_39_05.tif
srtm_39_06.tif europe_south_srtmV4_cgiar_LL.tif

# 2.7GB:
gdalinfo europe_south_srtmV4_cgiar_LL.tif

r.in.gdal ...

hope this helps,
Markus


More information about the grass-dev mailing list