Next Article in Journal
Effect of Incidence Angle on Temperature Measurement of Solar Panel with Unmanned Aerial Vehicle-Based Thermal Infrared Camera
Next Article in Special Issue
Innovative Rotating SAR Mode for 3D Imaging of Buildings
Previous Article in Journal
Temperature Effects in AMSR2 Soil Moisture Products and Development of a Removal Method Using Data at Ascending and Descending Overpasses
Previous Article in Special Issue
Synthetic Aperture Ladar Motion Compensation Method Based on Symmetric Triangle Linear Frequency Modulation Continuous Wave Segmented Interference
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Iterative Adaptive Based Multi-Polarimetric SAR Tomography of the Forested Areas

1
College of Electronic and Information Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 211106, China
2
Key Laboratory of Radar Imaging and Microwave Photonics, Ministry of Education, Nanjing University of Aeronautics and Astronautics, Nanjing 211106, China
3
Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100094, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(9), 1605; https://doi.org/10.3390/rs16091605
Submission received: 9 March 2024 / Revised: 27 April 2024 / Accepted: 29 April 2024 / Published: 30 April 2024
(This article belongs to the Special Issue Advances in Synthetic Aperture Radar Data Processing and Application)

Abstract

:
Synthetic aperture radar tomography (TomoSAR) is an extension of synthetic aperture radar (SAR) imaging. It introduces the synthetic aperture principle into the elevation direction to achieve three-dimensional (3-D) reconstruction of the observed target. Compressive sensing (CS) is a favorable technology for sparse elevation recovery. However, for the non-sparse elevation distribution of the forested areas, if CS is selected to reconstruct it, it is necessary to utilize some orthogonal bases to first represent the elevation reflectivity sparsely. The iterative adaptive approach (IAA) is a non-parametric algorithm that enables super-resolution reconstruction with minimal snapshots, eliminates the need for hyperparameter optimization, and requires fewer iterations. This paper introduces IAA to tomographicinversion of the forested areas and proposes a novel multi-polarimetric-channel joint 3-D imaging method. The proposed method relies on the characteristics of the consistent support of the elevation distribution of different polarimetric channels and uses the L 2 -norm to constrain the IAA-based 3-D reconstruction of each polarimetric channel. Compared with typical spectral estimation (SE)-based algorithms, the proposed method suppresses the elevation sidelobes and ambiguity and, hence, improves the quality of the recovered 3-D image. Compared with the wavelet-based CS algorithm, it reduces computational cost and avoids the influence of orthogonal basis selection. In addition, in comparison to the IAA, it demonstrates greater accuracy in identifying the support of the elevation distribution in forested areas. Experimental results based on BioSAR 2008 data are used to validate the proposed method.

1. Introduction

