Abstract

Low frequency oscillation is an important attribute of human brain activity, and the amplitude of low frequency fluctuation (ALFF) is an effective method to reflect the characteristics of low frequency oscillation, which has been widely used in the treatment of brain diseases and other fields. However, due to the low accuracy of the current analysis methods for low frequency signal extraction of ALFF, we propose the Fourier-based synchrosqueezing transform (FSST), which is often used in the field of signal processing to extract the ALFF of the low frequency power spectrum of the whole-time dimension. The low frequency characteristics of the extracted signal are compared with those of FSST and fast Fourier transform (FFT) through the resting-state data. It is clear that the signal extracted by FSST has more low frequency characteristics, which is significantly different from FFT.

1. Introduction

In recent years, with the rapid development of big data, especially artificial intelligence technology, intelligent methods have become increasingly popular in the operation and management optimization of complex large-scale systems, which often appear in ecological environments, communication networks, informatics, biology, and other disciplines. Extensive research has focused on medicine and health such as deep neural network, neuroscience, and brain science related purposes, data collection, accurate analysis, monitoring, and connecting available medical resources and healthcare services [16].

The brain is comprised of 1011 neurons and generally 1015 interconnections among them can handle massive information. Their shape and size are diverse, and their structures are distinct. The components of a neuron can be categorized as cell body, dendrite, and axon. Signals received from other neurons are reckoned into synapses and sent to the cell body. Once the signal entering the cell body surpasses the sustaining threshold, the neuron is burned, and the signal is transmitted to other neurons through axons [7].

The amplitude of low frequency fluctuation (ALFF) is a good measurement index of brain activation signal, which is calculated to obtain the power spectrum of brain voxel signal time series using time-frequency (TF) analysis method [811]. It is considered able to directly reflect the spontaneous synchronous changes of neural activity in a resting state [12, 13] and to some extent reflect the interaction and neural network connection among the relevant brain regions, which has been verified in the research of visual stimulus differences caused by eye opening and closing. Furthermore, ALFF can serve as a starting point for understanding brain diseases and can help reveal many pathological mechanisms, such as some mental diseases [14, 15]. Therefore, improving the analytical accuracy of ALFF has a positive effect on the improvement of brain science research and the treatment of mental diseases [1619].

Many signals such as audio signals and medical data (electrocardiogram and thoracic motion signals) can be modeled as superposition of amplitude-modulation and frequency-modulation (AM-FM) modes. Linear techniques such as wavelet transform and fast Fourier transform (FFT) are widely utilized to represent the characteristics of the signal in TF domain [2024]. However, they all have the same limitation; that is, they cannot locate a signal with arbitrary accuracy in time and frequency. Even though an improved TF representation method called “reassignment” method (RM) was proposed, its main problem is that the reassigned transformation is no longer reversible and does not allow for the mode reconstruction. Moreover, it is difficult to achieve compatibility between the time and frequency resolutions because the frequency resolution depends on the length of the analysis window, and the analytical window length needs to be fixed according to the frequency. In particular, if the period of a signal is not constant, such as the sample signal for the object being relatively active, the analytical precision decreases [25, 26]. Therefore, FFT cannot analyze the characteristics of the signal at a certain time accurately, resulting in a large amount of signal information losses, which is unsuitable for brain science analysis particularly.

In the analysis of audio signal [27, 28], another phase-based technique called “synchrosqueezing transform” (SST) was proposed. Its purpose is similar to that of RM, which is to sharpen the time-scale representation given by wavelet transform, and has the advantage of allowing for mode retrieval. Based on the principle of wavelet-based SST, an extension of SST to the TF representation given by short-time Fourier transform called Fourier-based synchrosqueezing transform (FSST) was proposed [29], and it is proved that this method is robust to small bounded perturbations and noise [30]. Moreover, some researches proposed improving the existing FSST by using higher-order amplitude and phase approximation to calculate more accurate instantaneous frequency estimates of the modal components of the signal. It is achieving a perfect concentration and reconstruction in a wider range of AM-FM modes, and most of the real signals are composed of strongly modulated AM-FM modes, such as chirps involved in radar [31], speech processing [32], or gravitational waves [9, 11, 33]. As a consequence, this technique provides a highly concentrated TF representation for a wide variety of AM-FM multicomponent signals and enables the reconstruction of their modes with high accuracy; we hypothesized that FSST can be used in the analysis of time series of brain science and to improve the analytical accuracy.

