Файл: Ersoy O.K. Diffraction, Fourier optics, and imaging (Wiley, 2006)(ISBN 0471238163)(427s) PEo .pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 28.06.2024

Просмотров: 923

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

308

COMPUTERIZED IMAGING TECHNIQUES I: SYNTHETIC APERTURE RADAR

Figure 17.2. A SAR image generated in September, 1995 [Courtesy of Center for Remote Imaging, Sensing and Processing, National University of Singapore (CRISP)].

systems carried on the space shuttles are similar SAR imaging systems. Examples of airborne SAR systems carried out with airplanes are the E-3 AWACS (Airborne Warning and Control System), used in the Persian Gulf region to detect and track maritime and airborne targets, and the E-8C Joint STARS (Surveillance Target Attack Radar System), which was used during the Gulf War to detect and locate ground targets.

An example of a SAR image is shown in Figure 17.2. This image of South Greenland was acquired on February 16, 2006 by Envisat’s Medium Resolution Imaging Spectrometer (MERIS) [ESA].

17.3RANGE RESOLUTION

In SAR as well as other types of radar, range resolution is obtained by using a pulse of EM wave. The range resolution has to do with ambiguity of the received signal due to overlap of the received pulse from closely spaced objects.

In addition to nearby objects, there are noise problems, such as random fluctuations due to interfering EM signals, atmospheric effects, and thermal variations in electronic components. Hence, it is necessary to increase signal-to- noise ratio (SNR) as well as to achieve large range resolution.

The distance R to a single object reflecting the pulse is tc=2, where t is the interval of time between sending and receiving the pulse, and c is the speed of light,

CHOICE OF PULSE WAVEFORM

309

3 108 m/sec. Suppose the pulse duration is T seconds. Then, the delay between two objects must be at least T seconds so that there is no overlap between the two pulse echoes. This means the objects must be separated by cT=2 meters (if MKS units are used). Reducing T results in better range resolution. However, high pulse energy is also required for good detection, and short pulses mean lower energy in practice. In order to avoid this problem, matched filtering discussed in Section 17.5 is often used to convert a pulse of long duration to a pulse of short duration at the receiver. In this way, the received echoes are sharpened, and the overall system possesses the range resolution of a short pulse. The peak transmitter power is also greatly reduced for a constant average power. In such systems, matched filtering is used both for pulse compression as well as detection by SNR optimization. This is further discussed in Section 17.5.

17.4CHOICE OF PULSE WAVEFORM

The shape of a pulse is significant in order to differentiate nearby objects. Suppose that pðtÞ is the pulse signal of duration T, which is nonzero for 0tT. The returned pulse from one object can be written as

p1ðtÞ ¼ 1pðt 1Þ

ð17:4-1Þ

where 1 is the attenuation constant, and 1 is the time delay. The returned pulse from a second object can be written as

p2ðtÞ ¼ 2pðt 2Þ

ð17:4-2Þ

The shape of the pulse should be optimized such that p1ðtÞ is as dissimilar from

p2

ðtÞ for 1 ¼6 2 as possible.

 

 

The most often used measure of similarity between two waveforms p1ðtÞ and

p2

ðtÞ is the Euclidian distance given by

 

 

D2 ¼ ð ½p1ðtÞ p2ðtÞ&2dt

ð17:4-3Þ

D2 can be written as

ð ð ð

D2 ¼ 21 p2ðt 1Þdt þ 22 p2ðt 2Þdt 2 1 2 pðt 1Þpðt 2Þdt

ð17:4-4Þ

The first two terms on the right-hand side above are proportional to the pulse energy, which can be separately controlled by scaling. Hence, only the last term is


310

COMPUTERIZED IMAGING TECHNIQUES I: SYNTHETIC APERTURE RADAR

significant for optimization. It should be minimized for

1 ¼6

2

in order to maximize

D2. The integral in the last term is rewritten as

 

 

Rð 1; 2Þ ¼ ð pðt 1Þpðt 2Þdt

 

ð17:4-5Þ

which is the same as

 

 

 

 

Rð Þ ¼ ð pðtÞpðt þ Þdt

 

 

ð17:4-6Þ

where equals 1 2 or 2 1. It is observed that Rð Þ is the autocorrelation of pðtÞ.

A linear frequency modulated (linear fm) signal, also called a chirp signal has the property of very sharp autocorrelation which is close to zero for 0. It can be written as

xðtÞ ¼ A cosð2pðft þ gt2ÞÞ

ð17:4-7Þ

or more generally as

xðtÞ ¼ ej2pðftþgt2Þ

ð17:4-8Þ

The larger g signifies larger variation of instantaneous frequency fi, which is the derivative of the phase:

fi ¼ f þ 2gt

ð17:4-9Þ

It is observed that fi varies linearly with t. The autocorrelation function of xðtÞ can be shown to be

ð

ð17:4-10Þ

Rð Þ ¼ ej2p ðf þg Þ ej4pg tdt

Suppose that a pulse centered at T0 and of duration T is expressed as

