Next Article in Journal
High-Resolution ISAR Imaging Based on Plug-and-Play 2D ADMM-Net
Next Article in Special Issue
Mining Cross-Domain Structure Affinity for Refined Building Segmentation in Weakly Supervised Constraints
Previous Article in Journal
Multitask Learning-Based for SAR Image Superpixel Generation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Learning a Fully Connected U-Net for Spectrum Reconstruction of Fourier Transform Imaging Spectrometers

1
Xi’an Institute of Optics and Precision Mechanics, Chinese Academy of Sciences, Xi’an 710119, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Pilot National Laboratory for Marine Science and Technology, Qingdao 266237, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(4), 900; https://doi.org/10.3390/rs14040900
Submission received: 25 December 2021 / Revised: 6 February 2022 / Accepted: 10 February 2022 / Published: 14 February 2022

Abstract

:
Fourier transform imaging spectrometers (FTISs) are widely used in global hyperspectral remote sensing due to the advantages of high stability, high throughput, and high spectral resolution. Spectrum reconstruction (SpecR) is a classic problem of FTISs determining the acquired data quality and application potential. However, the state-of-the-art SpecR algorithms were restricted by the length of maximum optical path difference (MOPD) of FTISs and apodization processing, resulting in a decrease in spectral resolution; thus, the applications of FTISs were limited. In this study, a deep learning SpecR method, which directly learned an end-to-end mapping between the interference/spectrum information with limited MOPD and without apodization processing, was proposed. The mapping was represented as a fully connected U-Net (FCUN) that takes the interference fringes as the input and outputs the highly precise spectral curves. We trained the proposed FCUN model using the real spectra and simulated pulse spectra, as well as the corresponding simulated interference curves, and achieved good results. Additionally, the performance of the proposed FCUN on real interference and spectral datasets was explored. The FCUN could obtain similar spectral values compared with the state-of-the-art fast Fourier transform (FFT)-based method with only 150 and 200 points in the interferograms. The proposed method could be able to enhance the resolution of the reconstructed spectra in the case of insufficient MOPD. Moreover, the FCUN performed well in visual quality using noisy interferograms and gained nearly 70% to 80% relative improvement over FFT for the coefficient of mean relative error (MRE). All the results based on simulated and real satellite datasets showed that the reconstructed spectra of the FCUN were more consistent with the ideal spectrum compared with that of the traditional method, with higher PSNR and lower values of spectral angle (SA) and relative spectral quadratic error (RQE).

1. Introduction

Imaging interferometers, also called Fourier transform imaging spectrometers (FTISs), have several superiorities, including high throughput, high wavenumber accuracy, high resolution, lower weight, and smaller size [1,2,3]. Due to these advantages, FTISs have been widely used in remote sensing of the moon, earth, and space [4,5]. With rapid developments in technology, various types of interferometric spectrometers that can obtain tens or even hundreds of hyperspectral data [6] have been established, such as the temporary modulated imaging interferometer proposed by Michelson, the spatial interferometric imaging spectrometer, and the spatially-temporally modulated interferometric imaging spectrometer [7].
Although there are currently a great variety of FTISs, the key part of the FTIS is the interferometer [8,9,10,11]. Figure 1 shows the optical layout of a type of FTIS, a large aperture static imaging spectrometry (LASIS) system. In theory and practice, the interferometer divides the light from the detection object into two coherent beams that will interfere on the optical sensors. By altering the optical path difference (OPD) of two beams of light, we will acquire a set of interference fringes [12]. To obtain the spectral information, the spectrum reconstruction (SpecR) is used to transform interference fringes into spectral information, and the SpecR algorithm determines the quality and the further application of the spectral data [13,14].
The fast Fourier transform (FFT) plays a vital role in the SpecR of FTISs, which correlates spectral profiles with interference fringes [15]. However, the spectrum recovered by FFT using the original interferogram may result in many side lobes. This is mainly because the interferogram cannot be scanned to infinity as the theoretical calculation of FFT, and then the value of interference fringes, out of the scope, is not recorded, which would cause spectrum leakage. The side lobes would reduce the accuracy of the spectral amplitude value. To solve the spectrum leakage of SpecR, many apodization functions (also called window functions) [16,17,18], such as triangle, Hamming, Hanning, and Happ–Genzel algorithms, have been applied to smooth down an interference curve to zero at the end of the recorded region. Even though apodization algorithms can reduce spectrum leakage, the full width at half maximum (FWHM) of the spectrum becomes large, resulting in a decrease in the spectral resolution. Thus, the peaks and valleys are hard to identify in the reconstructed spectrum. In addition, the spectral resolution of the FFT relies on the length of the interference data. However, long interference data cannot be transmitted due to bandwidth limitations in practical remote sensing applications and insufficient sampling points to record interference signals. When the signal is not long enough, it is difficult for FFT-based methods to recover the high precision spectrum [19]. The multiple signal classification (MUSIC) algorithm, which is maturely used in modern spectrum estimation, was introduced in the SpecR process to improve the resolution. Jian et al. [20] introduced the MUSIC algorithm and the autoregressive (AR) [21,22] model for better performance in the spectral recovery of the pulse signal. These algorithms all achieved good spectrum reconstruction in terms of resolution, however, they were not suited for polychromatic SpecR, particularly for low signal-to-noise ratio data. Hence, currently most of the methods for SpecR are still based on FFT [1,3,4].
Recently, deep learning has been well developed in the fields of image processing [23], such as image restoration [24], super-resolution [25], object detection [26], image segmentation [27], image classification [28], and remote sensing [29]. More closely related to our study task, the multilayer fully connected network (FCN) was applied for high and low-level computer vision problems. Although the FCN has the disadvantages of a large network and many parameters when processing image signals, its fully connected structure can well represent the relationship between the interference curve and spectrum information. The key process of SpecR, the Fourier transform, can be regarded as the linear FCN. To date, deep learning techniques have been successfully used for many studies of FTIS, such as medical image classification [30], digital staining [31], and diagnosis of abnormal thyroid function [32]. However, to the best of our knowledge, deep learning techniques have not been used in the SpecR of FTISs.
The transformation of the interference curve to a spectral curve using the Fourier transform can be defined as an end-to-end model, and it is possible for deep learning to achieve significant improvements in computational accuracy for the end-to-end model [33]. Motivated by this fact, we designed the fully connected U-Net (named FCUN) as an end-to-end mapping between interference fringes and spectral profiles. The architecture can extract and fuse multiscale features. The FCUN was fundamentally different from existing SpecR methods without needing too many intermediate steps, such as apodization or autoregression for recovering the spectrum. Furthermore, the intermediate steps were formulated as fully connected layers, which were involved in the optimization. In the proposed method, the entire SpecR pipeline was fully obtained through learning, with little pre/post-processing.
The main contributions of this research are listed as follows:
  • A fully connected U-Net was proposed for spectrum reconstruction of FTISs. The network directly learned an end-to-end mapping between the interference/spectrum information.
  • We proved that deep learning could improve the performance of spectrum reconstruction of FTISs, and could achieve good quality spectra. This guided the design of the algorithm of SpecR.
  • The proposed FCUN method can recover spectra that exceed the physical spectral resolution of the FTISs, within the limited MOPD of interference fringes.
  • Two datasets of SpecR were built for the first time and could be applied for further studies.

2. Theoretical Foundation

The traditional data processing of FTISs and limitations were introduced in the first part. Then, the relationship between FTIS data processing and the deep learning method was illustrated. The architecture of a typical U-Net was also described.

2.1. Traditional Data Processing of FTIS

