Générer un interférogramme (Outils Image Analyst)

Résumé

Crée un interférogramme en différenciant la phase des entrées SLC (Single Look Complex) secondaires et celle des entrées SLC de référence. Cette opération est réalisée par une corrélation croisée des paires de données radar complexes en entrée.

L’outil peut également être utilisé pour supprimer la contribution de phase de la topographie en utilisant un modèle numérique d’élévation (MNE) en entrée.

Le format de l’interférogramme en sortie sera également complexe. Un interférogramme complexe est composé des éléments suivants :

Phase—Représente la différence entre les deux images SLC dans la phase de signal radar. Ces informations de phase sont sensibles à la déformation de la surface, à la topographie et aux effets atmosphériques, ce qui en fait le composant essentiel de l’analyse interférométrique.
Amplitude—Représente l’intensité combinée des signaux radar provenant des deux images SLC. Ces informations d’amplitude représentent l’intensité combinée du signal radar de deux jeux de données SLC. Les zones avec une forte amplitude correspondent souvent à des retours radar plus stables et cohérents, qui peuvent indiquer une cohérence renforcée et une meilleure qualité des données.

Utilisation

Cet outil crée un interférogramme dans un format de données complexe. L’interférogramme est le produit principal créé dans des processus InSAR et DInSAR.
L’entrée doit être une paire de données radar complexes de la même mission et avec les mêmes mode d’acquisition et direction. La paire doit également être issue du même suivi d’orbite. Les données SLC Sentinel-1A et Sentinel-1B peuvent être appariées, car ces types de données correspondent tous les deux à la mission Sentinel-1.
La paire de données radar en entrée doit être coregistrée. Exécutez l’outil Apply Coregistration (Appliquer la coregistration) pour procéder à la coregistration de la paire de données radar en entrée avant d’exécuter l’outil.
Les bandes en copolarisation sont plus communément utilisées que les bandes en polarisation croisée pour optimiser les résultats. Les bandes en polarisation croisée ne sont généralement pas utilisées pour les processus InSAR en raison de leur faible cohérence inhérente.

Spécifiez la bande en copolarisation pour le paramètre Polarization Bands (Bandes de polarisation) afin d’améliorer la cohérence et d’optimiser le traitement en aval.

Les bandes VV et HH sont des exemples de bande en copolarisation.
La sortie de cet outil est affichée à l’aide d’une combinaison de couleurs d’échelle bathymétrique.

Les considérations relatives à la ligne de base pour les applications InSAR sont présentées dans le tableau ci-dessous :


Type de ligne de base	Définition	Ligne de base plus petite	Ligne de base plus grande
Ligne de base perpendiculaire	Composante de la séparation physique entre les deux positions d’acquisition satellite perpendiculaires à la ligne de visée (LOS).	Conserve la cohérence de phase, ce qui est idéal pour la surveillance de la déformation de surface. La sensibilité aux différences d’élévation est moindre, ce qui la rend moins idéale pour la génération de MNE.	Augmente la sensibilité à la topographie, ce qui est idéal pour la génération de MNE.
Ligne de base temporelle	Séparation temporelle entre les deux acquisitions.	Plus grande cohérence dans les zones dynamiques et capture les modifications rapides (les glissements de terrain et les séismes, par exemple).	Permet de détecter les lentes déformations cumulatives telles que les mouvements tectoniques, bien que cela comporte un risque de décorrélation.

Lignes de base perpendiculaires et temporelles

Lorsque le paramètre Remove topographic phase (Supprimer la phase topographique) est activé, l’outil supprime le signal de phase topographique de l’interférogramme en sortie à l’aide d’un MNE. Il s’agit d’une étape de traitement nécessaire pour les processus DInSAR, qui permet d’isoler le signal de déformation.

Pour obtenir des résultats optimaux, le MNE en entrée doit avoir la même résolution que les données radar en entrée.

Etiqueter	Explication	Type de données
Données radar de référence en entrée	Données radar complexes de référence en entrée.	Raster Dataset; Raster Layer
Données radar secondaires en entrée	Données radar complexes secondaires en entrée.	Raster Dataset; Raster Layer
Données radar en sortie	Données radar de l’interférogramme en sortie.	Raster Dataset
Bandes de polarisation (Facultatif)	Bandes de polarisation qui seront interférées. La première bande est sélectionnée par défaut.	String
Supprimer la phase topographique (Facultatif)	Indique si le signal de phase topographique sera supprimé. Activé—L’interférogramme est créé après suppression du signal de phase topographique. La sortie sera un interférogramme avec la phase topographique supprimée. Lorsque cette option est utilisée, la valeur du paramètre DEM Raster (Raster MNE) doit être indiquée. Décoché—Le signal de phase topographique ne sera pas supprimé de l’interférogramme. La sortie sera un interférogramme avec la phase orbitale supprimée. Il s’agit de l’option par défaut.	Boolean
Raster MNE (Facultatif)	Raster MNE qui est utilisé pour estimer le signal de phase résultant du terrain. Ce paramètre est requis lorsque le paramètre Remove topographic phase (Supprimer la phase topographique) est activé.	Mosaic Layer; Raster Layer
Appliquer une correction géodésique (Facultatif)	Indique si le système de référence verticale du MNE en entrée est transformé en hauteur ellipsoïdale. La plupart des jeux de données d’altitude étant référencés par rapport à la hauteur orthométrique au niveau de la mer, une correction est nécessaire dans ce cas pour effectuer la conversion en hauteur ellipsoïdale. Activé—Une correction de géoïde sera effectuée pour convertir la hauteur orthométrique en hauteur ellipsoïdale (en fonction de la géoïde EGM96). Il s’agit de l’option par défaut. Décoché—Aucune correction de géoïde ne sera effectuée. Utilisez cette option uniquement si le MNE est exprimé en hauteur ellipsoïdale.	Boolean
Interférogramme aplati en sortie (Facultatif)	Données radar de l’interférogramme aplati en sortie. Cette sortie compte uniquement pour la correction de la méthode planaire.	Raster Dataset

