IDW (3D Analyst)

Summary

Interpolates a raster surface from points using an inverse distance weighted (IDW) technique.

Learn more about how IDW works

Usage

  • The output value for a cell using inverse distance weighting (IDW) is limited to the range of the values used to interpolate. Because IDW is a weighted distance average, the average cannot be greater than the highest or less than the lowest input. Therefore, it cannot create ridges or valleys if these extremes have not already been sampled (Watson and Philip 1985).

  • The best results from IDW are obtained when sampling is sufficiently dense with regard to the local variation you are attempting to simulate. If the sampling of input points is sparse or uneven, the results may not sufficiently represent the desired surface (Watson and Philip 1985).

  • The influence of an input point on an interpolated value is isotropic. Since the influence of an input point on an interpolated value is distance related, IDW is not ridge preserving (Philip and Watson 1982).

  • The Output cell size parameter can be defined by a numeric value or obtained from an existing raster dataset. If the cell size hasn’t been explicitly specified as the parameter value, it is derived from the Cell Size environment if it has been specified. If the parameter cell size or the environment cell size have not been specified, but the Snap Raster environment has been set, the cell size of the snap raster is used. If nothing is specified, the cell size is calculated from the shorter of the width or height of the extent divided by 250 in which the extent is in the output coordinate system specified in the environment.

  • If the cell size is specified using a numeric value, the tool will use it directly for the output raster.

    If the cell size is specified using a raster dataset, the parameter will show the path of the raster dataset instead of the cell size value. The cell size of that raster dataset will be used directly in the analysis, provided the spatial reference of the dataset is the same as the output spatial reference. If the spatial reference of the dataset is different than the output spatial reference, it will be projected based on the specified Cell Size Projection Method value.

  • Some input datasets may have several points with the same x,y coordinates. If the values of the points at the common location are the same, they are considered duplicates and have no effect on the output. If the values are different, they are considered coincident points.

    The various interpolation tools may handle this data condition differently. For example, in some cases, the first coincident point encountered is used for the calculation; in other cases, the last point encountered is used. This may cause some locations in the output raster to have different values than what you might expect. The solution is to prepare your data by removing these coincident points. The Collect Events tool in the Spatial Statistics toolbox is useful for identifying any coincident points in your data.

  • The barriers option is used to specify the location of linear features known to interrupt the surface continuity. These features do not have z-values. Cliffs, faults, and embankments are typical examples of barriers. Barriers limit the selected set of the input sample points used to interpolate output z-values to those samples on the same side of the barrier as the current processing cell. Separation by a barrier is determined by line-of-sight analysis between each pair of points. This means that topological separation is not required for two points to be excluded from each other's region of influence. Input sample points that lie exactly on the barrier line will be included in the selected sample set for both sides of the barrier.

  • Barrier features are input as polyline features. IDW only uses the x,y coordinates for the linear feature; therefore, it is not necessary to provide z-values for the left and right sides of the barrier. Any z-values provided will be ignored.

  • Using barriers will significantly extend the processing time.

  • This tool has a limit of approximately 45 million input points. If your input feature class contains more than 45 million points, the tool may fail to create a result. You can avoid this limit by interpolating your study area in several pieces, making sure there is some overlap in the edges, then mosaicking the results to create a single large raster dataset. Alternatively, you can use a terrain dataset to store and visualize points and surfaces comprised of billions of measurement points.

  • The input feature data must contain at least one valid field.

  • For data formats that support Null values, such as file geodatabase feature classes, a Null value will be ignored when used as input.

  • References:

    Philip, G. M., and D. F. Watson. "A Precise Method for Determining Contoured Surfaces." Australian Petroleum Exploration Association Journal 22: 205–212. 1982.

    Watson, D. F., and G. M. Philip. "A Refinement of Inverse Distance Weighted Interpolation." Geoprocessing 2:315–327. 1985.

Parameters

LabelExplanationData Type
Input point features

The input point features containing the z-values to be interpolated into a surface raster.

Feature Layer
Z value field

The field that holds a height or magnitude value for each point.

This can be a numeric field or the Shape field if the input point features contain z-values.

Field
Output raster

The output interpolated surface raster.

It is always a floating-point raster.

Raster Dataset
Output cell size
(Optional)

The cell size of the output raster that will be created.

This parameter can be defined by a numeric value or obtained from an existing raster dataset. If the cell size hasn't been explicitly specified as the parameter value, the environment cell size value will be used if specified; otherwise, additional rules will be used to calculate it from the other inputs. See the usage section for more detail.

Analysis Cell Size
Power
(Optional)

The exponent of distance.

Controls the significance of surrounding points on the interpolated value. A higher power results in less influence from distant points. It can be any real number greater than 0, but the most reasonable results will be obtained using values from 0.5 to 3. The default is 2.

Double
Search radius
(Optional)

