[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