DEM to GRASS raster
Cliff Behrens
cliff at bellcore.com
Fri Feb 25 09:29:10 EST 1994
Try this...
1) dd if=diskfilename ibs=10240 | \
m.dmaUSGSread top=1 bottom=400 left=1 right=500 \
output=dem logfile=log
Of course you need to substitute the right parms for your area.
For example, bottom and right might be replaced by 1201 for all
the data. (I left the blocking parm as above.)
2) m.rot90 input=dem output=dem.r90
rows=400 cols=500 bpc=2
3) Create a region (in UTM coordinates) with g.region, then
r.in.ll input=dem.r90 output=dem.utm
dimension=500,400 bpc=2 res=3,3
corner=sw,37:13N,103:45W spheroid=wgs72
in the command above you have to give lat-long for a corner of your
dem using the wgs72 spheroid. The r.in.ll command should create all
the necessary support files for you. But if you run into problems,
just run r.support.
Best of luck,
Cliff Behrens
More information about the grass-user
mailing list