96

In this section the ISAR signal model for a monostatic configuration will be presented.

In this section the ISAR signal model for a monostatic configuration will be presented.

Figure 3.1 ISAR geometry

Consider the typical ISAR geometry shown in Figure 3.1 in which the transmitter and the receiver are colocated in the origin of the Cartesian reference system *T _{ξ}*. As stated in Chapter 1, the signal in the frequency/slow time domain at the output of the matched filter can be expressed as

3.1
$$S\left(f,n\right)=C\cdot \text{rect}\left(\frac{n}{N}\right){\left|{S}_{T}\left(f,n\right)\right|}^{2}{\int}_{V}{f}^{\prime}\left(\text{x}\right){e}^{-j2\pi f\tau \left(x,n\right)}d\text{x}$$

where in conventional active radar ${\left|{S}_{T}\left(f\right)\right|}^{2}\simeq \text{rect}\left(\frac{f-{f}_{0}}{B}\right)$
and $\text{rect}\left(\frac{n}{N}\right)=u\left(n+\frac{N}{2}\right)-u\left(n-\frac{N}{2}-1\right)$
defines the signal support in the frequency/slow time domain. The function
3.2
$$\tau \left(\text{x,}n\right)=\frac{2R\left(\text{x,}n\right)}{c}$$

where By assuming that the radar target distance is much greater than the size of the target (*straight iso-range* approximation), the term *R*(x, *n*) can be approximated as follows:

Figure 3.2 Rotation along the LoS

3.3
$$R\left(\text{x,}n\right)\simeq {R}_{0}\left(n\right)+{\text{x}}^{T}\cdot {\text{i}}_{LoS}\left(n\right)$$

where It is worth pointing out that Equation (3.3) holds true only when considering the target as a rigid body. In this case, its motion with respect to the radar can be considered as the superimposition of the translational motion of a reference point and a rotation vector applied to this reference point. The rotation vector takes into account both the rotation due to the translation and the proper target rotation. It is worth pointing out that, since the radar can measure only the distance variation along the LoS, any rotation around the *y*_{2} axis does not produce any effect. Then, the only component to be considered is **Ω**(*n*), as it contributes to the synthetic aperture formation. This result is evident by observing Figure 3.2 in which the case in which the target undergoes a rotation around the LoS is represented. In such a case, the distance between the radar and each scatterer on the target does not change, *R*(x, *t*_{1}) = *R*(x, *t*_{2}). Then, such a rotation does not produce any effect from the radar point of view.

Let us assume that the effective rotation vector is constant during the overall observation time, i.e., **Ω**(*n*) ≈ **Ω** for |*nT _{R}*| <

3.4
$$\dot{y}\left(n\right)={\Omega}_{T}\times y\left(n\right)$$

with the initial condition y(0) = x.
The resulting closed form solution is

3.5
$$y\left(n\right)=a+b\mathrm{cos}\left({\Omega}_{T}n{T}_{R}\right)+\frac{c}{{\Omega}_{T}}\mathrm{sin}\left({\Omega}_{T}n{T}_{R}\right)$$

where
$$\begin{array}{ll}\text{a}\hfill & =\frac{{\Omega}_{\text{T}}\cdot \text{x}}{{\Omega}_{T}^{2}}{\mathbf{\Omega}}_{T}\hfill \\ \text{b}\hfill & =\text{x}-\frac{{\Omega}_{\text{T}}\cdot \text{x}}{{\Omega}_{T}^{2}}{\mathbf{\Omega}}_{T}\hfill \\ \text{c}\hfill & ={\mathbf{\Omega}}_{T}\times \text{xand}{\Omega}_{T}=\left|{\mathbf{\Omega}}_{T}\right|\hfill \\ {\Omega}_{T}\hfill & =\left|{\mathbf{\Omega}}_{T}\right|\hfill \end{array}$$

with the above-mentioned choice of the reference system,
Equation (3.5) can be reasonably approximated by its first order Taylor series around *t* = 0 as follows:

3.6
$$y\left(n\right)\approx a+b+c\left(n{T}_{R}\right)=\text{x}+\text{c}\left(n{T}_{R}\right)$$

resulting in:
3.7
$${y}_{2}\left(\text{x},n\right)={x}_{2}+\Omega n{T}_{R}{x}_{1}$$

where Ω = |By substituting Equation (3.6) in Equation (3.1) and by defining $W\left(f,n\right)=C\cdot \text{rect}\left(\frac{f-{f}_{0}}{B}\right)$

3.8
$$S\left(f,n\right)=W\left(f,n\right){\displaystyle \underset{V}{\int}{f}^{\prime}\left(\text{x}\right)}{e}^{-j4\pi \frac{f}{c}\left({R}_{0}\left(n\right)+{x}_{2}+\Omega {x}_{1}n{T}_{R}\right)}d\text{x}$$

The same result can be obtained in a different way considering that the LoS unit vector in the *n ^{th}* sweep in the

3.9
$${\text{i}}_{LoS}\left(n\right)=\left[\begin{array}{c}\mathrm{sin}\left(\Omega n{T}_{R}\right)\\ \mathrm{cos}\left(\Omega n{T}_{R}\right)\\ 0\end{array}\right]$$

