Next Article in Journal
Insect Target Classes Discerned from Entomological Radar Data
Next Article in Special Issue
Multi-Label Remote Sensing Image Classification with Latent Semantic Dependencies
Previous Article in Journal
Phenological Characteristics of Global Ecosystems Based on Optical, Fluorescence, and Microwave Remote Sensing
Previous Article in Special Issue
Remote Sensing Image Ship Detection under Complex Sea Conditions Based on Deep Semantic Segmentation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of the Remapping Algorithms for the Advanced Technology Microwave Sounder (ATMS)

Earth System Science Interdisciplinary Center (ESSIC), University of Maryland, College Park, MD 20740, USA
*
Author to whom correspondence should be addressed.
Submission received: 10 January 2020 / Revised: 13 February 2020 / Accepted: 14 February 2020 / Published: 18 February 2020
(This article belongs to the Special Issue Computer Vision and Machine Learning Application on Earth Observation)

Abstract

:
One of the limitations in using spaceborne, microwave radiometer data for atmospheric remote sensing is the nonuniform spatial resolution. Remapping algorithms can be applied to the data to ameliorate this limitation. In this paper, two remapping algorithms, the Backus–Gilbert inversion (BGI) technique and the filter algorithm (AFA), widely used in the operational data preprocessing of the Advanced Technology Microwave Sounder (ATMS), are investigated. The algorithms are compared using simulations and actual ATMS data. Results show that both algorithms can effectively enhance or degrade the resolution of the data. The BGI has a higher remapping accuracy than the AFA. It outperforms the AFA by producing less bias around coastlines and hurricane centers where the signal changes sharply. It shows no obvious bias around the scan ends where the AFA has a noticeable positive bias in the resolution-enhanced image. However, the BGI achieves the resolution enhancement at the expense of increasing the noise by 0.5 K. The use of the antenna pattern instead of the point spread function in the algorithm causes the persistent bias found in the AFA-remapped image, leading not only to an inaccurate antenna temperature expression but also to the neglect of the geometric deformation of the along-scan field-of-views.

Graphical Abstract

1. Introduction

Microwave radiometers have wide applications in atmospheric remote sensing and provide essential inputs to numerical weather prediction models [1,2]. However, the application of these spaceborne, multispectral measurements from multiple sensors is often plagued with the problem of nonuniform spatial resolution caused by the limited size of the satellite instrument antenna and the frequency-dependent microwave emission from the earth–atmosphere system. The mismatch in resolution becomes a critical issue when observations from different frequencies and different sensors are combined to retrieve geophysical parameters like wind speed, cloud liquid water, total precipitation water, and snow depth [3,4,5,6]. It could also lead to uncertainty in the cross-calibration between similar satellite radiometer instruments [2,4,7,8].
To solve this problem, much effort has been made to manipulate the resolution of measurements from various sources, e.g., the Backus–Gilbert inversion (BGI) method [9,10,11,12,13,14,15,16], various filter algorithms [17,18,19], the truncated singular value decomposition method [20,21], and the scatterometer image-reconstruction algorithm [13,22,23]. Recently, many deep-learning-based super-resolution technologies have been proposed to improve the image quality of the natural environment [24,25,26,27,28,29]. Since the development by Hu et al. [30] of a deconvolution algorithm based on a convolutional neural network, the deep learning technique has been applied to microwave radiometer data to enhance their resolution [31,32,33]. However, among these algorithms, only a few manage to match the intrinsic antenna pattern of the observations. These algorithms are the focus of this study. These remapping techniques take advantage of the extra information provided by the overlaps between adjacent fields-of-view (FOVs) to produce observations with the expected antenna pattern. The result is an estimate of the antenna temperature (Ta) distribution the sensor would have measured, given a prescribed antenna pattern that is different from the actual one.
Currently available remapping algorithms for microwave sensors can be categorized as spatial-domain-based methods and frequency-domain-based methods. The BGI method represents the former. First developed by Backus and Gilbert [34,35] for solving integration equations, Stogryn [9] then applied the method to enhance the spatial resolution or reduce observations from spaceborne microwave instruments. Since then, the BGI method has been widely used as a remapping algorithm for multi-frequency microwave radiometers, such as the Special Sensor Microwave/Imager, the Microwave Radiation Imager, and the Advanced Technology Microwave Sounder (ATMS), among other instruments [11,12,15,16]. This method has been implemented in the ATMS resampling algorithm to produce ATMS brightness temperatures (Tbs) at each Cross-track Infrared Sounder FOV [36]. The filter algorithm represents the frequency-domain-based technique. Its theory is similar to frequency-domain-based deconvolution algorithms [17,18]. They both remove the impact of the original antenna pattern through deconvolution in the frequency domain. The filter algorithm for remapping does one more step of adding the influence of a target antenna pattern to the datasets in the same way. Hu et al. [33] used a similar method to generate high-resolution and low-resolution data pairs for training a deep learning residual network. Zou and Tian [6] used this algorithm to process ATMS data for hurricane warm-core retrievals. The filter algorithm was also applied in the ATOVS and AVHRR Pre-processing Package (AAPP) to remap ATMS data to Advanced Microwave Sounding Unit-A (AMSU-A)-like FOVs [37]. Hereafter, the filter algorithm will be called the AAPP Filter Algorithm (AFA). Although the BGI and the AFA have long been adopted in operational ATMS data pre-processing, and the remapped dataset is widely used by the scientific community, their strengths and weaknesses have never been compared.
This paper is devoted to the analysis and evaluation of these two remapping algorithms. For this purpose, the algorithms are applied to produce AMSU-A-like data from ATMS observations using both simulated and actual ATMS datasets. This paper is organized as follows. In Section 2, after a brief introduction to ATMS and AMSU-A instrument characteristics, descriptions of the BGI and AFA algorithms follow. Section 3 shows the point spread function (PSF) reconstructed by the BGI algorithm, which helps to understand some of the simulation results presented in this section. A comparison of the algorithms using actual ATMS data is also provided. Section 4 presents a discussion, and conclusions are made in Section 5.

2. Materials and Methods

2.1. ATMS and AMSU-A/MHS Instruments and Scan Geometries

