Document Type : Research Paper
Authors
1 department of electrical and electronic engineering- Malek Ashtar university of technology
2 Department of electrical and electronic engineering- malek ashtar university of technology
3 department of electrical and electronic engineering- malek ashtar university of technology
Abstract
Keywords
The synthetic aperture radar (SAR) technology is a major development of the modern active coherent microwave remote sensing technology. Compared with the optical system, SAR can not only acheive highe resolution 2-D images, but also has the characteristics of all weather and all time with the application in the field of remote sensing, explorations, agriculture, military, etc.[1-3]. This coherent procesing also known as the compression of frequency modulated pulse allows more energy to be transmitted and also better range resolution to be achieved[4]. SARs are available in different application modes which depends on the type of operation, target condition and image resolution[5]. Fig.1 shows a low squinted stripmap SAR with its simpled receiver model as one of the common types of airborn SARs with high imaging resolution. Due to its versatile applications, airborns with unusual flight characteristics can be equipped with this technique like missiles[6-9].
Fig. 1. The low squinted (α= 0 to 6 Deg) missile-borne stripmap SAR flight geometery and receiver model
The most common type of puls compression waveform is Linear Frequency Modulation (LFM) [6-9]. The match filter response of LFM has a side-lobe level of -13dB, which can be improved by some methods but at costs reducing SNR and wider mainlobe, which make it improper to use in SAR image formation [10]. Non-linear LFM (NLFM) is another pulse compression technique which can provides finer resolution, better SNR and good interference mitigation [5, 11-17].It has a spectrum weighting function inherently in their modulation function, which could shape the power spectral density (PSD) of signal due to its instantaneous frequency as a non-linear function of time and result in a pure match filter giving low side-lobes[18]. The NLFM technique has been used such as discrete and continous signal modulation [11-12]. Continuous NLFM waveforms with the help of phase modulation were devised in a way that provides very smaller sidelobe level, as well as the ability to be designed on hardware structures in contrast to discrete signal modulation [13-17]. In [13], the NLFM is designed based on piecewise linear functions as a method of instantaneous frequency function to generate the NLFM while no quality analysis is done for the NLFM. The radar ambiguity function (AF) as a complete quality assessment statement of the waveform provides an indication of the limitations and utility of different waveforms while assess the properties of the transmitted waveform regarding to their target resolution, measurement accuracy, ambiguity and responses to the backscattering [11, 14, 20]. In [13, 15-17], the NLFM is just used in an image formation algorithm (IFA) and no further AF quality analysis on the waveform design procedure is done. In this study, NLFM waveform design consideration based on AF quality analysis is firstly presented. Then, in comparison to [12-17] two different NLFM waveforms are not only simulated and evaluated based on AF quality analysis, but also were directly implemented to the SAR IFA for further quality assessment metrics extraction. The mentioned quality assessment techniques on the basis of the 2D image have not been done yet [11-17]. The qualitative AF analysis of the optimal SAR continuous NLFM waveform is done via ambiguity diagrams extraction [19], while the quantitative metric assessment is done via the calculation of mean square errors (MSE), signal-to-noise ratio (SNR), peak SNR (PSNR), PSD, normalized cross correlation (NCC) and integrated side lobe ratio (ISLR) [20-21]. It should be noted that the effective antenna pointing angle on the basis of ISLR value is also extracred.
This paper is organized as follows. Firstly, both the LFM and the continuous NLFM waveforms are modeled based on correlated AF analysis. Section II, also introduces the proposed optimum continuous NLFM waveform selection based on AF quality analysis as well as the raw data generation. Section III, utilizes a range-Doppler algorithm (RDA) to generate the images of the proposed optimum cosine NLFM (CNLFM), while other images is derived. Section IV, contains complete quality metrics assessment techniques for evaluation of CNLFM waveform application in 2D SAR images while system parameters is being changed. SectionV, presents conclusions on the simulation results and proposed topics for further researches.
On the basis of AF criterion of choice, this section includes the LFM and the NLFM waveforms simulations.
The pulse compression method is the best technique to increase range resolution meanwhile SNR remains suitable in the receive process. The general form of LFM waveform is:
(1)
Where, T represents the radar pulse width, is the carrier frequency, is the chirp rate or the slope of the frequency alteration and is the signal amplitude modulation term, while the quadratic phase modulation presented as (2):
(2)
The LFM waveform spectrum with its chirp rate alteration is shown in Fig. 2. If the LFM waveform is designed in a way that its frequency alters asymmetrically within the bandwidth, the quadratic phase modulation will be changed as (3) which is called continuous NLFM:
(3)
Where defines the non-linearity behavior of quadratic phase modulation of the waveform in shape of sine, cosine, tangent and etc. [22]. The designing of the NLFM waveform shapes the PSD to approach a required weighting function and cause the waveform does not have the SNR loss in contrast to LFM [11]. Fig.3 shows how the frequency alters with the time while the sine, the cosine and conventional linear functions are used. So, the NLFM has better detection rate characteristics than LFM and supposed to be more accurate in range determination for SAR application [23].
Fig. 2. SAR LFM waveform presentation [B=50MHz and T=10 and chirp rate=5 1012]
Fig. 3. Time-Frequency presentation of LFM and NLFM waveforms upchirp and downchirp state
The aim of the next section is to synthesize the LFM, continuous SNLFM and continuous CNLFM waveforms based on AF diagrams extraction.
Radar waveform AF is the output response of a single point scatterer at all possible combination of ranges and velocities. The AF of the SAR waveform can be defined in terms of cross-correlation of Doppler-shifted version of the transmitted waveform, , in the delay time and Doppler frequency plane, ( , as Eq. (4):
(4)
Where, is the delay time and is the Doppler frequency shift. According to Eq. (1) and Eq. (4), the AF of the LFM pulse can finally be written as Eq. (5):
(5)
This relation represents both the LFM ambiguity diagram and the LFM ambiguity contour levels which are shown in Fig. 4. According to Eq. (1), Eq. (4) and Eq. (5) the LFM waveform AF diagrams are derived
Fig. 4. 3D-view of ambiguity diagram and its pertinent ambiguity contour level
as Fig. 5. The same AF analysis including, ambiguity diagram, ambiguity contour levels, AF spectrogram and AF frequency profile for both of the continuous SNLFM and the continuous CNLFM are carried out, which are shown in Fig.6. It is deduced that different frequency chirp rates of the NLFM waveforms will lead to different ambiguity contour levels while presenting the similar ambiguity diagram is expected, because of similar inherent specifications of NLFM in comparison to LFM counterpart. The results are shown in Fig.5 (a), 5(b) and 6(a), 6(b).
According to Fig. 6 (c) , (d) and 6(e), (f) , the AF contour levels of CNLFM not only differs from the SNLFM but also is slightly more compact than LFM waveform which results in less spectral expansion of adjacent pixels and helps to improve the image quality. The NLFM waveform shapes the PSD to approach a required weighting function with lower sidelobes, while the CNLFM and SNLFM frequency profile has greater main-lobe level in comparison to LFM. On the other hand, CNLFM has distinctive lower side lobe level ratio in comparison to the SNLFM which makes it suitable for SAR image formation application.The results are clearly shown in Fig. 6(g), (h), (i), (j).
On the basis of AF diagram concepts, the requirements for ambiguity diagrams and exact resolution accuracy of SAR based on NLFM waveforms generation are not possible to achieve simultaneously. Hence, according to AF qualitative analysis, the continuous CNLFM with its compact contour levels and lower side-lobes is more proper for implementation than LFM and SNLFM.
According to flight geometry and SAR receiver model in Fig. 1, the received baseband echo signal of the point target in the distance fast time and azimuth slow time domain can be stated as Eq. (6):
(6)
Where, is the range time, is the azimuth time, is the received back scattered signal from the target, is the received signal amplitude, is the received pulse envelope, is the transmitted pulse envelope in azimuth, is the carrier frequency, is the chirp rate and is the speed of the light. The presented analysis is on the supposition that the platform stability is very high and is unaffected by the
(a) (b)
(c) (d)
Fig. 5. LFM waveform AF simulation results (a) LFM waveform ambiguity diagram (b) LFM waveform ambiguity contour levels (c)LFM waveform AF frequency profile (d)LFM waveform AF spectrogram for (Pulse width=1 , Chirp rate=4 1013 down chirp)
(a) (b)
(c) (d)
(e) (f)
(g) (h)
(i) (j)
(k) (l)
Fig. 6. NLFM wave forms AF simulation results (a) CNLFM waveform ambiguity diagram (b) SNLFM waveform ambiguity diagram (c) CNLFM waveform ambiguity contour levels (d) SNLFM waveform ambiguity contour levels (e) LFM- CNLFM waveform ambiguity contour levels (f) LFM- SNLFM waveform ambiguity contour levels (g) CNLFM waveform AF frequency profile (h) SNLFM waveform AF frequency profile (i) CNLFM –LFM AF frequency profile difference (j) SNLFM-LFM waveform AF frequency profile difference (k) CNLFM waveform AF spectrogram (l) SNLFM waveform AF spectrogram for (Pulse width = 1 , Chirp rate = 4 1013)
radiometric motions. According to Fig. 1, in the case of straight flight path with broadside radiation, the slant range distance between the platform and a point target can be expressed as Eq. (7) [11]:
(7)
Where, v is the velocity of platform, R0 is slant range at the center of the aperture, is the squint angle measured between R0 and zero Doppler plane, θ is the equivalent squint angle, and is the point scatterer location over the target terrain. Hence, to investigate the performance of the proposed continuous CNLFM waveform, the SAR simulations are done based on system parameters in Table. I.
Table I. SAR system parameters
Parameter |
Value |
Unit |
Modulation type |
CNLFM |
- |
Carrier frequency |
15 |
GHz |
Radar transmitter bandwidth |
50 |
MHz |
Chirp rate |
5 1012 |
Hz/Sec |
Pulse width |
10 |
|
Azimuth sampling rate |
2 |
KHz |
Range sampling rate |
60 |
KHz |
Initial slant range |
9452.2 |
Meter |
Minimum slant range |
9152 |
Meter |
Maximum slant range |
10947 |
Meter |
Platform velocity |
800 |
(a) (b)
(c) (d)
Fig. 7. NLFM difference with LFM (a) real (CNLFM-LFM) (b) image (CNLFM-LFM) (c) real (SNLFM-LFM) (d) images (SNLFM-LFM) for (Pulse width = 10 , Chirp rate = 50 1012 and Bandwidth 50MHz)
The magnitude differences of SAR chirp waveform for both of the continuous CNLFM and the continuous SNLFM with the LFM waveform on the basis of Table. I are shown in Fig.7. The non-linear inherent characteristics of NLFM modulation is the reason of differences in Fig. 7 (a) and 7(c) and 7(b), 7(d).
Fig. 8. RDA IFA for the proposed missile-borne SAR
This section simulates the resulted images of received backscattered signal of single point target on the basis of LFM and NLFMs waveforms with the help of RDA [11]. All the simulations are done based on the presumption of similar supersonic flight trajectory with the velocity of 800(m/sec). According to the flight geometry in Fig.1 and with the help of Eq. (6), the modified RDA processing procedure is proposed as Fig.8. It should be noted that the changes were made to the RDA because of different input waveforms.
The RF input, , firstly down converted into an intermediate frequency (IF) signal, , and then directly transformed in to baseband, , with the help of Eq. (6). Fig.1, Fig. 9(a) and 9(b) depict the non-linear characteristic of the frequency modulation in both the slow time and the fast time of the generated raw data, , while continuous CNLFM and the SNLFM waveforms are being used.
Range compression is performed with a fast convolution when the data are in the azimuth time domain ( ). In the other words, a range FFT is performed ( ) on Eq. (6) followed by a range match filter multiply , and finally a range IFTT, to complete the range compression. This is done as equation (8):
(8)
Where, is the generated raw data, is the IFFT of the window , is the range FFT of , is the match filter function and is the range compressed spectral
(a)
(b)
Fig. 9. 3D-view of Generated raw data (a) continuous SNLFM waveform (b) continuous CNLFM waveform
(a)
(b)
Fig. 10. 2D-view of CNLFM raw data (a) before range compression (b) after range compression
equation. For the simplicity the value of is supposed to unit in order to compare the result of LFM, continuous CNLFM and SNLFM implementation in the final image. For a rectangular window or tapered window, is sinc function or a sinc-like function with lower sidelobes. Fig. 10 (a) and (b) show both the CNLFM raw data and CNLFM raw data after range compression.
An azimuth FFT transforms the data into the Range-Doppler domain ( ) as Eq. (9):
(9)
Where, is an azimuth FFT transform of range compressed data in the Doppler domain. The azimuth beam pattern is transformed into , with its shape preserved. The second phase term is the azimuth modulation which also has LFM and NLFM characteristic in , while the first exponential term carries the inherent phase information of the target. The value of is due to center frequency and center time, .
RCMC is the range time and frequency dependent and it is performed in the Range-Doppler domain. The amount of RCM to correct is given by ( ) which is a function of azimuth frequency and represents the target displacement as a function of azimuth frequency. Note that is also a function of that is the range variant. The total procedure of RCMC is:
(10)
Where, is the amount of RCM, is the linear phase term of RCM, is the corrected signal with RCMC and is the corrected signal in both range time and azimuth time.
Since the data after RCMC are in the Range-Doppler domain, it is convenient and efficient to implement the azimuth match filter in this domain as a function of slant range, , and azimuth frequency, . This step can be performed as the frequency domain match filter, , which is multiplied with the range gate as (11):
(11)
Where, is the frequency domain match filter, is the function of , is the azimuth compression after RCMC.
This step, transforms the data back to the time domain, resulting in the complex image as Eq. (12):
(12)
Where, is the amplitude of the azimuth impulse response, a sinc function similar to, . The first exponential shows the target phase due to its range position, , and the second is the linear phase term due to doppler center frequency. Accordingly, the final results of RDA using LFM, continuous SNLFM and the continuous CNLFM waveforms are shown in Fig.11 and Fig.12.
The values and the values in simulation results are azimuth and range sampling rate. According to spectral expansion simulation results in Fig. 11(b), 12(c) and 12(d) it is evident that waveform plays a key role in the pixel overlapping on the basis of non-linearity behavior of quadratic phase modulation which directly affects the reconstructed image. Similarly, the phase contour levels as shown in Fig. 11(c), 12(e) and 12(f) also have the same results while being analyzed on the basis of inherent non-linear modulation characteristics. The difference between Fig. 11(c), 12(e) and 12(f) is because of inherent differences of the AF contour levels which are existied in Fig. 5(b), 6(c) and 6(d). Likewise, frequency profiles are derived as Fig. 11(d), 12(g) and 12(h). It is obvious that Fig. 11(d) is quite different from 12(g) and 12(h), which is the result of the AF frequency profile in Fig. 5(c), 6(e) and 6(f). It should be noted that the reason of similarity between Fig. 12(g) and 12(h) is because of negligible difference in Fig. 5(c) and 6(e). Hence, the total procedure of the AF quality analysis of continuous NLFM waveforms shows that CNLFM waveform has more effects on the image contrast in comparison to the calculated frequency profile. On the other hand, in order to find the differences
between SNLFM and CNLFM, quantitative quality metrics would be a key which will be done on the RDA images in next section.
This section introduces and implements full-reference objective image quality assessment metrics on reference resulted LFM image and reconstructed NLFM images which are shown in Fig.11(a), 12(a), and 12(b).
The full-reference objective image quality metrics includes the SNR, PSNR, MSE, PSD, NCC and the ISLR. The SNR is calculated as the ratio of the mean value of the image and the standard deviation of the noise within the image. The higher the SNR value, the sharper the reconstructed NLFM images will be.
(a) (b)
(c) (d)
Fig. 11. final RDA image of single point scattrer with LFM waveform (a) 3D-view of RDA after azimuth compression (b) 2D-view of target spectral expansion (c) 2D-view of scatterer contour levels (d) received frequency profile of target
The PSNR is the ratio between the maximum possible power of the image and the power of corrupting noise that affects the quality of the image with the help of MSE as a scaled definition of power by the number of pixels in the images. There is an inverse relationship between PSNR and MSE. So, the higher the PSNR value indicates the higher quality of the reconstructed NLFM images which introduces the lower value of errors.
The PSD as a major quality assessment metric is also measures the strength of the energy as a function of frequency which shows at which frequencies variations are strong and at which frequencies variations are weak. The PSD is the best way to study the Non-linear modulation characteristics of the implemented NLFM waveforms in the SAR applications.
The NCC is used as a criterion to evaluate the degree of similarity between the LFM image as the reference waveform and the reconstructed CNLFM and SNLFM images which are proposed to replace the LFM. In addition to the qualitative comparison analysis based on the AF diagrams and aforementioned objective image quality metrics, ISLR is another important metric in SAR image quality assessment techniques. The ISLR measurement is done in a square area with 10 cells width and a center of maximum
amplitude. It should be noted that, the simulations of continuous NLFM waveforms and LFM are carried out in the noise-free target terrain situation. However, the noisy condition does not affect the entity of the objective image quality assessment process on how the simulation is running, but the result values would
(a) (b)
(c) (d)
(e) (f)
(g) (h)
Fig. 12. (a) CNLFM 3D-view of RDA after azimuth compression (b) SNLFM 3D-view of RDA after azimuth compression (c) CNLFM 2D-view of target spectral expansion (d) SNLFM 2D-view of target spectral expansion (e) CNLFM 2D-view of scatterer contour levels (f) SNLFM 2D-view of scatterer contour levels (g) CNLFM received processed frequency profile (h) SNLFM received processed frequency profile
be different. In this paper it is basically focused on the NLFM waveforms applications and their effects on the quality while, the noisy channel is not opened.
According to Section IV.A, the pixel difference-based measurements results are listed in Table.II, while, the input reference LFM image has an SNR of 52.14dB. On the other hand, about 0.8dB difference in the SNR value between CNLFM and SNLFM confirms that CNLFM affects the image quality better than SNLFM. It should be noted that a noise free target terrain usually has SNR values higher than 50dB and the simulations verify the process of objective quality metrics extraction..
The PSNR calculation also illustrates that both the CNLFM and SNLFM waveforms have the same value because they are testing in noise-free situation and 23.008dB is acceptable. The MSE value is close to zero in such medium. The PSNR and MSE values are listed in Table.III.
In the case of dealing with the similar PSNR value for objective quality assessment, structural similarity index map (SSIM) is suitable alternative to demonstrate the performance of non-linear waveforms implementation in comparison to LFM. According to the presence of the single point scatterer, the SSIM difference of CNLFM and SNLFM reconstructed images is presented in Fig.13. As shown in Fig.13, it is deduced that the CNLFM image reconstruction is more efficient in comparison to the SNLFM. By applying 2D discrete Fourier transform on the reconstructed CNLFM and SNLFM images, their PSD is calculated. The PSD difference between the reference LFM image and the reconstructed CNLFM and SNLFM images, with the PSD difference between CNLFM and SNLFM images are also presented in Fig.14. Results in Fig.14, verifies that CNLFM has better PSD which means that energy distribution as a function of frequency are more efficient in comparison to SNLFM waveforms at the same frequencies. According to Fig.15, the NCC simulation results of CNLFM, SNLFM, an NCC difference about 0.4dB verifies that CNLFM has better effects in the image reconstruction based on cosine modulation. According to, Table. I and results in Fig .11(b) and Fig. 12(a), (b), the ISLR values for the LFM, continuous SNLFM and the continuous CNLFM are calculated as Table. IV. It shows that the ISLR value as an objective quality metric which is used specifically for SAR applications, verifies that the CNLFM has better performance in comparison to other waveforms. Besides the AF quality analysis and objective image quality assessment techniques, the impacts of SAR system parameters alteration on the ISLR while CNLFM is being used will be investigated here. Fig.16 shows the ISLR alteration while the SAR CNLFM transmitter bandwidth is being changed. It shows that the increase in in bandwidth up to 3 times will cause degradation of about 0.5dB in the ISLR value. Hence, the continuous CNLFM waveform is also robust to bandwidth alteration, which is the main parameter in SAR system resolutions. The antenna squint angle alteration is another concept in the analysis of SAR image objective image quality metrics based on ISLR value. Fig.17 shows that the increase in the antenna squint angle more than 6 degrees is not effective and for both the CNLFM waveform parameters an RDA must be modified to be adopted in the new squint angle. This alteration in the squint angle shows that a boundary of 6 degrees (from 0 to 6) is the effective antenna pointing angle while continuous CNLFM waveform is being used. It also shows that the ISLR is alternating between -10.46 dB to -16dB, which depicts 5dB robustness to the squint angle alteration. As
Wave form |
LFM |
CNLFM |
SNLFM |
SNR [dB] |
52.14 |
51.27 |
50.40 |
Wave form |
PSNR |
MSE |
CNLFM |
23.008 dB |
0.01 |
SNLFM |
23.008 dB |
0.01 |
(a) (b)
Fig. 13. SSIM difference of CNLFM and SNLFM (a) SSIM amplitude difference (b) 2D SSIM difference
(a) (b)
(c) (d)
Fig. 14. PSD presentation (a) reconstructed CNLFM (b) reconstructed SNLFM (c) PSD difference between LFM , CNLFM and SNLFM (d) PSD difference between reconstructed CNLFM and SNLFM images
mentioned in Section.II, on the basis of AF contour levels, the NLFM waveforms affect the azimuth spectral expansion and weighing it. According to Fig.18, it is concluded that the CNLFM is no more robust against the squint alteration for θ= . So, not only the waveform systematic parameters must be changed, but also the RDA must be modified indeed. It should be noted that the increase in the squint angle while using linear or non-linear waveforms will always have a destructive effect on the ISLR value.
(a) (b)
(c)
Fig. 15. NCC presentation (a) NCC of reconstructed CNLFM (b) NCC of reconstructed SNLFM (c) NCC difference contour between reconstructed CNLFM and SNLFM images
Wave form |
ISLR |
LFM |
-10.395 dB |
CNLFM |
-10.467 dB |
SNLFM |
-10.272 dB |
Fig. 16. ISLR alteration while radar transmitter bandwidth is being changed
The ISLR sensitivity to the center frequency alteration while other system parameters remain unchanged is also shown in Fig.19. It shows the irregular deviation in the ISLR value when the center frequency alters. More precisely, the impact of this deviation in the image quality index can be investigated in equations (10), (11), (12) and image formation design procedures. In other words, the main steps of the IFA is directly depending on the center frequency and sampling rates. On the other hand, the flight geometry and radar system parameters mutually depending together, in which modifying one of them must be accompanied with the correction of the other one. Otherwise, the system will not respond correctly. So,
Fig. 17. Proposed CNLFM waveform image, ISLR alteration while antenna squint angle is being changed
(a)
(b)
Fig. 18. Proposed CNLFM spectral expansion while squint angle increased up to 14deg (a) 3D-view (b) 2D-view
Fig. 19. ISLR alteration while center frequency is being changed based on proposed CNLFM waveform
changing the center frequency without correcting the other system parameters and IFA steps will lead to such response.
Type of change |
Total effects |
Squint angle deviation |
Direct effects on the ISLR value |
Transmitter bandwidth |
Direct effects on the range resolution |
Synthetic aperture observation time |
Direct effects on the azimuth resolution |
Hence, the continuous CNLFM waveform performance in comparison to the continuous SNLFM is discovered more robust and introduces itself as a robust alternative to LFM counterpart. Table.V. summarizes the results of the SAR system parameters alteration in the ISLR index qualitatively.
This article compares the characteristics of LFM, continuous SNLFM and continuous CNLFM waveforms based on their AF quality analysis and subsequently, the optimum continuous CNLFM waveform is selected for further analysis. The impact of CNLFM implementation was investigated in comparison to LFM and SNLFM with the help of objective image quality assessment metrics analysis. All the simulations were carried out based on the similar flight trajectory and system parameters. On the other hand, complete objective image quality assessment metrics were derived and compared to each other. The SNR, PSNR, MSE, PSD, NCC and ISLR value were extracted and compared based on using different
waveforms. All the results validate that CNLFM waveform is not only a suitable alternative for LFM but also is more robust than SNLFM and even some cases to LFM. As verified above, the PSD of CNLFM in comparison to SNLFM and LFM make it more suitable for using in electronic warfare environments.
The ISLR comparison results show that the continuous CNLFM waveform not only reduces the level of the side-lobes but also it increases the main-lobe level up to 1.2 times and makes the IFA more robust against spectral expansion and Doppler deviations in the image. The continuous CNLFM waveform robustness procedure is evaluated by changing the antenna squint angle and bandwidth while the ISLR value of the image is measured. The results presented in tables, indicate that the increase in bandwidth up
to 3 times will have a degradation of 0.5 dB in the ISLR, while increasing the squint angle up to 6 degrees will reduce the ISLR value about 5dB.
Future activities in the context of this article might involve using new the waveforms, increasing the squint angle and modifying the image formation algorithm, using the method of time-frequency processing in the IFA processing.