# Multiband Microwave Imaging Analysis of Ionosphere and Troposphere Refraction for Spaceborne SAR

Authors: Fan Zhang; Guojun Li; Wei Li; Wei Hu
Journal: International Journal of Antennas and Propagation
 (2014)

---

## Abstract

Ionosphere has different stratification at the different height. Troposphere has different refractivity at the different height. When microwave signals transmit through the ionosphere and the troposphere, the real propagation path is not an ideal straight line, but a slightly curved straight line. For the synthetic aperture radar (SAR) system, the actual distance errors will result in phase errors, which impact range section and azimuth section of SAR raw data. Consequently, the imaging precision has been decreased by imprecise slant range history. In this paper, we simulate the propagation path between satellite and the target according to Snell's law and analyze how the ionospheric and tropospheric refraction impact the spaceborne SAR imaging performance at L-band and X-band. The simulation results show that the two refraction effects should be compensated in low frequency band for better image focusing performance.


---

## Body

## 1. Introduction

As the microwave signals transmit between SAR and targets, the imaging results are inevitably affected by the ionosphere and troposphere. The ionizing is at the height of 70 km to 500 km from the ground. Low frequency signal often suffers ionosphere influences severely, such as dispersion, Faraday rotation, refraction, rotation of polarization caused by irregularities, and so on [1–4]. For example, the ionized ionosphere induces Faraday rotation (FR) that affects radar polarization and causes signal path delays, and all of these can be ignored at high frequencies (10 GHz) [5]. These effects all have strong relationship with the TEC (total electron content). Observations show that the electron density of the ionosphere presents a structure of stratification. Generally, we call the electron concentration distribution with height as the electron density profile. In recent years, there has been an increasing interest in the analysis of the impact of propagation link on low frequency and wide-bandwidth spaceborne SAR imaging [5]. Chen et al. put forward an operation mode for topside ionospheric sounding based on spaceborne high frequency SAR (HF-SAR) [6]. Theoretical analysis indicates that the azimuth resolution can be improved to tens of meters. Liu et al. develop a numerical model to investigate the SAR image degradation caused by an inhomogeneous ionosphere [7]. Both horizontal and vertical structures of the ionosphere are considered in this model. Although they do not study a method to compensate for the ionospheric effects on SAR, the numerical model they have developed is useful for testing different mitigation techniques. And they develop two numerical models about the ionosphere acting on the propagation characteristics of the spaceborne SAR [8]. One of the models simulates the ray-bending effects on the SAR system. Their results show substantial image shift due to this effect atP-band or lower frequency band. Wang et al. prensent the analytical study and numerical simulation to investigate the ionospheric effects on SAR imaging with Gaussian white noise atP-band [9]. They use the IRI 2001 model to generate the electron density profile and use the two-parameter spectrum to generate the horizontal electron density profile, which are also used in this paper.Troposphere is under the altitude of about 12 km. It is the closest layer to the ground and a mixture of kinds of gases (nitrogen, oxygen, carbon dioxide, and so on) and water vapor. The dielectric properties of the troposphere vary with time and space. So the microwave signals transmit in the troposphere differently from in the vacuum [10]. The effect of troposphere on SAR imaging is usually ignored for that the early moderate resolution imaging would not be affected by the relative minor tropospheric path error. Consequently, there are few thorough analysis and validation related to the troposphere in SAR field. Currently, the spaceborne SAR imaging resolution has reached half a meter level. So the troposphere environmental factor has become one of the factors that cannot be ignored on spaceborne SAR imaging with high resolution. Sun et al. introduce the impact of atmospheric refraction on the high resolution airborne SAR which mainly reflects on azimuth resolution [10]. The changes of atmosphere refractive index in troposphere with the Hopfield refractive model are analyzed, and the calculation of the apparent range between radar and ground targets is given. The algorithm mentioned in paper has certain reference value to the performance and imaging analysis of the high resolution airborne SAR system which is level of centimeter. The refraction effects are addressed with assuming a horizontally nonhomogeneous troposphere [11]. The impact of the range refractivity on a terrestrial radio link is investigated in terms of radio horizon range. By acquiring simultaneous weather radar data over the test site, it is possible to flag affected SAR images and exclude them from the procedure.Through the above introduction, the index of refraction will vary with the height both in ionosphere and troposphere [12]. According to Snell’s law, the propagation of radio waves in the multilayer medium will be bent and not straight as we expected. The path error will cause many problems in SAR imaging. Therefore, we will simulate and analyze the impact of the multiband microwave propagation errors caused by the ionosphere and troposphere on SAR imaging in the following contents.In this paper, we focus on the propagation path’s bending caused by ionosphere and troposphere refraction, which are seldom synthetically considered before. First, we discuss the effects of ionospheric and tropospheric refraction on SAR imaging. In Section2, we introduce the background of the ionospheric and tropospheric refractive model. The structure of ionosphere can be analyzed from two aspects: vertical profile and horizontal profile. In Section 3, we use the path tracing method to calculate the actual distance according to Snell’s law. In Section 4, we perform SAR imaging simulation with the path errors and analyze the imaging performance. In Section 5, we draw a conclusion and put forward some plans for future research.