The ATMS and AMSU-A/Microwave Humidity Sounder (MHS) are both cross-track scanning microwave sensors onboard the NOAA and MetOp series of polar-orbiting satellites, providing measurements of atmospheric thermal emission in the microwave region between 23 and 183 GHz. As a successor to and a combination of AMSU-A/MHS, the ATMS inherits most of its channels from its predecessor. ATMS channels 1–16 at frequencies ranging from 23 to 88 GHz are close to the AMSU-A channels, except that channel 4 at 51.76 GHz is a newly added ATMS channel, providing temperature information about the lower troposphere. ATMS samples have beam widths of 5.2° at 23 GHz and 31 GHz and 2.2° at the temperature-sounding channels. The beam widths of AMSU-A channels are all 3.3°. The ATMS channels located near the 183 GHz water vapor absorption line are similar to the MHS channels, except that channel 19 at 183.31 ± 4.5 GHz and channel 21 at 183.31 ± 1.8 GHz are new additions to the ATMS for better profiling atmospheric moisture. All channels of the two sensors operate with the same beam width of 1.1°. Therefore, by unifying the beam widths of AMSU-A and ATMS, their observations can be combined for studies of geophysical parameters and global climate changes.
The two sensors have different scan geometries. The ATMS has 96 scene resolution cells sampled at an interval of 8/3 s for each scan to cover 52.725° scan angles on both sides of the subsatellite path. By contrast, the AMSU-A scans the earth scene within ±49.5° of the nadir direction and has 30 FOVs. The larger scan angle of the ATMS creates a 2600 km swath width, which is larger than its predecessor AMSU-A (2343 km) and leaves almost no gaps near the equator, providing nearly full coverage of the low latitudes. The integration time for all ATMS channels is 18 ms, while that for AMSU-A channels 1–2 and 3–15 are 165 ms and 158 ms, respectively. Due to the shorter integration time, leading to a higher sampling rate, significant overlaps occur between the neighboring FOVs and between the neighboring scan lines of ATMS channels 1–16, but there is no overlap for AMSU-A FOVs. The oversampling feature of the ATMS allows for remapping the raw ATMS datasets to the observations with the AMSU-A-like antenna pattern. Because of the shorter integration time, the ATMS noise equivalent differential temperature (NEDT) is also larger than the corresponding AMSU-A channels. The comparison of the two sensors is well-documented in [38].
Note that due to the cross-scanning manner of the sensors and the curvature of the earth’s surface, the relative geometries of the data samples change over the scan. Figure 1a–b shows the ATMS channel-1 FOV (5.2° beam width) and the expected AMSU-A-like FOV (3.3° beam width) at nadir and at the scan edge, respectively. The beam width narrows when remapping from the ATMS channel-1 FOV to the AMSU-A-like FOV, referred to as resolution enhancement. Figure 1c–d shows the ATMS channel-3 FOV (2.2° beam width) and the AMSU-A-like FOV at nadir and scan ends, respectively. The beam width widens when remapping from the ATMS channel-3 FOV to the AMSU-A-like FOV, referred to as resolution degradation. Figure 1 illustrates that the relative geometries of the samples change significantly along the scan. Because the remapping algorithms are highly dependent on the overlaps between the raw antenna pattern and the expected one, the geometric deformation of the FOVs over the scan has a very important effect on the remapping algorithms [17,30,31,32].

2.2. ATMS Remapping Algorithms

This section describes the principles behind the BGI and AFA remapping algorithms. Since the beam widths of the ATMS channels 1 and 3 are broad compared with the sampling distance, the scan-motion smearing caused by the continuous scan mode of the ATMS is neglected to simplify the problem [37]. Based on this assumption, the Ta at location ρ 0 can be expressed as the convolution of the scene Tb with the antenna gain function G (normalized so that 4 π G   d Ω = 1 ):
T a ( ρ 0 ) = T b ( ρ ) · G ( ρ 0 , ρ ) d A
Note that the above integration requires that the variables Ta, Tb, and G are distributed on the same grids. Therefore, to be consistent with measurements, the antenna pattern needs to be projected onto the earth’s surface to generate the PSF.

2.2.1. Backus–Gilbert Inversion Algorithm

The BGI algorithm is a matrix inversion method for solving integral equations. In this algorithm, a set of optimal coefficients a i j is derived for generating Ta with the expected antenna pattern at the location ρ 0 ( T a t a r g e t ( ρ 0 ) ) as a linear combination of n × n numbers of adjacent original Tas at locations ρ i j ( T a s o u r c e ( ρ i j ) ):
T a t a r g e t ( ρ 0 ) = i = 1 n j = 1 n a i j T a s o u r c e ( ρ i j )
For definiteness, consider a n × n window stretching in the along-track and cross-track directions with T a t a r g e t at the center. Figure 1 shows the configuration. With the original gain function G s o u r c e , the variable T a s o u r c e ( ρ i j ) in Equation (2) can be expressed as Equation (1):
T a t a r g e t ( ρ 0 ) = ( i = 1 n j = 1 n a i j G s o u r c e ( ρ i j ,   ρ ) ) T b ( ρ ) d A
Comparing Equation (3) to Equation (1), the set of coefficients a i j makes the bracketed term in Equation (3) appear like the target gain function G t a r g e t ( ρ 0 ,   ρ ) . To obtain a tradeoff between the resolution enhancement and the noise suppression, optimal coefficients a i j are derived by minimizing the following objective function:
Q = Q 0 cos γ + e 2 w sin γ
w h e r e   Q 0 = ( i = 1 n j = 1 n a i j G s o u r c e ( ρ i j ,   ρ ) G t a r g e t ( ρ 0 ,   ρ ) ) 2 d A
e 2 = ( T r m s ) 2 i = 1 n j = 1 n a i j 2
where T r m s is the channel noise in terms of NEDT, w is a scale factor used to make the terms on the right side of Equation (4a) dimensionally and numerically compatible, and γ is a tradeoff factor that places degrees of emphasis on either the resolution enhancement or the noise suppression in the estimates of T a t a r g e t . Adjustable parameters are the window size n, w, and γ . These parameters have already been optimized in [16], where w is set to 0.001, n is set to 3 × 3 for resolution enhancement and 5 × 5 for resolution degradation, and γ is tuned for each FOV position to make the sum of the coefficients equal to 1. Because the resolution degradation does not amplify the noise as the resolution enhancement does, γ is set to zero in this case. These setups are used in this study. For reader completeness, the authors in [9,13] discuss some different criteria for selecting the parameters. Figure 1 shows the overlaps among the target FOV and the original FOVs within the selected calculation window.
The BGI algorithm in this study closely follows the one developed by Yang et al. [16], in which the antenna pattern is projected onto the earth’s surface, and the FOV deformation over the scan is taken into account by deriving a set of coefficients for each FOV position. One improvement made here is that the real antenna pattern, instead of the assumption of a Gaussian beam, is used in this algorithm. During a prelaunch ground test, antenna gain functions for 22 channels were measured at a resolution of 0.01° in an elevation angle and at intervals of 45° in azimuth angle at FOVs 1, 48, and 96 [39]. The antenna patterns at certain zenith and azimuthal angles at a certain FOV position are interpolated from the measurements through bilinear interpolation. Because the real AMSU-A antenna pattern is not available, a one-dimensional Gaussian function is adopted to approximate its gain function.

