USGS DEMs to GRASS raster maps

In an E-mail message, Bill Brown wrote:

Something like this should work:

dd if=usgsdem-w ibs=10240 | m.dmaUSGSread top=1 bottom=1201 \
left=1 right=1201 output=usgsdem.dat logfile=log

m.rot90 input=usgsdem.dat output=usgsdem.dem rows=1201 cols=1201 bpc=2

r.support
(create header & color files - get header info from log)

You will need to import the rotated file into Grass before you run support,
via the r.in.ll command as shown below:

r.in.ll input=usgsdem.dem output=YourDEM dimension=1201,1201 bpc=2 res=3,3 \
corner=sw,00:00N,00:00W spheroid=wgs72

Subsitute the actual latlon coords for where I showed zeros. Any 1=degree
paper map or USGS map atlas will have the corner coords.

Lastly you can run support. The east and west halves patch up easily
with r.patch.

--
-------------------------

Paul Loechl p-loechl@cecer.army.mil
USA CERL (217)352-6511x7668