In the present study, we propose and attempt the application of FSST to brain wave analysis and compare the results obtained by FFT and FSST analysis using data with no task resting state. Based on learning and improving the previous ALFF algorithm, the ALFF in the overall time domain is recalculated. According to the comparison with traditional ALFF values, we demonstrate that FSST is a highly precise method, which yields marked improvements in the reconstruction of TF signal. This technique is effective for analyzing and visualizing ALFF values and changes in brain waves.

2. Materials and Methods

The FFT as a Fourier-related transform has been commonly used for analyzing organic signals, such as the brain wave signal. The Fourier transform can be represented as follows:where t is the analysis window length (f1 = 1/t, n is a positive integer). Equation (1) is solved for determining the Fourier coefficients. As the Fourier transform is used for analyzing a completely periodic signal in an analysis window t, the calculated frequencies (nf1) depend on the window length t, and errors frequently occur in the analysis of nonharmonic signal frequencies. Moreover, if changes in the frequency of the sample signal for the object are relatively active over a short duration, a decrease in the length of the analysis window to increase the time resolution also decreases the frequency resolution; the resulting analytical precision decreases.

As we all know, it is not appropriate to describe the nonstationary signal with time-domain or frequency-domain representation. Therefore, the short-time Fourier transform (STFT) is introduced and defined as follows:where is the complex conjugate of . Then, the spectrum corresponds to , and the original signal f can be extracted from the STFT under the condition that does not disappear and is continuously at 0:

If f is analytic, which means , then ; the integral in (3) only happens on R+.

Synchrosqueezing transform can separate and adjust the signals of different modes by redistributing coefficients in series or frequency to avoid the uncertainty of linear transformation. It can analyze and compress the low frequency amplitude, obtain more accurate high-power signals, and find active points in brain regions. The signal can be expressed clearly while maintaining the reversibility of the signal.

The Fourier-based synchrosqueezing transform (FSST) technique is a newly proposed analytical method that has some advantages, as evidenced by the numerical processing of synthetic and observed gravitational wave signals [34, 35]. Time t and frequency are predicted by sharpening the blurred STFT representation using the following instantaneous frequency of FSST:where ∂t is the partial derivative with respect to t. Indeed, is redistributed to a new location () using the following synchrosqueezing operator:where γ is the threshold. By integrating near the corresponding ridge , the kth mode is reconstructed approximately, since the coefficients are redistributed along the frequency axis, and FSST maintains the nature of causality:where is an estimate of and the parameter d enables compensating for the imprecise approximation of and the error made by estimating the IF by means of . In addition, the approximation must be calculated before retrieving mode . As a consequence, a commonly used technique that is based on ridge extraction predicting and K is known [33]. Therefore, FSST as a highly precise frequency analysis method can be used to analyze brain waves and improve the analytical precision in the present study.

2.1. Experiment

Forty-seven participants in their 20s–40s with no mental illness history participated in the present study of resting state. All participants provided written consent to participate, and ethical approval was obtained from the institutional ethics committee [36, 37]. For avoiding the influence of outside equipment, the experiment was conducted in an anechoic chamber with lights off.

In functional imaging, echo planar imaging and fast gradient echo image were obtained by three-dimensional T1-weighted magnetization with the relevant information shown in Table 1, respectively [3840].

All of forty-seven participants’ closed-eye sample data were used in the resting experiment. Each condition consists of 240 functional volumes; in terms of data selection, we discard the original 10 point volumes and reserve 120 volumes.

The improvement of the algorithm is mainly reflected in the calculation and selection of power spectrum; in terms of analytical conditions, as the brain activity is the most obvious and active in this piece of data, FFT and FSST both select the average value in the time series with the frequency band of 0.01–0.1 Hz and the data length of 20–420 s in the total length of 9 minutes. The analysis window length is 2ABS (data)/sample length. Furthermore, FFT and FSST algorithms are used to extract ALFF values and conduct comparative analysis; the coefficient values of the corresponding ALFF are extracted as follows (Table 2).