2.2.2. AAPP Filter Algorithm

The filter algorithm implemented in the AAPP manipulates the antenna pattern in the frequency domain [37]. Based on the convolution theorem, the expression of Ta in the spatial domain can be simplified to a multiplication in the frequency domain. Thus, the Fourier-space equivalent of Equation (1) can be written as
T a ^ ( ξ ) = T b ^ ( ξ ) · G ^ ( ξ )
where ξ represents frequency, and the superscript ^ indicates the Fourier transform.
In the AAPP, the ATMS antenna gain function is approximated by a one-dimensional Gaussian function and its Fourier transform can be easily obtained. Let G ^ s o u r c e ( ξ ) and G ^ t a r g e t ( ξ ) be the Fourier transforms of the original and expected gain functions, respectively. In the frequency domain, the AFA can easily manipulate the beam width by first dividing the original antenna pattern and then multiplying the target one:
T a ^ t a r g e t ( ξ ) = T a ^ s o u r c e ( ξ ) · G ^ t a r g e t ( ξ ) G ^ s o u r c e ( ξ )
The second term on the right side of the equation is called the manipulating term. For resolution enhancement, to suppress the noise level, a cutoff factor c is added to attenuate the high-frequency components in the signal. With the cutoff factor c, the manipulating term takes the following form:
G ^ t a r g e t ( ξ ) G ^ o r i g i n a l ( ξ ) · exp ( ( ln G ^ t a r g e t ( ξ ) ) 2 ln 2 ( ln c ) 2 )
In the AAPP, c is set to 0.4 for the remapping from beam width 5.2° to 3.3°, based on experimental results.
After the beam width is manipulated in the frequency domain, T a ^ t a r g e t ( ξ ) is transformed back to the spatial domain. The whole remapping process is now complete. This study uses the remapping codes implemented in the latest version of the AAPP, i.e., AAPP v8.5, released in December 2019. It includes the software atms_beamwidth.F, modify_beamwidth.F, and the parameter file atms_beamwidth.dat.
Note that unlike the BGI in which Equation (1) is strictly carried out on rectangular grids on the earth’s surface, the AFA works out this integration in the antenna-viewing coordinate system, assuming that the antenna pattern is a one-dimensional Gaussian function, and the observations are distributed evenly with the sampling distance 1.11°. Since Equation (1) uses the antenna pattern instead of its projection on Earth, not considered in this algorithm is the along-scan FOV deformation. In terms of computational efficiency, the BGI is superior to the AFA. The BGI coefficients need to be determined only once, given the radiometer configuration, and the inversion procedure reduces to a weighted sum of the measurements Equation (2). The Fourier transform and inverse Fourier transform need to be applied repeatedly in the AFA when processing multiple satellite images.

3. Results

3.1. Reconstructed PSF by the BGI

Following the procedure described in Section 2.2, BGI coefficients are generated for the ATMS onboard the NOAA-20 satellite. As shown in the objective function, i.e., Equation (4b), the reconstructed PSF can be expressed as the sum of the integration of the original PSF weighted by its corresponding coefficient over the overlap between the original and target FOVs. For the AFA, in which the antenna patterns are not projected onto the earth’s surface, reconstructed PSFs are not available.
To illustrate the BGI’s remapping capability, reconstructed PSFs are computed for resolution enhancement and degradation and compared with the original and target PSFs. Figure 2 shows the resolution enhancement from ATMS channel 1 to an AMSU-A-like channel. The top three panels are the original, synthetic, and target PSFs. To clarify the comparison, all PSFs are displayed within the 2.5 times beam width of the target FOV. The comparison shows that the reconstructed PSF is narrower than the original PSF and closer to the target PSF. Fast Fourier transforms convert these PSFs to the frequency domain (Figure 2d–f). The BGI enhances high-frequency components of the PSF to some extent. Figure 3 shows the resolution degradation from ATMS channel 3 to an AMSU-A-like channel. The BGI suppresses the high-frequency components in the PSF (Figure 3d–f) and makes the synthetic PSF match the target one well (Figure 3a–c).
Based on the original, synthetic, and target PSFs in the space domain, the beam widths of the response functions before and after correction can be estimated. To do so, the half-power points of the PSFs at nadir, i.e., the points with a value of 0.5 in Figure 2a–c and Figure 3a–c, are extracted. These points are then fitted with a circle formula (Figure 4). From the diameter of the circle and the altitude of the satellite (824 km for NOAA-20), the beam widths can be estimated (Table 1). Results show that the enhancement procedure narrows the beam width by 15%. The correction for resolution degradation is almost perfect, as expected.
The comparison between the PSFs before and after applying the BGI correction demonstrates that the BGI can achieve some improvement for the resolution enhancement and can near-perfectly match the PSFs for the resolution degradation. This can be explained by the noise penalty term in the objective function of Equation (4a). The resolution enhancement is essentially a procedure to amplify the high-frequency components in a signal. It is inevitable that the noise of the signal is amplified as well. The resolution has to be sacrificed to some extent to avoid introducing any drastic increases in the instrument noise. The compromise between resolution and noise is made through minimizing Equation (4a) with a properly selected tradeoff factor γ . That is why the BGI cannot perfectly match the PSFs and why the resolution enhancement it can achieve is limited. On the contrary, for resolution degradation, the high-frequency components in the signal are suppressed, as is the noise. The noise penalty term in Equation (4a) is no longer needed. Therefore, the BGI can near-perfectly match the synthetic PSF to the target one. Note that the AFA also has its own strategy to suppress the noise, i.e., adding a low-pass filter to the manipulating term in Equation (7). This may explain the phenomena observed in the simulations in the following section.

3.2. Simulation Comparison between the BGI and the AFA

While the BGI and AFA have a long-term application in ATMS data pre-processing, and the remapped datasets are widely used by the scientific community, their strengths and weaknesses have never been compared. Direct validation of the remapped fields using real ATMS data is nearly impossible because the responses of the different channels to a given scene are different. However, simulated data are an excellent substitute [12,16].
In this study, Ta observed with a certain antenna pattern are simulated by the Community Radiative Transfer Model (CRTM) for the atmosphere and surface geophysical parameters provided by the Global Forecast System (GFS) with a 0.12° grid resolution. The specific procedure is as follows. Through a geolocation model incorporating the ATMS instrument scan geometry, the PSF at each FOV is generated by projecting the real antenna pattern onto the earth’s surface. The PSF within the main beam width is divided into sub-beam intervals at a 3-km resolution. Then the scene Tb at each projected point is calculated from the GFS geophysical parameters using the CRTM. Equation (1) derives Ta from the convolution of Tb with the PSF. Finally, Gaussian white noise is added to the NEDT of each channel (0.22 K and 0.32 K for ATMS channels 1 and 3, respectively).

