Вычислить индексы SAR (Image Analyst)

Доступно с лицензией Image Analyst.

Краткая информация

Вычисляет различные индексы SAR для данных радара с синтезированной апертурой (SAR), такие как радиолокационный индекс растительности (RVI), радиолокационный индекс деградации леса (RFDI) и индекс структуры полога (CSI).

Формулы, используемые для этих индексов, зависят от поляризаций, доступных во входном наборе данных радара.

Использование

  • Входные данные радара должны быть в линейных единицах измерения. Если данные указаны в децибелах (дБ), используйте инструмент Преобразовать единицы измерения SAR для преобразования данных радара в линейные единицы.

  • Для входных данных SAR с четырехкратной поляризацией, инструмент предоставит три возможных варианта формулы для индекса RVI. RVI требует поляризации HH, HV; VV, VH или HH, HV, VH, VV.

    Уравнения RVI

    Инструмент предоставит два возможных варианта индекса RFDI. RFDI требует поляризации HH, HV или VV, VH.

    Уравнения RFDI

    Инструмент предоставляет один вариант для индекса CSI. CSI требует поляризации HH, VV.

    Уравнение CSI

  • Для RFDI и RVI города, пустыни и водоемы могут иметь значения, соответствующие бесплодным, деградирующим и обезлесенным ландшафтам, поскольку этим регионам присуща слабая растительность или вообще ее отсутствие. Хотя большинство значений будут находиться в диапазоне от 0 до 1, могут быть выбросы, возникающие в результате двойного отражения или объемного рассеяния.

Параметры

ПодписьОписаниеТип данных
Входные радарные данные

Входные радарные данные.

Raster Dataset; Raster Layer
Выходной растр

Выходной растр индекса SAR.

Raster Dataset
Индекс
(Дополнительный)

Задает вычисляемый индекс SAR.

  • Вегетационный индекс радара (RVI)Будет использован Вегетационный индекс радара. RVI - это отношение обратного рассеяния перекрестной поляризации к общему обратному рассеянию от всех поляризаций. Значения находятся в диапазоне от 0 до 1. Значения RVI, близкие к 0, указывают на бесплодные ландшафты, в то время как большие значения указывают на ландшафты с растительностью. Это значение по умолчанию
  • Индекс деградации леса радара (RFDI)Будет использован Индекс деградации лесов радара. RFDI - это нормализованная разница между обратным рассеянием сополяризованного канала и канала с перекрестной поляризацией. Более низкие значения RFDI (менее 0,3) указывают на густой лес. Умеренные значения RFDI (от 0,4 до 0,6) соответствуют деградирующим лесам. Более высокие значения RFDI (более 0,6) указывают на обезлесенные ландшафты.
  • Индекс структуры крон (CSI)Будет использоваться Индекс структуры крон. CSI - это нормализованная разница сополяризованного обратного рассеяния (HH, VV). Значения находятся в диапазоне от -1 до +1, в котором кроны с преобладанием вертикальной структуры будут иметь значения CSI близкие к -1, в то время как кроны с преобладанием горизонтальной структуры будут иметь значения CSI близкие к 1. Эта опция поддерживается только в том случае, если входные данные радара содержат каналы HH и VV.
String
Каналы поляризации
(Дополнительный)

Задает каналы поляризации, которые будут использоваться при вычислении индекса.

Этот параметр поддерживается только в том случае, если значением параметра Входные данные радара является набор данных SAR с четырехкратной поляризацией, а значением параметра Индекс является Индекс растительности радара (RVI) или Индекс деградации лесов радара (RFDI).

  • HH, HVПри вычислении индекса будут использоваться каналы горизонтальный-горизонтальный и горизонтальный-вертикальный (двойная поляризация). Это значение по умолчанию
  • VV, VHПри вычислении индекса будут использоваться каналы вертикальный-вертикальный и вертикальный-горизонтальный (двойная поляризация).
  • HH, HV, VH, VVПри вычислении индекса будут использоваться каналы горизонтальный-горизонтальный, горизонтальный-вертикальный, вертикальный-горизонтальный и вертикальный-вертикальный (четырехкратная поляризация).