Derived from Fourier transform spectroscopy theory, the interferogram could be acquired by FTIS from the spectrum of the detected target [20]. The interference fringes of the interferogram can be idealized as:
I ( Δ ) = + B ( σ ) e j 2 π σ Δ d σ
where B is the spectrum, I is the interferogram, Δ is the path difference, and σ is the wavenumber. The reconstructed spectrum of the detected target is defined as:
B ( σ ) = + I ( Δ ) e j 2 π σ Δ d Δ
Equations (1) and (2) are the main relationship between the interferogram and the spectrum, and the fundamental theory of Fourier transform imaging spectrometers (FTISs).
For actual FTISs, the spectral range and the maximum optical path difference (MOPD) are limited [16], so the true relationship of the interferogram and the reconstructed spectrum are listed as follows:
I ( Δ ) = σ 1 σ 2 B ( σ ) e j 2 π σ Δ d σ = + T B ( σ ) B ( σ ) e j 2 π σ Δ d σ
B ( σ ) = L + L I ( Δ ) e j 2 π σ Δ d Δ = + T I ( Δ ) I ( Δ ) e j 2 π σ Δ d Δ
where σ 1 and σ 2 represent the range of wavenumbers and L is the maximum optical path difference (MOPD), which determines the FWHM of the distinguishable wavenumber of the instrument [9]. The FWHM of the distinguishable wavenumber is defined as:
δ σ = 1 / 2 L
The spectral resolution of the instrument is given by:
R = σ / δ σ
T I is the truncation function, and the Fourier transform of the truncation function is the instrument line shape. Since the nature of the Fourier transform determines the Gibbs phenomenon of the instrument linear function, the apodization function is applied to suppress the side lobes of Gibbs oscillation. Figure 2 shows the reconstructed spectrum of the laser, the blue line is the spectrum recovered by the original FFT, and the red line is the spectrum reconstructed using triangle apodization and FFT. It is clearly shown that triangle apodization can greatly suppress side lobes; however, the apodization method also brings a decrease in spectral resolution. The traditional SpecR processing of the interferogram [16] mainly includes preprocessing, phase correction, apodization, and FFT. Therefore, the spectral resolution of the hyperspectral data reconstructed by traditional SpecR methods has poor performance.

2.2. Relationship of FTIS Data Processing and Deep Learning Method

The FFT in traditional FTIS data processing can be viewed as a simplified, fully connected, single-layer network, as shown in Figure 2. We considered an interference curve as a vector with a length of l*1. In the FFT-based method, every point of the l*1 vector multiplied a coefficient on the corresponding band and then summed to obtain the spectral information of the band. As described in Equation (4), we can obtain spectral information of different bands by changing σ . As illustrated in Figure 3a, the FFT solver can be viewed as a simple linear mapping operator and was applied to transform the interference curve to the spectral curve (length is s*1). Figure 3b shows that the multilayer FCN and FFT transform have the same connection mode, but there was a nonlinear activation function in the multilayer FCN. It is possible to add nonlinear factors to increase the expression ability of the model. Significantly, his can acquire compelling results in numerous studies [25,34].
The above discussion shows that the FFT-based SpecR method can be viewed as a kind of FCN (without nonlinear mapping). According to Equation (4), the linear structure will produce side lobes with a limited optical path difference (OPD) [6]. In contrast, the multilayer FCN can be viewed as a nonlinear mapping with a stronger fitting ability than the linear model, and may produce a high-resolution spectrum without side lobes by the nonlinear transform. Motivated by the analogy, the deep learning method (fully connected U-Net) has been introduced to reconstruct the spectrum of FTIS.

2.3. The Architecture of a Typical U-Net

A U-Net is a U-shaped network structure which is widely used in the fields of image segmentation, image classification, image restoration, etc. The U-Net can extract multiscale features by encoder/decoder operating and can fuse features well by a concatenating procedure. The diagram (Figure 4a) of the original U-Net and an example of biomedical image segmentation (Figure 4b) are from Ronneberger’s work [27]. As shown in Figure 4b, the U-Net can identify the small separation borders of the adjacent touching cells and the precise contouring of cells, which indicates that the network can extract multiscale information and can fuse the extracted information well. Therefore, we used this structure for spectrum reconstruction (SpecR).
The process of the original U-Net for biomedical image segmentation is described in pseudocode as shown in Algorithm 1.
Algorithm 1: biomedical image segmentation network U-Net
Input: the pairs of training data, biomedical images Ib and the corresponding segmentation maps Ms
Output: the segmentation Ms
Step 1: load the training data and test data pairs (Ib, Ms)
Step 2: perform initialization of the network model, and randomly generate parameters θu of the U-Net
Step 3: training processing
For epoch in max number of Epochs
  Randomly load (Ib (n), Ms (n)) (n = [1, 2, …, N]), N is the batch size
  Calculate Ms’, and use the cross entropy to evaluate the consistency of segmentation results and true values
  Use gradient descent method to update the parameters θu, and save θu with the best segmentation result of training data
End
Step 4: the trained parameters θu are used to get the segmentation results of test datasets

3. Proposed Fully Connected U-Net for Spectrum Reconstruction

Our goal is to build a deep learning model that transforms interference curve I into spectral curve F(I), which is as similar as possible to the ground truth high spectral resolution spectrum B. In this work, we chose U-Net as the backbone network because the architecture contains the encoder/decoder paths that can extract multiscale features [27], and the concatenated structure of U-Net can fuse features well [33]. These features, with sufficient information, were the basis for accurate high-resolution spectrum reconstruction.

3.1. The Basic Architecture

Figure 5 and Table 1 show the schematic diagram and detailed parameters of the FCUN for SpecR. First, two individual fully connected (FC) layers were applied to extract the detailed features from the interference input. As a result, at each layer, the number of extracted features was expanded, which could be viewed as an increase in spectral resolution by extending the length of the interference curve. Then, six fully connected (FC) downsampling layers were used to calculate high-level features in growing numbers on a coarser scale. These features were connected with computed high-resolution features later using upsampling processing. To overcome the problem of gradient explosion and vanishing in the training processing, the residual connection was adopted to fuse upsampling/downsampling features. This type of residual connection in the U-Net structure can best exploit the potential of pairwise operations of different lengths by increasing the number of layers of upsampling/downsampling. There were six levels connected in the proposed network in this paper, and the length of the operational features of each level was half and double that of the previous level in the downsampling and upsampling layers, respectively. In the final layer, a linear fully connected structure was applied to obtain the spectral result. For interference curve of length J as input and K bands of spectrum to be estimated, the model returned predictions as K*1 vectors.
Table 2 shows the detailed structure of FC layers. The linear (in_len, out_len) indicates a linear transform to the input vector. In this paper, in_len indicates the size of the input sample and out_len stands for the size of the output sample. In the FC upsampling/downsampling layers and individual FC layer, the linear transform was followed by ReLU activation. It is possible to add more fully connected layers and ReLU activations to improve the nonlinearity. However, this can increase the complexity of the model, and more training time and data are needed. Hence, we used a single fully connected layer in the upsampling/downsampling layers because it has already produced compelling results.

3.2. Loss Function

