ESCOLA NAVAL
RADAR DE ABERTURA SINTÉTICA
APLICAÇÕES EM ISAR
José M. B. Dias
Instituto Superior Técnico
Instituto de Telecomunicações
Portugal
IST/IT
1
ESCOLA NAVAL
Sumário
•
•
•
•
•
O Radar Lateral
O Conceito de Abertura Sintética
Formação de Imagem SAR
Exemplos
Aplicações ISAR
IST/IT
2
ESCOLA NAVAL
Radar
• Sistema activo - opera 24 horas por dia
• Depende pouco das condições meteorológicas
• Gera imagens de microondas
IST/IT
3
ESCOLA NAVAL
Resolução: 1m X 1m
Dimensão: 630m x 630m
Pixeis:
IST/IT
940m x 940m
4
ESCOLA NAVAL
IST/IT
Banda
Espectral
Frequência (GHz)
Comprimento de Onda
(cm)
P
0.3 – 1.0
100 –30
L
1.0 – 2.0
30 - 15
S
2.0 – 4.0
15 – 7.5
C
4.0 – 8.0
7.5 – 3.75
X
8.0 – 12.5
3.75 – 2.4
Ku
12.5 – 18.0
2.4 – 1.67
K
18.0 – 26.5
1.67 – 1.1
Ka
26.5 – 40.0
1.1 – 0.75
Comprimento de onda e frequência das microondas.
5
ESCOLA NAVAL
Esquematização da aquisição de informação através do sistema SLAR.
IST/IT
6
ESCOLA NAVAL
Modos de operação SAR
IST/IT
7
ESCOLA NAVAL - Geometria de um Sistema de Abertura Sintética
Resolução
x 
c/2 1
B sin 
ERS-1
y 

D
r
DRA-Airborne SAR
  5.6cm
  3.1cm
B  15MHz
B  100MHz
D  10m
D  1.8m
  23º
  23º
r  850Km
r  14Km
x  25m
x  3.8m
y  0.24km
y  4.8km
!!!!!!
IST/IT
8
ESCOLA NAVAL - Compressão de Impulso (Resolução Lateral)
Ideia:
• Enviar sinais de longa duração
• Comprimi-los na recepção
IST/IT
9
ESCOLA NAVAL
- Compressão de Impulso (cont.)
Modulação Linear de Frequência (chirp)
[Dick, 53],[Darlington, 54]
p(t )  rect (t / T) exp(j / 2t 2 )
h(t )  p* (t )
| H ( ) |
(2 ) / 
T
2
( )  
c
2
IST/IT

2
p(t )  p* (t )
2
T
10
ESCOLA NAVAL
p(t )  p* (t )
IST/IT
- Compressão de Impulso (exemplo)
p(t )  p* (t )
11
ESCOLA NAVAL
- Resolução Azimutal
•Sinal Transmitido
p(t )
•Sinal Recebido
p(t  2r / c)w( y  y0 ) exp j 2kr 
Atraso
Ganho Fase de
da antena propagação
k
IST/IT

c

12
2

ESCOLA NAVAL - Resolução Azimutal (cont.)
Aproximação de Fresnel
( y  y0 ) 2
r  r0  ( y  y0 )  r0 
2r0
2
2
y  y0  v(t  t0 )
Sinal Recebido
s(t )  p(t  2r / c)w(v(t  t0 )) exp(t )
Chirp
• Duração: T 
(t ) 
IST/IT
4 r0

2 v 2 (t  t0 ) 2

r0
L
v
4 v 2
• Chirp rate:  
r0
13
ESCOLA NAVAL - Resolução Azimutal (cont.)
s(t) é um chirp
r
t  2  0
T 2vL
y  vt 
r0
2L
Mas
L

D
y 
r0
D
2
Resolução:
• Independente de r0
• Proporcional a D
IST/IT
14
ESCOLA NAVAL - Processador de SAR (aproximaçaõ de Fresnel)
Processador básico de SAR em modo Strip-Mapping
Desmodulação
de quadratura
Chirp rate:
Compressão
em distância
Compressão
em azimute
IST/IT

4 v 2

r0
Duração:
T
T
r0
Dv
15
ESCOLA NAVAL - 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)
IST/IT
2- Spaceborn com movimento
3-Airborn sem movimento
16
ESCOLA NAVAL
- Metodologia 2D
Recontrução de Frente de Onda
s(u, t )   f ( x' , y ' )a( x' , y 'u ) p(t  2r / c)dx' dy'
A
r  x'2 ( y'u ) 2
TF(t ) [s(u, t )]


s (u ,  )  P( )  f ( x' , y ' )a ( x' , y 'u ,  ) exp  j 2k x'2 ( y 'u ) 2 dx' dy'
A
k
IST/IT