## 2. Atmosphere Structure

As Figure1 shows, troposphere is at the height of 0 to 12 km. The layer between 60 km and 90 km is the D region and from 90 km to 140 km is the E region. The F region is at the height above 140 km. In Figure 1, curve ① represents the profile of vertical electron density; curve ② represents the profile of horizontal electron density. Because the sketch map is just used to indicate the atmosphere structure, the vertical and horizontal electron density profile are rough. We can see that the electron concentration varies with different heights and different horizontal distribution, which will lead to the change of the ionosphere refraction. Moreover, there are kinds of gases in the troposphere, and the tropospheric index of refraction will vary with different heights as well. So there must be path errors in the propagation of microwaves.Figure 1
Sketch map of ionosphere and troposphere.


### 2.1. Basic Structure of Ionosphere

The basic structure of ionosphere can be discussed from two aspects: vertical TEC profile and horizontal TEC profile.


#### 2.1.1. Vertical TEC Profile

According to Appleton-Hatreek, we can get the simplified formula about the refractive index of ionosphere [13, 14]:
(1)
n
=
1
-
80.6
N
e
f
2
,
where N
e is the electron concentration and f is the wave frequency in Hz.In this paper, we use the engineering ionospheric electron density profile model IRI 2001 to calculateN
e [15, 16]. The model is expressed as below [17]:
(2)
N
e
(
h
)
=
{
N
m
E
[
1
-
(
h
m
E
-
h
y
m
E
)
2
]
,
h
m
E
-
y
m
E
≤
h
≤
h
m
E
N
j
-
N
m
E
h
j
-
h
m
E
h
+
(
N
m
E
)
h
j
-
N
j
(
h
m
E
)
h
j
-
h
m
E
,
h
m
E
<
h
≤
h
j
N
m
F
2
[
1
-
(
h
m
F
2
-
h
y
m
F
2
)
2
]
,
h
j
<
h
<
h
m
F
2
N
m
F
2
exp
⁡
[
1
2
(
1
-
h
-
h
m
F
2
H
-
e
-
(
h
-
h
m
F
2
)
/
H
)
]
,
h
m
F
2
<
h
≤
1000
km
,
where N
m is the largest electronic concentration of each layer, N
j is for F
1 layer, h is the height corresponding to the largest electronic concentration of each layer, and h
j is for F
1 layer. y
m is the half-thickness for each layer, f
0 is the critical frequency for each layer, and f
j is for F
1 layer. h
m
E is 115 km; y
m
E is 20 km [18].The vertical electron concentration is shown in Figure2.Figure 2
Vertical electron density profile model.


#### 2.1.2. Horizontal TEC Profile

As shown in Figure3, the horizontal TEC varies with the different latitudes and longitudes. Since the spaceborne SAR aperture size is at least tens of kilometers, the horizontal TEC gradient will impact the SAR imaging performance.Figure 3
Horizontal electron density profile model.In this paper, we use the two-parameter spectrum to generate the horizontal profile of electron concentration [19, 20]. V
(
k
) can be expressed as [7]
(3)
V
(
k
)
=
{
π
A
v
1
-
1
1
(
k
2
+
k
0
2
)
v
1
-
1
-
π
A
(
k
b
2
+
k
0
2
)
v
1
-
1
(
1
v
1
-
1
-
1
v
2
-
1
)
,
k
≤
k
b
,
π
A
(
k
b
2
+
k
0
2
)
v
2
-
v
1
v
2
-
1
1
(
k
2
+
k
0
2
)
v
2
-
1
,
k
>
k
b
,
where k
b
=
2
π
/
L
b, k
0
=
2
π
/
L
0, k
b is the interrupt wave number, L
0 = 10 km, L
b = 500 m, and 2
v
1 = 3.5, 2
v
2 = 5.5.To generate the horizontal profile of electron concentration, a sequence of normally distributed random numbers is used in the phase of the profile spectrum. Using the discrete Fourier transform, the relationship between the electron density profile and the profile spectrum is described in the following equations:(4)
f
(
x
)
=
1
L
∑
n
=
N
/
2
N
/
2
-
1
F
(
k
n
)
exp
⁡
(
i
k
n
x
)
,
where
(5)
F
(
k
n
)
=
2
π
L
V
(
k
n
)
{
1
2
[
N
(
0,1
)
+
i
N
(
0,1
)
]
,
n
≠
0
,
N
2
N
(
0,1
)
,
n
=
0
,
N
2
and K
n
=
2
π
n
/
L, L is a total length of a profile. N
(
0,1
) denotes a sequence of normally distributed numbers in [
0,1
] with zero mean and unity standard deviation. The horizontal electron concentration is shown in Figure 4.Figure 4
The propagation path of microwave signals.


### 2.2. Model for Tropospheric Refractivity

As for troposphere, we often use the statistical model of the refractive index. In this paper, we use the Hopfield model [12], which includes mainly two parts. The first one is the so-called wet part, which is related to the concentration of water vapor. The other part is the so-called dry part, which is related to the gaseous nature of the lower part of the atmosphere. The dry parts and the wet parts are in form of function of fourth power [21]. 								(6)
N
(
h
)
=
N
d
+
N
w
(7)
N
d
=
{
N
0
d
(
H
d
-
h
H
d
-
h
0
)
4
,
h
<
H
d
0
,
h
≥
H
d
(8)
N
w
=
{
N
0
w
(
H
w
-
h
H
w
-
h
0
)
4
,
h
<
H
w
0
,
h
≥
H
w
,
where the subscripts d and w represent dry and wet parts; N
0
d and N
0
w are the ground atmospheric refractivity; H
d and H
w are the height where the dry and wet section of atmospheric refractivity decay to 0; h
0 is the height of ground.Then, we can calculate the ground atmospheric refractivity as below:(9)
N
0
d
=
77.6
P
0
(
t
s
+
273.15
)
N
0
w
=
3.37
×
1
0
3
R
H
×
a
(
t
s
+
237.15
)
2
exp
⁡
(
b
t
s
c
+
t
s
)
.Then, we can calculate the refractive index of troposphere by (6).


## 2.1. Basic Structure of Ionosphere

The basic structure of ionosphere can be discussed from two aspects: vertical TEC profile and horizontal TEC profile.


### 2.1.1. Vertical TEC Profile

According to Appleton-Hatreek, we can get the simplified formula about the refractive index of ionosphere [13, 14]:
(1)
n
=
1
-
80.6
N
e
f
2
,
where N
e is the electron concentration and f is the wave frequency in Hz.In this paper, we use the engineering ionospheric electron density profile model IRI 2001 to calculateN
e [15, 16]. The model is expressed as below [17]:
(2)
N
e
(
h
)
=
{
N
m
E
[
1
-
(
h
m
E
-
h
y
m
E
)
2
]
,
h
m
E
-
y
m
E
≤
h
≤
h
m
E
N
j
-
N
m
E
h
j
-
h
m
E
h
+
(
N
m
E
)
h
j
-
N
j
(
h
m
E
)
h
j
-
h
m
E
,
h
m
E
<
h
≤
h
j
N
m
F
2
[
1
-
(
h
m
F
2
-
h
y
m
F
2
)
2
]
,
h
j
<
h
<
h
m
F
2
N
m
F
2
exp
⁡
[
1
2
(
1
-
h
-
h
m
F
2
H
-
e
-
(
h
-
h
m
F
2
)
/
H
)
]
,
h
m
F
2
<
h
≤
1000
km
,
where N
m is the largest electronic concentration of each layer, N
j is for F
1 layer, h is the height corresponding to the largest electronic concentration of each layer, and h
j is for F
1 layer. y
m is the half-thickness for each layer, f
0 is the critical frequency for each layer, and f
j is for F
1 layer. h
m
E is 115 km; y
m
E is 20 km [18].The vertical electron concentration is shown in Figure2.Figure 2
Vertical electron density profile model.


### 2.1.2. Horizontal TEC Profile

As shown in Figure3, the horizontal TEC varies with the different latitudes and longitudes. Since the spaceborne SAR aperture size is at least tens of kilometers, the horizontal TEC gradient will impact the SAR imaging performance.Figure 3
Horizontal electron density profile model.In this paper, we use the two-parameter spectrum to generate the horizontal profile of electron concentration [19, 20]. V
(
k
) can be expressed as [7]
(3)
V
(
k
)
=
{
π
A
v
1
-
1
1
(
k
2
+
k
0
2
)
v
1
-
1
-
π
A
(
k
b
2
+
k
0
2
)
v
1
-
1
(
1
v
1
-
1
-
1
v
2
-
1
)
,
k
≤
k
b
,
π
A
(
k
b
2
+
k
0
2
)
v
2
-
v
1
v
2
-
1
1
(
k
2
+
k
0
2
)
v
2
-
1
,
k
>
k
b
,
where k
b
=
2
π
/
L
b, k
0
=
2
π
/
L
0, k
b is the interrupt wave number, L
0 = 10 km, L
b = 500 m, and 2
v
1 = 3.5, 2
v
2 = 5.5.To generate the horizontal profile of electron concentration, a sequence of normally distributed random numbers is used in the phase of the profile spectrum. Using the discrete Fourier transform, the relationship between the electron density profile and the profile spectrum is described in the following equations:(4)
f
(
x
)
=
1
L
∑
n
=
N
/
2
N
/
2
-
1
F
(
k
n
)
exp
⁡
(
i
k
n
x
)
,
where
(5)
F
(
k
n
)
=
2
π
L
V
(
k
n
)
{
1
2
[
N
(
0,1
)
+
i
N
(
0,1
)
]
,
n
≠
0
,
N
2
N
(
0,1
)
,
n
=
0
,
N
2
and K
n
=
2
π
n
/
L, L is a total length of a profile. N
(
0,1
) denotes a sequence of normally distributed numbers in [
0,1
] with zero mean and unity standard deviation. The horizontal electron concentration is shown in Figure 4.Figure 4
The propagation path of microwave signals.


## 2.1.1. Vertical TEC Profile

According to Appleton-Hatreek, we can get the simplified formula about the refractive index of ionosphere [13, 14]:
(1)
n
=
1
-
80.6
N
e
f
2
,
where N
e is the electron concentration and f is the wave frequency in Hz.In this paper, we use the engineering ionospheric electron density profile model IRI 2001 to calculateN
e [15, 16]. The model is expressed as below [17]:
(2)
N
e
(
h
)
=
{
N
m
E
[
1
-
(
h
m
E
-
h
y
m
E
)
2
]
,
h
m
E
-
y
m
E
≤
h
≤
h
m
E
N
j
-
N
m
E
h
j
-
h
m
E
h
+
(
N
m
E
)
h
j
-
N
j
(
h
m
E
)
h
j
-
h
m
E
,
h
m
E
<
h
≤
h
j
N
m
F
2
[
1
-
(
h
m
F
2
-
h
y
m
F
2
)
2
]
,
h
j
<
h
<
h
m
F
2
N
m
F
2
exp
⁡
[
1
2
(
1
-
h
-
h
m
F
2
H
-
e
-
(
h
-
h
m
F
2
)
/
H
)
]
,
h
m
F
2
<
h
≤
1000
km
,
where N
m is the largest electronic concentration of each layer, N
j is for F
1 layer, h is the height corresponding to the largest electronic concentration of each layer, and h
j is for F
1 layer. y
m is the half-thickness for each layer, f
0 is the critical frequency for each layer, and f
j is for F
1 layer. h
m
E is 115 km; y
m
E is 20 km [18].The vertical electron concentration is shown in Figure2.Figure 2
Vertical electron density profile model.


## 2.1.2. Horizontal TEC Profile

As shown in Figure3, the horizontal TEC varies with the different latitudes and longitudes. Since the spaceborne SAR aperture size is at least tens of kilometers, the horizontal TEC gradient will impact the SAR imaging performance.Figure 3
Horizontal electron density profile model.In this paper, we use the two-parameter spectrum to generate the horizontal profile of electron concentration [19, 20]. V
(
k
) can be expressed as [7]
(3)
V
(
k
)
=
{
π
A
v
1
-
1
1
(
k
2
+
k
0
2
)
v
1
-
1
-
π
A
(
k
b
2
+
k
0
2
)
v
1
-
1
(
1
v
1
-
1
-
1
v
2
-
1
)
,
k
≤
k
b
,
π
A
(
k
b
2
+
k
0
2
)
v
2
-
v
1
v
2
-
1
1
(
k
2
+
k
0
2
)
v
2
-
1
,
k
>
k
b
,
where k
b
=
2
π
/
L
b, k
0
=
2
π
/
L
0, k
b is the interrupt wave number, L
0 = 10 km, L
b = 500 m, and 2
v
1 = 3.5, 2
v
2 = 5.5.To generate the horizontal profile of electron concentration, a sequence of normally distributed random numbers is used in the phase of the profile spectrum. Using the discrete Fourier transform, the relationship between the electron density profile and the profile spectrum is described in the following equations:(4)
f
(
x
)
=
1
L
∑
n
=
N
/
2
N
/
2
-
1
F
(
k
n
)
exp
⁡
(
i
k
n
x
)
,
where
(5)
F
(
k
n
)
=
2
π
L
V
(
k
n
)
{
1
2
[
N
(
0,1
)
+
i
N
(
0,1
)
]
,
n
≠
0
,
N
2
N
(
0,1
)
,
n
=
0
,
N
2
and K
n
=
2
π
n
/
L, L is a total length of a profile. N
(
0,1
) denotes a sequence of normally distributed numbers in [
0,1
] with zero mean and unity standard deviation. The horizontal electron concentration is shown in Figure 4.Figure 4
The propagation path of microwave signals.


## 2.2. Model for Tropospheric Refractivity

As for troposphere, we often use the statistical model of the refractive index. In this paper, we use the Hopfield model [12], which includes mainly two parts. The first one is the so-called wet part, which is related to the concentration of water vapor. The other part is the so-called dry part, which is related to the gaseous nature of the lower part of the atmosphere. The dry parts and the wet parts are in form of function of fourth power [21]. 								(6)
N
(
h
)
=
N
d
+
N
w
(7)
N
d
=
{
N
0
d
(
H
d
-
h
H
d
-
h
0
)
4
,
h
<
H
d
0
,
h
≥
H
d
(8)
N
w
=
{
N
0
w
(
H
w
-
h
H
w
-
h
0
)
4
,
h
<
H
w
0
,
h
≥
H
w
,
where the subscripts d and w represent dry and wet parts; N
0
d and N
0
w are the ground atmospheric refractivity; H
d and H
w are the height where the dry and wet section of atmospheric refractivity decay to 0; h
0 is the height of ground.Then, we can calculate the ground atmospheric refractivity as below:(9)
N
0
d
=
77.6
P
0
(
t
s
+
273.15
)
N
0
w
=
3.37
×
1
0
3
R
H
×
a
(
t
s
+
237.15
)
2
exp
⁡
(
b
t
s
c
+
t
s
)
.Then, we can calculate the refractive index of troposphere by (6).


## 3. Calculation of the Path Error

When the microwave signals transmit through the ionosphere and troposphere, the path will become bent because of refraction as Figure5 shows. We can calculate the actual distance according to Snell’s law. Assuming that the radar situated at a height of h, radar’s depression angle is θ
0, the refraction index at the height of h is n
0. Then, the actual distance of the radar and the target is
(10)
R

=

∫
0
h
n
2
(
z
)
n
2
(
z
)
-
n
0
2
cos
⁡
2
θ
0
d
z
.Figure 5
Path tracing course.In (10), n
(
z
) is the refractive index of space, including ionosphere and troposphere. Due to the fact that different models are used, we independently calculate the actual distance in ionosphere and troposphere refraction. Furthermore, the ionospheric refraction includes two parts: vertical TEC section and horizontal TEC section.Accordingly, the horizontal distance between radar and the target isy
0:
(11)
y
0
=
∫
0
h
n
0
cos
⁡
θ
0
n
2
(
z
)
-
n
0
2
cos
⁡
2
θ
0
d
z
.According the equations above, we use the method of path tracing to calculate the actual distance of radar and the target. First, we substitutey
0 into (11) as a known quantity; thus we can get radar’s depression angle θ
0. Then, we substitute the θ
0 into (10), and get the actual distance of radar and the target. The flowchart of the course of path tracing is expressed in Figure 6.Figure 6
SAR imaging geometry.SAR imaging geometry is shown as Figure7. Spaceborne SAR is at the height of h and moves at a speed of v
a along the x-axis direction; the instant time is t
m. When it comes to the synthetic aperture center, t
m = 0. The shadow is the footprint, the horizontal width of the beam is  β, and the vertical width is  φ. At the t
m moment, the actual distance between SAR and the target is  R
(
t
m
).Figure 7
An example of horizontal TEC map over China and adjacent regions.If the ionospheric and tropospheric environment can be ignored, the actual distance is equal to the theoretical distance between SAR and the targets, which can be calculated as [22]
(12)
R
(
t
m
)
=
h
2
+
r
0
2
+
(
v
a
t
m
)
2
.So the path error is(13)
Δ
R
(
t
m
)
=
R
(
t
m
,
n
(
z
)
)
-
R
(
t
m
)
.The phase error caused by path error is(14)
Δ
φ
(
t
m
)
=
-
4
π
Δ
R
(
t
m
)
λ
.Assuming that the total synthetic aperture time isT
a, the time of range section is τ and the radar signal is chirp signal. We can get the echo that the radar receives from the target as below:
(15)
s
(
t
m
,
τ
)
=
σ
w
a
2
(
θ
)
rect
(
τ
-
2
R
(
t
m
)
/
c
T
a
)
·
exp
⁡
(
-
j
4
π
λ
R
(
t
m
)
)
·
exp
⁡
(
j
π
K
(
τ
-
2
R
(
t
m
)
c
)
2
)
,
where t is the azimuth slow time,   τ  is the range fast time,   K  is chirp rate, w
a
(
θ
) is the azimuth antenna illumination envelop,   σ  is the backscatter coefficient,   c is the speed of light, and  R is the distance from antenna phase center to the target.The echo signals that the radar receives from the target are related to the actual distance between the radar and the target. So the path error will lead to the delay and the offset error in range section.On the other hand, we can do quadratic Taylor expansion on the actual distanceR
(
t
m
,
n
(
z
)
) around t
c
=
r
0
/
v
a:
(16)
R
(
t
m
,
n
(
z
)
)
≈
R
(
(
t
m
-
t
c
)
,
n
(
z
)
)
+
Δ
t
*
R
′
(
(
t
m
-
t
c
)
,
n
(
z
)
)
+
0.5
Δ
t
2
*
R
c
′′
(
(
t
m
-
t
c
)
,
n
(
z
)
)
2
.In (16), there is a quadratic term. There must be quadratic phase error in azimuth section when the distance errors are considered. Obviously, the quadratic phase error will lead to the elevated sidelobe energy, thus affecting the SAR image quality. In summary, the path error will lead to two problems in SAR imaging: the offset error in range section and quadratic phase error in azimuth section, which can be illustrated as Figure 8.Figure 8
Sketch map of the offset error imaging in range section and quadratic phase error in azimuth section.


## 4. Multiband Spaceborne SAR Imaging Simulation with Refraction Path Error

To analyze the refraction effect of atmosphere on spaceborne SAR imaging, the propagation path errors should be introduced into echo generation process in every transmitted pulse. On the other hand, the multiband responses with atmospheric refraction effects are also the key topic, which can not only guide the course of spaceborne SAR system design, but also assist the in-orbit microwave imaging and analysis. Therefore, the multiband spaceborne SAR imaging simulation is carried out for atmospheric refraction error analysis and modeling.According to classical SAR simulation process, the imaging simulation with path error involves four steps, including refractive index calculation, propagation path error calculation, raw data simulation, and imaging analysis, as shown in Figure9. The simulation method is as follows.(i)
According to the description of Section2, we can calculate the different height’s refractive index in ionosphere and troposphere by the general engineering model. Then we change the microwave frequency to get the different refractive index distribution along the vertical height, which can be used in actual path calculation.
(ii)
Based on the path tracing method, the microwave propagation path under Snell’s law is calculated with considering the ionosphere and troposphere refraction. The specific algorithm is introduced in Section3. On this basis, we simulate the satellite orbit and calculate the path error at each pulse transmitting time and then get the error curve along the azimuth direction. To compare the multiband characteristic, it is assumed that the multiband simulation should be executed with same orbit and radar parameters. After multifrequency replacement, the multiband path error can be simulated.
(iii)
The spaceborne SAR raw data simulation includes four steps: target deployment, orbit calculation, path error injection, and raw data calculation. In order to measure the image quality, point targets are located in the simulated scene. Moreover, the time-domain algorithm that can easily introduce system errors is selected for spaceborne SAR raw data calculation. Through the simulation processing, the multiband SAR raw data can be generated for imaging analysis.
(iv)
The imaging analysis includes image processing and image quality assessment. First, the raw data should be processed by classical Chirp-Scaling imaging algorithm. Second, the point target data should be cut from SAR image for extracting the profile of range direction and azimuth direction. Finally, the SAR imaging index, such as resolution, expansion ratio, peak sidelobe ratio, and integral sidelobe ratio, can be calculated to analyze the multiband atmosphere refraction effects.Figure 9
Process of the simulation.


## 5. Simulation Result and Discussion

Simulations are carried out forL-band and X-band, which are applied for different remote sensing purpose and less discussed. The simulation parameters are shown in Table 1 and experimental conditions are described below:(1)
point target array is placed in the middle of simulated image;
(2)
in the course of refractive index calculation, the ionosphere height is 100~1000 km, the troposphere height is 0~12 km, and the vertical divided layer scale of ionosphere and troposphere is, respectively, 0.1 km and 0.001 km;
(3)
in horizontal electron density profile model, the average peak value is set to 1012 el/m3 and the perturbation is set to 10%. The simulated horizontal distance of electron density profile is 10 km with sampling interval of 100 m.Table 1
Simulation parameters.
Parameters
Value
Semimajor axis
7000 km
Eccentricity
0.0
Orbital inclination
95 degrees
RAAN
150 degrees
Argument of perigee
210 degrees
PRF
2000 Hz
Band width
60 MHz
Pulse width
30 us
Sampling rate
66.6 MHzFigure10 shows the microwave propagation path errors of L-band and X-band in a synthetic aperture time. With the change of slant range, the path error also varies in quadratic curve. If we do not take the horizontal TEC into consideration, the curve will be smooth. However, the curve shows a zigzag line by considering the refraction. From the results, the path error of L-band is greater than X-band, but the two bands have the similar fluctuation. Because ionosphere is higher than troposphere, the ionosphere contributes the majority of path errors for longer propagation path. It is seen that the higher frequency band radar works with the greater the refraction effects and worse image quality. Figure 11 depicts, respectively, the point target response profile of range direction and azimuth direction in L-band and X-band. The azimuth profile in L-band has a serious quadratic phase error, which leads to the sidelobe uplift. The azimuth profile in X-band is close to the ideal curve. Otherwise, both the two range profiles are very standard. We can see that the path errors will lead to the decline in image quality indicators of azimuth section. Tables 2 and 3 provide a quantitative proof of the conclusion.Table 2
Imaging indicators errors ofL-band and X-band compared to the theoretical indicators in range section.
Indicators
L
X
Resolution (m)
0.0327
0.0081
Expansion ratio (dB)
0.071
0.0017
Peak sidelobe ratio (dB)
0.1072
−0.0625
Integral sidelobe ratio (dB)
0.2509
−0.1326Table 3
Imaging indicators errors ofL-band to X-band compared to the theoretical indicators in azimuth section.
Indicators
L
X
Resolution (m)
1.9672
0.0479
Expansion ratio (dB)
0.4837
0.0118
Peak sidelobe ratio (dB)
−1.6258
−0.8447
Integral sidelobe ratio (dB)
4.6385
−0.5886The path error ofL-band and X-band.
(a)
L-band,f = 1.5 GHz
(b)
X-band,f = 9.6 GHzPoint target assessment curve atL-band and X-band.
(a)
L-band,f = 1.5 GHz
(b)
X-band,f = 9.6 GHz
(c)
Theoretical imageAs shown in Table2, compared with theoretical indicators, the resolution errors at L-band are 0.0246 m bigger than X-band in range direction and 1.9193 m more in azimuth section. Otherwise, the expansion ratio, the peak sidelobe ratio, and the integral sidelobe ratio are the same situations that the effects in L-band are worse than X-band, and the indicators in azimuth direction are worse than range direction. As shown in Table 4, the pixel positions of the point target at L-band and X-band are different from the theoretical position. It can be seen that the lower frequency imaging is seriously affected by the propagation.Table 4
Position of the point target in range section.
Theoretical position
L-band
X-band
(2048.00, 512.00)
(2031.00, 512.00)
(2047.00, 512.00)According to the analysis of Section3, the path errors also affect the position of point target. Therefore, we measure the position of two band SAR images in pixel. As can be seen from Table 4, the range direction offset is 17 pixels and 1 pixel in L-band and X-band compared with theoretical position in the image. The image offset is proportional to the path error.The simulation results show that path error caused by ionospheric and tropospheric refraction affects both range section and azimuth section imaging indicators especially on azimuth section, such as resolution, expansion ratio, peak sidelobe ratio, and integral sidelobe ratio. And the lower the frequency is, the more serious the effects are. Table4 shows that there is offset error imaging in range section. With the distance between the satellite and target increasing, the imaging results will get worse. The accurate estimation of ionosphere and troposphere model will assist the autofocus compensation to improve spaceborne SAR image quality and reduce the dependence on external calibration measurement.


## 6. Conclusions

In this paper, we use the engineering practical ionospheric electron density profile model to get the electron concentration of different heights. After that, we can get the ionospheric refractive index of different heights. As for troposphere, we use Hopfield model to get the tropospheric refractive index. The actual distance is calculated by path tracing method. Compared to the propagation path in the vacuum, the real path must be bent. The path error can lead to phase error which is the key to the spaceborne SAR imaging. After researching related literatures, we find that people care a lot about the effect of troposphere and ionosphere on SAR imaging independently, but few materials considering both of them. So the study of the effects of two kinds of refractions on the spaceborne SAR imaging performance is essential. Other imaging indicators like image quality analyzing still needs to be further studied. In the future, we will analyze the correction of refractive effects on spaceborne SAR imaging.



---

*Source: 913056-2014-05-20.xml*