CN104188649B - Ensure that linearly synthesizes a kind of method of real-time in multiple spot physiology pyroelectric monitor - Google Patents
Ensure that linearly synthesizes a kind of method of real-time in multiple spot physiology pyroelectric monitor Download PDFInfo
- Publication number
- CN104188649B CN104188649B CN201410471105.5A CN201410471105A CN104188649B CN 104188649 B CN104188649 B CN 104188649B CN 201410471105 A CN201410471105 A CN 201410471105A CN 104188649 B CN104188649 B CN 104188649B
- Authority
- CN
- China
- Prior art keywords
- monitoring
- linearly
- multiple spot
- signal
- real
- 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.)
- Active
Links
Abstract
The method that linearly synthesizes real-time is ensured in a kind of multiple spot physiology pyroelectric monitor, including step, (1) multiple measuring points of body surface are arranged as needed, connect monitoring system, prepare monitoring, (2) premonitoring records one section of multiple spot electro-physiological signals, with periodic component analytic approach or Independent component analysis or other method, solve corresponding optimization problem, obtain corresponding signal optimum linearity mix vector, (3) formal monitoring is entered, the signal optimum linearity mix vector determined with (2nd) step, multiple spot electro-physiological signals to each section of new monitoring are linearly synthesized, to synthesize the purified signal of needs, it is characterized in that：The signal optimum linearity mix vector that (2nd) step premonitoring stage was determined, it has been directly used in the linearly synthesis in (3rd) the formal monitoring process of step, and (3rd) step formally monitors stage no longer solving-optimizing problem, so as to ensure that the whole formal monitoring stage linearly synthesizes the real-time for the treatment of.
Description
Technical field
The application is related to ensure the method that linearly synthesizes real-time in a kind of multiple spot physiology pyroelectric monitor.
Dynamic monitoring electrocardiogram, dynamic monitoring electroencephalogram, dynamic monitoring electromyogram etc., the prison of these Human Physiology electric signals
Application is surveyed, is monitored in multiple points of body surface.The electro-physiological signals monitored by body surface each point are substantially all each
Plant the mixed signal of electro-physiological signals；For example, be often mixed with the EEG signals of head detection eye electricity, electrocardio and other interference and
Noise, from her abdominal detection fetal rhythm electric signal in be mixed with very big parent electrocardio signal, also have respiratory wave, myoelectricity and other
Interference and noise etc..People carry out linear combination by the multichannel electro-physiological signals in body surface multiple spot monitoring, can synthesize
More clean simple electro-physiological signals, greatly improve signal to noise ratio.
Multiple spot electro-physiological signals are carried out with linear synthesis and is related to the determination of optimum linearity mix vector.It is so-called to determine optimal line
Property mix vector, as the term suggests, seek to solve the optimization problem under a given object function, and often non-convex is more
Extremal optimization problem, amount of calculation is very big.In multiple spot physiology pyroelectric monitor, if each the time segment signal to monitoring first is asked
One optimization problem of solution, determines corresponding optimum linearity mix vector, and then the time segment signal is linearly closed accordingly again
Into treatment, then real-time is difficult to ensure that.At present, the various linear composition algorithms of multiple spot electro-physiological signals of research, are all using every
Individual time segment signal determines the optimum linearity mix vector of oneself, is suitable to offline Non real-time processing；In actual multiple spot physiology electric
In monitoring system, there is not yet linearly composition algorithm is applied to real-time processing.
The application seeks to propose a solution, so that linearly synthesis treatment can be applied to multiple spot physiology electric prison
In examining system, real-time processing is carried out.
Background technology
Multiple spot electro-physiological signals are linearly synthesized, there are periodic component analytic approach (Periodic Component
Analysis), Independent component analysis (independent component analysis, ICA) etc..These methods, it is determined that
Optimization object function during optimum linearity mix vector is not quite similar, but as it was previously stated, be all to solve for optimization problem and it is final really
Surely for the optimum linearity mix vector of specific synthesis application.Such as, have from M M roads mixed signal of point monitoring of pregnant woman's stomach wall
xiN () n=1 ..., N, i=1 ..., M, N are the sampling number for monitoring the time period, now synthesize one with this M roads linearly
The clean parent electrocardio signal mECG in road and Fetal ECG signal fECG clean all the way, no matter then how, finally all
An optimum linearity mix vector W for synthesis maternal ecg is determinedmWith an optimum linearity group for synthesis FECG
Resultant vector Wf.If note Wm=[wm1, wm2..., wmM]Wf=[wf1, wf2..., wfM], then have：
Preceding to have addressed, the algorithm that existing multiple spot electro-physiological signals linearly synthesize is determined using each time segment signal
The optimum linearity mix vector of oneself, is adapted for offline Non real-time processing；In actual multiple spot physiology electric real-time monitoring system
In, due to computationally intensive, there is not yet these algorithms are applied to real-time processing.The real-time of system is improved, can be with
By the processing speed for improving computer hardware, it is also possible to which specific optimized algorithm is improved, but at present apparently, can be pre-
The future seen, these approach all cannot solve problem at all.Must look for another way and propose a plan, to ensure multiple spot electro-physiological signals
The real-time of linear synthesis, can be applied in physiology electric real-time monitoring system.
Bibliography：
[1] Mohammed Assam Ouali, and Kheireddine Chafaa, " Separation of
Composite maternal ECG using SVD decomposition, " Computer Applications
Technology (ICCAT), 2013International Conference on, pp.1-4, Jan.2013.
[2] A.Van Oosterom, " Spatial filtering of the fetal electrocardiogram, "
J.Perinat.Med, vol.14, pp.411-419,1986.
[3] Zhenwei Shi, and Changshui Zhang, " Semi-blind source extraction for
fetal electrocardiogram extraction by combining non-Gaussianity and time-
Correlation, " Neurocomputing, vol.70, pp.1547-1581,2007.
[4] Reze Sameni, Christian Jutten, and Mohammad B.Shamsollahi,
“Multichannel Electrocardiogram Decomposition Using Periodic Component
Analysis, " IEEE Bio.Med.Eng., vol.55, No.8, pp.1935-1940, Aug.2008.
[5] Thato Tsalaile, Reza Sameni, Saeid Sanei, et al., " Sequential blind
Source extraction for quasi-periodic signals with time-varying period, " IEEE
Bio.Med.Eng., vol.56, No.3, pp.646-655, Mar.2009.
[6] Aapo Hyvarinen and Erkki Oja, " A fast fixed-point algorithm for
Independent component analysis, " Neural Comput.vol.9, pp.1483-1492,1997.
[7] Jean-Francois Cardoso, " Source Separation Using Higher Order
Moments, " ICASSP, proceedings, vol.4, May 1989.
[8] Wei Lu, and Jagath C.Rajapakse, " Approach and Applications of
Constrained ICA, " IEEE TNeuralNetwor., vol.16, No.1, pp.203-212, Jan.2005.
[9] Wei Lu, and Jagath C.Rajapakse, " ICA with Reference, " Neurocomputing,
Vol.69, pp.2244-2257,2006.
The content of the invention
Goal of the invention.
The existing linear composition algorithm of multiple spot electro-physiological signals, is all based on each monitoring time segment signal and determines this in itself
The optimum linearity mix vector of segment signal, is suitable to carry out offline Non real-time processing to the multiple spot electro-physiological signals monitored.This hair
Bright purpose is exactly to look for another way to propose a kind of scheme that can ensure the linear real-time of synthesis treatment of multiple spot electro-physiological signals, with
Meet demand of the actual physiology electric monitoring system to linearly synthesis treatment real-time.
Technical scheme.
Ensure the method that linearly synthesizes real-time in a kind of multiple spot physiology pyroelectric monitor, including step, (1) as needed
Multiple measuring points of body surface are arranged, monitoring system are connected, is prepared monitoring, (2) premonitoring records one section of multiple spot electro-physiological signals,
With periodic component analytic approach or Independent component analysis, corresponding optimization problem is solved, obtain corresponding signal optimum linearity
Mix vector, (3) newly monitor into formal monitoring, the signal optimum linearity mix vector determined with (2nd) step to each section
Multiple spot electro-physiological signals are linearly synthesized, to synthesize the purified signal of needs, it is characterized in that：(2nd) step premonitoring rank
The signal optimum linearity mix vector that section determines, has been directly used in the linearly synthesis in (3rd) the formal monitoring process of step, and
(3rd) step formally monitors stage no longer solving-optimizing problem, so as to ensure that the whole formal monitoring stage linearly synthesizes treatment
Real-time.The FB(flow block) of method is as shown in Figure 1.
The scheme for ensureing linearly synthesis real-time above is not arbitrarily conjecture, but grinding based on some basic laws
Study carefully and recognize what is proposed.These rules include：(1) the linear side of synthesis of the multi-point signal such as periodic component analysis, independent component analysis
Method, can be used for the blind source separating of multiple spot mixed signal, and each relatively pure source signal is isolated by linear synthesis, these
Method is confirmed from theory into action；(2) in theory, the linear synthesis side such as periodic component analysis, independent component analysis
Optimum linearity mix vector in method, is substantially to produce point to come to the characteristic of each monitoring point bang path of body surface by each signal source
Determine, be exactly by the physiology electric such as heart, brain, eyes generating source to each monitoring point of body surface for multiple spot physiology pyroelectric monitor
Bang path characteristic determine；(3) in multiple spot physiology electric monitoring process, the physiology electric generating source such as heart, brain, eyes
To each monitoring point of body surface position is usually changeless, bang path is metastable.Therefore, it is the premonitoring stage is true
The linear synthesis that fixed optimum linearity mix vector is directly used in formal monitoring process below is feasible from reason.Hair
A person of good sense's this conclusion also by lot of experiment validation, correlative theses will be delivered.
Beneficial effect.
First, application scheme obviously determines this section of optimum linearity mix vector with existing using the currently monitored signal segment
Method it is different, therefore, with novelty.Second, the optimal composite vector that the premonitoring stage determines is directly used in below
Formal monitoring process, makes formal monitoring process not further relate to solving-optimizing problem, and only merely multi-point signal is carried out linearly to add
Power summation, amount of calculation is greatly reduced, and can utterly ensure the real-time of linearly synthesis, therefore, this is a substantive progress.
Accompanying drawing 2 is online disclosed pregnant woman's stomach wall pyroelectric monitor data, and Fig. 2 (a) (b) (c) corresponds to 3 stomach wall measuring points respectively
Electro-physiological signals data (0-10 seconds section), they can all regard the mixed signal of parent electrocardio, fetal electrocardiogram and noise etc. as.
For the sake of simplicity, here as a example by synthesizing clean parent electrocardio source signal：Inventor is using 0-10 seconds section multi-point signal as pre-
Monitoring phase data, by independent component analysis method determine the optimum linearity mix vector of the segment data for Wm=[-
0.4748,0.0523, -0.8786], this section of maternal ecg result such as accompanying drawing 3 (a) with this vector synthesis is shown, it is clear that ratio
Former stomach wall electricity mixed signal is more pure；Then, by the premonitoring stage determine optimum linearity mix vector Wm=[- 0.4748,
0.0523, -0.8786] the linear synthesis of later each segment data is directly used in, equally obtain successfully, accompanying drawing 3 (b)-(h) is respectively
Show 20-30 seconds section, 70-80 seconds section, 120-130 seconds section, 170-180 seconds section, 220-230 seconds section, 270-280 seconds section with
And the 300-310 seconds composite result of section, it is obtained for the parent electrocardio for comparing " clean ".This example illustrates having for the present invention program
Effect property.
Brief description of the drawings
Fig. 1, ensures that linearly synthesizes the block diagram of Real-Time Performance in multiple spot physiology pyroelectric monitor of the present invention.
Fig. 2,3 physiology electric Monitoring Datas of premonitoring stage body surface in embodiment
Fig. 3, the whole monitoring some linear composite result schematic diagrames of stage signal in embodiment
Specific embodiment
Using online disclosed database MIT Non-Invasive Fetal Electrocardiogram Database
Data carry out example as multiple spot monitoring data.Include 55 groups of data in the database altogether, sample rate is 1kHz, every group of data
Include that 2 road chests lead signal and 3 or 4 road stomach wall electric signals.For convenience, the number that suite number is " ecgca748 " is only taken out
According to 3 stomach wall electric signals be specifically described, see shown in accompanying drawing (2)；They are all mixed signals, include parent electrocardio,
The noise such as fetal electrocardiogram and myoelectricity, in addition to parent electrocardio, other compositions are weaker；In the present embodiment, intend by the linear of them
Synthesis is obtained as a example by the parent electrocardio of " clean " all the way.
This group of total monitoring length of data is 318 seconds, preceding to have addressed, the step of by this programme, with its data of first 10 seconds
Used as premonitoring phase data, follow-up data (was marked altogether using every 10 seconds as formal monitoring phase data as a signal segment
30 sections).In the premonitoring stage, using independent component analysis method determine the optimum linearity for obtaining the segment data parent electrocardio combine to
Measure is Wm=[- 0.4748,0.0523, -0.8786]；Then, the optimum linearity mix vector Wm=for the premonitoring stage being determined
[- 0.4748,0.0523, -0.8786] is directly used in the linear synthesis of later data.Accompanying drawing 3 (b)-(h) respectively illustrates
20-30 seconds section, 70-80 seconds section, 120-130 seconds section, 170-180 seconds section, 220-230 seconds section, 270-280 seconds section and 300-
310 seconds composite results of section, are all the parent electrocardios for comparing " clean ", and synthesis is achieved successfully, illustrates the feasibility of method.
Claims (1)
1. the method that linearly synthesizes real-time ensured in a kind of multiple spot physiology pyroelectric monitor, including step, (1) cloth as needed
Multiple measuring points of body surface are put, monitoring system has been connected, is prepared monitoring, (2) premonitoring has recorded one section of multiple spot electro-physiological signals, uses
Periodic component analytic approach or Independent component analysis, solve corresponding optimization problem, obtain corresponding signal optimum linearity group
Resultant vector, (3), into formal monitoring, the signal optimum linearity mix vector determined with (2nd) step, what each section was newly monitored is more
Point electro-physiological signals are linearly synthesized, to synthesize the purified signal of needs, it is characterized in that：In (2nd) step premonitoring stage
The signal optimum linearity mix vector of determination, the linearly being directly used in (3rd) the formal monitoring process of step synthesizes, and the
(3) step formally monitors stage no longer solving-optimizing problem, so as to ensure that the whole formal monitoring stage linearly synthesizes the reality for the treatment of
Shi Xing.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410471105.5A CN104188649B (en) | 2014-09-15 | 2014-09-15 | Ensure that linearly synthesizes a kind of method of real-time in multiple spot physiology pyroelectric monitor |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410471105.5A CN104188649B (en) | 2014-09-15 | 2014-09-15 | Ensure that linearly synthesizes a kind of method of real-time in multiple spot physiology pyroelectric monitor |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104188649A CN104188649A (en) | 2014-12-10 |
| CN104188649B true CN104188649B (en) | 2017-06-06 |
Family
ID=52074157
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410471105.5A Active CN104188649B (en) | 2014-09-15 | 2014-09-15 | Ensure that linearly synthesizes a kind of method of real-time in multiple spot physiology pyroelectric monitor |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104188649B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109770889B (en) * | 2017-11-15 | 2022-03-11 | 深圳市理邦精密仪器股份有限公司 | Electrocardiogram data section selection method and device |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2003003905A2 (en) * | 2001-07-05 | 2003-01-16 | Softmax, Inc. | System and method for separating cardiac signals |
| WO2005044101A1 (en) * | 2003-10-31 | 2005-05-19 | The Board Of Trustees Of The University Of Illinois | Separation of one or more fetal heart component signals from heart signal information obtained from a pregnant female |
| CN101596108A (en) * | 2009-06-19 | 2009-12-09 | 南京大学 | The nonlinear separation and extract methods of fetal electrocardiogram |
| CN104027105A (en) * | 2014-04-23 | 2014-09-10 | 河南科技大学 | Novel maternal and fetal electrocardiogram separation method |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7333850B2 (en) * | 2004-05-28 | 2008-02-19 | University Of Florida Research Foundation, Inc. | Maternal-fetal monitoring system |
| CN1951320A (en) * | 2006-11-16 | 2007-04-25 | 上海交通大学 | Abnormal electrocardiogram recognition method based on ultra-complete characteristics |
| CN201005696Y (en) * | 2007-07-17 | 2008-01-16 | 中国人民解放军南京军区南京总医院 | Brain inducing electrical signal tester based on separate component analysis technique |
| CN101488189B (en) * | 2009-02-04 | 2012-01-18 | 天津大学 | Brain-electrical signal processing method based on isolated component automatic clustering process |
| CN102008300A (en) * | 2010-12-10 | 2011-04-13 | 吉林大学 | Wearable multiple physiological parameter recording device |
| CN103431859B (en) * | 2012-01-10 | 2014-12-10 | 西安交通大学 | Experimental method for determining brain load in multitask visual cognition |
| CN102835955B (en) * | 2012-09-08 | 2014-02-26 | 北京工业大学 | Method of automatically removing ocular artifacts from electroencephalogram signal without setting threshold value |
| CN104688220B (en) * | 2015-01-28 | 2017-04-19 | 西安交通大学 | Method for removing ocular artifacts in EEG signals |
-
2014
- 2014-09-15 CN CN201410471105.5A patent/CN104188649B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2003003905A2 (en) * | 2001-07-05 | 2003-01-16 | Softmax, Inc. | System and method for separating cardiac signals |
| WO2005044101A1 (en) * | 2003-10-31 | 2005-05-19 | The Board Of Trustees Of The University Of Illinois | Separation of one or more fetal heart component signals from heart signal information obtained from a pregnant female |
| CN101596108A (en) * | 2009-06-19 | 2009-12-09 | 南京大学 | The nonlinear separation and extract methods of fetal electrocardiogram |
| CN104027105A (en) * | 2014-04-23 | 2014-09-10 | 河南科技大学 | Novel maternal and fetal electrocardiogram separation method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104188649A (en) | 2014-12-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Dammers et al. | Integration of amplitude and phase statistics for complete artifact removal in independent components of neuromagnetic recordings | |
| Zong et al. | Hilbert-Huang transform based physiological signals analysis for emotion recognition | |
| Gonuguntla et al. | Event-related functional network identification: application to EEG classification | |
| Kahankova et al. | Non-adaptive methods of fetal ECG signal processing | |
| Veluvolu et al. | Adaptive estimation of EEG-rhythms for optimal band identification in BCI | |
| JP6423657B2 (en) | EEG signal analysis result display device and EEG signal analysis result display method | |
| Wang et al. | Feature extraction of brain-computer interface based on improved multivariate adaptive autoregressive models | |
| Bai et al. | The offline feature extraction of four-class motor imagery EEG based on ICA and Wavelet-CSP | |
| Gonuguntla et al. | Phase synchrony in subject-specific reactive band of EEG for classification of motor imagery tasks | |
| CN106236080A (en) | Based on the removing method of myoelectricity noise in multichannel EEG signals | |
| CN103932701A (en) | Individual brain function mapping method based on electrocorticogram high-frequency Gamma nerve oscillation | |
| Zou et al. | Automatic EEG artifact removal based on ICA and Hierarchical Clustering | |
| Khalighinejad et al. | NAPLib: An open source toolbox for real-time and offline Neural Acoustic Processing | |
| Liu et al. | An efficient and robust muscle artifact removal method for few-channel EEG | |
| Miao et al. | Automated CCA-MWF algorithm for unsupervised identification and removal of EOG artifacts from EEG | |
| Chang et al. | Inter-trial analysis of post-movement beta activities in EEG signals using multivariate empirical mode decomposition | |
| CN104188649B (en) | Ensure that linearly synthesizes a kind of method of real-time in multiple spot physiology pyroelectric monitor | |
| Hoole et al. | Autism, EEG and brain electromagnetics research | |
| Haritopoulos et al. | Foetal PQRST extraction from ECG recordings using cyclostationarity-based source separation method | |
| Ali et al. | Time-frequency coherence of multichannel EEG signals: Synchrosqueezing transform based analysis | |
| Haritopoulos et al. | Cyclostationarity-based estimation of the foetus subspace dimension from ECG recordings | |
| CN107450730B (en) | Low-speed eye movement identification method and system based on convolution mixed model | |
| Zhou et al. | Fetal Electrocardiogram Extraction and Performance Analysis. | |
| Gao et al. | Multivariate empirical mode decomposition and multiscale entropy analysis of EEG signals from SSVEP-based BCI system | |
| Ahmed et al. | Emotion recognition based on correlation between left and right frontal EEG assymetry |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |