最优插值 (Image Analyst)

获得 Image Analyst 许可后可用。

摘要

通过统计合并多个来源的数据生成输出栅格。 该工具可用于将背景数据(如模型输出)与观测点数据(如点测量)合并,进行插值操作。

了解有关“最优插值”工作原理的详细信息。

使用情况

  • 该工具接受背景数据集和观测点数据集作为输入。 输入背景栅格参数值通常是来自模型输出的格网栅格,输入观测点数据参数值通常是包含点观测的要素或轨迹数据集。

  • 要合并数据集,需根据其相对精度为背景数据集和观测点数据集分配权重,以最小化分析误差方差。 相对精度是根据背景误差方差观测点误差方差背景误差相关长度参数值确定的,用于确定背景测量之间的背景误差相关性。

  • 通常情况下,背景误差方差观测点误差方差参数值是从背景和观测点数据估计得到的全局常数。 或者,您可以提供另一个具有相同行和列、像元大小和范围的栅格,以表示背景误差方差。 您还可以提供观测点数据中表示观测点误差方差的字段。

  • 背景误差相关长度参数值用于估计观测点对输出的影响。 它以输入背景数据的空间参考单位表示。 较大的值对远离输出像元的点具有较高的影响力。

  • 假设观测点误差相关性为 0。

参数

标注说明数据类型
输入背景栅格

输入背景栅格,也称为背景字段。

Raster Dataset; Raster Layer; Image Service
输入观测点数据

用于插值的输入点要素。

Feature Layer; Trajectory Layer
观测点字段

包含将用于插值的观测点值的字段。

String
背景误差方差

背景测量的误差方差。

输入可以是单个值或误差方差栅格。 如果提供单个值,该值将用作所有背景测量的误差方差。 如果提供误差方差栅格,则背景数据中的每个像元将从相应的背景误差方差栅格获取其误差方差。 误差方差栅格必须具有与背景数据相同的像元大小和范围。

Double; Raster Dataset; Raster Layer; Image Service
观测点误差方差

观测点的误差方差。 输入可以是单个值或观测点数据中的字段。 如果提供单个值,该值将用作所有观测点的误差方差。 如果提供观测点数据中的字段,则字段中的值将用作每个相应观测点的误差方差。

Double; String
背景误差相关长度
(可选)

背景测量之间的相关长度。 默认值为输入背景栅格参数值的三倍像元大小。

Double

返回值

标注说明数据类型
输出栅格

输出多维栅格数据集。

Raster Dataset

OptimalInterpolation(in_background_raster, in_obs_data, obs_field, background_error_var, obs_error_var, {background_error_corr_length})
名称说明数据类型
in_background_raster

输入背景栅格,也称为背景字段。

Raster Dataset; Raster Layer; Image Service
in_obs_data

用于插值的输入点要素。

Feature Layer; Trajectory Layer
obs_field

包含将用于插值的观测点值的字段。

String
background_error_var

背景测量的误差方差。

输入可以是单个值或误差方差栅格。 如果提供单个值,该值将用作所有背景测量的误差方差。 如果提供误差方差栅格,则背景数据中的每个像元将从相应的背景误差方差栅格获取其误差方差。 误差方差栅格必须具有与背景数据相同的像元大小和范围。

Double; Raster Dataset; Raster Layer; Image Service
obs_error_var

观测点的误差方差。 输入可以是单个值或观测点数据中的字段。 如果提供单个值,该值将用作所有观测点的误差方差。 如果提供观测点数据中的字段,则字段中的值将用作每个相应观测点的误差方差。

Double; String
background_error_corr_length
(可选)

背景测量之间的相关长度。 默认值为 in_background_raster 参数值的三倍像元大小。

Double

返回值

名称说明数据类型
out_raster

输出多维栅格数据集。

Raster Dataset

代码示例

OptimalInterpolation 示例 1(独立脚本)

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

# Import system modules 
import arcpy from arcpy import * 

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

# Set local variables 
in_background_raster = r"C:\data\global_ssh.tif" 
in_obs_data = r"C:\data\jason_3_ssh.shp" 
obs_field = "SSH" 
background_error_var = 0.1 
obs_error_var = 0.2 
background_error_corr_length = 100 # km 

# Execute 
output = arcpy.ia.OptimalInterpolation(in_background_raster, in_obs_data, obs_field, background_error_var, obs_error_var, background_error_corr_length) 
Output.save(r”c:\output\OI_ssh.crf”)

相关主题