A Modified Stripmap SAR Processing for Vector Velocity Compensation Using the Cross-Correlation Estimation Method

Article information

J. Electromagn. Eng. Sci. 2019;19(3):159-165
Publication date (electronic) : 2019 July 31
doi : https://doi.org/10.26866/jees.2019.19.3.159
1School of Electrical Engineering, Korea Advanced Institute of Science and Technology, Daejeon, Korea
2Satellite Technology Research Center, Korea Advanced Institute of Science and Technology, Daejeon, Korea
*Corresponding Author: Seong-Ook Park (e-mail: soparky@kaist.ac.kr)
Received 2018 October 11; Revised 2018 December 15; Accepted 2019 May 07.

Abstract

This paper proposes a modified stripmap synthetic aperture radar (SAR) signal processing algorithm with an X-band SAR system, using the practical measurement results of automobile-based SAR (Auto-SAR). The quality of the image is degraded due to an unexpected direction change or anisotropy of the target position in radar image measurement. To solve the quality problem, signal processing is required for the vector velocity of each range in the azimuth direction. An X-band chirp pulse system was implemented and optimized by a signal processing algorithm suitable for high resolution. The stripmap SAR images are produced in five places. The validity of the proposed algorithm is verified by comparing impulse response function analysis and experimental images.

I. Introduction

Synthetic aperture radar (SAR) imaging technology is one of the most promising technologies utilized in various observation areas such as terrain, resource, and target searches. With the increasing demand for imaging technology, SAR systems have been developed to achieve high resolution and wide swath for accurate target detection, regardless of weather or environmental changes. To obtain high-resolution images, the SAR system can be improved in two ways. The first is a budget parameter optimization of the transmission and reception systems to reduce the ambiguity of radar signals and various interference problems. The second is to use diverse SAR signal processing algorithms for high-resolution image extraction. However, considering the cost and physical volume of the SAR system, an effective way to obtain high-resolution images is with diverse signal processing techniques.

There are numerous signal processing algorithms such as the range Doppler algorithm (RDA), chirp scaling algorithm (CSA), omega-K algorithm (ω-k), back projection algorithm (BPA). The RDA is widely used for efficient and stable signal processing because it is simple and easy to use. However, when the RDA is used for image processing, the final image is not focused due to random environmental variables. To cope with environmental variables and acquire high-resolution images, precise estimations of the Doppler centroid frequency and chirp rate are required in the algorithm. Many researchers have proposed algorithms to correct the received signal by removing the unexpected signals, using the magnitude or phase information of the chirp signal based on average cross-correlation coefficient (ACCC), multi-look cross correlation (MLCC), and multi-look beat frequency (MLBF) methods [1, 2]. However, the Doppler parameters extracted through previous algorithms still have uncertainty because of topographic change of interest area and interference of unexpected signal. Therefore, many researchers have attempted to increase the accuracy of Doppler parameters in various SAR experiments [36]. This paper proposes a novel method to solve the uncertainty by considering uncorrected Doppler parameters. In conventional SAR experiments, signal processing was performed by estimating the constant velocity with the lowest ambiguity error. This conventional method can obtain a reasonable image quality over the entire area if the area of interest has a nearly constant height over whole region. However, if the difference of altitude in the area and the direction of measurement is changeable, a high-quality SAR image cannot be obtained with only a constant velocity. Based on the changes of velocity in the area of interest, the quality of a SAR image can be improved using cross-correlation theory, designed for SAR. The estimated velocity is used to acquire a high resolution SAR image without complicated processing techniques. To evaluate the proposed method, the raw data obtained from practical experiments are used. Finally, high-resolution SAR images are constructed, applying the proposed algorithm to various environmental changes of automobile-based SAR (Auto-SAR).

II. SAR Signal Processing Analysis of the RDA

The RDA is a method of compressing raw data using a matched filter in the range direction and azimuth direction, respectively. As shown in Fig. 1, the RDA is typically used in stripmap standard mode with chirp pulse signals. The basic signal processing procedure consists of range and azimuth compressions and range cell migration, as shown in Fig. 2. Doppler centroid estimation can be applied to obtain high-resolution images in the SAR ground experiment. In case of the ground experiment, the image difference caused by the Doppler centroid variation is insignificant [1]. Therefore, it is possible to exclude the Doppler centroid estimation in a mathematical mechanism for SAR signal processing. The chirp pulse signals in the range and azimuth directions received from the SAR system are as follows, respectively.

Fig. 1

Stripmap SAR mode for chirp pulse radar.

Fig. 2

Conventional RDA.

(1) Srange(τ,η)=Arωr[τ-2R(η)c]ej4πf0R(η)cejπKr[τ-2R(η)c]2
(2) Sazimuth(τ,η)=Aaωa(η-ηc)ej2πffηc(τ)ejπKa(τ)η2

