Next Article in Journal
Spectral Mixture Analysis as a Unified Framework for the Remote Sensing of Evapotranspiration
Next Article in Special Issue
Geometric, Environmental and Hardware Error Sources of a Ground-Based Interferometric Real-Aperture FMCW Radar System
Previous Article in Journal
Climate Data Records from Meteosat First Generation Part I: Simulation of Accurate Top-of-Atmosphere Spectral Radiance over Pseudo-Invariant Calibration Sites for the Retrieval of the In-Flight Visible Spectral Response
Previous Article in Special Issue
Cross-Pol Transponder with Frequency Shifter for Bistatic Ground-Based Synthetic Aperture Radar
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Compressive Sensing for Ground Based Synthetic Aperture Radar

Department of Information Engineering, University of Florence, via Santa Marta, 3, 50139 Firenze, Italy
*
Author to whom correspondence should be addressed.
Remote Sens. 2018, 10(12), 1960; https://doi.org/10.3390/rs10121960
Submission received: 10 October 2018 / Revised: 30 November 2018 / Accepted: 2 December 2018 / Published: 5 December 2018

Abstract

:
Compressive sensing (CS) is a recent technique that promises to dramatically speed up the radar acquisition. Previous works have already tested CS for ground-based synthetic aperture radar (GBSAR) performing preliminary simulations or carrying out measurements in controlled environments. The aim of this article is a systematic study on the effective applicability of CS for GBSAR with data acquired in real scenarios: an urban environment (a seven-storey building), an open-pit mine, and a natural slope (a glacier in the Italian Alps). The authors tested the most popular sets of orthogonal functions (the so-called ‘basis’) and three different recovery methods (l1-minimization, l2-minimization, orthogonal pursuit matching). They found that Haar wavelets as orthogonal basis is a reasonable choice in most scenarios. Furthermore, they found that, for any tested basis and recovery method, the quality of images is very poor with less than 30% of data. They also found that the peak signal–noise ratio (PSNR) of the recovered images increases linearly of 2.4 dB for each 10% increase of data.

Graphical Abstract

1. Introduction

Ground-based synthetic aperture radar (GBSAR) systems are popular remote sensing instruments for short-term [1] and long-term [2] monitoring of landslides, glaciers [3], and mines [4] as well as for detecting small displacements of bridges [5] and dams [6].
The working principle of GBSAR is based on dense spatial sampling (steps smaller than a quarter of wavelength) along the synthetic aperture (usually a linear mechanical guide).
Compressive sensing (CS) is a recent sampling paradigm [7,8] which asserts one can recover certain signals from far fewer samples or measurements than traditional methods use. Its basic idea relies on the ‘sparsity’ of the signals of interest (the radar signals typically have this property [9]), and the incoherence of the sensing modality. The latter property is obtained through random sampling.
Generally speaking, CS can be applied in frequency domain (for reducing the number of transmitted frequencies) or in spatial domain (for reducing the number of spatial steps along the synthetic aperture).
In 2010, Huang et al. [10] applied CS to step-frequency through-wall radar that scanned along a horizontal mechanical guide. In 2011, Karlina and Sato [11] proposed to apply CS to GBSAR and they performed simulations, as well as Zonno in 2014 [12]. Yigit [13] extended CS to GB-SAR/ISAR acquisition modality and performed measurements in controlled environment. The interferometric properties of CS-GBSAR have been demonstrated by Giordano el al. [14] in 2015.
At the state of the art, CS techniques appear very promising for GBSAR applications, but all the cited works are relative to radar acquisitions in controlled environments (anechoic chamber or short-range experiments with controlled targets). Therefore, the crucial question is if these techniques are effectively applicable with data acquired in real scenarios, like GBSAR images of large buildings, landslides, open-pit mines, and glaciers.
Furthermore, a few specific questions about the CS methods have to be addressed. The first one is how to choose the orthogonal set of functions (the basis) that are used for recovering the signal. The possible bases are countless. The possible recovery algorithms are many as well. Therefore, a general selection criterion is necessary both for basis selection and recovery method. Another question is: what is an acceptable compression with real data. The claims declared in works that use simulated data or acquisitions performed in very cooperative scenarios are probably too optimistic.
Finally, we noticed that in previous papers [12,14] the ϕ matrix (see next chapter for his definition) is built by under-sampling the complete set of data acquired along a single scan. It means that each line could use a different set of positions along the scan. In this article, the ϕ matrix is built using the same reduced set of data for all the lines. This could reduce the incoherence of sampling (and so the quality of reconstruction), but it is a solution that could be easily implemented with a physical array or with a MIMO (multiple input multiple output), such as in [15,16,17].

2. Materials and Methods