3.2.1. Evaluation of Remapping Capabilities

NOAA-20 ATMS observations with original and target antenna patterns were simulated for Hurricane Dorian at 18:00 UTC 31 August 2019 to investigate the remapping capabilities of these algorithms. ATMS data used in this study are obtained from the NOAA Comprehensively Large Array-data Stewardship System website: https://www.bou.class.noaa.gov/saa/products/welcome;jsessionid=429B2AC656CF5F31AF12161E3B819B39. Channel noise is added to the simulation of the original observations. The BGI and AFA are then applied to the original measurements to generate the remapped results. The remapped and target data are compared to make a quantitative evaluation of the remapping algorithms.
Figure 5 presents the simulation results for ATMS channel 1. Figure 5a–b are the ATMS channel-1 datasets simulated with 5.2° and 3.3° beam widths. These two datasets are regarded as the raw observations and the simulated truth, respectively. Figure 5c–d are the observations with the expected 3.3° beam width remapped by the BGI and the AFA. Figure 5e–f presents the bias between the remapped data and the simulated truth. Both algorithms provide better definitions of the hurricane center and the islands than the non-enhanced image. However, both produce large biases around coastlines. Specifically, negative biases appear over land and positive biases over water, indicating that the resolution enhancement for the sharp change in signal is insufficient. This is due to the fact that the beam width of the synthetic PSF is larger than that of the expected one, which is consistent with Figure 4a and Table 1. The remapped PSF therefore has higher gain from the ocean (lower Tbs) than expected when the beam center is over land and higher gain from land (higher Tbs) when the beam center is over water. As discussed in Section 3.1, the resolution enhancement has to be sacrificed to suppress the amplification of noise and that leads to the imperfect matching of the PSFs and consequently the insufficient enhancement for the sharp change in signal. However, the issue of the insufficient enhancement appears less serious in the BGI image than in the AFA image (see the west coastline of south Florida). But the noise in the BGI image is larger than that in the AFA image. This indicates that the BGI places more emphasis on resolution matching and less on noise suppression than the AFA.
The most distinct difference between the AFA and BGI remapping error images (Figure 5e–f) is a noticeable bias, as high as 2.5 K, found towards the ends of the scan lines in the AFA image. The AFA image also shows a larger bias around the hurricane center than shown in the BGI image. Focusing on the north coastline of Cuba, the bias from the AFA increases with scan angle, a pattern not obvious in the BGI bias image. This demonstrates that the remapping accuracy of the AFA deteriorates as the scan angle increases. The BGI does not have this issue. That the AFA transformation uses the antenna pattern instead of the PSF and that the geometric deformation of the along-scan FOVs is neglected may explain this. Section 4 discusses this further.
Figure 6 presents the simulation results for ATMS channel 3. The BGI image approximates the simulated truth well, consistent with the fact that the reconstructed PSF for resolution degradation perfectly matches the target one (Figure 3). The AFA image still has some bias around the coastlines, with positive biases over land and negative biases over water, showing that the smoothness over the sharp signal change is insufficient. The use of antenna gain instead of the PSF in the AFA may explain this. Both algorithms reduce the observation noise significantly. This is because the degradation process reduces the resolution and the noise of these channels substantially.
The bias and standard deviation of the error for both original and remapped observations at different FOV positions are calculated to investigate the performance of the algorithms along the scan lines (Figure 7). The most significant reduction in the bias and standard error occurs at FOVs that contain coastlines. For the hurricane center over the ocean observed around nadir, the bias and standard error reduce after remapping. From FOVs 60–70, where no sharp change in Ta appears, the resolution enhancement amplifies the standard error, and the resolution degradation alleviates the standard error, illustrating the noise changes caused by the remapping process. The bias at the scan end of channel 1 caused by the AFA is clearly seen, compared to the observations at other FOVs over the ocean.
PSF mismatch and the noise variation are both responsible for the remapping errors observed in Figure 5 and Figure 6. To quantitatively evaluate the noise change arising from the remapping algorithms and its contribution to the total error, the original channel NEDT, the noise component, and the root-mean-square error (RMSE) of the raw and remapped data relative to the simulated truth were calculated (Table 2). For the BGI, the noise component after remapping was derived from the original NEDT and the BGI coefficients through Equation (4c). For the AFA, without such an explicit expression, the standard error of FOVs 60–70 was regarded as the noise component. This assumption is reasonable because Tbs are distributed evenly within these FOVs, so the standard error is contributed by noise rather than a PSF mismatch.
The numbers in Table 2 show that both algorithms reduce the RMSE dramatically after remapping. Specifically, for resolution enhancement, both methods reduce the RMSE to a similar level, i.e., ~60% of the RMSE of the non-enhanced image. The BGI amplifies the noise component by 0.5 K, while the AFA maintains it at its original level. Given that the RMSE is the total contribution of the PSF mismatch and the noise component, this indicates that the BGI does a better PSF matching than the AFA. For resolution degradation, the BGI and the AFA reduce the RMSE to 11% and 16% of the original level, respectively. Both algorithms reduce the noise component by 0.2 K. This also indicates that the BGI matches PSFs better than the AFA.

3.2.2. The Impact of Antenna Pattern Irregularity

Apart from the way of incorporating the antenna gain function, the BGI and the AFA also differ in the choice of the antenna pattern. In the AFA, a Gaussian function approximates the ATMS antenna pattern while in the BGI, adopted is the real one measured during the prelaunch test. For a comparison with the PSFs projected from the real antenna pattern (Figure 2a and Figure 3a), Figure 8 shows the PSFs of the Gaussian-shaped antenna pattern. The irregularities of the real antenna patterns are clear. Note that real antenna patterns generate the Ta fields in the experiments, making the simulations more accurate. This section discusses the impact of the irregularity of the antenna pattern on remapping.
Chosen to carry out this sensitivity study is the BGI algorithm because of its flexibility in switching the antenna pattern. A new set of BGI coefficients generated with the Gaussian-shaped antenna pattern used in the AFA is applied to the raw observations. Figure 9 shows the difference between the remapped data and the simulated truth. Compared with Figure 5e and Figure 6e, approximating the ATMS antenna pattern by a Gaussian function leads to some bias around the hurricane center and coastlines. The RMSE of the enhanced and degraded images increases to 1.65 K and 0.15 K, respectively. As mentioned in Section 2, the remapping algorithms are highly dependent on the overlaps among the raw PSFs and the expected ones. Therefore, the accuracy of the PSF has a direct impact on the accuracy of remapping results. Using a Gaussian function to approximate the real antenna pattern leads to errors in remapping the coefficients and, consequently, to errors in the reconstructed images.
In summary, both the BGI and the AFA are capable of enhancing and degrading the resolution of data. The BGI matches the antenna pattern more accurately than the AFA, especially for resolution degradation. It produces less bias around coastlines, the hurricane center, and scan ends than the AFA. However, the BGI improves the resolution at the expense of increasing the noise level by 0.5 K, while the AFA manages to keep the noise at its original level. Both algorithms are able to reduce the noise level to a large extent in resolution degradation. Approximating the antenna pattern with a Gaussian function could lead to some bias in the remapped results.

