[GRASSLIST:3255] KrigingProblem

[stef]

Dear GrassList,
I am trying to run predictions on a site map where I
don't have data. I constructed the sites tables with
attributes which I assigned such as the "%80" on the
first line (i.e. 89.7|50.85|#1 %80) and so on for the
rest of the rows, after deleting two colums of
"%..."(some large numbers created automatically when I
used s.in.ascii). I was able to obtain a nice
variogram from this, but the Kriging will not work as
I expect (i.e. I don't get any predictions of the
attributes of the site locations which do not have
this attributes).

Here is my command file:

# gstat command file, Linux version 2.4.3 
data(VarioRafWork2): 'RafWork', x=1, y=2, v=1,
min=20, max=40, radius=50;
variogram(VarioRafWork2): -48.4928 Nug(0) + 1065.59
Gau(29.992);
data(): 'PredictSites2', x=1, y=2; 
blocksize: dx=20, dy=20;
set output=  'RafWorkTest5.out';

And here is the message upon Kriging:

GRASS:~ > gstat VarioRafWork3.cmd
gstat: Linux version 2.4.3 (04 January 2004)
Copyright (C) 1992, 2004 Edzer J. Pebesma
using Marsaglia's random number generator
Warning: gstat: projection is 0 (is this a problem?)
name|RafWork
desc|s.in.ascii sites=RafWork input=RafWork.csv fs=
GRASS site list RafWork: 0 cat, 2 dim, 0 str, 1 dbl.
data(VarioRafWork2): gstat/grass: 224 sites read
successfully.
 RafWork     (GRASS site list)
attribute:             col[1]     [x:] x_1        : [ 
       0,       390]
n:                        224     [y:] y_2        : [ 
       0,       344]
sample mean:          56.1295     sample std.:        
 30.3768
Warning: gstat: projection is 0 (is this a problem?)
name|PredictSites2
desc|s.in.ascii sites=PredictSites2
input=PredictSites2.csv fs=
GRASS site list PredictSites2: 0 cat, 2 dim, 0 str, 1
dbl.
gstat/grass: 30 sites read successfully.
data():         PredictSites2     (GRASS site list)
attribute:             col[1]     [x:] x_1        : [ 
       0,       390]
n:                         30     [y:] y_2        : [ 
       0,       344]
sample mean:                1     sample std.:        
       0
[using ordinary kriging]
100% done

Thanks,
Stef


	
		
__________________________________
Do you Yahoo!?
Yahoo! Photos: High-quality 4x6 digital prints for 25¢
http://photos.yahoo.com/ph/print_splash