In the realm of earth observation, optical, light detection and ranging (LiDAR), and synthetic aperture radar (SAR) stand as the primary tools for remote sensing. While optical images can obtain forest canopy data, they are insufficient in providing the crucial vertical structural parameters [1]. LiDAR, with its high precision and ability to penetrate forest canopies, is capable of acquiring vertical structural information in dense forests. However, the high cost of data acquisition limits its application [2,3]. In contrast, SAR offers capabilities for all-day and all-weather operations, as well as the ability to penetrate vegetation. Long-waveband SAR can acquire internal structural information about forests [4,5]. Therefore, SAR has become an important tool for extracting forest structure information.
Traditional SAR imaging projects the scattering characteristics of a three-dimensional (3-D) target onto an azimuth-range plane to obtain a two-dimensional (2-D) image of the observed scene. However, due to the geometric characteristics of side-looking imaging, this projection leads to issues such as layover and shadows, adversely affecting the identification of targets in SAR images. The existing 3-D imaging techniques for forested areas, such as interferometric SAR (InSAR) [6] and polarimetric InSAR (PolInSAR) [7], are constrained by single-baseline observation, leading to limited visibility of the forest’s vertical structure. Furthermore, multi-baseline (MB) InSAR 3-D imaging technology has become a research hotspot. Exploiting the MB acquisitions with slightly different incidence angles, SAR tomography (TomoSAR) extends the aperture synthetic principle into the elevation direction, which is perpendicular to the azimuth-range plane. In lower frequency bands, such as L-band and P-band, forest structure features can be obtained through TomoSAR 3-D imaging. The longer wavelength of the P-band makes it highly sensitive to weak volume scattering, enabling effective ground identification due to its superior penetration capability [8], which is suitable for forest biomass monitoring, e.g., ESA’s BIOMASS mission [9]. In contrast, the L-band is highly sensitive to strong volume scattering, which assists in canopy identification [10,11]. Therefore, different frequency bands emphasize distinct aspects of forest scattering information. Then, the 3-D focused SAR image is obtained by reconstructing the elevation reflectivity function using spectral estimation (SE)- or compressive sensing (CS)-based methods [12,13,14].
SE-based algorithms can be categorized as non-parametric, including beamforming (BF), adaptive beamforming (Capon), and parametric methods like multiple-signal classification (MUSIC), singular-value decomposition (SVD), and truncated SVD (TSVD). In 2000, Reigber and Moreira initially demonstrated the theory of airborne TomoSAR and retrieved the elevation information of targets using the BF algorithm [12]. After introducing self-interference cancellation by weighting the steering vectors according to the covariance matrix, Capon can obtain the elevation image with higher resolution and fewer sidelobes than BF in TomoSAR [15]. In 2003, an SVD-based TomoSAR imaging method was proposed by Fornaro et al. and was further extended to TSVD [16] and SVD-Wiener [17]. Compared to BF, SVD-based inversion has a better performance in sidelobe suppression and slight super-resolving imaging. MUSIC is a model-based SE algorithm introduced to TomoSAR in 2002 [18,19]. In general, it has a better resolution and sidelobe suppression effect in TomoSAR than BF and Capon. Although the above SE-based algorithms can achieve commendable elevation resolution, the development of tomography is still constrained by uneven sampling, limited elevation apertures, and a limited number of baselines. Consequently, these techniques often encounter problems related to scatterer power and position estimation errors, as well as limitations in super-resolving capability. As an important development in sparse signal processing, CS was proposed by Donoho et al. in 2006 [20,21]. It can recover the sparse signal from far fewer samples than the number required by Shannon–Nyquist sampling theorem [22,23]. In urban environments, where the elevation distribution of artificial buildings tends to be sparse, CS is a favorable technique for recovering complex reflectivity functions along the elevation direction. Nevertheless, due to the typically non-sparse elevation distribution within the forested areas, recovering the reflectivity function through CS-based algorithms is challenging. To solve this problem, Aguilera et al. [24] analyzed the scattering mechanisms (SMs) in forested areas, introduced the wavelet basis [25] to the elevation sparse representation, and obtained the high-resolution image by solving an L 1 -norm regularization problem. In 2016, Li et al. proposed a CS-based fully polarimetric TomoSAR inversion method and achieved the 3-D imaging of forested areas based on the framework of CS [26]. In 2020, Bi et al. proposed a wavelet-based L 1 / 2 -regularized CS-TomoSAR imaging method, conducting experiments within the forested areas of Northern Sweden. It has been proven that this method can improve the quality of elevation recovered [27]. Additionally, Cazarra et al. compared BF, Capon, and CS algorithms for 3-D reconstruction of forests using L-band data over Traunstein. It has been proven that the CS-based algorithm can better reconstruct the forest elevation reflectivity, especially in cases with a lower number of acquisitions and complicated environments [28]. In 2022, Cazarra et al. focused on optimizing different wavenumber distributions to achieve a higher vertical resolution of forest structure in TomoSAR [11]. As an extension of TomoSAR, polarimetric SAR tomography can obtain abundant target information and obtain the 3-D structure of the same target under different SMs [29,30]. Similar to TomoSAR, after sparsely representing the elevation distribution of different polarimetric channels through a wavelet basis, CS also can achieve high-resolution 3-D imaging of the non-sparse forested areas [24,31,32,33]. However, the high computational cost is a serious issue with the wavelet-based CS method. The iterative adaptive approach (IAA) is a user parameter-free weighted least-square non-parametric algorithm. It enables super-resolution imaging without the need for hyperparameter optimization. In 2010, Yardibi et al. introduced IAA for array processing, extending to sparse results by using the Bayesian information criterion [34]. Then, Roberts et al. applied IAA to MIMO radar imaging [35]. In 2016, Campo et al. introduced a modified non-parametric IAA for amplitude and phase estimation in TomoSAR imaging. This approach bypasses the preprocessing step in sum of Kronecker product (SKP) decomposition. However, this method is only suitable for cases with few snapshots [36]. Therefore, Peng et al. proposed the non-parametric iterative approach based on maximum likelihood estimation in 2018 and applied it to forest TomoSAR imaging, performing effectively with a large number of snapshots [37]. In 2021, Feng et al. introduced an imaging method combining IAA and the generalized likelihood ratio test to achieve superior-elevation super-resolution in HoloSAR and when applied specifically to sparse scenes [38]. Compared to CS, IAA requires fewer iterations for the signal recovery, thereby reducing the computational cost. Its application in 3-D imaging of forested areas eliminates the need to construct an orthogonal basis for sparse representation of elevation, ensuring that TomoSAR imaging is not affected by these bases.
In this paper, leveraging the consistency of the target elevation support in different polarimetric channels, a novel IAA-based multi-polarimetric TomoSAR imaging method for forested areas is proposed. The proposed method first establishes a polarimetric TomoSAR imaging model for the forests. Then, an IAA-based multi-polarimetric channel joint TomoSAR imaging method is introduced to achieve the recovery of elevation reflectivity of the forested areas. By utilizing the L 2 -norm to constrain the multi-polarimetric results, the proposed method enhances the 3-D reconstruction accuracy, especially in identifying the support of the elevation distribution. The effectiveness of the proposed method is validated through comparisons with the SE-based algorithms [30], CS-based algorithm [28], and classical IAA [38]. The digital surface model (DSM) measured by LiDAR is chosen as the reference canopy height in this paper [2,3].
The rest of this paper is organized as follows. Section 2 provides an introduction to the TomoSAR imaging model as well as the polarimetric TomoSAR imaging mechanism. In Section 3, the proposed method for the TomoSAR inversion of forested areas is demonstrated in detail. Section 4 introduces the complete data-preprocessing process of the used BioSAR dataset. Section 5 presents the experimental results and performance analysis based on the real data. Section 6 discusses the proposed method in the paper and outlines future research. Conclusions are drawn in Section 7 with several useful remarks.

2. Polarimetric SAR Tomography Model

This section first introduces the TomoSAR imaging model of forested areas based on covariance matrices. Then, a multi-polarimetric imaging model is presented by considering the scattering information of all polarimetric channels. Finally, the reflection mechanisms of the forested areas are analyzed briefly.

2.1. Imaging Model

