Discrete Sibson interpolation

This tool can be used to interpolate a regular grid raster from a ShapeFile of Point ShapeType using a discrete Sibson (natural neighbour) interpolation method. The user inputs one or more ShapeFile datasets and the field within the file's attribute table to be interpolated (if the input file is of a 'Z' type, then the 'Z' value can also be used in place of an attribute). The name of the output raster file must also be specified. The grid cell size and spatial extent of the output grid is determined by specifying either the Cell Size or an existing Base Raster File. One of these two options must be specified. If the user inputs a desired Cell Size, the output grid will have a spatial extent approximately matching that of the input ShapeFile. If the user instead inputs an existing base file raster, the output image will have the same grid cell resolution and extent as the base file. The user must also specify the discretization factor (5-20). The algorithm is extremely computationally intenstive and may take a long while to complete.

Grid cells that have no points within the search radius assigned the NoData value in the output image. The output raster is of the float data type and continuous data scale.

See Also:

Scripting:

The following is an example of a Python script using this tool:

wd = pluginHost.getWorkingDirectory()
# Specify the input file name and attribute # separated by a semicolon inputData = wd + "input.shp";"HEIGHT"
useZValues = "false" # Only used with 'Z' ShapeType
outputFile = wd + "output.shp"
gridRes = "1.5"
baseFile = "not specified"
discretizationFactor = "5"
args = [inputData, useZValues, outputFile, gridRes, baseFile, maxDistance]
pluginHost.runPlugin("DiscreteSibsonInterpolation", args, False)

This is a Groovy script also using this tool:

def wd = pluginHost.getWorkingDirectory()
// Specify the input file name and attribute // separated by a semicolon def inputData = wd + "input.shp";""
def useZValues = "true" // Used with 'Z' ShapeType
def outputFile = wd + "output.shp"
def gridRes = "not specified"
def baseFile = wd + "existingRaster.dep"
def discretizationFactor = "5"
String[] args = [inputData, useZValues, outputFile, gridRes, baseFile, discretizationFactor]
pluginHost.runPlugin("DiscreteSibsonInterpolation", args, false)

Credits: