计算变换误差 (编辑)

摘要

根据要用于空间数据变换的已知控制点之间的输入链接的坐标来计算残差和均方根误差 (RMSE)。

使用情况

  • 变换要素工具相比,优先使用此工具来确定控制点是否适合预期变换,尤其是在转换大量要素之前。

  • 输入链接要素是代表相应控制点之间的变换链接的线。

  • 变换操作可在直角或平面坐标系中执行。 建议使用投影坐标系 (PCS)。 将地理坐标系 (GCS) 与经纬度值一同使用可导致变形或造成计算错误。

  • 误差将通过以下三种变换方法之一进行计算:仿射、相似和投影。 每种方法都至少需要某一数量的变换链接。 有关变换公式在内的详细信息,请参阅转换要素

    • AFFINE - 至少需要 3 个变换链接
    • PROJECTIVE - 至少需要 4 个变换链接
    • SIMILARITY - 至少需要 2 个变换链接
  • 变换结果取决于输入链接的质量。 链接应起始于已知的源位置,结束于相应的目标位置,也就是所谓的控制点。 控制点建立得越好,变换的结果越准确。 链接的开始和结束位置的坐标用于派生用于变换方程的参数,这些参数是源控制点和目标控制点之间的最佳拟合,正如变换要素中所述。 即使使用变换参数来变换实际的源控制点,变换后的位置也不会与目标控制点位置完全匹配。 这就是所谓的残差,每个变换链接都会生成残差。 输入链接残差将写入包含以下字段的指定输出表格:

    • Orig_FID- 输入链接要素 ID
    • X_Source- 源的 x 坐标或来自链接的结束位置
    • Y_Source- 链接的源或起始端位置的 y 坐标
    • X_Destination- 链接的目标或结束端位置的 x 坐标
    • Y_Destination- 目标的 y 坐标或到链接的结束位置
    • Residual_Error- 变换位置的残差

    RMSE 也称为均方根偏差 (RMSD),将在残差的基础上计算得出,用于指示派生变换的常规适宜性。 在处理消息中写出 RMSE 值。 它还是可在脚本或模型工作流中使用的派生输出参数。

    变换要素主题提供了有关残差和 RMSE 计算的详细信息。 您必须根据所掌握的输入要素及控制点的定位精度相关知识来确定可接受的 RMSE 值。 如果 RMSE 值过高,请查看残差,并放弃或替换残差过高的链接。

参数

标注说明数据类型
输入链接要素

空间变换的已知控制点的输入链接要素。

Feature Layer
输出链接表
(可选)

包含输入链接要素 ID 及其残差的输出表。 输入链接残差将写入包含以下字段的指定输出表格:

  • Orig_FID- 输入链接要素 ID
  • X_Source- 源的 x 坐标或来自链接的结束位置
  • Y_Source- 链接的源或起始端位置的 y 坐标
  • X_Destination- 链接的目标或结束端位置的 x 坐标
  • Y_Destination- 目标的 y 坐标或到链接的结束位置
  • Residual_Error- 变换位置的残差
Table
方法
(可选)

指定将用于转换输入要素坐标的变换方法。

  • 仿射变换至少需要 3 个变换链接。 这是默认设置。
  • 投影变换至少需要 4 个变换链接。
  • 相似变换至少需要 2 个变换链接。
String

派生输出

标注说明数据类型
RMSE

报告均方根误差 (RMSE) 值。

Double

arcpy.management.CalculateTransformationErrors(in_link_features, {out_link_table}, {method})
名称说明数据类型
in_link_features

空间变换的已知控制点的输入链接要素。

Feature Layer
out_link_table
(可选)

包含输入链接要素 ID 及其残差的输出表。 输入链接残差将写入包含以下字段的指定输出表格:

  • Orig_FID- 输入链接要素 ID
  • X_Source- 源的 x 坐标或来自链接的结束位置
  • Y_Source- 链接的源或起始端位置的 y 坐标
  • X_Destination- 链接的目标或结束端位置的 x 坐标
  • Y_Destination- 目标的 y 坐标或到链接的结束位置
  • Residual_Error- 变换位置的残差
Table
method
(可选)

指定将用于转换输入要素坐标的变换方法。

  • AFFINE至少需要 3 个变换链接。 这是默认设置。
  • PROJECTIVE至少需要 4 个变换链接。
  • SIMILARITY至少需要 2 个变换链接。
String

派生输出

名称说明数据类型
out_rmse

报告均方根误差 (RMSE) 值。

Double

代码示例

CalculateTransformationErrors 示例 1(Python 窗口)

以下 Python 窗口脚本演示了如何在即时模式下使用 CalculateTransformationErrors 函数。

import arcpy
arcpy.env.workspace = "C:/data/Countries.gdb"
arcpy.edit.CalculateTransformationErrors("control_Links", "output_Table", "AFFINE")
CalculateTransformationErrors 示例 2(独立脚本)

以下独立脚本演示了如何在脚本环境中应用 CalculateTransformationErrors 函数。

import arcpy
import os

# All input data is in country.gdb and output will also go to this gdb
arcpy.env.workspace = os.path.join(os.getcwd(), "country.gdb")

in_links_feats = "link_features"
out_link_table = "output_table"

# Transformation method
method = "SIMILARITY"

result = arcpy.edit.CalculateTransformationErrors(in_links_feats, out_link_table, method)

# Get the transformation error
error = float(result.getOutput(1))

# If error is less than 12.234, run Transform Features
if error < 20.0:
    # Make a copy of the input features 
    arcpy.management.CopyFeatures(in_links_feats, "in_links_copy")
    arcpy.edit.TransformFeatures("in_links_copy", in_links_feats, method, "out_link_table")
else:
    print("Transformation error {} is too high".format(error))