For single-channel SAR, as shown in Figure 1, using the complex image data acquired by the MB observation for the same area with slightly different incidence angles, TomoSAR can synthesize an aperture along the elevation direction s, which is perpendicular to the azimuth-slant range x r plane, where y is the range direction and Δ b is the length of the elevation aperture, and, hence, obtain the 3-D focused SAR image. Let B = b 1 , b 2 , , b M denote the elevation aperture distribution; for the mth SAR acquisition with b m , at a specific azimuth-range cell x 0 , r 0 , the focused measurement g m x 0 , r 0 can be expressed as [12,14]
g m x 0 , r 0 = Δ s γ s exp j 2 π ξ m s d s
where ξ m = 2 b m 2 b m λ r λ r is the spatial (elevation) frequency, γ s is the complex reflectivity function along the elevation s, Δ s is the elevation scope, and λ and r are the wavelength and slant range, respectively. After discretizing γ s along s by s l l = 1 , 2 , , L , the imaging model in (1) can be approximated by
g m δ s · l = 1 L γ s l exp j 2 π ξ m s l
where L is the number of discrete indices in the elevation, and the constant δ s = Δ s / L 1 is the discretization interval. Let g = g 1 , g 1 , , g M T and γ = [ γ s 1 , γ s 1 , , γ s L ] T denote the total data of all baselines and the discrete complex reflectivity function vector at x 0 , r 0 . One can rewrite (2) as
g M × 1 = Φ M × L γ L × 1
where Φ C M × L is the mapping matrix according to the TomoSAR imaging geometry, with Φ m , l = e j k m s l , where k m = 4 π λ r b m is the vertical wavenumber. Polarimetric SAR can record the horizontal and vertical polarization information of microwave signals. Therefore, multi-polarimetric SAR can provide richer terrain information compared to traditional SAR. Different polarization modes emphasize different aspects of forest structure representation, making it a valuable tool for obtaining comprehensive structure information. Then, the multi-polarimetric covariance matrix of the focused measurement can be expressed as [24]
C ch = Φ ch · diag P ch · Φ ch H
where C ch C HH , C HV , C VV ( HH , HV , VV denote different polarimetric channels), Φ ch Φ HH , Φ HV , Φ VV , P ch P HH , P HV , P VV , and, in each ch , C = E gg H , diag P R L × L , with P = γ 1 2 , γ 2 2 , . . . , γ L 2 as its main diagonal and zeros in the off-diagonal elements. Assuming that the target’s backscatter structure approximates similar azimuth angles, and due to the approximate azimuth angles of different polarization channels, the observation matrices of different polarization channels are considered to be uniform. Since the observation matrix of all polarimetric channels is consistent, thus Φ HH = Φ HV = Φ VV = Φ . The theoretical elevation resolution ρ s can be calculated as [14]
ρ s = λ r 2 Δ b
where λ is the wavelength.

2.2. SM Analysis

To achieve the expected recovered result, an analysis of the radar scattering property of the forested area is first performed. As discussed in [30,39,40], there are four dominated SMs in a forested area, i.e., ground backscattering, trunk–ground scattering, canopy–ground scattering, and canopy backscattering (see Figure 2a). Canopy backscattering is a volumetric SM whose phase center is determined by the height of the canopy and is distributed above the ground along the elevation direction. As discussed in [31], although the other three SMs have different scattering characteristics, their phase scattering centers are all located on the ground plane and distributed along the elevation. Actually, the above SM model has been comprehensively validated using different datasets [30,41,42]. According to this analysis, the elevation backscattered power is typically attributed to contributions from both the ground and canopy, as depicted in Figure 2b.

3. Imaging Method

IAA can achieve super-resolution imaging with fewer iterations, eliminating the need for hyperparameter optimization. Consequently, the IAA-based TomoSAR imaging method recovers the elevation signals using only a few 2-D single-look complex (SLC) images. Polarimetric data offer more information for the 3-D recovery of forested areas. The supports of multiple polarimetric channels are consistent; only the scattering intensities are different (see Figure 3). Figure 3 shows the schematic diagram of backscattering coefficients of HH, HV, and VV polarimetric channels. Different colors represent different signal values, whereas white corresponds to a value of 0. The horizontal direction indicates the height, with s 1 and s 2 representing two scattering centers. It can be observed that the height positions of the scattering centers for different polarimetric channels are consistent. According to this property, a novel IAA-based multi-polarimetric TomoSAR imaging method is proposed and used for the high-precision recovery of forested areas. It introduces the L 2 -norm to constrain the polarimetric data to improve the reconstruction accuracy of the target elevation position. In the calculation of the elevation scattering distribution for each pixel unit, we firstly initialize the model based on Equation (4). Then, the IAA algorithm is used to compute the initial scattering distribution along the elevation direction for each polarization channel. Once the elevation scattering distributions for all polarization channels are obtained, an L 2 -norm constraint is utilized to further refine the distribution. Finally, the refined elevation scattering distributions, after meeting the specified criteria, are outputted. The detailed iterative procedure of the proposed method based on the model in Equation (4) is summarized as follows:
(i)
Initialization:
Elevation power distribution:
p ( 0 ) = diag Φ H · C HH + C HV + C VV · Φ ;
Noise vector: d ( 0 ) = 0 M × 1 ;
Identity matrix: V M × M ;
Maximum number of iterations I max ;
Error parameter ε ;
(ii)
Iteration:
While i < I max and Residual > ε ,
R = Φ · diag p ( i ) · Φ H + diag d ( i )
P l , ch = Φ l H · R 1 · C ch · R 1 H · Φ l Φ l H · R 1 · Φ l 2
D m , ch = V m H · R 1 · C ch · R 1 H · V m V m H · R 1 · V m 2
p ( i + 1 ) = P ch 2
d ( i + 1 ) = D ch 2
Residual = p ( i + 1 ) p ( i ) 2
i = i + 1
(iii)
Output:
Reconstructed elevation reflectivity function p ( i + 1 ) .
In each iteration, the L 2 -norm constraint is applied to multiple polarimetric channels, which, ultimately, results in a reconstructed joint elevation power distribution between all channels.

