Aplicar aplanamiento radiométrico del terreno (Image Analyst)

Disponible con licencia de Image Analyst.

Resumen

Corrige distorsiones radiométricas debidas a la topografía en los datos de radar de apertura sintética (SAR) de entrada.

Debido a la visión lateral de los sensores SAR, las entidades orientadas hacia el sensor aparecen artificialmente más brillantes y las entidades orientadas hacia el lado opuesto del sensor aparecen artificialmente más oscuras. El aplanamiento radiométrico del terreno normaliza los valores de retrodispersión de modo que las variaciones de los valores se deban a las propiedades de dispersión de la superficie.

El aplanamiento radiométrico del terreno es necesario para obtener una retrodispersión significativa que se pueda relacionar directamente con las propiedades de dispersión de la superficie de las entidades de una imagen de SAR sobre cualquier terreno.

Uso

  • Los datos de SAR de entrada deben estar calibrados como beta cero.

  • Utilice la herramienta Aplicar calibración radiométrica para calibrar los datos de SAR como beta cero.

  • Si el DEM de entrada no abarca todo el dataset de SAR, la herramienta generará valores NoData para los píxeles fuera de la extensión de DEM correspondientes a salidas de gamma cero, sigma cero, área de dispersión y distorsión geométrica. Para la salida de máscara de distorsión geométrica, la herramienta generará valores de undetermined para los píxeles fuera de la extensión de DEM.

  • El DEM de entrada debe estar en el sistema de coordenadas geográficas WGS 1984 (EPSG:4326).

Parámetros

EtiquetaExplicaciónTipo de datos
Datos de radar de entrada

Los datos de radar de entrada.

Los datos se deben calibrar radiométricamente como beta cero.

Raster Dataset; Raster Layer
Datos de radar de salida

Los datos de radar con aplanamiento radiométrico del terreno.

Raster Dataset
Ráster DEM

El DEM de entrada.

El DEM se utilizará para estimar el área iluminada local y el ángulo de incidencia local.

Mosaic Layer; Raster Layer
Aplicar corrección geoide
(Opcional)

Especifica si el sistema de referencia vertical del DEM de entrada se transformará a altura elipsoidal. La mayoría de los datasets de elevación hacen referencia a altura ortométrica a nivel del mar, por lo que es necesario realizar una corrección en estos casos para convertirla a altura elipsoidal.

  • Activada: se realizará una corrección geoide para convertir la altura ortométrica en altura elipsoidal (basándose en el geoide EGM96). Esta es la opción predeterminada.
  • Desactivada: no se realizará ninguna corrección geoide. Utilice esta opción solamente si el DEM ya está expresado en altura elipsoidal.

Boolean
Bandas de polarización
(Opcional)

Las bandas de polarización cuyo terreno se aplanará radiométricamente.

La primera banda está seleccionada de forma predeterminada.

String
Tipo de calibración
(Opcional)

Especifica si la salida tendrá un terreno aplanado con sigma cero o gamma cero.

  • Gamma ceroLa retrodispersión beta cero se corregirá utilizando un cómputo exacto de un área usando un DEM. Esta es la opción predeterminada.
  • Sigma ceroLa retrodispersión beta cero se corregirá utilizando el área de unidad de un plano que sea tangente localmente al DEM.
String
Área de dispersión de salida
(Opcional)

El dataset de radar de área de dispersión.

Raster Dataset
Distorsión geométrica de salida
(Opcional)

El dataset de radar de distorsión geométrica de 4 bandas. La primera banda es la pendiente del terreno, la segunda banda es el ángulo de visión, la tercera banda es la proporción de escorzo y la cuarta banda es el ángulo de incidencia local.

Raster Dataset
Máscara de distorsión geométrica de salida
(Opcional)

El dataset de radar de máscara de distorsión geométrica de 1 banda. Los píxeles se clasifican mediante seis valores únicos, uno para cada tipo de distorsión:

  • Indeterminada Rojo, Verde, Azul = 160, 177, 186: valor de 0
  • Escorzo Rojo, Verde, Azul = 0, 154, 222: valor de 1
  • Alargamiento Rojo, Verde, Azul = 175, 88, 186: valor de 2
  • Sombra Rojo, Verde, Azul = 255, 198, 30: valor de 3
  • Inversión por relieve Rojo, Verde, Azul = 0, 205, 108: valor de 4
  • Inversión por relieve y sombra Rojo, Verde, Azul = 166, 118, 29: valor de 5

Raster Dataset