Compressive sensing theory states a sparse signal can be recovered by fewer samples than what is required by the Nyquist theorem. With reference to Figure 1, a standard GBSAR transmits and receives at N positions along a linear mechanical guide. The Nyquist theorem requires that the spatial step has to be smaller than a quarter of wavelength (λ/4) for omnidirectional antennas (this constraint is a bit more relaxed for directional antennas, but it is not essential in the discussion that follows). Let E be the vector of the N samples acquired according the Nyquist theorem. A subset of M random positions is selected between the N positions. These M positions are the effective points where the radar-head performs the measurements that will be used for recovering the complete set of N samples.
The next step is to define M variables (ym) as random linear combinations of the N samples using only the M positions as sketched in Figure 2.
The vector y of the M variables can be written as
y = φ E
with ϕ matrix M × N as defined in Figure 2. It should be noted that the vector y is a measurable quantity, as it can be calculated from the NA measurements.
A radar signal is often sparse in a suitable set of functions. As an example, an ideal scenario with a single point-target in far field appears as a single peak in the domain of spatial frequencies. Inspired by this idea, we can write the vector E as
E = ψ b
where
ψ = I F F T ( [ 1 0 0 0 1 0 0 0 0 0 0 0 1 ] )
and where b is the vector of the coefficients of spatial Fast Fourier Transform (FFT) of the signal. If the signal is sparse, many coefficient bn are almost null.
The FFT or the DCT (discrete cosine transform) should be able to effectively represent ideal point targets, but they could be a too rough approximation for more realistic distributed targets. Wavelets transforms give better results with distributed targets [18]. The possible complete sets of wavelets are countless, so we have limited this study to the most popular sets [18,19,20,21] specifically: Haar wavelets, Daubechies wavelets, biorthogonal wavelet, coiflets, symlets, LeGall wavelets, and discrete Meyer wavelets.
By substituting Equaiton (2), in Equaiton (1), the vector y can be expressed as
y = ϑ b
with θ = ϕΨ. This matrix is the so-called ‘dictionary’. Equaiton (4) is a linear system with many more variables than equations. It is an ill-posed problem, but it can be solved with a suitable recovery method. The most popular are: l1-minimization [22], l2-minimization [23], and orthogonal matching pursuit [24]. Using one of these methods, the vector b is recovered, so the E vector of N samples is obtained by Equation (2).
The procedure described above has to be repeated for each frequency for obtaining the matrix Ek,i, with k-index ranging from 1 to Nf (number of frequencies) and i-index ranging from 1 to N.
The next step is to focus the matrix Ek,i using a back-propagation algorithm that takes into account of the phase history of each contribution relative to one specific position and one frequency.
With reference to Figure 3, the image value I(x,y) in a generic image point P(x,y) can be calculated as
I ( x , y ) = k = 1 N f i = 1 N E k , i e 1 4 π c f k R i ( x , y )
where c is the speed of light. This formula is correct but computationally very inefficient. More efficient algorithms make use of FFT and interpolation (see for example: [25]).
Since I(x,y) is a complex number, it provides the phase information too. This can be exploited for generating differential interferograms. Displacement maps can be obtained from differential interferograms using the well-known relationship [5]
Δ r ( x , y ) = λ 4 π Δ φ ( x , y )
where Δr(x,y) is the displacement in the point P(x,y), Δϕ(x,y) is the differential phase in the point P(x,y), and λ is the wavelength at the central frequency.

3. Results

3.1. Simulations

The CS algorithms have been preliminary tested with the simulation of a single target at 150 m distance in front of radar. The used parameters were: initial frequency f1 = 9.915 GHz, final frequency f2 = 10.075 GHz, number of frequencies Nf = 801, length of the mechanical scan L = 1 m, number of samples (according to Nyquist theorem) N = 100.
The quality of the recovered images were estimated with the peak signal–noise ratio (PSNR) calculated as
P S N R = 10 log 10 ( max ( | I i , j | 2 ) 1 m n i = 1 m j = 1 n | I i , j R i , j | 2 )
where Ri,j is the reference GBSAR image obtained from the data sampled according to Nyquist theorem without application of CS techniques; Ii,j is the image obtained using the CS techniques. Both images are complex matrix m × n. As rule of thumb, images with PSNR lower than 20–25 dB appear of low quality, with features hardly recognizable [26].
The used basis was the Haar wavelets, the recovery method was l2–minimization (this choice of Haar wavelets will be motivated in the experimental section, l2–minimization has been selected because it works even with a number of variable smaller than the number of equations, nevertheless the choice of both basis and recovery method is not essential for the discussion that follows). The PSNR has been calculated for an image 200 × 200 in polar coordinates. The range limits were 100 m and 200 m. The azimuth limits were −30 deg and +30 deg.
Figure 4 shows the plots of PSNR versus M/N. Each PSNR was calculated as the average of 25 random iterations and the standard deviation of each value is reported in the plot (blue line). It is interesting to note that for M/N > 0.3, PSNR increases linearly of 2.4 dB for each 10% increase of M/N. This is a notable result that we found also true for simulations of more complex scenarios (several targets), using different measurement parameters, different dictionaries, and different recovery methods, and in experimental data too. Note that PSNR does not tend to infinity for M = N, but it stops at about 37dB. The reason of it is related to the limited size of the φ matrix. As it is known, the CS reconstruction relies on the incoherence of both random sensing and the random coefficients of the matrix φ. When M = N, the sensing is no longer incoherent, but the incoherence of the ϕ coefficients persists and limits the quality of reconstruction. The φ coefficients are uniformly distributed random numbers between 0 and 1. When M = N, the number of elements is N × N = 100 × 100 = 10,000, the sum of 10,000 uniformly distributed random numbers between 0 and 1 gives exactly 37.0 dB, that is just the maximum PSNR obtained in Figure 4.
The noise can affect both PSNR and displacement retrieved by interferometry. In order to verify/evaluate it, we have repeated the simulations reported above with controlled Gaussian noise. Figure 5 shows the PSNR versus M/N calculated with different SNR values (the SNR reported in the plots is relative to one single measurement specified by one single frequency and one single position along the mechanical guide). Each PSNR value of the plots in Figure 6 is calculated as the average of 12 random iterations. As expected, by reducing the SNR the PSNR decreases. When the SNR of a single measurement is >0, the PSNR is practically equal to the PSNR without noise calculated in the previous chapter.
As mentioned above, the noise can give an error in the displacement retrieved by interferometry. In order to verify/evaluate it, we have estimated the phase error (Δϕ) and we have converted it in displacement error (ΔR) using the well-known relationship [5]:
Δ R = λ 4 π Δ φ
Figure 6 shows the estimated error in function of SNR. Each PSNR value is calculated as the average of 40 random iterations. For NA/N = 1 the error plot in function of SNR overlaps perfectly the plot of the error calculated without CS that, in log scale, is a straight line with slope −1. As expected, for M/N = 0.5 the error increases of 3dB and for M/N = 0.1 the error increases of 10 dB.
When M/N increases, the number of samples increases too. So the error should decrease. This intuitive idea has been verified with a simulation. Figure 7 shows how the error decreases increasing M/N, for two different values of the SNR. Both the plots are approximately straight lines with slope −1.
Generally speaking, the interferometric error could not depend only on Gaussian noise (i.e., on thermal noise). In that case, calculating the error in function of M/N does not make much sense, as it could be hard to obtain a straight line. On the contrary, the differential error (i.e., the error with respect to the image obtained using the CS with 100% of the samples) is expected to have a linear behavior (in log scale). Figure 8 shows the differential displacement in function of M/N for three different SNR. It is evident that the plots can be approximated as straight lines with a slope of −2.