3.3. Assessing the BGI and the AFA Using Actual Data

Having examined simulated images, we now consider actual data in a qualitative assessment of algorithm performance. Remapping the real measurements is similar to remapping the simulated dataset. Figure 10 shows the raw observations and remapping results of NOAA-20 ATMS channels 1 and 3 for Hurricane Dorian at 18:00 UTC 31 August 2019. As predicted by the simulation results, details of the hurricane center structure and the islands stand out more clearly in the enhanced images (Figure 10a,c,e), becoming obscure in the degraded images (Figure 10b,d,f). Although not explicitly shown in the figure, the AFA-enhanced image has a 1.5 K positive deviation from the BGI image at the edges of the scan lines. In the enhancement case, the BGI noticeably amplifies the instrument noise, well maintained by the AFA. In the degradation case, both algorithms manage to reduce the noise. These findings are consistent with the conclusions drawn from the simulations in the previous section.
The white lines in Figure 10a,b show a cross-section of Ta selected to show the impact of remapping around coastlines. Figure 11 shows a plot of the observations before and after the remapping process for channels 1 and 3. FOVs 5–10 are over the Gulf of Mexico, FOVs 20–25 are over the Atlantic Ocean, and FOV 15 is over Lake Okeechobee. For channel 1, after enhancement, more sharply delineated are the Atlantic and Gulf coasts, shown by an increase in the gradient across the coastlines. For channel 3, after degradation, decreased is the gradient across the coastlines, and smoothed out is Lake Okeechobee. For both cases, the BGI seems to have a slightly stronger effect in remapping than the AFA. This is expected because the simulations in Section 3.2 show that the BGI has a higher remapping accuracy than the AFA, and the issues of insufficiency in resolution enhancement and degradation around coastlines are less serious in the BGI than in the AFA.

4. Discussion

The application of the two remapping algorithms to the simulated and actual ATMS datasets shows that the remapping accuracy of the BGI is higher than that of the AFA. Compared to the BGI, the AFA-enhanced image has a larger bias around coastlines and the hurricane center where sharp changes in the signal occur. The remapping accuracy of the AFA deteriorates as the scan angle increases, producing more than a 1.5 K positive bias towards the scan ends. As for the degraded image, the bias around coastlines still exists in the AFA image but not in the BGI image. An explanation of the persistence of this bias found in AFA remapped images is that the AFA performs the remapping with the antenna pattern in the antenna-viewing coordinate system instead of projected onto the earth’s surface and neglects the geometric deformation of the along-scan FOVs. Equation (1) shows that remapping algorithms should be carried out in the image domain, the same domain in which the observations lie in. Therefore, what is needed is the transformation of the antenna pattern to a PSF before application in the remapping process. Figure 1 shows that the relative geometry of the along-scan FOV changes dramatically due to the ATMS scan mode and the curvature of the earth’s surface. Since remapping algorithms highly depend on the overlap among FOVs, space-variant PSFs should be used in the algorithms. The BGI incorporates these factors into its process by projecting the antenna pattern onto the earth’s surface and producing coefficients for each FOV position. This strategy makes its remapping more accurate than the AFA remapping.
Some filter-based deconvolution algorithms similar to the AFA have already incorporated the along-scan deformation of FOVs using space-variant PSFs [17,30,31,32]. In these algorithms, filtering is completed along each column of the image with a PSF at that specific scan position, since the relative geometry changes of the data in the along-track direction almost stay the same. The scan is usually divided into a few regions, and the filtering is carried out in each region with the same PSF to reduce the computation time. We believe that the remapping accuracy of the AFA can be improved by following the steps of these deconvolution algorithms. To avoid the bias caused by the irregularity of the antenna pattern, the real antenna pattern measured in the prelaunch ground test, instead of an approximation to a Gaussian function, should also be utilized in the remapping algorithms.

5. Conclusions

The applications of spaceborne, microwave radiometer data in atmospheric remote sensing are often hampered by the nonuniform spatial resolutions available from various sensors and at different frequencies. Much effort has been made to develop remapping algorithms unifying the resolutions to address this problem. In this paper, presented was a comparison of two typical methods, the Backus–Gilbert inversion (BGI) and a filter algorithm (AFA), remapping in spatial and frequency domains, respectively. They both have long been applied in Advanced Technology Microwave Sounder (ATMS) operational data preprocessing. An investigation of their strengths and weaknesses is thus necessary.
In this study, compared were the BGI and AFA algorithms via simulations and actual ATMS data. Both algorithms can produce Advanced Microwave Sounding Unit-A (AMSU-A)-like ATMS observations. The BGI matches the antenna pattern more accurately than the AFA. Specifically, for resolution enhancement, i.e., remapping ATMS channel-1 data with a 5.2° field-of-view (FOV) to an AMSU-A-like 3.3° FOV, both algorithms are capable of enhancing the resolution to some extent. The BGI produces less bias than the AFA around coastlines and hurricane centers where the antenna temperature changes sharply. In addition, the AFA has a noticeable bias towards the ends of scan lines. However, the BGI achieves the resolution enhancement at the expense of increasing the noise level by 0.5 K, while the AFA maintains the noise at its original level. For resolution degradation, i.e., remapping ATMS channel-3 observations with a 2.2° FOV to an AMSU-A-like 3.3° FOV, the BGI shows no obvious bias, while the AFA still has some bias around coastlines. Both algorithms suppress the noise level by 0.2 K.
The use of the antenna pattern instead of a point spread function in the AFA algorithm causes the persistent bias found in the AFA-remapped image. This not only leads to an inaccurate antenna temperature expression but also leads to the neglect of the geometric deformation of the along-scan FOVs. Neglecting the relative geometry change will make the calculation of the overlaps among FOVs inaccurate, which has an important effect on the remapping algorithms. Approximating the antenna pattern with a Gaussian function also contributes toward some bias in the remapped results. Therefore, to improve the remapping accuracy of the AFA, recommended is using the real antenna pattern and its projection onto the earth’s surface in the remapping procedure. The amplification of noise in resolution enhancement by the BGI, which outperforms the AFA by producing less bias around areas of large signal changes, is inevitable. A possible solution might be passing the enhanced image through a carefully designed low-pass filter. Further fine-tuning of the parameters in both algorithms to adapt them to specific Earth surface characteristics and FOV positions may also help improve their remapping qualities. Our future work will explore these topics to further improve the ATMS remapping accuracy.