Learning the end-to-end mapping function F(I) for SpecR requires the estimation of the parameters of the FCUN. This is accomplished by minimizing the loss between the reconstructed spectral curves and the corresponding ground truth spectral curves. Here, the parameters are denoted by θ, the reconstructed spectrum was defined as F(I; θ) or B, and the ground truth spectral curves are represented by B. Note that a good setting of the loss function is the key point to solving the parameter estimation. In general, the spectral angle (SA) [3] and the relative spectral quadratic error (RQE) [4] are two efficient indicators for quantitative evaluation of the error of the reconstructed spectrum. They are defined by:
S A ( B ( n ) , B ( n ) ) = cos 1 ( k = 1 K B k ( n ) B k ( n ) / ( k = 1 K ( B k ( n ) ) 2 k = 1 K ( B k ( n ) ) 2 ) )
R Q E ( B ( n ) , B ( n ) ) = k = 1 K ( B k ( n ) B k ( n ) ) 2 / k = 1 K B k ( n )
where B ( n ) stands for the ground truth of the spectrum, B ( n ) is the reconstructed spectrum, and K is the number of spectral bands. Then, the loss functions based on the above indicators are defined as:
l S A = 1 N n = 1 N S A ( B ( n ) , B ( n ) )
l R Q E = n = 1 N R Q E ( B ( n ) , B ( n ) )
where N is the total number of examples for the SpecR model training. Accordingly, we defined the final loss function as a linear combination of the two above components in Equations (9) and (10).
l F = l S A + τ l R E Q
Using l F above as the final loss function favors a high degree of consistency of shape and absolute value between the reconstructed spectrum and ground truth. Here, the value of τ was set as 0.5.

3.3. The Process of Generating SpecR Network FCUN

Our SpecR method based on FCUN is realized by training process. First, randomly select the training data and test data from the dataset; the paired interference curves and the corresponding spectral curves are used for training. Second, perform initialization of the network model by setting hyperparameters; the specific values are shown in Section 4.3, and randomly generate parameters θ of the FCUN. Third, train the network and update the parameters θ; after several iterations of parameter updates, the parameter θ is finally obtained.
The process of generating a SpecR network FCUN is illustrated in pseudocode as shown in Algorithm 2.
Algorithm 2: Spectrum reconstruction network FCUN
Input: the pairs of training data, interference curves and the corresponding spectral curves (I, B)
Output: the reconstructed spectral curves described as F(I; θ) or B
Step 1: load the training data and test data pairs (I, B)
Step 2: perform initialization of the network model, and randomly generate parameters θ of the FCUN
Step 3: training processing
For epoch in max number of Epochs
  Randomly load (I(n), B(n)) (n = [1, 2, …, N]), N is the batch size
  Calculate B ,   and   use   the   function   l F to evaluate the loss
  Update parameters θ, and save θ with the best reconstruction result of training data
End
Step 4: the trained parameters θ are applied to reconstruct the spectral curves B’ of test datasets

4. Experiments and Results

4.1. Datasets

There is no dedicated dataset for spectrum reconstruction of FTIS, especially for testing deep learning methods. We evaluated our SpecR network FCUN on two spectrum reconstruction datasets that we built. Both datasets were divided into two parts: interference curves and spectral curves.
In the first dataset, we simulated spectrum reconstruction datasets according to Equation (3). The real collected spectral data and simulated pulse spectrum were used to generate the corresponding interference curves, making the data sample diverse. We produced 13.7% asymmetric interferograms (almost single-sided) for calculations, as asymmetric structures can retain more spectral information at a fixed detector length for interference [19]. The unit OPD that indicates the OPD of two adjacent pixels was set as 206.96 nanometers (nm), which is the same as the visible and near infrared (VNIR) hyperspectral imager (HSI), a kind of FTIS on Chinese HJ-2(A/B) satellites that launched on 27 September 2020, and the MOPD was 45,738.16 nm in the simulation. Note that simulated interferograms were ideally sampled with one sample exactly located, and there was no inconsistency of the detector response and phase error. The spectra ranged from 455.06 nm to 898.73 nm with a resolution of 3.18 nm and 202 bands in total. The simulated pulsed spectrum with 202 bands was generated by adding Gaussian noise on the full zero spectrum. The two examples of simulated data were shown in Figure 6 and Figure 7.
The second dataset was constructed using current operational satellite images acquired by visible and near-infrared FTIS for testing the proposed FCUN for real data. The dataset included 640,000 samples of different locations on the Earth’s surface. One sample contained an interference curve where the length was 256 and a corresponding spectral curve with 100 bands equally spaced range from 450 nm to 920 nm. All interference curves were corrected by preprocessing and phase correction [19]. Considering tests of the spectral resolution of the reconstructed spectrum by different compared SpecR methods, the long side of corrected asymmetric interference curves was cut to different lengths for testing. Because different interference curve lengths indicate different MOPDs, the smaller MOPD resulted in lower spectral resolution according to Equations (5) and (6). If an algorithm can recover interference curves with a small MOPD into the high-resolution spectrum, the method can perform well in SpecR for FTISs. There is an example of real data shown in Figure 8.

4.2. Methods for Comparison

To test the proposed method, we compared the proposed FCUN with the most effective spectrum reconstruction (SpecR) method, the state-of-the-art FFT-based method [19], and the multitaper power spectral density estimate (MPSDE) method proposed by Thomson [35]. The PyTorch package was used for training and testing implementations [36].

4.3. Implementation Details

The proposed SpecR model was trained on a computer with an AMD Ryzen 9 4900HS CPU, 16 GB memory, and an RTX2060 MaxQ GPU with 6 GB GPU memory. During the training, the batch size was set as 2048, and the Adam algorithm was used for parameter optimization. The parameters for the learning rate, β1, β2, and epsilon were 0.001, 0.9, 0.999, and 10−9, respectively. The training and test data were normalized to avoid the problem of overfitting and further strengthen the training, and dropout was set to 0.5 in the training process. Our network was implemented on the PyTorch framework using the GPU, with 2000 training cycles on the first dataset and 1000 training cycles on the second dataset, because the training cycles were enough to produce a stable network. Finally, we saved the model with minimum training loss. The detailed parameters are listed in the Table 3. In our experiments, the length of the input interference curve for the first dataset was 256, and the lengths for the second data set were 200 and 150, respectively; therefore, we trained three FCUN models for spectrum reconstruction.

4.4. Evaluation

The spectral angle (SA) [3] in Equation (7), the relative spectral quadratic error (RQE) [4] in Equation (8), and peak signal to noise ratio (PSNR) [25] were used to evaluate the spectral shape, magnitude error, and similarity between the reconstructed spectrum and the ground truth of the two test datasets, respectively.

4.4.1. Evaluation of the First Dataset

