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
Schuyler Fishman
Schuylerfish@earthlink.net
www.mindlikesky.com
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: 3566376.00000 4179541.00000 0.02640
566396.00000 4179541.00000 0.02788
566416.00000 4179541.00000 0.02962
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@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: -99991 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: 3566376.00000 4179541.00000 0.02640
566396.00000 4179541.00000 0.02788
566416.00000 4179541.00000 0.02962
Schuyler Fishman
Schuylerfish@earthlink.net
www.mindlikesky.com