4. Date Preprocessing

In general, the data format typically obtained is SLC image, which contains irrelevant information and cannot be directly used for 3-D imaging. Additionally, data suitable for TomoSAR imaging must ensure coherence between images. Therefore, the preprocessing is essential, as it provides high-quality data for the following 3-D reconstruction. This section introduces the dataset used in this paper and outlines the overall data-processing workflow, with detailed explanations of the data-preprocessing procedures.

4.1. BioSAR 2008 Dataset

The BioSAR 2008 dataset was acquired by the E-SAR sensor of the German Aerospace Agency (DLR) in the northern forests of Sweden [43]. BioSAR provides data in both the P-band and L-band, with the P-band data specifically suited for forest biomass monitoring [8]. This paper aims to obtain the 3-D scattering structure and canopy height. Therefore, the BioSAR 2008 airborne multi-polarimetric L-band dataset is used to verify the proposed method. It covers terrain exhibiting a significant elevation variation, ranging from 100 to 400 m. The parameters of this campaign are listed in Table 1. The elevation aperture position is shown in Figure 4.

4.2. Preprocessing

This section delineates the specific steps involved in the 3-D imaging of the forested areas, as depicted in Figure 5. Firstly, multiple SLC images of the same observed scene at different incidence angles are inputted to the proposed method. Then, a series of data-preprocessing steps are performed on these datasets. Finally, we utilize the proposed method for tomographic inversion to obtain the 3-D scattering information of the observed target. This preprocessing [44] part consists of five key steps—data registration, interference analysis, phase flattening, topographic phase removal, and filtering and sampling—as follows:
(i) 
Data registration. This requires one scene to be designated as the master image and the rest of the scenes to be aligned with the master one as slave images, ensuring that the elevation direction of multiple 2-D complex image data in each pixel cell is consistent;
(ii) 
Interference analysis. This validates the coherence and phase information between images to ensure their suitability for 3-D imaging;
(iii) 
Phase flattening. Due to the wide coverage of the surveillance region, spanning over 2000 m in slant range, this step aims to rectify phase discrepancies caused by slant range;
(iv) 
Topographic phase removal. Due to significant terrain variations in the surveillance region, it is crucial to eliminate phase changes caused by terrain alterations for obtaining information about surface forests. High-precision digital elevation model (DEM) data are necessary for this step. DEM is estimated from the laser mapping of Krycklan, with the ground level subtracted. It is presented in a grid size of 0.5 m × 0.5 m, using the UTM Zone 34N geographic datum;
(v) 
Filtering and sampling. This is employed to eliminate the noise and, hence, improve the quality of 2-D SAR images.
After the above preprocessing, the amplitude–phase results of the HH, HV, and VV polarimetric channels in the experimental scenario (see Figure 6) are shown in Figure 7. The lower left corner of Figure 7 shows the interference phase between pairs of six images, while the upper right corner displays their corresponding coherence results. The coherence results are represented as grayscale images, where lighter shades indicate higher coherence. Subsequently, 3-D imaging is conducted on the preprocessed data. Initially, the preprocessed data are inputted, and a specialized TomoSAR imaging model is developed for the forested areas to effectively mitigate coherent speckles along the elevation direction. Then, the proposed method is applied to derive the height scattering distribution. Finally, a 3-D point cloud of the forested area is obtained.

5. Experimental Results