Defines which of the input points will be used to interpolate the value for each cell in the output raster.

There are two options: Variable and Fixed. Variable is the default.

  • Variable

    Uses a variable search radius in order to find a specified number of input sample points for the interpolation.

    • Number of points—An integer value specifying the number of nearest input sample points to be used to perform interpolation. The default is 12 points.
    • Maximum distance—Specifies the distance, in map units, by which to limit the search for the nearest input sample points. The default value is the length of the extent's diagonal.
  • Fixed

    Uses a specified fixed distance within which all input points will be used for the interpolation.

    • Distance—Specifies the distance as a radius within which input sample points will be used to perform the interpolation.

      The value of the radius is expressed in map units. The default radius is five times the cell size of the output raster.

    • Minimum number of points—An integer defining the minimum number of points to be used for interpolation. The default value is 0.

      If the required number of points is not found within the specified distance, the search distance will be increased until the specified minimum number of points is found.

      When the search radius needs to be increased it is done so until the Minimum number of points fall within that radius, or the extent of the radius crosses the lower (southern) and/or upper (northern) extent of the output raster. NoData is assigned to all locations that do not satisfy the above condition.

Radius
Input barrier polyline features
(Optional)

Polyline features to be used as a break or limit in searching for the input sample points.

Feature Layer

arcpy.ddd.Idw(in_point_features, z_field, out_raster, {cell_size}, {power}, {search_radius}, {in_barrier_polyline_features})
NameExplanationData Type
in_point_features

The input point features containing the z-values to be interpolated into a surface raster.

Feature Layer
z_field

The field that holds a height or magnitude value for each point.

This can be a numeric field or the Shape field if the input point features contain z-values.

Field
out_raster

The output interpolated surface raster.

It is always a floating-point raster.

Raster Dataset
cell_size
(Optional)

The cell size of the output raster that will be created.

This parameter can be defined by a numeric value or obtained from an existing raster dataset. If the cell size hasn't been explicitly specified as the parameter value, the environment cell size value will be used if specified; otherwise, additional rules will be used to calculate it from the other inputs. See the usage section for more detail.

Analysis Cell Size
power
(Optional)

The exponent of distance.

Controls the significance of surrounding points on the interpolated value. A higher power results in less influence from distant points. It can be any real number greater than 0, but the most reasonable results will be obtained using values from 0.5 to 3. The default is 2.

Double
search_radius
(Optional)

Defines which of the input points will be used to interpolate the value for each cell in the output raster.

There are two ways to specify the searching neighborhood: Variable and Fixed.

Variable uses a variable search radius in order to find a specified number of input sample points for the interpolation. Fixed uses a specified fixed distance within which all input points will be used. Variable is the default.

The syntax for these parameters are:

  • Variable, number_of_points, maximum_distance, where
    • number_of_points—An integer value specifying the number of nearest input sample points to be used to perform interpolation. The default is 12 points.
    • maximum_distance—Specifies the distance, in map units, by which to limit the search for the nearest input sample points. The default value is the length of the extent's diagonal.
  • Fixed, distance, minimum_number_of_points, where
    • distance—Specifies the distance as a radius within which input sample points will be used to perform the interpolation. The value of the radius is expressed in map units. The default radius is five times the cell size of the output raster.
    • minimum_number_of_points—An integer defining the minimum number of points to be used for interpolation. The default value is 0.

      If the required number of points is not found within the specified distance, the search distance will be increased until the specified minimum number of points is found.

      When the search radius needs to be increased it is done so until the minimum_number_of_points fall within that radius, or the extent of the radius crosses the lower (southern) and/or upper (northern) extent of the output raster. NoData is assigned to all locations that do not satisfy the above condition.

Radius
in_barrier_polyline_features
(Optional)

Polyline features to be used as a break or limit in searching for the input sample points.

Feature Layer

Code sample

IDW example 1 (Python window)

This example inputs a point shapefile and interpolates the output surface as a TIFF raster.

import arcpy
from arcpy import env  
env.workspace = "C:/data"
arcpy.Idw_3d("ozone_pts.shp", "ozone", "C:/output/idwout.tif", 2000, 2, 10)
IDW example 2 (stand-alone script)

This example inputs a point shapefile and interpolates the output surface as a Grid raster.

# Name: IDW_3d_Ex_02.py
# Description: Interpolate a series of point features onto a
#    rectangular raster using Inverse Distance Weighting (IDW).
# Requirements: 3D Analyst Extension

# Import system modules
import arcpy
from arcpy import env

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

# Set local variables
inPointFeatures = "ca_ozone_pts.shp"
zField = "ozone"
outRaster = "C:/output/idwout01"
cellSize = 2000.0
power = 2
searchRadius = 150000

# Execute IDW
arcpy.ddd.Idw(inPointFeatures, zField, outRaster, cellSize, 
             power, searchRadius)

Related topics