Abstract

The study focused on the application value of ultrasound Doppler fetal heart rate detection algorithm based on short-time Fourier transform (STFT) in the diagnosis of twin-to-twin selective intrauterine growth restriction (sIUGR) and the correlation between twin-to-twin sIUGR and sFas/sFasL levels of umbilical cord blood. The normalized method was introduced into the STFT algorithm to optimize it to detect the fetal instantaneous heart rate. 82 pregnant women with twin pregnancies were selected as the research subjects and they were divided into the restricted group (41 cases) and the control group (41 cases) according to whether the fetus had selective growth restriction. The two groups were compared for the differences in the fetal mortality, complication rate, and sFas/sFasL expression levels. The results showed that the STFT-based ultrasonic Doppler fetal heart rate detection algorithm could ensure the quality of the fetal heart rate signal and had high resolution at the 200–400 Hz characteristic frequency band and that the accuracy in distinguishing S1 and S2 was 5.8% higher than that of the traditional autocorrelation algorithm. The proportion of abnormal fetal heart rate in the restricted group was significantly higher than that in the control group (), birth weight was significantly lower than that in the control group (), and fetal mortality was significantly higher than that in the control group (). There was no statistical difference in the incidence of complications between the two groups (). In restricted group, the content of sFas in cord blood was (3326.54 ± 317.42) pg/mL and that in the control group was (2003.29 ± 196.45) pg/mL. The content of sFas in cord blood of the restricted group was significantly higher than that of the control group (). In the restricted group, the content of sFasL in cord blood was (382.52 ± 36.17) pg/mL, and that in the control group was (180.84 ± 16.20) pg/mL. The content of sFasL in cord blood of the restricted group was significantly higher than that of the control group (). It was concluded that the STFT-based ultrasound Doppler fetal heart rate monitoring is beneficial to early diagnosis and timely intervention of twin-to-twin sIUGR.

1. Introduction

With the continuous development of assisted reproductive technology, the incidence of twin-to-twin pregnancy has shown a significant upward trend, from 1/80 of natural pregnancy to about 20% [1]. The intrauterine growth environment of twin-to-twin pregnancy is more complicated than a single pregnancy, and the incidence of premature delivery, intrauterine growth restriction, perinatal mortality, and congenital turbulence is significantly increased [2]. Fetal growth restriction (FGR) is one of the important complications of the perinatal period, with an incidence rate of 2.75%∼15.53%, and the domestic average FGR incidence rate is 6.39% [3]. FGR has a certain impact on fetal development and later physical and intellectual development, and its perinatal mortality rate is about 5 times that of normal fetuses [4], so it is of great significance to explore the pathogenesis of FGR. There are many and complex causes of FGR. The Fas/FasL system is mainly responsible for inducing cell apoptosis. The current research results show that Fas and FasL are expressed in placental trophoblasts and decidua cells of full-term parturients. Fas and FasL expression changes after pregnancy, which has a certain protective effect on the fetus. Studies have analyzed the Fas/FasL expression in the placenta of FGR patients and normal pregnant women. The results found that the Fas expression of FGR patients was significantly higher than that of normal pregnant women, while FasL expression was significantly lower than that of normal pregnant women, indicating that Fas and FasL expression levels are associated with the occurrence of FGR [5]. However, there are few studies on the correlation between Fas/FasL and twin-to-twin selective intrauterine growth restriction (sIUGR).

Fetal heart rate monitoring is a common clinical method to ensure the safety of mothers and fetuses. The traditional fetal heart rate monitoring uses an autocorrelation algorithm to calculate the fetus’s instantaneous heart rate, but it is susceptible to interference from maternal and fetal activities, sensors, and the external environment, so it is difficult to accurately measure the fetus’s instantaneous heart rate [6]. Short-time Fourier transform (STFT) is a commonly used time-frequency analysis method, which is characterized by convenient analysis and is used for frequency positioning of sound signals [7]. Some researchers have applied STFT to the processing of fetal heart sound signals by ultrasonic Doppler. The results showed that this algorithm can overcome the weakness of fetal heart sound easily disturbed by noise, realize the real-time calculation of instantaneous heart rate of fetal heart sound, and improve the accuracy of fetal heart rate calculation. However, STFT cannot achieve high resolution in both frequency domain and time domain, so it needs to be further optimized.

Above, there are currently few studies on the correlation between Fas/FasL and twin-to-twin SIUGR, and the resolution of the STFT algorithm needs to be further improved. In this study, normalized image processing method was introduced based on STFT algorithm to improve the resolution of fetal heart rate by ultrasound Doppler. Meanwhile, 82 twin-to-twin pregnant women were taken as research subjects to explore the correlation between twin-to-twin sIUGR and the level of sFas/sFasL in cord blood, expected to provide a reference for the diagnosis and treatment of twin-to-twin sIUGR.

