GEOSTATS: Variowin - GSTAT conversion??
I am a newcomer to this list, indeed to geostatistics, but I couldn't
find this question in the FAQ. My colleagues and I have several hundred
dates of both potential evapotranspiration and precipitation data as
well as soil profiles that we are attempting to krige. We purchased
several intro-to-geostatistics-texts, the VARIOWIN book/software, and
have created variogram models for each of our datasets. We plan to krige
our data using GSTAT 2.0. However, we're having a couple difficulties
with the conversion. The following is a sample variowin model
The variogram model (for the same data) in VARIOWIN looks like the
Measure of spatial continuity: Variogram
Direction: 50 - Tolerance: 30 - MaxBW.: N/A
Direction: 0 - Tolerance: 30 - MaxBW.: N/A
Gamma(h): 1.058 + 3.404 Gauss.18 (h)
Dir.(1): 50 | anis.(1): 1000
A priori covariance: 4.5763e+00 A priori correlation: 1e+00
MODEL PARAMETERS TO BE INCLUDED IN A GSLIB PARAMETER FILE
1 1.058 :nst, nugget
3 10.392304 3.404 40 1000 :it, aa, cc, ang1, anis
When I translate this model to use in a GSTAT command file, I have used
the following (contents of jun97pet.cmd):
data(jun97pet):'jun97pet.eas', x=1, y=2, v=3;
variogram(jun97pet): 1.058 Nug (0) + 3.404 Gau (18,40,1000);
When the above command file is run in GSTAT, it crashes with the
following report to the screen:
gstat 2.0i (February 1999) for Win32
Copyright (C) 1992, 1999 Edzer J. Pebesma
parsing anis(40,0,0,1000,1), error on 1000
gstat: variable outside valid range: this value should be in <0..1]
Is this because GSTAT is expecting an anisotropy correction as a ratio
between 0..1 and we
are supplying something else? How do we calculate the correct value from
our model (above)?
Also, we think that there are directional issues between the two
softwares VARIOWIN and GSTAT. We understand that Zero degrees in
Variowin is 90 degrees in GSTAT, 5 degrees in Variowin is 85 in GSTAT,
and so on. Are we correct in our assumption?
Is there anything missing from our GSTAT command file that we can
correctly gain from our VARIOWIN model to aid in our kriging, and ought
to use? Thanks in advance for any help you can offer.
GIS & Remote Sensing Laboratory
Woods Hole Research Center
PO Box 296, 13 Church Street
Woods Hole, MA 02543
(508)540-9900, ext 118
*To post a message to the list, send it to ai-geostats@....
*As a general service to list users, please remember to post a summary
of any useful responses to your questions.
*To unsubscribe, send email to majordomo@... with no subject and
"unsubscribe ai-geostats" in the message body.
DO NOT SEND Subscribe/Unsubscribe requests to the list!