David,
I get the same error on your binary file. But the converted file from
gmt grid work fine. So, the problem is in your 'grdreformat'.
I put the new binary file on:
www.talha-mar.com.br/programas/triangle2.zip
output of grdinfo:
triangle.grd: Title:
triangle.grd: Command: triangulate -JX25i/40i median.xyz -I10 -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: 10 units: user_x_unit
nx: 1394
triangle.grd: y_min: 5.32243e+06 y_max: 5.34661e+06 y_inc: 10 units:
user_y_unit ny: 2419
triangle.grd: z_min: -11.8154 z_max: 47.3911 units: user_z_unit
triangle.grd: scale_factor: 1 add_offset: 0
output of r.info:
Type of Map: raster Number of Categories: 255
Data Type: FCELL
Rows: 2419
Columns: 1394
Total Cells: 3372086
Projection: UTM (zone 0)
N: 5346610 S: 5322430 Res: 9.99586606
E: 548020 W: 534090 Res: 9.9928264
Range of data: min = -11.815351 max = 47.391102
Data Source:
Data Description: generated by r.in.bin
On Sat, 2003-07-26 at 19:09, David Finlayson wrote:
I placed the files on our web server:
http://www.ocean.washington.edu/data/pugetsound/grass/triangle.zip
Thanks for your help with this.
David
Sandro Klippel wrote:
> Perhaps you right. I'm not know cygwin, but it worked pretty well on
> linux (grdreformat-r.in.bin).
>
> I can test it for you. Send me your arquives (both .grd and .bin).
>
> On Sat, 2003-07-26 at 16:40, David Finlayson wrote:
>
>>I removed the -R and got the same result. I wonder if this is a cygwin
>>problem?
>>
>>Triangle.grd works fine in GMT plots. Something is falling appart
>>during either the reformat step or during import.
>>
>>Unfortuantley, I don't have a linux box with grass installed to verify that.
>>
>>David
>>
>>Sandro Klippel wrote:
>>
>>>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
>>
>>--
>>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@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