2. Materials and Methods

2.1. Research Subjects and the Grouping

82 pregnant women who were diagnosed as having twin-to-twin selective fetuses in the obstetrics department of our hospital from February 2017 to December 2020 were selected as the research subjects, aged 25–39 years old, with an average age of (28.33 ± 4.53) years. According to whether the fetus has selective growth restriction, they were divided into the restricted group (41 cases) and control group (41 cases). Inclusion criteria for this study: (I) patients not accompanied by structural or chromosomal abnormalities; (II) all pregnancies were of natural fertilization. The exclusion criteria: (I) the presence of typical twin-to-twin transfusion syndrome; (II) those with fetal malformations during the prenatal examination. The study has been approved by the ethics committee of the hospital, and all subjects included in the study had signed an informed consent form.

2.2. Establishment of Optimized STFT Algorithm

The STFT establishes a bridge from the time domain to the frequency domain [8], and it is a commonly used time-frequency analysis method. It has the characteristics of strong matching ability, stable extraction of features, and high matching accuracy [9]. Assuming that the analyzed signal is , , and the analysis window function is , then, for nonstationary signals , the STFT algorithm can be expressed as follows:where represents the analysis window shifted from time 0 to time t. By windowing and truncating the original signal through , a short-term signal is obtained.

Assuming that the duration of the analysis window is , the time of time-domain windowing and truncation is . Then, the STFT at can be expressed as follows:where is the spectrum of the analysis window; and represents the spectrum of the analysis window shifted from to a; by windowing and truncating the signal spectrum through , the short-frequency signal is calculated as follows:

Fourier transform is performed on the continuous signal and windowing is performed on the result of the continuous Fourier transform. Then, the inverse Fourier transform algorithm is expressed as follows:

According to the algorithm, STFT can be expressed as follows:

To completely reconstruct the signal, the STFT should meet the following conditions:

When , the STFT complete reconstruction condition can be expressed as follows:

To realize discrete STFT, it is necessary to sample the STFT in the digital frequency domain and the corresponding STFT can be expressed as follows:where represents the analysis window.

The frequency variable is continuous and the frequency variable needs to be discretized [10]. Then, the discrete STFT of the digital signal can be expressed as follows:

Then, the discrete STFT can be expressed as follows:

For nonstationary signals , , the window sequence is . Then, the sequence STFT algorithm transformation can be expressed as follows:where represents a short-time sequence, and the inverse transformation of at time n can be expressed as follows:

Suppose , let , and then the sequence STFT inverse transformation can be expressed as follows:

During the processing of the ultrasonic Doppler signal, the spectral shape and performance of the window function have a significant impact on the spectrum leakage and fence effect [11]. The rectangular window belongs to the zeroth power window of the time variable [12] and the discretized rectangular window of length N can be expressed as follows:

Then, its discretization form can be expressed as follows:

The triangular window is the first-order form of the power window [13] and the discretized rectangular window of length N (even number) is as follows:

Then, its discretization form can be expressed as follows:

The time domain of the discrete Hanning window of length N is expressed as follows:

The time-domain discretization form of the Hanning window can be expressed as follows:

The time-domain expression of the discrete Hamming window of length N is as follows:

The time-domain discretization form of Hamming window can be expressed as follows:where is the spectral function of the rectangular window.

Normalized cross-correlation (NCC) is one of the robust matching methods. It has the characteristics of strong antibackground noise, strong adaptability to image distortion, and high matching accuracy [14]. In the study, the normalization algorithm is introduced into STFT to optimize it. The normalized correlation coefficient between the feature template M and the two-dimensional time-frequency diagram S in the two-dimensional time-frequency plane can be expressed as follows:where is the mean value of the template image, means the mean value of the area of with as the starting point, and the template size is .

The one-dimensional vectors of the template M and the matched image area can be expressed as follows:where is the template one-dimensional vector and is the one-dimensional vector of the image area to be matched. Then, the template row vector can be expressed as , the image area vector to be matched can be expressed as , and the normalized cross-correlation coefficient can be further expressed as follows:where represents the inner product of and ; and and represent the 2 norms of the sum.

In this study, the STFT algorithm is optimized by introducing the normalized algorithm. The flowchart of the fetal heart rate detection algorithm is shown in Figure 1.

2.3. Algorithm Operating Environment and sIUGR Evaluation Criteria

The optimized STFT algorithm to detect fetal instantaneous heart rate is programmed in MATLAB 7.11.0 and the operating environment is Windows XP. The sampling frequency of heart sounds is qq025 Hz, 14-bit A/D conversion is adopted, and the sampling data is transmitted to the upper computer through serial communication.