In the experiments of the first dataset, 640,000 paired examples of interference curves and spectral curves were used for training and 480,000 paired examples were used for testing. The results of SA, RQE, and PSNR of the first test dataset, summarized in Table 4, indicated that the proposed FCUN improved the accuracy of the reconstructed spectrum obtained by simulation. Note that the proposed FCUN yielded the lowest mean SA score in the experiments, which indicated that the shape of the reconstructed spectrum was more consistent with the reference value. The proposed FCUN also gained the best RQE score, and thus, the absolute value of the reconstructed spectra showed better agreement with the ideal reference spectra. In addition, the FCUN achieved the highest mean PSNR score, suggesting that the reconstructed results have relatively better perceptual quality. The specific reconstructed examples from the test dataset are shown in Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13.
Figure 9a showed the typical recovered spectra of the wall. It is evident from the plot that the spectrum recovered by FFT and MPSDE had broader oxygen absorption peaks at 761.9 nm than FCUN. It can be proven that the FCUN achieved better performance in the spectral resolution. In addition, the relative errors [37] (Figure 9b) of spectra reconstructed by FCUN, FFT, and MPSDE to the reference spectrum were 0.65%, 5.60%, and 37.47%, respectively. Therefore, the FCUN had a greater potential to obtain accurate intensity values for spectral recovery.
The recovered spectrum of the tree was shown in Figure 10a, and the intensity value of the proposed FCUN nearly overlapped with the reference spectrum with a relative error of 1.24% (Figure 10b). The FFT produced a spectrum with little deviation from the actual reference with a relative error of 12.07% (Figure 10b), and the intensity value of the spectrum recovered by the MPSDE was quite different from the reference spectrum with a relative error of 64.34% (Figure 10b). Thus, the proposed FCUN also had a better performance for SpecR of the tree.
Pulse spectrum reconstruction was also a vital part for measuring the goodness of the algorithms. Figure 11 shows that the peak positions of the spectra recovered by all compared methods were consistent with the reference values. This result indicated that the three methods can perform well in wavelength accuracy. The spectra reconstructed by FCUN (the purple dash-dot line) and MPSD (the orange line) have few side lobes, while the spectrum of FFT (the red line) has side lobes that are approximately 3% of the high peak. This meant that the FFT-based SpecR method did not perform well in terms of the spectral accuracy of the relative value. Noted that the recovered spectral FWHMs of FCUN, MPSDE, and FFT were 2.22 nm, 4.00 nm, and 5.33 nm, respectively. A smaller FWHM value indicated a higher spectral resolution. Hence, the FCUN and MPSDE performed relatively well in terms of spectral resolution and spectral accuracy for one pulse spectrum.
Figure 12 shows the reconstructed spectra of multiple pulses with scattered peaks, which were obtained by the FFT, MPSDE, and FCUN. All methods accomplished the task of spectral reconstruction of the multipulse interference curve in their own way and obtained precise spectral locations of three peaks. As shown by the red line, the reconstructed spectrum obtained through FFT has more side lobes than the other two methods, which indicated that the intensity value of each band would be influenced by the adjacent bands. This would affect the accuracy of subsequent hyperspectral quantification analysis. As shown by the orange line, even though the spectrum recovered by MPSDE had fewer side lobes than the FFT-based method, the intensity value at 719.9 nm was nearly one-third of the spectrum of FFT (the red line), which would alter the spectral properties of intensity. Compared with FFT and MPSDE, the proposed FCUN had a better performance in spectral resolution, because the spectrum of the FCUN (the purple dash-dot line) with a very narrow FWHM is more consistent with the reference spectrum (the blue line) than the other two methods.
The reconstructed results of the multipulse spectrum with very close pulse peaks using the three compared methods are displayed in Figure 13. The FFT and MPSDE can only reconstruct two peaks in the blue dotted box, while the reference spectrum had three closed peaks, which indicated that the resolving power of FFT and MPSDE was not enough to distinguish very close peaks in the spectrum. The recovered spectrum of the FCUN had three peaks in the blue dotted box with accurate pulse positions. Although the intensity of the middle peak in the 642.7 nm band was different from the reference spectrum, FCUN still performed better than FFT and MPSDE.
Based on the comparison of the reconstructed pulse spectra of Figure 11, Figure 12, Figure 13, it was clear that the FCUN achieved better performance than the other two methods.

4.4.2. Evaluation of the Second Dataset

To assess the effectiveness of the proposed method on real data, 440,000 samples were extracted from the second dataset as training samples and 200,000 as test samples, with interference curves of lengths 200 and 150, respectively. The SA and RQE of the reconstructed results of the test samples using real interference curves with different data lengths were shown in Table 5. Because the MPSDE does not perform well in the continuous spectrum as shown in Table 4, we only compared the state-of-the-art FFT and the proposed FCUN. From the table, the spectra reconstructed by FCUN had a lower value in both quantitative evaluation indicators than the FFT with the same data length. Thus, the FCUN had better performance than the FFT in the real SpecR of FTIS. Additionally, with the decrease of the data length ranges from 200 to 150, the SA of FFT increased 8.6%, the RQE of FFT increased 16.9%, and the SA and RQE of FCUN also increased 15.1% and 18.4%, respectively. This also indicated that a larger length of the interference curve can produce a better spectrum. The specific reconstructed examples from the test dataset are shown in Figure 14 and Figure 15.
The reconstructed spectra of nudation in the test dataset and the relative errors to the reference spectrum are shown in Figure 14. The proposed FCUN method provided recovered spectra that were more consistent with the ideal reference spectrum than the spectra produced by FFT (Figure 14a). The consistency rankings of details (the red rectangle box in Figure 14a) with the reference spectrum were the FCUN with a data length of 200, the FCUN with a data length of 150, the FFT with a data length of 200, and the FFT with a data length of 150. In addition, the relative errors (Figure 14b) of the compared methods with different data lengths were consistent with the detailed comparison of Figure 14a. Moreover, the compared methods provided results of peaky spectrum in Figure 15a in agreement with the results in Figure 14a; the reconstructed spectra of FCUN were in closer agreement with the ideal reference spectrum with the minimal relative error (Figure 15b). Hence, the FCUN had better performance than the FFT on the real dataset, and data with longer interference length gave better results in both methods.

4.4.3. Evaluation of the Noise Data

Our proposed method also had good generalization performance on simulated noisy interferograms. Here, different signal-to-noise ratio (SNR) data were generated by adding different levels of Gaussian noise on the original interferograms according to the following Equation:
S N R = 20 log 10 ( m e a n ( I o r g ) / M S E )
where Iorg denotes the interference curve with DC component, MSE is the variance of the noise, and mean(*) denotes the mean value. According to the different SNRs, we calculated the corresponding MSE and superimposed Gaussian noise with mean value of 0 and variance value of MSE on interferograms.
Figure 16 shows the true color images of the reconstructed hyperspectral data cubes of the first dataset for different signal-to-noise ratio (SNR) [8] interferograms. The outputs of the FFT contained a large amount of noise even though the SNR is 50 dB (Figure 16b). While the FCUN exhibited superior visual performance when the SNR was 40 dB (Figure 16e), there was no obvious noise on the images (Figure 16d,e). Moreover, the details of the sky in the red rectangular box are in the lower left corner of each subfigure, and compared to the FFT method, the proposed FCUN method had less noise in the subimage and the sky was more uniform. Most remarkably, the mean relative error (MRE) [37] of FFT with SNRs of 40 dB and 50 dB were 31.97% and 11.93%, respectively, while the MRE of FCUN with SNRs of 40 dB and 50 dB were 5.28% and 2.68% (Figure 17). The FCUN gained nearly 80% relative improvement over FFT for MRE. This gain is due to the learned representations of the FCUN.
Figure 18 shows the true color building images of the reconstructed hyperspectral data cubes of the first dataset for different SNR interferograms. From Figure 18d, the outputs of the FCUN had a similar visual quality to that of the reference image (Figure 18a). While the visual quality of the recovered results of FFT was poor, some details of the area corresponding to the red rectangular box were lost (Figure 18b,c). As shown in Figure 19, the proposed approach had the lowest relative error for the interferograms with an SNR of 50. The FFT had higher relative errors than the FCUN. This indicated that the nonlinear method had better performance in SpecR than the linear approach.
In the actual imaging system of FTIS, the data will be corrupted by Gaussian (readout) noise and Poisson (shot) noise. In the following experiment, we accounted for both noises. The original data were generated by adding Poisson and Gaussian noise according to light conditions and characteristics of the FTIS on Chinese HJ-2(A/B) satellites. For the FTIS, each photon received by the detector produces a voltage of 4.2 millivolt; under the default imaging conditions, the system saturation digital number (DN) value was 4095 when the system voltage was 2V. So for the FTIS system, one DN value corresponds to 116 photons. In this experiment, we simulated Poisson noise for three different luminance conditions (low luminance, medium luminance, and high luminance). The average DN of the interferograms are about 500, 1000, and 2000, corresponding to the photon numbers of 58000, 116000, and 232000, respectively. The simulation of Gaussian noise was based on variance of dark current of the FTIS, which is equivalent to 97 photons.
Figure 20, Figure 21 and Figure 22 show the true color images of the reconstructed hyperspectral data cubes for the three different luminances. From the three figures, the outputs of the FCUN have a similar visual quality to that of the reference image. The reconstructed results of the FFT had noise-induced color distortion, and some details were obscured by noise. As shown in Figure 23, the relative error of FCUN is smaller than FFT, the corresponding mean relative errors of FFT and FCUN were 3.33% and 0.87%, and the FCUN gained nearly 73.9% relative improvement than the FFT for SpecR of mixed noise (Gaussian and Poisson) data.
Through quantitative and qualitative comparisons, the proposed FCUN can accurately reconstruct the spectral information of Fourier transform imaging spectrometers (FTISs) and has better stability than the compared methods.

