## 1.

## Introduction

Ground moving target indication plays an important role in the synthetic aperture radar (SAR) image processing community. For a moving point scatterer, the image is displaced in an azimuth dimension and then defocused in the range dimension due to its range velocity component. In addition, it is smeared in the azimuth dimension due to its azimuth velocity component. In the range Doppler domain of the returns from a moving target, the range velocity component leads to a nonzero Doppler centroid and the azimuth velocity component leads to a biased Doppler modulation rate with respect to that of the background. Researchers have exploited the affected phase information to identify the moving targets in a complex-valued SAR image.^{1}2.3.4.^{–}^{5}

For a moving target with a nonzero range velocity component, it is possible for it to be detected by computing the Doppler centroid in the returns via pulse-by-pulse, and the corresponding algorithms, which are relatively mature.^{6}7.^{–}^{8} For a target moving in the azimuth direction only, many methods have been proposed, such as shear-average based autofocusing,^{2}^{,}^{3} information extraction from the antenna beam pattern,^{4} and image stack technology.^{5} Most of them are based on using phase estimation or Doppler filter-bank techniques. However, the traditional digital algorithms mainly concentrate on operations such as Fourier transformation and phase estimation. These algorithms are so complex that they are usually computationally expensive to implement. Furthermore, they are not robust enough to deal with clutters or interference.

As an optical cylindrical lens acts as a natural Fourier transform operator, it is used in many early SAR processors.^{9} In this paper, we will present an optical moving target indicator (MTI) to detect moving targets with nonzero azimuth velocities in a complex-valued SAR image, and show how to estimate their azimuth velocities.

Using our approach, two refocused SAR images are generated using two refocusing filters. In one refocused SAR image, a moving target will be in better focus, while in the other, it will be defocused. By comparing the sharpness of the two refocused SAR images, it is possible to identify the moving target.

Our optical MTI architecture uses two groups of refocusing subsystems. Each of the subsystems consists of two cylindrical lenses and one phase-only matching filter. The majority of this paper is a brief analytical description of the detection and azimuth velocity estimation methods that we used. The proposed architecture has been verified by both computer simulations and field data experiments. The results show that the optical processor is effective and practicable.

The paper is organized as follows. The principles for moving target detection and azimuth velocity estimation in the spatial domain of a given SAR image are reviewed. We then propose and discuss an optical processor using the detection and estimation principles. Both simulated and field data are used to verify the proposed scheme. Finally, we state our conclusion of the proposed concept.

## 2.

## Principles

## 2.1.

### Moving Target Detection in Space Domain

A strip-map mode SAR is considered in the following analysis. For a stationary point scatterer, its azimuth phase history has the form:^{10}

## (2)

$${s}_{a}(x)=\mathrm{exp}[-j\pi {K}_{a}{(1-{\epsilon}_{a})}^{2}{x}^{2}],\phantom{\rule[-0.0ex]{2em}{0.0ex}}|x-{x}_{0}|\le {l}_{s},$$Taking the Fourier transform of ${s}_{a}(x)$ according to the stationary phase theory and ignoring the amplitude constant and the phase constant, one obtains

## (3)

$${S}_{a}(u)=\mathrm{exp}\left[j\frac{{u}^{2}}{4\pi {(1-{\epsilon}_{a})}^{2}{K}_{a}}\right]\mathrm{exp}(-ju{x}_{0}),\phantom{\rule[-0.0ex]{2em}{0.0ex}}|u|\le \frac{2\pi}{{l}_{a}},$$In the traditional range Doppler algorithm,^{11} if ${S}_{a}(u)$ is filtered by

## (6)

$${S}_{f}(u)=\mathrm{exp}\left(j\frac{{\epsilon}_{a}{u}^{2}}{2\pi {K}_{a}}\right)\mathrm{exp}(-ju{x}_{0}),\phantom{\rule[-0.0ex]{2em}{0.0ex}}|u|\le \frac{2\pi}{{l}_{a}}.$$Taking the inverse Fourier transform of the above in Eq. (6), the focused signal can be shown to be

where## (9)

$${s}_{r}(x)={\mathsf{F}}_{u\to x}^{-1}\left[\mathrm{exp}\right(j\frac{{\epsilon}_{a}{u}^{2}}{2\pi {K}_{a}}\left)\right],\phantom{\rule[-0.0ex]{2em}{0.0ex}}|u|\le \frac{2\pi}{{l}_{a}},$$Equation (10) indicates that a residual quadratic phase term appears in the refocused signal. This residual term affects the focusing effectiveness.