3.2. Experimental Test Site “Campus”

The CS techniques have been tested in a controlled experimental test site at the University of Florence (named “Campus”). A single metallic target, provided with a micrometric linear translator, has been positioned in front of the radar equipment at 150 m distance in a large flat garden, as shown in the pictures in Figure 9.
The measurement parameters were the same used in the simulations. The acquired experimental data have been used for testing the most popular bases and recovery methods. Specifically we used discrete cosine transform (dct), fast Fourier transform (fft), the Haar wavelts (haar), the Daubechies wavelets from order 2 to order 10 (db2–db10), seven biorthogonal wavelets (bior1.3, bior2.2, bior2.6, bior3.1, bior3.5, bior3.9, bior5.5), the “coiflets” from order 1 to 5 (coif1–coif5), the “simlets” from order 2 to 8 (sim2–sim8), LeGall wavelets (legal5.3), the discrete Meyer wavelets (dmey). As recovery methods, we used l1-minimization (L1), l2-minimization (L2), and orthogonal matching pursuit (OMP).
Table 1 reports the obtained PSNR for M/N = 0.5 and M/N = 0.33. The PSNR has been calculated for an image 200 × 200 in polar coordinates. The range limits were 100 m and 200 m. The azimuth limits were −30 deg and +30 deg.
In the table, we have highlighted the highest (in bold) and the lowest (in red bold) values. Using 50% of data, PSNR varies between 26.01 dB (Haar wavelets, OMP) and 21.89 dB (FFT, OMP). Using 30% of data, PSNR varies between 20.03 dB (Haar wavelets, OMP) and 15.87 dB (Bior3.1 wavelets, OMP).
In simulated data, we have noted that PSNR increases linearly of 2.4 dB for each 10% increase of M/N. This relationship is perfectly confirmed also for these experimental data as shown in Figure 10.
For evaluating the interferometric error we acquired 100 measurements in 96 h. By considering one pixel associated to a stable point of one building at about 200 m in the background, we evaluated the standard deviation of the displacement retrieved by interferometry. Figure 11 shows the obtained values varying M/N. The error is rather constant. It means that thermal (Gaussian) noise is not the main source of error. The effective contribution of thermal noise is more evident in the plot of the differential error (Figure 12) that confirms the error (in log scale) linearly decreases with slope −2, increasing the M/N (in log scale), as predicted by simulation.
The metallic target in the middle of the grass was provided with a micrometric positioner able to translate it with a nominal accuracy of 0.1 mm. So, we performed 20 radar acquisitions by moving the target of 1 mm forward the radar after each acquisition. Figure 13 shows the cumulative displacement detected by interferometry using 100% of data without CS, using 100% of data with CS, using 50% of data with CS, and using 30% of data with CS. The four plots are practically overlapping. Figure 14 shows the plots of the differences. In the worst case (CS with 30% of data) the error is lower than 25 μm.

3.3. A Seven-Storey Building

In order to test the application of CS techniques in an operative scenario representative of an urban environment, we performed radar measurements using as target the seven-storey building shown in Figure 15.
The measurement parameters were: initial frequency f1 = 9.915 GHz, final frequency f2 = 10.075 GHz, number of frequencies Nf = 801, length of the mechanical scan L = 1.8 m, number of samples (according to Nyquist theorem) N = 180. The obtained radar image (without application of CS techniques) is shown in Figure 16.
As well as for simulated data, and for the experimental data acquired at the “Campus” test site, we verified that PSNR increases linearly of 2.4 dB for each 10% increase of M/N and this relationship has been perfectly confirmed also for these experimental data (we do not report the graph for sake of brevity).
Table 2 reports the obtained PSNR for M/N = 0.5 and M/N = 0.3. We have highlighted the highest (in bold) and the lowest (in red bold) values. Using 50% of data, PSNR varies from 22.78 dB (db2 wavelets, OMP) to 8.19 dB (dct, OMP). Using 30% of data, the PSNR ranges between 19.12 dB (db5 wavelets or sym2 wavelets, OMP) and 1.80 dB (fft, OMP).

