CN108013872A - System for maternal fetus rhythm of the heart - Google Patents

System for maternal fetus rhythm of the heart Download PDF

Info

Publication number
CN108013872A
CN108013872A CN201810021545.9A CN201810021545A CN108013872A CN 108013872 A CN108013872 A CN 108013872A CN 201810021545 A CN201810021545 A CN 201810021545A CN 108013872 A CN108013872 A CN 108013872A
Authority
CN
China
Prior art keywords
module
signal
mrow
mother
parent
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201810021545.9A
Other languages
Chinese (zh)
Inventor
王妍
陈坤
赵扬玉
陈练
王汉斌
张嘉宾
张珏
方竞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Peking University
Peking University Third Hospital
Original Assignee
Peking University
Peking University Third Hospital
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Peking University, Peking University Third Hospital filed Critical Peking University
Priority to CN201810021545.9A priority Critical patent/CN108013872A/en
Publication of CN108013872A publication Critical patent/CN108013872A/en
Pending legal-status Critical Current

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
    • A61B5/0245—Detecting, measuring or recording pulse rate or heart rate by using sensing means generating electric signals, i.e. ECG signals
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02411—Detecting, measuring or recording pulse rate or heart rate of foetuses
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/43—Detecting, measuring or recording for evaluating the reproductive systems
    • A61B5/4306—Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
    • A61B5/4343—Pregnancy and labour monitoring, e.g. for labour onset detection
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Cardiology (AREA)
  • Signal Processing (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Gynecology & Obstetrics (AREA)
  • Pregnancy & Childbirth (AREA)
  • Reproductive Health (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

The invention discloses a kind of system for maternal fetus rhythm of the heart, including signal acquisition module, signal pre-processing module, signal quality assessment module, mother's electrocardio extraction module, mother's heart rate detection module, mother's electrocardio suppression module, Fetal ECG extraction module, fetal heart frequency detection module and result output module.Wherein：Signal acquisition module gathers her abdominal electric signal, signal pre-processing module removes power frequency, baseline drift and myoelectricity interference, signal quality assessment module rejects low-quality signal, mother's electrocardio extraction module and mother's heart rate detection module calculate mother's heart rate, mother's electrocardio suppression module, Fetal ECG extraction module and fetal heart frequency detection module calculate fetal heart frequency, and as a result output module exports mother's heart rate and fetal heart frequency at the same time.The system rejects the low-quality signal in monitoring process and fully suppresses the noise signal before detection Fetal Heart Rate, can accurately extract fetal heart frequency with less cardiac diagnosis lead, meet the needs of clinical long-range fetal heart monitoring.

Description

System for maternal fetus rhythm of the heart
Technical field
It is more particularly to a kind of to be for maternal fetus rhythm of the heart the invention belongs to electrocardiogram (ECG) data processing technology field System.
Background technology
Electronic fatal monitoring (EFM) is it is currently understood that fetal in utero situation and the important method of fetus reserve capabillity.Utilize Fetal Heart Rate Monitoring curve judges fetal in utero whether there is hypoxic-ischemic state, and the active situation of fetal central nervous system, is Clinical judgment Perinatal morbidity and guiding treatment provide important evidence.Fetal rhythm is supervised when being divided into antenatal and production according to the time of application Shield.Antepartum monitoring is mostly non-stress test (NST), refers to monitor i.e. antenatal anodinia before childbirth, stimulates feelings without external load Under condition, observation and the record of Fetal Heart Rate are carried out to fetus, to understand fetus reserve capabillity；Monitoring during labor is uterine contraction irritant test (CST), as there is the fetal heart monitoring of uterine contraction, meet change for understanding placenta transient anoxic when the uterine contraction, measure fetus Reserve function.
Existing scheme mainly has following：The country is widely used that Doppler electronic fetal heart monitoring, it is necessary to instrument active Send ultrasonic signal and act on fetus, extract Fetal Heart Rate in the different reflected signals from the various tissues of fetus to ultrasonic wave, this side Method is sensitive to lie, and pregnant woman need to use fixed posture in monitoring process, sitting posture or lies on one's side, and is not suitable for carrying out long-time prison Survey；Another scheme is that the electrocardiosignal gathered from her abdominal surface by multi-channel electrode extracts Fetal Heart Rate.This method It is adapted to long-term monitoring Fetal Heart Rate, but has power frequency, breathing, mother's electrocardio from the frequent aliasing of electric signal of her abdominal acquisition surface With the interference of the various noises such as myoelectricity.
Chinese patent invention 201410164558.3 proposes a kind of female fetal electrocardiogram separation method, gathers the electrocardio letter of multichannel Number, determined through SVD decomposition algorithms in multiline message after active ingredient, it is female to multichannel using probability independent component analysis (ICA) Tire signal mixed signal carries out isolated Fetal ECG, parent electrocardio and its interference, then shows Fetal ECG and parent Electrocardio.The advantages of program is：Clear model, it is readily appreciated that.But the program has the drawback that：Probability independent component analysis is managed Observation variable number is required to be more than or equal to independent variable number on, the program does not pre-process known noise fully, right Electrocardiogram acquisition number of active lanes requires height, can paste many electrodes in her abdominal and (lead her abdominal mixing in the embodiment for six to imitate True signal) skin can be caused uncomfortable and increase monitoring cost.
The content of the invention
It is an object of the invention to the deficiency for above-mentioned prior art, there is provided a kind of for maternal fetus rhythm of the heart System, the system only gather triple channel parent belly electric signal, can be under conditions of low signal-to-noise ratio, and stabilization extracts mother exactly Heart rate and fetal heart frequency, meet the needs of clinical long-range fetal heart monitoring.
To achieve these goals, the system for maternal fetus rhythm of the heart of the invention, including：
Signal acquisition module (1), gathers triple channel parent belly electric signal, and collection signal is amplified, AD digital-to-analogues Conversion；
Signal pre-processing module (2), locates the parent belly electric signal that signal acquisition module (1) collects in advance respectively Reason, removes industrial frequency noise, baseline drift and myoelectricity interference；
Signal quality assessment module (3), quality evaluation is carried out to the signal after signal pre-processing module (2) processing, and is carried Useful signal is taken to be output to parent electrocardio extraction module (4) and parent electrocardio suppression module (6)；
Parent electrocardio extraction module (4), adopts effective parent belly electric signal after signal quality assessment module (3) processing Separated with independent composition analysis algorithm；
Maternal heart rate detection module (5), choose parent electrocardio extraction module (4) the clean parent electrocardio signal isolated into Row R ripples detect to obtain maternal heart rate time series and are output to parent electrocardio suppression module (6) and result output module (9)；
Parent electrocardio suppression module (6), effective parent belly electric signal after received signal quality evaluation module (3) processing The maternal heart rate time series transmitted with maternal heart rate detection module (5), removes each passage parent electrocardio component；
Fetal ECG extraction module (7), receives the parent belly electric signal after parent electrocardio suppression module (6) processing and adopts Separated with independent composition analysis algorithm；
Fetal heart frequency detection module (8), choose Fetal ECG extraction module (7) separation after clean Fetal ECG signal into Row R ripples detect to obtain fetal heart frequency time series and are output to result output module (9)；
As a result output module (9), maternal heart rate detection module (5) and fetal heart frequency detection module (8) module are transmitted Heart rate time sequence exports at the same time.
Preferably, window of the signal pre-processing module (2) by 10~60s of duration, is moved with each sliding window The mode of 50% window width chooses the parent belly electric signal that signal acquisition module (1) collects, and can so export current prison in real time Survey the result of calculation of data.
Preferably, the signal pre-processing module (2) removes industrial frequency noise and baseline drift using Bezier Wave filter, because the phase response of Bessel filter is almost linear, can effectively avoid phase distortion.
Locally thrown for phase space preferably, the signal pre-processing module (2) removes the method that myoelectricity interference uses Shadow method, since myoelectricity and the cardiac electrical frequency spectrum of belly have a large amount of aliasings, the in-band frequency that traditional filtering algorithm suppresses echo signal can It can change the dynamics of former time series, dynamics can effectively be kept using phase space partial projection method, specifically Step is：
1) phase space reconfiguration：According to C-C algorithms, calculate and remove industrial frequency noise and baseline drift postabdomen electric signal phase space One-dimensional belly electric signal, is reconstructed into the matrix of m dimensions by the optimal embedding dimension m and optimum delay time τ of reconstruct；
2) neighborhood is determined：The optimal radius of neighbourhood ε of partial projection is calculated using recursive analysis, determines neighborhood；
3) non-orthogonal projection and bob-weight structure：In neighborhood, the covariance square of the neighborhood square formed by calculating partial vector Battle array, estimates different noise component(s)s, calculates correction value, and iteration updates the value of observation station；By bob-weight structure, by m dimension matrix conversions For the belly electrocardiosignal after denoising.
Preferably, described signal quality assessment module (3) specific algorithm is based on mother and fetus mixing electrocardio work( Rate Spectral structure, since long-term monitoring parent belly electric signal electrode delamination that may be present, the big, movement of Skin Resistance etc. are disturbed, The valuation mistake of heart rate is caused, so as to cause the false alarm of patient monitor, has seriously affected clinician to true critical alarm Response and monitoring effect, and mother and fetus mixing electrocardio QRS complex energy be concentrated mainly on about centered on 10Hz Width be about that the ratio that power spectral density (PSD) value in the frequency band accounts for total PSD values can be used as and judge in the frequency band of 10Hz The reference index of QRS wave quality.Therefore, define the signal quality index (sSQI) based on Power Spectrum Distribution be frequency f be 5~ 14Hz frequency bands account for the ratio of 3~49Hz frequency bands, according to formula:
Wherein, PSD (k, f) is power spectral density, and k is current segment signal sequence number, and f is frequency, and threshold value T is adjusted by testing Whole empirical value, scope are 0.4~0.7, and when the signal sSQI for having a more than passage is labeled as 1, then this period signal is to have Imitate signal；If three passage sSQI are marked and are, the segment signal is judged for exception and exports alarm signal.
Preferably, the parent electrocardio extraction module (4) and Fetal ECG extraction module (7), independent component analysis is calculated Method is fast independent component analysis algorithm (FastICA).Correct application for FastICA, it is necessary to meet claimed below：Signal Source is separate, and signal source non-gaussian (most one), the quantity of the constant hybrid matrix of instantaneous linear and measuring signal is equal Or the quantity more than source.In practical situations, last hypothesis is not fully met, because industrial frequency noise caused by alternating current, Myoelectricity interference caused by baseline drift caused by respiratory movement and contraction of muscle adds the quantity in source and causes one linearly The hybrid matrix of change, but mother's electrocardio is that triple channel belly electric signal is most strong and more common independent source, Signal Pretreatment mould Block (2) effectively avoids phase distortion when industrial frequency noise, baseline drift is removed, and is remained while myoelectricity interference is removed The dynamics of mother's electrocardio and Fetal ECG, then can be readily separated out with FastICA；Similarly, Signal Pretreatment mould Block (2) and parent electrocardio suppression module (6) main purpose are to increase the signal-to-noise ratio of Fetal ECG in three passages, are improved The accuracy of FastICA isolating fetal electrocardio components.
Preferably, in the parent electrocardio suppression module (6), concretely comprise the following steps：
1) the average value RR of mother's cardiac RR intervals is calculatedmeanWith R ripple position sequences：By maternal heart rate detection module (5) Mother's heart rate time sequence transmitted changes into RR interval series, and calculating average value is RRmean, then the R ripple position sequences that add up to obtain；
2) the QRS wave matrix of each passage is constructed respectively：To each passage after signal quality assessment module (3) processing The each mother's ecg-r wave position of effective parent belly electric signal preceding 0.3*RRmeanWith rear 0.7*RRmeanIt is a length of during selection RRmeanSignal segment for column vector construct QRS wave matrix X, wherein " * " is multiplication；
3) SVD decomposition is carried out to the matrix X of each passage respectively：According to formula X=U Σ VT, herein U for it is left it is unusual to Amount, Σ is 0 except cornerwise element, and the element on diagonal arranges from big to small, is known as singular value, VT (transposition of V) For right singular vector；Mother's electrocardio component can be represented by rule of thumb taking preceding 3~5 diagonal element reconstruction signals M of Σ；
4) respective M signals are individually subtracted in each passage parent belly electric signal and suppress mother's electrocardio component.
The invention has the advantages that：
The system has a clear superiority in long-range maternal fetus rhythm of the heart task, and the cardiac diagnosis lead long-range of triple channel is adopted Collection, signal quality from lie influence；Existing method requires her abdominal electrical signal collection number of active lanes high and pre- The myoelectricity interference of paroxysmal is not considered during processing or only removes noise with low-pass filtering, causes the dynamics of signal special Sexually revise and then increase later stage separated difficulty.The system is abundant before Fetal ECG extraction module (7) carries out Signal separator Suppress industrial frequency noise, baseline drift, myoelectricity interference and mother's electrocardio component and be effectively retained the dynamics spy of Fetal ECG signal Property, thus can use triple channel cardiac diagnosis lead long-range collection her abdominal electric signal, it is more stables, exactly extraction the fetus heart Rate, improves the comfort level of monitoring process；The system signal quality assessment module reject long-range monitoring electrocardiogram (ECG) data there may be Electrode delamination, the low-quality signal caused by reason such as the big, movement of Skin Resistance so that the result of long-range fetal heart monitoring more may be used Lean on.
Brief description of the drawings
The structure diagram of Fig. 1 present invention；
Fig. 2 triple channel parent belly electric signals；
Fig. 3 parent belly electric signals after signal pre-processing module denoising；
Embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, with reference to the accompanying drawings and embodiments, it is right The present invention is further elaborated.
As shown in Figure 1, the embodiment of the present invention provides a kind of system for maternal fetus rhythm of the heart, including：
1) signal acquisition module (1), threeway is picked up using wet gel electrode slice in palace bottom, pubic symphysis and palace bottom lower left side Road her abdominal electrocardiosignal, and collection signal is amplified and AD digital-to-analogue conversions obtain signal X, as shown in Figure 2；
2) signal pre-processing module (2), receive the triple channel parent belly electrocardiosignal that signal acquisition module (1) collects Sliding window pretreatment is carried out respectively, is removed industrial frequency noise, baseline drift and myoelectricity interference, is comprised the following steps that：
(a) by the window of duration 10s, in a manner of each sliding window moves 5s window widths choose signal acquisition module (1) and adopt The parent belly electric signal collected；
(b) baseline drift is filtered out, is removed using the Bezier high-pass filter that cutoff frequency is 3Hz since respiratory movement is drawn The ECG baseline drift risen；
(c) 50Hz industrial frequency noises are filtered out, removing alternating current using the Bezier bandstop filter that frequency is 48~52Hz causes Industrial frequency noise；
(d) myoelectricity interference is filtered out, phase-space reconstruction delay parameter t is 20ms in the present embodiment, and embedded dimension is 2, according to The geometric properties of chaos attractor, the local tangent space of manifold where CHAOTIC INTERFERENCE is projected to attractor by regional area, then Inverse transformation is to time domain space.Sub-step is as follows:
(d-1) phase space reconfiguration.According to Takens embedding theorems, using C-C algorithms, suitable time delay and embedding is selected Enter dimension, phase space reconstruction so that noisy fetal rhythm electric signal expands to m by One-dimension Time Series and ties up, and is sufficiently spread out attractor Structure.Original fetal rhythm electric signal is the time series that one group of long degree is N, is denoted as, x (i), i=1,2 ..., N, and Embedded dimensions are M, when time delay is τ, obtains the space vector of one group of m dimension, Y1,Y2,...,Yk,...YLMiddle L=N- (m-1) τ, then:
Yi=[x (i), x (i+ τ), x (i+2 τ) ..., x (i+ (m-1) τ)] (3)
(d-2) determine neighborhood, the optimal radius of neighbourhood of partial projection is calculated using recursive analysis.
Calculate recurrence plot Ri,j：
Ri,j=Θ (ε-| | Yi-Yj| |), i, j=1,2 ..., N
In formula, ε is the radius of neighbourhood, and Θ () is Heaviside functions；
Calculate recurrence plot diagonal (not including leading diagonal) length hi:
In formula, t is equal to hiAverage add three times standard deviation；
Definition, Np(ε) is equal to and meets conditionCornerwise number,
Definition,The corresponding ε of β (ε) minimalization is the optimal radius of neighbourhood.
(d-3) non-orthogonal projection.In neighborhoodIt is interior, calculate partial vectorThe neighborhood square A of compositioni, wherein,And non-orthogonal projection is carried out, calculate covariance matrix Ci=(RAi)T(RAi), wherein, T is transposition computing；R is M ties up diagonal square matrix, R11And RmmLarger, preferable R should be taken as far as possible11And Rmm100 are taken, remaining is 1；
(d-4) determine local neighborhood noise subspace, and calculate correction value.Estimate different noise component(s)s, composition is non-just Hand over projection matrixWherein eqFor covariance matrix CiFeature vector.Computed correction θi：
(d-5) update the data a little, determine observation station YiThe updated value of iteration once：
So as to obtain the correction value of the point in whole phase space, the point of whole phase space correspondence position is averaged, bob-weight Structure can obtain the sequence after noise reduction of iteration
(d-6) iteration, repeat the above steps d-2- step d-5, preferably 3-5 rear stopping of iteration.
Signal X1 is obtained after above-mentioned steps are handled, as shown in Figure 3.
3) specific algorithm of signal quality assessment module (3) is to be based on ECG signal Power Spectrum Distribution, the energy of QRS complex waves It is about power spectral density (PSD) value in the frequency band in the frequency band of 10Hz to measure the width being concentrated mainly on about centered on 10Hz The ratio for accounting for total PSD values can be as the reference index for judging QRS wave quality.We calculate frequency f and are accounted for for 5~14Hz frequency bands The ratio of 3~49Hz frequency bands simultaneously defines the signal quality index (sSQI) based on Power Spectrum Distribution and is:
Wherein, PSD (k, f) is power density, and threshold value 0.6 is the empirical value by testing adjustment, when have a passage with On signal sSQI be labeled as 1, then this period signal is useful signal.
4) mother's electrocardio extraction module (4), effective parent belly electrocardio after received signal quality evaluation module (3) processing Signal carries out blind source signal separation.This embodiment blind source signal separation algorithm is specifically FastlCA.Object function selection kurtosis pole Greatly, the step of algorithm is as follows：
(1) centralization is carried out to the Three-channel data X1 of observation, it is 0 to make its average；
(2) albefaction, X1 → Z are carried out to data；
(3) the component number for needing to estimate is set as 3, if value iterations p ← 1；
(4) an initial weight vector W is randomly selectedp；
(5) W is madep=E { Zg (WpTZ)}-E{g'(Wp TZ)}W
Wherein, E [] is mean operation；G () is nonlinear function, takes g (y)=tanh (y)；
(6)(7) W is madep=Wp/||Wp||；
(8) W is worked aspWhen not restraining, then the 5th step is returned；
(9) make p=p+1, if p≤3, return to the 4th step；
(10) the separated result of algorithm
5) mother's heart rate detection module (5) receiving module (4) output Y and by each channel signal basis signal matter in Y Measure evaluation module (3) and calculate SDR, SDR maximums are corresponded to clean parent electrocardio signal carries out R ripple detections, and the heart rate is calculated Method is Pan-Tompkins, it is by analyzing slope, amplitude and width information identification QRS complex waves；The sampling of the present embodiment signal Under conditions of rate is 250Hz, 30 width of window selection are the most suitable (120ms), identify R ripple position sequences, difference obtains RR sequences, between 60 divided by RR the phase obtain mother's heart rate sequence and be output to parent electrocardio suppression module (6) and result output module (9)；
6) mother's electrocardio suppression module (6), effective parent belly electrocardio after received signal quality evaluation module (3) processing The maternal heart rate time series that signal and maternal heart rate detection module (5) transmit, removes the electrocardio component of wherein parent, specific step Suddenly it is：
(a) the average value RR of mother's cardiac RR intervals is calculatedmeanWith R ripple position sequences：By maternal heart rate detection module (5) Mother's heart rate time sequence transmitted changes into RR interval series, and calculating average value is RRmean, then the R ripple position sequences that add up to obtain.
(b) the QRS wave matrix of each passage is constructed respectively：To each passage after signal quality assessment module (3) processing The each mother's ecg-r wave position of effective parent belly electric signal preceding 0.3*RRmeanWith rear 0.7*RRmeanIt is a length of during selection RRmeanSignal segment for column vector construct QRS wave matrix X, wherein " * " is multiplication；
(c) SVD decomposition is carried out to the matrix X of each passage respectively：According to formula X=U Σ VT, herein U for it is left it is unusual to Amount, Σ is 0 except cornerwise element, and the element on diagonal arranges from big to small, is known as singular value, VT(transposition of V) For right singular vector；Mother's electrocardio component can be represented by rule of thumb taking preceding 3 diagonal element reconstruction signals M of Σ；
(d) respective M signals are individually subtracted in each passage parent belly electric signal and suppress mother's electrocardio component, be left Triple channel electric signal is X2.
7) fetus extraction ECG module (7) carries out X2 blind source separating using FastICA algorithms, and algorithm steps are female with step Close electrocardio extraction module (4), algorithm separating resulting are X3.
8) fetal heart frequency detection module (8), receive fetus extraction ECG module (7) and export and comment X3 basis signal quality To estimate module (3) and calculate sSQI, the passage for being to mark carries out R ripple detections respectively, and the algorithm of heart rate is Pan-Tompkins, Under conditions of sample rate is 250Hz, 15 width of window selection are the most suitable (60ms), identify R ripple position sequences, difference RR interval series are obtained, 60 divided by RR interval series obtain Fetal Heart Rate sequence；According to Fetal Heart Rate in normal range (NR) (110- The Fetal Heart Rate sequence that 160bpn) accounting is most or variance is minimum is output to result output module (9)；
9) what result output module (9), reception mother's heart rate detection module (5) and fetal heart frequency detection module (8) transmitted Heart rate sequence, while export result.
It should be noted that in the case where there is no conflict, the feature in embodiment and embodiment in the present invention can phase Mutually combination.Although present invention has been a degree of description, it will be apparent that, do not departing from the bar of the spirit and scope of the present invention Under part, the appropriate change of each condition can be carried out.It is appreciated that the invention is not restricted to the embodiment, and be attributed to right and want The scope asked, it includes the equivalent substitution of each factor.