SIUGR diagnostic criteria [15]: the estimated weight of a fetus of twins is lower than the 10% of the fetus at normal gestational weeks and the twin-to-twin weight difference is more than 25%.

2.4. sFas/sFasL Level in Umbilical Cord Blood

In the second trimester of pregnancy, the two groups of pregnant women, 2 mL of umbilical cord blood was drawn under the guidance of B ultrasound, anticoagulated using heparin, and centrifuged to obtain plasma, which was stored at −20°C. The sFas/sFasL levels in the cord blood of the two groups of pregnant women were detected according to the instructions of the ELISA test kit (R&D System, USA).

2.5. Statistical Methods

The experimental data were processed using SPSS 19.0 statistical software, and the count data were tested by χ2, and indicates that the difference was statistically significant.

3. Results

3.1. The Optimized STFT Algorithm

Figure 2 shows the initial fetal heart sound data collected by Doppler ultrasound. This one-dimensional data cannot effectively display all the features of fetal heart sounds while the two-dimensional signal can, and it has a high resolution at 200∼400 Hz, so in the follow-up research, the frequency band of 200∼400 Hz was used.

After the original fetal heart rate signal map was normalized, the normalized correlation curve of the feature template and the target image is shown in Figure 3. It was noted that the relationship numbers were all above 0.86 at different times.

3.2. Fetal Heart Rate Detected by the Optimized STFT Algorithm

The fetal heart sound signals collected by Doppler ultrasound (Figure 4) were extracted by the autocorrelation algorithm and the optimized STFT algorithm, respectively. For the fetal heart signal with a small amplitude, autocorrelation cannot accurately locate it (Figure 5), only showing an abnormal fetal heart rate. The instantaneous heart rate extracted by the optimized STFT algorithm can accurately locate the fetal heart sound signal with a small amplitude and can correctly display the instantaneous fetal heart rate waveform (Figure 6).

The segmentation result of the fetal heart sound signal by the optimized STFT algorithm is shown in Figure 7. It was noted that the STFT optimization algorithm can accurately segment the S1 and S2 sounds of the fetal heart sound.

3.3. Fetal Heart Sound Detection Results by the Optimized STFT Algorithm

Figure 8 shows the detection results of clinically collected fetal heart sounds. For fetuses at 30–40 weeks of gestation, the autocorrelation algorithm [16] and optimized STFT algorithm were used to detect fetal heart sounds. The results showed that the detection error rate of the autocorrelation algorithm at 32 weeks and 38 weeks was 4% and 2%, respectively; those of the optimized STFT were 1% and 0, respectively, and the total detection error rate was reduced by 5%.

3.4. Comparison of Basic Clinical Data of Two Groups of Pregnant Women

The restricted group and the control group were compared for the age, gestational age at delivery, and birth weight of the fetus (Figure 9). The birth weight of the fetus in the restricted group was significantly lower than that of the control group (). There was no statistical difference in gestational weeks ().

3.5. Analysis of Abnormal Results of Fetal Heart Rate Detection in the Two Groups

Figure 10 shows the abnormal fetal heart rate detection in the restricted group and the control group. There were 30 cases of the abnormal fetal heart rate in the restricted group, accounting for 73.17%; there were 17 cases of abnormal fetal heart rate in the control group, accounting for 41.46%. The proportion of fetal heart rate abnormalities in the restricted group was significantly higher than that in the control group and there was a statistical difference between the two ().

3.6. Comparison of Perinatal Outcomes between the Two Groups of Patients

Figure 11 shows the perinatal outcomes of the restricted group and the control group. There were 22 deaths and 17 survivors in the restricted group, accounting for 53.66% and 41.46%, respectively. There were 10 deaths and 32 survivors in the control group, accounting for 24.39% and 78.09%, respectively. The survival rate and mortality of the two groups were significantly different (). There was no statistical difference in the incidence of complications between the two groups ().

3.7. Analysis of sFas/sFasL Expression Levels in Cord Blood of the Two Groups

Figure 12 shows the sFas and sFasL expression levels. The sFas content of the restricted group was (3326.54 ± 317.42) pg/mL and the sFas content of the control group was (2003.29 ± 196.45) pg/mL. The sFas content of the restricted group was significantly higher than that of the control group (). The sFasL content of the restricted group was (382.52 ± 36.17) pg/mL and the sFasL content of the control group was (150.84 ± 16.20) pg/mL. The sFasL content of the restricted group was significantly higher than that of the control group ().

4. Discussion