3.4. An Open-Pit Copper Mine

The most popular use of the GBSAR systems is the monitoring of open-pit mines. Therefore, we tested the application of the CS techniques to a representative case study of a copper mine in South America. The measurement parameters were: initial frequency f1 = 17.1 GHz, final frequency f2 = 17.3 GHz, number of frequencies Nf = 5333, length of the mechanical scan L = 1.275 m, number of samples (according to Nyquist theorem) N = 256. The obtained radar image (without application of CS techniques) is shown in Figure 17.
As well as for previous data, we verified that PSNR increases linearly of 2.4 dB for each 10% increase of M/N and this relationship has been perfectly confirmed also for these experimental data (we do not report the graph for sake of brevity).
Table 3 reports the obtained PSNR for M/N = 0.5 and M/N = 0.3. Using 50% of data, PSNR varies from 29.95 dB (Haar wavelets, OMP) to 17.25 dB (FFT, L2). Using 30%, the PSNR ranges between 29.44 dB (Haar wavelets, OMP) and 9.33 dB (FFT, OMP).
Figure 18 shows an interferogram between two images (focused without application of CS techniques) taken at different time. It shows a clear movement on the left side of the slope. In order to test the capability of the CS techniques to reconstruct interferograms, we applied them to 50% of data (Figure 19) and to 33% of data (Figure 20) using Haar wavelets and OMP.
In order to evaluate the interferometric error, we have processed 41 measurements acquired in 4 h and 20 min. By considering one pixel associated to one stable point (without visible fringes), we have evaluated the standard deviation of the displacement retrieved by interferometry. The point is labeled with letter A in Figure 17. Figure 21 shows the error (calculated as standard deviation) of the displacement at the stable point A varying M/N. The error is rather constant. It means that thermal (Gaussian) noise is not the main source of error. The effective contribution of thermal noise is more evident in the plot of the differential error (Figure 22). It confirms that the error linearly decreases (in log scale) with slope −2, increasing the M/N (in log scale), as predicted by simulation.
In order to verify the capability of CS to detect a cumulative displacement, we have taken into account the point B in Figure 18 where several fringes are evident in the interferograms (see Figure 19). Figure 23 shows the plot obtained without CS and the plot obtained with CS applied to 100% of data are perfectly overlapping. The CS plots with 50% and 30% of data tend to underestimate the displacement.

3.5. A Glacier

As test-site representative of a natural scenario we considered the “Belvedere” glacier in the Italian Alps. The data was acquired during a measurement campaign in 2006, whose details are reported in [3,27,28]. The measurement parameters were: f1 = 5.97 GHz, f2 = 5.99 GHz, Nf = 801, L = 1.71 m, N = 141. The obtained radar image (without application of CS techniques) is shown in Figure 24.
As well as for previous data, we verified that PSNR increases linearly of 2.4 dB for each 10% increase of M/N and this relationship has been perfectly confirmed also for these experimental data (we do not report the graph for sake of brevity).
Table 4 reports the obtained PSNR for M/N = 0.5 and M/N = 0.3. The PSNR has been calculated for an image 200 × 200 in polar coordinates. The range limits were 350 m and 1000 m. The azimuth limits were −25 deg and +25 deg. Using 50% of data, PSNR varies from 21.16 dB (Haar wavelets, OMP) to −3.50 dB (DCT, L2). Using 30%, the PSNR ranges between 16.52 dB (Haar wavelets, OMP) and −5.09 dB (FFT, L1).
In order to verify the capability of CS to detect a cumulative displacement, we have taken into account a pixel inside the sliding body of the glacier. Figure 25 shows the plots obtained without CS and with CS applied to 100, 50, and 30% of data. The plots are almost perfectly overlapping.
Table 5 resumes the parameters of each measurement campaign and the highest value of PSNR obtained using 50% of data.

4. Conclusions

Authors found that CS techniques are able to reconstruct GBSAR images of good quality (PSNR > 20–30 dB) with almost 30–40% of data. Furthermore, they found that for more than 30% of data, the PSNR of the recovered image increases linearly of 2.4 dB for each 10% increase of data for all the tested bases and recovery methods. Nevertheless, the quality of reconstruction depends on the selected basis and recovery method. In a test site with just one single strong target, the best image has been obtained with FFT and OMP. In an urban scenario the best choice has been Daubechies wavelets and OMP. For an open-pit mine the highest quality image has been reconstructed with Haar wavelets and OMP. For a natural slope (a glacier) the optimal choice has been again Haar wavelets and OMP. In all cases, Haar wavelets with OMP gave good results (even if not always the best). Therefore, we can state that Haar wavelets with OMP is a reasonable choice for most practical cases. Furthermore, we found that CS reconstruction does not worsen significantly the interferometric error. The plots of the cumulative displacement recovered with 100, 50, and 30% of data overlaps almost completely the plot obtained without application of CS techniques. All these findings constitute the basic scientific starting point for designing optimal sparse array or multiple-input multiple-output (MIMO) radar able to exploit the CS techniques.

Author Contributions