Author Contributions

Methodology, investigation, formal analysis, validation, software, data curation, writing—original draft, J.Z.; conceptualization, supervision, funding acquisition, writing—review and editing, H.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by NOAA (grant nos. NA14NES4320003 and NA19NES4320002; Cooperative Institute for Satellite Earth System Studies-CISESS).

Acknowledgments

We gratefully acknowledge Fanglin Yang at NOAA’s Environmental Modeling Center for providing the high-resolution geophysical parameters of Hurricane Dorian simulated from the Global Forecast System. We thank M.C. Cribb at the University of Maryland for her assistance in editing this article. We are thankful for the insightful comments from three anonymous reviewers.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ulaby, F.T.; Moore, R.K.; Fung, A.K. Microwave Remote Sensing: Active and Passive, Volume I: Microwave Remote Sensing Fundamentals and Radiometry; Artech House: Norwood, MA, USA, 1981. [Google Scholar]
  2. Yang, H.; Weng, F.; Lv, L.; Lu, N.; Liu, G.; Bai, M.; Qian, Q.; He, J.; Xu, H. The FengYun-3 microwave radiation imager on-orbit verification. IEEE Trans. Geosci. Remote Sens. 2011, 49, 4552–4560. [Google Scholar] [CrossRef]
  3. Goodberlet, M.A.; Swift, C.T.; Wilkerson, J.C. Ocean surface wind speed measurements of the Special Sensor Microwave/Imager (SSM/I). IEEE Trans. Geosci. Remote Sens. 1990, 28, 823–828. [Google Scholar] [CrossRef]
  4. Yang, H.; Zou, X.; Li, X.; You, R. Environmental data records from FengYun-3B microwave radiation imager. IEEE Trans. Geosci. Remote Sens. 2012, 50, 4986–4993. [Google Scholar] [CrossRef]
  5. Li, X.; Zhao, K.; Wu, L.; Zheng, X.; Jiang, T. Spatiotemporal analysis of snow depth inversion based on the FengYun-3B microwave radiation imager: A case study in Heilongjiang Province, China. J. Appl. Remote Sens. 2014, 8, 084692. [Google Scholar] [CrossRef]
  6. Zou, X.; Tian, X. Hurricane Warm-Core Retrievals from AMSU-A and Remapped ATMS Measurements with Rain Contamination Eliminated. J. Geophys. Res. Atmos. 2018, 123, 10815–10829. [Google Scholar] [CrossRef]
  7. Yang, S.; Weng, F.; Yan, B.; Sun, N.; Goldberg, M. Special Sensor Microwave Imager (SSM/I) intersensor calibration using a simultaneous conical overpass technique. J. Appl. Meteorol. Climatol. 2011, 50, 77–95. [Google Scholar] [CrossRef]
  8. Wu, S.; Chen, J. Instrument performance and cross calibration of FY-3C MWRI. In Proceedings of the 2016 IEEE International Geoscience and Remote Sensing Symposium (IGARSS), Beijing, China, 10–15 July 2016; pp. 388–391. [Google Scholar]
  9. Stogryn, A. Estimates of brightness temperatures from scanning radiometer data. IEEE Trans. Antennas Propag. 1978, AP-26, 720–726. [Google Scholar] [CrossRef]
  10. Poe, G.A. Optimum Interpolation of Imaging Microwave Radiometer Data. IEEE Trans. Geosci. Remote Sens. 1990, 28, 800–810. [Google Scholar] [CrossRef] [Green Version]
  11. Farrar, M.R.; Smith, E.A. Spatial resolution enhancement of terrestrial features using deconvolved SSM/I microwave brightness temperatures. IEEE Trans. Geosci. Remote Sens. 1992, 30, 349–354. [Google Scholar] [CrossRef]
  12. Robinson, W.D.; Kummerow, C.; Olson, W.S. A technique for enhancing and matching the resolution of microwave measurements from the SSM/I instrument. IEEE Trans. Geosci. Remote Sens. 1992, 30, 419–427. [Google Scholar] [CrossRef]
  13. Long, D.G.; Daum, D.L. Spatial resolution enhancement of SSM/I data. IEEE Trans. Geosci. Remote Sens. 1998, 36, 407–417. [Google Scholar] [CrossRef] [Green Version]
  14. Wang, Y.; Shi, J.; Jiang, L.; Du, J.; Tian, B. The development of an algorithm to enhance and match the resolution of satellite measurements from AMSR-E. Sci. China Earth Sci. 2011, 54, 410–419. [Google Scholar] [CrossRef]
  15. Yang, H.; Shang, J.; Lu, L.; He, J.; Xu, H. Study of channel resolution matching of spaceborne microwave radiometer and its application in MWRI of FY-3 satellite. Aerosp. Shanghai 2012, 29, 23–51. [Google Scholar]
  16. Yang, H.; Zou, X. Optimal ATMS Remapping algorithm for Climate Research. IEEE Trans. Geosci. Remote Sens. 2014, 52, 7290–7296. [Google Scholar] [CrossRef]
  17. Sethmann, R.; Burns, B.A.; Heygster, G.C. Spatial resolution improvement of SSM/I data with image restoration techniques. IEEE Trans. Geosci. Remote Sens. 1994, 32, 1144–1151. [Google Scholar] [CrossRef]
  18. Piles, M.; Camps, A.; Vall-llossera, M.; Talone, M. Spatial-resolution enhancement of SMOS data: A deconvolution-based approach. IEEE Trans. Geosci. Remote Sens. 2009, 47, 2182–2192. [Google Scholar] [CrossRef]
  19. Liu, D.; Liu, K.; Lv, C.; Miao, J. Resolution enhancement of passive microwave images from geostationary Earth orbit via a projective sphere coordinate system. J. Appl. Remote Sens. 2014, 8, 083656. [Google Scholar] [CrossRef] [Green Version]
  20. Migliaccio, M.; Gambardella, A. Microwave Radiometer Spatial Resolution Enhancement. IEEE Trans. Geosci. Remote Sens. 2005, 43, 1159–1169. [Google Scholar] [CrossRef]
  21. Lenti, F.; Nunziata, F.; Estatico, C.; Migliaccio, M. On the spatial resolution enhancement of microwave radiometer data in Banach spaces. IEEE Trans. Geosci. Remote Sens. 2014, 52, 1834–1842. [Google Scholar] [CrossRef]
  22. Long, D.G.; Hardin, P.J.; Whiting, P.T. Resolution enhancement of spaceborne scatterometer data. IEEE Trans. Geosci. Remote Sens. 1993, 31, 700–715. [Google Scholar] [CrossRef] [Green Version]
  23. Early, D.S.; Long, D.G. Image reconstruction and enhanced resolution imaging from irregular samples. IEEE Trans. Geosci. Remote Sens. 2001, 39, 291–302. [Google Scholar] [CrossRef] [Green Version]
  24. Dong, C.; Chen, C.L.; He, K.; Tang, X. Learning a Deep Convolutional Network for Image Super-Resolution; Springer International Publishing: Berlin, Germany, 2014; pp. 184–199. [Google Scholar]
  25. Kim, J.; Lee, J.K.; Lee, K.M. Deeply-recursive convolutional network for image super-resolution. In Proceedings of the 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Las Vegas, NV, USA, 27–30 June 2016; pp. 1637–1645. [Google Scholar]
  26. Shi, W.; Caballero, J.; Huszár, F.; Totz, J.; Aitken, A.P.; Bishop, R.; Rueckert, D.; Wang, Z. Real-time single image and video super-resolution using an efficient sub-pixel convolutional neural network. In Proceedings of the 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Las Vegas, NV, USA, 27–30 June 2016; pp. 1874–1883. [Google Scholar]
  27. Schuler, C.J.; Burger, H.C.; Harmeling, S.; Schölkopf, B. A machine learning approach for non-blind image deconvolution. In Proceedings of the 2013 IEEE Conference on Computer Vision and Pattern Recognition, Portland, OR, USA, 23–28 June 2013; pp. 1067–1074. [Google Scholar]
  28. Burger, H.C.; Schuler, C.J.; Harmeling, S. Image denoising: Can plain neural networks compete with BM3D? In Proceedings of the 2012 IEEE Conference on Computer Vision and Pattern Recognition, Providence, RI, USA, 16–21 June 2012; pp. 2392–2399. [Google Scholar]
  29. Zhang, K.; Zuo, W.; Gu, S.; Zhang, L. 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; pp. 2808–2817. [Google Scholar]
  30. Hu, W.; Zhang, W.; Chen, S.; Lv, X.; An, D.; Lighart, L. A deconvolution technology of microwave radiometer data using convolutional neural networks. Remote Sens. 2018, 10, 275. [Google Scholar] [CrossRef] [Green Version]
  31. Hu, W.; Li, Y.; Zhang, W.; Chen, S.; Lv, X.; Lighart, L. Spatial resolution enhancement of satellite microwave radiometer data with deep residual convolutional neural network. Remote Sens. 2019, 11, 771. [Google Scholar] [CrossRef] [Green Version]
  32. Li, Y.; Hu, W.; Chen, S.; Zhang, W.; Guo, R.; He, J.; Ligthart, L. Spatial resolution matching of microwave radiometer data with convolutional neural network. Remote Sens. 2019, 11, 2432. [Google Scholar] [CrossRef] [Green Version]
  33. Hu, T.; Zheng, F.; Li, W.; Hu, W.; Tao, R. Microwave Radiometer Data Superresolution Using Image Degradation and Residual Network. IEEE Trans. Geosci. Remote Sens. 2019, 57, 8954–8967. [Google Scholar] [CrossRef]
  34. Backus, G.E.; Gilbert, J.F. Numerical applications of a formalism for geophysical inverse problems. Geophys. J. R. Astron. Soc. 1967, 13, 247–276. [Google Scholar] [CrossRef] [Green Version]
  35. Backus, G.E.; Gilbert, J.F. Resolving power of gross Earth data. Geophys. J. R. Astron. Soc. 1968, 16, 169–205. [Google Scholar] [CrossRef] [Green Version]
  36. Joint Polar Satellite System (JPSS) Advanced Technology Microwave Sounder (ATMS) SDR Calibration Algorithm Theoretical Basis Document (ATBD); Center for Satellite Applications and Research: College Park, MD, USA, 2013. Available online: https://www.star.nesdis.noaa.gov/jpss/documents/ATBD/D0001-M01-S01-001_JPSS_ATBD_ATMS-SDR_A.pdf (accessed on 13 February 2020).
  37. Atkinson, N.C. Annex to AAPP Scientific Documentation: Pre-Processing of ATMS and CrIS; Document (NWPSAF-MO-UD-027); EUMESAT and the Met Office: Devon, UK, 2011.
  38. Weng, F.; Zou, X.; Wang, X.; Yang, S.; Goldberg, M.D. Introduction to Suomi national polar-orbiting partnership advanced technology microwave sounder for numerical weather prediction and tropical cyclone applications. J. Geophys. Res. Atmos. 2012, 117, D19112. [Google Scholar] [CrossRef]
  39. Joint Polar Satellite System (JPSS) Advanced Technology Microwave Sounder (ATMS) Calibration Data Book; Document JPSS1 ATMS P/N 1362460-1 S/N 303; Northrop Grumman: Azusa, CA, USA, 2017.
