CN104997518A - Motion interference offset method and system for blood oxygen measuring equipment - Google Patents
Motion interference offset method and system for blood oxygen measuring equipment Download PDFInfo
- Publication number
- CN104997518A CN104997518A CN201510364159.6A CN201510364159A CN104997518A CN 104997518 A CN104997518 A CN 104997518A CN 201510364159 A CN201510364159 A CN 201510364159A CN 104997518 A CN104997518 A CN 104997518A
- Authority
- CN
- China
- Prior art keywords
- hongguang
- power data
- receiving light
- light power
- infrared light
- 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.)
- Granted
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/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/1455—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 using optical sensors, e.g. spectral photometrical oximeters
- A61B5/14551—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 using optical sensors, e.g. spectral photometrical oximeters 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/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
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
-
- 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/7221—Determining signal validity, reliability or quality
Abstract
The invention provides a motion interference offset method and system for blood oxygen measuring equipment. The method comprises following steps: acquiring red light reception intensity data and infrared light reception intensity data measured by the blood oxygen measuring equipment in a set period; calculating shape similarity of the red light reception intensity data and the infrared light reception intensity data; when the shape similarity is larger than a set threshold, calculating red light interference parameters according to the red light reception intensity data and preset interference reference data sequences and eliminating the motion interference of the red light reception intensity data according to the red light interference parameters, and calculating infrared light interference parameters according to the infrared light reception intensity data and preset interference reference data sequences and eliminating the motion interference of the infrared light reception intensity data according to the infrared light interference parameters; and calculating blood oxygen values according to the infrared light reception intensity data and the red light reception intensity data with the motion interference being eliminated. By employing the method mentioned above, the motion interference can be accurately determined, the delay time is short, and the motion interference can be further eliminated.
Description
Technical field
The present invention relates to Oximetry technical field, particularly relate to motion artifacts counteracting method and the system of blood oxygen measurement device.
Background technology
Blood oxygen saturation (SaO2) is that the capacity of the HbO2 Oxyhemoglobin (HbO2) combined by oxygen in blood accounts for whole combinative hemoglobin (Hb, hemoglobin) percentage ratio of capacity, the i.e. concentration of blood oxygen in blood, it is the important physiological parameter of breath cycle.
Current measuring method adopts finger cot type photoelectric sensor, during measurement, only need by sensor sleeve on finger, utilize finger as the transparent vessel of splendid attire hemoglobin, the HONGGUANG of use wavelength 660nm and the near infrared light of 940nm, as injecting light source, measure the light conductive strength by tissue bed, calculate hemoglobin concentration and blood oxygen saturation, instrument can show human body blood oxygen saturation, provides a kind of not damaged Oximetry instrument continuously for clinical.But when human body is kept in motion, its blood flow and muscular states can constantly change, and the light penetration region of blood oxygen measurement device also constantly can produce displacement, thus affect result accuracy.Differentiate from Oximetry result, when measurement blood oxygen levels result Rapid Variable Design within a period of time, then think motion artifacts, now SHIELDING CALCULATION result, and motion artifacts of reporting to the police.
Current this method, owing to being judge according to result stability, motion artifacts judges to have larger time delay, meanwhile, only can report to the police, and can not measure result preferably when there is motion artifacts.
Summary of the invention
Based on this, be necessary to provide a kind of and judge that motion artifacts comparatively accurately, time delay is less and better can eliminate the motion artifacts counteracting method of the blood oxygen measurement device of motion artifacts and the motion artifacts bucking-out system of blood oxygen measurement device.
A motion artifacts counteracting method for blood oxygen measurement device, comprises the following steps:
Obtain HONGGUANG receiving light power data and infrared light receiving light power data that in setting-up time section, blood oxygen measurement device is measured;
Calculate the plesiomorphism degree of described HONGGUANG receiving light power data and infrared light receiving light power data, when described plesiomorphism degree is greater than setting threshold value, calculate HONGGUANG interference parameter according to described HONGGUANG receiving light power data and default interference reference data sequence, and eliminate the motion artifacts in HONGGUANG receiving light power data according to described HONGGUANG interference parameter;
Calculate infrared light interference parameter according to infrared light receiving light power data and described interference reference data sequence, and eliminate the motion artifacts in infrared light receiving light power data according to described infrared light interference parameter;
The infrared light receiving light power data after eliminating motion artifacts and HONGGUANG receiving light power data are utilized to calculate blood oxygen levels.
A motion artifacts bucking-out system for blood oxygen measurement device, comprising:
Acquisition module, for obtaining the HONGGUANG receiving light power data and infrared light receiving light power data that in setting-up time section, blood oxygen measurement device is measured;
Eliminate HONGGUANG interference module, for calculating the plesiomorphism degree of described HONGGUANG receiving light power data and infrared light receiving light power data, when described plesiomorphism degree is greater than setting threshold value, calculate HONGGUANG interference parameter according to described HONGGUANG receiving light power data and default interference reference data sequence, and eliminate the motion artifacts in HONGGUANG receiving light power data according to described HONGGUANG interference parameter;
Eliminating infrared light interference module, for calculating infrared light interference parameter according to infrared light receiving light power data and described interference reference data sequence, and eliminating the motion artifacts in infrared light receiving light power data according to described infrared light interference parameter;
Blood oxygen levels computing module, calculates blood oxygen levels for utilizing the infrared light receiving light power data after eliminating motion artifacts and HONGGUANG receiving light power data.
The motion artifacts counteracting method of above-mentioned blood oxygen measurement device and the motion artifacts bucking-out system of blood oxygen measurement device, by calculating the plesiomorphism degree of described HONGGUANG receiving light power data and infrared light receiving light power data, when plesiomorphism degree is greater than setting threshold value, namely be judged to there occurs motion artifacts, judge that motion artifacts is comparatively accurate, time delay is less; Then by calculating HONGGUANG interference parameter and infrared light interference parameter, motion artifacts is eliminated according to two parameters and default interference reference data sequence, then the infrared light receiving light power data after eliminating motion artifacts and HONGGUANG receiving light power data are utilized to calculate blood oxygen levels, thus the motion artifacts can eliminated in blood oxygen measurement device measurement blood oxygen levels, improve measurement accuracy.
Accompanying drawing explanation
Fig. 1 is the motion artifacts counteracting method flow chart of an embodiment blood oxygen measurement device;
Fig. 2 is an embodiment blood oxygen measurement device instrumentation plan;
Fig. 3 is the HONGGUANG of embodiment when not moving interference and infrared light oscillogram;
Fig. 4 is the HONGGUANG of embodiment when moving interference and infrared light oscillogram;
Fig. 5 is that an embodiment eliminates the HONGGUANG after motion artifacts and infrared light oscillogram;
Fig. 6 is the motion artifacts bucking-out system structural representation of an embodiment blood oxygen measurement device.
Detailed description of the invention
Be described in detail below in conjunction with the detailed description of the invention of accompanying drawing to the motion artifacts counteracting method of blood oxygen measurement device of the present invention and the motion artifacts bucking-out system of blood oxygen measurement device.
Refer to Fig. 1, Fig. 1 is the motion artifacts counteracting method flow chart of an embodiment blood oxygen measurement device.
A motion artifacts counteracting method for blood oxygen measurement device, comprises the following steps:
Step S101: obtain HONGGUANG receiving light power data and infrared light receiving light power data that in setting-up time section, blood oxygen measurement device is measured;
In step S101, setting-up time section can be 8s, 10s, 12s etc., sets according to demand.
Step S103: the plesiomorphism degree calculating described HONGGUANG receiving light power data and infrared light receiving light power data, when described plesiomorphism degree is greater than setting threshold value, calculate HONGGUANG interference parameter according to described HONGGUANG receiving light power data and default interference reference data sequence, and eliminate the motion artifacts in HONGGUANG receiving light power data according to described HONGGUANG interference parameter;
In step s 103, the plesiomorphism degree calculating described HONGGUANG receiving light power data and infrared light receiving light power data is due to when not moving interference, the plesiomorphism degree of HONGGUANG receiving light power data and infrared light receiving light power data is higher, and when sending motion artifacts, the plesiomorphism degree of HONGGUANG receiving light power data and infrared light receiving light power data is lower, can judge motion artifacts situation more accurately and compared with low delay by the plesiomorphism degree calculating HONGGUANG receiving light power data and infrared light receiving light power data.
In one embodiment, the step of the plesiomorphism degree of described calculating described HONGGUANG receiving light power data and infrared light receiving light power data can comprise:
The first cost function value is calculated according to HONGGUANG receiving light power data and infrared light receiving light power data, and with the minima of the first cost function value for plesiomorphism degree; Wherein, the formula calculating the first cost function value is:
d
1＝(E
1-b*(D
1-a)-c)
2+(E
2-b*(D
2-a)-c)
2+…+(E
n-b*(D
n-a)-c)
2；
Wherein, d
1represent the first cost function value, E
1, E
2, E
nrepresent infrared light receiving light power data, D
1, D
2, D
nrepresent HONGGUANG receiving light power data, a represents the baseline value of HONGGUANG, and c represents the baseline value of infrared light, and b represents the ratio of the peak valley difference of infrared light and the peak valley difference of HONGGUANG, and n represents data amount check.
When trying to achieve the minima of the first cost function value, the minima of the first cost function value is plesiomorphism degree, meanwhile, namely a, b, c are now similar parameter G.a, G.b, G.c, represent the similarity degree of HONGGUANG receiving light power data and infrared light receiving light power data.
In one embodiment, describedly calculate in the step of HONGGUANG interference parameter according to described HONGGUANG receiving light power data and the interference reference data sequence preset, the interference reference data sequence preset is:
L
1＝E
1-G.b*(D
1-G.a)-G.c；
L
2＝E
2-G.b*(D
2-G.a)-G.c；
……；
L
n＝E
n-G.b*(D
n-G.a)-G.c；
Wherein, G.a, G.b, G.c represent and calculate the numerical value that when plesiomorphism is spent, a, b, c are corresponding.
The interference reference data sequence preset and interfering signal have strong correlation, can be obtained by the mode of statistical measurement.
Further, in one embodiment, the described step according to described HONGGUANG receiving light power data and default interference reference data sequence calculating HONGGUANG interference parameter can comprise:
Calculate the second cost functional minimum value according to described HONGGUANG receiving light power data and default interference reference data sequence, calculate HONGGUANG interference parameter according to the second cost functional minimum value; Wherein, the second cost function is:
d
2＝(D
1-M.b*(L
1-M.a)-M.c)
2+(D
2-M.b*(L
2-M.a)-M.c)
2+…+(D
n-M.b*(L
n-M.a)-M.c)
2；
Wherein, d
2represent the second cost function value, D
1, D
2, D
nrepresent HONGGUANG receiving light power data, L
1, L
2, L
nrepresent interference reference data sequence, M.a, M.b, M.c represent HONGGUANG interference parameter.
The value of n is larger, and the degree of accuracy that follow-up blood oxygen levels is measured is higher, but the complexity of calculating is also higher.
In one embodiment, described according in the motion artifacts step in described HONGGUANG interference parameter elimination HONGGUANG receiving light power data, eliminating HONGGUANG interference formula is:
J
1＝D
1-M.b*(L
1-M.a)-M.c；
J
2＝D
2-M.b*(L
2-M.a)-M.c；
……；
J
n＝D
n-M.b*(L
n-M.a)-M.c；
Wherein, J
1, J
2, J
nrepresent the HONGGUANG receiving light power data after eliminating motion artifacts.
Step S105: calculate infrared light interference parameter according to infrared light receiving light power data and described interference reference data sequence, and eliminate the motion artifacts in infrared light receiving light power data according to described infrared light interference parameter;
In one embodiment, the described step according to described infrared light receiving light power data and described interference reference data sequence calculating infrared light interference parameter can comprise:
Calculate third generation valency functional minimum value according to described infrared light receiving light power data and default interference reference data sequence, calculate infrared light interference parameter according to third generation valency functional minimum value; Wherein, third generation valency function is:
D
3=(E
1-N.b* (L
1-N.a)-N.c)
2+ (E
2-N.b* (L
2-N.a)-N.c)
2+ ... + (E
n-N.b* (L
n-N.a)-N.c)
2; Wherein, d
3represent third generation valency functional value, E
1, E
2, E
nrepresent infrared light receiving light power data, L
1, L
2, L
nrepresent interference reference data sequence, N.a, N.b, N.c represent infrared light interference parameter.
Further, in one embodiment, described according in the motion artifacts step in described infrared light interference parameter elimination infrared light receiving light power data, eliminating infrared light interference formula is:
K
1＝E
1-N.b*(L
1-N.a)-N.c；
K
2＝E
2-N.b*(L
2-N.a)-N.c；
……；
K
n＝E
n-N.b*(L
n-N.a)-N.c；
Wherein, K
1, K
2, K
nrepresent the infrared light receiving light power data after eliminating motion artifacts.
Step S107: utilize the infrared light receiving light power data after eliminating motion artifacts and HONGGUANG receiving light power data to calculate blood oxygen levels.
The motion artifacts counteracting method of above-mentioned blood oxygen measurement device, by calculating the plesiomorphism degree of described HONGGUANG receiving light power data and infrared light receiving light power data, when plesiomorphism degree is greater than setting threshold value, namely there occurs motion artifacts, judge that motion artifacts is comparatively accurate, time delay is less; Then by calculating HONGGUANG interference parameter and infrared light interference parameter, motion artifacts is eliminated according to two parameters and default interference reference data sequence, recycling eliminates the infrared light receiving light power data after motion artifacts and HONGGUANG receiving light power data calculate blood oxygen levels, thus the motion artifacts can effectively eliminated in blood oxygen measurement device measurement blood oxygen levels, improve measurement accuracy.
In order to further describe the motion artifacts counteracting method of blood oxygen measurement device of the present invention in detail, be described below in conjunction with embody rule example.
Refer to Fig. 2, Fig. 2 is an embodiment blood oxygen measurement device instrumentation plan.
Current measuring method adopts finger cot type photoelectric sensor, during measurement, only need by sensor sleeve on finger, as shown in Figure 2, utilize finger as the transparent vessel of splendid attire hemoglobin, the HONGGUANG of use wavelength 660nm and the near infrared light of 940nm are as injecting light source, by ballistic device to finger red-emitting and infrared light, HONGGUANG receiving light power data and infrared light receiving light power data are obtained by receiving device, and the light conductive strength measured by tissue bed, calculate hemoglobin concentration and blood oxygen saturation, instrument can show human body blood oxygen saturation.
Refer to Fig. 3, Fig. 3 is the HONGGUANG of embodiment when not moving interference and infrared light oscillogram.
During without motion artifacts, the optical signal that blood oxygen measurement device measures, the receiving light power of HONGGUANG and infrared light morphologically closely.As shown in Figure 3, upper strata curve is the receiving light power of HONGGUANG, and bottom plot line is the receiving light power of infrared light.
Refer to Fig. 4, Fig. 4 is the HONGGUANG of embodiment when moving interference and infrared light oscillogram.
During motion artifacts, the optical signal that blood oxygen measurement device measures, the receiving light power of HONGGUANG and infrared light morphologically differs greatly.As shown in Figure 4, upper strata curve is the receiving light power of HONGGUANG, and bottom plot line is the receiving light power of infrared light.
The method of Oximetry will be described below, first predetermined movement interference differentiates that duration is A, and namely each one group of data that is A with nearest duration differentiate whether motion artifacts, and A can be 10s.
Note equipment Optical Sampling rate is B, differentiates whether motion artifacts with C=A*B data in the nearest time period namely at every turn.The sequence of nearest C data of note HONGGUANG receiving light power is D, and the sequence of nearest C data of infrared light receiving light power is E.Use plesiomorphism degree algorithm, try to achieve the receiving light power morphologically similarity degree F of HONGGUANG and infrared light, and similar parameter G:
Plesiomorphism degree algorithm is as follows:
D
1=(E
1-b* (D
1-a)-c)
2+ (E
2-b* (D
2-a)-c)
2+ ... + (E
n-b* (D
n-a)-c)
2, wherein, d
1represent the first cost function value, E
1, E
2, E
nrepresent infrared light receiving light power data, D
1, D
2, D
nrepresent HONGGUANG receiving light power data, a represents the baseline value of HONGGUANG, and c represents the baseline value of infrared light, and b represents the ratio of the peak valley difference of infrared light and the peak valley difference of HONGGUANG, and n represents data amount check, the desirable any real number of a, b, c, n=C.
Then try to achieve and can make the first cost function value d
1the value of a, b, c that value is minimum, and d
1minima, d
1minima be similarity degree F in form.Make d obtain a, b, c of minima, be similar parameter G, hereafter will be described as G.a, G.b, G.c respectively.
Predetermined movement interference decision threshold H, if F>H, is then judged as motion artifacts state; If F≤H, be then judged as non-athletic disturbance state, H value is less then more responsive to motion artifacts.If be judged as motion artifacts state, then use motion artifacts cancellation algorithms, by G, D, E, nearest C the data sequence of trying to achieve the HONGGUANG receiving light power after offsetting motion artifacts is J, and nearest C the data sequence of offsetting the infrared light receiving light power after motion artifacts is K.
G used herein can be the similar parameter G that non-motion artifacts state stores, and also can be the G that data J, K after using elimination motion artifacts hereafter try to achieve again through calculating.
Motion artifacts cancellation algorithms is as follows:
d
2＝(D
1-M.b*(L
1-M.a)-M.c)
2+(D
2-M.b*(L
2-M.a)-M.c)
2+…+(D
n-M.b*(L
n-M.a)-M.c)
2；
Wherein, d
2represent the second cost function value, D
1, D
2, D
nrepresent HONGGUANG receiving light power data, L
1, L
2, L
nrepresent interference reference data sequence, M.a, M.b, M.c represent HONGGUANG interference parameter.
Interference reference data sequence L is:
L
1＝E
1-G.b*(D
1-G.a)-G.c；
L
2＝E
2-G.b*(D
2-G.a)-G.c；
……；L
n＝E
n-G.b*(D
n-G.a)-G.c；
Wherein, G.a, G.b, G.c represent and calculate the numerical value that when plesiomorphism is spent, a, b, c are corresponding.
According to L, D, try to achieve HONGGUANG interference parameter M, specifically,
Try to achieve and can make the second cost function value d
2the value of a, b, c that value is minimum, a, b, c value is now HONGGUANG interference parameter M, hereafter will be described as M.a, M.b, M.c respectively.
The HONGGUANG receiving light power data J eliminated after motion artifacts is:
J
1＝D
1-M.b*(L
1-M.a)-M.c；
J
2＝D
2-M.b*(L
2-M.a)-M.c；
……；
J
n＝D
n-M.b*(L
n-M.a)-M.c；
According to L, E, try to achieve infrared light interference parameter N, specifically, design third generation valency function: d
3=(E
1-N.b* (L
1-N.a)-N.c)
2+ (E
2-N.b* (L
2-N.a)-N.c)
2+ ... + (E
n-N.b* (L
n-N.a)-N.c)
2; Wherein, d
3represent third generation valency functional value, E
1, E
2, E
nrepresent infrared light receiving light power data, L
1, L
2, L
nrepresent interference reference data sequence, N.a, N.b, N.c represent infrared light interference parameter.
Try to achieve and can make cost function value d
3the value of a, b, c that value is minimum, this, b, c value is infrared light interference parameter N, hereafter will be described as N.a, N.b, N.c respectively.
The Infrared Data sequence K offset after interference is:
K
1＝E
1-N.b*(L
1-N.a)-N.c；
K
2＝E
2-N.b*(L
2-N.a)-N.c；
……；
K
n＝E
n-N.b*(L
n-N.a)-N.c；
Wherein, K
1, K
2, K
nrepresent the infrared light receiving light power data after eliminating motion artifacts.
Use J to replace D, K to replace E, carry out follow-up blood oxygen levels and calculate, the impact of motion artifacts can be eliminated.
Fig. 5 is the effect eliminating motion artifacts, and visible elimination disturbs front and back oscillogram difference larger.
Refer to Fig. 6, Fig. 6 is the motion artifacts bucking-out system structural representation of an embodiment blood oxygen measurement device.
A motion artifacts bucking-out system for blood oxygen measurement device, comprising:
Acquisition module 210, for obtaining the HONGGUANG receiving light power data and infrared light receiving light power data that in setting-up time section, blood oxygen measurement device is measured;
Eliminate HONGGUANG interference module 230, for calculating the plesiomorphism degree of described HONGGUANG receiving light power data and infrared light receiving light power data, when described plesiomorphism degree is greater than setting threshold value, calculate HONGGUANG interference parameter according to described HONGGUANG receiving light power data and default interference reference data sequence, and eliminate the motion artifacts in HONGGUANG receiving light power data according to described HONGGUANG interference parameter;
Eliminating infrared light interference module 250, for calculating infrared light interference parameter according to infrared light receiving light power data and described interference reference data sequence, and eliminating the motion artifacts in infrared light receiving light power data according to described infrared light interference parameter.
Blood oxygen levels computing module 270, calculates blood oxygen levels for utilizing the infrared light receiving light power data after eliminating motion artifacts and HONGGUANG receiving light power data.
The motion artifacts bucking-out system of above-mentioned blood oxygen measurement device, setting-up time section can be 8s, 10s, 12s etc., set according to demand, the plesiomorphism degree calculating described HONGGUANG receiving light power data and infrared light receiving light power data is due to when not moving interference, the plesiomorphism degree of HONGGUANG receiving light power data and infrared light receiving light power data is higher, and when sending motion artifacts, the plesiomorphism degree of HONGGUANG receiving light power data and infrared light receiving light power data is lower, motion artifacts situation can be judged more accurately and compared with low delay by the plesiomorphism degree calculating HONGGUANG receiving light power data and infrared light receiving light power data.
In one embodiment, the process that described elimination HONGGUANG interference module 230 performs the plesiomorphism degree of described calculating described HONGGUANG receiving light power data and infrared light receiving light power data can be further used for:
The first cost function value is calculated according to HONGGUANG receiving light power data and infrared light receiving light power data, and with the minima of the first cost function value for plesiomorphism degree; Wherein, the formula calculating the first cost function value is:
d
1＝(E
1-b*(D
1-a)-c)
2+(E
2-b*(D
2-a)-c)
2+…+(E
n-b*(D
n-a)-c)
2；
Wherein, d
1represent the first cost function value, E
1, E
2, E
nrepresent infrared light receiving light power data, D
1, D
2, D
nrepresent HONGGUANG receiving light power data, a represents the baseline value of HONGGUANG, and c represents the baseline value of infrared light, and b represents the ratio of the peak valley difference of infrared light and the peak valley difference of HONGGUANG, and n represents data amount check.
When trying to achieve the minima of the first cost function value, the minima of the first cost function value is plesiomorphism degree, meanwhile, namely a, b, c are now similar parameter G.a, G.b, G.c, represent the similarity degree of HONGGUANG receiving light power data and infrared light receiving light power data.
In one embodiment, described elimination HONGGUANG interference module 230 performs and describedly calculates in the process of HONGGUANG interference parameter according to described HONGGUANG receiving light power data and the interference reference data sequence preset, and the interference reference data sequence preset is:
L
1＝E
1-G.b*(D
1-G.a)-G.c；
L
2＝E
2-G.b*(D
2-G.a)-G.c；
……；
L
n＝E
n-G.b*(D
n-G.a)-G.c；
Wherein, G.a, G.b, G.c represent and calculate the numerical value that when plesiomorphism is spent, a, b, c are corresponding.
The interference reference data sequence preset and interfering signal have strong correlation, can be obtained by the mode of statistical measurement.
Further, in one embodiment, described elimination HONGGUANG interference module 230 performs the described process calculating HONGGUANG interference parameter according to described HONGGUANG receiving light power data and the interference reference data sequence preset and can be further used for:
Calculate the second cost functional minimum value according to described HONGGUANG receiving light power data and default interference reference data sequence, calculate HONGGUANG interference parameter according to the second cost functional minimum value; Wherein, the second cost function is:
d
2＝(D
1-M.b*(L
1-M.a)-M.c)
2+(D
2-M.b*(L
2-M.a)-M.c)
2+…+(D
n-M.b*(L
n-M.a)-M.c)
2；
Wherein, d
2represent the second cost function value, D
1, D
2, D
nrepresent HONGGUANG receiving light power data, L
1, L
2, L
nrepresent interference reference data sequence, M.a, M.b, M.c represent HONGGUANG interference parameter.
The value of n is larger, and the degree of accuracy that follow-up blood oxygen levels is measured is higher, but the complexity of calculating is also higher.
In one embodiment, described elimination HONGGUANG interference module 230 performs described according in the motion artifacts process in described HONGGUANG interference parameter elimination HONGGUANG receiving light power data, and eliminating HONGGUANG interference formula is:
J
1=D
1-M.b* (L
1-M.a)-M.c; J
2=D
2-M.b* (L
2-M.a)-M.c; J
n=D
n-M.b* (L
n-M.a)-M.c; Wherein, J
1, J
2, J
nrepresent the HONGGUANG receiving light power data after eliminating motion artifacts.
In one embodiment, described elimination infrared light interference module 250 performs the described process calculating infrared light interference parameter according to described infrared light receiving light power data and described interference reference data sequence and can be further used for:
Calculate third generation valency functional minimum value according to described infrared light receiving light power data and default interference reference data sequence, calculate infrared light interference parameter according to third generation valency functional minimum value; Wherein, third generation valency function is:
D
3=(E
1-N.b* (L
1-N.a)-N.c)
2+ (E
2-N.b* (L
2-N.a)-N.c)
2+ ... + (E
n-N.b* (L
n-N.a)-N.c)
2; Wherein, d
3represent third generation valency functional value, E
1, E
2, E
nrepresent infrared light receiving light power data, L
1, L
2, L
nrepresent interference reference data sequence, N.a, N.b, N.c represent infrared light interference parameter.
Further, in one embodiment, described elimination infrared light interference module 250 performs described according in the motion artifacts process in described infrared light interference parameter elimination infrared light receiving light power data, and eliminating infrared light interference formula is:
K
1=E
1-N.b* (L
1-N.a)-N.c; K
2=E
2-N.b* (L
2-N.a)-N.c; K
n=E
n-N.b* (L
n-N.a)-N.c; Wherein, K
1, K
2, K
nrepresent the infrared light receiving light power data after eliminating motion artifacts.
The motion artifacts bucking-out system of above-mentioned blood oxygen measurement device, by calculating the plesiomorphism degree of described HONGGUANG receiving light power data and infrared light receiving light power data, when plesiomorphism degree is greater than setting threshold value, namely there occurs motion artifacts, judge that motion artifacts is comparatively accurate, time delay is less; Then by calculating HONGGUANG interference parameter and infrared light interference parameter, motion artifacts is eliminated according to two parameters and default interference reference data sequence, the infrared light receiving light power data after eliminating motion artifacts and HONGGUANG receiving light power data are utilized to calculate blood oxygen levels, thus the motion artifacts can eliminated in blood oxygen measurement device measurement, improve measurement accuracy.
The motion artifacts bucking-out system of blood oxygen measurement device of the present invention and the motion artifacts counteracting method one_to_one corresponding of blood oxygen measurement device of the present invention, the technical characteristic of setting forth in the embodiment of the motion artifacts counteracting method of above-mentioned blood oxygen measurement device and beneficial effect thereof are all applicable to, in the embodiment of motion artifacts bucking-out system of blood oxygen measurement device, hereby state.
The above embodiment only have expressed several embodiment of the present invention, and it describes comparatively concrete and detailed, but therefore can not be interpreted as the restriction to the scope of the claims of the present invention.It should be pointed out that for the person of ordinary skill of the art, without departing from the inventive concept of the premise, can also make some distortion and improvement, these all belong to protection scope of the present invention.Therefore, the protection domain of patent of the present invention should be as the criterion with claims.
Claims (10)
1. a motion artifacts counteracting method for blood oxygen measurement device, is characterized in that, comprise the following steps:
Obtain HONGGUANG receiving light power data and infrared light receiving light power data that in setting-up time section, blood oxygen measurement device is measured;
Calculate the plesiomorphism degree of described HONGGUANG receiving light power data and infrared light receiving light power data, when described plesiomorphism degree is greater than setting threshold value, calculate HONGGUANG interference parameter according to described HONGGUANG receiving light power data and default interference reference data sequence, and eliminate the motion artifacts in HONGGUANG receiving light power data according to described HONGGUANG interference parameter;
Calculate infrared light interference parameter according to infrared light receiving light power data and described interference reference data sequence, and eliminate the motion artifacts in infrared light receiving light power data according to described infrared light interference parameter;
The infrared light receiving light power data after eliminating motion artifacts and HONGGUANG receiving light power data are utilized to calculate blood oxygen levels.
2. the motion artifacts counteracting method of blood oxygen measurement device according to claim 1, is characterized in that, the step of the plesiomorphism degree of described calculating described HONGGUANG receiving light power data and infrared light receiving light power data comprises:
The first cost function value is calculated according to HONGGUANG receiving light power data and infrared light receiving light power data, and with the minima of the first cost function value for plesiomorphism degree; Wherein, the formula calculating the first cost function value is:
d
1＝(E
1-b*(D
1-a)-c)
2+(E
2-b*(D
2-a)-c)
2+…+(E
n-b*(D
n-a)-c)
2；
Wherein, d
1represent the first cost function value, E
1, E
2, E
nrepresent infrared light receiving light power data, D
1, D
2, D
nrepresent HONGGUANG receiving light power data, a represents the baseline value of HONGGUANG, and c represents the baseline value of infrared light, and b represents the ratio of the peak valley difference of infrared light and the peak valley difference of HONGGUANG, and n represents data amount check.
3. the motion artifacts counteracting method of blood oxygen measurement device according to claim 2, it is characterized in that, describedly calculate in the step of HONGGUANG interference parameter according to described HONGGUANG receiving light power data and the interference reference data sequence preset, the interference reference data sequence preset is:
L
1＝E
1-G.b*(D
1-G.a)-G.c；
L
2＝E
2-G.b*(D
2-G.a)-G.c；
……；
L
n＝E
n-G.b*(D
n-G.a)-G.c；
Wherein, G.a, G.b, G.c represent and calculate the numerical value that when plesiomorphism is spent, a, b, c are corresponding.
4. the motion artifacts counteracting method of blood oxygen measurement device according to claim 3, is characterized in that, the described step according to described HONGGUANG receiving light power data and default interference reference data sequence calculating HONGGUANG interference parameter comprises:
Calculate the second cost functional minimum value according to described HONGGUANG receiving light power data and default interference reference data sequence, calculate HONGGUANG interference parameter according to the second cost functional minimum value; Wherein, the second cost function is:
d
2＝(D
1-M.b*(L
1-M.a)-M.c)
2+(D
2-M.b*(L
2-M.a)-M.c)
2+…+(D
n-M.b*(L
n-M.a)-M.c)
2；
Wherein, d
2represent the second cost function value, D
1, D
2, D
nrepresent HONGGUANG receiving light power data, L
1, L
2, L
nrepresent interference reference data sequence, M.a, M.b, M.c represent HONGGUANG interference parameter.
5. the motion artifacts counteracting method of blood oxygen measurement device according to claim 4, is characterized in that, described according in the motion artifacts step in described HONGGUANG interference parameter elimination HONGGUANG receiving light power data, eliminating HONGGUANG interference formula is:
J
1＝D
1-M.b*(L
1-M.a)-M.c；
J
2＝D
2-M.b*(L
2-M.a)-M.c；
……；
J
n＝D
n-M.b*(L
n-M.a)-M.c；
Wherein, J
1, J
2, J
nrepresent the HONGGUANG receiving light power data after eliminating motion artifacts.
6. the motion artifacts counteracting method of blood oxygen measurement device according to claim 3, is characterized in that, the described step according to described infrared light receiving light power data and described interference reference data sequence calculating infrared light interference parameter comprises:
Calculate third generation valency functional minimum value according to described infrared light receiving light power data and default interference reference data sequence, calculate infrared light interference parameter according to third generation valency functional minimum value; Wherein, third generation valency function is:
D
3=(E
1-N.b* (L
1-N.a)-N.c)
2+ (E
2-N.b* (L
2-N.a)-N.c)
2+ ... + (E
n-N.b* (L
n-N.a)-N.c)
2; Wherein, d
3represent third generation valency functional value, E
1, E
2, E
nrepresent infrared light receiving light power data, L
1, L
2, L
nrepresent interference reference data sequence, N.a, N.b, N.c represent infrared light interference parameter.
7. the motion artifacts counteracting method of blood oxygen measurement device according to claim 6, is characterized in that, described according in the motion artifacts step in described infrared light interference parameter elimination infrared light receiving light power data, eliminating infrared light interference formula is:
K
1＝E
1-N.b*(L
1-N.a)-N.c；
K
2＝E
2-N.b*(L
2-N.a)-N.c；
……；
K
n＝E
n-N.b*(L
n-N.a)-N.c；
Wherein, K
1, K
2, K
nrepresent the infrared light receiving light power data after eliminating motion artifacts.
8. a motion artifacts bucking-out system for blood oxygen measurement device, is characterized in that, comprising:
Acquisition module, for obtaining the HONGGUANG receiving light power data and infrared light receiving light power data that in setting-up time section, blood oxygen measurement device is measured;
Eliminate HONGGUANG interference module, for calculating the plesiomorphism degree of described HONGGUANG receiving light power data and infrared light receiving light power data, when described plesiomorphism degree is greater than setting threshold value, calculate HONGGUANG interference parameter according to described HONGGUANG receiving light power data and default interference reference data sequence, and eliminate the motion artifacts in HONGGUANG receiving light power data according to described HONGGUANG interference parameter;
Eliminating infrared light interference module, for calculating infrared light interference parameter according to infrared light receiving light power data and described interference reference data sequence, and eliminating the motion artifacts in infrared light receiving light power data according to described infrared light interference parameter;
Blood oxygen levels computing module, calculates blood oxygen levels for utilizing the infrared light receiving light power data after eliminating motion artifacts and HONGGUANG receiving light power data.
9. the motion artifacts bucking-out system of blood oxygen measurement device according to claim 8, it is characterized in that, the process that described elimination HONGGUANG interference module performs the plesiomorphism degree of described calculating described HONGGUANG receiving light power data and infrared light receiving light power data is further used for:
The first cost function value is calculated according to HONGGUANG receiving light power data and infrared light receiving light power data, and with the minima of the first cost function value for plesiomorphism degree; Wherein, the formula calculating the first cost function value is:
d
1＝(E
1-b*(D
1-a)-c)
2+(E
2-b*(D
2-a)-c)
2+…+(E
n-b*(D
n-a)-c)
2；
Wherein, d
1represent the first cost function value, E
1, E
2, E
nrepresent infrared light receiving light power data, D
1, D
2, D
nrepresent HONGGUANG receiving light power data, a represents the baseline value of HONGGUANG, and c represents the baseline value of infrared light, and b represents the ratio of the peak valley difference of infrared light and the peak valley difference of HONGGUANG, and n represents data amount check.
10. the motion artifacts bucking-out system of blood oxygen measurement device according to claim 9, it is characterized in that, described elimination HONGGUANG interference module execution is described to be calculated in the process of HONGGUANG interference parameter according to described HONGGUANG receiving light power data and default interference reference data sequence, and the interference reference data sequence preset is:
L
1＝E
1-G.b*(D
1-G.a)-G.c；
L
2＝E
2-G.b*(D
2-G.a)-G.c；
……；
L
n＝E
n-G.b*(D
n-G.a)-G.c；
Wherein, G.a, G.b, G.c represent and calculate the numerical value that when plesiomorphism is spent, a, b, c are corresponding.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510364159.6A CN104997518B (en) | 2015-06-25 | 2015-06-25 | The motion artifacts counteracting method and system of blood oxygen measurement device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510364159.6A CN104997518B (en) | 2015-06-25 | 2015-06-25 | The motion artifacts counteracting method and system of blood oxygen measurement device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104997518A true CN104997518A (en) | 2015-10-28 |
| CN104997518B CN104997518B (en) | 2019-03-26 |
Family
ID=54370586
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510364159.6A Active CN104997518B (en) | 2015-06-25 | 2015-06-25 | The motion artifacts counteracting method and system of blood oxygen measurement device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104997518B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107595296A (en) * | 2017-09-19 | 2018-01-19 | 广东乐心医疗电子股份有限公司 | Blood oxygen saturation detection method and blood oxygen detection device for eliminating motion interference |
| CN109009159A (en) * | 2018-07-24 | 2018-12-18 | 华中科技大学 | A kind of anti motion interference accelerating circuit suitable for wearable blood oxygen chip |
| CN112043287A (en) * | 2020-09-30 | 2020-12-08 | 重庆大学 | Noninvasive monitoring method and monitoring device for cerebral blood oxygen |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100061596A1 (en) * | 2008-09-05 | 2010-03-11 | Varian Medical Systems Technologies, Inc. | Video-Based Breathing Monitoring Without Fiducial Tracking |
| CN101836863A (en) * | 2009-03-19 | 2010-09-22 | 深圳迈瑞生物医疗电子股份有限公司 | Method and system for monitoring patients by utilizing two channels |
-
2015
- 2015-06-25 CN CN201510364159.6A patent/CN104997518B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100061596A1 (en) * | 2008-09-05 | 2010-03-11 | Varian Medical Systems Technologies, Inc. | Video-Based Breathing Monitoring Without Fiducial Tracking |
| CN101836863A (en) * | 2009-03-19 | 2010-09-22 | 深圳迈瑞生物医疗电子股份有限公司 | Method and system for monitoring patients by utilizing two channels |
Non-Patent Citations (1)
| Title |
|---|
| 龚渝顺等: "动态环境血氧饱和度监测的运动干扰分离自适应对消方法", 《航天医学与医学工程》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107595296A (en) * | 2017-09-19 | 2018-01-19 | 广东乐心医疗电子股份有限公司 | Blood oxygen saturation detection method and blood oxygen detection device for eliminating motion interference |
| CN107595296B (en) * | 2017-09-19 | 2020-06-05 | 广东乐心医疗电子股份有限公司 | Blood oxygen saturation detection method and blood oxygen detection device for eliminating motion interference |
| CN109009159A (en) * | 2018-07-24 | 2018-12-18 | 华中科技大学 | A kind of anti motion interference accelerating circuit suitable for wearable blood oxygen chip |
| CN109009159B (en) * | 2018-07-24 | 2020-06-02 | 华中科技大学 | Anti-motion interference accelerating circuit suitable for wearable blood oxygen chip |
| CN112043287A (en) * | 2020-09-30 | 2020-12-08 | 重庆大学 | Noninvasive monitoring method and monitoring device for cerebral blood oxygen |
| CN112043287B (en) * | 2020-09-30 | 2021-07-20 | 重庆大学 | Noninvasive monitoring method and monitoring device for cerebral blood oxygen |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104997518B (en) | 2019-03-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107157492B (en) | Embedded human physiological information noninvasive detection system and data processing method | |
| CN106413534B (en) | Continuous blood pressure measuring device, measuring model establishing method and system | |
| CN101991410B (en) | Pulse rate searching and calculating method | |
| AU2013204274B2 (en) | Cardiac activation time detection | |
| CN104173043A (en) | Electrocardiogram (ECG) data analysis method suitable for mobile platform | |
| CN105997043B (en) | A kind of pulse frequency extracting method based on wrist wearable device | |
| KR20160094265A (en) | Computing device and method for analyzing biophysiological rate features thereof | |
| Datta et al. | Blood pressure estimation from photoplethysmogram using latent parameters | |
| US10357164B2 (en) | Method and device for non-invasive blood pressure measurement | |
| CN112545472B (en) | PPG signal quality evaluation method, device, equipment and storage medium | |
| CN104224197A (en) | Calculation method for arterial oxygen saturation by using characteristic quantity of dual-wavelength photoelectric volume pulse wave | |
| CN112233813A (en) | Non-contact non-invasive heart rate and respiration measurement method and system based on PPG | |
| CN103263271A (en) | Non-contact automatic blood oxygen saturation degree measurement system and measurement method | |
| Sun et al. | PPG signal motion artifacts correction algorithm based on feature estimation | |
| CN104665803A (en) | Atrial fibrillation detecting system based on intelligent platform | |
| CN108937957A (en) | Detection method, device and detection device | |
| CN104997518A (en) | Motion interference offset method and system for blood oxygen measuring equipment | |
| CN111166306A (en) | Physiological signal acquisition method, computer device and storage medium | |
| CN104434312B (en) | Custodial care facility and its physiological parameter processing method and system | |
| CN115177215A (en) | Neurovascular coupling analysis method and device for peripheral nerves and electronic equipment | |
| Philip et al. | Continous monitoring of blood glucose using photophlythesmograph signal | |
| CN208511016U (en) | A kind of vital sign comprehensive detection analysis system | |
| CN204618246U (en) | A kind of non-contact vein blood oxygen saturation imaging device | |
| CN103876731B (en) | A kind of Fetal ECG signal extracting device and method | |
| CN104887246B (en) | Blood oxygen measuring method and blood oxygen measuring device |
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 |