Claims (6)

  1. A kind of 1. system for maternal fetus rhythm of the heart, it is characterised in that including：
    Signal acquisition module (1), gathers triple channel parent belly electric signal, and collection signal is amplified, AD digital-to-analogue conversions；
    Signal pre-processing module (2), pre-processes the parent belly electric signal that signal acquisition module (1) collects respectively, Remove industrial frequency noise, baseline drift and myoelectricity interference；
    Signal quality assessment module (3), quality evaluation is carried out to the signal after signal pre-processing module (2) processing, and extraction has Signal output is imitated to parent electrocardio extraction module (4) and parent electrocardio suppression module (6)；
    Parent electrocardio extraction module (4), to effective parent belly electric signal after signal quality assessment module (3) processing using only Vertical constituent analysis algorithm is separated；
    Maternal heart rate detection module (5), chooses the clean parent electrocardio signal that parent electrocardio extraction module (4) is isolated and carries out R Ripple detects to obtain maternal heart rate time series and is output to parent electrocardio suppression module (6) and result output module (9)；
    Parent electrocardio suppression module (6), effective parent belly electric signal and mother after received signal quality evaluation module (3) processing The maternal heart rate time series that body heart rate detection module (5) transmits, removes each passage parent electrocardio component；
    Fetal ECG extraction module (7), receives the parent belly electric signal after parent electrocardio suppression module (6) processing and using only Vertical constituent analysis algorithm is separated；
    Fetal heart frequency detection module (8), chooses the clean Fetal ECG signal after Fetal ECG extraction module (7) separation and carries out R Ripple detects to obtain fetal heart frequency time series and is output to result output module (9)；
    As a result output module (9), the heart rate that maternal heart rate detection module (5) and fetal heart frequency detection module (8) module are transmitted Time series exports at the same time.
  2. 2. signal pre-processing module (2) according to claim 1, its feature is in the window by 10~60s of duration, with every The mode that secondary sliding window moves 50% window width chooses the parent belly electric signal that signal acquisition module (1) collects.
  3. 3. signal pre-processing module (2) according to claim 1, it is characterised in that remove industrial frequency noise and baseline drift Using Bessel filter.
  4. 4. signal pre-processing module (2) according to claim 1, it is characterised in that myoelectricity interference is removed to each passage The method used concretely comprises the following steps for phase space partial projection method：
    1) phase space reconfiguration：According to C-C algorithms, calculate and remove industrial frequency noise and baseline drift postabdomen electric signal phase space reconfiguration Optimal embedding dimension m and optimum delay time τ, by one-dimensional belly electric signal be reconstructed into m dimension matrix；
    2) neighborhood is determined：The optimal radius of neighbourhood ε of partial projection is calculated using recursive analysis, determines neighborhood；
    3) non-orthogonal projection and bob-weight structure：In neighborhood, the covariance matrix of the neighborhood square formed by calculating partial vector, estimates Different noise component(s)s is counted, calculates correction value, iteration updates the value of observation station；By bob-weight structure, m dimension matrixes are converted into denoising Belly electrocardiosignal afterwards.
  5. 5. signal quality assessment module (3) according to claim 1, it is characterised in that definition is based on Power Spectrum Distribution Signal quality index (sSQI) is that signal frequency f is the ratio that 5~14Hz frequency bands account for 3~49Hz frequency bands, according to formula:
    <mrow> <mi>S</mi> <mi>D</mi> <mi>R</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>f</mi> <mo>=</mo> <mn>5</mn> </mrow> <mn>14</mn> </munderover> <mi>P</mi> <mi>S</mi> <mi>D</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>,</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>/</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>f</mi> <mo>=</mo> <mn>3</mn> </mrow> <mn>49</mn> </munderover> <mi>P</mi> <mi>S</mi> <mi>D</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>,</mo> <mi>f</mi> <mo>)</mo> </mrow> </mrow>
    <mrow> <mi>s</mi> <mi>S</mi> <mi>Q</mi> <mi>I</mi> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> <mi> </mi> <mi>S</mi> <mi>D</mi> <mi>R</mi> <mo>&amp;GreaterEqual;</mo> <mi>T</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> <mi> </mi> <mi>S</mi> <mi>D</mi> <mi>R</mi> <mo>&lt;</mo> <mi>T</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
    Wherein, PSD (k, f) is power spectral density, and k is current segment signal sequence number, and f is frequency, and threshold value T is adjusted by testing Empirical value, scope are 0.4~0.7, and when the signal sSQI for having a more than passage is labeled as 1, then this period signal is effectively letter Number；If three passage sSQI are marked and are, the segment signal is judged for exception and exports abnormal signal standby signal.
  6. 6. parent electrocardio suppression module (6) according to claim 1, it is characterised in that concretely comprise the following steps：
    1) the average value RR of mother's cardiac RR intervals is calculatedmeanWith R ripple position sequences：Maternal heart rate detection module (5) is transmitted Mother's heart rate time sequence changes into RR interval series, and calculating average value is RRmean, then the R ripple position sequences that add up to obtain；
    2) the QRS wave matrix of each passage is constructed respectively：It is effective after signal quality assessment module (3) processing to each passage The preceding 0.3*RR of each mother's ecg-r wave position of parent belly electric signalmeanWith rear 0.7*RRmeanA length of RR during selectionmeanLetter Number fragment constructs QRS wave matrix X for column vector, wherein " * " is multiplication；
    3) SVD decomposition is carried out to the matrix X of each passage respectively：According to formula X=U Σ VT, U is left singular vector herein, and Σ is removed Cornerwise element is all 0, and the element on diagonal arranges from big to small, is known as singular value, VT(transposition of V) is right unusual Vector；Mother's electrocardio component can be represented by rule of thumb taking preceding 3~5 diagonal element reconstruction signals M of Σ；
    4) respective M signals are individually subtracted in each passage parent belly electric signal and suppress mother's electrocardio component.