String

ComputeSARIndices(in_radar_data, out_raster, {index}, {polarization_bands})
ИмяОписаниеТип данных
in_radar_data

Входные радарные данные.

Raster Dataset; Raster Layer
out_raster

Выходной растр индекса SAR.

Raster Dataset
index
(Дополнительный)

Задает вычисляемый индекс SAR.

  • RVIБудет использован Вегетационный индекс радара. RVI - это отношение обратного рассеяния перекрестной поляризации к общему обратному рассеянию от всех поляризаций. Значения находятся в диапазоне от 0 до 1. Значения RVI, близкие к 0, указывают на бесплодные ландшафты, в то время как большие значения указывают на ландшафты с растительностью. Это значение по умолчанию
  • RFDIБудет использован Индекс деградации лесов радара. RFDI - это нормализованная разница между обратным рассеянием сополяризованного канала и канала с перекрестной поляризацией. Более низкие значения RFDI (менее 0,3) указывают на густой лес. Умеренные значения RFDI (от 0,4 до 0,6) соответствуют деградирующим лесам. Более высокие значения RFDI (более 0,6) указывают на обезлесенные ландшафты.
  • CSIБудет использоваться Индекс структуры крон. CSI - это нормализованная разница сополяризованного обратного рассеяния (HH, VV). Значения находятся в диапазоне от -1 до +1, в котором кроны с преобладанием вертикальной структуры будут иметь значения CSI близкие к -1, в то время как кроны с преобладанием горизонтальной структуры будут иметь значения CSI близкие к 1. Эта опция поддерживается только в том случае, если входные данные радара содержат каналы HH и VV.
String
polarization_bands
(Дополнительный)

Задает каналы поляризации, которые будут использоваться при вычислении индекса.

Этот параметр поддерживается только в том случае, если значением параметра in_radar_data является набор данных SAR с четырехкратной поляризацией, а значением параметра index является RVI или RFDI.

  • HH_HVПри вычислении индекса будут использоваться каналы горизонтальный-горизонтальный и горизонтальный-вертикальный (двойная поляризация). Это значение по умолчанию
  • VV_VHПри вычислении индекса будут использоваться каналы вертикальный-вертикальный и вертикальный-горизонтальный (двойная поляризация).
  • HH_HV_VH_VVПри вычислении индекса будут использоваться каналы горизонтальный-горизонтальный, горизонтальный-вертикальный, вертикальный-горизонтальный и вертикальный-вертикальный (четырехкратная поляризация).
String

Пример кода

ComputeSARIndices, пример 1 (окно Python)

В этом примере RVI вычисляется на основе набора данных SAR с четырехкратной поляризацией.

import arcpy
arcpy.env.workspace = r"C:\Data\SAR"
    
out = arcpy.ia.ComputeSARIndices(  
    "Quad-Polarization_manifest_CalB0_TNR_RTFG0_Dspk_GTC.crf", "RVI",  
    "HH, HV, VH, VV") 
out.save("Quad-Polarization_manifest_CalB0_TNR_RTFG0_Dspk_GTC_RVI.crf")
ComputeSARIndices, пример 2 (автономный скрипт)

В этом примере RVI вычисляется на основе набора данных SAR с четырехкратной поляризацией.

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

# Set local variables
arcpy.env.workspace = r"C:\Data\SAR" 
in_radar_data = "Quad-Polarization_manifest_CalB0_TNR_RTFG0_Dspk_GTC.crf" 
out_raster = "Quad-Polarization_manifest_CalB0_TNR_RTFG0_Dspk_GTC_RVI.crf"
index = "RVI" 
polarization_bands = "HH, HV, VH, VV" 

# Execute  
out = arcpy.ia.ComputeSARIndices(
            in_radar_data, out_raster, index, polarization_bands) 
out.save(out_raster)

Связанные разделы