The panoramic polarimetric SAR image of the surveillance area of the BioSAR 2008 campaign is presented in Figure 6, where the x-axis corresponds to the azimuth direction and the y-axis denotes the range direction. The yellow-highlighted azimuthal slices in Figure 6 denote the chosen experimental segment. Slice 1 consists of both ground and forest parts. Slice 2 reveals intricate scattering details within the forest. In order to validate the proposed method, the reconstructed results of three SE-based algorithms [12,15,18,19], a CS-based algorithm [27,28], and IAA [38] are used for a fair comparison. Considering the usage of the multi-polarimetric data in the proposed method, Figure 8 presents a comparison of the incoherent sum of the results for all polarimetric channels. The white line in these figures represents the processed LiDAR DSM data, utilized as the reference height. The original reference data are the DSM estimated from the laser mapping of Krycklan, with the ground level subtracted. Information is presented in a grid size of 0.5 m × 0.5 m, using the UTM Zone 34N geographic datum. To maintain the accuracy of the experiments, the original reference data underwent the same filtering and sampling as applied to the experimental data in the paper, thereby excluding the influence of other factors. Figure 9 depicts the corresponding results for Slice 2.
From the comparison of slices, it can be seen that the proposed method has better resolving ability than the three SE-based algorithms. Simultaneously, it demonstrates a more accurate recognition of the canopy height. For the SE-based algorithms, it is shown that MUSIC stands out for its precise description of the canopy of the forests. However, it suffers from considerable sidelobe artifacts and lacks accuracy in identifying scattering intensity as well as ground scattering center. Conversely, the proposed method effectively suppresses these sidelobes and can acquire more accurate canopy and ground intensity information. For the CS-TomoSAR imaging, the wavelet-based L 1 -norm regularization method for the scene recovery is employed. It introduces Daubechies Symmlet wavelets [24] as the sparse basis to represent the elevation distribution. To use this method for TomoSAR inversion, it is necessary to solve the optimization problem of wavelet coefficients and then use the orthogonal basis and wavelet coefficients to obtain the elevation reflection function. Although it provides excellent resolution capabilities, this method produces anomalous artifacts during the elevation reconstruction, resulting in suboptimal results. Moreover, it is computationally expensive. Furthermore, compared to the IAA, as shown in the results for Slice 1, it is found that the proposed method improves the suppression of error information around the 5250 m position. Similar advantages also can be seen in the far-range areas of the recovered images from Slice 2. In the complicated elevation scattering distribution within the forested areas, the proposed method shows commendable performance, effectively mitigating the sidelobes, especially in far-range areas.
To further demonstrate the effectiveness of the proposed method, root mean square error (RMSE) of the altitude estimation error and computational time are used to quantitatively compare the performance of different methods. In this paper, the canopy height reconstructed by TomoSAR is represented by the vertical position of the scattering centers in the canopy, as shown in Figure 2b. Table 2 lists the RMSE for Slice 1 and Slice 2, calculated using different methods, along with the average computation time per azimuth-range cell. From Table 2, it is seen that the SE-based algorithms, i.e., BF, Capon, and MUSIC, have lower computation times in multi-polarimetric channels. Among these, MUSIC is able to estimate the canopy height most accurately, i.e., it has the smallest RMSE value. Even in Slice 2, the MUSIC algorithm in multi-polarimetric channels exhibits the lowest RMSE. However, its performance in identifying the elevation distribution of forests is still slightly inferior to the proposed method, particularly in ground scattering centers and scattering intensity identification. Compared to the wavelet-based L 1 algorithm, the proposed method not only significantly enhances computational efficiency, but also notably reduces the RMSE. Furthermore, compared to IAA, it still provides a better estimate of height, but the calculation time is slightly longer. The results in Table 2 are consistent with the above experimental results and analysis. It is shown that the proposed method can achieve higher accuracy in height estimation, albeit with slightly longer computation times. However, compared with the wavelet-based L 1 algorithm, its computational cost is still considerable. Figure 10 shows the 3-D point cloud of the whole forested area reconstructed by the proposed method. It vividly demonstrates the height and distribution of trees in densely forested areas and flat terrain.

6. Discussion

The proposed method in this paper combines multi-channel data. While the primary aim of this paper is to acquire 3-D height information of the forest canopy, it overlooks the differences between various polarization channels. However, each polarization channel emphasizes different aspects of acquiring vertical structural information about forests. Therefore, future research will consider preserving the differences between polarization channels while obtaining a more comprehensive understanding of the vertical structure of the forested areas. Additionally, applying P-band data in forest biomass estimation is also a key focus of future research efforts.

7. Conclusions

This paper focuses on acquiring high-precision 3-D scattering information in an extensive forested area. It demonstrates a comprehensive procedural approach that starts with multiple SLC data and ultimately generates a 3-D point cloud representation of the forest. The multi-polarimetric data enriches the information base, which is essential for capturing 3-D scattering details. In this paper, a novel IAA-based multi-polarimetric TomoSAR imaging method for forested areas is proposed. It first establishes a polarimetric TomoSAR imaging model for the forests. Then, an IAA-based multi-polarimetric channel joint TomoSAR imaging method is introduced to achieve the recovery of elevation reflectivity. By using the L 2 -norm to constrain the multi-polarimetric results, the proposed method achieves high-precision and 3-D reconstruction accuracy, especially in identifying the support of the elevation distribution. Compared with SE-based algorithms, it suppresses the elevation sidelobes and ambiguity dramatically. Compared with the wavelet-based CS algorithm, it has reduced computational cost and avoids the influence of orthogonal-basis selection. In addition, in comparison to the IAA algorithm, it demonstrates greater accuracy in identifying the support of the elevation distribution in the forested areas. Experimental results based on real BioSAR 2008 data validate the proposed method.

Author Contributions

S.J. and H.B. conceived the article. S.J., J.Z., W.H. and H.B. processed the BioSAR data and performed related experiments. S.J., H.B. and Q.G. participated in the writing of this article. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Foundation of China under Grant 62271248 and 61901213, in part by the Natural Science Foundation of Jiangsu Province under Grant BK20230090, and in part by the Aeronautical Science Foundation of China under Grant 201920052001.

Data Availability Statement

The E-SAR data was provided by the European Space Agency (ESA) under the BIOSAR 2008 campaign. Available online: https://earth.esa.int/eogateway/campaigns/biosar-2 (accessed on 26 April 2024).

Acknowledgments

The authors would like to thank Dragon 3 Project (ID10609) and Chen Erxue for providing the BioSAR 2008 dataset.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
LiDARLight Detection And Ranging
SARSynthetic Aperture Radar
3-DThree-Dimensional
2-DTwo-Dimensional
InSARInterferometric SAR
PolInSARPolarimetric InSAR
MBMulti-Baseline
TomoSARSAR Tomography
SESpectral Estimation
CSCompressive Sensing
BFBeamforming
CaponAdaptive Beamforming
MUSICMultiple Signal Classification
SMsScattering Mechanisms
IAAIterative Adaptive Approach
SKPSum of Kronecker Product
SLCSingle-Look Complex
DSMDigital Surface Model
DEMDigital Elevation Model
RMSERoot Mean Square Error