τ and η are the time, ω r and ω a are the window functions each chirp pulse signal, A r and A a are the amplitude constant values, f 0 and f η c are the center frequency of chip pulse signal, and K r and K a are the chirp rate of the signals. The first and second parameters mentioned above are for the range and azimuth direction, respectively. R(η) in (1) is the instantaneous slant range. Eqs. (1) and (2) are approximated as (3) through the range and azimuth compressions and range cell migration. After the range cell migration of SAR data, the performance of resolution for the matched filter in the azimuth direction can be determined according to K a and f η c . The final SAR image signal, S final (τ,η), is as follows.

(3) Sfinal(τ,η)=A0pr(τ-2R0c)pa(η)e(-j4πf0R0c)ej2πfηcη

Therefore, it is more important to estimate the exact matched filter in the azimuth direction than the range direction because the errors of K a and f η c are vulnerable by the SAR measurement environment, such as the actual speed, changing rate in the distance to the target, and motion vibrating conditions. However, the matched filter in the range direction has a low error rate because the raw data are compressed using the system parameters or loop-back chirp pulse signal. Therefore, our proposed algorithm estimates the accurate velocity in each range, which is an important parameter for the azimuth matched filter.

III. Proposed SAR Processing Algorithm

Based on the RDA, the modified signal process of the Doppler parameter is proposed in this section. To estimate the azimuth matched filter that fits with the actual signal, the exact Doppler centroid, f η c (τ) and Doppler chirp rate, K a (τ) are required under variable terrain conditions [79]. The Doppler parameters are as follows.

(4) fηc(τ)=2Vssin(θs)λ
(5) Ka(τ)2Vs2cos3θsλR0(τ)

θ s is the squint angle for measuring the target, λ is the wavelength of the center frequency in the chirp pulse signal, and R 0(τ) is the slant range along to the range time. The Doppler parameters are addressed as a function of the vehicle velocity, V s . Therefore, the main reason for the uncorrected azimuth matched filter is the misestimated velocity. This velocity leads to the low quality of the SAR image. The velocity estimation of each range can lead to more the high-resolution quality and precise target classification than the conventional results of the SAR processing. To find the accurate velocity V s , the proposed algorithm incorporates cross-correlation technique into the estimation of the Doppler parameter [10, 11]. If S(τ,η) is a 2D signal constructed in the range and azimuth, the cross-correlation function is defined as follows.

(6) (sLook1*sLook2)(τn,η)=def-sLook1*(τn)sLook2(τn,η+ω)dω         (n=1,2,3,4,,maximum range cell)

τ n is the slow time and, η is the fast time. Based on (6), the chirp pulse frequency of the range-Doppler domain along the range is divided into two looks (Look1, Look2), having V s + ΔV error in phase information. A matched filter of assumed V s multiplies each look, and then inverse.

FFT is used to generate the images of each look with a phase error. The cross-correlation values between two images are compared until reaching the smallest correlation value, changing V s . Range cell migration correction (RCMC) and azimuth compression are performed, using the V s defined from finishing the cross-correlation estimation in each range cell with phase variance compensation. In the conventional SAR process, each step of RDA uses its own technology to compensate the phase error for every range cell. It is time-consuming and complicated, resulting in low efficiency. However, in the proposed technique, all steps are performed at once without any other process using the estimated V s . Unlike airborne- and spaceborne-SAR, which are sensitive to environmental changes, the proposed algorithm is suitable for Auto-SAR to make real-time processing fast and easy, resulting in high efficiency.

From the estimated results, the defined V s can vary over a range, such as a vector. Therefore, the V s is defined as the vector velocity. The most suitable V s depends on the direction of the target. Before the practical experiment, impulse response function (IRF) analyses using the proposed algorithm are performed in Fig. 3. Fig. 3(a) and (b) show the analysis structure in 2D and the Doppler chip pulse signal according to the range, respectively. Fig. 3(c) shows that estimating the correct V s improves the peak value and the 3-dB resolution of IRF of the target in the SAR image. Fig. 4 shows over whole signal processing of the proposed algorithm by organizing the process described in this section.

Fig. 3

Information of the proposed estimation: (a) the method from the raw data structure for vector velocity, (b) the chirp pulse in the azimuth along to the range, and (c) the difference in the ideal IRF due to the velocity variation.

Fig. 4

Proposed RDA for high-resolution stripmap SAR.

IV. Experiment Results and Performance

To evaluate the performance of the proposed algorithm, the X-band SAR system was implemented and the practical experiment was conducted. Table 1 shows the specifications of the X-band SAR system.

X-band SAR system specifications

As shown in Fig. 5, two antennas are installed as the transmitter and receiver, respectively. The stop-and-go method is applied to measure the practical SAR experiment. To prove the validity of the proposed SAR algorithm, five different regions are selected. In particular, two curved terrains are included in the practical experiment regions. Fig. 6 shows the practical measurement locations, which are taken while driving through an express highway. Scenes #1 and #2 are places to drive through curved loads, Scene #3 is a lower elevation terraced region as the range increases along the distance. Scenes #4 and #5 are rice fields measured in both directions at a height of 60 m. Applying vector velocity estimation, the quality of the SAR image is significantly improved over conventional RDA in all regions, as shown in Fig. 7. Fig. 8 shows the difference before and after applying the vector velocity to Scenes #1 and #2, respectively. As shown in the images in Fig. 8, we can see the performance of the algorithm through the focusing improvement in practical experiments. In addition, as confirmed in Figs. 7 and 8, the proposed technique can be applied to various terrain and experimental conditions in Auto-SAR. In Fig. 9, analysis of the corner reflector in Scene #5 shows that the vector estimation improves the 3-dB resolution of the IRF and increases the peak energy of the main lobe. Through the estimated velocity is closer to the accurate velocity, the energy balance of the sinc-function is addressed more clearly. These results suggest that the desired resolution and high-quality SAR image can be obtained by the correction of the vector velocity.

Fig. 5

The picture of the SAR measurement setup in an automobile.

Fig. 6