4.4.4. Evaluation of the Structure of FCUN

We then examined the sensitivity of the network to different layers and structures. In previous experiments, we set the proposed FCUN with two individual FC layers, six FC downsampling/upsampling layers, and one linear transform layer. Here, we change the number of layers of the FCUN to obtain two other SpecR networks with different sizes. First, we cut the two individual FC layers, and only use four downsampling/upsampling layers for SpecR as a lightweight FCUN. Second, we cut FC downsampling/upsampling layers from six to three of the proposed FCUN to evaluate the performance of layers in U-Net. All the experiments were done using the first dataset, and the hyperparameters were consistent with those in Table 3. The quantitative evaluation of the different structures for testing the dataset was listed in Table 6.
The ranking of reconstruction accuracy of the compared structures is: the proposed FCUN, the FCUN with 3 layers cut, and lightweight FCUN, in order. This suggested that a reasonably larger size of the FCUN could extract richer spectral information, which would lead to better reconstructed results of FTIS. It is important to emphasize that the feature vector upsampled in the beginning for two layers is important for quality improvement of SpecR. However, the size of the model’s parameters and deployment speed also increases with a larger network size. Therefore, the choice of the network parameters should always be a compromise between speed and performance based on the constraint of computing power.

5. Discussion

Some recent research has shown that deep learning techniques can improve the performance of FTISs in many fields [32,33,34] because of the ability of feature extraction and nonlinear transformation. However, as far as we know, the models of deep learning have not been used in the spectrum reconstruction (SpecR) of FTISs. In this study, we first developed a nonlinear model based on deep learning for spectrum reconstruction (SpecR). The above experiments supported our hypothesis that the nonlinear deep learning model can achieve good performance in spectral resolution (Figure 11) and intensity value (Figure 8 and Figure 9) for SpecR. This accords with previous studies that deep learning can improve the performance of FTISs. The SA and RQE of the proposed FCUN are lower than the state-of-the-art methods, and the FCUN can also achieve better visual quality of the reconstructed images (Figure 16, Figure 18, Figure 20, Figure 21, Figure 22). A possible explanation for the better performance of the FCUN might be that the U-Net structure with residual connection is suitabl for extracting and fusing multiscale features of interference curve. Compared with the other network architectures, the FCUN also worked better in SpecR, which proves the effectiveness of the proposed FCUN model to some extent.
The reconstructed result of the multipulse spectrum with very close pulse peaks may be unsatisfactory because of the limited spectral resolution of FCUN. It is possible that the training data is not abundant enough or the structure of FCUN needs further improvement. It could be hypothesized that the performance of the deep learning method for SpecR can be further promoted by increasing the network layers and using a more efficient network structure, such as attention and transformer [38].
Despite promising results of the FCUN for most testing data, further work on analyzing network structure and building up abundant datasets for spectrum reconstruction should be conducted.

6. Conclusions

In this paper, we proposed a novel fully connected U-Net (FCUN) for spectrum reconstruction (SpecR) of Fourier transform imaging spectrometers (FTISs). We first showed that the traditional fast Fourier transform (FFT) SpecR methods can be reformulated into a deep learning algorithm. The presented approach, FCUN, learned an end-to-end mapping between interference curves and reconstructed spectrum, with little pre/post-processing before the model training, and outperformed the competing state-of-the-art methods for testing simulated and real datasets. It is also important to note that the proposed FCUN easily reconstructs a spectrum with higher resolution than the other algorithms, although it produced some distortion when there were many pulses close to each other in the spectrum. Moreover, the FCUN achieved excellent SpecR results for noisy interferograms. We speculate that better performance can be further obtained by exploring more effective structures of the network, larger training datasets, and different training strategies. In addition, our method has the advantage of stability and simplicity and could be further used for SpecR of a wide range of FTISs.

Author Contributions