In this study, the normalized image processing method was introduced into the STFT algorithm to optimize it to detect the fetal instantaneous heart rate. It was found that the cross-correlation coefficients were all above 0.86. The normalized cross-correlation algorithm is to calculate the correlation between vectors according to the row/column. A higher correlation coefficient indicates the correlation between the two images is better [17]. In the study, the correlation coefficients were all above 0.86, indicating that the normalized image processing method had a significant advantage in instantaneous heart rate extraction. The normal heart mainly includes four heart sounds S1∼S4. In most cases, only two heart sounds S1 and S2 can be monitored by Doppler ultrasound [18]. In a cardiac cycle, the time characteristics of heart sounds can be fully expressed by the two heart sounds in systole S1 and diastolic S2 [19]. For normal fetuses, the normal heart rate is 120 to 160 beats/min [20]. The optimized STFT algorithm in this study can accurately segment the S1 and S2 fetal heart sounds. The detection error rates of the traditional autocorrelation algorithm at 32 weeks and 38 weeks were 4% and 2%, respectively. The detection error rates of the optimized STFT algorithm were 1% and 0, respectively, and the total detection error rate was reduced by 5%. These results showed that the optimized STFT algorithm had a significant advantage in instantaneous heart rate detection. It may be because the introduction of the normalized cross-correlation algorithm strengthens the accurate positioning performance of the fetal heart sound signal with a smaller amplitude and ultimately improves the calculation accuracy of fetal heart rates.

Fetal growth restriction is one of the important complications in the perinatal period. sIUGR can occur at any period of pregnancy and is one of the common complications in twin-to-twin pregnancy, mostly caused by congenital abnormalities. The results of this study found that the birth weight of the fetus in the restricted group was significantly lower than that in the control group (), indicating that the weight of the fetus in the restricted group was reduced. It may be related to the expected nutritional abundance, congenital abnormal diseases, and poor fetal development [21]. At the same time, the proportion of abnormal fetal heart rate in the restricted group was 73.17%, the proportion of abnormal fetal heart rate in the control group was 41.46%, and the proportion of abnormal fetal heart rate in the restricted group was significantly higher than that in the control group (), indicating that abnormal fetal heart rate can provide a reference for the diagnosis of twin-to-twin sIUGR. The abnormal fetal heart rate in the control group may be caused by the abnormal fetal blood supply and oxygen supply [22]. In the restricted group, there were 22 deaths and 17 survivors, accounting for 53.66% and 41.46%, respectively. In the control group, there were 10 deaths and 32 survivors, accounting for 24.39% and 78.09%, respectively. The survival rate and mortality of the two groups were significantly different (). There was no statistical difference in the incidence of complications between the two groups (). It showed that selective fetal restriction increased the mortality and survival rate of the fetus, so early ultrasound Doppler monitoring of the fetal heart rate is of great significance. The results in the study found that the sFas content of the restricted group was (3326.54 ± 317.42) pg/mL and the sFas content of the control group was (2003.29 ± 196.45) pg/mL. The sFas content of the restricted group was significantly higher than that of the control group (). The sFasL content of the restricted group was (382.52 ± 36.17) pg/mL and the sFasL content of the control group was (180.84 ± 16.20) pg/mL. The sFasL content of the restricted group was significantly higher than that of the control group (). sFas is a transmembrane glycoprotein, widely distributed in various tissues and organs of the human body, and plays an important role in cell apoptosis [23]. sFasL is the ligand of sFas, a member of the tumor necrosis factor superfamily, and is mainly expressed in activated lymphocytes [24]. The sFas/sFasL system can play a role in immune tolerance by mediating cell apoptosis, which is related to placental apoptosis [25]. The research of Karowicz et al. [26] pointed out that the expression of Fas in maternal blood of pregnant women with pregnancy-induced hypertension combined with sIUGR was significantly increased, the level of Fas of the fetus at birth was also significantly increased, and the expression of Fas in maternal blood was higher than that in cord blood (), consistent with the results of this study. Above, the changes in sFas/sFasL levels are correlated with the occurrence of sIUGR.

5. Conclusion

In the study, the normalized method was introduced into the STFT algorithm to optimize it to detect the fetal instantaneous heart rate and analyze the differences in the fetal mortality, complication rate, and cord blood sFas/sFasL expression levels between the restricted group and the control group. It was found that the optimized STFT in this study elevated the accuracy of fetal heart rate detection. The levels of sFas and sFasL had a certain correlation with twin-to-twin sIUGR. However, some limitations in the study should be noted. The sample size is small, the research subjects of this study have geographical limitations, which will reduce the power of the study. In the follow-up, an expanded sample size is necessary to strengthen the findings of the study. In conclusion, prenatal STFT-based ultrasound Doppler monitoring of fetal heart rate is conducive to early diagnosis and timely intervention in twin-to-twin sIUGR. The results that cell apoptosis increases in twin-to-twin sIUGR fetuses provide a reference for the diagnosis and treatment of twin-to-twin sIUGR.

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 no conflicts of interest.