p t

Þ ¼

rect

t T0

 

x t

Þ

ð

17:4-11

Þ

ð

 

T

ð

 

The autocorrelation function of pðtÞ is given by

R X ¼ ej2p Ð ½f þ2gðT0þT2Þ&tri ðT j jÞ sin c pg2 ðT j jÞ


THE MATCHED FILTER

311

The width of the main lobe of this function is approximately equal to 1=gT. The rectangular windowing function is often replaced by a more smooth function such as a Gaussian function, as discussed in Section 14.2. Because of the properties discussed above, a chirp signal within a finite duration window is often the pulse waveform chosen for good range resolution.

17.5THE MATCHED FILTER

One basic problem addressed by matched filtering is to decide whether a signal of a given form exists in the presence of noise. In the classical case, the filter is also constrained to be a LTI system. Matched filters are used in many other applications as well, such as pulse compression as discussed in the next section, and image reconstruction as discussed in Sections 17.8–17.10.

Suppose the input consists of the deterministic signal xðtÞ plus noise NðtÞ. Let the corresponding outputs from the linear system used be x0ðtÞ and N0ðtÞ, respectively. This is shown in Figure 17.3. The criterion of optimality to be used to determine whether xðtÞ is present is the maximum SNR at the system output. It is shown below that the LTI filter that maximizes the SNR is the matched filter.

If n0ðtÞ is assumed to be a sample function of a wide-sense stationary (WSS) process, the SNR at time T0 can be defined as

SNR

¼

jx0ðT0Þj2

ð

17:5-1

Þ

E½N02ðT0Þ&

 

 

 

The output signal x0ðtÞ is given by

 

 

 

 

 

 

1

 

 

 

 

 

x0ðTÞ ¼

ð

Xðf ÞHðf Þej2pfT0 df

ð17:5-2Þ

 

1

 

 

 

 

The output average noise power is given by

 

1

 

 

E½N02ðtÞ& ¼

ð

jHðf Þj2SN ðf Þdf

ð17:5-3Þ

 

1

 

 

where SN ðf Þ is the spectral density of NðtÞ.

Figure 17.3. The matched filter as a LTI filter to optimally reduce noise.


312

COMPUTERIZED IMAGING TECHNIQUES I: SYNTHETIC APERTURE RADAR

 

Now, the SNR can be written as

 

 

 

 

 

 

2

 

 

 

 

 

 

1

X

ð

f

Þ

H f ej2pfT0

 

 

 

 

 

 

ð

 

 

ð Þ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

¼

 

1

 

 

 

 

2

 

ð

Þ

 

SNR

 

 

 

 

 

 

 

 

17:5-4

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

ð

jHðf Þj SN ðf Þdf

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

To optimize the SNR, the Schwarz inequality can be used. If A(f) and B(f) are two possibly complex functions of f, the Schwarz inequality is given by

 

 

1

Aðf ÞBðf Þdf

2

 

1 jAðf Þj2df

1 jBðf Þj2df

 

ð17:5-5Þ

 

 

ð

 

 

 

 

 

ð

 

 

 

 

ð

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

with equality iff

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Aðf Þ ¼ CB ðf Þ

 

 

 

 

 

 

 

 

ð17:5-6Þ

C being an arbitrary real constant. Let

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Aðf Þ ¼ SN ðf ÞHðf Þ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 fT

 

 

 

 

 

 

 

 

 

 

 

17:5-7

 

 

 

 

 

 

 

 

 

 

pj p

0

 

 

 

 

 

 

 

 

 

ð

Þ

 

 

 

 

 

B f

Þ ¼

Xðf Þe

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ð

 

SN ðf Þ

 

 

 

 

 

 

 

 

 

 

 

 

 

Then, the Schwarz inequality gives

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

2

1

SN ðf ÞjHðf Þj2df 32

1

 

X

f

2

 

3

 

 

 

Xðf ÞHðf Þej2pfT df

ð

ð

 

j ð

Þj

 

df

ð17:5-8Þ

S

f

 

ð

 

 

 

 

 

 

 

 

 

 

 

 

N

ð Þ

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

54

 

 

 

 

 

5

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

or

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

SNR

 

1

jXðf Þj2

df

 

 

 

 

 

 

 

 

 

17:5-9

 

 

 

 

 

 

 

ð

 

 

 

 

 

 

 

 

ð

Þ

 

 

 

 

 

 

 

 

SN ðf Þ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

The SNR is maximized and becomes equal to the right-hand side of Eq. (17.5-9) when the equality holds according to Eq. (17.5-6), in other words, when

H f

Þ ¼

H

f

Þ ¼

C

X ðf Þ

e j2pfT0

ð

17:5-10

Þ

ð

 

optð

S

f

Þ

 

 

 

 

 

 

 

 

 

N ð

 

 

 

 

The filter whose transfer function is given by Eq. (17.5-10) is called the matched filter. It is observed that Hoptðf Þ is proportional to the complex conjugate of the FT of