M.P. conceived and designed the experiments; L.M. performed the experiments; N.R. analyzed the data; M.P. wrote the paper.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Luzi, G.; Pieraccini, M.; Mecatti, D.; Noferini, L.; Macaluso, G.; Galgaro, A.; Atzeni, C. Advances in groundbased microwave interferometry for landslide survey: A case study. Int. J. Remote Sens. 2006, 27, 2331–2350. [Google Scholar] [CrossRef]
  2. Noferini, L.; Pieraccini, T.T.M.; Mecatti, D.; Macaluso, G.; Luzi, G.; Atzeni, C. Analysis of ground-based SAR data with diverse temporal baselines. IEEE Trans. Geosci. Remote Sens. 2008, 46, 1614–1623. [Google Scholar] [CrossRef]
  3. Noferini, L.; Mecatti, D.; Macaluso, G.; Pieraccini, M.; Atzeni, C. Monitoring of Belvedere Glacier using a wide angle GB-SAR interferometer. J. Appl. Geophys. 2009, 68, 289–293. [Google Scholar] [CrossRef]
  4. Liu, B.; Ge, D.; Li, M.; Zhang, L.; Wang, Y.; Zhang, X. Using GB-SAR technique to monitor displacement of open pit slope. In Proceedings of the International Geoscience and Remote Sensing Symposium (IGARSS), Beijing, China, 10–15 July 2016; pp. 5986–5989. [Google Scholar]
  5. Pieraccini, M.; Tarchi, D.; Rudolf, H.; Leva, D.; Luzi, G.; Bartoli, G.; Atzeni, C. Structural static testing by interferometric synthetic radar. NDT E Int. 2000, 33, 565–570. [Google Scholar] [CrossRef]
  6. Jenkins, W.; Rosenblad, B.; Gomez, F.; Le-garsky, J.; Loehr, E. Deformation measurements of earth dams using a ground based inter-ferometric radar. In Proceedings of the 2012 ASDSO Annual Conference on Dam Safety, Denver, CO, USA, 16–20 September 2012. [Google Scholar]
  7. Candès, E.J.; Wakin, M.B. An introduction to compressive sampling. IEEE Signal Process. Mag. 2008, 25, 21–30. [Google Scholar] [CrossRef]
  8. Massa, A.; Rocca, P.; Oliveri, G. Compressive sensing in electromagnetics-a review. IEEE Antennas Propag. Mag. 2015, 57, 224–238. [Google Scholar] [CrossRef]
  9. Hadi, M.A.; Alshebeili, S.; Jamil, K.; El-Samie, F.E.A. Compressive sensing applied to radar systems: An overview. Signal Image Video Process. 2015, 9, 25–39. [Google Scholar] [CrossRef]
  10. Huang, Q.; Qu, L.; Wu, B.; Fang, G. UWB through-wall imaging based on compressive sensing. IEEE Trans. Geosci. Remote Sens. 2010, 48, 1408–1415. [Google Scholar] [CrossRef]
  11. Karlina, R.; Sato, M. Compressive sensing applied to imaging by ground-based polarimetric SAR. In Proceedings of the 2011 IEEE International Geoscience and Remote Sensing Symposium (IGARSS 2011), Vancouver, BC, Canada, 24–29 July 2011. [Google Scholar]
  12. Zonno, M. GBSAR data focusing based on compressive sensing. In Proceedings of the 10th European Conference on Synthetic Aperture Radar EUSAR 2014, Berlin, Germany, 3–5 June 2014. [Google Scholar]
  13. Yiğit, E. Compressed sensing for millimeter-wave ground based SAR/ISAR imaging. J. Infrared Millim. Terahertz Waves 2014, 35, 932–948. [Google Scholar] [CrossRef]
  14. Giordano, R.; Guccione, P.; Cifarelli, G.; Mascolo, L.; Nico, G. Focusing SAR images by compressive sensing: Study of interferometric properties. In Proceedings of the International Geoscience and Remote Sensing Symposium (IGARSS), Milan, Italy, 26–31 July 2015; pp. 5352–5355. [Google Scholar]
  15. Feng, W.; Yi, L.; Sato, M. Near range radar imaging based on block sparsity and cross-correlation fusion algorithm. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 2079–2089. [Google Scholar] [CrossRef]
  16. Tarchi, D.; Oliveri, F.; Sammartino, P.F. MIMO radar and ground-based SAR imaging systems: Equivalent approaches for remote sensing. IEEE Trans. Geosci. Remote Sens. 2013, 51, 425–435. [Google Scholar] [CrossRef]
  17. Hu, C.; Wang, J.; Tian, W.; Zeng, T.; Wang, R. Design and Imaging of Ground-Based Multiple-Input Multiple-Output Synthetic Aperture Radar (MIMO SAR) with Non-Collinear Arrays. Sensors 2017, 17, 598. [Google Scholar] [CrossRef] [PubMed]
  18. He, L.; Carin, L. Exploiting structure in wavelet-based Bayesian compressive sensing. IEEE Trans. Signal Process. 2009, 57, 3488–3497. [Google Scholar]
  19. Chui, C.K. An Introduction to Wavelets; Academic Press: San Diego, CA, USA, 1992. [Google Scholar]
  20. Walker, J.S. A Primer on Wavelets and Their Scientific Applications; CRC Press: Boca Raton, FL, USA, 1999. [Google Scholar]
  21. Wickerhauser, M.V. Adapted Wavelet Analysis from Theory to Software; IEEE Press: Piscataway, NJ, USA, 1994. [Google Scholar]
  22. Yang, A.Y.; Sastry, S.S.; Ganesh, A.; Ma, Y. Fast ℓ1-minimization algorithms and an application in robust face recognition: A review. In Proceedings of the 17th IEEE International Conference on Image Processing (ICIP), Hong Kong, China, 26–29 September 2010; pp. 1849–1852. [Google Scholar]
  23. Baraniuk, R.G. Compressive sensing [lecture notes]. IEEE Signal Process. Mag. 2007, 24, 118–121. [Google Scholar] [CrossRef]
  24. Tropp, J.A.; Gilbert, A.C. Signal recovery from random measurements via orthogonal matching pursuit. IEEE Trans. Inf. Theory 2007, 53, 4655–4666. [Google Scholar] [CrossRef]
  25. Pieraccini, M.; Miccinesi, L. ArcSAR: Theory, simulations, and experimental verification. IEEE Trans. Microw. Theory Tech. 2017, 65, 293–301. [Google Scholar] [CrossRef]
  26. Salomon, D. Data Compression: The Complete Reference; Springer: New York, NY, USA, 2004. [Google Scholar]
  27. Pieraccini, M.; Noferini, L.; Mecatti, D.; Macaluso, G.; Luzi, G.; Atzeni, C. Digital elevation models by a GBSAR interferometer for monitoring glaciers: The case study of Belvedere Glacier. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, IGARSS 2008, Boston, MA, USA, 7–11 July 2008; Volume 4, p. IV-1061. [Google Scholar]
  28. Pieraccini, M.; Luzi, G.; Mecatti, D.; Noferini, L.; Atzeni, C. Ground-based SAR for short and long term monitoring of unstable slopes. In Proceedings of the 3rd European Radar Conference, EuRAD 2006, Manchester, UK, 13–15 September 2006; pp. 92–95. [Google Scholar]
