从时空点进行插值 (Image Analyst)

获得 Image Analyst 许可后可用。

摘要

将时间点数据插入到多维栅格中。

使用情况

• 可以使用时间聚合参数，按时间对点进行分组。

• 时间段内的所有点都将用于插值，以在输出多维栅格中创建切片。

• 该工具的结果是一个多维栅格，其中包含跨时间和空间插值的变量。

• 该工具支持以下插值方法：

  • 反距离权重法 - 根据点距离倒数的 2 次方加权平均值来计算输出。 搜索半径是输出像元大小的 20 倍。
  • 三角测量法 - 使用根据输入点构造的不规则三角网络的线性插值来计算输出。
  • 均值 - 使用输入点的均值计算输出。
  • 中值 - 使用输入点的中值计算输出。
  • 自然邻域法 - 使用根据输入点构造的不规则三角网络的自然邻域插值来计算输出。
  • 最邻近法 - 使用距输出像元中心最近的点来计算输出。
  • 二次 - 使用二次函数来计算输出：

    

    其中 x,y 为测量点的位置，t 为测量时间。 此函数适合输出像元内的所有测量，并且使用迭代最小二乘法来估计系数。 根据给定输出像元的中心坐标和日期的二次函数来估计每个像元的输出值。

参数

对话框

标注	说明	数据类型
输入时态数据集	输入点图层、轨迹图层或轨迹数据集。	Trajectory Layer; Feature Layer; Mosaic Dataset; Mosaic Layer
变量字段	包含变量值的字段。	String
时间字段	包含时间值的字段。	String
时间聚合 (可选)	指定输出多维栅格的时间聚合。 插值算法在这些时间段内使用所有可用数据来计算输出切片。 每天—数据值将聚合到每日时间步长。 这是默认设置。 每周—数据值将聚合到每周时间步长。 每月—数据值将聚合到每月时间步长。 季度—数据值将聚合到季度时间步长。 每年—数据值将聚合为年度时间步长。	String
像元大小 (可选)	输出像元大小。 默认情况下，像元大小为输入点要素范围的宽度或高度中的较短值除以 250。	Double
插值方法	指定将使用的插值方法。 反距离权重法—将使用反距离权重插值法。 三角测量法—将使用三角测量插值法。 均值—将使用均值插值法。 中值—将使用中值插值法。 自然邻域法—将使用自然邻域插值法。 最邻近法—将使用最邻近插值法。 二次—将使用二次插值法。	String

返回值

标注	说明	数据类型
输出栅格	输出多维栅格数据集。	Raster Dataset

Python

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

名称	说明	数据类型
in_dataset	输入点图层、轨迹图层或轨迹数据集。	Trajectory Layer; Feature Layer; Mosaic Dataset; Mosaic Layer
variable_field	包含变量值的字段。	String
time_field	包含时间值的字段。	String
temporal_aggregation (可选)	指定输出多维栅格的时间聚合。 插值算法在这些时间段内使用所有可用数据来计算输出切片。 DAILY—数据值将聚合到每日时间步长。 这是默认设置。 WEEKLY—数据值将聚合到每周时间步长。 MONTHLY—数据值将聚合到每月时间步长。 QUARTERLY—数据值将聚合到季度时间步长。 YEARLY—数据值将聚合为年度时间步长。	String
cell_size (可选)	输出像元大小。 默认情况下，像元大小为输入点要素范围的宽度或高度中的较短值除以 250。	Double
interpolation_method	指定将使用的插值方法。 IDW—将使用反距离权重插值法。 TRIANGULATION—将使用三角测量插值法。 MEAN—将使用均值插值法。 MEDIAN—将使用中值插值法。 NATURAL_NEIGHBOR—将使用自然邻域插值法。 NEAREST_NEIGHBOR—将使用最邻近插值法。 QUADRATIC—将使用二次插值法。	String

返回值

名称	说明	数据类型
out_raster	输出多维栅格数据集。	Raster Dataset

代码示例

InterpolateFromSpatioTemporalPoints 示例 1（独立脚本）

此示例将高程点数据插入到多维栅格中。

# 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)

环境

像元大小, 当前工作空间, 范围, 地理变换, 掩膜, 输出坐标系, 并行处理因子, 捕捉栅格