Kriging by Indication
This window presents the Kriging by Indication step which aims at obtaining a regular grid of values acquired from point sample data. This module was implemented in SPRING based on the indicator kriging using the "ik3d" code from GSLIB -
Geostatistical Software Library ( Deutsch and Journel, 1992). This module relies on
measures
of spatial variability according to regular grids, of spatial attributes, continuous or categorical. Along with an attributes map, an uncertainty map is generated with an estimation of a regular grid representation.
Executing Kriging by Indication:
- Select on the "Control Panel" a numeric IL that contains point samples;
- on the SPRING main menu choose Analysis - Geostatistics - Kriging by Indication...; Two windows will be presented: Model Parameters and Kriging by Indication. The window Model Parameters can also be activated by pressing the button Models/Probability... inside the window Kriging by Indication. Follow the instructions described at Models/Probability... to fill out the fields;
- on top of the window "Kriging by Indication" you will see the variable correspondent to the active IL;
- the structural parameters can be changed at any time by clicking the button Models/Probability...;
- define the attribute nature by choosing Continuous or Categorical on the field Variable;
- choose between Ordinary or Simple at the field Kriging Type;
- Define an option for Kriging by Indication: Full IK or Median. If the Option is Median, type the median value at the field Threshold;
- the button Indirect Data... must be activated in case you want to choose a file with auxiliary data of attributes that are correlated with the attribute being studied. These indirect data must be in a file with Geo-EAS ( Deutsch e Journel, 1992 ) format, X and Y coordinates on the columns 1 and 2 respectively, and indirect attribute values already transformed by indication, placed on the column 3.
- Press the Bounding Box... button in case you want to define a new area for the grid to be generated. The default value is the active Information Layer;
- The fields X Res. and Y Res. contain predefined values for the active IL. Change these values in case you want to use different resolution. Notice that the values defined for the bounding box and resolution determine the size of the grid and can not be greater than 1000 rows or 1000 columns;
- The values Minimum and Maximum at the number of points in the searching area correspond to the minimum and maximum number of points at the searching ellipsoid. The default values are 4 and 16 but you can change them if you wish.
- Next, define the radius and orientation for the Searching Ellipsoid. The fields R.min, R.max and Angle have predefined values for the Isotropic case: R.min and R.max (in meters) are equivalent to the reach of the isotropic variogram and Angle is equal to zero. If there is anisotropy, the parameters must be adjusted accordingly ( Deutsch e Journel, 1992);
- As outputs of Kriging by Indication we have two Information Layers with regular grid representations, one with attribute values and another with estimation uncertainties. To define these ILs it is necessary to choose a Category, type in a name for the IL, IL values (different from the active IL), choose the type of value and uncertainty. For continuous attributes you should choose for Value mean or median and for Uncertainty confidence intervals based on standard deviations or quartiles. For categorical attributes, the distribution mode is chosen for Value and the maximum probability or Shannon entropy for Uncertainty . The uncertainty Information Layer has a _Inc suffix added to the IL values field.
- Click on Apply.
[Procedures]
[Concepts] [Exploratory
Analysis ] [
Semivariogram Analysis
]
[Semivariogram Modeling] [Model
Validation ] [Kriging]
See also :
SPRING - Spatial Analysis
Theoretical Models for adjusting experimental variogram
|