c
17
ESCOLA NAVAL - Reconstrução de frente de onda (cont)
TF(u ) [s(u, )]
(Aproximação de fase estacionária)
|
{z
}
2
s(ku ,  )  P( ) A(ku ,  )  f ( x' , y' ) exp  j 2k 4k 2  ku x' jk y y' dx' dy'


A
2

 K x  4k 2  k u


 K y  ku
s(ku , )  P( ) A(ku , )F (kx , k y )
IST/IT
Transformada de
Fourier de f(x,y)
18
ESCOLA NAVAL - Reconstrução de frente de onda (cont)
s(ku , )  P( ) A(ku , )F (kx , k y )
Estimador de F
Fˆ (k x , k y ) 
s ( ku ,  )
P( ) A(ku ,  )
• Filtro Inverso
•Filtro Adaptado
Fˆ (k x , k y )  s(ku ,  ) P ( ) A (ku ,  )
IST/IT
19
ESCOLA NAVAL - Reconstrução de frente de onda (cont)
Resposta Impulsional do filtro Adaptado
h(k x , k y )  P ( ) A(ku ,  )
A(2k sin  (u), )  a( x' , u,  )
2
| A (ku ,  ) |2
| P ( ) |2
 kM
kM
kM  2k sin( A / 2)
h(t )

2
T
kM
t
IST/IT
ku


D

k A 2
u
20
ESCOLA NAVAL
Modo de Operação ISAR
• Alvo move-se
• Sensor estático
ou móvel
Complexidade adicional
•Parametros de trajectória
IST/IT
21
ESCOLA NAVAL - Exemplo de ISAR
GAIR – GROUND-TO-AIR IMAGING RADAR PROCESSOR
Funções
•Seguimento de
aeronaves a baixa
altitude
• Adquirir dados ISAR
• Processar os dados
adquiridos e gerar imagens
de alta resolução
IST/IT
22
ESCOLA NAVAL
Parâmetros de operação
IST/IT
23
ESCOLA NAVAL
Alvo
IST/IT
Aspect angle –40º
Grazing angle –23º
24
ESCOLA NAVAL
IST/IT
25
ESCOLA NAVAL
-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
26
ESCOLA NAVAL
Interferometria SAR
x1  ρ1e i 1
i 2
x2  ρ2e
S1
_
B
R1
S2
2π
i 
2 Ri  ψi
λ
h  f 1( 12 ) .
R2
h(x)
1  2 
IST/IT
4

R1  R2   ( 1  2 ) .
27
Absolute
ESCOLA
NAVALPhase estimation in InSAR (Interferometric SAR)
InSAR Problem: Estimate 2- 1 from signal read by
28
s1 IST/IT
and s2
ESCOLA NAVAL
Dificuldades
 É apenas possível determinar módulo 2 e i 1 é
uma função periódica de período 2. 
 Os dados radar estão imersos em ruído (por
exemplo, ruído térmico)
1   2
 O termo
não é nulo devido a descorrelação
dos tipos geométrico, temporal e espacial.
IST/IT
29
ESCOLA NAVAL
Geometria InSAR
2 R 1 B  sin ( 0   )  B2 2 R 1 B cos ( 0   ) 
R 

.
2R
2R
 


ΔR 0
R
S2
_
h
S2
S1
_
B



'
S1


R1

H
b
a
IST/IT
R2

'

B



R 0  B sin ( 0  )

 

R h  B cos (0  ) .


B||
B
h
c
Rt
h
  
h
R 1   sin .
B
 4 
30
ESCOLA NAVAL
Registo
I2 (N,N)
I
Rˆ 12 ( k, l)   I1 (i  k, j  l)I*2 (i, j)
2
I1
I1 (N,N)
y
i, jJ
y
i+1
9
8
7
6
5
4
3
2
1
I2 (k,l)
I
i
2
I2 (1,1)
j
1 2 34 5 6 78 9
j+1
x
( kˆ , ˆl )  arg max Rˆ ( k, l).
k, l
I1 (5,5)
I1 (1,1)
a
x
b
N
ˆk  1  kˆ
i
N J i 1
IST/IT
N
ˆl  1  ˆl
i
N J i 1
31
ESCOLA NAVAL
Exemplo
Imagem de amplitude
IST/IT
Disparidades da imagem reamostrada
32
ESCOLA NAVAL
Interferograma
Coerência
N
*
 x i1 x i 2
ˆ ML 
i 1
N
 x i1
i 1
IST/IT
2
N
  x i2
2
i 1
33
ESCOLA NAVAL
Histogramas da coerência
Após registo grosso
IST/IT
Após registo fino
34
ESCOLA NAVAL
IST/IT
Results
35
ESCOLA NAVAL Results
Estimated surface (15:15:120)+2
} truth surface (15:15:120)
IST/IT
36
ESCOLA NAVAL
Results
Original
IST/IT
Estimated
37
ESCOLA NAVAL
Outro Exemplo
IST/IT
38
ESCOLA NAVAL
IST/IT
39