3. Results and Discussion

The experimental data of each participant were extracted twice. We divided these two groups of data into Test1 and Test2; FFT and FSST are used to analyze the two groups of data. Results of participant 32 as an example are shown in Figure 1 since it is nonspecific, which are Test1 FSST, Test1 FFT, Test2 FSST, and Test2 FFT, respectively. Each group of images shows 12 cross sections of the brain, and the corresponding position of each cross section is shown in the longitudinal view of brain area in the lower right corner of each group of images.

Although the high amplitude of FSST and FFT is only about 3000 in the resting state, the difference between FSST and FFT is obvious. From the point of view of images, FSST detects the results that cannot be resolved by FFT and their mean coefficients, which proves that FSST is more accurate and sensitive than FFT in data analysis and image extraction and highlighting. In particular, FSST highlights the superior frontal fornix, rectus, and other areas, which are valuable for the scientific research of brain activity and mental disease. Therefore, the difference between FSST and FFT is large, whose accuracy is also higher.

Since the amplitude of the low frequency band is used in this experiment, the average sum of these values should be carried out in FFT analysis of these multiple detection points, and then the frequency power spectrum is obtained by judging the ratio of the low frequency band (0.01–0.1 Hz) to each frequency band; the obtained frequency power spectrum is analyzed to determine whether these points conform to the low frequency amplitude attribute. To further analyze the significance of these signals, we introduce the mean correlation coefficient to correlate the extracted signals and calculate the average value, as shown in the figure below.

We can clearly see that, in both the mean coefficient and the low frequency ratio of the two groups of data, the multiple detection parts of FSST are indeed better than those of FFT, as shown in Figure 2. Since the significance of the low frequency coefficient can reflect the significance of the ALFF index and the difference of the low frequency amplitude is the reason for the difference of the low frequency coefficient, the low frequency amplitude value targeted by this experiment can show the characteristics of ALFF. Therefore, it is also proved that FSST is clearly superior to FFT in the analysis of ALFF value, which has stronger low frequency characteristics and more obvious ability to extract low frequency characteristics.

Figure 3 shows the area with value less than 0.05 after the examination. Set 18 adjacent voxels (edge connection, the default value of SPM) and voxel number of 200 as the cluster size parameter, perform cluster reading on these regions, and the region with voxel greater than 70 is selected from the cluster for display, since the number of regions is too large. The positive value is the region with significant difference between FSST and FFT. Negative values represent areas where there is a significant difference between FFT and FSST. As shown in Figure 2, the FSST and FFT also show significant differences in Figure 3.

Based on the analysis and the performance of the actual effect, the significance of the algorithm improvement is proved. The performance of ALFF index has been greatly improved which promotes the application of ALFF and its value.

4. Conclusions

At present, the research on human brain occupational plasticity has been developed rapidly, but the research on ALFF needs to be further improved. Because of the complexity and variety of EEG signals, extracting the information of the signals in time series can effectively reduce the signal ignorance caused by frequency band concentration.

Based on the previous research, this paper proposes the FSST algorithm to improve the resolution and enhance the analytical accuracy of ALFF, aiming to improve the visualization effect of ALFF and the extraction effect of low frequency signals. From the experimental results, it is proved that there is a significant difference between the overall analytical effects of FSST and FFT in resting state, and the extraction precision of FSST is significantly higher, which can reflect the TF structure characteristics of brain waves clearly.

In addition, FSST can detect the changes of brain dynamics in time series. After the processing, the TF map of the signal can be recovered to a greater extent. Based on the accuracy of the recovered signal and the characteristics of low frequency attributes, it is proved that FSST algorithm is feasible for the detection and analysis of ALFF in the time domain. In the future, we will verify the applicability of FSST to ALFF in task-related experiments.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research was supported by the Natural Science Foundation of the Jiangsu Higher Education Institutions of China under Grants 19KJB160001 and 19KJB520004, the Lianyungang city Haiyan project under Grant 2019-QD-004, the MOE Key Laboratory of TianQin Project, Sun Yat-sen University, and the Open Fund Project of Jiangsu Institute of Marine Resources Development under Grant JSIMR202018. The authors would like to thank Qin Yu for his contribution at the early phase of this work.