[STATSGRASS] grass, gstat, and 3d maps
Javier
andresjavier.garcia at wanadoo.es
Wed Apr 24 04:53:06 EDT 2002
Hi Edzer;
I've tried to save you a little of work. I send you the sites file, and in
the file "paraEdzer.txt", there are all the posibilities I've tried in the
command file, followed by the corresponding answers.
------
Javier
----- Original Message -----
From: Edzer J. Pebesma <e.pebesma at geog.uu.nl>
To: Javier <andresjavier.garcia at wanadoo.es>
Cc: <statsgrass at grass.itc.it>
Sent: Tuesday, April 23, 2002 9:24 AM
Subject: Re: [STATSGRASS] grass, gstat, and 3d maps
> I wrote grass support in gstat such that also the z-dimension was
> used; if the z-dimension is in column 3 of the site list, you should use:
>
> data(something): 'sites', z=3;
>
> i.e., x- and y-coordinates are assumed to be in column 1 and 2.
>
> It is not unlikely that you are the first one to test this feature, so
> please bear with us. Could you send me the offending command file, the
> sites data (or part of it), and the output gstat gives you?
> --
> Edzer
>
-------------- next part --------------
If I don't specify any coordinate or data, the first set of data is read, but only 2d are considered: data(traz1):"trazas1";
gstat: Linux version 2.3.6 (28 January 2002)
Copyright (C) 1992, 2001 Edzer J. Pebesma
name|trazas1
desc|s.in.ascii sites=trazas1
GRASS site list trazas1: 0 cat, 3 dim, 0 str, 2 dbl.
Warning: gstat/grass: specify z coordinate if you want to use it
data(traz1): gstat/grass: 547 sites read successfully.
trazas1 (GRASS site list)
attribute: col[1] [x:] x_1 : [ 686000, 711000]
n: 547 [y:] y_2 : [ 4.185e+06, 4.205e+06]
sample mean: 1.4757 sample std.: 0.621017
[starting interactive mode]
--------------------------------------------------------------------------------------
Just the same if I specify: data(traz1):"trazas1",v=1;
gstat: Linux version 2.3.6 (28 January 2002)
Copyright (C) 1992, 2001 Edzer J. Pebesma
name|trazas1
desc|s.in.ascii sites=trazas1
GRASS site list trazas1: 0 cat, 3 dim, 0 str, 2 dbl.
Warning: gstat/grass: specify z coordinate if you want to use it
data(traz1): gstat/grass: 547 sites read successfully.
trazas1 (GRASS site list)
attribute: col[1] [x:] x_1 : [ 686000, 711000]
n: 547 [y:] y_2 : [ 4.185e+06, 4.205e+06]
sample mean: 1.4757 sample std.: 0.621017
[starting interactive mode]
press return to continue...
-------------------
In this way, I can specify the second set of data, but still are considered 2d: data(traz1):"trazas1",v=2;
Copyright (C) 1992, 2001 Edzer J. Pebesma
name|trazas1
desc|s.in.ascii sites=trazas1
GRASS site list trazas1: 0 cat, 3 dim, 0 str, 2 dbl.
Warning: gstat/grass: specify z coordinate if you want to use it
data(traz1): gstat/grass: 547 sites read successfully.
trazas1 (GRASS site list)
attribute: col[2] [x:] x_1 : [ 686000, 711000]
n: 547 [y:] y_2 : [ 4.185e+06, 4.205e+06]
sample mean: 67.3692 sample std.: 36.3099
[starting interactive mode]
press return to continue...
------------------------------------------------------------------------------------------
If I try to specify the first set of data in this way: data(traz1):"trazas1",v=4;
gstat: Linux version 2.3.6 (28 January 2002)
Copyright (C) 1992, 2001 Edzer J. Pebesma
name|trazas1
desc|s.in.ascii sites=trazas1
GRASS site list trazas1: 0 cat, 3 dim, 0 str, 2 dbl.
Warning: gstat/grass doubles in site list are ignored, check data v field!
Warning: gstat/grass: specify z coordinate if you want to use it
gstat: value not allowed for: data v field not found in sites list
--------------------------------------------------------------
If I just specify the z column (as z=3), at last z column is read but now no set of data is established: data(traz1):"trazas1",z=3;
gstat: Linux version 2.3.6 (28 January 2002)
Copyright (C) 1992, 2001 Edzer J. Pebesma
name|trazas1
desc|s.in.ascii sites=trazas1
GRASS site list trazas1: 0 cat, 3 dim, 0 str, 2 dbl.
data(traz1): using coordinate column 3 as z-coordinategstat/grass: 547 sites read successfully.
trazas1 (GRASS site list)
attribute: col[1]
n: 547
no v attribute set for data(traz1)
gstat: variable not set: trazas1.cmd
------------------------------------------------------------------------------------
Just the same if I specify: data(traz1):"trazas1",z=3,v=1;
gstat: Linux version 2.3.6 (28 January 2002)
Copyright (C) 1992, 2001 Edzer J. Pebesma
name|trazas1
desc|s.in.ascii sites=trazas1
GRASS site list trazas1: 0 cat, 3 dim, 0 str, 2 dbl.
data(traz1): using coordinate column 3 as z-coordinategstat/grass: 547 sites read successfully.
trazas1 (GRASS site list)
attribute: col[1]
n: 547
no v attribute set for data(traz1)
gstat: variable not set: trazas1.cmd
----------------------------------------------------------------------------------------
A different message when I try: data(traz1):"trazas1",z=3,v=4;
gstat: Linux version 2.3.6 (28 January 2002)
Copyright (C) 1992, 2001 Edzer J. Pebesma
name|trazas1
desc|s.in.ascii sites=trazas1
GRASS site list trazas1: 0 cat, 3 dim, 0 str, 2 dbl.
Warning: gstat/grass doubles in site list are ignored, check data v field!
gstat: value not allowed for: data v field not found in sites list
----------------------------------
I have also tried: data(traz1):"trazas1",x=1,y=2,z=3,v=4;
gstat: Linux version 2.3.6 (28 January 2002)
Copyright (C) 1992, 2001 Edzer J. Pebesma
name|trazas1
desc|s.in.ascii sites=trazas1
GRASS site list trazas1: 0 cat, 3 dim, 0 str, 2 dbl.
Warning: gstat/grass doubles in site list are ignored, check data v field!
gstat: value not allowed for: data v field not found in sites list
---------------------------------------------
And: data(traz1):"trazas1",x=1,y=2,z=3,v=5;
gstat: Linux version 2.3.6 (28 January 2002)
Copyright (C) 1992, 2001 Edzer J. Pebesma
name|trazas1
desc|s.in.ascii sites=trazas1
GRASS site list trazas1: 0 cat, 3 dim, 0 str, 2 dbl.
Warning: gstat/grass doubles in site list are ignored, check data v field!
gstat: value not allowed for: data v field not found in sites list
-------------- next part --------------
A non-text attachment was scrubbed...
Name: trazas1
Type: application/octet-stream
Size: 21524 bytes
Desc: not available
Url : http://lists.osgeo.org/pipermail/grass-stats/attachments/20020424/b7aaf13b/trazas1.obj
More information about the grass-stats
mailing list