Conversion of GTOPO30

[Jacques Bouchard]

Hi Stephan,

> I have to import gtopo30 to my GRASS, but unfortunatly my region falls
> into two tiles of USGS data and the resolution differs:
> 
> g.region:
> 
> CURRENT REGION: N=60S  S=80S  RES=0:06  ROWS=200
>                 E=45E  W=45W  RES=0:20  COLS=270
> PROJECTION: 3 (Latitude-Longitude)
> ZONE:       0
>
> Could your r.in.bin be a solution? If so: I am using an IRIX system,
> so I would have to compile it, or does anybody have the SGI/IRIX binary
> and is willing to send it to me.

Yes r.in.bin is a solution. Here is how you can do:

r.in.bin input=W060S60/W060S60.DEM output=w060s60 north=-60 south=-90
east=0 west=-60 r=3600 c=7200 subst=-9999,0 bytes=2 -s

r.in.bin input=W000S60/W000S60.DEM output=w000s60 north=-60 south=-90
east=60 west=0 r=3600 c=7200 subst=-9999,0 bytes=2 -s

r.mapcalc elevation.dem=w060s60+w000s60

g.remove rast=w060s60,w000s60

You will need some space for the temporary raster maps,
but the final raster map will be quite small:

-rw-r--r--   1 root     root        52757 Jul 27 11:17 elevation.dem
-rw-r--r--   1 root     root     15082536 Jul 27 11:17 w000s60
-rw-r--r--   1 root     root     10727262 Jul 27 11:15 w060s60

The r.in.bin binary I have attached is for SGI IRIX 6.2.

Bye

Jacques
-------------- next part --------------
A non-text attachment was scrubbed...
Name: r.in.bin.gz
Type: application/x-gzip
Size: 42153 bytes
Desc: not available
Url : http://lists.osgeo.org/pipermail/grass-user/attachments/19980727/ba435be4/r.in.bin.gz