[GRASSLIST:3166] grass to gmt (2)

marino vetuschi zuccolini zucco at dipteris.unige.it
Mon Feb 18 03:55:51 EST 2002


Sorry, I wrote a partial procedure.
The good one is

1) r.out.ascii in=filein out=fileout null=NaN

2) sed '1,6 d ; s/ /\
/g' fileout > file_1_column
with this we convert the out ascii in 1-column file

3) xyz2grd file_1_column -Gfileout.grd -Rx/x/x/x -Ix[m|c]/x[m|c] 
-D/=/=/=/=/= -ZTLa -F -V

-R identifies the region, then x/x/x/x are west/east/south/north in 
units as meters, decimal degrees
-I dentifies x_increment and y_increment (the resolution) in meters 
or minutes[m]/seconds[c] if in lat_long of the grid file
-D give the values for 
x_units/y_units/z_units/scale/offset/title/remark. If untouched fix 
to =.
-ZTLa read 1-column file from TopLeft in ASCII format
-F force pixel registration
-V in verbose mode

For instances we can plot the .grd file imposing a vector clip file 
($FILES/new_buff2)  to masquerade a map:

4) psclip $FILES/new_buff2 $frame -Jm$centr_merid/$centr_par/$scala 
-M -K -O -: >> $image_ps
grdcontour $FILES/fileout.grd -Jm$centr_merid/$centr_par/$scala 
-C$sigla.cpt -O -K $V >> $image_ps
psclip -C  -K -O >> $image_ps


marino



More information about the grass-user mailing list