Técnicas de interferometría diferencial SAR para la obtención de mapas de deformación del terreno y control de edificios

Óscar Mora

Description

Resumen
En este artículo se presentan diversas técnicas para la obtención de mapas de deformación del terreno usando datos SAR (Radar de Apertura Sintética). El primero de los algoritmos estima las componentes lineal y no lineal del desplazamiento, el error del Mapa de Elevaciónd el Terreno (MET) utilizado para eliminar el término topográfico de los interferogramas, y los artefactos atmosféricos. La entrada de datos se compone de un conjunto de Interferogramas Diferenciales de baja resolución espacial y sus correspondientes mapas de coherencia. El algoritmo no tiene restricciones teóricas respecto al mínimo número de interferogramas requerido y las dimensiones de las líneas de base. Sin embargo, los mejores resultados se obtienen con grandes conjuntos de datos de línea de base corta. El algoritmo selecciona aquellos píxeles que presentan una alta coherencia, evitando aquellos afectados por decorrelación temporal. La segunda técnica presentada se encarga de estimar la deformación relativa respecto al terreno de estructuras, tales como edificios. Dicha técnica se basa en el cálculo del residuo obtenido al restar los INterferogramas Diferenciales de baja y alta resolución. Utilizando esta información es posible obtener la evolución temporal de la deformación de puntos coherentes. Los algoritmos presentados han sido validados como datos SAR de los satélites ERS.
Abstract
In this paper, advanced techiques for the generation of Terrain deformation maps using SAR (Synthetic Aperture Radar) data are presented. One of the algorithms estimates the linear and non-linear components of the displacement, the error of the Digitial Elevation Model (DEM) used to cancel the topographic term from the interferograms, and the atmospheric artifacts. The input data is a set of spatial low resolution Differential Interferograms and their associated coherence images. The algorithm has no theoretical restrictions with respect to the minimum number of required interferograms and their baseline dimensions. However, better results are obtained with large datasets of short baseline interferograms. The algorithm takes advantage of those pixels presenting a good coherence level in the whole set of interferograms, avoiding the rest affected by temporal decorrelation. Another algorithm has been developed to estimate the deformation respect to the ground of structures, such as buildings. The main point is the calculation of the phase residue obtained by subtracting equivalent low and high resolution versions of Differential Interferograms. Using this information it is possible to obtain the deformation evolution of high coherent points, like buildings. THese algorithms have been tested with ERS SAR data.

Sujets

• SAR
• Radar
• Interferometría
• Subsidencia
• Deformación

Informations

• Publié par: erevistas
• Publié le: 01 janvier 2004
• Nombre de lectures: 23
• Langue: Español

Extrait

