Calcular coherencia (Image Analyst)

Disponible con licencia de Image Analyst.

Resumen

Calcula la similitud entre los datos de radar complejos de entrada de referencia y secundarios.

La salida es un ráster de coherencia con un rango de valores de 0 a 1 en el que 0 indica ausencia de coherencia y 1 indica coherencia perfecta. Un valor igual o superior a 0,3 se considera un buen valor de coherencia.

Uso

  • Las entradas deben ser un par de datasets de radar complejos corregistrados.

  • El ráster de coherencia de salida puede utilizarse para la detección de cambios coherentes. El ráster de coherencia de salida puede utilizarse en flujos de trabajo InSAR y DInSAR para evaluar la fiabilidad de los píxeles del interferograma.

    En las aplicaciones InSAR y DInSAR, los datos de radar de entrada deben tener la misma geometría de satélite.

  • En el panel Geoprocesamiento, los parámetros Tamaño de la ventana de rango y Tamaño de la ventana acimutal incluyen información sobre el tamaño estimado necesario para crear un píxel cuadrado aproximado para el valor del parámetro Datos de radar de entrada.

  • En la tabla siguiente se muestran las consideraciones de línea base de las aplicaciones InSAR:

    Tipo de línea baseDefiniciónLínea base más pequeñaLínea base más grande

    Línea base perpendicular

    El componente de la separación física entre las dos posiciones de adquisición del satélite que es perpendicular a la línea de visión (LOS).

    Mantiene la coherencia de fase, lo que resulta ideal para monitorizar la deformación de la superficie. Se reduce la sensibilidad a las diferencias de elevación, por lo que es menos idóneo para la generación de DEM.

    Aumenta la sensibilidad a la topografía, por lo que es perfecto para la generación de DEM.

    Línea base temporal

    El tiempo transcurrido entre las dos adquisiciones.

    Mejor coherencia en áreas dinámicas y captura cambios rápidos como deslizamientos de tierra y terremotos.

    Adecuado para detectar deformaciones lentas y acumulativas, como el movimiento tectónico, aunque conlleva riesgo de decorrelación.

    Líneas de base perpendiculares y temporales

Parámetros

EtiquetaExplicaciónTipo de datos
Datos del radar de referencia de entrada

Datos del radar de referencia de entrada complejos.

Raster Dataset; Raster Layer
Datos de radar secundarios de entrada

Datos de radar secundarios de entrada complejos.

Raster Dataset; Raster Layer
Datos de radar de salida

Datos de radar de coherencia de salida.

Raster Dataset
Bandas de polarización
(Opcional)

Se corregirán las bandas de polarización.

La primera banda está seleccionada de forma predeterminada.

String
Tamaño de la ventana de rango
(Opcional)

Tamaño de la ventana de rango en píxeles.

El valor predeterminado es 10.

Long
Tamaño de la ventana acimutal
(Opcional)

Tamaño de la ventana acimutal en píxeles.

El valor predeterminado es el número mínimo de píxeles necesarios para crear una ventana cuadrada aproximada. Por ejemplo, si el valor del parámetro Tamaño de la ventana de rango es 10, el valor predeterminado será 3.

Long

ComputeCoherence(in_reference_radar_data, in_secondary_radar_data, out_radar_data, {polarization_bands}, {range_window_size}, {azimuth_window_size})
NombreExplicaciónTipo de datos
in_reference_radar_data

Datos del radar de referencia de entrada complejos.

Raster Dataset; Raster Layer
in_secondary_radar_data

Datos de radar secundarios de entrada complejos.

Raster Dataset; Raster Layer
out_radar_data

Datos de radar de coherencia de salida.

Raster Dataset
polarization_bands
[polarization_bands,...]
(Opcional)

Se corregirán las bandas de polarización.

La primera banda está seleccionada de forma predeterminada.

String
range_window_size
(Opcional)

Tamaño de la ventana de rango en píxeles.

El valor predeterminado es 10.

Long
azimuth_window_size
(Opcional)

Tamaño de la ventana acimutal en píxeles.

El valor predeterminado es el número mínimo de píxeles necesarios para crear una ventana cuadrada aproximada. Por ejemplo, si el valor del parámetro range_window_size es 10, el valor predeterminado será 3.

Long

Muestra de código

Ejemplo 1 de ComputeCoherence (ventana de Python)

Este ejemplo calcula la similitud entre las imágenes SAR de referencia y secundarias utilizando un tamaño de ventana de rango y acimutal de 10 píxeles.

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")
Ejemplo 2 de ComputeCoherence example 2 (secuencia de comandos independiente)

Este ejemplo calcula la similitud entre las imágenes SAR de referencia y secundarias utilizando un tamaño de ventana de rango y acimutal de 10 píxeles.

# 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")

Información de licenciamiento

  • Basic: Requiere Image Analyst
  • Standard: Requiere Image Analyst
  • Advanced: Requiere Image Analyst

Temas relacionados