Summary
The SearchNeighborhoodSmooth class can be used to define the search neighborhood for IDW, Local Polynomial Interpolation and Radial Basis Functions (only when the INVERSE_MULTIQUADRIC_FUNCTION keyword is used). The smooth search neighborhood class accepts inputs for a minor and major axis, the size of the search ellipse, the angle of the search ellipse, and a smoothing factor.
Syntax
SearchNeighborhoodSmooth ({majorSemiaxis}, {minorSemiaxis}, {angle}, {smoothFactor})| Parameter | Explanation | Data Type |
majorSemiaxis | The distance, in map units, specifying the length of the major semi axis of the ellipse within which data is selected from. | Double |
minorSemiaxis | The distance, in map units, specifying the length of the minor semi axis of the ellipse within which data is selected from. | Double |
angle | The angle of the search ellipse. | Double |
smoothFactor | Determines how much smoothing will be performed. 0 is no smoothing; 1 is the maximum amount of smoothing. | Double |
Properties
| Property | Explanation | Data Type |
| angle (Read and Write) | The angle of the search ellipse. | Double |
| majorSemiaxis (Read and Write) | The distance, in map units, specifying the length of the major semi axis of the ellipse within which data is selected. | Double |
| minorSemiaxis (Read and Write) | The distance, in map units, specifying the length of the minor semi axis of the ellipse within which data is selected. | Double |
| nbrType (Read Only) | The neighborhood type: Smooth or Standard. | String |
| smoothFactor (Read and Write) | Determines how much smoothing will be performed: 0 is no smoothing, and 1 is the maximum amount of smoothing. | Double |
Code sample
SearchNeighborhoodSmooth with Local Polynomial Interpolation to produce an output raster.
import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.LocalPolynomialInterpolation_ga(
"ca_ozone_pts", "OZONE", "outLPI", "C:/gapyexamples/output/lpiout", "2000",
"2", arcpy.SearchNeighborhoodSmooth(300000, 300000, 0, 0.5), "QUARTIC",
"", "", "", "", "PREDICTION")SearchNeighborhoodSmooth with Local Polynomial Interpolation to produce an output raster.
# Name: LocalPolynomialInterpolation_Example_02.py
# Description: Local Polynomial interpolation fits many polynomials, each
# within specified overlapping neighborhoods.
# Requirements: Geostatistical Analyst Extension
# Import system modules
import arcpy
# Set environment settings
arcpy.env.workspace = "C:/gapyexamples/data"
# Set local variables
inPointFeatures = "ca_ozone_pts.shp"
zField = "ozone"
outLayer = "outLPI"
outRaster = "C:/gapyexamples/output/lpiout"
cellSize = 2000.0
power = 2
kernelFunction = "QUARTIC"
bandwidth = ""
useConNumber = ""
conNumber = ""
weightField = ""
outSurface = "PREDICTION"
# Set variables for search neighborhood
majSemiaxis = 300000
minSemiaxis = 300000
angle = 0
smoothFactor = 0.5
searchNeighbourhood = arcpy.SearchNeighborhoodSmooth(majSemiaxis, minSemiaxis,
angle, smoothFactor)
# Execute LocalPolynomialInterpolation
arcpy.LocalPolynomialInterpolation_ga(inPointFeatures, zField, outLayer, outRaster,
cellSize, power, searchNeighbourhood,
kernelFunction, bandwidth, useConNumber,
conNumber, weightField, outSurface)