ラベル | 説明 | データ タイプ |
入力レーダー データ | 入力レーダー データ。 | Raster Dataset; Raster Layer |
出力ラスター | 出力 SAR 指数ラスター。 | Raster Dataset |
インデックス (オプション) | 計算される SAR 指数を指定します。
| String |
偏波バンド (オプション) | 指数計算に使用する偏波バンドを指定します。 このパラメーターは、[入力レーダー データ] パラメーター値がクアッド偏波 SAR データセットで、[指数] パラメーター値が [レーダー植生指数 (RVI)] または [レーダー森林劣化指数 (RFDI)] である場合にのみサポートされます。
| String |
Image Analyst ライセンスで利用できます。
サマリー
レーダー植生指数 (RVI)、レーダー森林劣化指数 (RFDI)、樹冠構造指数 (CSI) など、SAR (合成開口レーダー) データのさまざまな SAR 指数を計算します。
これらの指数で使用される式は、入力レーダー データセットで使用できる偏波によって異なります。
使用法
入力レーダー データには距離単位を使用する必要があります。 データがデシベル (dB) 単位の場合、[SAR 単位の変換 (Convert SAR Units)] ツールを使用してレーダー データを距離単位に変換します。
このツールは、クアッド偏波の入力 SAR データについては、RVI 指数に対して 3 つの式の選択肢を提供します。 RVI には、HH, HV または VV, VH または HH, HV, VH, VV 偏波が必要です。
このツールは、RFDI 指数に対して 2 つの選択肢を提供します。 RFDI には、HH, HV または VV, VH 偏波が必要です。
このツールは、CSI 指数に対して 1 つの選択肢を提供します。 CSI には HH, VV 偏波が必要です。
RFDI と RVI の場合、都市、砂漠、水域は本質的に植生が少ないか、植生が存在しないため、不毛、荒廃、森林破壊された地形に対応する値を持つことがあります。 ほとんどの値は 0 ~ 1 の範囲内に収まりますが、2 回反射散乱または体積散乱による外れ値が生じることがあります。
パラメーター
ComputeSARIndices(in_radar_data, out_raster, {index}, {polarization_bands})
名前 | 説明 | データ タイプ |
in_radar_data | 入力レーダー データ。 | Raster Dataset; Raster Layer |
out_raster | 出力 SAR 指数ラスター。 | Raster Dataset |
index (オプション) | 計算される SAR 指数を指定します。
| String |
polarization_bands (オプション) | 指数計算に使用する偏波バンドを指定します。 このパラメーターは、in_radar_data パラメーター値がクアッド偏波 SAR データセットで、index パラメーター値が RVI または RFDI である場合にのみサポートされます。
| String |
コードのサンプル
この例では、クアッド偏波 SAR データセットから RVI を計算します。
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")
この例では、クアッド偏波 SAR データセットから RVI を計算します。
# 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)