Figure 1. Advanced Technology Microwave Sounder (ATMS) original (black) and expected (red) 3-dB footprints on the earth’s surface. The top panels show the resolution enhancement from the ATMS channel-1 5.2° FOV to the AMSU-A-like 3.3° FOV at (a) nadir and at (b) the scan edge. The bottom panels show the resolution degradation from the ATMS channel-3 2.2° FOV to the AMSU-A-like 3.3° FOV at (c) nadir and at (d) the scan edge.
Figure 1. Advanced Technology Microwave Sounder (ATMS) original (black) and expected (red) 3-dB footprints on the earth’s surface. The top panels show the resolution enhancement from the ATMS channel-1 5.2° FOV to the AMSU-A-like 3.3° FOV at (a) nadir and at (b) the scan edge. The bottom panels show the resolution degradation from the ATMS channel-3 2.2° FOV to the AMSU-A-like 3.3° FOV at (c) nadir and at (d) the scan edge.
Remotesensing 12 00672 g001
Figure 2. The Backus–Gilbert inversion (BGI) reconstruction from ATMS channel 1 to an AMSU-A point spread function (PSF) at nadir. (a) The original ATMS channel-1 PSF, (b) the reconstructed AMSU-A-like PSF, and (c) the target AMSU-A PSF. (df) are the Fourier transforms of the PSFs (ac). All response functions have been normalized to 1.0 to facilitate comparisons.
Figure 2. The Backus–Gilbert inversion (BGI) reconstruction from ATMS channel 1 to an AMSU-A point spread function (PSF) at nadir. (a) The original ATMS channel-1 PSF, (b) the reconstructed AMSU-A-like PSF, and (c) the target AMSU-A PSF. (df) are the Fourier transforms of the PSFs (ac). All response functions have been normalized to 1.0 to facilitate comparisons.
Remotesensing 12 00672 g002
Figure 3. The BGI reconstruction from ATMS channel 3 to an AMSU-A PSF at nadir. (a) the original ATMS channel-3 PSF, (b) the reconstructed AMSU-A-like PSF, and (c) the target AMSU-A PSF. (d–f) are the Fourier transforms of the PSFs (ac). All response functions have been normalized to 1.0 to facilitate comparisons.
Figure 3. The BGI reconstruction from ATMS channel 3 to an AMSU-A PSF at nadir. (a) the original ATMS channel-3 PSF, (b) the reconstructed AMSU-A-like PSF, and (c) the target AMSU-A PSF. (d–f) are the Fourier transforms of the PSFs (ac). All response functions have been normalized to 1.0 to facilitate comparisons.
Remotesensing 12 00672 g003
Figure 4. Half-power points and their circular curve fitting for (a) resolution enhancement and (b) degradation. The black, red, and blue markers and circular curves represent the half-power points and their fitting curves for the original, synthetic, and target PSFs, respectively.
Figure 4. Half-power points and their circular curve fitting for (a) resolution enhancement and (b) degradation. The black, red, and blue markers and circular curves represent the half-power points and their fitting curves for the original, synthetic, and target PSFs, respectively.
Remotesensing 12 00672 g004
Figure 5. Simulation results of ATMS channel-1 resolution enhancement from the original 5.2° FOV to the expected 3.3° FOV. Simulated (a) original and (b) expected observations; remapped observations by (c) the BGI and (d) the AAPP Filter Algorithm (AFA); the difference between the remapping results of (e) the BGI and (f) the AFA and the simulated truth (remapping minus truth). Units are K.
Figure 5. Simulation results of ATMS channel-1 resolution enhancement from the original 5.2° FOV to the expected 3.3° FOV. Simulated (a) original and (b) expected observations; remapped observations by (c) the BGI and (d) the AAPP Filter Algorithm (AFA); the difference between the remapping results of (e) the BGI and (f) the AFA and the simulated truth (remapping minus truth). Units are K.
Remotesensing 12 00672 g005
Figure 6. Simulation results of ATMS channel-3 resolution degradation from the original 2.2° FOV to the expected 3.3° FOV. Simulated (a) original and (b) expected observations; remapped observations by (c) the BGI and (d) the AFA; the difference between the remapping results of (e) the BGI and (f) the AFA and the simulated truth (remapping minus truth). Units are K.
Figure 6. Simulation results of ATMS channel-3 resolution degradation from the original 2.2° FOV to the expected 3.3° FOV. Simulated (a) original and (b) expected observations; remapped observations by (c) the BGI and (d) the AFA; the difference between the remapping results of (e) the BGI and (f) the AFA and the simulated truth (remapping minus truth). Units are K.
Remotesensing 12 00672 g006
Figure 7. Bias and standard error for ATMS observations before and after the remapping process for ATMS (a) channel 1 and (b) channel 3. Black curves and error bars represent the original observations. Red and blue curves and error bars represent the remapped observations by the BGI and the AFA, respectively.
Figure 7. Bias and standard error for ATMS observations before and after the remapping process for ATMS (a) channel 1 and (b) channel 3. Black curves and error bars represent the original observations. Red and blue curves and error bars represent the remapped observations by the BGI and the AFA, respectively.
Remotesensing 12 00672 g007
Figure 8. ATMS (a) channel-1 and (b) channel-3 PSFs projected from Gaussian-shaped antenna patterns.
Figure 8. ATMS (a) channel-1 and (b) channel-3 PSFs projected from Gaussian-shaped antenna patterns.
Remotesensing 12 00672 g008
Figure 9. Simulation results of resolution enhancement and degradation by the BGI with the coefficients produced by the Gaussian-shaped antenna pattern for ATMS channels (a) 1 and (b) 3. Units are K.
Figure 9. Simulation results of resolution enhancement and degradation by the BGI with the coefficients produced by the Gaussian-shaped antenna pattern for ATMS channels (a) 1 and (b) 3. Units are K.
Remotesensing 12 00672 g009
Figure 10. Application of the BGI and the AFA to real ATMS observations for resolution enhancement (left) and degradation (right). (a) Raw observations of ATMS channel 1, and (c) BGI and (e) AFA remapping results; (b) raw observations of ATMS channel 3, and (d) BGI and (f) AFA remapping results. The white lines in (a,b) indicate the location of the Ta cross-sections shown in Figure 11. Units are K.
Figure 10. Application of the BGI and the AFA to real ATMS observations for resolution enhancement (left) and degradation (right). (a) Raw observations of ATMS channel 1, and (c) BGI and (e) AFA remapping results; (b) raw observations of ATMS channel 3, and (d) BGI and (f) AFA remapping results. The white lines in (a,b) indicate the location of the Ta cross-sections shown in Figure 11. Units are K.
Remotesensing 12 00672 g010
Figure 11. NOAA-20 ATMS antenna temperature cross-section through the Florida peninsula for (a) channel-1 resolution enhancement and (b) channel-3 resolution degradation. The white lines in the top panels of Figure 10 show the location of the cross-section.
Figure 11. NOAA-20 ATMS antenna temperature cross-section through the Florida peninsula for (a) channel-1 resolution enhancement and (b) channel-3 resolution degradation. The white lines in the top panels of Figure 10 show the location of the cross-section.
Remotesensing 12 00672 g011
Table 1. The beam widths of the response functions before and after correction.
Table 1. The beam widths of the response functions before and after correction.
Beam Width (deg)
OriginalSyntheticTarget
Enhancement5.44.63.4
Degradation2.33.33.3
Table 2. Root-mean-square errors (RMSEs) and noise components for the raw and remapped observations relative to the simulated truth.
Table 2. Root-mean-square errors (RMSEs) and noise components for the raw and remapped observations relative to the simulated truth.
ChannelNEDTAlgorithmRMSE (K)Noise Component (K)
10.22none2.61
BGI1.590.69
AFA1.540.24
30.32none1.14
BGI0.130.09
AFA0.200.11

Share and Cite

MDPI and ACS Style

Zhou, J.; Yang, H. Comparison of the Remapping Algorithms for the Advanced Technology Microwave Sounder (ATMS). Remote Sens. 2020, 12, 672. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12040672

AMA Style

Zhou J, Yang H. Comparison of the Remapping Algorithms for the Advanced Technology Microwave Sounder (ATMS). Remote Sensing. 2020; 12(4):672. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12040672

Chicago/Turabian Style

Zhou, Jun, and Hu Yang. 2020. "Comparison of the Remapping Algorithms for the Advanced Technology Microwave Sounder (ATMS)" Remote Sensing 12, no. 4: 672. https://0-doi-org.brum.beds.ac.uk/10.3390/rs12040672

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