[GRASSLIST:822] Re: Problems importing a GMT grd file

[Sandro Klippel]

David,

you should not put the '-R' option without specify boundaries. You only
need this to select partial areas from original grid.

try,

grdreformat triangle.grd=0 triangle.bin=1 -V (without '-R')


Sandro.

On Sat, 2003-07-26 at 03:53, David Finlayson wrote:
> I tried the grdreformat line as you suggested, now r.in.bin gives the 
> following error:
> 
> North must be larger than South
> 
> Has the header been corrupted in the conversion?  Does r.in.bin accept 
> scientific (exponential) notation in the header?  The Northing 
> coordinates are given in exponential notiation in triangle.grd as you 
> can see below in my first post.
> 
> This is a UTM Zone 10 grid created from XYZ data with the following 
> script (Note the UTM coordinates and non-geographic projection):
> 
> echo " *** Filter the data and convert to grid *** "
> gmtconvert 02022_utm.xyz -bo > utm.b
> blockmedian utm.b -bi -I10 -R534090/548020/5322430/5346610 -V > median.xyz
> 
> echo " *** Build Triangle Grid *** "
> triangulate median.xyz -I10 -E -Gtriangle.grd -R -JX25i/40i -V -M > 
> lines.dat
> 
> echo " *** Convert Grid to proper type for GRASS *** "
> grdreformat triangle.grd=0 triangle.bin=1 -V -R
> 
> Then in GRASS:
> 
> r.in.bin -fh in=triangle.bin out=shoals.gmt
> 
> Thanks for your help!
> 
> David
> 
> Sandro Klippel wrote:
> 
> > Try use 
> > 
> > 
> >>grdreformat triangle.grd=0 triangle.bin=1
> > 
> > 
> > and then
> > 
> > 
> >>r.in.bin -hf in=triangle.bin out=shoals.gmt
> > 
> > 
> > On Fri, 2003-07-25 at 15:18, David Finlayson wrote:
> > 
> >>I have created a GMT grid file from point data using GMT's triangulate 
> >>command.  When I tried to import the grd with r.in.bin, I received an 
> >>error saying there was an invalid row value. Is there something I am 
> >>doing wrong?
> >>
> >>Here are some details:
> >>
> >>GRASS:~/shoals > gridinfo triangle.grd
> >>triangle.grd: Title:
> >>triangle.grd: Command: triangulate -JX25i/40i median.xyz -I2 -E 
> >>-Gtriangle.grd -R534090/548020/5322430/5346610 -V -M
> >>triangle.grd: Remark:
> >>triangle.grd: Normal node registration used
> >>triangle.grd: grdfile format # 0
> >>triangle.grd: x_min: 534090 x_max: 548020 x_inc: 2 units: user_x_unit 
> >>nx: 6966
> >>triangle.grd: y_min: 5.32243e+06 y_max: 5.34661e+06 y_inc: 2 units: 
> >>user_y_unit ny: 12091
> >>triangle.grd: z_min: -11.945 z_max: 48.91 units: user_z_unit
> >>triangle.grd: scale_factor: 1 add_offset: 0
> >>
> >>GRASS:~/shoals > r.in.bin -hf in=triangle.grd out=shoals.gmt
> >>Illegal row value
> >>
> >>I am using GRASS 5.0.2 (April 2003) and GMT 3.4.3 both on Cygwin/XP
> >>
> >>Thanks
> >>
> >>
> >>-- 
> >>David Finlayson
> >>School of Oceanography
> >>Box 357940
> >>University of Washington
> >>Seattle, WA  98195-7940
> >>USA
> >>
> >>Office: Marine Sciences Building, Room 112
> >>Phone: (206) 616-9407
> >>Web: http://students.washington.edu/dfinlays
> 
> -- 
> David Finlayson
> School of Oceanography
> Box 357940
> University of Washington
> Seattle, WA  98195-7940
> USA
> 
> Office: Marine Sciences Building, Room 112
> Phone: (206) 616-9407
> Web: http://students.washington.edu/dfinlays
-- 
Sandro Klippel <sandro.klippel at talha-mar.com.br>

talha-mar projetos ambientais
Conselheiro d'Avila 190 
Porto Alegre,RS 91040-450
+55 (51) 3344-7451

w w w . t a l h a - m a r . c o m . b r