1.一種基于分段濾波處理的子宮肌電強(qiáng)度信息提取方法,其特征在于,包括以下步驟:
步驟A1、采集孕產(chǎn)婦的子宮肌電信號,記為原始子宮肌電信號;
步驟A2、識別上述子宮肌電信號的爆發(fā)波段與非爆發(fā)波段;
步驟A3、對所述子宮肌電信號不同信號波段進(jìn)行不同方式的濾波處理,即將子宮肌電信號劃分為爆發(fā)波段與非爆發(fā)波段,并分段進(jìn)行不同的濾波處理,得到濾波處理后的子宮肌電信號消除子宮肌電信號中的脈沖性噪聲干擾;
步驟A4、利用包絡(luò)提取算法提取所述子宮肌電強(qiáng)度信息。
2.根據(jù)權(quán)利要求1所述基于分段濾波處理的子宮肌電強(qiáng)度信息提取方法,其特征在于,所述步驟A1中,所述子宮肌電信號是通過貼附在孕產(chǎn)婦腹部的Ag-AgCl表面電極采集得到。
3.根據(jù)權(quán)利要求1所述基于分段濾波處理的子宮肌電強(qiáng)度信息提取方法,其特征在于,所述步驟A2中,識別爆發(fā)波段與非爆發(fā)波段的方法,包括以下步驟:
步驟B1、對所述子宮肌電信號進(jìn)行預(yù)處理;
步驟B2、由預(yù)處理后子宮肌電信號特點將其分為爆發(fā)波段與非爆發(fā)波段。
4.根據(jù)權(quán)利要求3所述基于分段濾波處理的子宮肌電強(qiáng)度信息提取方法,其特征在于,所述步驟B1中,對于采集到的原始子宮肌電信號預(yù)處理包括如下步驟:
步驟C1、對采集到的子宮肌電信號進(jìn)行降采樣;
步驟C2、計算上述降采樣后信號中待求樣本點所在窗口的過零率。
5.根據(jù)權(quán)利要求3所述基于分段濾波處理的子宮肌電強(qiáng)度信息提取方法,其特征在于,所述步驟B2中,子宮肌電信號爆發(fā)波與非爆發(fā)波段的判斷方法包括:
步驟D1、設(shè)定子宮肌電信號爆發(fā)波段過零率閾值;
步驟D2、比較權(quán)利要求3中計算待求樣本點所在窗口的過零率,若其過零率大于等于所述爆發(fā)波過零率設(shè)定閾值,則判斷該點處于爆發(fā)波段;反之處于非爆發(fā)波段。
6.根據(jù)權(quán)利要求2所述的基于分段濾波處理的子宮肌電強(qiáng)度信息提取方法,其特征在于,所述子宮肌電信號爆發(fā)波段與非爆發(fā)波段的判斷方法還包括:
步驟E1、計算子宮肌電信號基線;
步驟E2、通過判斷子宮肌電信號波形與子宮肌電信號基線之間的位置關(guān)系得出爆發(fā)波段與非爆發(fā)波段信號,當(dāng)子宮肌電信號波形處于該信號基線上方時,則認(rèn)為該信號處于爆發(fā)波段;反之處于非爆發(fā)波段。
7.根據(jù)權(quán)利要求1所述的基于分段濾波處理的子宮肌電強(qiáng)度信息提取方法,其特征在于,所述步驟A3中,對子宮肌電信號不同信號波段進(jìn)行不同方式的濾波處理,方法如下:
在相同濾波函數(shù)的條件下,對子宮肌電信號的爆發(fā)波段與非爆發(fā)波段進(jìn)行不同長度的濾波窗口處理,即實現(xiàn)對子宮肌電信號進(jìn)行分段濾波處理,得到濾波處理后的子宮肌電信號。
8.根據(jù)權(quán)利要求1所述的基于分段濾波處理的子宮肌電強(qiáng)度信息提取方法,其特征在于,所述步驟A3中,對子宮肌電信號不同信號波段進(jìn)行不同方式的濾波處理,方法還包括:
對子宮肌電信號的爆發(fā)波段與非爆發(fā)波段分別使用不同長度的濾波窗口與不同濾波函數(shù)的處理,得到濾波處理后的子宮肌電信號。
9.根據(jù)權(quán)利要求1所述的基于分段濾波處理的子宮肌電強(qiáng)度信息提取方法,其特征在于,步驟A4中,使用包絡(luò)提取算法從分段濾波處理后的子宮肌電信號中提取子宮肌電強(qiáng)度信號,具體為:
步驟F1、使用包絡(luò)提取算法提取子宮肌電強(qiáng)度信號;
步驟F2、為了抵消包絡(luò)提取過程中可能出現(xiàn)的頻譜泄露,對所述子宮肌電信號包絡(luò)進(jìn)行漢寧窗處理;
所述漢寧窗計算公式如下:
其中xi為輸入信號,即第i點的子宮肌電強(qiáng)度信號幅值,yi為當(dāng)前點對應(yīng)的輸出信號值,N為窗口長度。
10.根據(jù)權(quán)利要求9所述的基于分段濾波處理的子宮肌電強(qiáng)度信息提取方法,其特征在于,所述步驟F1中,采用RMS法提取子宮肌電強(qiáng)度信號,在計算子宮肌電強(qiáng)度信號時需要選擇合適的窗口長度LRMS,且在信號的不同位置,其窗口長度選擇也會存在差異。