GenerateInterferogram(in_reference_radar_data, in_secondary_radar_data, out_radar_data, {polarization_bands}, {remove_topographic_phase}, {in_dem_raster}, {geoid}, {out_flattened_interferogram})

Nom	Explication	Type de données
in_reference_radar_data	Données radar complexes de référence en entrée.	Raster Dataset; Raster Layer
in_secondary_radar_data	Données radar complexes secondaires en entrée.	Raster Dataset; Raster Layer
out_radar_data	Données radar de l’interférogramme en sortie.	Raster Dataset
polarization_bands (Facultatif)	Bandes de polarisation qui seront interférées. La première bande est sélectionnée par défaut.	String
remove_topographic_phase (Facultatif)	Indique si le signal de phase topographique sera supprimé. `TOPO`—L’interférogramme est créé après suppression du signal de phase topographique. La sortie sera un interférogramme avec la phase topographique supprimée. Lorsque cette option est utilisée, la valeur du paramètre `in_dem_raster` doit être indiquée. `NONE`—Le signal de phase topographique ne sera pas supprimé de l’interférogramme. La sortie sera un interférogramme avec la phase orbitale supprimée. Il s’agit de l’option par défaut.	Boolean
in_dem_raster (Facultatif)	Raster MNE qui est utilisé pour estimer le signal de phase résultant du terrain. Ce paramètre est requis lorsque la valeur du paramètre `remove_topographic_phase` est `TOPO`.	Mosaic Layer; Raster Layer
geoid (Facultatif)	Indique si le système de référence verticale du MNE en entrée est transformé en hauteur ellipsoïdale. La plupart des jeux de données d’altitude étant référencés par rapport à la hauteur orthométrique au niveau de la mer, une correction est nécessaire dans ce cas pour effectuer la conversion en hauteur ellipsoïdale. `GEOID`—Une correction de géoïde sera effectuée pour convertir la hauteur orthométrique en hauteur ellipsoïdale (en fonction de la géoïde EGM96). Il s’agit de l’option par défaut. `NONE`—Aucune correction de géoïde ne sera effectuée. Utilisez cette option uniquement si le MNE est exprimé en hauteur ellipsoïdale.	Boolean
out_flattened_interferogram (Facultatif)	Données radar de l’interférogramme aplati en sortie. Cette sortie compte uniquement pour la correction de la méthode planaire.	Raster Dataset

Exemple de code

Exemple 1 d’utilisation de la fonction GenerateInterferogram (fenêtre )

Cet exemple génère un interférogramme en sortie.

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

out_radar_data = arcpy.ia.GenerateInterferogram("Reference.crf",
    "Secondary_Coreg.crf", "VV", "TOPO", "dem.tif", "GEOID",
    "flattened_IFG.crf")
outRadar.save(out_radar_data)

Exemple 2 d’utilisation de la fonction GenerateInterferogram (script autonome)

Cet exemple génère un interférogramme en sortie.

# Import system modules
import arcpy
from arcpy.ia import *

# Set local variables
in_ref=r"C:/Data/SAR/Reference.crf"
in_second=r"C:/Data/SAR/Secondary_Coreg.crf"
polarization="VV"
in_DEM=r"C:/Data/SAR/DEM.tif"
out_flat_ifg=r"C:/Data/SAR/flat_IFG.crf"

#Create output
outRadar = arcpy.ia.GenerateInterferogram(in_ref, in_second,
    polarization, "TOPO", in_DEM, "GEOID", out_flat_ifg)
outRadar.save(out_radar)

Environnements

Alignement des cellules, Taille de cellule, Compression, Espace de travail courant, Étendue, Transformations géographiques, ID GPU, Apparier la variable multidimensionnelle, NoData, Système de coordonnées en sortie, Facteur de traitement parallèle, Type de processeur, Pyramide, Statistiques raster, Méthode de rééchantillonnage, Espace de travail temporaire, Raster de capture, Taille de tuile, Union des dimensions

Informations de licence