Table of Contents
This chapter discuss how to run GSTAT from within Matlab.
GSTAT can be run non-interactively using by parsing a parameter file. It is by reading and writing this printer file that mGstat interfaces with GSTAT. Consider the following GSTAT parameter file
# # Local simple point kriging on a mask map # data(ln_zinc): 'zinc.eas', x=1, y=2, v=3, log, min=20, max=40, radius=1000, sk_mean=5.9; variogram(ln_zinc): 0.0554 Nug(0) + 0.581 Sph(900); mask: 'mask_map'; predictions(ln_zinc): 'lzn_skpr'; variances(ln_zinc): 'lzn_skvr';
using read_gstat_par to read the parameter into the data structure G, gives
>> G=read_gstat_par(par) G = mgstat: [1x1 struct] data: {[1x1 struct]} variogram: {[1x1 struct]} mask: {[1x1 struct]} predictions: {[1x1 struct]} variances: {[1x1 struct]} >> G.mgstat ans = parfile: '/home/tmh/RESEARCH/PROGRAMMING/gstat-2.4.4/cmd/ex05.cmd' comment: {'#' '#Local simple point kriging on a mask map' '#'} >> G.data{1} ans = data: 'ln_zinc' file: 'zinc.eas' x: 1 y: 2 v: 3 log: '' min: 20 max: 40 radius: 1000 sk_mean: 5.9000 >> G.predictions{1} ans = data: 'ln_zinc' file: 'lzn_skpr'
The mgstat
field to the G
structure, is specific to mGstat and stores the comments of the parameter file, and the location on disk of the parameter file.
The rest of the fields of G
refers to lines in the GSTAT parameter file.
Take for example the data
field as listed above. The field is called data, because this is the identifier (first string) in the parameter file. Filename and data identifiers are always specified using 'data' and 'file' fields. Several options can be specified for the 'data'. This are simply listed as fields of G.data{1}
. If an option does not supply a value (as the 'log' option) it simply refers to an empty string.
In case several data type are listed in the parameter file, they are read into separate structures, as G.data{1}
, G.data{2}
,...
© 2004-2010 Thomas Mejer Hansen. | This site is hosted by |