[GRASSLIST:5990] Re: problem with r.in.ascii

[Schuyler Fishman]

Thanks for your reply, I really appreciate the help everyone has given 
me on this list.
The problem seems to be that when I run s.to.rast, I am unsure of what 
options to
choose, and when I run it with default, I get a segmentation error.

Options - questions

input= site file
output= rater file
size= value
Number of cells to surround site cell
If my model output is a 20m grid and my location resolution is 1m 
should this size be 20?

title= obvious
field= value
Attribute field type to use for operation
options: dim,decimal,cat
default: decimal

I am not sure what to use here, because I imported a floating point 
value and when I do d.info for sites it appears there are two 
categories - the second one is the value I need.

findex= value
Attribute field number to use for operation
default: 1

So should this be 1 or 2?

string= value
String attribute number to use for description
default: 1

Is this the difference between cat and cat description? I am confused.

Results of s.info
----------------------------------------------------------

SITES FILENAME: Brite.89s at PERMANENT
--------------

Header Information:
------------------
         name          Brite.89s
         description   s.in.ascii sites=Brite.89s input=Brite.89s 
fs=space

Number of DIMENSIONS:         2
--------------------
                  - - MIN - -     - - MAX - -
         dim  1  566376.000000   567356.000000   Easting
         dim  2  4179541.000000  4180521.000000  Northing

Type of CATEGORY information: CELL_TYPE
----------------------------
                  - - MIN - -     - - MAX - -
                            1            2500

Number of DOUBLE attributes:  1
---------------------------
                  - - MIN - -     - - Q1 - -      - - MED - -     - - Q3 
- -
         dbl  1             0        0.101395        0.172225        
0.325565
      36.1833

Number of STRING attributes:  0
---------------------------

TOTAL SITES COUNTED: 2500
----------------------------------------------------------


Results of d.s.info
566512.8(E) 4180380(N)
Brite.89s in PERMANENT  566516|4180381 2108  0.0808

566858.2(E) 4180100.6(N)
Brite.89s in PERMANENT  566856|4180101 1425  2.34583
It is the second column that I want.

Thanks in advance for your help

On Wednesday, April 9, 2003, at 03:47  AM, H Bowman wrote:

> Problem is, r.in.ascii doesn't take x,y,z as input. It takes blocks of
> numbers.
>
> For x,y,z data, use s.in.ascii and s.to.rast and see how it goes.
>
>
> Alternatively, get good with Matlab/Octave and convert in into the form
> r.in.ascii does like:
>
> EXAMPLE (from the man page)
> The following is a sample input file to r.in.ascii:
>
> north: 4299000.00
> south: 4247000.00
> east: 528000.00
> west: 500000.00
> rows: 10
> cols: 15
> null: -9999
>
> 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
> 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
> 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
> 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
> 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
> 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
> 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
> 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
> 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
> 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
>
>
>
>> Ok - now I revise my question, the output is x,y,z and the grid
>> spacing (20x20m) is of lower resolution than the location (1x1m).
>> How do I create a raster layer from this data? Does anyone have
>> any quick suggestions or scripts to get the data into the right
>> format?
>>
>> The grid is 500 rows x 500 columbs
>
> [...]
>
>> Here is a sample header from the file
>>
>> south: 4179541
>> east: 567376
>> west: 566376
>> rows: 2500
>> cols: 3
>>
>> 566376.00000 4179541.00000 0.02640
>> 566396.00000 4179541.00000 0.02788
>> 566416.00000 4179541.00000 0.02962
>
>
Schuyler Fishman
Schuylerfish at earthlink.net
www.mindlikesky.com


-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: text/enriched
Size: 5009 bytes
Desc: not available
Url : http://lists.osgeo.org/pipermail/grass-user/attachments/20030409/e1e2b2b1/attachment.bin