This manuscript was written by T.C.; T.C., X.S. and H.L. presented the main idea of this work; T.C., G.Z. and X.F. wrote the program code and completed all experiments; S.W., J.L. and Y.W. revised the paper structure. C.Z. provided the simulated data. X.S., X.L. and S.L. supervised the work and contributed to the experiments and discussions. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Natural Science Foundation of Shaanxi Province (Grant NO.2020JQ-436), the Open Research Fund of National Earth Observation Data Center (Grant NO.NODAOP2021007), the National Science Basic Research Foundation of Shaanxi Province (Grant NO.2021JQ-324), the National Natural Science Foundation of China under Grants (Grant NO. 41706217, 42176182), the National Science Basic Research Foundation of Shaanxi Province (Grant NO.2019JQ-929), and Youth Innovation Promotion Association CAS.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the editors and all of the anonymous reviewers for their constructive comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xiangli, B.; Cai, Q.; Du, S. Large aperture spatial heterodyne imaging spectrometer: Principle and experimental results. Opt. Commun. 2015, 357, 148–155. [Google Scholar] [CrossRef]
  2. Wang, S.; Li, L.B.; Pi, H.F. Research of Spectrum Signal-to-Noise Ratio of Large Aperture Static Imaging Spectrometer. Spectrosc. Spectr. Anal. 2014, 3, 851–856. [Google Scholar]
  3. Su, L.; Yuan, Y.; Xiangli, B.; Huang, F.; Cao, J.; Li, L.; Zhou, S. Spectrum Reconstruction Method for Airborne Temporally–Spatially Modulated Fourier Transform Imaging Spectrometers. IEEE Trans. Geosci. Remote Sens. 2014, 52, 3720–3728. [Google Scholar] [CrossRef]
  4. Zhang, G.; Shi, D.; Wang, S.; Yu, T.; Hu, B. Data correction techniques for the airborne large-aperture static image spectrometer based on image registration. J. Appl. Remote Sens. 2015, 9, 095088. [Google Scholar] [CrossRef]
  5. Mueller, D.; Ferrão, M.F.; Marder, L.; Costa, A.D.; Rosana, D. Fourier Transform Infrared Spectroscopy (FTIR) and Multivariate Analysis for Identification of Different Vegetable Oils Used in Biodiesel Production. Sensors 2013, 13, 4258–4271. [Google Scholar] [CrossRef] [Green Version]
  6. Amenabar, I.; Poly, S.; Goikoetxea, M.; Nuansing, W.; Lasch, P.; Hillenbrand, R. Hyperspectral infrared nanoimaging of organic samples based on Fourier transform infrared nanospectroscopy. Nat. Commun. 2017, 8, 14402. [Google Scholar] [CrossRef] [Green Version]
  7. Zhang, C. Complete throughput analysis of the static four-phase-divided wind imaging interferometer. Opt. Eng. 2012, 51, 046202. [Google Scholar] [CrossRef]
  8. Ferrec, Y. Noise sources in imaging static Fourier transform spectrometers. Opt. Eng. 2012, 51, 111716. [Google Scholar] [CrossRef]
  9. Huang, F.; Yuan, Y.; Li, J.; Cao, J. Spectrum reconstruction using relative-deviation-based kernel regression in temporally and spatially modulated Fourier transform imaging spectrometer. Appl. Opt. 2015, 54, 6737–6743. [Google Scholar] [CrossRef]
  10. Soncco, D.-C.; Barbanson, C.; Nikolova, M.; Almansa, A.; Ferrec, Y. Fast and Accurate Multiplicative Decomposition for Fringe Removal in Interferometric Images. IEEE Trans. Comput. Imaging 2017, 3, 187–201. [Google Scholar] [CrossRef] [Green Version]
  11. Ren, W.; Zhang, C.; Mu, T.; Dai, H. Spectrum reconstruction based on the constrained optimal linear inverse methods. Opt. Lett. 2012, 37, 2580–2582. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, G.; Wang, S.; Li, L.; Hu, X.; Hu, B.; Wang, S. Interference data correction methods for lunar observation with a large-aperture static imaging spectrometer. Appl. Opt. 2016, 55, 8770. [Google Scholar] [CrossRef] [PubMed]
  13. Cao, J.; Yuan, Y.; Su, L.; Zhu, C.; Yan, Q. Online Correction Method for the Registration Error between TSMFTIS Detector and Interferogram. Sensors 2020, 20, 1195. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Cai, Q.; Xiangli, B.; Huang, M.; Han, W.; Pei, L.; Bu, M. Prototype development and field measurements of high etendue spatial heterodyne imaging spectrometer. Opt. Commun. 2018, 410, 403–409. [Google Scholar] [CrossRef]
  15. Mu, T.; Zhang, C.; Zhao, B. Analysis of a moderate resolution Fourier transform imaging spectrometer. Opt. Commun. 2009, 282, 1699–1705. [Google Scholar] [CrossRef]
  16. Zhang, W.; Wen, D.; Song, Z. Spectrum reconstruction in interference spectrometer based on sparse Fourier transform. Optik 2018, 154, 157–164. [Google Scholar] [CrossRef]
  17. Markowski, K.; Jdrzejewski, K.; Osuch, T. Optimization of group delay response of (apodized) tapered fiber Bragg grating by shaping taper transition and apodization window. Proc. SPIE Int. Soc. Opt. Eng. 2016, 100311, 100311G. [Google Scholar]
  18. Stankwitz, H.C.; Dallaire, R.J.; Fienup, J.R. Nonlinear apodization for sidelobe control in SAR imagery. IEEE Trans. Aerosp. Electron. Syst. 1995, 31, 267–279. [Google Scholar] [CrossRef]
  19. Zhao, X.; Lu, J.; Gong, A.; Xiao, Z.; Kang, Q.; Li, Q. Operational spectrum reconstruction of data from the Fourier transform hyperspectral imager onboard HJ-1A satellite. Chin. Sci. Bull. 2010, 55, 1808–1812. [Google Scholar] [CrossRef]
  20. Jian, X.; Zhang, C.; Zhao, B.; Zhu, B. The application of MUSIC algorithm in spectrum reconstruction and interferogram processing. Opt. Commun. 2008, 281, 2424–2428. [Google Scholar] [CrossRef]
  21. Zhang, W.; Wen, D.; Song, Z.; Wei, X.; Liu, G.; Li, Z. High Resolution and Fast Processing of Spectral Reconstruction in Fourier Transform Imaging Spectroscopy. Sensors 2018, 18, 4159. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Zhang, C.; Jian, X. Wide-spectrum reconstruction method for a birefringence interference imaging spectrometer. Opt. Lett. 2010, 35, 366–368. [Google Scholar] [CrossRef]
  23. Panboonyuen, T.; Jitkajornwanich, K.; Lawawirojwong, S.; Srestasathiern, P.; Vateekul, P. Transformer-Based Decoder Designs for Semantic Segmentation on Remotely Sensed Images. Remote Sens. 2021, 13, 5100. [Google Scholar] [CrossRef]
  24. Kai, Z.; Zuo, W.; Gu, S.; Lei, Z. Learning Deep CNN Denoiser Prior for Image Restoration. In Proceedings of the 2017 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Honolulu, HI, USA, 21–26 July 2017. [Google Scholar]
  25. Chao, D.; Chen, C.L.; He, K.; Tang, X. Learning a Deep Convolutional Network for Image Super-Resolution. In Proceedings of the Computer Vision—ECCV 2014: 13th European Conference, Zurich, Switzerland, 6–12 September 2014. [Google Scholar]
  26. Yu, B.; Chen, F.; Wang, Y.; Wang, N.; Yang, X.; Ma, P.; Zhou, C.; Zhang, Y. Res2-Unet+, a Practical Oil Tank Detection Network for Large-Scale High Spatial Resolution Images. Remote Sens. 2021, 13, 4740. [Google Scholar] [CrossRef]
  27. Ronneberger, O.; Fischer, P.; Brox, T. U-Net: Convolutional Networks for Biomedical Image Segmentation. In Medical Image Computing and Computer-Assisted Intervention—MICCAI 2015; Springer: Cham, Switzerland, 2015. [Google Scholar]
  28. Zhou, Z.; Crawford, M.M. A Batch-Mode Regularized Multimetric Active Learning Framework for Classification of Hyperspectral Images. IEEE Trans. Geosci. Remote Sens. 2017, 55, 6594–6609. [Google Scholar]
  29. Torres, D.L.; Turnes, J.N.; Soto Vega, P.J.; Feitosa, R.Q.; Silva, D.E.; Marcato Junior, J.; Almeida, C. Deforestation Detection with Fully Convolutional Networks in the Amazon Forest from Landsat-8 and Sentinel-2 Images. Remote Sens. 2021, 13, 5084. [Google Scholar] [CrossRef]
  30. Berisha, S.; Lotfollahi, M.; Jahanipour, J.; Gurcan, I.; Walsh, M.; Bhargava, R.; Van Nguyen, H.; Mayerich, D. Deep learning for FTIR histology: Leveraging spatial and spectral features with convolutional neural networks. Analyst 2019, 144, 1642–1653. [Google Scholar] [CrossRef]
  31. Lotfollahi, M.; Berisha, S.; Daeinejad, D.; Mayerich, D. Digital Staining of High-Definition Fourier Transform Infrared (FT-IR) Images Using Deep Learning. Appl. Spectrosc. 2019, 73, 556–564. [Google Scholar] [CrossRef]
  32. Yue, F.; Chen, C.; Yan, Z.; Chen, C.; Guo, Z.; Zhang, Z.; Chen, Z.; Zhang, F.; Lv, X. Fourier transform infrared spectroscopy combined with deep learning and data enhancement for quick diagnosis of abnormal thyroid function. Photodiagnosis Photodyn. Ther. 2020, 32, 101923. [Google Scholar] [CrossRef]
  33. Han, Y.; Ye, J.C. Framing U-Net via Deep Convolutional Framelets: Application to Sparse-View CT. IEEE Trans. Med. Imaging 2018, 37, 1418–1429. [Google Scholar] [CrossRef] [Green Version]
  34. Yang, Z.; Diao, C.; Li, B. A Robust Hybrid Deep Learning Model for Spatiotemporal Image Fusion. Remote Sens. 2021, 13, 5005. [Google Scholar] [CrossRef]
  35. Percival, D.B.; Walden, A.T. Spectral Analysis for Physical Applications; Multitaper and Conventional Univariate Techniques; Multitaper Spectral Estimation; Cambridge University Press: London, UK, 1993; pp. 331–377. [Google Scholar] [CrossRef]
  36. Rauber, J.; Brendel, W.; Bethge, M. Foolbox v0.8.0: A Python toolbox to benchmark the robustness of machine learning models. arXiv 2017, arXiv:1707.04131. [Google Scholar]
  37. Moustakides, G.V.; Poor, H.V. On the relative error probabilities of linear multiuser detectors. IEEE Trans. Inf. Theory 2001, 47, 450–456. [Google Scholar] [CrossRef]
  38. Li, T.; Song, H.; Zhang, K.; Liu, Q. Recurrent Reverse Attention Guided Residual Learning for Saliency Object Detection. Neurocomputing 2020, 389, 170–178. [Google Scholar] [CrossRef]
