SAR Instrument Principles and
Processing
J.A. Johannessen and F. Collard
with support from
B. Chapron, R. Romeiser and W. Alpers
Content
• SAR antenna
• Frequency domain – EM spectrum
• Incidence angles and Ranges
• Surface roughness and backscatter (NRCS)
• Bragg scattering and image formation
• Spatial Resolution
– Range (across track)
– Azimuth (along track)
• Processing from raw data to SAR image
• Range Doppler signal
ENVISAT ASAR Antenna
AATSR
MIPAS
SCIAMACHY
MERIS
ASAR antenna
MWR
10 m x 1.33 m
GOMOS
DORIS
RA-2
LRR
ASAR
Synthetic Aperture Radar
• Side-looking
• Active antenna that
transmits/receives electromagnetic
radiation in VV, HH, VH, HV pol
• Records both signal amplitude and
phase
• Works both day and night
• Can “see” through clouds
7 October 2011 Morten Wergeland Hansen 4
Electromagnetic Spectrum
l= c * T = c * 1/f
SAR f-range
Angles and Ranges
Surface Roughness and SAR backscatter
The ocean surface roughness is influenced by wind
and waves, currents, surface slicks and sea ice
and is often different in open ocean versus coastal or
ice covered regions due to fetch effects
The surface roughness is the source for the backscatter of
the SAR signal.
The signal that arrives at the antenna is registered both in
Amplitude and Phase.
Surface Roughness and SAR backscatter
Surface Roughness
&
Surface Motion
Wind, Waves
SAR
and Current, Sea
NRCS & Doppler
ice, natural film, Shift
oil spill
Surface Roughness and SAR backscatter
ASAR, 2.10.2011
Smooth surface (light
wind)
Enhanced
backscatter
Rough surface (stronger wind)
Bragg Scattering
At oblique incidence angles, the SAR
backscattering arising from the sea surface is
caused by surface waves of the order of the
radar wavelength.
These waves are called “ Bragg waves“. They obey
the “ Bragg resonance condition“:
lB = lr /2 sin
lB
where lB = Bragg wavelength, lr = radar wavelength, and
= incidence angle
Bragg Scattering and Radar Cross-section
SYNTHETIC The radar backscatter is primarily determined by the Bragg
APERTURE scattering (determined by the11surface roughness): l = 2 lB sin
RADAR for incidence angles in the range of 20 to 50 degrees
Smooth vs. Rough Surfaces
Specular
Backscatter
from rough
surface b >
from smooth Bragg
surface a
lB
SAR is a transmitting-receiving
From [Link] instrument where
SO431: Environmental Remote Sensing 13.10.2003
Pr = (Pt / 4pR2) G ( / 4pR2) A
P (r=receive, t=transmit), R = range distance, G = antenna gain, = radar cross section, A = antenna area,
= 10 log10 ( /A) is defined as radar backscatter (=function of surface roughness)
Image Formation of Surface Roughness
z A short pulse is emitted by the antenna and then the
amplitude and phase of the backscattered signal is
recorded as a function of time.
Antenna This is repeated over again while the platform is
moving, and a 2-dimensional image is thus generated.
Azimuth
Range
0
x
0() = 0 (surface roughness) = y
0(wind, waves, breaking, current)
Spatial Resolution
The SAR spatial resolution is independent of the platform
height. This is remarkable and unique for SAR
instruments.
In range Xr = c t/2sin (where c: speed of light, t: pulse length, incidence angle)
In azimuth Xa = D/2 (D = antenna length)
A frequency modulation - chirp of the pulse in used in range. In azimuth a synthetic aperture
principle is used whereby the motion of the platform induces a frequency modulation. In
result a very long antenna is synthesized by the motion of the platform. For ASAR the length
of the synthetic antenna is around 20 km!
Consequently, in both directions, the signals are frequency modulated. This modulation is
sensed in the Doppler shift.
Range Resolution
Increase of the range resolution
We have noticed already that the range resolution Xr = c t/2 sinq of a RAR or SAR is
independent of the platform height.
However, it is technically not possible to generate a radar pulse that has a
length of only a few meters.
Radar engineers use a long pulse with a (linearly) modulated frequency - called
a Chirp. With this technique it is possible to increase the range resolution
Range Resolution
Pulse is not frequency modulated
L = ct Target 1 Target 2
L/2
L=pulse length L/2 = separation of the targets
The two targets separated by L/2 can only be resolved when the pulse length ct is equal to or
smaller smaller than L/2.
Range resolution: Xr = ct/2sinq
Range Resolution
Improvement of the range resolution by
using a frequency modulated pulse
S2(t)
S1(t)
T1 T2
Backscattered signal S(t), is the
sum of the backcattered signals
from target 1 and target 2.
S(t) = S1(t) +S2(t)
Range Resolution
Improvement of the range resolution by
using a frequency modulated pulse
Uref S(t)
C(t)
(values are normalized)
The positions of the two targets show up in the correlation function c(t) as two
separate peaks. The minimum distance of these peaks is given by Xr = ct/2
Range Resolution
A frequency modulated pulse can resolve targets
which are separated by less than L/2.
This is achieved by cross-correlating the
backscatterd pulse s(t) = s1(t) + s2(t) with a
reference signal uref (t), which is the
complex conjugate of the emitted signal:
Azimuth Resolution
Deriving fine azimuth resolution (courtesy of Prof. Werner Alpers)
Acoustic anology of a
antenna
SAR
frequency
Azimuth Resolution
Synthetic aperture radar
V
V = platform velocity
antenna
T = integration time beam
L = length of the synthetic antenna
target L
The target is for T seconds (T = L/V) in the antenna beam
Azimuth Resolution
Synthetic aperture radar principle (more after Alpers)
V antenna
target V
fd = + Vq/ l fd = Doppler shift
Azimuth Resolution
Synthetic aperture radar principle - 2
V antenna
target V
fd = 0
Azimuth Resolution
Synthetic aperture radar principle - 3
V antenna
target V
fd = - Vq/ l
Azimuth Resolution
Change of the Doppler shift fd across the aperture
antenna
V= velocity of the target
through the antenna beam
q/2
Vr=Vsinq/2
=V q/2 Vr
(approx.) V q/2
fd = +2Vr/l (= (2Vr/c) f ) =+Vq/ l
Change of Doppler shift across the aperture = fd -(- fd)=2fd = 4Vr/ l= 2Vq/ l
2fd = B is called the Azimuthal Bandwidth of the SAR
Azimuth Resolution
Application of the general rule in signal processing that, if
an electrical system has a bandwidth B, then it can resolve a signal that has a time length
of t=1/B, to SAR:
The time interval that can be resolved is
t=1/B=1/2fd= l /2Vq =D/2V (because of q = l/D).
The spatial interval in flight direction that can be
resolved =azimuthal resolution = Xa = Vt =D/2.
Thus, the unique azimuthal resolution of a SAR is
independent of range R and is proportional to the
antenna length D
Schematic Summary of the SAR Image azimuth resolution
Processing
In the SAR processor these frequency modulated signals
are used to improve the resolutions in range Xr and
azimuth Xa. This is called range compression and
azimuth compression, respectively.
Xa = D/2 Xr = c t/2sin = c/2Bsin
Therefore the SAR processor consists essentially of 2
correlators, one for range and one for azimuth.
Processing
STAGES IN SAR IMAGE COMPRESSION
(after Robinson, 2004)
Processing
Raw SAR data
Range compressed data
Range + azimuth
compressed data
= SAR image
SAR imaging coordinate system
The illuminated area can
be referenced to a coordinate
system of concentric circles
(equidistances) and coaxial
hyperbolas (equi-Doppler).
Each point in the image plane
can be uniquely identified by
its time delay and Doppler
shift.
Speckle Noise and Removal
Averaging
Rhine plume
reduces
front
noise
signature,
subimage size = 6 km × 6 km
Image Formation and Doppler anomalies
z
APPROACH: SAR Imaging of Roughness and Doppler Shift
Antenna
0() = 0 (surface roughness) =
Azimuth
0(wind, waves, breaking, current
Range
0
x VD sinVD
y
VD = VD (wind driven scatterer, waves, breaking, current)
Doppler Centroid Anomaly
f Dca = f Dc - f Dp
f Dc : estimated Doppler centroid
frequency shift
Power f Dp : predicted Doppler shift
Chapron et al. (2003, 2005)
Doppler centroid Time
Frequency
Target moving
Target at rest
Time
7 October 2011 Morten Wergeland Hansen 33
Doppler Processing
Envisat ASAR scene off the South African coast
14 September 2010, 21:15 UTC
Doppler centroid
Observed Predicted anomaly
Land
7 October 2011 Morten Wergeland Hansen 34
Azimuth Bias Correction
Before correction After correction
7 October 2011 Morten Wergeland Hansen 35
Range Bias Correction
Measurements over
land, where the
Doppler shift should be
zero
Ideal antenna pointing
Incidence angle, q [degrees]
7 October 2011 Morten Wergeland Hansen 36
Range Bias Correction
Before correction After correction
Land
7 October 2011 Morten Wergeland Hansen 37
Retrieval Accuracy
24.5 cm/s
Incidence angle, q [degrees]
7 October 2011 Morten Wergeland Hansen 38
SAR backscatter and Doppler shift anomaly
The ocean surface roughness is the sources of the SAR
backscatter.
The SAR therefore sense the near surface wind and waves,
currents, surface slicks and sea ice.
The SAR is also registering the Doppler centroid anomaly.
The challenge is to convert and partition the SAR
backscatter signal and Doppler centroid anomaly reliable
estimate of wind, waves and current.
Outlook
• Sentinel-1A will be
launched in early 2014
• The range Doppler shift is
planned to become a
standard product, with
– Significantly better
accuracies
– Improved capability to
monitor the temporal and
spatial variability of the
ocean surface circulation
from SAR
7 October 2011 Morten Wergeland Hansen 40