Revista de Teledetección. 2004. 21: 51-54
Técnicas de interferometría diferencial SAR para la
obtención de mapas de deformación del terreno y
control de edificios
Óscar Mora*, **, Jordi Mallorquí*, Ricardo Lanari*** y Antoni Broquetas*
[e-mail]
* Universitat Politècnica de Catalunya (UPC). Dept. Teoría del Senyal i Comunicacions (TSC).
Campus Nord UPC D3-211.
Jordi Girona 1-3, 08034 Barcelona, Spain
** Institut Cartogràfic de Catalunya (ICC). Unitat de Teledetecció.
Parc de Montjuïc. 08038 Bar
*** Instituto per il Rilevamento Elettromagnetico dell´Ambiente (IREA). Consiglio Nazionale delle
Ricerche (CNR). Via Diocleziano, 328 - 80124 Napoli, Italia
RESUMEN ABSTRACT
En este artículo se presentan diversas técnicas In this paper, advanced techiques for the genera-
para la obtención de mapas de deformación del tion of Terrain deformation maps using SAR (Synt-
terreno usando datos SAR (Radar de Apertura Sinté- hetic Aperture Radar) data are presented. One of the
tica). El primero de los algoritmos estima las com- algorithms estimates the linear and non-linear com-
ponentes lineal y no lineal del desplazamiento, el ponents of the displacement, the error of the Digi-
error del Mapa de Elevaciónd el Terreno (MET) uti- tial Elevation Model (DEM) used to cancel the
lizado para eliminar el término topográfico de los topographic term from the interferograms, and the
interferogramas, y los artefactos atmosféricos. La atmospheric artifacts. The input data is a set of spa-
entrada de datos se compone de un conjunto de tial low resolution Differential Interferograms and
Interferogramas Diferenciales de baja resolución their associated coherence images. The algorithm
espacial y sus correspondientes mapas de coheren- has no theoretical restrictions with respect to the
cia. El algoritmo no tiene restricciones teóricas res- minimum number of required interferograms and
pecto al mínimo número de interferogramas requeri- their baseline dimensions. However, better results
do y las dimensiones de las líneas de base. Sin are obtained with large datasets of short baseline
embargo, los mejores resultados se obtienen con interferograms. The algorithm takes advantage of
grandes conjuntos de datos de línea de base corta. El those pixels presenting a good coherence level in
algoritmo selecciona aquellos píxeles que presentan the whole set of interferograms, avoiding the rest
una alta coherencia, evitando aquellos afectados por affected by temporal decorrelation. Another algo-
decorrelación temporal. La segunda técnica presen- rithm has been developed to estimate the deforma-
tada se encarga de estimar la deformación relativa tion respect to the ground of structures, such as
respecto al terreno de estructuras, tales como edifi- buildings. The main point is the calculation of the
cios. Dicha técnica se basa en el cálculo del residuo phase residue obtained by subtracting equivalent
obtenido al restar los INterferogramas Diferenciales low and high resolution versions of Differential
de baja y alta resolución. Utilizando esta informa- Interferograms. Using this information it is possible
ción es posible obtener la evolución temporal de la to obtain the deformation evolution of high cohe-
deformación de puntos coherentes. Los algoritmos rent points, like buildings. THese algorithms have
presentados han sido validados como datos SAR de been tested with ERS SAR data.
los satélites ERS.
PALABRAS CLAVE: SAR, radar, interferometría,
subsidencia, deformación.
N.º 21- Junio 2004 51O. Mora, J. Mallorquí, R.o Lanari y A. Broquetas
Una vez estimada la componente lineal de laINTRODUCCIÓN
deformación, la técnica presentada propone una
solución para el cálculo del término no lineal (MoraAl generar un interferograma combinando dos
2002a). De este modo se puede analizar la evolu-imágenes SAR, los incrementos de fase entre pixe-
ción temporal del movimiento. Las dos etapas prin-les pueden expresarse como:
cipales del algoritmo propuesto son las siguientes:
DY = DY + DY + DY +int flat topo mov • Cálculo del residuo de fase resultante de restar
+ DY + DY (1)atmos noise a las fases originales el modelo lineal obtenido
anteriormente. Este residuo estará compuesto
donde DY es el término de tierra plana, DY seflat topo por el movimiento no lineal, los artefactos
corresponde con la fase topográfica, DY es lamov atmosféricos y otros términos de ruido.
componente debida al desplazamiento del terreno en • Separación de los términos de movimiento no-
la dirección slant-range, DY es la fase relaciona- lineal y atmósfera, mediante un filtrado espacio-atmos
da con los artefactos atmosféricos, y DY incluye temporal.noise
diversos factores de decorrelación y ruido térmico.
Finalmente, la combinación de las componentesLa interferometría diferencial clásica trata de eli-
lineal y no lineal del movimiento resultan en unaminar las componentes conocidas de la ecuación (1)
completa descripción de la deformación del terrenopara obtener la fase relacionada con el desplaza-
durante el periodo de tiempo de observación.miento del terreno. Esto se puede conseguir gene-
rando la fase topográfica con la ayuda de un MET
Movimientos de edificiosexterno. Sin embargo, estos interferogramas dife-
renciales suelen estar seriamente afectados por arte- La técnica propuesta anteriormente está pensada
factos atmosféricos, dificultando la obtención de para el control de deformaciones que ocurren a gran
resultados precisos. escala, por ejemplo mediante el uso de pixeles de
100 m x 100 m. Sin embargo, no es apropiada para
el estudio de deformaciones muy localizadas que
Movimientos de áreas extensas afecten a edificios aislados. Por lo tanto, se propo-
ne un nuevo algoritmo que extiende las capacidades
La primera parte de la técnica propuesta conside- del método anterior utilizando un conjunto de inter-
ra un modelo lineal de deformación del terreno ferogramas diferenciales de línea de base corta para
(Ferretti 2000). Esto significa que se asume que la minimizar los efectos de la decorrelación espacial
parte más importante del movimiento es un vector (Mora 2002b). El proceso es el siguiente:
de velocidad constante. El modelo lineal puede ser
• Se generan dos conjuntos de interferogramasajustado a los datos utilizando un conjunto de inter-
diferenciales del área de interés; uno de ellosferogramas diferenciales con diferentes intervalos
con alta resolución espacial (un sólo look) y eltemporales. El proceso completo es el siguiente
otro a baja resolución (multi-look).(Mora 2002a):
• Obtención del residuo resultante de restar las
fases de baja resolución a las de alta. Este resi-• Selección de los pixeles basada en su coheren-
duo estará compuesto por el movimiento relati-cia. Este paso es necesario para eliminar aque-
vo de las estructuras del tamaño del orden delllos puntos seriamente afectados por ruido.
pixel utilizado en la alta resolución con respec-• Triangulación entre los pixeles seleccionados.
to al terreno sobre el que están instaladas.Esto permite trabajar con los incrementos de
• Cálculo de la componente lineal de movimientofase entre puntos cercanos, minimizando los
de forma similar al algoritmo de áreas extensas.efectos atmosféricos (Hanssen 1998).
• Estimación de la componente no-lineal median-• Ajuste de los datos a un modelo lineal de defor-
te el estudio temporal de los residuos obtenidosmación. También se tiene en cuenta la estimación
de restar el modelo lineal a los datos originales.del error en altura del MET utilizado para la eli-
minación de la componente topográfica de la fase. Debe destacarse que la información obtenida con
•Integración de los resultados obtenidos para las dos técnicas (áreas extensas y edificios) es abso-
cada par de puntos con el fin de estimar la velo- lutamente complementaria, y se deben combinar para
cidad de deformación absoluta. realizar un estudio completo de la zona de interés.
52 N.º 21 - Junio 2004Técnicas de interferometría diferencial SAR para la obtención de mapas de deformación del terreno...
RESULTADOS EXPERIMENTALES
Se han seleccionado dos zonas urbanas para vali-
dar los métodos presentados. La primera contiene
dos poblaciones rodeadas de bosque y la segunda
corresponde a una ciudad con escasa vegetación.
En la primera zona se ha aplicado el algoritmo de
áreas extensas. La Figura 1 muestra la deformación
lineal y el error del MET obtenidos para una zona
de 10 Km x 16 Km usando un conjunto de 28 inter-
Figura 2. Detalles de dos áreas de interés. Puntos sin
ferogramas de línea de base corta. Los resultados deformación (A y D) y puntos afectados por movimiento (B
muestran que prácticamente todos los pixeles selec- y C).
cionados no presentan deformación. Sin embargo,
se ha detectado una pequeña zona en la parte supe-
rior de la imagen con problemas de subsidencia