Figure 1. The optical layout of LASIS.
Figure 1. The optical layout of LASIS.
Remotesensing 14 00900 g001
Figure 2. Reconstructed spectrum using laser. The blue line is the reconstructed spectrum of original interference curve, the red line is the reconstructed spectrum of interference curve by apodization processing.
Figure 2. Reconstructed spectrum using laser. The blue line is the reconstructed spectrum of original interference curve, the red line is the reconstructed spectrum of interference curve by apodization processing.
Remotesensing 14 00900 g002
Figure 3. An illustration of the relationship of FFT and deep learning for FTIS data processing: (a) The relationship of FFT and linear mapping; (b) The relationship of multilayer FCN and FFT transform.
Figure 3. An illustration of the relationship of FFT and deep learning for FTIS data processing: (a) The relationship of FFT and linear mapping; (b) The relationship of multilayer FCN and FFT transform.
Remotesensing 14 00900 g003
Figure 4. The diagram of original U-Net and segmentation result [27]: (a) The diagram of the original U-Net for biomedical image segmentation; (b) An example of segmentation results of touching cells.
Figure 4. The diagram of original U-Net and segmentation result [27]: (a) The diagram of the original U-Net for biomedical image segmentation; (b) An example of segmentation results of touching cells.
Remotesensing 14 00900 g004
Figure 5. The architecture of fully connected U-Net (FCUN) for spectrum reconstruction (SpecR). The green arrow is the fully connected downsampling layer, the red arrow is the fully connected upsampling layer, the gray arrow is the linear transform, the residual connection is indicated by gray arrows and plus signs. The architecture of FCUN contains six residual connections.
Figure 5. The architecture of fully connected U-Net (FCUN) for spectrum reconstruction (SpecR). The green arrow is the fully connected downsampling layer, the red arrow is the fully connected upsampling layer, the gray arrow is the linear transform, the residual connection is indicated by gray arrows and plus signs. The architecture of FCUN contains six residual connections.
Remotesensing 14 00900 g005
Figure 6. Simulated data of sky: (a) Original spectrum; (b) Simulated interference curve.
Figure 6. Simulated data of sky: (a) Original spectrum; (b) Simulated interference curve.
Remotesensing 14 00900 g006
Figure 7. Simulated data of pulse spectroscopy: (a) Original spectrum; (b) Simulated interference curve.
Figure 7. Simulated data of pulse spectroscopy: (a) Original spectrum; (b) Simulated interference curve.
Remotesensing 14 00900 g007
Figure 8. Real interference curve and spectrum: (a) Real interference curve; (b) The corresponding spectrum.
Figure 8. Real interference curve and spectrum: (a) Real interference curve; (b) The corresponding spectrum.
Remotesensing 14 00900 g008
Figure 9. (a) Reconstructed spectrum of wall by FFT, MPSDE, and FCUN; (b) Relative error to the reference spectrum. The red line is the result of FFT, the orange line is the result of MPSDE, and the purple dash-dot line is the result of the proposed FCUN.
Figure 9. (a) Reconstructed spectrum of wall by FFT, MPSDE, and FCUN; (b) Relative error to the reference spectrum. The red line is the result of FFT, the orange line is the result of MPSDE, and the purple dash-dot line is the result of the proposed FCUN.
Remotesensing 14 00900 g009
Figure 10. (a) Reconstructed spectrum of tree by FFT, MPSDE and FCUN; (b) Relative error to the reference spectrum. The red line is the result of FFT, the orange line is the result of MPSDE, and the purple dash-dot line is the result of the proposed FCUN.
Figure 10. (a) Reconstructed spectrum of tree by FFT, MPSDE and FCUN; (b) Relative error to the reference spectrum. The red line is the result of FFT, the orange line is the result of MPSDE, and the purple dash-dot line is the result of the proposed FCUN.
Remotesensing 14 00900 g010
Figure 11. Single pulse spectrum reconstruction using different methods. The blue line is the reference spectrum, the red line is the spectrum reconstructed by FFT, the orange line is the spectrum reconstructed by MPSDE, and the purple dash-dot line is the spectrum reconstructed by FCUN.
Figure 11. Single pulse spectrum reconstruction using different methods. The blue line is the reference spectrum, the red line is the spectrum reconstructed by FFT, the orange line is the spectrum reconstructed by MPSDE, and the purple dash-dot line is the spectrum reconstructed by FCUN.
Remotesensing 14 00900 g011
Figure 12. Multipulse spectrum reconstruction using different methods. The blue line is the reference spectrum, the red line is the spectrum reconstructed by FFT, the orange line is the spectrum reconstructed by MPSDE, and the purple dash-dot line is the spectrum reconstructed by FCUN.
Figure 12. Multipulse spectrum reconstruction using different methods. The blue line is the reference spectrum, the red line is the spectrum reconstructed by FFT, the orange line is the spectrum reconstructed by MPSDE, and the purple dash-dot line is the spectrum reconstructed by FCUN.
Remotesensing 14 00900 g012
Figure 13. Reconstructed results of the multipulse spectrum with very close pulse peaks using different methods. The blue line is the reference spectrum, the red line is the spectrum reconstructed by FFT, the orange line is the spectrum reconstructed by MPSDE, and the purple dash-dot line is the spectrum reconstructed by FCUN.
Figure 13. Reconstructed results of the multipulse spectrum with very close pulse peaks using different methods. The blue line is the reference spectrum, the red line is the spectrum reconstructed by FFT, the orange line is the spectrum reconstructed by MPSDE, and the purple dash-dot line is the spectrum reconstructed by FCUN.
Remotesensing 14 00900 g013
Figure 14. Reconstructed results and relative errors of nudation: (a) Reconstructed spectrum of nudation; (b) Relative errors to the reference spectrum. The blue line is the reference spectrum, the red line is the result of FFT using the interference curve of length 150, the orange line is the result of MPSDE using the interference curve of length 200, the purple dash-dot line is the result of FCUN using the interference curve of length 150, and the light green dash line is the result of FCUN using the interference curve of length 200.
Figure 14. Reconstructed results and relative errors of nudation: (a) Reconstructed spectrum of nudation; (b) Relative errors to the reference spectrum. The blue line is the reference spectrum, the red line is the result of FFT using the interference curve of length 150, the orange line is the result of MPSDE using the interference curve of length 200, the purple dash-dot line is the result of FCUN using the interference curve of length 150, and the light green dash line is the result of FCUN using the interference curve of length 200.
Remotesensing 14 00900 g014
Figure 15. Reconstructed results and relative errors of peaky spectrum: (a) Reconstructed spectrum of peaky spectrum; (b) Relative error to the reference spectrum. The blue line is the reference spectrum, the red line is the result of FFT using the interference curve of length 150, the orange line is the result of MPSDE using the interference curve of length 200, the purple dash-dot line is the result of FCUN using the interference curve of length 150, and the light green dash line is the result of FCUN using the interference curve of length 200.
Figure 15. Reconstructed results and relative errors of peaky spectrum: (a) Reconstructed spectrum of peaky spectrum; (b) Relative error to the reference spectrum. The blue line is the reference spectrum, the red line is the result of FFT using the interference curve of length 150, the orange line is the result of MPSDE using the interference curve of length 200, the purple dash-dot line is the result of FCUN using the interference curve of length 150, and the light green dash line is the result of FCUN using the interference curve of length 200.
Remotesensing 14 00900 g015
Figure 16. Reconstruction results of buildings and trees with different signal-to-noise ratio (SNR) data: (a) Reference data; (b) FFT (SNR is 50 dB); (c) FFT (SNR is 40 dB); (d) FCUN (SNR is 50 dB); (e) FCUN (SNR is 40 dB).
Figure 16. Reconstruction results of buildings and trees with different signal-to-noise ratio (SNR) data: (a) Reference data; (b) FFT (SNR is 50 dB); (c) FFT (SNR is 40 dB); (d) FCUN (SNR is 50 dB); (e) FCUN (SNR is 40 dB).
Remotesensing 14 00900 g016
Figure 17. Relative error of the reconstruction results for different SNRs of Figure 16.
Figure 17. Relative error of the reconstruction results for different SNRs of Figure 16.
Remotesensing 14 00900 g017
Figure 18. Reconstruction results of buildings with different signal-to-noise ratio (SNR) data: (a) Reference data; (b) FFT (SNR is 50 dB); (c) FFT (SNR is 40 dB); (d) FCUN (SNR is 50 dB); (e) FCUN (SNR is 40 dB).
Figure 18. Reconstruction results of buildings with different signal-to-noise ratio (SNR) data: (a) Reference data; (b) FFT (SNR is 50 dB); (c) FFT (SNR is 40 dB); (d) FCUN (SNR is 50 dB); (e) FCUN (SNR is 40 dB).
Remotesensing 14 00900 g018
Figure 19. Relative error of the reconstruction results for different SNRs of Figure 18.
Figure 19. Relative error of the reconstruction results for different SNRs of Figure 18.
Remotesensing 14 00900 g019
Figure 20. Reconstruction results for high luminance data with Poisson and Gaussian noise: (a) Reference data; (b) Result of the FFT method; (c) Result of the proposed FCUN.
Figure 20. Reconstruction results for high luminance data with Poisson and Gaussian noise: (a) Reference data; (b) Result of the FFT method; (c) Result of the proposed FCUN.
Remotesensing 14 00900 g020
Figure 21. Reconstruction results for medium luminance data with Poisson and Gaussian noise: (a) Reference data; (b) Result of the FFT method; (c) Result of the proposed FCUN.
Figure 21. Reconstruction results for medium luminance data with Poisson and Gaussian noise: (a) Reference data; (b) Result of the FFT method; (c) Result of the proposed FCUN.
Remotesensing 14 00900 g021
Figure 22. Reconstruction results for low luminance data with Poisson and Gaussian noise: (a) Reference data; (b) Result of the FFT method; (c) Result of the proposed FCUN.
Figure 22. Reconstruction results for low luminance data with Poisson and Gaussian noise: (a) Reference data; (b) Result of the FFT method; (c) Result of the proposed FCUN.
Remotesensing 14 00900 g022
Figure 23. Relative error of the reconstruction results by adding Gaussian and Poisson noises. The red line is relative error of the FCUN, and the blue line is relative error of FFT.
Figure 23. Relative error of the reconstruction results by adding Gaussian and Poisson noises. The red line is relative error of the FCUN, and the blue line is relative error of FFT.
Remotesensing 14 00900 g023
Table 1. Detailed parameters of the FCUN.
Table 1. Detailed parameters of the FCUN.
LayerInput SizeOutput SizeLayerInput SizeOutput Size
Individual FC-1(J,1)(512,1)FC upsampling-1(16,1)(32,1)
Individual FC-2(512,1)(1024,1)FC upsampling-2(32,1)(64,1)
FC downsampling-1(1024,1)(512,1)FC upsampling-3(64,1)(128,1)
FC downsampling-2(512,1)(256,1)FC upsampling-4(128,1)(256,1)
FC downsampling-3(256,1)(128,1)FC upsampling-5(256,1)(512,1)
FC downsampling-4(128,1)(64,1)FC upsampling-6(512,1)(1024,1)
FC downsampling-5(64,1)(32,1)Linear transform(1024,1)(K,1)
FC downsampling-6(32,1)(16,1)
Table 2. Fully connected block diagram of the base architecture.
Table 2. Fully connected block diagram of the base architecture.
LayerOperationInputOutput
FC upsamplingLinear(n, 2n), Relu()n*1 vector(2n)*1 vector
FC downsamplingLinear(n, n/2), Relu()n*1 vector(n/2)*1 vector
Individual FC Linear(n, m), Relu()n*1 vectorm*1 vector
Table 3. Hyper-parameters of FCUN are used in this paper.
Table 3. Hyper-parameters of FCUN are used in this paper.
Parameter NameParameter Setting
Batch size2048
OptimizerAdam
Initial learning rate0.001
β10.9
β20.999
Epsilon10−9
Epochs2000/1000
Dropout0.5
Table 4. Quantitative evaluation of the tested methods.
Table 4. Quantitative evaluation of the tested methods.
AlgorithmSARQEPSNR (dB)
FFT0.06550.005629.61
MPSDE1.11690.272016.07
FCUN0.3227 × 10−60.0672 × 10−638.43
Table 5. Quantitative evaluation of the tested methods for interference curves of different data lengths.
Table 5. Quantitative evaluation of the tested methods for interference curves of different data lengths.
AlgorithmData LengthSARQEPSNR (dB)
FFT1500.15130.018024.19
2000.13930.015424.81
FCUN1500.1379 × 10−60.1189 × 10−629.61
2000.1198 × 10−60.1004 × 10−634.55
Table 6. Quantitative evaluation of the different network structure.
Table 6. Quantitative evaluation of the different network structure.
AlgorithmSARQEParameters Size (Mb)
The proposed FCUN0.3227 × 10−60.0672 × 10−68.65
Lightweight FCUN (The first)0.8527 × 10−60.3209 × 10−60.54
FCUN with 3 layers cut (The second)0.4298 × 10−60.0740 × 10−68.56
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, T.; Su, X.; Li, H.; Li, S.; Liu, J.; Zhang, G.; Feng, X.; Wang, S.; Liu, X.; Wang, Y.; et al. Learning a Fully Connected U-Net for Spectrum Reconstruction of Fourier Transform Imaging Spectrometers. Remote Sens. 2022, 14, 900. https://doi.org/10.3390/rs14040900

AMA Style

Chen T, Su X, Li H, Li S, Liu J, Zhang G, Feng X, Wang S, Liu X, Wang Y, et al. Learning a Fully Connected U-Net for Spectrum Reconstruction of Fourier Transform Imaging Spectrometers. Remote Sensing. 2022; 14(4):900. https://doi.org/10.3390/rs14040900

Chicago/Turabian Style

Chen, Tieqiao, Xiuqin Su, Haiwei Li, Siyuan Li, Jia Liu, Geng Zhang, Xiangpeng Feng, Shuang Wang, Xuebin Liu, Yihao Wang, and et al. 2022. "Learning a Fully Connected U-Net for Spectrum Reconstruction of Fourier Transform Imaging Spectrometers" Remote Sensing 14, no. 4: 900. https://doi.org/10.3390/rs14040900

APA Style

Chen, T., Su, X., Li, H., Li, S., Liu, J., Zhang, G., Feng, X., Wang, S., Liu, X., Wang, Y., & Zou, C. (2022). Learning a Fully Connected U-Net for Spectrum Reconstruction of Fourier Transform Imaging Spectrometers. Remote Sensing, 14(4), 900. https://doi.org/10.3390/rs14040900

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