Local Polynomial Interpolation

Fits the specified order (zero, first, second, third, and so on) polynomial, each within specified overlapping neighborhoods, to produce an output surface

Interpolation

Contributor(s)
Initial contribute: 2019-05-09

Authorship

  • ESRI
    • :  
    View
    Is authorship not correct? Feed back

    Classification(s)

    ● Method-focused categoriesData-perspectiveGeoinformation analysis

    Detailed Description

    English {{currentDetailLanguage}} English

    Summary

    Fits the specified order (zero, first, second, third, and so on) polynomial, each within specified overlapping neighborhoods, to produce an output surface.

    Usage

    • Use Local Polynomial Interpolation when your dataset exhibits short-range variation.

    • Global Polynomial Interpolation is useful for creating smooth surfaces and identifying long-range trends in the dataset. However, in earth sciences the variable of interest usually has short-range variation in addition to long-range trend. When the dataset exhibits short-range variation, Local Polynomial Interpolation maps can capture the short-range variation.

    Syntax

    LocalPolynomialInterpolation_ga (in_features, z_field, {out_ga_layer}, {out_raster}, {cell_size}, {power}, {search_neighborhood}, {kernel_function}, {use_condition_number}, {bandwidth}, {condition_number}, {weight_field}, {output_type})
    Parameter Explanation Data Type
    in_features

    The input point features containing the z-values to be interpolated.

    		 Feature Layer
    z_field

    Field that holds a height or magnitude value for each point. This can be a numeric field or the Shape field if the input features contain z-values or m-values.

    		 Field
    out_ga_layer
    (Optional)

    The geostatistical layer produced. This layer is required output only if no output raster is requested.

    		 Geostatistical Layer
    out_raster
    (Optional)

    The output raster. This raster is required output only if no output geostatistical layer is requested.

    		 Raster Dataset
    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
    power
    (Optional)

    The order of the polynomial.

    		 Long
    search_neighborhood
    (Optional)

    Defines which surrounding points will be used to control the output. Standard is the default.

    This is a Search Neighborhood class SearchNeighborhoodStandard,SearchNeighborhoodSmooth, SearchNeighborhoodStandardCircular, and SearchNeighborhoodSmoothCircular.

    Standard

    • Major semiaxis—The major semiaxis value of the searching neighborhood.
    • Minor semiaxis—The minor semiaxis value of the searching neighborhood.
    • Angle—The angle of rotation for the axis (circle) or semimajor axis (ellipse) of the moving window.
    • Maximum neighbors—The maximum number of neighbors that will be used to estimate the value at the unknown location.
    • Minimum neighbors—The minimum number of neighbors that will be used to estimate the value at the unknown location.
    • Sector type—The geometry of the neighborhood.
      • One sector—Single ellipse.
      • Four sectors—Ellipse divided into four sectors.
      • Four sectors shifted—Ellipse divided into four sectors and shifted 45 degrees.
      • Eight sectors—Ellipse divided into eight sectors.

    Smooth

    • Major semiaxis—The major semiaxis value of the searching neighborhood.
    • Minor semiaxis—The minor semiaxis value of the searching neighborhood.
    • Angle—The angle of rotation for the axis (circle) or semimajor axis (ellipse) of the moving window.
    • Smoothing factor—The Smooth Interpolation option creates an outer ellipse and an inner ellipse at a distance equal to the Major Semiaxis multiplied by the Smoothing factor. The points that fall outside the smallest ellipse but inside the largest ellipse are weighted using a sigmoidal function with a value between zero and one.

    StandardCircular

    • Radius—The length of the radius of the search circle.
    • Angle—The angle of rotation for the axis (circle) or semimajor axis (ellipse) of the moving window.
    • Maximum neighbors—The maximum number of neighbors that will be used to estimate the value at the unknown location.
    • Minimum neighbors—The minimum number of neighbors that will be used to estimate the value at the unknown location.
    • Sector type—The geometry of the neighborhood.
      • One sector—Single ellipse.
      • Four sectors—Ellipse divided into four sectors.
      • Four sectors shifted—Ellipse divided into four sectors and shifted 45 degrees.
      • Eight sectors—Ellipse divided into eight sectors.

    SmoothCircular

    • Radius—The length of the radius of the search circle.
    • Smoothing factor—The Smooth Interpolation option creates an outer ellipse and an inner ellipse at a distance equal to the Major Semiaxis multiplied by the Smoothing factor. The points that fall outside the smallest ellipse but inside the largest ellipse are weighted using a sigmoidal function with a value between zero and one.
    		 Geostatistical Search Neighborhood
    kernel_function
    (Optional)

    The kernel function used in the simulation.

    • EXPONENTIAL — The function grows or decays proportionally.
    • GAUSSIAN — Bell-shaped function that falls off quickly toward plus/minus infinity.
    • QUARTIC — Fourth-order polynomial function.
    • EPANECHNIKOV — A discontinuous parabolic function.
    • POLYNOMIAL5 — Fifth-order polynomial function.
    • CONSTANT —An indicator function.
    		 String
    use_condition_number
    (Optional)

    Option to control the creation of prediction and prediction standard errors where the predictions are unstable. This option is only available for polynomials of order 1, 2 and 3.

    • NO_USE_CONDITION_NUMBER —Prediction and prediction standard errors can be created where the predictions are unstable. This is the default.
    • USE_CONDITION_NUMBER —Prediction and prediction standard errors will not be created where the predictions are unstable.
    		 Boolean
    bandwidth
    (Optional)

    Used to specify the maximum distance at which data points are used for prediction. With increasing bandwidth, prediction bias increases and prediction variance decreases.

    		 Double
    condition_number
    (Optional)

    Every invertible square matrix has a condition number that indicates how inaccurate the solution to the linear equations can be with a small change in the matrix coefficients (it can be due to imprecise data). If the condition number is large, a small change in the matrix coefficients results in a large change in the solution vector.

    		 Double
    weight_field
    (Optional)

    Used to emphasize an observation. The larger the weight, the more impact it has on the prediction. For coincident observations, assign the largest weight to the most reliable measurement.

    		 Field
    output_type
    (Optional)

    Surface type to store the interpolation results.

    • PREDICTION —Prediction surfaces are produced from the interpolated values.
    • PREDICTION_STANDARD_ERROR —Standard Error surfaces are produced from the standard errors of the interpolated values.
    • CONDITION_NUMBER —The Spatial condition number surface indicates variation in the numerical model stability and provides additional information on the prediction uncertainty as the prediction standard error surface is created, assuming that the model is correct, which implies that there is no reason for model instability.
    		 String

    Code Sample

    LocalPolynomialInterpolation example 1 (Python window)

    Interpolate point features onto a rectangular 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")
    LocalPolynomialInterpolation example 2 (stand-alone script)

    Interpolate point features onto a rectangular 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)


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