CN201810021545.9A 2018-01-10 2018-01-10 System for maternal fetus rhythm of the heart Pending CN108013872A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810021545.9A CN108013872A (en) 2018-01-10 2018-01-10 System for maternal fetus rhythm of the heart

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810021545.9A CN108013872A (en) 2018-01-10 2018-01-10 System for maternal fetus rhythm of the heart

Publications (1)

Publication Number Publication Date
CN108013872A true CN108013872A (en) 2018-05-11

Family

ID=62071410

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810021545.9A Pending CN108013872A (en) 2018-01-10 2018-01-10 System for maternal fetus rhythm of the heart

Country Status (1)

Country Link
CN (1) CN108013872A (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108836315A (en) * 2018-07-23 2018-11-20 智联时空(北京)科技有限公司 Intelligent waistband for Fetal ECG monitoring
CN109009083A (en) * 2018-07-19 2018-12-18 电子科技大学 A kind of Fetal ECG extracting method and device merging small echo and fastICA
CN109199375A (en) * 2018-11-30 2019-01-15 东南大学 A kind of noninvasive Fetal ECG detection device and ecg signal data processing method
CN111265204A (en) * 2019-04-30 2020-06-12 索思（苏州）医疗科技有限公司 Algorithm for extracting fetal heart rate from mother abdomen mixed ECG signal
CN112826513A (en) * 2021-01-05 2021-05-25 华中科技大学 Fetal heart rate detection system based on deep learning and specificity correction on FECG
CN113907765A (en) * 2021-10-10 2022-01-11 北京工业大学 Noninvasive fetal electrocardiosignal quality evaluation method
CN114027852A (en) * 2021-11-11 2022-02-11 浙江智柔科技有限公司 Device and method for analyzing conditions of intrauterine child
CN117442212A (en) * 2023-12-25 2024-01-26 科普云医疗软件(深圳)有限公司 Intelligent monitoring method for obstetrical nursing

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101972145A (en) * 2010-10-12 2011-02-16 华南理工大学 Fetus electrocardio blind separation method based on relative sparsity of time domain of source signal
CN103720469A (en) * 2014-01-02 2014-04-16 山东大学 Wearable type dynamitic maternal fetus electrocardio-detecting device
US20140350421A1 (en) * 2008-11-21 2014-11-27 Massachusetts Institute Of Technology Extraction of fetal cardiac signals
CN105266800A (en) * 2015-12-02 2016-01-27 广东工业大学 Fetal electrocardiogram blind separation method based on low signal-to-noise ratio
CN105411577A (en) * 2015-12-30 2016-03-23 深圳先进技术研究院 Method and system for separating fetal ECG (electrocardiogram)
CN105530857A (en) * 2013-09-09 2016-04-27 皇家飞利浦有限公司 Fetal heart rate extraction from maternal abdominal ECG recordings
US9579055B1 (en) * 2008-10-17 2017-02-28 Orbital Research Inc. Apparatus for non-invasive fetal biosignal acquisition

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9579055B1 (en) * 2008-10-17 2017-02-28 Orbital Research Inc. Apparatus for non-invasive fetal biosignal acquisition
US20140350421A1 (en) * 2008-11-21 2014-11-27 Massachusetts Institute Of Technology Extraction of fetal cardiac signals
CN101972145A (en) * 2010-10-12 2011-02-16 华南理工大学 Fetus electrocardio blind separation method based on relative sparsity of time domain of source signal
CN105530857A (en) * 2013-09-09 2016-04-27 皇家飞利浦有限公司 Fetal heart rate extraction from maternal abdominal ECG recordings
CN103720469A (en) * 2014-01-02 2014-04-16 山东大学 Wearable type dynamitic maternal fetus electrocardio-detecting device
CN105266800A (en) * 2015-12-02 2016-01-27 广东工业大学 Fetal electrocardiogram blind separation method based on low signal-to-noise ratio
CN105411577A (en) * 2015-12-30 2016-03-23 深圳先进技术研究院 Method and system for separating fetal ECG (electrocardiogram)

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
李桥: "重症监护病人心电导联信号质量评估", 《山东大学学报（医学版）》 *
王君: "基于神经网络的混沌时间序列预测", 《中国优秀硕士学位论文全文数据库信息科技辑》 *
王赛红: "基于盲分离的胎心音心率检测算法与实现", 《中国优秀硕士学位论文全文数据库医药卫生科技辑》 *
田文龙: "无创胎儿心电信号提取方法的研究", 《中国优秀硕士学位论文全文数据库医药卫生科技辑》 *
蔡坤: "胎儿心电信号的盲分离研究", 《中国博士学位论文全文数据库信息科技辑》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109009083A (en) * 2018-07-19 2018-12-18 电子科技大学 A kind of Fetal ECG extracting method and device merging small echo and fastICA
CN108836315A (en) * 2018-07-23 2018-11-20 智联时空(北京)科技有限公司 Intelligent waistband for Fetal ECG monitoring
CN109199375A (en) * 2018-11-30 2019-01-15 东南大学 A kind of noninvasive Fetal ECG detection device and ecg signal data processing method
CN109199375B (en) * 2018-11-30 2021-11-02 东南大学 Noninvasive fetal electrocardiogram detection device and electrocardiogram signal data processing method
CN111265204A (en) * 2019-04-30 2020-06-12 索思（苏州）医疗科技有限公司 Algorithm for extracting fetal heart rate from mother abdomen mixed ECG signal
CN112826513A (en) * 2021-01-05 2021-05-25 华中科技大学 Fetal heart rate detection system based on deep learning and specificity correction on FECG
CN113907765A (en) * 2021-10-10 2022-01-11 北京工业大学 Noninvasive fetal electrocardiosignal quality evaluation method
CN113907765B (en) * 2021-10-10 2024-02-23 北京工业大学 Noninvasive fetal electrocardiosignal quality assessment method
CN114027852A (en) * 2021-11-11 2022-02-11 浙江智柔科技有限公司 Device and method for analyzing conditions of intrauterine child
CN117442212A (en) * 2023-12-25 2024-01-26 科普云医疗软件(深圳)有限公司 Intelligent monitoring method for obstetrical nursing
CN117442212B (en) * 2023-12-25 2024-03-12 科普云医疗软件(深圳)有限公司 Intelligent monitoring method for obstetrical nursing

Similar Documents

Publication Publication Date Title
CN108013872A (en) System for maternal fetus rhythm of the heart
CN102160787B (en) Time-frequency-transformation-based blind extraction method of fetal electrocardiography
CN103083013B (en) Electrocardio signal QRS complex wave detection method based on morphology and wavelet transform
CN101972145B (en) Fetus electrocardio blind separation method based on relative sparsity of time domain of source signal
Vullings et al. Dynamic segmentation and linear prediction for maternal ECG removal in antenatal abdominal recordings
CN103263262B (en) System and method for measuring heart rate of fetus
CN106889981B (en) A kind of intelligent terminal for being used to extract fetal heart frequency
CN103610460B (en) A kind of Fetal ECG method for extracting signal based on self adaptation FLANN wave filter
CN102626310A (en) Electrocardiogram signal feature detection algorithm based on wavelet transformation lifting and approximate envelope improving
CN108175384A (en) Method and device based on uterine myoelectricity signal identification uterine contraction
CN103070683A (en) Sleep breathing mode identification method and device based on bioelectrical impedance
CN104887220A (en) Method and system for extracting fetus electrocardiosignals from abdominal wall electrocardiosignals
CN104688220A (en) Method for removing ocular artifacts in EEG signals
CN109199375A (en) A kind of noninvasive Fetal ECG detection device and ecg signal data processing method
CN108577834A (en) A method of it is detected automatically for phase spike between epilepsy
CN107411736A (en) Fetal ECG signal detection system
CN104473631A (en) Fetal electrocardiogram instantaneous heart rate recognition method and system based on non-negative blind separation
CN111265204A (en) Algorithm for extracting fetal heart rate from mother abdomen mixed ECG signal
CN102258368A (en) Time-domain sparsity linear aliasing blind separation model discrimination method in fetal electrocardiogram detection
CN110353704A (en) Mood assessments method and apparatus based on wearable ECG monitoring
CN106889987A (en) Uterine myoelectricity strength information extracting method based on region filtering treatment
CN113907765B (en) Noninvasive fetal electrocardiosignal quality assessment method
Karvounis et al. Detection of fetal heart rate through 3-D phase space analysis from multivariate abdominal recordings
Di Maria et al. An algorithm for the analysis of fetal ECGs from 4-channel non-invasive abdominal recordings
CN105310688A (en) Fetal ECG characteristic signal extraction method based on nonnegative blind separation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20180511