ApplyRadiometricTerrainFlattening(in_radar_data, out_radar_data, in_dem_raster, {geoid}, {polarization_bands}, {calibration_type}, {out_scattering_area}, {out_geometric_distortion}, {out_geometric_distortion_mask})
NombreExplicaciónTipo de datos
in_radar_data

Los datos de radar de entrada.

Los datos se deben calibrar radiométricamente como beta cero.

Raster Dataset; Raster Layer
out_radar_data

Los datos de radar con aplanamiento radiométrico del terreno.

Raster Dataset
in_dem_raster

El DEM de entrada.

El DEM se utilizará para estimar el área iluminada local y el ángulo de incidencia local.

Mosaic Layer; Raster Layer
geoid
(Opcional)

Especifica si el sistema de referencia vertical del DEM de entrada se transformará a altura elipsoidal. La mayoría de los datasets de elevación hacen referencia a altura ortométrica a nivel del mar, por lo que es necesario realizar una corrección en estos casos para convertirla a altura elipsoidal.

  • GEOIDSe realizará una corrección geoide para convertir la altura ortométrica en altura elipsoidal (basándose en el geoide EGM96). Esta es la opción predeterminada.
  • NONENo se realizará ninguna corrección geoide. Utilice esta opción solamente si el DEM ya está expresado en altura elipsoidal.
Boolean
polarization_bands
[polarization_bands,...]
(Opcional)

Las bandas de polarización cuyo terreno se aplanará radiométricamente.

La primera banda está seleccionada de forma predeterminada.

String
calibration_type
(Opcional)

Especifica si la salida tendrá un terreno aplanado con sigma cero o gamma cero.

  • GAMMA_NOUGHTLa retrodispersión beta cero se corregirá utilizando un cómputo exacto de un área usando un DEM. Esta es la opción predeterminada.
  • SIGMA_NOUGHTLa retrodispersión beta cero se corregirá utilizando el área de unidad de un plano que sea tangente localmente al DEM.
String
out_scattering_area
(Opcional)

El dataset de radar de área de dispersión.

Raster Dataset
out_geometric_distortion
(Opcional)

El dataset de radar de distorsión geométrica de 4 bandas. La primera banda es la pendiente del terreno, la segunda banda es el ángulo de visión, la tercera banda es la proporción de escorzo y la cuarta banda es el ángulo de incidencia local.

Raster Dataset
out_geometric_distortion_mask
(Opcional)

El dataset de radar de máscara de distorsión geométrica de 1 banda. Los píxeles se clasifican mediante seis valores únicos, uno para cada tipo de distorsión:

  • Indeterminada Rojo, Verde, Azul = 160, 177, 186: valor de 0
  • Escorzo Rojo, Verde, Azul = 0, 154, 222: valor de 1
  • Alargamiento Rojo, Verde, Azul = 175, 88, 186: valor de 2
  • Sombra Rojo, Verde, Azul = 255, 198, 30: valor de 3
  • Inversión por relieve Rojo, Verde, Azul = 0, 205, 108: valor de 4
  • Inversión por relieve y sombra Rojo, Verde, Azul = 166, 118, 29: valor de 5

Raster Dataset

Muestra de código

Ejemplo 1 de ApplyRadiometricTerrainFlattening (ventana de Python)

En este ejemplo, se corrige un dataset de radar polarizado cruzado utilizando un DEM y parámetros gamma cero.

import arcpy
arcpy.env.workspace = "D:\Data\SAR\S1\20181014"
outRadar = arcpy.ia.ApplyRadiometricTerrainFlattening(
     "IW_manifest_TNR_CalB0_Dspk.crf", "D:\Data\DEM\dem.tif", "GEOID", 
     "VH;VV", "GAMMA_NOUGHT")
outRadar.save("IW_manifest_TNR_CalB0_Dspk_RTFG0.crf")
Ejemplo 2 de ApplyRadiometricTerrainFlattening (script independiente)

En este ejemplo, se corrige un dataset de radar polarizado cruzado utilizando un DEM y parámetros gamma cero.

# 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 = "D:\Data\SAR\S1\20181014\IW_manifest_TNR_CalB0_Dspk.crf"
out_radar = "D:\Data\SAR\S1\20181014\IW_manifest_TNR_CalB0_Dspk_RTFG0.crf"
in_dem_raster = "D:\Data\DEM\dem.tif"
ApplyGeoid = "GEOID" 
polarization = "VH;VV"
calibration_type = "GAMMA_NOUGHT"

# Execute 
outRadar = arcpy.ia.ApplyRadiometricTerrainFlattening(
     in_radar, in_dem_raster, ApplyGeoid, polarization, calibration_type)
outRadar.save(out_radar)

Información de licenciamiento

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

Temas relacionados