Figure 1. Sampling along the linear mechanical guide of a GBSAR.
Figure 1. Sampling along the linear mechanical guide of a GBSAR.
Remotesensing 10 01960 g001
Figure 2. Definition of the M variables ym. The number different from zero in the matrix are selected randomly between 0 and 1.
Figure 2. Definition of the M variables ym. The number different from zero in the matrix are selected randomly between 0 and 1.
Remotesensing 10 01960 g002
Figure 3. Focusing geometry.
Figure 3. Focusing geometry.
Remotesensing 10 01960 g003
Figure 4. PSNR in function of M/N (simulated data).
Figure 4. PSNR in function of M/N (simulated data).
Remotesensing 10 01960 g004
Figure 5. PSNR in function of M/N and SNR (simulated data).
Figure 5. PSNR in function of M/N and SNR (simulated data).
Remotesensing 10 01960 g005
Figure 6. Displacement error in function of SNR for different M/N (simulated data). The dark line is the displacement error calculated without CS technique.
Figure 6. Displacement error in function of SNR for different M/N (simulated data). The dark line is the displacement error calculated without CS technique.
Remotesensing 10 01960 g006
Figure 7. Displacement error in function of M/N for different SNR (simulated data).
Figure 7. Displacement error in function of M/N for different SNR (simulated data).
Remotesensing 10 01960 g007
Figure 8. Differential displacement error in function of M/N for different SNR (simulated data).
Figure 8. Differential displacement error in function of M/N for different SNR (simulated data).
Remotesensing 10 01960 g008
Figure 9. Experimental test site “Campus”.
Figure 9. Experimental test site “Campus”.
Remotesensing 10 01960 g009
Figure 10. PSNR in function of M/N (data acquired at the experimental test site “Campus”).
Figure 10. PSNR in function of M/N (data acquired at the experimental test site “Campus”).
Remotesensing 10 01960 g010
Figure 11. Displacement error in function of M/N for a stable point in the test site “Campus”.
Figure 11. Displacement error in function of M/N for a stable point in the test site “Campus”.
Remotesensing 10 01960 g011
Figure 12. Differential displacement error in function of M/N for a stable point in the test site “Campus”.
Figure 12. Differential displacement error in function of M/N for a stable point in the test site “Campus”.
Remotesensing 10 01960 g012
Figure 13. Cumulative measured displacement of the controlled metallic target at the test site “Campus”.
Figure 13. Cumulative measured displacement of the controlled metallic target at the test site “Campus”.
Remotesensing 10 01960 g013
Figure 14. Differences between displacements detect through 100% CS, 50% CS, 30% and displacements detected using 100% of data without application of CS techniques (test site “Campus”).
Figure 14. Differences between displacements detect through 100% CS, 50% CS, 30% and displacements detected using 100% of data without application of CS techniques (test site “Campus”).
Remotesensing 10 01960 g014
Figure 15. Aerial picture of the seven-storey building used as urban test site.
Figure 15. Aerial picture of the seven-storey building used as urban test site.
Remotesensing 10 01960 g015
Figure 16. Radar image (without application of CS techniques) of the building shown in Figure 15.
Figure 16. Radar image (without application of CS techniques) of the building shown in Figure 15.
Remotesensing 10 01960 g016
Figure 17. Radar image (without application of CS techniques) of a copper mine in South America.
Figure 17. Radar image (without application of CS techniques) of a copper mine in South America.
Remotesensing 10 01960 g017
Figure 18. Radar image (without application of CS techniques) of a copper mine in South America.
Figure 18. Radar image (without application of CS techniques) of a copper mine in South America.
Remotesensing 10 01960 g018
Figure 19. Radar image obtained with 50% of data and CS.
Figure 19. Radar image obtained with 50% of data and CS.
Remotesensing 10 01960 g019
Figure 20. Radar image obtained with 33% of data and CS.
Figure 20. Radar image obtained with 33% of data and CS.
Remotesensing 10 01960 g020
Figure 21. Displacement error in function of M/N for a stable point in the open pit mine.
Figure 21. Displacement error in function of M/N for a stable point in the open pit mine.
Remotesensing 10 01960 g021
Figure 22. Differential displacement error in function of M/N for a stable point in the open pit mine.
Figure 22. Differential displacement error in function of M/N for a stable point in the open pit mine.
Remotesensing 10 01960 g022
Figure 23. Cumulative measured displacement of the point B in the open pit mine.
Figure 23. Cumulative measured displacement of the point B in the open pit mine.
Remotesensing 10 01960 g023
Figure 24. Radar image (without application of CS techniques) of a glacier in the Italian Alps.
Figure 24. Radar image (without application of CS techniques) of a glacier in the Italian Alps.
Remotesensing 10 01960 g024
Figure 25. Cumulative measured displacement of a pixel inside the sliding body the glacier.
Figure 25. Cumulative measured displacement of a pixel inside the sliding body the glacier.
Remotesensing 10 01960 g025
Table 1. PSNR calculated with different bases and recovery methods for experimental data acquired at the test site “Campus”.
Table 1. PSNR calculated with different bases and recovery methods for experimental data acquired at the test site “Campus”.
50%33%
L1 (dB)L2 (dB)OMP (dB)L1 (dB)L2 (dB)OMP (dB)
Dct22.9222.8922.2918.7719.0219.62
Fft22.6122.7321.8917.0017.1617.45
Haar24.8324.8326.0118.3018.3020.03
Db224.3224.2924.3218.1418.1318.13
Db324.3224.3324.3617.6817.6618.01
Db423.9723.9723.6718.2218.2317.98
Db523.7123.7323.1017.4917.4817.35
Db623.8723.9022.7417.7617.7717.98
Db723.2423.2522.6017.4117.4516.11
Db823.2423.2322.4817.7117.7017.82
Db923.2623.2721.4617.7117.7217.82
Db1023.0423.0723.4817.7417.7418.52
Bior1.324.6524.6425.0018.1618.1619.51
Bior2.224.4524.4523.5317.0016.9916.96
Bior2.624.2224.2123.3717.6517.6517.37
Bior3.123.5123.5122.6617.3917.3915.87
Bior3.524.2924.2823.5017.9917.9917.21
Bior3.924.0024.0423.6717.6017.6017.67
Bior5.523.9723.9823.6218.0418.0718.23
Coif124.0424.0923.5217.7117.7217.25
Coif224.2024.2124.0416.8416.8417.26
Coif323.7223.7223.2717.3917.3917.36
Coif423.4723.4722.7417.4917.5117.52
Coif523.2823.1822.7917.0517.0416.90
Sym224.3224.2924.3218.1418.1318.13
Sym324.3224.3324.3617.6817.6618.01
Sym424.4224.4424.0817.2917.2817.04
Sym523.7723.7723.6217.6317.6417.54
Sym624.2424.1523.7817.6417.6416.94
Sym724.2624.2724.1818.0218.0317.68
Sym824.1024.1123.5317.9417.9417.04
Legall5.324.1724.1724.4918.3118.3018.51
Dmey23.6023.6023.5617.9317.9417.96
Table 2. PSNR calculated with different bases and recovery methods for experimental data acquired at the seven-storey building.
Table 2. PSNR calculated with different bases and recovery methods for experimental data acquired at the seven-storey building.
50%33%
L1 (dB)L2 (dB)OMP (dB)L1 (dB)L2 (dB)OMP (dB)
Dct16.4716.598.1611.0410.855.93
Fft15.6215.3010.729.339.981.80
Haar22.7322.7321.5919.0019.0018.42
Db221.5021.5222.7818.3318.3218.85
Db321.6021.6222.7218.1318.1318.86
Db421.5521.5820.9917.8417.8317.85
Db521.0121.0221.6017.9117.9319.12
Db620.0920.9021.3917.9117.9118.98
Db720.4820.5220.6717.7917.7518.49
Db819.1619.1719.1514.7014.7016.17
Db919.1119.1417.8117.0517.0416.32
Db1019.8219.8319.2717.0317.0517.69
Bior1.322.2622.2721.4918.7818.7917.58
Bior2.221.9121.9222.4218.1318.1318.27
Bior2.621.3121.3222.3617.6717.6818.58
Bior3.121.2721.3120.5620.1020.1019.16
Bior3.520.1220.1219.6317.4917.5417.33
Bior3.920.7920.8019.4117.7717.8016.45
Bior5.522.1322.1422.7617.8017.8319.11
Coif121.7621.7621.7418.3318.3417.96
Coif221.9821.9621.9017.8517.8517.86
Coif321.7721.7721.7718.2918.3018.74
Coif420.2120.2121.8817.8717.8918.89
Coif521.3821.3821.6718.0018.0018.66
Sym221.5021.5222.7818.3318.3218.84
Sym321.6121.6122.7218.1318.1318.86
Sym421.4521.4421.1818.1518.2518.00
Sym521.3621.3820.1817.8417.8118.76
Sym621.1321.1220.9617.4717.4418.04
Sym721.7921.8021.5217.2217.2218.34
Sym821.4721.4720.0416.9016.9317.86
Legall5.322.0322.0422.6418.0618.0618.47
Dmey20.5520.5919.1617.4617.4616.95
Table 3. PSNR calculated with different bases and recovery methods for experimental data acquired at a copper mine in South America.
Table 3. PSNR calculated with different bases and recovery methods for experimental data acquired at a copper mine in South America.
50%33%
L1 (dB)L2 (dB)OMP (dB)L1 (dB)L2 (dB)OMP (dB)
Dct17.5418.2317.3615.8416.3115.90
Fft17.4117.2515.2611.6912.4410.05
Haar29.3929.3929.9525.8125.8129.44
Db228.2928.3028.2526.2926.3027.69
Db325.3425.3727.8523.9423.9426.22
Db424.5424.5323.1422.4822.4822.94
Db525.6925.7023.6821.3821.3722.74
Db625.2425.2426.3022.2622.2622.49
Db725.8825.8924.1421.4821.4923.01
Db824.0824.1021.7020.8220.8219.29
Db924.0424.0519.9821.2821.2916.25
Db1023.6823.6818.8222.2522.2618.27
Bior1.328.7928.7927.4525.4225.4227.47
Bior2.227.6027.5926.6322.7322.7218.68
Bior2.626.4126.4126.2822.6222.6019.09
Bior3.126.4326.4422.5624.3824.3814.96
Bior3.525.5425.5521.4122.8522.8617.30
Bior3.925.1825.1920.4922.7023.1718.03
Bior5.525.2425.2427.2623.3323.3526.33
Coif127.1327.1426.5821.3221.3320.95
Coif224.8324.8424.0120.7420.7521.29
Coif325.0525.0623.7322.3722.3721.18
Coif426.1026.0925.5622.9722.9722.85
Coif526.0526.0525.9422.8022.8025.93
Sym228.2928.3028.2526.2926.3027.69
Sym325.3425.3527.8523.9423.9426.22
Sym425.4725.4625.6321.9421.9419.47
Sym527.6227.6325.5423.0123.0624.35
Sym625.5325.5425.0822.4522.4518.56
Sym724.9524.9626.1024.4824.4825.98
Sym825.5025.5424.2622.1822.2118.73
Legall5.328.7828.7829.3926.7826.7827.58
Dmey25.1525.1621.7222.3122.3118.04
Table 4. PSNR calculated with different bases and recovery methods for experimental data acquired at a glacier in the Italian Alps.
Table 4. PSNR calculated with different bases and recovery methods for experimental data acquired at a glacier in the Italian Alps.
50%33%
L1 (dB)L2 (dB)OMP (dB)L1 (dB)L2 (dB)OMP (dB)
Dct−2.26−3.50−2.771.341.390.21
Fft−1.671.85−0.92−5.09−4.93−4.06
Haar20.5320.5321.1615.0115.0216.52
Db213.5313.5319.278.218.2614.09
Db39.549.5418.835.645.7514.02
Db44.054.1013.058.758.758.56
Db55.565.6911.147.177.1910.53
Db62.892.8717.513.793.8711.77
Db711.2611.2410.678.688.675.58
Db82.622.648.831.781.770.92
Db96.696.7410.255.535.52−3.23
Db102.322.333.29−2.96−2.902.16
Bior1.313.9613.9510.2711.4511.418.57
Bior2.214.3714.3919.179.939.9213.99
Bior2.69.849.8618.156.917.0113.03
Bior3.16.656.645.167.507.468.46
Bior3.54.814.812.937.217.2110.65
Bior3.93.313.27-3.016.086.537.56
Bior5.510.5510.5214.624.114.0914.02
Coif110.9811.1119.7210.5810.6114.27
Coif212.1712.1619.667.987.9914.55
Coif39.489.4818.726.566.5610.16
Coif411.6511.6514.465.765.7511.57
Coif55.875.8811.369.499.493.42
Sym213.5313.5319.288.218.2614.09
Sym39.549.5418.835.645.7514.02
Sym411.1411.1917.3910.1010.1014.03
Sym59.529.519.784.834.513.63
Sym612.4812.4912.947.377.377.37
Sym79.359.3716.862.442.4511.40
Sym813.6413.5912.177.007.0011.34
Legall5.311.7611.7915.759.439.4014.62
Dmey6.836.8410.725.605.6111.78
Table 5. Parameters of each measurement campaign and the highest value of PSNR obtained using 50% of data.
Table 5. Parameters of each measurement campaign and the highest value of PSNR obtained using 50% of data.
CampusBuildingMineGlacier
f1 (GHz)9.9159.91517.15.97
f2 (GHz)10.07510.07517.35.99
Nf8018015333801
L (m)11.81.2751.71
N100180256141
PSNR max 50% (dB)26.0122.7829.9521.16

Share and Cite

MDPI and ACS Style

Pieraccini, M.; Rojhani, N.; Miccinesi, L. Compressive Sensing for Ground Based Synthetic Aperture Radar. Remote Sens. 2018, 10, 1960. https://doi.org/10.3390/rs10121960

AMA Style

Pieraccini M, Rojhani N, Miccinesi L. Compressive Sensing for Ground Based Synthetic Aperture Radar. Remote Sensing. 2018; 10(12):1960. https://doi.org/10.3390/rs10121960

Chicago/Turabian Style

Pieraccini, Massimiliano, Neda Rojhani, and Lapo Miccinesi. 2018. "Compressive Sensing for Ground Based Synthetic Aperture Radar" Remote Sensing 10, no. 12: 1960. https://doi.org/10.3390/rs10121960

APA Style

Pieraccini, M., Rojhani, N., & Miccinesi, L. (2018). Compressive Sensing for Ground Based Synthetic Aperture Radar. Remote Sensing, 10(12), 1960. https://doi.org/10.3390/rs10121960

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop