コヒーレンスの計算 (Image Analyst)

Image Analyst ライセンスで利用できます。

サマリー

参照レーダー データとセカンダリ入力複素数レーダー データの類似度を計算します。

値の範囲は 0 ~ 1 で、0 はコヒーレンスなし、1 は完全なコヒーレンスを表すコヒーレンス ラスターが出力されます。 0.3 以上の値であれば、コヒーレンス値は良好とみなされます。

使用法

  • 入力は、コレジストレーションされた複素数レーダー データセットのペアでなくてはなりません。

  • 出力コヒーレンス ラスターは、一貫性の変化の検出に使用できます。 出力コヒーレンス ラスターを InSAR ワークフローと DInSAR ワークフローで使用し、インターフェログラム ピクセルの信頼性を評価できます。

    InSAR アプリケーションと DInSAR アプリケーションでは、入力レーダー データの衛星ジオメトリは同じである必要があります。

  • [ジオプロセシング] ウィンドウの [レンジのウィンドウ サイズ] パラメーターと [アジマスのウィンドウ サイズ] パラメーターには、[入力レーダー データ] パラメーター値の近似正方形ピクセルを作成するために必要な、推定サイズに関する情報が含まれています。

パラメーター

ラベル説明データ タイプ
入力参照レーダー データ

入力参照複素数レーダー データ。

Raster Dataset; Raster Layer
入力セカンダリ レーダー データ

入力セカンダリ複素数レーダー データ。

Raster Dataset; Raster Layer
出力レーダー データ

出力コヒーレンス レーダー データ。

Raster Dataset
偏波バンド
(オプション)

補正する偏波バンド。

デフォルトでは、最初のバンドが選択されています。

String
レンジのウィンドウ サイズ
(オプション)

レンジのウィンドウ サイズ (ピクセル単位)。

デフォルト値は 10 です。

Long
アジマスのウィンドウ サイズ
(オプション)

アジマスのウィンドウ サイズ (ピクセル単位)。

デフォルト値は、近似正方形ウィンドウを作成するために必要なピクセルの最小数です。 たとえば、[レンジのウィンドウ サイズ] パラメーター値が 10 の場合、デフォルトは 3 です。

Long

ComputeCoherence(in_reference_radar_data, in_secondary_radar_data, out_radar_data, {polarization_bands}, {range_window_size}, {azimuth_window_size})
名前説明データ タイプ
in_reference_radar_data

入力参照複素数レーダー データ。

Raster Dataset; Raster Layer
in_secondary_radar_data

入力セカンダリ複素数レーダー データ。

Raster Dataset; Raster Layer
out_radar_data

出力コヒーレンス レーダー データ。

Raster Dataset
polarization_bands
[polarization_bands,...]
(オプション)

補正する偏波バンド。

デフォルトでは、最初のバンドが選択されています。

String
range_window_size
(オプション)

レンジのウィンドウ サイズ (ピクセル単位)。

デフォルト値は 10 です。

Long
azimuth_window_size
(オプション)

アジマスのウィンドウ サイズ (ピクセル単位)。

デフォルト値は、近似正方形ウィンドウを作成するために必要なピクセルの最小数です。 たとえば、range_window_size パラメーター値が 10 の場合、デフォルトは 3 です。

Long

コードのサンプル

ComputeCoherence の例 1 (Python ウィンドウ)

この例では、レンジとアジマスのウィンドウ サイズを 10 ピクセルとして、参照画像とセカンダリ SAR 画像の類似度を計算します。

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

outRadar = arcpy.ia.ComputeCoherence("Reference_SAR.crf",
    "Secondary_SAR_Coreg.crf", "VV", 10, 10) 
outRadar.save("Reference_Secondary_SAR_Coh.crf")
ComputeCoherence の例 2 (スタンドアロン スクリプト)

この例では、範囲と方位角のウィンドウ サイズを 10 ピクセルとして、参照画像とセカンダリ SAR 画像の類似度を計算します。

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

# Set local variables 
in_reference_radar_data = r"C:\SAR\Reference_SAR.crf" 
in_secondary_radar_data = r"C:\SAR\Secondary_SAR_Coreg.crf" 
polarization_bands ="VV" 
range_window_size = 10 
azimuth_window_size=10 

# Execute  
outRadar = arcpy.ia.ComputeCoherence(in_reference_radar_data, 
    in_secondary_radar_data, polarization_bands, range_window_size, 
    azimuth_window_size)  
outRadar.save(r"C:\SAR\Reference_Secondary_SAR_Coh.crf")

ライセンス情報

  • Basic: 次のものが必要 Image Analyst
  • Standard: 次のものが必要 Image Analyst
  • Advanced: 次のものが必要 Image Analyst

関連トピック