430 likes | 616 Views
INTERFERÓMETRIA EM RADARES DE ABERTURA SINTÉTICA. José M. B. Dias Instituto Superior Técnico Instituto de Telecomunicações Portugal. Sumário. O Conceito de Abertura Sintética Formação de Imagem SAR Princípio de Interferometria SAR (InSAR) Geração de Interferogramas
E N D
INTERFERÓMETRIAEM RADARES DE ABERTURA SINTÉTICA José M. B. Dias Instituto Superior Técnico Instituto de Telecomunicações Portugal IST/IT
Sumário • O Conceito de Abertura Sintética • Formação de Imagem SAR • Princípio de Interferometria SAR (InSAR) • Geração de Interferogramas • Estimação de Fase Absoluta • Determinação de Mapas Digitais do Terreno IST/IT
- Geometria de um Sistema de Abertura Sintética Resolução ERS-1 DRA-Airborne SAR !!!!!! IST/IT
- Compressão de Impulso (Resolução Lateral) Ideia: • Enviar sinais de longa duração • Comprimi-los na recepção IST/IT
- Compressão de Impulso (cont.) Modulação Linear de Frequência (chirp) [Dick, 53],[Darlington, 54] IST/IT
- Resolução Azimutal • Sinal Transmitido • Sinal Recebido Atraso Ganho da antena Fase de propagação IST/IT
- Resolução Azimutal (cont.) Aproximação de Fresnel Sinal Recebido Chirp • Duração: • Shirp rate: IST/IT
- Resolução Azimutal (cont.) s(t) é um chirp Mas • Resolução: • Independente de r0 • Proporcional a D IST/IT
- Processador de SAR (aproximaçaõ de Fresnel) Processador básico de SAR em modo Strip-Mapping Desmodulação de quadratura Chirp rate: Duração: Compressão em distância Compressão em azimute IST/IT
- Presença de “range walk” e “range curvature” 3 2 1 1- Spaceborn sem movimento Dificuldade Encontrar a assinatura de um alvo no plano (r,y) 2- Spaceborn com movimento 3-Airborn sem movimento IST/IT
- Metodologia 2D Recontrução de Frente de Onda IST/IT
- Reconstrução de frente de onda (cont) (Aproximação de fase estacionária) Transformada de Fourier de f(x,y) IST/IT
- Reconstrução de frente de onda (cont) • Estimador de F • Filtro Inverso • Filtro Adaptado IST/IT
- Reconstrução de frente de onda (cont) Resposta Impulsional do filtro Adaptado IST/IT
-Técnicas de Formação de Imagem • Tópicos de Investigação em SAR • Técnicas de Formação de Imagem • Diversidade na polarização • Pertubações na trajectória do sensor • Detecção/Identificação de alvos em movimento • Estimação de reflectividade • Estimação de mapas de altura-InSAR IST/IT
Interferómetria SAR IST/IT
Absolute Phase estimation in InSAR (Interferometric SAR) InSAR Problem: Estimate f2- f1 from signal read by s1 and s2 IST/IT
Dificuldades • É apenas possível determinar módulo 2p ( é • uma função periódica de período 2p. ) • Os dados radar estão imersos em ruído (por • exemplo, ruído térmico) • O termo não é nulo devido a descorrelação • dos tipos geométrico, temporal e espacial. IST/IT
Decorrelation Factors • Spatial decorrelation (misregistration) • Temporal decorrelation (scatterer displacements) • Geometric Decorrelation (finite resolution) • Electronic Noise IST/IT
Geometria InSAR IST/IT
Registo IST/IT
Exemplo Imagem de amplitude Disparidades da imagem reamostrada IST/IT
Coerência Interferograma IST/IT
Histogramas da coerência Após registo fino Após registo grosso IST/IT
Classical Approaches • Phase Unwrapping • Bayesian and Regularization • Take a two-step approach • Step 1: Principal phase values are inferred from noisy InSAR image pairs • Step 2: Phase is unwrapped by determining the 2p multiples Phase Unwrapping IST/IT
Phase Unwrapping Methods • Path following • Phase is unwrapped through selected image paths [Goldstein et al., 88], [Lim et al., 95], [Xu et al., 95] • Minimum norm • Cast the unwrapping problem into the minimization of a Lp norm [Flyn et al., 97], [Costantini, 98], [Ghiglia et al., 94], [Ghiglia et al., 98]. IST/IT
Path Following Methods - If then + and IST/IT
Minimum Norm Methods Define Discrete optimization problem Exact solutions only for p = 1 [Flyn et al., 97], [Costantini, 98] Network programming algorithms IST/IT
Bayesian and Regularization Enforce smoothness on the solution by including prior and/or regularization terms • Nonlinear Filtering [Leitão et al., 97], [Dias et al., 98] • Quadratic regularization [Marroquin et al., 95] IST/IT
Interferogram Original Some results using classical algorithms Smoothed Interferogram IST/IT
Some results (cont.) Goldstein (smooth) Flyn (smooth) Goldstein (nonsmooth) Flyn (nonsmooth) IST/IT
Proposed Bayesian Approach Prior: Compound Gauss Markov Random Field [Jeng et al., 91] IST/IT
ZpM - Algorithm IST/IT
Results IST/IT
Results Posterior density Phase error histogram IST/IT
Results IST/IT
Results Ramp with 0.1rad/pixel IST/IT
Results IST/IT
Results Estimated surface (15:15:120)+2 } truth surface (15:15:120) IST/IT
Results Estimated Original IST/IT
Outro Exemplo IST/IT