Inverted EEMD: a robust method to identify narrow absorption features form spectral data and cubes

Zhenzhen He,Guang-Xing Li
2024-10-19
Abstract:Extracting information from complex data is a challenge shared by multiple frontiers of modern astrophysical research. Among those, analyzing spectra cubes, where the emission is mapped in the position-position-velocity space is a difficult task given the vast amount of information contained within. The cubes often contain a superposition of emissions and absorptions, where extracting absorption signatures is often necessary. One example is the extraction of narrow absorption structures in HI 21 cm emission spectra. These HI self-absorption (HISA) clouds trace the cold HI gas in interstellar space. We introduce an automatic and robust method called the \emph{inverted EEMD} method to extract narrow features from spectral cubes. Our method is based on the EEMD method, an established method to decompose 1d signals. The method is robust and parameter-free, making it useful in analyzing spectral cubes containing localized absorption signals of different types. The inverted-EEMD method is suitable for the analysis of spectral cubes where it can produce a cube containing the absorption signal and one containing the unabsorbed signal, where cold clouds can be identified as coherent regions in the absorption map. A Python implementation of the method is available at \url{<a class="link-external link-https" href="https://github.com/zhenzhen-research/inverted_eemd_map" rel="external noopener nofollow">this https URL</a>}.
Instrumentation and Methods for Astrophysics,Astrophysics of Galaxies
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to extract narrow absorption features from complex spectral data, especially to identify narrow self - absorption features (HISA) from the HI 21 - centimeter emission spectral line. Specifically, the paper proposes a new method - Inverted Ensemble Empirical Mode Decomposition (Inverted EEMD) to automatically and robustly separate these narrow absorption features from the spectral cube. ### Main problem background 1. **Information extraction from complex data**: One of the main challenges in modern astrophysics research is to extract useful information from complex data. In particular, when analyzing spectral cubes, the emission and absorption signals are usually superimposed, making it very difficult to separate these signals. 2. **Importance of narrow absorption features**: Narrow absorption features, such as HI self - absorption (HISA) clouds, are crucial for tracing cold HI gas in interstellar space. These cold gas clouds play an important role in the galactic ecosystem, and understanding them helps to reveal the relationships between different gas phases and their circulation mechanisms. ### Limitations of existing methods - **Gaussian component assumption**: Some existing methods assume that the spectrum is composed of multiple positive and negative Gaussian components superimposed, and identify cold gas by fitting these components. However, this method has limitations in dynamic range. - **Second - derivative method**: Other methods attempt to use the second derivative on the frequency axis to extract absorption features, but this requires the data to have a high signal - to - noise ratio. - **Algorithm implementation not public**: The implementation codes of many existing methods are not public, which makes it difficult to evaluate their effectiveness and also hinders their wide application. ### Solution proposed in the paper The paper proposes a new method based on Ensemble Empirical Mode Decomposition (EEMD) - Inverted EEMD, for extracting narrow absorption features from the spectral cube. The main features of this method include: - **Robustness and parameter - free**: The Inverted EEMD method does not need to set additional parameters and can handle nonlinear features, and is especially suitable for complex and rapidly changing HI emission profiles. - **Automated processing**: This method can be automatically applied to each spectrum to generate two cubes containing the absorption signal and the unabsorbed signal. - **Strong applicability**: This method is suitable for analyzing different types of spectral cubes containing local absorption signals, and is especially suitable for identifying coherent regions of cold gas clouds. ### Method overview 1. **Spectral extraction**: Extract the spectrum \(I(v)\) at each position from the observed HI 21 - centimeter emission data. 2. **Spectral decomposition**: Apply the EEMD algorithm to decompose the spectrum \(I(v)\) into multiple Intrinsic Mode Functions (IMFs), and extract the absorption feature \(A(v)\) and the background signal \(B(v)\) from them. 3. **Cube reconstruction**: Reconstruct the absorption feature and the background signal at each position into two new cubes, which contain all the absorption signals and the unabsorbed signals respectively. Through the above steps, this method can effectively separate narrow absorption features from complex spectral cubes, providing strong support for further research on cold gas clouds.