Interpolate From Spatiotemporal Points (Image Analyst)

Available with Image Analyst license.

Summary

Interpolates temporal point data into a multidimensional raster.

Usage

  • The points are grouped by time using the Temporal Aggregation parameter.

  • All points within a temporal period are used for interpolation to create a slice in the output multidimensional raster.

  • The result of this tool is a multidimensional raster that contains the variable interpolated across time and space.

  • The tool supports the following interpolation methods:

    • IDW—The output is calculated from the average of points weighted by their inverse distance to the power of 2. The search radius is 20 times the output cell size.
    • Triangulation—The output is calculated using a linear interpolation from a triangulated irregular network constructed from the input points.
    • Mean—The output is calculated using the average of input points.
    • Median—The output is calculated using the median of input points.
    • Natural neighbor—The output is calculated using a natural neighbor interpolation from a triangulated irregular network constructed from the input points.
    • Nearest neighbor—The output is calculated using the nearest point from the center of the output cell.
    • Quadratic—The output is calculated using a quadratic function:

      Quadratic equation

      Where x,y are the location of the measurement's points, t is the time of measurements. The function is fit to all measurements within an output cell and the coefficients are estimated using an iterative least squares approach. The output value for each cell is estimated from the quadratic function given the center coordinate of the output cell and dates.

Parameters

LabelExplanationData Type
Input Temporal Dataset

The input point layer, trajectory layer, or trajectory dataset.

Trajectory Layer; Feature Layer; Mosaic Dataset; Mosaic Layer
Variable Field

A field containing variable values.

String
Time Field

A field containing time values.

String
Temporal Aggregation
(Optional)

Specifies the temporal aggregation of the output multidimensional raster. The interpolation algorithm uses all available data within these time periods to calculate the output slice.

  • DailyThe data values will be aggregated into daily time steps. This is the default.
  • WeeklyThe data values will be aggregated into weekly time steps.
  • MonthlyThe data values will be aggregated into monthly time steps.
  • QuarterlyThe data values will be aggregated into quarterly time steps.
  • YearlyThe data values will be aggregated into yearly time steps.
String
Cellsize
(Optional)

The output cell size. By default, the cell size will be the shorter of the width or the height of the input point feature extent, divided by 250.

Double
Interpolation Method

Specifies the interpolation method that will be used.

  • IDWInverse distance weighted interpolation will be used.
  • TriangulationTriangulation interpolation will be used.
  • MeanMean interpolation will be used.
  • Median Median interpolation will be used.
  • Natural NeighborNatural neighbor interpolation will be used.
  • Nearest NeighborNearest neighbor interpolation will be used.
  • QuadraticQuadratic interpolation will be used.
String

Return Value

LabelExplanationData Type
Output Raster

The output multidimensional raster dataset.

Raster Dataset

InterpolateFromSpatiotemporalPoints(in_dataset, variable_field, time_field, {temporal_aggregation}, {cell_size}, interpolation_method)
NameExplanationData Type
in_dataset

The input point layer, trajectory layer, or trajectory dataset.

Trajectory Layer; Feature Layer; Mosaic Dataset; Mosaic Layer
variable_field

A field containing variable values.

String
time_field

A field containing time values.

String
temporal_aggregation
(Optional)

Specifies the temporal aggregation of the output multidimensional raster. The interpolation algorithm uses all available data within these time periods to calculate the output slice.

  • DAILYThe data values will be aggregated into daily time steps. This is the default.
  • WEEKLYThe data values will be aggregated into weekly time steps.
  • MONTHLYThe data values will be aggregated into monthly time steps.
  • QUARTERLYThe data values will be aggregated into quarterly time steps.
  • YEARLYThe data values will be aggregated into yearly time steps.
String
cell_size
(Optional)

The output cell size. By default, the cell size will be the shorter of the width or the height of the input point feature extent, divided by 250.

Double
interpolation_method

Specifies the interpolation method that will be used.

  • IDWInverse distance weighted interpolation will be used.
  • TRIANGULATIONTriangulation interpolation will be used.
  • MEANMean interpolation will be used.
  • MEDIAN Median interpolation will be used.
  • NATURAL_NEIGHBORNatural neighbor interpolation will be used.
  • NEAREST_NEIGHBORNearest neighbor interpolation will be used.
  • QUADRATICQuadratic interpolation will be used.
String

Return Value

NameExplanationData Type
out_raster

The output multidimensional raster dataset.

Raster Dataset

Code sample

InterpolateFromSpatioTemporalPoints example 1 (stand-alone script)

This example interpolates elevation point data into a multidimensional raster.

# Import system modules
import arcpy
from arcpy import *

#Check out the ArcGIS Image Analyst extension license
arcpy.CheckOutExtension("ImageAnalyst")

# Set local variables
in_dataset = "icesat_trajectory"
out_raster = r"C:\temp\icesat_surface.crf"
variable_field = "elevation"
time_field = "Time"
temporal_aggregation = "Daily"
cell_size = 5000
interpolation_method = "Quadratic"

# Execute
interpolation_output = arcpy.ia.InterpolateFromSpatiotemporalPoints(in_dataset, out_raster, variable_field,
		       time_field, temporal_aggregation, cell_size, interpolation_method)

Related topics