应用几何地形校正 (Image Analyst)

获得 Image Analyst 许可后可用。

摘要

使用距离多普勒反向地理编码算法对输入合成孔径雷达 (SAR) 数据进行正射校正。

距离多普勒反向地理编码方法使用轨道状态向量计算每个 DEM 格网点的雷达距离和方位角指数。 如果未提供 DEM,则该工具将使用元数据中包含的连接点来执行距离多普勒地形校正。

使用情况

  • 在运行此工具之前,对输入 SAR 数据运行去斑工具。

  • 该工具可校正 SAR 影像中的几何畸变,准确定位影像并执行定量分析。

  • 该工具不支持将地理数据库用作输出位置。

  • 如果未指定 DEM,则该工具将创建一个近似 DEM,该 DEM 从元数据连接点进行插值。 仅对完全海洋雷达场景使用连接点方法;当雷达场景中包含陆地要素时,必须指定 DEM。

参数

标注说明数据类型
输入雷达数据

输入雷达数据。

Raster Dataset; Raster Layer
输出雷达数据

校正后的几何地形雷达数据。

Raster Dataset
极化波段
(可选)

将进行校正的极化波段。

默认情况下,第一个波段处于选中状态。

String
DEM 栅格
(可选)

输入 DEM。

如果未指定 DEM 或在指定 DEM 未覆盖的区域中,将创建从元数据连接点插值的近似 DEM。

仅对完全海洋雷达场景使用连接点方法;当雷达场景中包含陆地要素时,必须指定 DEM。

Raster Dataset; Raster Layer; Mosaic Layer
应用大地水准面校正
(可选)

指定是否将输入 DEM 的垂直参考系统转换为椭球体高度。 大多数高程数据集均参考海平面正高,因此在这些情况下,需要进行校正以将海平面正高转换为椭球体高度。

  • 选中 - 将进行大地水准面校正以将正高转换为椭球体高度(根据 EGM96 大地水准面)。 这是默认设置。
  • 未选中 - 不会进行大地水准面校正。 仅当以椭球体高度提供 DEM 时,才会使用此选项。

Boolean

ApplyGeometricTerrainCorrection(in_radar_data, out_radar_data, {polarization_bands}, {in_dem_raster}, {geoid})
名称说明数据类型
in_radar_data

输入雷达数据。

Raster Dataset; Raster Layer
out_radar_data

校正后的几何地形雷达数据。

Raster Dataset
polarization_bands
[polarization_bands,...]
(可选)

将进行校正的极化波段。

默认情况下,第一个波段处于选中状态。

String
in_dem_raster
(可选)

输入 DEM。

如果未指定 DEM 或在指定 DEM 未覆盖的区域中,将创建从元数据连接点插值的近似 DEM。

仅对完全海洋雷达场景使用连接点方法;当雷达场景中包含陆地要素时,必须指定 DEM。

Raster Dataset; Raster Layer; Mosaic Layer
geoid
(可选)

指定是否将输入 DEM 的垂直参考系统转换为椭球体高度。 大多数高程数据集均参考海平面正高,因此在这些情况下,需要进行校正以将海平面正高转换为椭球体高度。

  • GEOID将进行大地水准面校正以将正高转换为椭球体高度(根据 EGM96 大地水准面)。 这是默认设置。
  • NONE不会进行大地水准面校正。 仅当以椭球体高度提供 DEM 时,才会使用此选项。
Boolean

代码示例

ApplyGeometricTerrainCorrection 示例 1(Python 窗口)

此示例使用 DEM 对雷达数据进行正射校正。

import arcpy
arcpy.env.workspace = r"C:\Data\SAR"
outRadar = arcpy.ia.ApplyGeometricTerrainCorrection(
     "IW_manifest_TNR_CalB0_Dspk_RTFG0.crf", "VH;VV", r"C:\Data\DEM\dem10m.tif", 
     "GEOID") 
outRadar.save("IW_manifest_TNR_CalB0_Dspk_RTFG0_GTC.crf")
ApplyGeometricTerrainCorrection 示例 2(独立脚本)

此示例使用 DEM 对雷达数据进行正射校正。

# Import system modules and check out ArcGIS Image Analyst extension license
import arcpy
arcpy.CheckOutExtension("ImageAnalyst")
from arcpy.ia import *

# Set local variables
in_radar = r"C:\Data\SAR\IW_manifest_TNR_CalB0_Dspk_RTFG0.crf"
out_radar = r"C:\Data\SAR\IW_manifest_TNR_CalB0_Dspk_RTFG0_GTC.crf"
polarization = "VH"
in_dem_raster = r"C:\Data\DEM\dem10m.tif"
geoid_correction = "NONE"


# Execute 
outRadar = arcpy.ia.ApplyGeometricTerrainCorrection(in_radar, polarization, 
     in_dem_raster, geoid_correction)
outRadar.save(out_radar)

相关主题