As a consequence, the inner product in Equation (3.3) results in

3.10
$${y}_{2}\left({\text{x}}^{T},n\right)=\text{x}\cdot {\text{i}}_{LoS}\left(n\right)={x}_{2}\mathrm{cos}\left(\Omega n{T}_{R}\right)+{x}_{1}\mathrm{sin}\left(\Omega n{T}_{R}\right)$$

Figure 3.3 ISAR processing chain

Equation (3.1) can therefore be written as

3.11
$$S\left(f,n\right)=W\left(f,n\right){\displaystyle \underset{{x}_{1},{x}_{2}}{\iint}f\left({x}_{1},{x}_{2}\right){e}^{-j\frac{4\pi f}{c}\left[{R}_{0}\left(n\right)+{x}_{2}\mathrm{cos}\left(\Omega n{T}_{R}\right)+{x}_{1}\mathrm{sin}\left(\Omega n{T}_{R}\right)\right]}d{x}_{1}d{x}_{2}}$$

where $f\left({x}_{1},{x}_{2}\right)={\displaystyle \int {x}_{3}}{f}^{\prime}\left(\text{x}\right)d{x}_{3}$
is the projection of the target's reflectivity function onto the image projection plane (IPP).
By assuming the total aspect angle variation is small during the observation time, i.e., Δ*ϑ* = Ω*T _{obs}* ≪ 1, Equation (3.11) becomes Equation (3.8).

The steps to be performed to obtain an ISAR image are summarized in the functional block diagram shown in Figure 3.3. A time interval of the overall observation within which the effective rotation vector can be assumed constant time is selected first by means of the *time windowing selection*. This operation is important because when the radar observes the target for long time intervals, this assumption does not hold true and the model in Equation (3.8) cannot be applied. Through the time windowing process, a subset of data for which the effective rotation vector is constant is properly selected to form the ISAR image. After that, the motion compensation operation is performed. This operation aims at removing the phase term *R*_{0}(*n*) that accounts for the target radial motion. After motion compensation, the image formation step is performed via a range-Doppler (RD) approach. This operation provides the ISAR image in the delay time/Doppler coordinates (*τ, ν*). Finally, in order to better asses the target size, a scaling operation to spatial coordinates is needed.

All these processing steps will be extensively described in the following sections.

The ISAR image is formed via the RD approach. Consider the received signal expressed in Equation (3.11). Assuming that the phase term due to the translational motion *R*_{0}(*n*) is perfectly compensated during the ISAR image formation processing, the motion compensated signal is expressed as follows:

Figure 3.4 Spatial frequency domain

3.12
$${S}_{C}\left({X}_{1},{X}_{2}\right)=W\left({X}_{1},{X}_{2}\right){\displaystyle \underset{{x}_{1},{x}_{2}}{\iint}f\left({x}_{1},{x}_{2}\right){e}^{-j2\pi \left({x}_{1}{X}_{1}+{x}_{2}{X}_{2}\right)}d{x}_{1}d{x}_{2}}$$

