Calculer les indices SAR (Image Analyst)

Disponible avec une licence Image Analyst.

Synthèse

Calcule divers indices SAR pour les données SAR (radar à synthèse d’ouverture), tels que l’indice de végétation radar (RVI), l’indice de dégradation des forêts radar (RFDI) et l’indice de structure de la canopée (CSI).

Les formules utilisées pour ces indices dépendent des polarisations disponibles dans le jeu de données radar en entrée.

Utilisation

  • Les données radar en entrée doivent être en unités linéaires. Si les données sont en décibels (dB), utilisez l’outil Convertir les unités SAR pour convertir les données radar en unités linéaires.

  • Pour les données SAR en entrée à polarisation quadruple, l’outil fournit les trois options de formule possibles pour l’indice RVI. RVI requiert des polarisations HH, HV ou VV, VH ou HH, HV, VH, VV.

    Équations RVI

    L’outil fournit deux options possibles pour l’indice RFDI. RFDI requiert des polarisations HH, HV ou VV, VH.

    Équations RFDI

    L’outil fournit une option pour l’indice CSI. CSI requiert des polarisations HH, HV.

    Équation CSI

  • Pour les indices RFDI et RVI, les villes, déserts et plans d’eau peuvent avoir des valeurs correspondant à des paysages arides, dégradés et déboisés, car une végétation très limitée, voire inexistante, est inhérente à ces régions. Même si la plupart des valeurs sont comprises entre 0 et 1, des points aberrants peuvent résulter de la double réflexion ou de la diffusion volumétrique.

Paramètres

ÉtiquetteExplicationType de données
Données radar en entrée

Données radar en entrée.

Raster Dataset; Raster Layer
Raster en sortie

Raster d’indice SAR en sortie.

Raster Dataset
Index
(Facultatif)

Spécifie l’indice SAR qui est calculé.

  • Indice de végétation radar (RVI)Indice de végétation radar (RVI) qui est utilisé. Le RVI est le ratio de rétrodiffusion en polarisation croisée par rapport à la rétrodiffusion totale de toutes les polarisations. Les valeurs sont comprises entre 0 et 1. Des valeurs RVI proches de 0 désignent des paysages arides, tandis que des valeurs supérieures indiquent des paysages végétalisés. Il s’agit de l’option par défaut.
  • Indice de dégradation des forêts radar (RFDI)L’indice de dégradation des forêts radar (RFDI) est utilisé. Le RFDI est la différence normalisée entre la rétrodiffusion en copolarisation et la rétrodiffusion en polarisation croisée. Des valeurs RFDI inférieures (inférieures à 0,3) indiquent une forêt plus dense. Des valeurs RFDI modérées (comprises entre 0,4 et 0,6) correspondent à des forêts dégradées. Des valeurs RFDI plus élevées (supérieures à 0,6) indiquent des paysages déboisés.
  • Indice de structure de la canopée (CSI)L’indice de structure de la canopée (CSI) est utilisé. Le CSI est la différence normalisée de la rétrodiffusion en copolarisation (HH, VV). Les valeurs sont comprises entre -1 et +1. Les canopées dominées par des structures verticales ont des valeurs CSI proches de -1, tandis que celles dominées par des structures horizontales ont des valeurs CSI proches de 1. Cette option est prise en charge uniquement lorsque les données radar en entrée contiennent des bandes HH et VV.
String
Bandes de polarisation
(Facultatif)

Spécifie les bandes de polarisation qui sont utilisées dans le calcul des indices.

Ce paramètre est pris en charge uniquement lorsque la valeur du paramètre Données radar en entrée est un jeu de données SAR à polarisation quadruple et la valeur du paramètre Indice est Indice de végétation radar (RVI) ou Indice de dégradation des forêts radar (RFDI).

  • HH, HVLes bandes de type horizontal-horizontal et horizontal-vertical sont utilisées dans le calcul des indices (double polarisation). Il s’agit de l’option par défaut.
  • VV, VHLes bandes de type vertical-vertical et vertical-horizontal sont utilisées dans le calcul des indices (double polarisation).
  • HH, HV, VH, VVLes bandes de type horizontal-horizontal, horizontal-vertical, vertical-horizontal et vertical-vertical sont utilisées dans le calcul des indices (polarisation quadruple).
String

ComputeSARIndices(in_radar_data, out_raster, {index}, {polarization_bands})
NomExplicationType de données
in_radar_data

Données radar en entrée.

Raster Dataset; Raster Layer
out_raster

Raster d’indice SAR en sortie.

Raster Dataset
index
(Facultatif)

Spécifie l’indice SAR qui est calculé.

  • RVIIndice de végétation radar (RVI) qui est utilisé. Le RVI est le ratio de rétrodiffusion en polarisation croisée par rapport à la rétrodiffusion totale de toutes les polarisations. Les valeurs sont comprises entre 0 et 1. Des valeurs RVI proches de 0 désignent des paysages arides, tandis que des valeurs supérieures indiquent des paysages végétalisés. Il s’agit de l’option par défaut.
  • RFDIL’indice de dégradation des forêts radar (RFDI) est utilisé. Le RFDI est la différence normalisée entre la rétrodiffusion en copolarisation et la rétrodiffusion en polarisation croisée. Des valeurs RFDI inférieures (inférieures à 0,3) indiquent une forêt plus dense. Des valeurs RFDI modérées (comprises entre 0,4 et 0,6) correspondent à des forêts dégradées. Des valeurs RFDI plus élevées (supérieures à 0,6) indiquent des paysages déboisés.
  • CSIL’indice de structure de la canopée (CSI) est utilisé. Le CSI est la différence normalisée de la rétrodiffusion en copolarisation (HH, VV). Les valeurs sont comprises entre -1 et +1. Les canopées dominées par des structures verticales ont des valeurs CSI proches de -1, tandis que celles dominées par des structures horizontales ont des valeurs CSI proches de 1. Cette option est prise en charge uniquement lorsque les données radar en entrée contiennent des bandes HH et VV.
String
polarization_bands
(Facultatif)

Spécifie les bandes de polarisation qui sont utilisées dans le calcul des indices.

Ce paramètre est pris en charge uniquement lorsque la valeur du paramètre in_radar_data est un jeu de données SAR à polarisation quadruple et la valeur du paramètre index est RVI ou RFDI.

  • HH_HVLes bandes de type horizontal-horizontal et horizontal-vertical sont utilisées dans le calcul des indices (double polarisation). Il s’agit de l’option par défaut.
  • VV_VHLes bandes de type vertical-vertical et vertical-horizontal sont utilisées dans le calcul des indices (double polarisation).
  • HH_HV_VH_VVLes bandes de type horizontal-horizontal, horizontal-vertical, vertical-horizontal et vertical-vertical sont utilisées dans le calcul des indices (polarisation quadruple).
String

Exemple de code

Exemple 1 d’utilisation de l’outil ComputeSARIndices (fenêtre Python)

Cet exemple calcule le RVI à partir d’un jeu de données SAR à polarisation quadruple.

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")
Exemple 2 d’utilisation de l’outil ComputeSARIndices (script autonome)

Cet exemple calcule le RVI à partir d’un jeu de données SAR à polarisation quadruple.

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

Rubriques connexes