CN108670238A - The autonomous regulating power evaluation method of brain blood flow based on brain oxygen signal and blood pressure signal - Google Patents
The autonomous regulating power evaluation method of brain blood flow based on brain oxygen signal and blood pressure signal Download PDFInfo
- Publication number
- CN108670238A CN108670238A CN201810481779.1A CN201810481779A CN108670238A CN 108670238 A CN108670238 A CN 108670238A CN 201810481779 A CN201810481779 A CN 201810481779A CN 108670238 A CN108670238 A CN 108670238A
- Authority
- CN
- China
- Prior art keywords
- brain
- correlation
- local
- hemoglobin concentration
- small echo
- 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
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/0205—Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
-
- 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/021—Measuring pressure in heart or blood vessels
-
- 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/026—Measuring blood flow
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/14542—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring blood gases
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/14546—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring analytes not otherwise provided for, e.g. ions, cytochromes
-
- 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/7235—Details of waveform analysis
-
- 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/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/726—Details of waveform analysis characterised by using transforms using Wavelet transforms
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Optics & Photonics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Vascular Medicine (AREA)
- Pulmonology (AREA)
- Hematology (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
The autonomous regulating power evaluation method of the brain blood flow based on brain oxygen signal and blood pressure signal of the present invention detects body local brain tissue pathological changes, local organization hemoglobin concentration index and local brain tissue total hemoglobin concentration variable quantity first；And measure human body artery pressure；Then correlation and phase correlation between each parameter are analyzed using small echo crossbar transistion method；In the small echo correlation of different frequency range and Wavelet Phase correlation between ultimate analysis brain oxygen parameter and angiosthenia.The present invention is by measuring human brain cortex hemoglobin concentration, brain oxygen saturation and angiosthenia, and correlation between several measurement results, the autonomous regulating power of brain blood flow can be evaluated indirectly, it is preoperative to estimate patient death rate, the range and postoperative evaluation patient's recovery situation of the angiosthenia of different people are quantitatively determined in art.
Description
Technical field
The present invention relates to a kind of biomedical engineering technology fields, more particularly to a kind of to be based on brain oxygen signal and blood pressure signal
The autonomous regulating power evaluation method of brain blood flow.
Background technology
By measuring human brain cortex hemoglobin concentration, brain oxygen saturation and angiosthenia, and pass through ad hoc fashion meter
It counts stating correlation between several measurement results in, brain blood flow can be evaluated indirectly from main regulation (CA, Cerebral
Autoregulation) ability estimates patient death rate so as to preoperative, the angiosthenia of quantitative determination different people in art
Range and postoperative evaluation patient's recovery situation.
Based on brain blood flow from the assessment of main regulation, has relevant personnel and studied, the Mx that such as clinically uses,
PRx.But their use is there is significant limitation, what Mx parameters indicated be blood flow velocity in middle cerebral and angiosthenia it
Between correlation, and blood flow velocity in middle cerebral is obtained by transcranial Doppler technology measurement, and this measurement method can
The technology of operator is relied heavily on by property, and when bone window penetrates bad can not be checked, and its can only carry out it is short
The continuous measurement of time.What PRx was indicated is the relationship between intracranial pressure and angiosthenia, and intracranial pressure is needed when measuring by encephalic
Pressure detection survey meter probe is placed in encephalic, that is, it is invasive to measure, therefore significantly limits its scope of application.Brain oxygen is indirect
Reflect brain blood flow situation of change, calculate the correlation between brain oxygen and angiosthenia, may be implemented clinically it is noninvasive, continuous,
Monitoring in real time, and it is easy to operate.The existing related researcher of relationship between brain oxygen and angiosthenia is made at present certain
Research, but algorithm is incomplete, and accurate fixed and validity needs further to be verified, and does not promote also clinically.
Wavelet analysis is the frontier of a rapid development in current application mathematics and engineering discipline, by nearly spy in 30 years
Rope is studied, and important mathematical form system has built up.Compared with Fourier transformation, wavelet transformation is space (time) and frequency
The partial transformation of rate, thus information can be effectively extracted from signal.By the calculation functions such as flexible and translation can to function or
Signal carries out multiple dimensioned refinement analysis, solves the indeterminable many difficult problems of Fourier transformation.
The present invention proposes new research method and parameter with regard to the relationship between brain oxygen and angiosthenia, uses wavelet analysis side
Method analyzes parameter with regard to different frequency range on frequency domain, and different frequency ranges has different physiological significances, therefore makes this hair
Bright evaluation method is more targeted.
Invention content
The present invention provides a kind of autonomous regulating power evaluation method of the brain blood flow based on brain oxygen signal and blood pressure signal, solves
Problems of the prior art improve the precision and detection efficiency of the autonomous regulating power evaluation of brain blood flow.
In a first aspect, technical solution of the present invention is as follows：
A kind of autonomous regulating power evaluation method of brain blood flow based on brain oxygen signal and blood pressure signal of the present invention, including with
Lower step：
Step (1), detection body local brain tissue pathological changes TOI, local organization hemoglobin concentration index THI drawn games
Portion brain tissue total hemoglobin concentration variation delta CtHb；
Step (2), while carrying out step (1), measure human body artery pressure；
Step (3), using small echo crossbar transistion method, analyze the body local brain tissue pathological changes TOI, local organization
Phase between hemoglobin concentration index THI, local brain tissue total hemoglobin concentration variation delta CtHb and human body artery pressure
Closing property and phase correlation；
Small echo correlation between step (4), analysis brain oxygen parameter and angiosthenia in different frequency range is related to Wavelet Phase
Property.
Preferably, the step (3), including：
Continuous wavelet transform process data, on a specific time, each frequency is as a corresponding window
Mouthful, on the window, the Pan and Zoom of quadractically integrable function Ψ 0 (η) that mean value is 0 in time can to frequency into
Row partial analysis, Ψ0(η) is known as morther wavelet；
Continuous wavelet transform is expressed as：
S indicates that scale, δ t indicate sampling interval, xn(n=1,2,3 ... N) indicates time series；
The intersection small echo of two different time domain signals Xn and Yn are defined as：
WXY=WXWY*Formula (2)
* indicates complex conjugate in formula, definition | WXY| it is power spectrum, the intersection small echo W of two different time domain signals Xn and YnXY
It is the result of two signal Xn and Yn crossbar transistion；Its absolute value is power spectrum, wherein complex matrix | WXY| phase angle be two
Phase differences of a time series Xn and Yn on time-frequency domain:
Φ in formulai(s) be two time series in time tiPhase delay under scale s；
Step (4), including：
Small echo correlation (WCO) formula is：
Wavelet Phase correlation (WPCO) formula is：
Wherein,It is two groups of signal phase differences being averaged in the time domain.
Preferably, in step (1), using Near-infrared Brain oxygen detection detection body local brain tissue pathological changes TOI,
Local organization hemoglobin concentration index THI and local brain tissue total hemoglobin concentration variation delta CtHb；Specially：
Two electrode slices are attached to the forehead of human body, are measured 30-60 minutes, body local brain tissue pathological changes are extracted
TOI and part tissue hemoglobin concentration index THI and local brain tissue total hemoglobin concentration variation delta CtHb.
Second aspect, technical solution of the present invention are as follows：
A kind of autonomous regulating power evaluating apparatus of brain blood flow based on brain oxygen signal and blood pressure signal, including：
First detection unit, for detecting body local brain tissue pathological changes TOI, local organization hemoglobin concentration refers to
Number THI and local brain tissue total hemoglobin concentration variation delta CtHb；
Second detection unit, for while first detection unit is run, measuring human body artery pressure；
Analytic unit analyzes the body local brain tissue pathological changes TOI, office for using small echo crossbar transistion method
Portion tissue hemoglobin concentration index THI, local brain tissue total hemoglobin concentration variation delta CtHb and human body artery press it
Between correlation and phase correlation, and in the small echo correlation of different frequency range and small between analysis brain oxygen parameter and angiosthenia
Wave phase correlation.
Preferably, the analytic unit is used for：
Continuous wavelet transform process data, on a specific time, each frequency is as a corresponding window
Mouthful, on the window, quadractically integrable function Ψ that mean value is 00The Pan and Zoom of (η) in time can to frequency into
Row partial analysis, Ψ0(η) is known as morther wavelet；
Continuous wavelet transform is expressed as：
S indicates that scale, δ t indicate sampling interval, xn(n=1,2,3 ... N) indicates time series；
The intersection small echo of two different time domain signals Xn and Yn are defined as：
WXY=WXWY*Formula (2)
* indicates complex conjugate in formula, definition | WXY| it is power spectrum, the intersection small echo W of two different time domain signals Xn and YnXY
It is the result of two signal Xn and Yn crossbar transistion；Its absolute value is power spectrum, wherein complex matrix | WXY| phase angle be two
Phase differences of a time series Xn and Yn on time-frequency domain:
Φ in formulai(s) be two time series in time tiPhase delay under scale s；
Analyze small echo correlation and Wavelet Phase correlation in different frequency range between brain oxygen parameter and angiosthenia：
Small echo correlation (WCO) formula is：
Wavelet Phase correlation (WPCO) formula is：
Wherein,It is two groups of signal phase differences being averaged in the time domain.
Technical solution provided by the invention has the advantages that：
The present invention proposes new research method and parameter with regard to the relationship between brain oxygen and angiosthenia, uses wavelet analysis side
Method analyzes parameter with regard to different frequency range on frequency domain, and different frequency ranges has different physiological significances, therefore measures more
With specific aim.The present invention realizes noninvasive real-time monitoring of the brain blood flow from main regulation, to for it is preoperative, in art and postoperative disease
The monitoring of people's state provides reliable reference.
Compared with presently disclosed the relevant technologies, the features of the present invention and advantage are：First, the method for the present invention can be with
Reach noninvasive, real-time monitoring and evaluation, it is easy to operate；Second, algorithm of the present invention is novel, accuracy is high, and has
Specific aim.
Description of the drawings
Fig. 1 is that human body spontaneity blood vessel shakes schematic diagram；Wherein (a) is human vas schematic diagram；(b) it is indicated with (d)
It is the blood pressure concussion schematic diagram of big blood vessel under time domain and frequency domain；(c) what is indicated with (d) is the blood of capillary under time domain and frequency domain
Pressure concussion schematic diagram；
Fig. 2 is that the physiological significance of spontaneous concussion under human vas different frequency range compares schematic diagram；
Fig. 3 is the autonomous regulating power evaluation of a kind of brain blood flow based on brain oxygen signal and blood pressure signal in the embodiment of the present application
The flow chart of method；
Fig. 4 is the autonomous regulating power evaluation of a kind of brain blood flow based on brain oxygen signal and blood pressure signal in the embodiment of the present application
The structure chart of device.
Specific implementation mode
In order to understand the present invention in depth, the present invention is described in detail below in conjunction with the accompanying drawings.
The present invention obtains brain oxygen parameter by Near-infrared Brain oxygen monitor first and measures human body artery pressure, then by small
Wave analysis method analyzes the relationship between brain oxygen parameter and angiosthenia.Although having related researcher with regard between brain oxygen and blood pressure
Relationship studied, but their calculation is unfolded in the time domain.And it is needed sometimes in the fluctuation of angiosthenia
It carries out artificial intervention or prolonged measure has obtained the blood pressure fluctuated by a relatively large margin, and the accuracy measured also needs
Further verification.The Algorithm Analysis that the present invention uses is in blood pressure and brain oxygen with regard to being analyzed on different frequency range
It carries out under autonomous state, do not need human intervention or measures for a long time, therefore detect simplicity, unfavorable shadow will not be caused to human body
It rings, and data are more targeted.
As shown in Figure 1, under autonomous state, there are the concussions of varying strength in different frequency range pressures for human body artery blood pressure；
And the concussion of different frequency range has different physiological significances, as shown in Figure 2.
Blood flow (oxygen) signal of human body is non-stationary signal, has apparent time and frequency domain characteristics.If with traditional Fourier
Leaf transformation analysis method, which carries out analysis, has the shortcomings that resolution ratio is low, spectrum is rough.And wavelet analysis can be non-flat for analyzing
Steady signal, have signal high band temporal resolution precision it is higher, in the higher spy of signal low-frequency range frequency resolution precision
Point.But wavelet transformation can only discuss single time series signal, it is difficult to analyze influencing each other between more element sequence signals and when
Frequency correlation.Cross wavelet analysis can diagnose correlation between unlike signal, when ductility and phase structure.So this is of the invention
Analysis on frequency domain is carried out using cross wavelet analysis analysis method.
Embodiment one
As shown in figure 3, the autonomous regulating power evaluation method of the brain blood flow based on brain oxygen signal and blood pressure signal of the present invention,
Include the following steps：
Step (1)：Detect body local brain tissue pathological changes TOI, local organization hemoglobin concentration index THI drawn games
Portion brain tissue total hemoglobin concentration variation delta CtHb.
Two electrode slices are attached to the forehead of human body, are measured 30-60 minutes, Regional cerebral oxygen saturation TOI, office are measured
Portion tissue hemoglobin concentration index THI and local brain tissue total hemoglobin concentration variation delta CtHb.The brain group measured
It is the time series changed over time to knit oxygen saturation TOI, time series Xn in continuous wavelet transform in corresponding step (3)；It surveys
The local brain tissue total hemoglobin concentration variation delta CtHb measured is also the time series changed over time, also corresponding continuous
Time series Xn in wavelet transformation.
Step (2)：While measuring brain oxygen relevant parameter, human body artery pressure is measured.The arterial pressure signal measured is
The time series changed over time corresponds to time series Xn in continuous wavelet transform.
By step (1) and step (2), brain tissue pathological changes TOI, the local organization hemoglobin for extracting measurement are dense
It spends index THI, local brain tissue total hemoglobin concentration variation delta CtHb and human body artery and presses initial data.Extract completion
Afterwards, it is also necessary to it is filtered, the relevant treatments such as noise reduction, then carry out the small echo cross-correlation analysis of step (3).
Step (3)：Using small echo crossbar transistion method, the body local brain tissue pathological changes TOI, local organization are analyzed
Phase between hemoglobin concentration index THI, local brain tissue total hemoglobin concentration variation delta CtHb and human body artery pressure
Closing property and phase correlation.
Continuous wavelet transform process data, on a specific time, each frequency can regard as one it is corresponding
Window.On the window, quadractically integrable function Ψ that mean value is 00The Pan and Zoom of (η) in time can be to frequency
Carry out partial analysis, Ψ0(η) is known as morther wavelet.
Continuous wavelet transform is expressed as：
S indicates that scale, δ t indicate sampling interval, xn(n=1,2,3 ... N) indicates time series.
The intersection small echo W of two different time domain signals Xn and YnXYIt is defined as：
WXY=WXWY*Formula (2)
* indicates complex conjugate in formula；|WXY| it is power spectrum；The intersection small echo W of two different time domain signals Xn and YnXYIt is two
The crossbar transistion result of a signal；Its absolute value is power spectrum；Wherein power spectrum complex matrix | WXY| phase angle can be seen
Make the phase difference of two time serieses Xn and Yn on time-frequency domain:
Φ in formulai(s) be two time series in time tiPhase delay under scale s.
Step (4)：Analyze brain oxygen parameter (i.e.：TOI, THI and Δ CtHb) between angiosthenia in the small echo of different frequency range
Correlation and Wavelet Phase correlation.
Under normal conditions, at low-frequency range (about 0.005-0.145Hz), small echo correlation is lower, and brain blood flow is from main regulation
Function is relatively preferable.
Small echo correlation (WCO) formula is：
Wavelet Phase correlation (WPCO) formula is：
WhereinIt is two groups of signal phase differences being averaged in the time domain.
Wavelet Phase correlation has reacted the synchronism of phase between angiosthenia and brain oxygen parameter, and general Wavelet Phase is related
Property illustrates that synchronism is good between angiosthenia and brain oxygen parameter when higher, and CA regulating powers are poor；When relatively low, synchronism is poor,
CA regulating powers are relatively preferable.Related researcher has found the Pearson came between brain oxygen and angiosthenia when carrying out time-domain analysis
When correlation is less than 0.3-0.5, brain blood flow is relatively normal from main regulation, judges whether CA is relatively normal related based on wavelet analysis
The critical value of property needs further experiment to determine.
Figure of description human body spontaneous blood vessel shakes in schematic diagram 1, and (a) is human vas schematic diagram；(b) and (d) table
What is shown is the blood pressure concussion of big blood vessel under time domain and frequency domain；(c) what is indicated with (d) is the blood pressure of capillary under time domain and frequency domain
Concussion.
Since infant's brain development is perfect not enough, correlation between brain oxygen parameter and angiosthenia may relative at
Nian Renhui is relatively higher.And at high band (about 0.145-2Hz), fluctuation of blood pressure feelings have been reacted in the fluctuation of brain oxygen to a certain extent
Condition, therefore the correlation between brain oxygen and angiosthenia is higher, brain function is relatively preferable.Under normal circumstances, in high band normal person
The opposite CA impaired subjects of correlation meeting between brain oxygen parameter and angiosthenia are higher.Figure of description 2 is human vas different frequencies
The physiological significance spontaneously shaken under section.
The autonomous regulating power evaluation method of the brain blood flow based on brain oxygen signal and blood pressure signal of the present invention is with regard to brain oxygen and moves
Relationship between pulse pressure proposes new research method and parameter, using wavelet analysis method on frequency domain with regard to different frequency range to ginseng
Number is analyzed, and different frequency ranges has different physiological significances, therefore is measured more targeted.Meanwhile the present invention is real
Showed noninvasive real-time monitoring of the brain blood flow from main regulation, to for it is preoperative, in art and the monitoring of postoperative patient condition provide can
The reference leaned on.
Embodiment two
Based on same inventive concept, as shown in figure 4, present embodiments providing a kind of based on brain oxygen signal and blood pressure signal
The autonomous regulating power evaluating apparatus of brain blood flow, including：
First detection unit 401, it is dense for detecting body local brain tissue pathological changes TOI, local organization hemoglobin
Spend index THI and local brain tissue total hemoglobin concentration variation delta CtHb；
Second detection unit 402, for while first detection unit is run, measuring human body artery pressure；
Analytic unit 403, for use small echo crossbar transistion method, analyze the body local brain tissue pathological changes TOI,
Local organization hemoglobin concentration index THI, local brain tissue total hemoglobin concentration variation delta CtHb and human body artery pressure
Between correlation and phase correlation.
As a kind of optional embodiment, analytic unit 403 is specifically used for：
Continuous wavelet transform process data, on a specific time, each frequency is as a corresponding window
Mouthful, on the window, quadractically integrable function Ψ that mean value is 00The Pan and Zoom of (η) in time can to frequency into
Row partial analysis, Ψ0(η) is known as morther wavelet；
Continuous wavelet transform is expressed as：
S indicates that scale, δ t indicate sampling interval, xn(n=1,2,3 ... N) indicates time series；
The intersection small echo of two different time domain signals Xn and Yn are defined as：
WXY=WXWY*Formula (2)
* indicates complex conjugate in formula, definition | WXY| it is power spectrum, the intersection small echo W of two different time domain signals Xn and YnXY
It is the result of two signal Xn and Yn crossbar transistion；Its absolute value is power spectrum, wherein complex matrix | WXY| phase angle be two
Phase differences of a time series Xn and Yn on time-frequency domain:
Φ in formulai(s) be two time series in time tiPhase delay under scale s.
As a kind of optional embodiment, analytic unit 403 is additionally operable to：
Analyze small echo correlation and Wavelet Phase correlation in different frequency range between brain oxygen parameter and angiosthenia；
Small echo correlation (WCO) formula is：
Wavelet Phase correlation (WPCO) formula is：
Wherein,It is two groups of signal phase differences being averaged in the time domain.
As a kind of optional embodiment, first detection unit 401 is specifically used for：It is examined by Near-infrared Brain oxygen detection
Survey body local brain tissue pathological changes TOI, local organization hemoglobin concentration index THI and local brain tissue total hemoglobin
Two electrode slices specifically can be attached to the forehead of human body under the cooperation of worker, measure 30- by concentration variation delta CtHb
60 minutes, extract body local brain tissue pathological changes TOI and part tissue hemoglobin concentration index THI and local brain tissue
Total hemoglobin concentration variation delta CtHb.
The above embodiment is a preferred embodiment of the present invention, but embodiments of the present invention are not by above-described embodiment
Limitation, it is other it is any without departing from the spirit and principles of the present invention made by changes, modifications, substitutions, combinations, simplifications,
Equivalent substitute mode is should be, is included within the scope of the present invention.
Claims (5)
1. a kind of autonomous regulating power evaluation method of brain blood flow based on brain oxygen signal and blood pressure signal, which is characterized in that described
Method includes：
Step (1), detection body local brain tissue pathological changes TOI, local organization hemoglobin concentration index THI and local cerebral
Organize total hemoglobin concentration variation delta CtHb；
Step (2), while carrying out step (1), measure human body artery pressure；
Step (3), using small echo crossbar transistion method, it is blood red to analyze the body local brain tissue pathological changes TOI, local organization
Correlation between albumen concentration index THI, local brain tissue total hemoglobin concentration variation delta CtHb and human body artery pressure
And phase correlation；
In the small echo correlation and Wavelet Phase correlation of different frequency range between step (4), analysis brain oxygen parameter and angiosthenia.
2. the method as described in claim 1, which is characterized in that
The step (3), including：
Continuous wavelet transform process data, on a specific time, each frequency as a corresponding window,
On the window, the Pan and Zoom of quadractically integrable function Ψ 0 (η) that mean value is 0 in time can be to frequency carry out office
Portion analyzes, Ψ0(η) is known as morther wavelet；
Continuous wavelet transform is expressed as：
S indicates that scale, δ t indicate sampling interval, xn(n=1,2,3 ... N) indicates time series；
The intersection small echo of two different time domain signals Xn and Yn are defined as：
WXY=WXWY*Formula (2)
* indicates complex conjugate in formula, definition | WXY| it is power spectrum, the intersection small echo W of two different time domain signals Xn and YnXYIt is two
The result of signal Xn and Yn crossbar transistion；Its absolute value is power spectrum, wherein complex matrix | WXY| phase angle be two times
The phase difference of sequence X n and Yn on time-frequency domain:
Φ in formulai(s) be two time series in time tiPhase delay under scale s；
The step (4), including：
Small echo correlation (WCO) formula is：
Wavelet Phase correlation (WPCO) formula is：
Wherein,It is two groups of signal phase differences being averaged in the time domain.
3. the autonomous regulating power evaluation side of a kind of brain blood flow based on brain oxygen signal and blood pressure signal according to claim 1
Method, which is characterized in that in the step (1), body local brain tissue pathological changes are detected using Near-infrared Brain oxygen detection
TOI, local organization hemoglobin concentration index THI and local brain tissue total hemoglobin concentration variation delta CtHb；Specially：
Two electrode slices are attached to the forehead of human body, measure 30-60 minute, extraction body local brain tissue pathological changes TOI with
Local organization hemoglobin concentration index THI and local brain tissue total hemoglobin concentration variation delta CtHb.
4. a kind of autonomous regulating power evaluating apparatus of brain blood flow based on brain oxygen signal and blood pressure signal, which is characterized in that described
Device includes：
First detection unit, for detecting body local brain tissue pathological changes TOI, local organization hemoglobin concentration index
THI and local brain tissue total hemoglobin concentration variation delta CtHb；
Second detection unit, for while first detection unit is run, measuring human body artery pressure；
Analytic unit analyzes the body local brain tissue pathological changes TOI, local group for using small echo crossbar transistion method
It knits between hemoglobin concentration index THI, local brain tissue total hemoglobin concentration variation delta CtHb and human body artery pressure
In the small echo correlation and small echo phase of different frequency range between correlation and phase correlation, and analysis brain oxygen parameter and angiosthenia
Position correlation.
5. the autonomous regulating power evaluation dress of a kind of brain blood flow based on brain oxygen signal and blood pressure signal according to claim 4
It sets, which is characterized in that the analytic unit is used for：
Continuous wavelet transform process data, on a specific time, each frequency as a corresponding window,
On the window, quadractically integrable function Ψ that mean value is 00The Pan and Zoom of (η) in time can be to frequency carry out office
Portion analyzes, Ψ0(η) is known as morther wavelet；
Continuous wavelet transform is expressed as：
S indicates that scale, δ t indicate sampling interval, xn(n=1,2,3 ... N) indicates time series；
The intersection small echo of two different time domain signals Xn and Yn are defined as：
WXY=WXWY*Formula (2)
* indicates complex conjugate in formula, definition | WXY| it is power spectrum, the intersection small echo W of two different time domain signals Xn and YnXYIt is two
The result of signal Xn and Yn crossbar transistion；Its absolute value is power spectrum, wherein complex matrix | WXY| phase angle be two times
The phase difference of sequence X n and Yn on time-frequency domain:
Φ in formulai(s) be two time series in time tiPhase delay under scale s；
Analyze small echo correlation and Wavelet Phase correlation in different frequency range between brain oxygen parameter and angiosthenia：
Small echo correlation (WCO) formula is：
Wavelet Phase correlation (WPCO) formula is：
Wherein,It is two groups of signal phase differences being averaged in the time domain.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810481779.1A CN108670238A (en) | 2018-05-18 | 2018-05-18 | The autonomous regulating power evaluation method of brain blood flow based on brain oxygen signal and blood pressure signal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810481779.1A CN108670238A (en) | 2018-05-18 | 2018-05-18 | The autonomous regulating power evaluation method of brain blood flow based on brain oxygen signal and blood pressure signal |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108670238A true CN108670238A (en) | 2018-10-19 |
Family
ID=63806791
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810481779.1A Pending CN108670238A (en) | 2018-05-18 | 2018-05-18 | The autonomous regulating power evaluation method of brain blood flow based on brain oxygen signal and blood pressure signal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108670238A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109730668A (en) * | 2018-11-30 | 2019-05-10 | 上海浩聚医疗科技有限公司 | Signal processing method based on histogram estimation brain self-regulation ability |
| CN112206140A (en) * | 2020-09-03 | 2021-01-12 | 国家康复辅具研究中心 | Pneumatic circulation training method and system for heart and brain coupling guidance |
| CN112826457A (en) * | 2020-12-31 | 2021-05-25 | 苏州爱琴生物医疗电子有限公司 | Brain-based sign information detection device and method and related equipment |
| CN113384256A (en) * | 2021-05-24 | 2021-09-14 | 苏州爱琴生物医疗电子有限公司 | Brain blood flow autonomous regulation capacity assessment method |
| CN116269249A (en) * | 2023-02-14 | 2023-06-23 | 国家康复辅具研究中心 | Cerebral apoplexy risk prediction method and system |
-
2018
- 2018-05-18 CN CN201810481779.1A patent/CN108670238A/en active Pending
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109730668A (en) * | 2018-11-30 | 2019-05-10 | 上海浩聚医疗科技有限公司 | Signal processing method based on histogram estimation brain self-regulation ability |
| CN112206140A (en) * | 2020-09-03 | 2021-01-12 | 国家康复辅具研究中心 | Pneumatic circulation training method and system for heart and brain coupling guidance |
| CN112826457A (en) * | 2020-12-31 | 2021-05-25 | 苏州爱琴生物医疗电子有限公司 | Brain-based sign information detection device and method and related equipment |
| CN113384256A (en) * | 2021-05-24 | 2021-09-14 | 苏州爱琴生物医疗电子有限公司 | Brain blood flow autonomous regulation capacity assessment method |
| CN116269249A (en) * | 2023-02-14 | 2023-06-23 | 国家康复辅具研究中心 | Cerebral apoplexy risk prediction method and system |
| CN116269249B (en) * | 2023-02-14 | 2023-11-17 | 国家康复辅具研究中心 | Cerebral apoplexy risk prediction method and system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108670238A (en) | The autonomous regulating power evaluation method of brain blood flow based on brain oxygen signal and blood pressure signal | |
| CN112274126B (en) | Noninvasive continuous blood pressure detection method and device based on multiple pulse waves | |
| Aysin et al. | Effect of respiration in heart rate variability (HRV) analysis | |
| EP3358485A1 (en) | General noninvasive blood glucose prediction method based on timing analysis | |
| Tayel et al. | Poincaré plot for heart rate variability | |
| US20120016209A1 (en) | Fetal ecg monitoring | |
| CN106264499B (en) | Analysis method for quantifying interaction of heart-lung system | |
| Tadi et al. | Seismocardiography: Toward heart rate variability (HRV) estimation | |
| US20130102856A1 (en) | Non-invasive detection of fetal or maternal illness | |
| Babaeian et al. | Real time driver drowsiness detection using a logistic-regression-based machine learning algorithm | |
| US20090124914A1 (en) | Analysis system and a method for pulse diagnosis in chinese medicine | |
| CN105686810A (en) | A temperature measuring method and device | |
| CN112274120A (en) | Noninvasive arteriosclerosis detection method and device based on one-way pulse wave | |
| Li et al. | The correlation study of Cun, Guan and Chi position based on wrist pulse characteristics | |
| Johnson et al. | Detecting aortic valve-induced abnormal flow with seismocardiography and cardiac MRI | |
| EP1982646A2 (en) | Analysis system and a method for pulse diagnostic in chinese medicine | |
| Chen et al. | Cerebral blood flow autoregulation measurement via bioimpedance technology | |
| CN111000544B (en) | Method and system for constructing hybrid continuous blood pressure measurement model based on PPG waveform | |
| WO2020220403A1 (en) | Fatigue state detection method, apparatus and device, and storage medium | |
| Akşahin et al. | Obstructive sleep apnea classification with artificial neural network based on two synchronic hrv series | |
| CN113171061B (en) | Noninvasive vascular function assessment method, noninvasive vascular function assessment device, noninvasive vascular function assessment equipment and noninvasive vascular function assessment medium | |
| CN106618517B (en) | Method for evaluating artery elasticity function based on pulse wave of arm position change | |
| Kang et al. | Multi bio-signal based algorithm using EMD and FFT for stress analysis | |
| CN110897621B (en) | Physiological signal extraction method and system | |
| Addison | Identifying stable phase coupling associated with cerebral autoregulation using the synchrosqueezed cross-wavelet transform and low oscillation morlet wavelets |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20181019 |
|
| WD01 | Invention patent application deemed withdrawn after publication |