where (
3.13
$$\{\begin{array}{c}{X}_{1}\left(f,n\right)=\frac{2f}{c}\mathrm{sin}\left(\Omega n{T}_{R}\right)\\ {X}_{2}\left(f,n\right)=\frac{2f}{c}\mathrm{sin}\left(\Omega n{T}_{R}\right)\end{array}$$

From Equation (3.13) it is obvious that the motion compensated signal corresponds to the windowed Fourier transform (FT) of the projected target's reflectivity function. The term *W*(*X*_{1}, *X*_{2}) defines the region (a sector of annulus) in the spatial frequencies domain in which the signal is defined, see Figure 3.4.

As stated above, by assuming that the total aspect angle variation is small, i.e., Δ*ϑ* = Ω*T _{obs}* ≪ 1, the spatial frequencies can be approximated as follows:

3.14
$$\{\begin{array}{c}{X}_{1}\left(n\right)=\frac{2{f}_{0}}{c}\left(\Omega n{T}_{R}\right)\\ {X}_{2}\left(f\right)=\frac{2f}{c}\end{array}$$

It is worth noting that for *X*_{1} the frequency *f* has been substituted by the central frequency *f*_{0}, being the result of the approximation of the polar domain in Figure 3.4 with a rectangular window that intersects the angular sector at the coordinates ${X}_{2}=\frac{2{f}_{0}}{c}$

Figure 3.5 Spatial frequency domain with small aspect angle variation

By assuming that the target is composed by *K* point-like scatterers and neglecting all the interactions among the scatterers, the target reflectivity function can be expressed as follows:

3.15
$$f\left({x}_{1},{x}_{2}\right)={\displaystyle \sum _{k=1}^{K}{\sigma}_{k}\delta \left({x}_{1}-{x}_{{1}_{k}},{x}_{2}-{x}_{{2}_{k}}\right)}$$

where ${x}_{{1}_{k}}$
, ${x}_{{2}_{k}}$
and By substituting Equation (3.15) in Equation (3.12) the received signal after motion compensation results

3.16
$${S}_{C}\left({X}_{1},{X}_{2}\right)=W\left({X}_{1},{X}_{2}\right){\displaystyle \sum _{k=1}^{K}{\sigma}_{k}{e}^{-j2\pi \left({X}_{1}{x}_{{1}_{k}}+{X}_{2}{x}_{{2}_{k}}\right)}}$$

The ISAR image can be obtained by applying a 2D-FT to Equation (3.16) and can be expressed as

3.17
$$I\left({x}_{1},{x}_{2}\right)=w\left({x}_{1},{x}_{2}\right)\otimes \otimes f\left({x}_{1},{x}_{2}\right)={\displaystyle \sum _{k=1}^{K}{\sigma}_{k}w\left({x}_{1}-{x}_{{1}_{k}},{x}_{2}-{x}_{{2}_{k}}\right)}$$

By substituting Equation (3.14) in Equation (3.16), the received signal after motion compensation becomes:

3.18
$${S}_{C}\left(f,n\right)=W\left(f,n\right){\displaystyle \sum _{k=1}^{K}{\sigma}_{k}{e}^{-j2\pi \left(f{\tau}_{k}+n{T}_{R}{v}_{k}\right)}}$$

where the delay-time and the Doppler coordinates are defined as:
3.19
$$\begin{array}{c}\tau =\frac{2{x}_{2}}{c}\\ v=\frac{2{f}_{0}\Omega {x}_{1}}{c}\end{array}$$

The ISAR image in the delay time/Doppler domain, *I*(*τ, ν*) can then be obtained via a two-dimensional Fourier transform of Equation (3.18) as follows:

3.20
$$I\left(\tau ,v\right)=w\left(\tau ,v\right)\otimes \otimes {\displaystyle \sum _{k=1}^{K}{\sigma}_{k}\delta \left(\tau -{\tau}_{k},v-{v}_{k}\right)={\displaystyle \sum _{k=1}^{K}{\sigma}_{k}w}\left(\tau -{\tau}_{k},v-{v}_{k}\right)}$$

where
3.21
$$w\left(\tau ,v\right)=2\text{D-FT}\left[W\left(f,t\right)\right]=B{T}_{obs}\mathrm{sin}c\left(\tau B\right)\mathrm{sin}c\left(v{T}_{obs}\right){e}^{-j2\pi {f}_{0}\tau}$$

and 2D-FT[
3.22
$$\begin{array}{c}{\delta}_{\tau}=\frac{1}{B}\\ {\delta}_{v}=\frac{1}{{T}_{obs}}\end{array}$$

It is worth pointing out that Equation (3.20) represents the ISAR image in the delay time/Doppler domain. To obtain the target ISAR image in the spatial coordinates (range/cross-range) a scaling operation is needed. This scaling operation can be performed by exploiting Equation (3.19) as follows:

3.23
$$\{\begin{array}{l}{x}_{1}=\frac{c}{2{f}_{0}\Omega}v\hfill \\ {x}_{2}=\frac{c}{2}\tau \hfill \end{array}$$

and, as a consequence, the spatial resolutions can be expressed as:
3.24
$$\{\begin{array}{l}{\delta}_{{x}_{1}}=\frac{c}{2{f}_{0}\Omega {T}_{obs}}\hfill \\ {\delta}_{{x}_{2}}=\frac{c}{2B}\hfill \end{array}$$

Equation (3.23) shows that the scaling operation from the delay time to the range coordinate is straightforward while this is not true for the scaling operation from the Doppler frequency to the cross-range coordinates. In fact, this coordinates conversion involves the knowledge of the modulus of the effective rotation vector, Ω, which depends on the motion of the non-cooperative moving target and, as a consequence, cannot be assumed known in the ISAR scenario.

As can be noted from Equation (3.24), the resolution becomes finer as the aspect angle variation increases. This is in contrast with the small aspect angle variation assumption which is mandatory for the applicability of the range Doppler image formation algorithm. This hypothesis determines an upper bound in the cross-range resolution that can be achieved with this algorithm.

Another important consideration concerns the definition of the IPP where the target is represented. As stated above, the reference system *T _{x}* has been chosen so that the

3.25
$$\Omega \left(n\right)={\text{i}}_{LoS}\times \left[{\Omega}_{T}\left(n\right)\times {\text{i}}_{Lo\text{S}}\right]$$

Therefore, the cross-range axis x_{1} is defined as:

3.26
$${\text{x}}_{1}=\Omega \times {\text{i}}_{LoS}$$

where the dependency on It is worth pointing out that the above-mentioned RD technique is based on the assumption that the Doppler frequency of each scatterer, relative to the focusing point, is constant during the CPI (coherent processing interval). This hypothesis usually holds true for low spatial resolution (of the order of meter) and in the case in which the target does not undergo fast maneuvers and/or is affected by significant oscillation such as roll, pitch and yaw. In the case in which very high resolutions are required, longer CPIs are needed and the Doppler frequency becomes time-varying. This can happen also when the target undergoes angular motions, as in the case of ships. In this case the first order Taylor approximation in Equation (3.6) and (3.7) is not valid. To solve this problem, the *time frequency transforms* (TFTs) is used instead of the RD approach to form ISAR images. TFTs are, in fact, suitable for the analysis of non-stationary signals. An example of a TFT-based technique is the range instantaneous Doppler (RID) approach described in [2].

As stated before, the applicability of the RD technique relies on the assumption that the effective rotation vector is constant during the observation time, i.e., **Ω**(*n*) *≈* **Ω**. However, the target's own motion may induce a non-uniform target rotation vector. In particular, maritime targets can undergo very complex motions (roll, pitch and yaw) induced by the sea surface. In order to minimize the target's rotation variations, the CPI can be controlled via a time-windowing approach.

Typically, a fixed window length can be set and a sliding window processing along the whole observation time is applied to obtain a sequence of ISAR images of the targets. Among the obtained frames, the more suitable for target classification and recognition is selected.

It is worth pointing out that, in order to obtain good classification and recognition performances, the most important requirement is a good level of focus, as the target details would then appear sharper than in defocused images. As stated above, the image cross-range resolution is inversely proportional to the total aspect angle variation and, as a consequence, to the considered CPI so that longer CPI involves a finer resolution in the cross-range dimension. On the other hand, longer CPIs increase the chance that the target rotation vector may not be considered constant and therefore producing a defocused image. It is quite obvious that a trade-off must be identified to obtain a well focused and a high resolved image at the same time.

An automatic time windowing technique is proposed in [9] and will be here recalled for the sake of clarity. Specifically, the position of the window along the whole observation time and its length are automatically selected in order to obtain the image with the highest focus. The concept can be explained referring to Figure 3.6, in which the data is refereed to be distributed along the time axis *n* (*n* denotes the sweep index at time *t* = *nT _{R}*). The considered temporal window is defined by two parameters, namely

The measure of the image focus is made through the image contrast (IC) which is assumed to be maximum when the window position (*n _{w}*) and the window length (

In particular, the optimal time window is obtained by maximizing the IC with respect to the pair (*n _{w}*,

Figure 3.6 Time windowing concept

Figure 3.7 Optimal time window estimator

3.27
$$\left({n}_{w,opt},{N}_{w,opt}\right)=\underset{\left({n}_{w},{N}_{w}\right)}{\mathrm{arg}\mathrm{max}}\left(IC\left({n}_{w},{N}_{w}\right)\right)$$

It should be noted that the variables (*n _{w}*,

- Contrast maximization with respect to
*n*for a given guess_{w}*N*. Let_{w,in}*n*be the solution of such maximization;_{w,opt} - Optimization with respect to
*N*with_{w}*n*=_{w}*n*._{w,opt}

The procedure is depicted in Figure 3.7 for the sake of clarity. The justification of this procedure is heuristic. It can be observed both in simulated and real ISAR data of several types of targets that the position of the optimal time-window is almost independent on the window length and, as a consequence, the IC peak position along the slow-time axis, *n*, does not change by changing *N _{w}*. This can be physically explained by considering that the target own motion is characterized by regular motions at given times and less regular motion at other times. For example, a ship usually undergoes regular pitch and roll motions due to the sea surface state. This regular motion may be disturbed by an incoming wave that generates complex motions which cause the rotation vector to change rapidly. For this reason, the estimation of

As stated in Section 3.2.1, the radial motion compensation is a fundamental step in forming well-focused ISAR images. This step is typically termed *image autofocus* and consists of removing the phase term $-\frac{4\pi f}{c}{R}_{0}\left(n\right)$

Several techniques have been developed by the ISAR research community, each of them showing pros and cons. Such techniques can be classified as parametric and non-parametric methods [1]. Parametric methods need a parametric model of the radar received signal while non-parametric techniques do not make any assumption. Some autofocus techniques are detailed in the next subsections.

The *image contrast based autofocus* (ICBA) is a parametric technique that is based on the obvious concept that a higher value of IC corresponds to a more focused image [10]. This autofocus technique aims at removing the term *R*_{0}(*n*) via an iterative estimation of the motion parameters. For a relatively small observation time, *T _{obs}* =

3.28
$${R}_{0}\left(n\right)={\displaystyle \sum _{l=0}^{L}\frac{1}{l!}}{\alpha}_{l}{\left(n{T}_{R}\right)}^{l}$$

in which The estimation of *R*_{0}(*n*) results in the estimation of the target motion parameters via the maximization of the IC with respect to *α* as:

3.29
$$\widehat{\alpha}=\underset{\alpha}{\mathrm{arg}\mathrm{max}}\left[IC\left(\alpha \right)\right]$$

where the IC is defined as follows:
3.30
$$IC\left(\alpha \right)=\frac{\sqrt{{A}_{v}\left[{\left[I\left(\tau ,v,\alpha \right)-{A}_{v}\left[I\left(\tau ,v,\alpha \right)\right]\right]}^{2}\right]}}{{A}_{v}\left[I\left(\tau ,v,\alpha \right)\right]}$$

and where
3.31
$$I\left(\tau ,v,\alpha \right)={\left|2\text{D-FT}\left[S\left(f,n\right){e}^{j\frac{4\pi f}{c}{R}_{0}\left(\alpha ,n\right)}\right]\right|}^{p}$$

Referring to Equation (3.31), *p* = 1 in the case of image amplitude and *p* = 2 in the case of image intensity. *R*_{0}(*α, n*) is the radial motion model evaluated with the *α* motion parameters.

A similar approach to the one proposed in Section 3.2.3.1 can be devised by substituting the IC with the *image entropy* (IE) expressed as follows:

3.32
$$IE\left(\alpha \right)=-\int \int \mathrm{ln}\left(\overline{I}\left(\tau ,v,\alpha \right)\right)\overline{I}\left(\tau ,v,\alpha \right)d\tau dv$$

where $\overline{I}\left(\tau ,\nu ,\alpha \right)=\frac{I\left(\tau ,\nu ,\alpha \right)}{{A}_{\upsilon}\left[I\left(\tau ,\nu ,\alpha \right)\right]}$
.
As for the IC, the IE is a good indicator of the image focus. Unlike the IC, the smaller IE value the better the image focus.

The *dominant scatterer autofocus* (DSA), also known as the *hot spot* (HS), is a non-parametric two stages technique. The first stage consists of a rough alignment before applying the phase compensation as a second step. This technique is derived from the studies in two different areas of research, namely the time delay estimation [3] and the adaptive beamforming [12]. A complete description of this algorithm can be found in [4] while a brief overview is reported here.

After measuring and storing the complex envelopes of the echo samples, high resolution range profiles in the *n ^{th}* sweep,

3.33
$${S}_{\alpha}\left(\tau ,n\right)=A\left(\tau ,n\right){e}^{j\varphi \left(\tau ,n\right)}$$

where Then, a search along the range coordinate is performed in order to find a dominant and stable scatterer. The range cell where such scatterer is located is the reference range denoted by *τ*_{0} and the corresponding range profile is given by

3.34
$${S}_{\alpha}\left({\tau}_{0},n\right)={S}_{0}\left(n\right)=A\left({\tau}_{0},n\right){e}^{j{\varphi}_{0}\left(n\right)}$$

The *τ*_{0} value is found by measuring the normalized echo variance in each range cell and is determined to be the range for which the variance value is minimum. This approach relies on the assumption that a dominant scatterer with a large radar cross-section exists and, therefore, the measured phase can be attributed to the phase generated by one point scatterer. The next step performs a phase compensation by using the phase history of *S*_{0}(*n*). For the range cell data corresponding to *τ*_{0} the result is expressed as follows:

3.35
$${S}_{0,c}\left(n\right)=A\left({\tau}_{0},n\right)\approx A$$

Applying the same operation to the other range cells, motion compensation is achieved and the result is expressed as

3.36
$${S}_{C}\left(\tau ,n\right)=A\left(\tau ,n\right){e}^{j\left[\varphi \left(\tau ,n\right)-{\varphi}_{0}\left(n\right)\right]}$$

This algorithm is also known as *minimum variance algorithm* (MVA) due to the criterion used to choose the reference range cell.

A more robust version of this algorithm, which combines the echoes from several range cells has been proposed in [4] and it is named *multiple scatterer algorithm* (MSA). It is based on an average operation performed on the phase differences of *M _{c}* reference scatterer (after phase unwrapping) in order to provide the phase correction. Typically, three reference range bins (

3.37
$${S}_{m}\left(n\right)=A\left({\tau}_{m},n\right){e}^{j{\varphi}_{m}\left(n\right)}$$

and let
3.38
$${\varphi}_{0}\left(n\right)=\frac{1}{{M}_{c}}{\displaystyle \sum _{m=1}^{{M}_{c}}{\varphi}_{m}\left(n\right)}$$

Phase compensation is then performed as expressed in Equation (3.35) and (3.36). The DSA algorithm is summarized in the functional block depicted in Figure 3.8.

The *Phase Gradient Autofocus* (PGA) can be seen as an extension of the DSA algorithm proposed in [4]. It also exploits the information that remains in those range cells discarded by the *minimum variance* (MV) range selection. The PGA algorithm replaces the phase compensation approach used in DSA with a solution based on the *maximum likelihood* (ML) estimator. The ML estimator theoretically uses the information contained in all the range cells. However, only the range cells in which the SNR is high enough are considered to improve the PGA performance.

Figure 3.8 DSA algorithm functional block diagram

Figure 3.9 PGA functional block diagram

The functional block diagram that summarizes the PGA algorithm is depicted in Figure 3.9.

A range Doppler image is formed and a rough range alignment is performed. Then, the peak value in each range cell, which is supposed to be the return from the dominant scatterer, is found along the Doppler dimension and centershifted and filtered in the Doppler domain (low-pass filter). Each range cell is then transformed back in the pulse domain via an inverse Fourier transform (IFT) to obtain the phase shifted and filtered time histories.

This operation corresponds to isolating the strongest scatterer contribution for each range cell and forcing it to zero-Doppler, which can be interpreted as a way to remove the scatterer radial motion.

Then the phase gradient is estimated by using the ML approach as described as follows.

Let the arbitrary *m ^{th}* range cell time history be represented as follows, where two consecutive samples are considered:

3.39
$$\begin{array}{c}{S}_{\alpha}\left({\tau}_{m},n-1\right)=A\left({\tau}_{m}\right){e}^{j\varphi \left(n-1\right)}+{S}_{n}\left({\tau}_{m},n-1\right)\\ {S}_{\alpha}\left({\tau}_{m},n\right)=A\left({\tau}_{m}\right){e}^{j\varphi \left(n\right)}+{S}_{n}\left({\tau}_{m},n\right)\end{array}$$

where τThe phase gradient between two consecutive samples can be defined as

3.40
$$\Delta \varphi =\varphi \left(n\right)=\varphi \left(n-1\right)$$

The derivation of the ML estimator of the phase gradient is given in [5]. The solution is expressed as follows

3.41
$$\Delta \widehat{\varphi}\left(n\right)=\angle \left({\displaystyle \sum _{m=1}^{{M}_{c}}{S}_{\alpha}\left({\tau}_{m},n\right){S}_{\alpha}^{*}\left({\tau}_{m},n-1\right)}\right)$$

where The phase correction term is then obtained by integrating the estimated phase gradient, as follows:

3.42
$$\begin{array}{c}\widehat{\varphi}\left(n\right)={\displaystyle {\sum}_{l=1}^{n}\Delta \widehat{\varphi}\left(l\right)}\\ \widehat{\varphi}\left(0\right)=0\end{array}$$

Finally, the phase correction term is used to compensate the phase error and obtain the compensated signal.

The ISAR processing generates a two-dimensional high-resolution image of a target in the delay-time/Doppler domain. In order to determine the size of the target it is, however, required to have a fully scaled image in the range/cross-range domain. The scaling along the range dimension can be performed by using the well-known relationship ${x}_{2}=\frac{c\tau}{2}$

, whereUnder the assumption that the target rotation vector is constant within the CPI, the chirp rate (which corresponds to the variation rate of the Doppler frequency of a scatterer) produced by the scattering centers is related to the modulus of the target rotation vector. Therefore, each scattering center carries information about the modulus of the target rotation vector through its chirp rate. To demonstrate this we will consider the radial motion compensated signal expressed in Equation (3.11).

The phase term in the integral can be approximated with a second order Taylor polynomial expansion, as follows:

3.43
$$\varphi \left(f,n;{x}_{1},{x}_{2}\right)=-j\frac{4\pi f}{c}\left({x}_{2}+{x}_{1}\Omega n{T}_{R}-\frac{1}{2}\cdot {x}_{2}{\Omega}^{2}{\left(n{T}_{R}\right)}^{2}\right)$$

so that the echo of an ideal scatterer located at $\left({x}_{{1}_{k}},{x}_{{2}_{k}}\right)$
, with finite reflectivity function $f\left({x}_{1},{x}_{2}\right)={\sigma}_{k}\delta \left({x}_{1}-{x}_{{1}_{k}},{x}_{2}-{x}_{{2}_{k}}\right)$
, can be written as follows:
3.44
$${S}_{C}\left(f,n;{x}_{{1}_{k}},{x}_{{2}_{k}}\right)=W\left(f,n\right){\sigma}_{k}{e}^{\left(-j\frac{4\pi f}{c}\left({x}_{{2}_{k}}+{x}_{{1}_{k}}\Omega n{T}_{R}-\frac{1}{2}\cdot {x}_{{2}_{k}}{\Omega}^{2}{\left(n{T}_{R}\right)}^{2}\right)\right)}$$

After applying an FT to the signal in Equation (3.44) along the frequency domain, the range compressed signal is obtained, as shown in Equation (3.45):

3.45
$$\begin{array}{c}{S}_{C}\left(\tau ,n;{x}_{{1}_{k}},{x}_{{2}_{k}}\right)=B{\sigma}_{k}\mathrm{sin}c\left(B\left(\tau -\frac{2}{c}{x}_{{2}_{k}}\right)\right)\text{rect}\left(\frac{n}{N}\right)\\ {e}^{-j\frac{4\pi {f}_{0}}{c}\left({x}_{{2}_{k}}+{x}_{{1}_{k}}\Omega n{T}_{R}-\frac{1}{2}{x}_{{2}_{k}}{\Omega}^{2}{\left(n{T}_{R}\right)}^{2}\right)}\end{array}$$

It should be noted that the phase term in Equation (3.45) represents the phase of a chirp signal, where ${\mu}_{k}=\frac{2{f}_{0}}{c}{x}_{{2}_{k}}{\Omega}^{2}$

is usually referred to as
3.46
$$\mu =\frac{2{f}_{0}}{c}{x}_{2}{\Omega}^{2}$$

In fact, as the scatterer range coordinate *x*_{2} can be readily obtained by measuring the delay-time *τ*, the effective rotation vector can be obtained by inverting Equation (3.46), as follows:

3.47
$$\Omega =\sqrt{\frac{c}{2{f}_{0}{x}_{2}}\mu}$$

In practice, a scatterer chirp rate, as well as its range, must be estimated from the received data. Therefore, the estimation of the effective rotation vector magnitude would, in general, be affected by an error. Different techniques for estimating target scattering center chirp rates making use of atomic decomposition [15], CLEAN technique [13, 8] or based on the IC method have been proposed in [7].

Specifically, the method proposed in [7] is considered here, as it has been proven to be the most effective. According to such a method, the signal components of the scattering centers are processed by means of a polynomial Fourier transform (PFT) in order to estimate the chirp rate. The use of such a method, namely local polynomial Fourier transform (LPFT), requires the solution of an optimization problem [14]. It has been largely proven in literature that the image contrast (IC) and the image entropy (IE) are good indexes to assess the quality of an image. Therefore, the proposed algorithm estimates the scatterers chirp rate that maximize the IC.

To make the estimation more accurate and robust, the chirp rates of a number of target scatterers can be measured together with their ranges. The problem of estimating the effective rotation vector magnitude is then transformed into a problem of estimating the slope of a straight line that fits the scatterplot generated by the set of range and chirp rate estimates. One way of solving this problem is to apply a least square error (LSE) approach [6]. The mathematical problem can be set as follows:

3.48
$${\widehat{\mu}}_{k}=a{x}_{{2}_{k}}+b+{\epsilon}_{k}$$

where $a=\frac{2{f}_{0}}{c}{\Omega}^{2}$
,
3.49
$$\left(\widehat{a},\widehat{b}\right)=\mathrm{arg}\underset{a,b}{\mathrm{min}}{\displaystyle \sum _{k=1}^{K}{\epsilon}_{k}^{2}}$$

3.50
$$\widehat{a}=\frac{K{\displaystyle {\sum}_{K=1}^{K}{\widehat{\mu}}_{k}{x}_{{2}_{k}}-{\displaystyle {\sum}_{k=1}^{K}{\widehat{\mu}}_{k}}{\displaystyle {\sum}_{k=1}^{K}{x}_{{2}_{k}}}}}{K{\displaystyle {\sum}_{k=1}^{K}{\left({x}_{{2}_{k}}\right)}^{2}-{\left({\displaystyle {\sum}_{k=1}^{K}{x}_{{2}_{k}}}\right)}^{2}}}$$

For the sake of clarity, the algorithm steps can be summarized as follows:

- Image segmentation
- Chirp rate estimation via LPFT
- Effective rotation vector estimation

In detail, after motion compensation, an ISAR image obtained by means of the RD technique can be modeled as the composition of several sinc-like functions centered at each scatterer location in the IPP, as shown in Equation (3.51).

3.51
$$I\left(\tau ,v\right)=C{\displaystyle \sum _{k=1}^{K}{\sigma}_{k}\mathrm{sin}c\left[B\left(\tau -\frac{2{x}_{{2}_{k}}}{c}\right)\right]\mathrm{sin}c\left[{T}_{obs}\left(v-\frac{2{f}_{0}\Omega {x}_{{1}_{k}}}{c}\right)\right]}$$

where
3.52
$${S}_{C}\left({\tau}_{k},n\right)={\sigma}_{k}{e}^{-j\frac{4\pi {f}_{0}}{c}{x}_{{2}_{x}}}{e}^{\left[-j\frac{4\pi {f}_{0}}{c}\left({x}_{{1}_{k}}\Omega n{T}_{R}-\frac{{x}_{{2}_{k}}}{2}{\Omega}^{2}{\left(n{T}_{R}\right)}^{2}\right)\right]}$$

By applying a second order LPFT [14] a complex three variables ISAR image can be obtained in the delay time/Doppler/chirp rate domain as follows:

3.53
$${S}_{C}\left({\tau}_{k},v,\mu \right)={\sigma}_{k}{e}^{-j\frac{4\pi {f}_{0}}{c}{x}_{{2}_{x}}}{\displaystyle \sum _{n}{e}^{\left[-j\pi \left(\left(v+\frac{2{f}_{0}}{c}{x}_{{1}_{k}}\Omega \right)n{T}_{R}+\frac{1}{2}\left(\mu -\frac{2{f}_{0}}{c}{x}_{{2}_{k}}{\Omega}^{2}\right){\left(n{T}_{R}\right)}^{2}\right)\right]}}$$

From Equation (3.53) it is possible to estimate the Doppler frequency, *ν _{k}*, and the chirp rate,

From the estimated value of *μ _{k}* and the scatterer location

In a realistic scenario, the signal to work with is fully digital so all the variables introduced in the signal modeling are discrete and can be expressed as:

3.54
$$\begin{array}{cc}f={f}_{0}+{m}_{f}{\Delta}_{f}& {m}_{f}=0,\cdots ,{M}_{f}-1\\ \tau =m{\Delta}_{\tau}& m=0,\dots ,M-1\\ v={n}_{v}\Delta v& {n}_{v}=0,\cdots ,{N}_{v}-1\end{array}$$

Since the relationship between (*f, n*) and (*τ, ν*) is given by a Fourier transformation, the sampling intervals in the image domain depend on the ISAR signal spectral occupancy, as follows:

3.55
$$\begin{array}{c}{\delta}_{\tau}=\frac{1}{M\Delta f}=\frac{1}{B}\\ {\delta}_{v}=\frac{1}{N{T}_{R}}=\frac{1}{{T}_{obs}}\end{array}$$

where
3.56
$$\begin{array}{c}{\delta}_{r}=\frac{c}{2B}\\ {\delta}_{cr}=\frac{c}{2{f}_{0}\Omega {T}_{obs}}\end{array}$$

Moreover, the above-mentioned sampling operation in the ISAR signal domain leads to an ambiguity in the ISAR image domain that must be taken into account in the radar design. Specifically, the non-ambiguity region in the cross-range and range dimension, i.e., ${\Delta}_{{x}_{1}}$

and ${\Delta}_{{x}_{2}}$ , can be expressed by considering the discrete spatial frequency domain shown in Figure 3.10, as follows:
3.57
$$\begin{array}{c}{\Delta}_{{x}_{1}}=\frac{c}{2{f}_{0}\Omega {T}_{R}}\\ {\Delta}_{{x}_{2}}=\frac{c}{2\Delta f}\end{array}$$

In order to avoid ambiguity, the relations in Equation (3.57) define the maximum sizes of the target. Any target bigger than the non-ambiguity region will be seen as folded within the image making target recognition more difficult. It is worth pointing out that ${\Delta}_{{x}_{1}}$

depends also on the modulus of the effective rotation vector, Ω, which, as stated above, cannot be a priori predicted leading to a more difficult design of an ISAR system.Figure 3.10 Non-ambiguity regions

Moreover, considering Equation (3.24) some remarks about spatial resolutions can be made. Specifically, the range resolution depends on the transmitted signal bandwidth and can be directly controlled in the radar design process, since it can be considered as a known parameter. On the other hand, it is evident that the spatial resolution in the cross-range dimension depends on the total aspect angle variation, Δ*ϑ* = *ΩT _{obs}*. This results in a parameter that is not controllable and cannot be selected during the radar design process, unless making some hypothesis about the target motion

The main concepts and theory behind ISAR imaging have been treated in this chapter. More specifically, the ISAR geometry and the received signal model have been defined. Then, mathematical details have been introduced in this chapter to support the main concepts behind ISAR technique and to provide the basis for implementing basic ISAR imaging algorithms. The concepts presented in this chapter mainly apply to monostatic active radar system. Nevertheless, the ISAR signal model formulation presented in this chapter is independent of the transmitted waveform; therefore, as it will be clarified in the following chapters, the same concepts apply to active and passive bistatic radar systems.

F. Berizzi, M. Martorella, B. Haywood, E. Dalle Mese, and S. Bruscoli. A survey on ISAR autofocusing techniques. Image Processing, 2004. ICIP ’04. 2004 International Conference on, 1:9–12, Oct. 2004.

F. Berizzi, E.D. Mese, M. Diani, and M. Martorella. High-resolution ISAR imaging of maneuvering targets by means of the range instantaneous doppler technique: modeling and performance analysis. Image Processing, IEEE Transactions on, 10(12):1880–1890, Dec 2001.

G. Clifford Carter. Time delay estimation for passive sonar signal processing. Acoustics, Speech and Signal Processing, IEEE Transactions on, 29(3):463–470, Jun 1981.

B. Haywood and R. J. Evans. Motion compensation for ISAR imaging. Proceedings of ASSPA 89, pages 113–117, 1989.

C.V.J. Jakowatz, D.E. Wahl, P.H. Eichel, D.C. Ghiglia, and P.A. Thompson. Spotlight-Mode Synthetic Aperture Radar: A Signal Processing Approach: A Signal Processing Approach. Springer, 2011.

Steven M. Kay. Fundamentals of Statistical Signal Processing: Estimation Theory. Signal Processing. 1993.

M. Martorella. Novel approach for ISAR image cross-range scaling. Aerospace and Electronic Systems, IEEE Transactions on, 44(1):281–294, 2008.

M. Martorella, N. Acito, and F. Berizzi. Statistical CLEAN technique for ISAR imaging. IEEE Tr. on Geoscience and Remote Sensing, 45(11):3552–3560, 2007.

M. Martorella and F. Berizzi. Time windowing for highly focused ISAR image reconstruction. Aerospace and Electronic Systems, IEEE Transactions on, 41(3):992–1007, 2005.

M. Martorella, F. Berizzi, and B. Haywood. Contrast maximisation based technique for 2-D ISAR autofocusing. Radar, Sonar and Navigation, IEE Proceedings, 152(4):253–262, 2005.

R. Gray Parker and Ronald L. Rardin. Discrete Optimization. Academic Press Professional, Inc., San Diego, CA, 1988.

B.D. Steinberg. Radar imaging from a distorted array: The radio camera algorithm and experiments. Antennas and Propagation, IEEE Transactions on, 29(5):740–748, Sept. 1981.

J. Tsao and B.D. Steinberg. Reduction of sidelobe and speckle artifacts in microwave imaging: The CLEAN technique. IEEE Tr. on Antennas and Propagation, 36:543–556, 1988.

Yongmei Wei and Guoan Bi. Efficient analysis of time-varying multicomponent signals with modified LPTFT. EURASIP J. Adv. Sig. Proc., 2005(8):1261–1268, 2005.

O.A. Yeste-Ojeda, J. Grajal, and G. Lopez-Risueno. Atomic decomposition for radar applications. Aerospace and Electronic Systems, IEEE Transactions on, 44(1):187–200, Jan. 2008.