The effect of the residual quadratic phase is shown in Fig. 1. In Fig. 1(a), the solid line represents the time frequency distribution of a stationary target, the dotted line and the dashed line represent that of two moving targets with the negative and positive azimuth velocity component, respectively.

It can be seen that the Doppler chirp rate of a moving target varies with its azimuth velocity component. In the traditional range Doppler imaging algorithm, the azimuth matched filter is optimally designed for the stationary background, so when it is used to compress the azimuth signal history of a moving target, a residual quadratic Doppler phase appears as shown in Fig. 1(b). In this figure, after being compressed in azimuth, it can be seen that the quadratic Doppler phase term of a stationary target is ideally removed after azimuth compression, which means that the stationary target can be focused optimally.

For the two moving targets, as the residual quadratic Doppler phase terms exist, the compressed signal of the target with ${\epsilon}_{a}>0$ and that of the target with ${\epsilon}_{a}<0$ have the pulse width of $\mathrm{\Delta}{t}_{1}$ and $\mathrm{\Delta}{t}_{2}$, respectively.

As the bandwidth of ${s}_{r}(x)$ is $4\pi /{l}_{a}$, its azimuth width in space domain satisfies

## (11)

$$2\pi \frac{{K}_{a}}{2|{\epsilon}_{a}|}\phantom{\rule[-0.0ex]{2em}{0.0ex}}\mathrm{\Delta}x=\frac{4\pi}{{l}_{a}},$$If the SAR image is refocused in azimuth with a pair of new refocusing filters described by

and where ${\epsilon}_{a}^{\prime}={V}_{a}^{\prime}/{V}_{p}$ with ${V}_{a}^{\prime}$ is a probing azimuth velocity, then two refocused SAR images will be obtained. In the two refocused SAR images, the background is smeared to the same extent, while for a moving target, the two refocused Doppler spectra will be## (16)

$${S}_{f1}(u)=\mathrm{exp}[j\frac{{u}^{2}}{2\pi {K}_{a}}({\epsilon}_{a}-{\epsilon}_{a}^{\prime})]\mathrm{exp}(-ju{x}_{0}),$$## (17)

$${S}_{f2}(u)=\mathrm{exp}[j\frac{{u}^{2}}{2\pi {K}_{a}}({\epsilon}_{a}+{\epsilon}_{a}^{\prime})]\mathrm{exp}(-ju{x}_{0}).$$Equations (16) and (17) illustrate that the moving target will be smeared differently in the two refocused images. In one image it is focused better, while in another image it is defocused worse. For a quadratic phase error, in the smeared image the energy tends to be spread uniformly over the distance of the smear for an un-weighted aperture.^{2} The moving target image will be in better focus when the refocused residual phase term is smaller, with the opposite also being true.

Different measurements can be used to describe the focusing quality, such as sharpness, entropy, and contrast. We adopt the sharpness function^{12} defined by

## 2.2.

### Azimuth Velocity Component Measurement

In Eq. (16), the refocused image is smeared in

## (20)

$$M({\epsilon}_{a},{\epsilon}_{a}^{\prime})=1+|{\epsilon}_{a}-{\epsilon}_{a}^{\prime}|{\mu}_{x}$$## (21)

$$M({\epsilon}_{a},-{\epsilon}_{a}^{\prime})=1+|{\epsilon}_{a}+{\epsilon}_{a}^{\prime}|{\mu}_{x}$$## (22)

$${D}_{c}({\epsilon}_{a},{\epsilon}_{a}^{\prime})=\frac{{|b|}^{4}}{M({\epsilon}_{a},{\epsilon}_{a}^{\prime})}-\frac{{|b|}^{4}}{M({\epsilon}_{a},-{\epsilon}_{a}^{\prime})}.$$Equation (22) will be analyzed on the condition that ${\epsilon}_{a}^{\prime}>0$. When ${\epsilon}_{a}>0$, in the case of $0<{\epsilon}_{a}^{\prime}\le {\epsilon}_{a}$,

## (23)

