Low Complexity Radio Frequency Interference Mitigation for Radio Astronomy Using Large Antenna Array

Zaid Bin Tariq,Teviet Creighton,Louis P. Dartez,Naofal Al-Dhahir,Murat Torlak
2024-03-07
Abstract:With the ongoing growth in radio communications, there is an increased contamination of radio astronomical source data, which hinders the study of celestial radio sources. In many cases, fast mitigation of strong radio frequency interference (RFI) is valuable for studying short lived radio transients so that the astronomers can perform detailed observations of celestial radio sources. The standard method to manually excise contaminated blocks in time and frequency makes the removed data useless for radio astronomy analyses. This motivates the need for better radio frequency interference (RFI) mitigation techniques for array of size M antennas. Although many solutions for mitigating strong RFI improves the quality of the final celestial source signal, many standard approaches require all the eigenvalues of the spatial covariance matrix ($\textbf{R} \in \mathbb{C}^{M \times M}$) of the received signal, which has $O(M^3)$ computation complexity for removing RFI of size $d$ where $\textit{d} \ll M$. In this work, we investigate two approaches for RFI mitigation, 1) the computationally efficient Lanczos method based on the Quadratic Mean to Arithmetic Mean (QMAM) approach using information from previously-collected data under similar radio-sky-conditions, and 2) an approach using a celestial source as a reference for RFI mitigation. QMAM uses the Lanczos method for finding the Rayleigh-Ritz values of the covariance matrix $\textbf{R}$, thus, reducing the computational complexity of the overall approach to $O(\textit{d}M^2)$. Our numerical results, using data from the radio observatory Long Wavelength Array (LWA-1), demonstrate the effectiveness of both proposed approaches to remove strong RFI, with the QMAM-based approach still being computationally efficient.
Signal Processing,Instrumentation and Methods for Astrophysics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to efficiently mitigate radio - frequency interference (RFI) in radio astronomical observations. With the development of radio communication technology, radio - frequency interference has caused increasing pollution to radio astronomical data, which has hindered astronomers' research on cosmic radio sources. Especially when studying short - lived radio transient sources, quickly mitigating strong RFI becomes particularly important because these transient sources require astronomers to take measures quickly for detailed observations. The existing standard methods are usually to manually cut out the contaminated time and frequency segments, but this will make the cut - out data unavailable for radio astronomical analysis. In addition, many existing RFI mitigation techniques can improve the quality of the final celestial signal, but usually need to calculate all eigenvalues of the spatial covariance matrix of the received signal. The computational complexity of this process is \(O(M^3)\), where \(M\) is the size of the antenna array, and \(d\) is the size of the RFI to be removed, and \(d \ll M\). This high computational complexity limits the application of these methods in real - time or near - real - time processing, especially in time - domain radio astronomy, where the need for fast processing is more urgent. Therefore, this paper proposes two new low - complexity RFI mitigation methods: 1. **Fast QMAM (Quadratic Mean to Arithmetic Mean) method based on the Lanczos method**: Using previously collected data under similar radio sky conditions, find the Rayleigh - Ritz values of the covariance matrix through the Lanczos method, thereby reducing the computational complexity of the overall method to \(O(dM^2)\). 2. **RFI mitigation method with celestial sources as references**: This method uses a known celestial source as a reference to mitigate RFI. These two methods aim to provide a faster and more efficient RFI mitigation tool, especially for transient source monitoring in time - domain radio astronomy. Through numerical experiments using data from the Long Wavelength Array (LWA - 1), the effectiveness of these two methods has been proven, especially in significantly improving the quality of the target signal while maintaining computational efficiency.