References

  1. Van Tricht, K.; Gobin, A.; Gilliams, S.; Piccard, I. Synergistic use of radar Sentinel-1 and optical Sentinel-2 imagery for crop mapping: A case study for Belgium. Remote Sens. 2018, 10, 1642. [Google Scholar] [CrossRef]
  2. Yan, Z.; Liu, R.; Cheng, L.; Zhou, X.; Ruan, X.; Xiao, Y. A concave hull methodology for calculating the crown volume of individual trees based on vehicle-borne LiDAR data. Remote Sens. 2019, 11, 623. [Google Scholar] [CrossRef]
  3. Tarsha Kurdi, F.; Lewandowicz, E.; Shan, J.; Gharineiat, Z. Three-dimensional modeling and visualization of single tree LiDAR point cloud using matrixial form. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2024, 17, 3010–3022. [Google Scholar] [CrossRef]
  4. Hirschmugl, M.; Deutscher, J.; Sobe, C.; Bouvet, A.; Mermoz, S.; Schardt, M. Use of SAR and optical time series for tropical forest disturbance mapping. Remote Sens. 2020, 12, 727. [Google Scholar] [CrossRef]
  5. Xu, Z.; Wang, Y. Radar satellite image time series analysis for high-resolution mapping of man-made forest change in Chongming Eco-Island. Remote Sens. 2020, 12, 3438. [Google Scholar] [CrossRef]
  6. Zebker, H.A. Goldstein, R.M. Topographic mapping from interferometric synthetic aperture radar observations. J. Geophys. 1986, 91, 4993–4999. [Google Scholar] [CrossRef]
  7. Xue, F.; Wang, X.; Xu, F.; Wang, R. Polarimetric SAR interferometry: A tutorial for analyzing system parameters. IEEE Geosci. Remote Sens. Mag. 2020, 8, 83–107. [Google Scholar] [CrossRef]
  8. Santi, E.; Paloscia, S.; Pettinato, S.; Cuozzo, G.; Padovano, A.; Notarnicola, C.; Albinet, C. Machine-Learning Applications for the Retrieval of Forest Biomass from Airborne P-Band SAR Data. Remote Sens. 2020, 12, 804. [Google Scholar] [CrossRef]
  9. Banda, F.; Giudici, D.; Le Toan, T.; Mariotti d’Alessandro, M.; Papathanassiou, K.; Quegan, S.; Riembauer, G.; Scipal, K.; Soja, M.; Tebaldini, S.; et al. The BIOMASS level 2 prototype processor: Design and experimental results of above-ground biomass estimation. Remote Sens. 2020, 12, 985. [Google Scholar] [CrossRef]
  10. Tello, M.; Cazcarra-Bes, V.; Pardini, M.; Papathanassiou, K. Forest structure characterization from SAR tomography at L-band. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2018, 11, 3402–3414. [Google Scholar] [CrossRef]
  11. Cazcarra-Bes, V.; Pardini, M.; Papathanassiou, K. Definition of tomographic SAR configurations for forest structure applications at L-band. IEEE Geosci. Remote Sens. Lett. 2022, 19, 1–5. [Google Scholar] [CrossRef]
  12. Reigber, A.; Moreira, A. First demonstration of airborne SAR tomography using multibaseline L-band data. IEEE Trans. Geosci. Remote Sens. 2000, 38, 2142–2152. [Google Scholar] [CrossRef]
  13. Fornaro, G.; Serafino, F.; Lombardini, F. Three-dimensional multipass SAR focusing: Experiments with long-term spaceborne data. IEEE Trans. Geosci. Remote Sens. 2005, 43, 702–714. [Google Scholar] [CrossRef]
  14. Zhu, X.; Bamler, R. Tomographic SAR inversion by L1-norm regularization-the compressive sensing approach. IEEE Trans. Geosci. Remote Sens. 2010, 48, 3839–3846. [Google Scholar] [CrossRef]
  15. Lombardini, F.; Gini, F.; Matteucci, P. Application of array processing techniques to multibaseline InSAR for layover solution. In Proceedings of the 2001 IEEE Radar Conference, Atlanta, GA, USA, 3 May 2001. [Google Scholar]
  16. Fornaro, G.; Serafino, F.; Soldovieri, F. Three-dimensional focusing with multipass SAR data. IEEE Trans. Geosci. Remote Sens. 2003, 41, 507–517. [Google Scholar] [CrossRef]
  17. Zhu, X.X.; Bamler, R. Very high resolution spaceborne SAR tomography in urban environment. IEEE Trans. Geosci. Remote Sens. 2010, 48, 4296–4308. [Google Scholar] [CrossRef]
  18. Gini, F.; Lombardini, F. Multilook APES for multibaseline SAR interferometry. IEEE Trans. Signal Process. 2002, 50, 1800–1803. [Google Scholar] [CrossRef]
  19. Gini, F.; Lombardini, F. Multibaseline cross-track SAR interferometry: A signal processing perspective. IEEE Aerosp. Electron. Syst. Mag. 2005, 20, 71–93. [Google Scholar] [CrossRef]
  20. Donoho, D.L. Compressed sensing. IEEE Trans. Inf. Theory 2006, 52, 1289–1306. [Google Scholar] [CrossRef]
  21. Candes, E.J.; Tao, T. Near-optimal signal recovery from random projections: Universal encoding strategies. IEEE Trans. Inf. Theory 2006, 52, 5406–5425. [Google Scholar] [CrossRef]
  22. Nyquist, H. Certain topics in telegraph transmission theory. Trans. Am. Inst. Electr. Eng. 1928, 47, 617–644. [Google Scholar] [CrossRef]
  23. Shannon, C.E. Communication in the presence of noise. Proc. IRE 1949, 37, 10–21. [Google Scholar] [CrossRef]
  24. Aguilera, E.; Nannini, M.; Reigber, A. Wavelet-based compressed sensing for SAR tomography of forested areas. IEEE Trans. Geosci. Remote Sens. 2013, 51, 5283–5295. [Google Scholar] [CrossRef]
  25. Strang, G.; Nguyen, T. Wavelet and Filter Banks; MIT: Wellesley, MA, USA, 1997. [Google Scholar]
  26. Li, X.; Liang, L.; Guo, H.; Huang, Y. Compressive sensing for multibaseline polarimetric SAR tomography of forested areas. IEEE Trans. Geosci. Remote Sens. 2016, 54, 153–166. [Google Scholar] [CrossRef]
  27. Bi, H.; Cheng, Y.; Zhu, D.Y.; Hong, W. Wavelet-based L1/2 regularization for CS-TomoSAR imaging of forested area. J. Syst. Eng. Electron. 2020, 31, 1160–1166. [Google Scholar]
  28. Cazcarra-Bes, V.; Pardini, M.; Tello, M.; Papathanassiou, K.P. Comparison of tomographic SAR reflectivity reconstruction algorithms for forest applications at L-band. IEEE Trans. Geosci. Remote Sens. 2020, 58, 147–164. [Google Scholar] [CrossRef]
  29. Tebaldini, S. Algebraic synthesis of forest scenarios from multibaseline polInSAR data. IEEE Trans. Geosci. Remote Sens. 2009, 47, 4132–4142. [Google Scholar] [CrossRef]
  30. Tebaldini, S. Single and multipolarimetric SAR tomography of forested areas: A parametric approach. IEEE Trans. Geosci. Remote Sens. 2010, 48, 2375–2387. [Google Scholar] [CrossRef]
  31. Aguilera, E.; Nannini, M.; Reigber, A. A data-adaptive compressed sensing approach to polarimetric SAR tomography of forested areas. IEEE Geosci. Remote Sens. Lett. 2013, 10, 543–547. [Google Scholar] [CrossRef]
  32. Aguilera, E.; Nannini, M.; Reigber, A. Multisignal compressed sensing for polarimetric SAR tomography. IEEE Geosci. Remote Sens. Lett. 2012, 9, 871–875. [Google Scholar] [CrossRef]
  33. Ponce, O.; Prats-Iraola, P.; Scheiber, R.; Reigber, A.; Moreira, A. First airborne demonstration of Holographic SAR tomography with fully polarimetric multicircular acquisitions at L-Band. IEEE Trans. Geosci. Remote Sens. 2016, 54, 6170–6196. [Google Scholar] [CrossRef]
  34. Yardibi, T.; Li, J.; Stoica, P.; Xue, M.; Baggeroer, B. Source localization and sensing: A nonparametric iterative adaptive approach based on weighted least squares. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 425–443. [Google Scholar] [CrossRef]
  35. Roberts, W.; Stoica, P.; Li, J.; Yardibi, T.; Sadjadi, F.A. Iterative adaptive approaches to MIMO radar imaging. IEEE J. Sel. Top. Signal Process. 2010, 4, 5–20. [Google Scholar] [CrossRef]
  36. Campo, G.D.M.d.; Reigber, A.; Shkvarko, Y.V. Resolution enhanced SAR tomography a nonparametric iterative adaptive approach. In Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Beijing, China, 10–15 July 2016. [Google Scholar]
  37. Peng, X.; Li, X.; Wang, C.; Fu, H.; Du, Y. A maximum likelihood based nonparametric iterative adaptive method of synthetic aperture radar tomography and its application for estimating underlying topography and forest height. Sensors 2018, 18, 2459. [Google Scholar] [CrossRef] [PubMed]
  38. Feng, D.; An, D.; Chen, L.; Huang, X. Holographic SAR tomography 3-D reconstruction based on iterative adaptive approach and generalized likelihood ratio test. IEEE Trans. Geosci. Remote Sens. 2021, 59, 305–315. [Google Scholar] [CrossRef]
  39. Lin, Y.; Sarabandi, K. Electromagnetic scattering model for a tree trunk above a tilted ground plane. IEEE Trans. Geosci. Remote Sens. 1995, 33, 1063–1070. [Google Scholar]
  40. Sarabandi, K. Scattering from dielectric structures above impedance surfaces and resistive sheets. IEEE Trans. Antennas Propag. 1992, 40, 67–78. [Google Scholar] [CrossRef]
  41. Nannini, M.; Scheiber, R.; Horn, R.; Moreira, A. First 3-D reconstructions of targets hidden beneath foliage by means of polarimetric SAR tomography. IEEE Geosci. Remote Sens. Lett. 2012, 9, 60–64. [Google Scholar] [CrossRef]
  42. Tebaldini, S.; Rocca, F. Multibaseline polarimetric SAR tomography of a boreal Forest at P- and L-Bands. IEEE Trans. Geosci. Remote Sens. 2012, 50, 232–246. [Google Scholar] [CrossRef]
  43. Hajnsek, I.; Scheiber, R.; Keller, M.; Horn, R.; Lee, S.; Ulander, L.M.H.; Gustavsson, A.; Sandberg, G.; Le Toan, T.; Tebaldini, S.; et al. BioSAR 2008: Final Report; Technical Report; ESA-ESTEC: Noordwijk, The Netherlands, 2009. [Google Scholar]
  44. Ferretti, A.; Monti-Guarnieri, A.; Prati, C.; Rocca, F.; Massonnet, D. InSAR Principles: Guidelines for SAR Interferometry Processing and Interpretation; ESA: Paris, France, 2007. [Google Scholar]
