Copyright © 2019 EarthSoft, Inc • Modified: 12 Nov 2019
The standard method for contouring involves the interpolation of measured values onto a rectangular mesh and contouring of the gridded data. Inverse distance, kriging, and natural neighbor are offered as interpolation methods. Contours are constructed using either the xfarbe routines (Algorithm 671 - FARB-E-2D: Fill Area with Bicubics on Rectangles, ACM Transactions on Mathematical Software, vol. 15, no. 1, March 1989) or a 2D version of the marching cubes algorithm. For contours generated on a vertical profile, the interpolation is performed in three dimensions to obtain the data values on a three-dimensional grid. In plan mode, the interpolation procedure does not consider the elevation or depth of a particular measurement in its interpolation onto the rectangular mesh corner points, although data queries can exclude source data measurements by depth or elevation. For example, to contour only shallow measurements, set the depth or elevation query to select only those measurements.
An alternative triangulation contouring procedure is available for plan view contours. This involves the generation of triangles, with corners at the location of the measured data points. The contours are constructed within each triangle based on linear interpolation and the values of the contoured data at the corner nodes. This method is frequently useful for data reported on a rectangular grid.
•View a training video on contouring using the inverse distance interpolation method in EnviroInsite here.
•View a training video on contouring using the kriging interpolation method in EnviroInsite here.
Confine the area for contouring to the data extent. Go to the Grid tab within the Contour dialog box and select the Draw Grid button or the Draw Polygon button. Draw a rectangle or shape that excludes areas beyond the extent of the data.
For contouring values that vary over several orders of magnitude, select Log Transform on the Options tab of the Contour dialog box. This interpolate on the log of the values rather than the actual values. The effect reduces the spread of the contours away from the peak values.
Use the kriging interpolation option. In many cases, the inverse distance interpolation will create a kind of bubbly effect around the data points. Natural neighbor depends on triangulation and can be problematic contouring to the edge of the data or beyond.
Adjust the kriging interpolation range values to a distance for which you want the measured values to have an impact on the contours. For example, if high interpolation values extend a great distance from a peak measurement location, reduce the range value to limit the influence of the peak measurement on the contours. Reducing the search radius may have the same impact as reducing the range, but may cause the results to be non-smooth.
If the direction of flow is known, you may want to elongate the persistence of interpolated values in the direction of flow. Do this by selecting the kriging interpolation option and adjusting the maximum range. Set the maximum range angle to be aligned with the direction of groundwater flow. The angle is specified in degrees clockwise relative to due north.
Creating a Contour Plot
Click Plot> Data from the main menu and select Contours. Modify the contour properties on the Query tab, Interval tab, Format tab, Options tab, Grid tab, Interpolation tab, Virtual Points tab, and EQuIS Query tab as desired. Click the OK button to save changes.
Cross-Section Contour Option
In contouring on a cross-section, multiple measurement points are assigned along the length of the specified measurement interval. The first value in the Points per Sample Interval is the minimum number of points that are assigned to a sample interval. Then the interval height per measurement point is calculated. If this exceeds the entered Max Interval Height per Point, the number of points to achieve that value is calculated and that number of points is distributed along the interval horizon.
This tab allows the user to define the grid on which the measured values will be interpolated prior to contouring. Increasing the number of cells may improve the faithfulness of the contoured result to the interpolated field, but increases the time to generate the contours.
This tab allows the user to select the interpolation scheme and the parameters of the interpolation method. The correct selection of interpolation parameters is critical to generate contours that accurately reflect the field data and our expectations of how the values vary between the measured data points. The default parameters are frequently adequate, although some improvement can be anticipated through trial and error. There is no single, objectively optimal set of interpolation parameters. Different methods and parameters work best for different data sets.
Method: Select the contouring method from inverse distance, kriging, natural neighbor, and triangulation.
Radius: The cutoff distance establishes the maximum distance from the interpolated point for which a measured value will be used in the interpolation function. Used by both kriging and inverse distance routines.
Default Value for Points Outside Search Radius: This is the default value for points within the grid that are interpolated because they lie outside the search radius or outside the user-defined polygon from the Grid tab.
Use Default Value: Specify whether or not to use the default value. If not selected, then the non-interpolated nodes will not be contoured.
The kriging routines are derived from the kt3d routine of the Geostatistical Software Library (GSLIB) authored by Clayton Deutsch and Andre Journel (www.gslib.com).
Variogram Model: Select from the spherical, exponential, gaussian, power, and hole variogram models using conventions of GSLIB.
Sill: Set the variogram model sill or c value.
Nugget: Set local scale variability.
Range a max: set maximum horizontal range.
Range a min: Set minimum horizontal range.
Range a vert: Set vertical range.
Angle: Rotation angle in degrees clockwise about 12 o'clock of axis defining the direction of maximum horizontal range. In some cases it may be advantageous to increase the a max value and to align this axis with the direction of ground-water flow.
Inverse Distance Parameters
Z Scale: The vertical distances are multiplied by this value prior to calculating the interpolation weight assigned to a measured value. Z Scale values greater than one will result in lower weights assigned to measured values that are not at the same elevation as the grid point.
Exponent: The measured values are assigned weights that are equal to one over the distance between the measured value and the grid point raised to some power. This exponent is the power used in that calculation.
Smooth Distance: The inverse distance method tends to cause a kind of bubbly surface with the bubbles coinciding with measurement points. This type of trend can be reduced by introducing a non-zero smoothing distance. The smoothing distance is effected by adding this distance to the separation distance between the measurement point and the grid point prior to calculating the interpolation weight.
Virtual Points Tab
Virtual points are used to control generating contours with sparse data. In those cases or in cases of water bodies that are hydraulically continuous with groundwater, it may be advantageous to create virtual measurement points that will control the resulting contours.