|
1.IntroductionGround 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–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–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.Principles2.1.Moving Target Detection in Space DomainA strip-map mode SAR is considered in the following analysis. For a stationary point scatterer, its azimuth phase history has the form:10 where is the azimuth coordinate, and with being the carrier wavelength and being the slant range to the stationary scatterer. If a scatterer is moving in azimuth direction with a constant velocity , then the azimuth signal becomes where is the azimuth position of the moving scatterer when the radar bore-sight is directed to the scatterer, is the synthetic aperture length, and with being the platform velocity.Taking the Fourier transform of according to the stationary phase theory and ignoring the amplitude constant and the phase constant, one obtains where is the antenna length. Usually , so the following approximation can be made:In the traditional range Doppler algorithm,11 if is filtered by then the focused version of isTaking the inverse Fourier transform of the above in Eq. (6), the focused signal can be shown to be where and the symbol means the convolution operation. In Eq. (9), denotes the inverse Fourier transform. According to the stationary phase theory, and neglecting the constant coefficient, takes the formEquation (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 and that of the target with have the pulse width of and , respectively. As the bandwidth of is , its azimuth width in space domain satisfies and thus where with being the theoretical azimuth resolution which is the effective space width of . The image of a moving scatterer will be smeared to the azimuth lengthIf the SAR image is refocused in azimuth with a pair of new refocusing filters described by and where with 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 andEquations (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 function12 defined by where is a patch sized by a resolution cell. The moving target can be indicated by comparing the sharpness between the two refocused SAR images: where and are the sharpness distribution image of the two SAR images refocused by and , respectively. During the sharpness computation procedure, the patches are overlapping so that a moving target with a typical size will correspond to a set of bright pixels in the sharpness difference image. The bright pixels will be kept when a median filter is used to remove isolated bright pixels caused by noise.2.2.Azimuth Velocity Component MeasurementIn Eq. (16), the refocused image is smeared in azimuth resolution cells, and in Eq. (17), the refocused image is smeared in azimuth resolution cells. For a moving scatterer having the intensity , its sharpness difference of the two refocused SAR images can be formulated byEquation (22) will be analyzed on the condition that . When , in the case of , And in the case of , The two cases show that (1) when , is a monotonic increasing function that reaches its maximum value at the point where , and (2) when , is a monotonic decreasing function that reaches its maximum value at the point where . In addition, it infinitely approaches zero with incremental increases of . The same analysis procedure can be repeated for the case of . The analysis results show that (1) when , is a monotonic decreasing function that reaches its minimum value at the point where , and (2) when , is a monotonic increasing function that reaches its minimum value at the point where . In addition, the sharpness difference curve infinitely approaches zero with the increment of . Figure 2 presents two sharpness difference curves versus for two moving targets with and , 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 in a reasonable range. 3.ImplementationThe optical MTI processor architecture is presented in Fig. 3. In this architecture, a given complex-valued SAR image named is separated into an amplitude part and a phase part in the computer, i.e., where is the phase angle of . Two spatial light modulators (SLM), and , are modulated by and , respectively. The modulated light is split up into two beams by a beam-splitter. One beam is transmitted to the cylindrical lens , the other part is reflected to the cylindrical lens . The lenses and transform the complex-valued SAR image into a range Doppler map. The phase-only SLM is modulated by the refocusing filter , and the filtered range Doppler map is transformed by the cylindrical lens 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 , , and except that is modulated by the refocusing filter . 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 . For a moving target, either a too small or a too large will lead to a small . To detect moving targets with different azimuth velocities, will vary. For example, it can be orderly selected from a set . After a moving target is detected, it is isolated from the given SAR image and sent to the SLMs and . By varying , a series of refocusing filter pairs are generated and displayed in and . 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.Results4.1.Computer SimulationWe 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 and . 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 , the returns are collected and processed, resulting in an SAR image embedded with three targets labeled by , , and shown in Fig. 4(a). From this figure, we can see that the stationary scatterer is optimally focused, while and are smeared in the azimuth direction because of their azimuth velocity components. The complex-valued SAR image is put in the SLMs and and is defocused by and . The SLMs and are modulated by and for (), 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 is defocused to the same extent, while in Fig. 4(b), is smeared worse and is focused better, and in Fig. 4(c), is focused better and is smeared worse. So the moving targets and 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 are presented in Fig. 5. In the estimation process, ranges from 0 to with regular increments of . The sharpness difference curve of is flat because a stationary target tends to be defocused to the same extent in the SLMs and . The sharpness difference curve of reaches its valley value at , and it indicates that the target’s azimuth velocity is . The sharpness difference curve of reaches its peak value at , and it indicates that the azimuth velocity of is . 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 DataA 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 in the same direction. The targets are labeled by and , 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 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 and are modulated by and with (), respectively. After being filtered by and , 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 are presented in Fig. 7. In the estimation process, ranges from 0 to with regular increments of . The two curves reach their valley values at approximately. It indicates that the speed of the vehicles is about . 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.ConclusionIn 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. AcknowledgmentsWe 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). ReferencesR. Raney,
“Synthetic aperture imaging radar and moving targets,”
IEEE Trans. Aerosp. Electron. Syst., AES-7
(3), 499
–505
(1971). http://dx.doi.org/10.1109/TAES.1971.310292 IEARAX 0018-9251 Google Scholar
J. Fienup,
“Detecting moving targets in SAR imagery by focusing,”
IEEE Trans. Aerosp. Electron. Syst., 37
(3), 794
–809
(2001). http://dx.doi.org/10.1109/7.953237 IEARAX 0018-9251 Google Scholar
J. WangX. Liu,
“Velocity estimation of moving targets using SAR,”
in Proc. IEEE IGARSS,
340
–343
(2011). Google Scholar
P. MarquesJ. Dias,
“Moving targets processing in SAR spatial domain,”
IEEE Trans. Aerosp. Electron. Syst., 43
(3), 864
–874
(2007). http://dx.doi.org/10.1109/TAES.2007.4383579 IEARAX 0018-9251 Google Scholar
D. Weihinget al.,
“Detection of along-track ground moving targets in high resolution spaceborne SAR images,”
ISPRS J. Photogramm. Remote Sens., 61
(3–4), 135
–140
(2006). http://dx.doi.org/10.1016/j.isprsjprs.2006.09.011 IRSEE9 0924-2716 Google Scholar
R. Bamler,
“Doppler frequency estimation and the Cramer-Rao bound,”
IEEE Trans. Geosci. Remote. Sens., 29
(3), 385
–390
(1991). http://dx.doi.org/10.1109/36.79429 IGRSD2 0196-2892 Google Scholar
I. CummingS. Li,
“Improved slope estimation for SAR Doppler ambiguity resolution,”
IEEE Trans. Geosci. Remote. Sens., 44
(3), 707
–718
(2006). http://dx.doi.org/10.1109/TGRS.2005.861925 IGRSD2 0196-2892 Google Scholar
F. McFadden,
“Precise pose estimation for synthetic aperture radar images of vehicles,”
Opt. Eng., 46
(10), 107201
(2007). http://dx.doi.org/10.1117/1.2793710 OPEGAR 0091-3286 Google Scholar
D. PsaltisK. Wagner,
“Real-time optical synthetic aperture radar (SAR) processor,”
Opt. Eng., 21
(5), 822
–828
(1982). http://dx.doi.org/10.1117/12.7972988 OPEGAR 0091-3286 Google Scholar
F. DickeyJ. Mason,
“Optical Synthetic-aperture radar processor architecture with quadratic phase-error correction,”
Opt. Lett., 15
(20), 1162
–1164
(1990). http://dx.doi.org/10.1364/OL.15.001162 OPLEDP 0146-9592 Google Scholar
B. Wang, Digital Signal Processing Techniques and Applications in Radar Image Processing, 254
–255 Wiley, Hoboken, NJ
(2008). Google Scholar
J. Fienup,
“SAR autofocus by maximizing sharpness,”
Opt. Lett., 25
(4), 221
–223
(2000). http://dx.doi.org/10.1364/OL.25.000221 OPLEDP 0146-9592 Google Scholar
BiographyYuan 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. |