Feed-Forward Neural Network Denoising Applied to Goldstone Solar System Radar Images
Abstract
:1. Introduction
2. GSSR Radar Images
3. Selected Datasets
4. Denoising Methodology
- Process BPC image from a number of radar files. The total number of radar files available differs from track to track as it depends on the round-trip-time to the object.
- Apply Z-score normalization to the BPC image, following Equation (1). Store original and values.
- Apply denoising technique.
- Denormalize denoised BPC image using original and values.
- Apply Z-score normalization to the denoised BPC image, following Equation (1) using denoised and values.
- Daubechies: The Daubechies wavelet transform is a generalization of the Haar transform, which is implemented as a succession of decompositions. Daubechies allow for a filter length (N) bigger than 2, corresponding to the Haar transform, which provides more localization and smoothing effect, see Figure 1 for some examples. Daubechies transform are described by the number of vanishing moments, i.e., N/2 and are based on the use of compactly supported orthonormal wavelets, which makes discrete wavelet analysis feasible.
- Coiflet: The Coiflet wavelet transform is also based on scaling functions, similar to Daubechies wavelet, that exhibit vanishing moments.
- PCA: The PCA is a statistical technique that consists of simplifying a dataset by reducing it to a dataset composed of a lower number of dimensions. It uses orthogonal properties to transform a set of observations of possibly correlated variables into a set of values of uncorrelated variables. The denoising application intends to remove the noise while keeping the signal information. Basically, PCA is used to find the principal components that describe the maximum variance of the matrix (in this case an image), eliminating therefore those that contain residual information, expected to be related to the noise of the image. The image is then formed back from the reduced number of principal components which contain most of the information, and the noise is reduced.
- BM3D: Block-matching 3D (BM3D) filtering is a two-stage non-locally collaborative filtering method in the transform domain. This method is based on selecting 2D image fragments, known as blocks, and creating 3D arrays with those blocks that share similarities, known as groups, i.e., block matching. The 3D groups are then transformed into another domain, and a collaborative filtering is applied, that by attenuating the noise, reveals more fine details shared by grouped blocks and at the same time preserves the essential unique features of each individual block. This filtering can be a simple thresholding filter or something in the lines of a Wiener filter. After filtering the inverse transform is applied, and all the grouped images are aggregated to reconstruct the image. It is known that when the noise level is high, the BM3D performance decreases and can create undesired artifacts.
Feed-Forward Neural Network (FFNN)
5. FFNN Applied to Selected NEA
5.1. FFNN Set-Up: Data and Net
5.2. FFNN Training and Validation
5.3. FFNN Testing
6. Discussion: The Relevance of Denoised Radar Images
- Preservation of the radar scattering function after denoising is applied.
- Preservation of echo pixels close to the noise level.
- Preservation of the center of mass of the echo after denoising is applied.
- Denoised images to produce similar shape models for the asteroid, as compared to original images.
7. Conclusions
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Acknowledgments
Conflicts of Interest
References
- Asteroid Radar History Website. Available online: https://echo.jpl.nasa.gov/asteroids/PDS.asteroid.radar.history.html (accessed on 18 October 2021).
- Press Release on 3 September 2021: Planetary Radar Observes 1000th Near-Earth Asteroid Since 1968. Available online: https://www.jpl.nasa.gov/news/planetary-radar-observes-1000th-near-earth-asteroid-since-1968 (accessed on 18 October 2021).
- National Research Council. Defending Planet Earth: Near-Earth-Object Surveys and Hazard Mitigation Strategies; The National Academies Press: Washington, DC, USA, 2010. [Google Scholar]
- Lee, C.; Jao, J.; Slade, M.A.; Rodriguez-Alvarez, N. Micro-Meteoroid and Orbital Debris Radar from Goldstone Radar Observations. J. Space Saf. Eng. 2020, 7, 242–248. [Google Scholar] [CrossRef]
- McGuire, M.L.; Burke, L.M.; Hack, K.J.; Strange, N.J.; McElrath, T.P.; Landau, D.F.; Lantoine, G.; Lopez, P.; McDonald, M.A. Potential cislunar and interplanetary proving ground excursion trajectory concepts. In Proceedings of the 2016 IEEE Aerospace Conference, Big Sky, MT, USA, 5–12 March 2016; pp. 1–11. [Google Scholar] [CrossRef]
- Rodriguez-Alvarez, N.; Jao, J.S.; Lee, C.G.; Slade, M.A.; Lazio, J.; Oudrhiri, K.; Andrews, K.S.; Snedeker, L.G.; Liou, R.R.; Stanchfield, K.A. The Improved Capabilities of the Goldstone Solar System Radar Observatory. IEEE Trans. Geosci. Remote Sens. 2021, 60, 5103415. [Google Scholar] [CrossRef]
- Schroeder, D.J. Astronomical Optics, 2nd ed.; Academic Press: San Diego, CA, USA, 1999; p. 278. ISBN 978-0-12-629810-9. [Google Scholar]
- Daubechies, I. Ten Lectures on Wavelets; SIAM: Philadelphia, PA, USA, 1992; ISBN 978-0-89871-274-2. [Google Scholar]
- Zhang, L.; Dong, W.; Zhang, D.; Shi, G. Two-stage Image Denoising by Principal Component Analysis with Local Pixel Grouping. Pattern Recognit. 2010, 43, 1531–1549. [Google Scholar] [CrossRef] [Green Version]
- Dabov, K.; Foi, A.; Katkovnik, V.; Egiazarian, K. Image Denoising by Sparse 3-D Transform-Domain Collaborative Filtering. IEEE Trans. Image Process. 2007, 16, 2080–2095. [Google Scholar] [CrossRef] [PubMed]
- Thakur, R.S.; Yadav, R.N.; Gupta, L. State-of-art analysis of image denoising methods using convolutional neural networks. IET Image Process. 2019, 13, 2367–2380. [Google Scholar] [CrossRef]
- Fan, L.; Zhang, F.; Fan, H.; Zhang, C. Brief review of image denoising techniques. Vis. Comput. Ind. Biomed. Art 2019, 2, 7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Vojtekova, A.; Lieu, M.; Valtchanov, I.; Altieri, B.; Old, L.; Chen, Q.; Hroch, F. Learning to denoise astronomical images with U-nets. Mon. Not. R. Astron. Soc. 2021, 503, 3204–3215. [Google Scholar] [CrossRef]
- Slade, M.A.; Benner, L.A.M.; Silva, A. Goldstone Solar System Radar Observatory: Earth-Based Planetary Mission Support and Unique Science Results. Proc. IEEE 2010, 99, 757–769. [Google Scholar] [CrossRef]
- Levenberg, K. A Method for the Solution of Certain Problems in Least Squares. Q. Appl. Math. 1944, 2, 164–168. [Google Scholar] [CrossRef] [Green Version]
- Marquardt, D. An Algorithm for Least-Squares Estimation of Nonlinear Parameters. J. Soc. Ind. Appl. Math. 1963, 11, 431–441. [Google Scholar] [CrossRef]
- Hudson, S. Three-dimensional reconstruction of asteroids from radar observations. Remote Sens. Rev. 1993, 8, 195–203. [Google Scholar] [CrossRef]
- Magri, C.; Ostro, S.J.; Scheeres, D.; Nolan, M.; Giorgini, J.D.; Benner, L.A.; Margot, J.-L. Radar observations and a physical model of Asteroid 1580 Betulia. Icarus 2007, 186, 152–177. [Google Scholar] [CrossRef]
- Naidu, S.P.; Margot, J.-L.; Busch, M.W.; Taylor, P.A.; Nolan, M.C.; Brozovic, M.; Benner, L.A.M.; Giorgini, J.D.; Magri, C. Radar imaging and physical characterization of near-Earth Asteroid (162421) 2000 ET70. Icarus 2013, 226, 323–335. [Google Scholar] [CrossRef] [Green Version]
Name | H | P [h] | D [m] | |
---|---|---|---|---|
2018 DH1 | 21.13 | 5 | unknown | 200 * |
(505657) 2014 SR339 | 18.54 | 8.71 | 0.068 | 971 ^ |
(441987) 2010 NY65 | 21.39 | 4.97 | 0.071 | 228 ^ |
2017 YE5 1 | 19.2 | 20.6 | dark | 860 * |
2018 EB | 21.9 | 3.16 | 0.07 * | 240 ^ |
2017 VR12 | 20.5 | 1.38 | bright | 250 * |
2017 WX12 | 22.1 | 16.23 | unknown | 100 * |
2016 AJ193 | 18.99 | unknown | 0.032 | 1374 + |
(444584) 2006 UK | 20.2 | 5.72 | unknown | 300 * |
Method | Image Alteration | |
---|---|---|
Daubechies 1 WT | ×4.53 | Yes-smoothing |
Daubechies 4 WT | ×3.19 | Yes-smoothing |
Coiflet 1 WT | ×3.32 | Yes-smoothing |
Coiflet 4 WT | ×1.58 | Yes-smoothing |
PCA | increased | Yes-shape |
BM3D | ×1.25 | Yes-shape |
FFNN | ×3.66 | Not evident |
Object | Processing Information |
---|---|
2018 DH1 | FFT length = 2048, frequency = 0.5 Hz, code length = 511, baud rate = 0.125 us |
2014 SR339 | FFT length = 2048, frequency = 1 Hz, code length = 255, baud rate = 0.5 us |
2010 NY65 | FFT length = 1024, frequency = 0.1 Hz, code length = 127, baud rate = 0.125 us |
2017 YE5 | FFT length = 2048, frequency = 0.025 Hz, code length = 255, baud rate = 0.25 us |
1988 XB | FFT length = 512, frequency = 0.5 Hz, code length = 127, baud rate = 0.5 us |
2016 AJ193 | FFT length = 254, frequency = 3 Hz, code length = 255, baud rate = 0.25 us |
2008 EB | FFT length = 256, frequency = 0.5 Hz, code length = 127, baud rate = 0.25 us |
2017 VR12 | FFT length = 1024, frequency = 1 Hz, code length = 511, baud rate = 0.125 us |
2017 WX12 | FFT length = 256, frequency = 0.04 Hz, code length = 127, baud rate = 0.25 us |
2006 UK | FFT length = 511, frequency = 0.25 Hz, code length = 255, baud rate = 0.125 us |
2000 ET70 | FFT length = 512, frequency = 1 Hz, code length = 255, baud rate = 0.25 us |
Asteroid | Noise Reduction | ||
---|---|---|---|
2016 AJ193 | 0.04 | 0.01 | ×4 |
2018 EB | 0.7 | 0.22 | ×3.2 |
2017 VR12 | 11.14 | 3.93 | ×2.83 |
2017 WX12 | 0.04 | 0.01 | ×4 |
2006 UK | 0.63 | 0.22 | ×2.86 |
Object | RMSE | Correlation Coefficient |
---|---|---|
2014 SR339 (Figure 12a) | 3.690% | 90.5% |
2006 UK (Figure 12b) | 1.670% | 72.9% |
2017 VR12 (Figure 12c) | 0.062% | 99.7% |
2008 EB (Figure 12d) | 0.610% | 84.9% |
2017 WX12 (Figure 12e) | 1.260% | 88.5% |
2016 AJ193 (Figure 12f) | 1.310% | 98.1% |
2010 NY65 (Figure 12g) | 0.061% | 98.2% |
2018 DH1 (Figure 12h) | 1.210% | 99.1% |
Object | Original image [Doppler, Delay] [Hz, Bin] | FFNN Denoised [Doppler, Delay] [Hz, Bin] | Difference [Doppler, Delay] [Hz, Bin] |
---|---|---|---|
2018 DH1 | [101, 477] | [101, 477] | [0, 0] |
2014 SR339 | [500, 267] | [500, 266] | [0, 1] |
2010 NY65 | [20.9, 250] | [20.9, 249] | [0, 1] |
2017 YE5 | [10.1, 253] | [10.1, 253] | [0, 0] |
1988 XB | [20, 256] | [20, 255] | [0, 1] |
2016 AJ193 | [98.05, 252] | [98.05, 251] | [0, 1] |
2008 EB | [10.25, 271] | [10.75, 271] | [−0.5, 0] |
2017 VR12 | [20.05, 328] | [20.05, 328] | [0, 0] |
2017 WX12 | [−0.28, 242] | [−0.24, 242] | [−0.04, 0] |
2006 UK | [−35.25, 239] | [−35.25, 239] | [0, 0] |
2000 ET70 | [−80.68, 290] | [−80.68, 287] | [0, 3] |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Rodriguez-Alvarez, N.; Jao, J.S.; Munoz-Martin, J.F.; Lee, C.G.; Oudrhiri, K. Feed-Forward Neural Network Denoising Applied to Goldstone Solar System Radar Images. Remote Sens. 2022, 14, 1643. https://doi.org/10.3390/rs14071643
Rodriguez-Alvarez N, Jao JS, Munoz-Martin JF, Lee CG, Oudrhiri K. Feed-Forward Neural Network Denoising Applied to Goldstone Solar System Radar Images. Remote Sensing. 2022; 14(7):1643. https://doi.org/10.3390/rs14071643
Chicago/Turabian StyleRodriguez-Alvarez, Nereida, Joseph S. Jao, Joan Francesc Munoz-Martin, Clement G. Lee, and Kamal Oudrhiri. 2022. "Feed-Forward Neural Network Denoising Applied to Goldstone Solar System Radar Images" Remote Sensing 14, no. 7: 1643. https://doi.org/10.3390/rs14071643
APA StyleRodriguez-Alvarez, N., Jao, J. S., Munoz-Martin, J. F., Lee, C. G., & Oudrhiri, K. (2022). Feed-Forward Neural Network Denoising Applied to Goldstone Solar System Radar Images. Remote Sensing, 14(7), 1643. https://doi.org/10.3390/rs14071643