从时空点进行插值 (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

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)

相关主题