# Execute LocalPolynomialInterpolation
arcpy.LocalPolynomialInterpolation_ga(inPointFeatures, zField, outLayer, outRaster,
                                      cellSize, power, searchNeighbourhood,
                                      kernelFunction, bandwidth, useConNumber,
                                      conNumber, weightField, outSurface)
    View Full Text

    模型元数据

    {{htmlJSON.HowtoCite}}

    ESRI (2019). Local Polynomial Interpolation, Model Item, OpenGMS, https://geomodeling.njnu.edu.cn/modelItem/4e1e1733-b16e-4898-9af0-800fed1f0f10
    {{htmlJSON.Copy}}

    Contributor(s)
    Initial contribute : 2019-05-09

    {{htmlJSON.CoContributor}}

    Authorship

  • ESRI
    • :  
    View
    Is authorship not correct? Feed back

    QR Code

    ×
  • {{curRelation.name}}
    {{curRelation.name}}
    {{curRelation.overview}}
    {{curRelation.author.join('; ')}}
    {{curRelation.journal}}





    {{htmlJSON.RelatedItems}}

    {{htmlJSON.LinkResourceFromRepositoryOrCreate}}{{htmlJSON.create}}.

    Drop the file here, orclick to upload.
    Select From My Space
    + add
    {{htmlJSON.Cancel}} {{htmlJSON.OK}}

    {{htmlJSON.authorshipSubmitted}}

    Use My Info Cancel Submit
    Cancel Submit     {{htmlJSON.Cancel}} {{htmlJSON.Submit}}
    {{htmlJSON.Localizations}} + {{htmlJSON.Add}}
    {{ item.label }} {{ item.value }}
    {{htmlJSON.ModelName}}:
    {{htmlJSON.Cancel}} {{htmlJSON.Submit}}     * 名称 别名 {{tag}} +
    系列名 * 版本号 * 目的 * 修改内容 * 创建/修改日期 * 作者
    * 摘要 详细描述
    {{tag}} + 添加关键字
    * 时间参考系
    * 空间参考系类型 * 空间参考系名称

    * 起始日期 终止日期 * 进展 * 开发者
    * 是否开源 * 访问方式 * 使用方式 * 开源协议 * 传输方式 * 获取地址 * 发布日期 * 发布者



    * 编号 * 目的 * 修改内容 * 创建/修改日期 * 作者




    * 时间分辨率 * 时间尺度 * 时间步长 * 时间范围 * 空间维度 * 格网类型 * 空间分辨率 * 空间尺度 * 空间范围
    {{tag}} +     * 类型
    图例

    * 名称 * 描述
    示例描述 * 名称 * 类型 * 值/链接 或 上传


    {{htmlJSON.Cancel}} {{htmlJSON.Submit}}
    Title Author Date Journal Volume(Issue) Pages Links Doi Operation
    {{htmlJSON.Cancel}} {{htmlJSON.Submit}}
    {{htmlJSON.EnterDoiClickSearch}}
    {{htmlJSON.Add}} {{htmlJSON.Cancel}}

    {{articleUploading.title}}

    Authors:  {{articleUploading.authors[0]}}, {{articleUploading.authors[1]}}, {{articleUploading.authors[2]}}, et al.

    Journal:   {{articleUploading.journal}}

    Date:   {{articleUploading.date}}

    Page range:   {{articleUploading.pageRange}}

    Link:   {{articleUploading.link}}

    DOI:   {{articleUploading.doi}}

    Yes, this is it Cancel

    The article {{articleUploading.title}} has been uploaded yet.

    OK     {{htmlJSON.Cancel}} {{htmlJSON.Confirm}}
    9ttHvUSYwmCT