Figure 1. TomoSAR imaging geometry.
Figure 1. TomoSAR imaging geometry.
Remotesensing 16 01605 g001
Figure 2. Scattering Scattering of a forested area. (a) Scattering mechanism, (b) Scattering distribution.
Figure 2. Scattering Scattering of a forested area. (a) Scattering mechanism, (b) Scattering distribution.
Remotesensing 16 01605 g002
Figure 3. Schematic diagram of backscattering coefficients of HH, HV, and VV polarimetric channels.
Figure 3. Schematic diagram of backscattering coefficients of HH, HV, and VV polarimetric channels.
Remotesensing 16 01605 g003
Figure 4. Elevation aperture position in the BioSAR 2008 dataset.
Figure 4. Elevation aperture position in the BioSAR 2008 dataset.
Remotesensing 16 01605 g004
Figure 5. Implementation process of the TomoSAR 3-D imaging of the forested areas based on the proposed method.
Figure 5. Implementation process of the TomoSAR 3-D imaging of the forested areas based on the proposed method.
Remotesensing 16 01605 g005
Figure 6. Polarimetric SAR image of the surveillance area (The yellow area numbers 1 and 2 respectively represent the two slices selected for the experiment).
Figure 6. Polarimetric SAR image of the surveillance area (The yellow area numbers 1 and 2 respectively represent the two slices selected for the experiment).
Remotesensing 16 01605 g006
Figure 7. Amplitude and phase results after data preprocessing for the (a) HH, (b) HV, and (c) VV polarimetric channels.
Figure 7. Amplitude and phase results after data preprocessing for the (a) HH, (b) HV, and (c) VV polarimetric channels.
Remotesensing 16 01605 g007aRemotesensing 16 01605 g007b
Figure 8. The incoherent sum of the results for all polarization channels (Slice 1). (a) BF. (b) Capon. (c) MUSIC. (d) Wavelet-based L 1 . (e) IAA. (f) The proposed method. The white line represents the LiDAR DSM.
Figure 8. The incoherent sum of the results for all polarization channels (Slice 1). (a) BF. (b) Capon. (c) MUSIC. (d) Wavelet-based L 1 . (e) IAA. (f) The proposed method. The white line represents the LiDAR DSM.
Remotesensing 16 01605 g008aRemotesensing 16 01605 g008b
Figure 9. The incoherent sum of the results for all polarization channels (Slice 2). (a) BF. (b) Capon. (c) MUSIC. (d) Wavelet-based L 1 . (e) IAA. (f) The proposed method. The white line represents the LiDAR DSM.
Figure 9. The incoherent sum of the results for all polarization channels (Slice 2). (a) BF. (b) Capon. (c) MUSIC. (d) Wavelet-based L 1 . (e) IAA. (f) The proposed method. The white line represents the LiDAR DSM.
Remotesensing 16 01605 g009aRemotesensing 16 01605 g009b
Figure 10. 3-D point cloud map of the entire surveillance region reconstructed by the proposed method.
Figure 10. 3-D point cloud map of the entire surveillance region reconstructed by the proposed method.
Remotesensing 16 01605 g010
Table 1. Parameters of BioSAR 2008 dataset.
Table 1. Parameters of BioSAR 2008 dataset.
Parameter    Value
Tracks    6
Radar center frequency    1.3 GHz
Center slant range    ≈4500 m
Slant range resolution    1.5 m
Azimuth resolution    1.6 m
Height resolution    6∼25 m (near range to far range)
Table 2. RMSE of the height estimation error and computational time of different TomoSAR imaging methods.
Table 2. RMSE of the height estimation error and computational time of different TomoSAR imaging methods.
AlgorithmRMSE [m]Time [s]
Slice 1Slice 2
BF (All channels)16.7113.080.003
Capon (All channels)9.1011.690.003
MUSIC (All channels)6.425.460.004
Wavelet-based L 1 (All channels)10.3112.950.080
IAA (All channels)4.936.220.007
The proposed method4.575.580.021
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jin, S.; Bi, H.; Guo, Q.; Zhang, J.; Hong, W. Iterative Adaptive Based Multi-Polarimetric SAR Tomography of the Forested Areas. Remote Sens. 2024, 16, 1605. https://doi.org/10.3390/rs16091605

AMA Style

Jin S, Bi H, Guo Q, Zhang J, Hong W. Iterative Adaptive Based Multi-Polarimetric SAR Tomography of the Forested Areas. Remote Sensing. 2024; 16(9):1605. https://doi.org/10.3390/rs16091605

Chicago/Turabian Style

Jin, Shuang, Hui Bi, Qian Guo, Jingjing Zhang, and Wen Hong. 2024. "Iterative Adaptive Based Multi-Polarimetric SAR Tomography of the Forested Areas" Remote Sensing 16, no. 9: 1605. https://doi.org/10.3390/rs16091605

APA Style

Jin, S., Bi, H., Guo, Q., Zhang, J., & Hong, W. (2024). Iterative Adaptive Based Multi-Polarimetric SAR Tomography of the Forested Areas. Remote Sensing, 16(9), 1605. https://doi.org/10.3390/rs16091605

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