$$\frac{\mathrm{d}{D}_{a}}{\mathrm{d}{\epsilon}_{a}^{\prime}}=\frac{{|b|}^{4}{\mu}_{x}}{{[1+{\mu}_{x}({\epsilon}_{a}-{\epsilon}_{a}^{\prime})]}^{2}}+\frac{{|b|}^{4}{\mu}_{x}}{{[1+{\mu}_{x}({\epsilon}_{a}+{\epsilon}_{a}^{\prime})]}^{2}}>0.$$And in the case of ${\epsilon}_{a}^{\prime}>{\epsilon}_{a}$,

## (24)

$$\frac{\mathrm{d}{D}_{a}}{\mathrm{d}{\epsilon}_{a}^{\prime}}=\frac{{|b|}^{4}{\mu}_{x}}{{[1+{\mu}_{x}({\epsilon}_{a}+{\epsilon}_{a}^{\prime})]}^{2}}-\frac{{|b|}^{4}{\mu}_{x}}{{[1+{\mu}_{x}({\epsilon}_{a}^{\prime}-{\epsilon}_{a})]}^{2}}<0.$$The two cases show that (1) when $0<{\epsilon}_{a}^{\prime}\le {\epsilon}_{a}$, ${D}_{c}({\epsilon}_{a},{\epsilon}_{a}^{\prime})$ is a monotonic increasing function that reaches its maximum value at the point where ${\epsilon}_{a}^{\prime}={\epsilon}_{a}$, and (2) when ${\epsilon}_{a}^{\prime}>{\epsilon}_{a}$, ${D}_{c}({\epsilon}_{a},{\epsilon}_{a}^{\prime})$ is a monotonic decreasing function that reaches its maximum value at the point where ${\epsilon}_{a}^{\prime}={\epsilon}_{a}$. In addition, it infinitely approaches zero with incremental increases of ${\epsilon}_{a}^{\prime}$.

The same analysis procedure can be repeated for the case of ${\epsilon}_{a}<0$. The analysis results show that (1) when $0<{\epsilon}_{a}^{\prime}\le -{\epsilon}_{a}$, ${D}_{c}({\epsilon}_{a},{\epsilon}_{a}^{\prime})$ is a monotonic decreasing function that reaches its minimum value at the point where ${\epsilon}_{a}^{\prime}=-{\epsilon}_{a}$, and (2) when ${\epsilon}_{a}^{\prime}>-{\epsilon}_{a}$, ${D}_{c}({\epsilon}_{a},{\epsilon}_{a}^{\prime})$ is a monotonic increasing function that reaches its minimum value at the point where ${\epsilon}_{a}^{\prime}=-{\epsilon}_{a}$. In addition, the sharpness difference curve infinitely approaches zero with the increment of ${\epsilon}_{a}^{\prime}$.

Figure 2 presents two sharpness difference curves versus ${\epsilon}_{a}^{\prime}$ for two moving targets with ${\epsilon}_{a}>0$ and ${\epsilon}_{a}<0$, respectively. From this diagram, we can see that the azimuth velocity of a moving target can be estimated from its sharpness difference curve at the peak or valley point. The peak indicates that the target’s azimuth velocity component is positive, while the valley indicates a negative value.

The discussion above provides us a chance to estimate the azimuth velocity of a detected moving target by varying ${\epsilon}_{a}^{\prime}$ in a reasonable range.

## 3.

## Implementation

The optical MTI processor architecture is presented in Fig. 3. In this architecture, a given complex-valued SAR image named $I(m,n)$ is separated into an amplitude part and a phase part in the computer, i.e.,

where $\phi (m,n)$ is the phase angle of $I(m,n)$. Two spatial light modulators (SLM), ${M}_{0}$ and ${M}_{1}$, are modulated by $|I(m,n)|$ and $\mathrm{exp}[j\phi (m,n)]$, respectively. The modulated light is split up into two beams by a beam-splitter. One beam is transmitted to the cylindrical lens ${L}_{2}$, the other part is reflected to the cylindrical lens ${L}_{4}$. The lenses ${L}_{2}$ and ${L}_{4}$ transform the complex-valued SAR image into a range Doppler map. The phase-only SLM ${M}_{2}$ is modulated by the refocusing filter ${H}_{1}(u)$, and the filtered range Doppler map is transformed by the cylindrical lens ${L}_{3}$ resulting in a refocused SAR image. The computer samples the refocused image by using a charge-coupled diode (CCD) camera. The same operations are performed at ${L}_{4}$, ${M}_{3}$, and ${L}_{5}$ except that ${M}_{3}$ is modulated by the refocusing filter ${H}_{2}(u)$. After the two refocused SAR images are sampled and stored, their sharpness difference is computed and the moving targets are indicated.Let us discuss the selection of ${\epsilon}_{a}^{\prime}$. For a moving target, either a too small or a too large ${\epsilon}_{a}^{\prime}$ will lead to a small ${D}_{c}$. To detect moving targets with different azimuth velocities, ${\epsilon}_{a}^{\prime}$ will vary. For example, it can be orderly selected from a set ${V}_{a}^{\prime}=\{5,10,15,20,25,30\}\text{\hspace{0.17em}}\text{\hspace{0.17em}}(\mathrm{m}/\mathrm{s})$.

After a moving target is detected, it is isolated from the given SAR image and sent to the SLMs ${M}_{0}$ and ${M}_{1}$. By varying ${\epsilon}_{a}^{\prime}$, a series of refocusing filter pairs are generated and displayed in ${M}_{2}$ and ${M}_{3}$. As a result, a sharpness difference curve is drawn in the computer. The target’s azimuth velocity can be estimated from this curve.

This scheme is feasible because the sharpness of clutter is erased from the sharpness difference image, and thus the moving target can be detected with greater ease. Furthermore, the optical architecture ensures high computational efficiency because complex Fourier transforms are implemented by lenses instead of digital computation.

## 4.

## Results

## 4.1.

### Computer Simulation

We can verify the proposed scheme by computer simulations combined with optical instruments. In one of the simulations, three point scatterers are used as the scene. The three scatterers are set up such that one is stationary, and the other two have azimuth velocity components of $-10$ and $10\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$. The radar waveform is a chirp pulse, its bandwidth is 400 MHz, the pulse width is 20 *μ*s, and the pulse repetition frequency is 2000 MHz. The radar carrier frequency is 10 GHz. While the radar platform is moving at the speed of $200\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$, the returns are collected and processed, resulting in an SAR image embedded with three targets labeled by ${T}_{1}$, ${T}_{2}$, and ${T}_{3}$ shown in Fig. 4(a). From this figure, we can see that the stationary scatterer ${T}_{1}$ is optimally focused, while ${T}_{2}$ and ${T}_{3}$ are smeared in the azimuth direction because of their azimuth velocity components. The complex-valued SAR image is put in the SLMs ${M}_{0}$ and ${M}_{1}$ and is defocused by ${M}_{2}$ and ${M}_{3}$. The SLMs ${M}_{2}$ and ${M}_{3}$ are modulated by ${H}_{1}(u)$ and ${H}_{2}(u)$ for ${\epsilon}_{a}^{\prime}=0.025$ (${V}_{a}^{\prime}=5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$), and the defocused SAR images are shown in Fig. 4(b) and 4(c), respectively. From the two defocused images, we see that the scatterer ${T}_{1}$ is defocused to the same extent, while in Fig. 4(b), ${T}_{2}$ is smeared worse and ${T}_{3}$ is focused better, and in Fig. 4(c), ${T}_{2}$ is focused better and ${T}_{3}$ is smeared worse. So the moving targets ${T}_{2}$ and ${T}_{3}$ can be indicated by comparing the sharpness difference of the two defocused SAR images.

The azimuth velocity components of the three detected scatterers are estimated according to Eq. (22) by using the optical processor. The measured sharpness difference curves varying with ${V}_{a}^{\prime}$ are presented in Fig. 5. In the estimation process, ${V}_{a}^{\prime}$ ranges from 0 to $20\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$ with regular increments of $0.5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$. The sharpness difference curve of ${T}_{1}$ is flat because a stationary target tends to be defocused to the same extent in the SLMs ${M}_{2}$ and ${M}_{3}$. The sharpness difference curve of ${T}_{2}$ reaches its valley value at ${V}_{a}^{\prime}=10\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$, and it indicates that the target’s azimuth velocity is $-10\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$. The sharpness difference curve of ${T}_{3}$ reaches its peak value at ${V}_{a}^{\prime}=10\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$, and it indicates that the azimuth velocity of ${T}_{3}$ is $10\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$. From Fig. 5 we see that the maximum sharpness difference values appear at the expected azimuth velocity points. The simulation results show that the proposed scheme is feasible and effective.

