提取水体 (Image Analyst)

摘要

使用合成孔径雷达 (SAR) 数据和 DEM 查找水体。

该工具使用输入雷达后向散射来确定哪些像素应被归类为水体,然后为水域创建面。 该工具还会为非水域创建面,并将其视为陆地区域。

插图

提取水体工具图示

使用情况

  • 可以使用应用辐射定标工具将输入激光雷达数据校准至 gamma-nought。 这将优化描绘和分类过程,尤其在大激光雷达场景中。

  • 雷达输入数据可能与输出要素类不一致。 如果输入雷达数据未正射校正,则该工具会使用 DEM 栅格参数转换输出要素类参数值。 如果未提供 DEM,则该工具将使用椭球面执行转换。 要对输出要素类参数值进行最佳转换,请为 DEM 栅格参数提供输入 DEM。 输入的 DEM 必须位于 WGS 84 (EPSG:4326) 地理坐标系中。

  • 此工具还可用于为其他 SAR 工具输入创建水体面或陆地面,例如检测亮色海洋对象检测暗色海洋区域

参数

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

输入雷达数据。

Raster Dataset; Raster Layer
输出要素类

输出面要素类将展示水体和陆地面。

Feature Class
最小面积
(可选)

提取为水体的最小面积。 默认值为 50,000 平方米。

Areal Unit
DEM 栅格
(可选)

输入 DEM。

如果输入雷达数据未经过正射校正,则将使用此 DEM 对其进行正射校正。

此 DEM 还将用于优化面构建。

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

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

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

Boolean

ExtractWater(in_radar_data, out_feature_class, {min_area}, {in_dem_raster}, {geoid})
名称说明数据类型
in_radar_data

输入雷达数据。

Raster Dataset; Raster Layer
out_feature_class

输出面要素类将展示水体和陆地面。

Feature Class
min_area
(可选)

提取为水体的最小面积。 默认值为 50,000 平方米。

Areal Unit
in_dem_raster
(可选)

输入 DEM。

如果输入雷达数据未经过正射校正,则将使用此 DEM 对其进行正射校正。

此 DEM 还将用于优化面构建。

Mosaic Layer; Raster Layer
geoid
(可选)

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

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

代码示例

ExtractWater 示例 1(Python 窗口)

此示例提取面积超过 1 平方公里的水体面。

import arcpy
arcpy.env.workspace = r"C:\Data\SAR"

arcpy.ia.ExtractWater("LowNoise_manifest_CalG0_TNR.crf", 
    "LowNoise_manifest_CalG0_TNR_Water.shp", "1 SquareKilometer", 
    "dem_COP30_ortho.tif", "GEOID")
ExtractWater 示例 2(独立脚本)

此示例提取面积超过 1 平方公里的水体面。

# 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:\SAR\Low Noise_manifest_CalG0_TNR.crf"
out_radar=r"C:\SAR\Low Noise_manifest_CalG0_TNR_Water.shp"
min_area="1 SquareKilometer"
in_dem_raster=r"C:\DEM\dem_COP30_ortho.tif"
geoid="GEOID"

# Execute 
arcpy.ia.ExtractWater(in_radar, out_radar, min_area, in_dem_raster, geoid)

相关主题