The measurement area for verifying proposed SAR (location of Scenes #1 and #2, Gyeo-nam bridge; Scene #3, Chonchon bridge; Scenes #4 and #5, Magoksa IC).

Fig. 7

Signal processing results by the proposed algorithm: (a) Scene #1, (b) Scene #2, (c) Scene #3, (d) Scene #4, and (e) Scene #5.

Fig. 8

The results of Scenes #1 and #2 when processed by the original RDA (a), (c) and the proposed RDA (b), (d), respectively.

Fig. 9

IRF analysis using the proposed algorithm: (a) the location and image of corner reflector and (b) the impulse response function of a corner reflector according to velocity.

V. Conclusion

The proposed algorithm is suitable for measuring the image in Auto-SAR. The recent technology has complex process to compensate for the phase error in each range cell for airborne-and spaceborne-SAR. It is inefficient in creating SAR images in Auto-SAR. However, the proposed technology is easy to handle and provides fast performance, allowing an immediate response to a variety of situations that can occur on the ground. Applying the precise estimation of the vector velocity can provide more accurate signal processing for high-quality SAR images. To prove the effectiveness of the proposed algorithm, high-quality results are acquired from practical experiments. In addition, the proposed technique of vector velocity is expected to be applied in various conditions, such as sub-aperture correction, spotlight mode, and autofocus.

Acknowledgments

This work was supported by Institute for information & communications Technology Promotion (IITP) grant funded by the Korea government (MSIT) (No. 2018-0-01658, Key Technologies Development for Next Generation Satellites).

References

1. Cumming IG, Wong FH. Digital Processing of Synthetic Aperture Radar Data London: Artech House Inc; 2005.
2. Curlander JC, McDonough RN. Synthetic Aperture Radar Systems and Signal Processing New York, NY: John Wiley & Sons Inc; 1991.
3. Liu B, Liu C, Wang Y. Slice convolution based slope estimation for SAR Doppler ambiguity resolver. In : Proceedings of 2013 IEEE International Geoscience and Remote Sensing Symposium. Melbourne, Australia; 2013; p. 2055–2058.
4. Yang J, Zhang Y, Kang X. A Doppler ambiguity tolerated algorithm for airborne SAR ground moving target imaging and motion parameters estimation. IEEE Geoscience and Remote Sensing Letters 12(12):2398–2402. 2015;
5. Liu Q, Wang Z, Xia W, Chen B. Ground moving target positioning with Doppler ambiguity and scene ambiguity in SAR. In : Proceedings of 2017 40th International Conference on Telecommunications and Signal Processing (TSP). Barcelona, Spain; 2017; p. 465–469.
6. Li W, Huang Y, Yang J, Wu J. Comparison of geometry-based Doppler ambiguity resolver in squint SAR. In : Proceedings of 2011 IEEE International Geoscience and Remote Sensing Symposium. Vancouver, Canada; 2011; p. 4443–4445.
7. Zheng H, Wang J, Liu X, Gao Y, Zhang L. Velocity estimation of the moving target for high-resolution wide-swath SAR systems. In : Proceedings of 2017 IEEE International Geoscience and Remote Sensing Symposium (IGARSS). Fort Worth, TX; 2017; p. 2054–2057.
8. Liu X, Xiong J, Huang Y, Wu J, Yang J, Pu W. SAR Doppler frequency rate estimation based on time-frequency rate distribution. In : Proceedings of 2009 2nd Asian-Pacific Conference on Synthetic Aperture Radar. Xian, China; 2009; p. 926–930.
9. Moreira A, Mittermayer J, Scheiber R. A SAR autofocus technique based on azimuth scaling. In : Proceedings of 1997 IEEE International Geoscience and Remote Sensing Symposium Proceedings Remote Sensing-A Scientific Vision for Sustainable Development. Singapore; 1997; p. 2028–2030.
10. Yue Y, Zhang X, Jun S. Motion compensation method for translational variant bistatic SAR using autocorrelation variance. In : Proceedings of IET International Radar Conference. Guilin, China; 2009;
11. Ouchi K, Hwang SI. Improvement of ship detection accuracy by SAR multi-look cross-correlation technique using adaptive CFAR. In : Proceedings of 2010 IEEE International Geoscience and Remote Sensing Symposium. Honolulu, HI; 2010; p. 3716–3719.

Biography

Chul-Ki Kim was born in Gangneung, Korea, in July 1989. He obtained his B.S. degree in Electronic Engineering from Soongsil University, Seoul, Korea, in 2014, his M.S. degree in Electrical Engineering from Korea Advanced Institute of Science and Technology, Daejeon, in 2016, and is currently the candidate-student, working toward a Ph.D. in Electrical Engineering at Korea Advanced Institute of Science and Technology (KAIST). His current research interests include synthetic aperture radar (SAR).

Jang-Soo Chae was born in Suncheon, Korea, in August 1959. He obtained his B.S. degree from In Ha University, Korea in 1987, his M.S. degree from Seoul National University, Seoul, Korea in 1989, his Ph.D. from POSTECH in 1992, and another Ph.D. from Ajou University, Suwon, in 2004. From March 1990 to August 1995, he was a Research Engineer with Korea Aerospace Research Institute, Daejeon, working with LEO satellite system engineering, and is currently working at Satellite Research Center in KAIST. He is engaging in research on the development of passive antenna for synthetic aperture radar (SAR) onboard a small satellite and dynamics and the control of flexible and space structure.

Jung-Su Lee was born in Busan, Korea in February 1978. He obtained his B.S. degree in Electronic Engineering from Korea University, Seoul, Korea in 2001 and his M.S. degree in Radio Sciences and Engineering, Korea University, Seoul, Korea in 2003. Since July 2003, he has been a Senior Research Engineer at the Satellite Technology Research Center, Korea Advanced Institute of Science and Technology, Daejeon, Korea. His research interests include satellite wireless communication and ranging systems.

Seong-Ook Park was born in Kyungpook, Korea, in December 1964. He obtained his B.S. degree from Kyungpook National University, Korea, in 1987, his M.S. degree from Korea Advanced Institute of Science and Technology, Daejeon, Korea, in 1989, and his Ph.D. degree from Arizona State University, Tempe, AZ, in 1997, all in electrical engineering. From March 1989 to August 1993, he was a Research Engineer with Korea Telecom, Daejeon, working with microwave systems and networks. He later joined the Telecommunication Research Center, Arizona State University, until September 1997. Since October 1997, he has been with the Information and Communications University, Daejeon, and currently is a Professor at the Korea Advanced Institute of Science and Technology. His research interests include mobile handset antenna and analytical and numerical techniques in the area of electromagnetics. Dr. Park is a member of Phi Kappa Phi.

Article information Continued

Fig. 1

Stripmap SAR mode for chirp pulse radar.

Fig. 2

Conventional RDA.

Fig. 3

Information of the proposed estimation: (a) the method from the raw data structure for vector velocity, (b) the chirp pulse in the azimuth along to the range, and (c) the difference in the ideal IRF due to the velocity variation.

Fig. 4

Proposed RDA for high-resolution stripmap SAR.

Fig. 5

The picture of the SAR measurement setup in an automobile.

Fig. 6

The measurement area for verifying proposed SAR (location of Scenes #1 and #2, Gyeo-nam bridge; Scene #3, Chonchon bridge; Scenes #4 and #5, Magoksa IC).

Fig. 7

Signal processing results by the proposed algorithm: (a) Scene #1, (b) Scene #2, (c) Scene #3, (d) Scene #4, and (e) Scene #5.

Fig. 8

The results of Scenes #1 and #2 when processed by the original RDA (a), (c) and the proposed RDA (b), (d), respectively.

Fig. 9

IRF analysis using the proposed algorithm: (a) the location and image of corner reflector and (b) the impulse response function of a corner reflector according to velocity.

Table 1

X-band SAR system specifications

Parameter Value
Center frequency X-band
Frequency bandwidth (MHz) 300
Pulse repetition frequency (Hz) 1,000
Measurement velocity (km/hr) 80