During our simulations, we find that the scheme provides high measurement precision because of the sharp and high sharpness difference curve for a large azimuth velocity, while for a small azimuth velocity, the measurement precision is poor because its sharpness difference curve is fairly flat and low.

## 4.2.

### Field Data

A set of field SAR data is used to verify the optical processor. Figure 6(a) shows an image of two vehicles moving along a runway at the speed of $5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$ in the same direction. The targets are labeled by ${T}_{1}$ and ${T}_{2}$, respectively. The radar waveform parameters are listed as follows. The carrier frequency is 9.6 GHz, the pulse width is 20 *μ*s, the bandwidth is 400 MHz, and the pulse repetition frequency is 1200 MHz. The platform moves at the speed of $218\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$ along the runway. From Fig. 6(a) we see that the images of the two vehicles are smeared in the azimuth direction.

In the target detection process, the SLMs ${M}_{2}$ and ${M}_{3}$ are modulated by ${H}_{1}(u)$ and ${H}_{2}(u)$ with ${\epsilon}_{a}^{\prime}=0.029$ (${V}_{a}^{\prime}=5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$), respectively. After being filtered by ${M}_{2}$ and ${M}_{3}$, two defocused SAR images are generated and shown in Fig. 6(b) and 6(c), respectively. It can be seen that the two vehicles are smeared even worse in Fig. 6(b) and focused better in Fig. 6(c). By computing the sharpness difference of the two defocused images, the moving vehicles can be indicated.

The azimuth velocity components of the two detected vehicles are measured according to Eq. (22) by using the proposed optical processor. The measured sharpness difference curves varying with ${V}_{a}^{\prime}$ are presented in Fig. 7. In the estimation process, ${V}_{a}^{\prime}$ ranges from 0 to $10\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$ with regular increments of $0.5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$. The two curves reach their valley values at ${V}_{a}^{\prime}=5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$ approximately. It indicates that the speed of the vehicles is about $-5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{m}/\mathrm{s}$.

The curves in Fig. 7 show a nonmonotonic behavior due to the radar imaging process. Uncontrollable measurement conditions, such as the vibration of antenna, nonstraight trajectory of the platform, external interference, and digitalization, will cause phase errors in the echoes during the signal collection process. So there are some subtle differences between the two defocused images of the background, which is mixed in the detected moving target. As a result, the sharpness difference curve of a real moving target appears uneven. However, the trends of the two curves in Fig. 7 are similar to that in Fig. 2. This experimental result shows that the proposed scheme can deal with the field data practically.

## 5.

## Conclusion

In summary, we have presented an optical MTI to detect moving targets and estimate their azimuth velocities in SAR images. In principle, it is robust against clutter and interference. The experimental results using both simulated and field data show that the proposed optical architecture is feasible, effective, and practical.

## Acknowledgments

We appreciate the support from the Shandong Young Scientists Award Foundation (Grant Number: BS2010DX021) and the National Natural Science Foundation of China (Grant Number: 61070175 and 61175053).

## References

## Biography

**Yuan Li** received her BSc degree from the Qufu Normal University, Shandong, China, in 1998, a MSc degree from University of Electronic Science and Technology of China, Sichuan, China, in 2001 and a PhD degree from the Chinese Academy of Engineering Physics (CAEP), Beijing, China, in 2007. From 2007 to 2010, she was an engineer with IEE of CAEP, Mianyang, China, where she worked on signal processing and ISAR systems. She is currently a teacher in Shandong Institute of Business and Technology, and a researcher in the Key Laboratory of Intelligent Information Processing in Universities of Shandong, Yantai, China. Her research interests include coherent optical signal processing, radar signal processing, and ground moving target indication.

**Gaohuan Lv** received his BSc degree from the Qufu Normal University, Shandong, China, in 1998, a MSc degree from the Chinese Academy of Engineering Physics, Beijing, China, in 2001 and a PhD degree from the Shanghai Jiao Tong University, Shanghai, China, in 2013. From 2001 to 2008, he was an engineer with IEE of CAEP, Mianyang, China, where he worked on control system and software engineering. His research interests include coherent optical signal processing, radar signal processing, digital image processing, and ground moving target indication.