Moving Window Kriging (Geostatisical Analyst)

License Level:BasicStandardAdvanced

Summary

Recalculates the Range, Nugget, and Partial Sill semivariogram parameters based on a smaller neighborhood, moving through all location points.

Usage

Syntax

GAMovingWindowKriging_ga (in_ga_model_source, in_datasets, in_locations, neighbors_max, out_featureclass, {cell_size}, {out_surface_grid})
ParameterExplanationData Type
in_ga_model_source

The geostatistical model source to be analyzed.

File; Geostatistical Layer
in_datasets
[in_dataset,...]

Input datasets displays two pieces:

  • Dataset—The path and the name of the input data.
  • Field—The name of the required field.
Geostatistical Value Table
in_locations

Point locations where predictions will be performed.

Feature Layer
neighbors_max

Number of neighbors to use in the moving window.

Long
out_featureclass

Feature class storing the results.

Feature Class
cell_size
(Optional)

The cell size at which the output raster will be created.

This value can be explicitly set under Raster Analysis from the Environment Settings. If not set, it is the shorter of the width or the height of the extent of the input point features, in the input spatial reference, divided by 250.

Analysis Cell Size
out_surface_grid
(Optional)

The prediction values in the output feature class are interpolated onto a raster using the Local polynomial interpolation method.

Raster Dataset

Code Sample

MovingWindowKriging example 1 (Python window)

Predict values at select point locations.

import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.GAMovingWindowKriging_ga(
    "C:/gapyexamples/data/kriging.lyr", "C:/gapyexamples/data/ca_ozone_pts.shp OZONE",
    "C:/gapyexamples/data/obs_pts.shp", "10", "C:/gapyexamples/output/outMWK", "", "")
MovingWindowKriging example 2 (stand-alone script)

Predict values at select point locations.

# Name: MovingWindowKriging_Example_02.py
# Description: The kriging model is automatically estimated for each neighborhood
#              as the kriging interpolation moves through all the location points.
# Requirements: Geostatistical Analyst Extension

# Import system modules
import arcpy

# Set environment settings
arcpy.env.workspace = "C:/gapyexamples/data"

# Set local variables
inLayer = "C:/gapyexamples/data/kriging.lyr"
inPoints = "C:/gapyexamples/data/ca_ozone_pts.shp OZONE"
obsPoints = "C:/gapyexamples/data/obs_pts.shp"
maxNeighbors = 10
outPoints = "C:/gapyexamples/output/outMWK"

# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")

# Execute MovingWindowKriging
arcpy.GAMovingWindowKriging_ga(inLayer, inPoints, obsPoints, maxNeighbors,
                               outPoints)

Environments

Related Topics

Licensing Information

ArcGIS for Desktop Basic: Requires Geostatistical Analyst
ArcGIS for Desktop Standard: Requires Geostatistical Analyst
ArcGIS for Desktop Advanced: Requires Geostatistical Analyst
6/24/2013