CN108784675A - Can eliminate or avoid offset interference and can accurate measurement physiological characteristic electronic device and method - Google Patents
Can eliminate or avoid offset interference and can accurate measurement physiological characteristic electronic device and method Download PDFInfo
- Publication number
- CN108784675A CN108784675A CN201710290128.XA CN201710290128A CN108784675A CN 108784675 A CN108784675 A CN 108784675A CN 201710290128 A CN201710290128 A CN 201710290128A CN 108784675 A CN108784675 A CN 108784675A
- Authority
- CN
- China
- Prior art keywords
- signal
- offset correction
- offset
- physiological characteristic
- correction amount
- 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
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000005259 measurement Methods 0.000 title claims description 4
- 238000012937 correction Methods 0.000 claims abstract description 106
- 239000000203 mixture Substances 0.000 claims abstract description 31
- 230000035479 physiological effects, processes and functions Effects 0.000 claims abstract description 31
- 238000004020 luminiscence type Methods 0.000 claims abstract description 14
- 230000004044 response Effects 0.000 claims abstract description 6
- 238000012360 testing method Methods 0.000 claims description 69
- 238000012545 processing Methods 0.000 claims description 22
- 238000001228 spectrum Methods 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 5
- 230000000873 masking effect Effects 0.000 claims description 4
- 235000013399 edible fruits Nutrition 0.000 claims description 2
- 230000009471 action Effects 0.000 description 20
- 230000003287 optical effect Effects 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 210000004369 blood Anatomy 0.000 description 3
- 239000008280 blood Substances 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000005611 electricity Effects 0.000 description 3
- 230000033001 locomotion Effects 0.000 description 3
- 210000000707 wrist Anatomy 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000001121 heart beat frequency Effects 0.000 description 2
- 238000013186 photoplethysmography Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 210000003743 erythrocyte Anatomy 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000011514 reflex Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
A kind of method measuring physiological characteristic, including：Using an at least luminescence unit, to emit an at least light, which includes to correspond to the first light composition of an one of first wave length and correspond to the second light composition of an one of second wave length, which is different from the second wave length；Using an image sensing circuit, in response in an at least light, sensing generates an at least physiology characteristic quantity and surveys signal；Signal is surveyed for an at least physiology characteristic quantity and carries out an offset correction operation, signal is measured to generate the physiological characteristic after at least one correction；And signal is measured to the physiological characteristic after at least one correction and is calculated, a physiology characteristic results are obtained with estimation.
Description
Technical field
The present disclosure generally relates to a kind of physiological characteristics to measure mechanism, espespecially a kind of electronic device measuring physiological characteristic and side
Method.
Background technology
In general, traditional photo-plethysmographic method (Photoplethysmography, PPG) detection technology uses more waves
Long detecting heartbeat signal, but it is highly prone to action artifact (motion artifact), optical sensing system itself and/or outside
The interference effect of noise, the estimation that those interference easily lead to PPG signals generate bias effect, cause heart beat detection mistake, even if
Traditional technology directly takes the PPG signals of two different wave lengths to do ratio proccessing, can mitigate the interference of partial noise slightly, however
Always the action artifact of the direct current offset caused by light reflection, optical sensing system itself and/or external noise can not be solved
Interference.
Invention content
Therefore, the first purpose of the present invention is to provide a kind of electronic device and method measuring physiological characteristic, to solve
Prior art problem encountered.
According to embodiments of the present invention, there is disclosed a kind of electronic devices measuring physiological characteristic.The electronic device includes one
Image sensing circuit, at least a luminescence unit and a processing circuit, at least a luminescence unit are to emit an at least light, this is extremely
Few light include correspond to the first light composition of an one of first wave length and corresponding to the second light of an one of second wave length at
Part, which is different from the second wave length.Image sensing circuit can be in response in an at least light, and sensing generates at least one
Physiological characteristic measures signal.Processing circuit is coupled to the image sensing circuit, and to be surveyed for an at least physiology characteristic quantity
Signal carries out an offset correction operation, measures signal to generate the physiological characteristic after at least one correction, and to an at least school
Physiological characteristic after just measures signal and is calculated, and a physiology characteristic results are obtained with estimation.
According to embodiments of the present invention, a kind of method measuring physiological characteristic is separately disclosed.This method includes：Use at least one
Luminescence unit, to emit an at least light, which includes to correspond to the first light composition of an one of first wave length
And corresponding to the second light composition of an one of second wave length, which is different from the second wave length；Use an image sensing
Circuit, in response in an at least light, sensing generates an at least physiology characteristic quantity and surveys signal；For an at least physiology characteristic quantity
It surveys signal and carries out an offset correction operation, signal is measured to generate the physiological characteristic after at least one correction；And to this at least one
Physiological characteristic after correction measures signal and is calculated, and a physiology characteristic results are obtained with estimation.
According to embodiments of the present invention, major advantage is to solve the problems, such as traditional technology, improves photo-plethysmographic method
Sensor resistance acts the ability of artifact, increases the accuracy rate of heart beat detection, the offset adjustment of multi-wavelength can be directed to, to improve the heart
Jump detecting is influenced by action artifact, such as the photo-plethysmographic method signal by taking two different wave lengths, removes respectively dynamic
Making the direct current offset of noise influences, then carries out ratio proccessing or act to make an uproar to two different wave length signals after removal direct current offset
Sound deducts skill, can be effectively improved the influence of action artifact.
Description of the drawings
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technology description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
The structure shown according to these attached drawings obtains other attached drawings.
Fig. 1 be one embodiment of the invention to measure physiological characteristic electronic device 100 schematic diagram.
Fig. 2 is the FFT spectrum diagrams for the example that the electronic device of Fig. 1 carries out the test result signal of offset correction.
Fig. 3 is the FFT spectrum signal for another example that the electronic device of Fig. 1 carries out the test result signal of offset correction
Figure.
Fig. 4 is the brief operation workflow block diagram of the processing circuit shown in Fig. 1.
Drawing reference numeral explanation：
| Label | Title | Label | Title |
| 100 | Electronic device | 120 | Processing circuit |
| 105 | Image sensing circuit | 1201 | Offset determining means |
| 110 | Luminescence unit | 1202 | Action noise corrects unit |
| 115 | Light source controller | 1203 | Physiological characteristic estimates unit |
The embodiments will be further described with reference to the accompanying drawings for the realization, the function and the advantages of the object of the present invention.
Specific implementation mode
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation describes, it is clear that described embodiment is only a part of the embodiment of the present invention, instead of all the embodiments.Base
Embodiment in the present invention, those of ordinary skill in the art obtained without creative efforts it is all its
His embodiment, shall fall within the protection scope of the present invention.
If it is to be appreciated that related in the embodiment of the present invention directionality instruction (such as up, down, left, right, before and after ...),
Then directionality instruction be only used for explaining relative position relation under a certain particular pose (as shown in the picture) between each component,
Motion conditions etc., if the particular pose changes, directionality instruction also correspondingly changes correspondingly.
If in addition, relating to the description of " first ", " second " etc. in the embodiment of the present invention, it is somebody's turn to do " first ", " second " etc.
Description be used for description purposes only, be not understood to indicate or imply its relative importance or implicitly indicate indicated skill
The quantity of art feature." first " is defined as a result, the feature of " second " can explicitly or implicitly include at least one spy
Sign.In addition, the technical solution between each embodiment can be combined with each other, but must be with those of ordinary skill in the art's energy
It is enough realize based on, when the knot that conflicting or cannot achieve when will be understood that this technical solution occurs in the combination of technical solution
Conjunction is not present, also not the present invention claims protection domain within.
The spirit of the present invention is to be reduced or avoided action artifact, optical sensing system itself and/or external noise to light
Bias effect caused by the estimation of plethysmography signal, action artifact is mainly started to be made because acting movement in user's generation
At the inaccuracy that physiological characteristic senses, such as when the wrist type electronic device that user can detect one physiological characteristic is worn on it
When wrist, the slight typing action of user can influence the wrist type electronic device completely causes the sensing of physiological characteristic inaccurate
Really, the offset caused by optical sensing system itself is mainly started due to having narrow difference when each optical sensing apparatus is in assembling
Cause the sensing of physiological characteristic inaccurate, in addition, being to enter when light enters skin reflex one of the reason of external noise
The external noise felt by sensor when sensor, the external noise will also cause the sensing of physiological characteristic inaccurate.Following
The electronic device of invention and corresponding method can be reduced or avoided action artifact, optical sensing system itself and/or outside and make an uproar
The bias effect caused by the estimation of photo-plethysmographic method signal is especially reduced or avoided to photo-plethysmographic method signal in sound
Estimation caused by direct current offset influence, by estimating and correcting the offset of an at least photo-plethysmographic method signal, enable follow-up
It can get when estimating physiological characteristic (such as palmic rate etc.) using the photo-plethysmographic method signal after amendment more accurate raw
Feature sensing result is managed, is not interfered by action artifact, optical sensing system itself and/or external noise.Described below
The embodiment of invention.
Please refer to Fig. 1, depicted is one embodiment of the invention measuring the signal of the electronic device 100 of physiological characteristic
Figure.Electronic device 100 include an image sensing circuit 105 (such as an image sensor array), an at least luminescence unit 110,
One light source controller 115 and a processing circuit 120 can be a portable electronic devices or a wearing formula device (example in implementation
Such as smartwatch or bracelet), processing circuit 120 include offset determining means 1201, action noise correction unit 1202 and
Physiological characteristic estimates unit 1203.
Such as LED controller of light source controller 115, to control an at least luminescence unit 110, this is extremely
A few such as light emitting diode of luminescence unit 110, to emit an at least light to user's human skin table of electronic device 100
Face enables the red blood cell composition of user's blood of human body absorb the light emitted, is reflected because of the absorptivity of different light difference
Image in the image sensing circuit 105, with obtain different light photo-plethysmographic method (Photoplethysmography,
PPG) signal, accordingly estimation calculate the physiological characteristic result (such as palmic rate etc.) of user.In implementation, at least one luminous list
Member 110 for example comprising a weight-light ligh-emitting diode unit or includes the light emitting diodes of the multiple and different wavelength of light of correspondence
(such as feux rouges and green light LED unit, the limitation of this non-this case), the light which emits
Line includes the light composition with multiple and different wavelength of light, and the light emitting diode of multiple different wave length is to respectively
In different time points independent transmission respectively the light composition with different wavelength of light to human skin surface；It is above-mentioned this at least one
The possible embodiment of luminescence unit 110 each falls within the scope of this case.
Image sensing circuit 105 generates different body of light products and retouches to sense with the light composition of different wavelength of light
Notation signal, in implementation, if in 110 implementation of luminescence unit only including weight-light ligh-emitting diode unit, accordingly image
Sensing circuit 105 can use a chromatic filter array (Color Filter Array, CFA) 1051, the chromatic filter number
Group 1051 includes multiple Color Channels (color channel) corresponding to different wavelength of light, is filtered out respectively with not
With the light composition of wavelength of light, such as the light composition of red light wavelength, green wavelength, blue light wavelength, image sense are filtered out respectively
Slowdown monitoring circuit 105 generates different body of light products further according to red light wavelength, the light composition of green wavelength, blue light wavelength, respectively imaging
Graphical method signal.In addition, if in 110 implementation of luminescence unit including corresponding different wavelength of light (such as first wave length
The green light of feux rouges and second wave length) multiple light emitting diodes, then accordingly image sensing circuit 105 without colour
Filter array, image sensing circuit 105 can be made in first time point according to the light produced by red light-emitting diode
The result signal that is imaged generates the photo-plethysmographic method signal of red light wavelength after user's blood of human body absorbs, and at second
Between point according to the light produced by green light LED, the result signal that is imaged after absorbed by user's blood of human body produces
The photo-plethysmographic method signal of raw green wavelength, can generate different light respectively according to the light composition of different wave length in this way
Plethysmography signal；Therefore, different in response to the design in the luminescence unit 110, chromatic filter array 1051 is optional
(optional), the limitation of non-present invention.
Above-mentioned measurement to photo-plethysmographic method signal can be used to estimate physiological characteristic, as physiology of the present invention is special
Sign measures signal.It should be noted that above-mentioned each photo-plethysmographic method signal may shift or not deviate, therefore in generation
After the physiological characteristic stated measures signal, in order to which action artifact, optical sensing system itself and/or external noise pair is reduced or avoided
Bias effect caused by the estimation of photo-plethysmographic method signal, processing circuit 120 receive an above-mentioned at least physiology characteristic quantity and survey
Signal simultaneously corrects the offset that an above-mentioned at least physiology characteristic quantity surveys signal, generates the physiological characteristic after correction and measures signal, then root
Signal is measured according to the physiological characteristic after correction to estimate a physiology characteristic results (such as palmic rate of user).Place described below
Manage the running of circuit 120.
In implementation, processing circuit 120 includes offset determining means 1201, action noise correction unit 1202 and physiology
Feature estimates unit 1203, and for a physiology characteristic quantity surveys signal (that is, PPG signals), offset determining means 1201 is used
To generate or estimate an offset correction amount of the PPG signals, action noise corrects unit 1202 to according to the offset correction amount
It corrects the offset of the PPG signals, generates the PPG signals after correction, physiological characteristic estimates unit 1203 according to the correction later
Rear PPG signals estimate such as user's palmic rate.The combination that software, hardware or software and hardware can be used in aforesaid operations comes in fact
It is existing, that is, the combination of software or hardware or software and hardware can be used to realize in above-mentioned processing circuit 120 and its unit 1201 to 1203,
For implemented in software, unit 1201 to 1203 is the software unit for being contained in a specific program code, and processing circuit 120 is for example
It is operated for a processor and to execute the specific program code to carry out above-mentioned offset correction operation and physiological characteristic estimation.
Hardware-implemented words, then said units different circuit units can be used to realize it；The realization of the above software/hardware combination is equal
It belongs to the scope of the present invention.
Offset determining means 1201 generates or the running of one of one PPG signals of estimation offset correction amount is as follows.Processing electricity
Road 120 receives such as three PPG signals, such as corresponds respectively to red light wavelength channel, green wavelength channel and blue light wavelength letter
The PPG signals in road (RGB channels, that is, different light compositions), the present embodiment, in order to accurately estimate DC offset, processing electricity
It (is, for example, the PPG of red light wavelength and green wavelength that road 120, which can be used and correspond respectively to two PPG signals of two different wave length channels,
The PPG signals of signal or blue light wavelength and green wavelength), physiological characteristic estimation operation is carried out, therefore only need to two PPG
Signal carries out offset correction operation, does not need to carry out offset correction to all PPG signals；So this is not the limit of the present invention
System, other embodiment also can carry out offset correction to single PPG signals or all PPG signals.
For example, the DC offset caused by the yield value for correction process variation or respective pixel, processing electricity
Road 120 is by first capturing completely black background video, and it is inclined to be considered as the direct current for the sensing value of generation corresponding to respective pixel at this time
Shifting amount, also that is, processing circuit 120 controls the equivalent environment in masking ambient light of electronic device 100, then in masking ambient light
When measure a PPG signals, the DC offset of the PPG signals under total darkness environment is estimated, as the offset for correcting the direct current offset
Correcting value corrects the DC offset caused by the yield value of process variation or respective pixel accordingly.
And for the direct current offset of corrective action artifact, it lies in and does not cover the action puppet that ambient light for example applies specific frequency
It being carried out when shadow, electronic device 100 is for example placed on the skin of user and is applied with specific frequency or the action artifact of interference,
Offset determining means 1201 can be by gradually applying or adjusting corresponding to multiple and different direct current offset test volumes to red light wavelength
A Color Channel, generate multiple and different PPG test result signals corresponding to red light wavelength, surveyed from those different PPG
Test result signal, the anti-DC offset correction amount for estimating the PPG signals for measuring red light wavelength, action noise correct unit 1202
The DC offset correction amount can be used to execute an offset correction to the PPG signals of red light wavelength, similarly, offset determines single
Member 1201 can by gradually applying or adjusting the Color Channel corresponding to multiple and different direct current offset test volumes to green wavelength,
Multiple and different PPG test results signals corresponding to green wavelength are generated, the PPG test result signal different from those is counter to push away
A DC offset correction amount of the PPG signals of green wavelength is estimated, and using the DC offset correction amount come to green wavelength
PPG signals execute an offset correction.
Above-mentioned application or adjust different direct current offset test volumes when, can by one fix step size (step size) come
Those different direct current offset test volumes are generated, or a step size can be adjusted by dynamic and surveyed to generate those different direct current offsets
Examination amount.Wherein, for fixing step size, adjusting range can for example be selected red caused by optical sensing system itself
The half of positive value of the half of the negative value of the preliminary direct current offset R_DC of optical wavelength to the preliminary direct current offset, that is,-R_DC/2
To R_DC/2, and using the twentieth size (that is, R_DC/20) of the adjusting range as fixed step, it is inclined to carry out direct current
The testing and debugging of shifting.Similarly, for the testing and debugging of green wavelength and blue light wavelength, aforesaid operations are also such as same as, it is no longer superfluous
It states.It should be noted that above-mentioned example explanation is only that reader is enabled more to be readily appreciated that embodiments of the present invention, it is not the limitation of the present invention.
In addition, can also count and analyze the Fourier transform function of multiple test result signals of one or more Color Channels
Frequency spectrum, Fourier transform function (such as Fast Fourier Transform FFT) can be used by each survey in offset determining means 1201
Test result signal is converted from time-domain to the frequency spectrum of frequency domain, according to a ceiling capacity of one of each test result signal frequency spectrum
Value and a ratio, generate a certain threshold, and according in each test result signal be more than the certain threshold it
One local extremum number of radio-frequency component, from multiple and different DC offsets corresponding to multiple and different test result signals
In, the anti-DC offset correction amount for estimating the PPG signals for measuring a particular light ray wavelength.When a fc-specific test FC result signal it
When one local extremum number is less than a certain number, offset determining means 1201 is chosen corresponding to the fc-specific test FC result signal
One of DC offset as the DC offset correction amount.
With reference to Fig. 2, the depicted electronic device for Fig. 1 carries out the FFT of an example of the test result signal of offset correction
The time frequency analysis schematic diagram of frequency spectrum.As shown in Fig. 2, X-axis represents the time, the energy intensity numerical value that Y axis represents FFT spectrum (represents
The signal energy intensity of certain frequency), horizontal dotted line represents certain threshold TH, and 1201 system of offset determining means, which first detects, to be obtained
The maximum energy value of the frequency spectrum of the test result signal is located at A points, using corresponding to A points energy value and a ratio (such as
1/4) TH is calculated, that is, the numerical value of TH removes 4 for the energy value corresponding to A points, then compares in the test result signal
One local extremum number of the radio-frequency component more than certain threshold TH may detect that local extremum number is 2, if the spy
Determine number and be set as 4, then the local extremum number of this test result signal is 2 and is less than certain number 4, therefore offset
Determining means 1201 can choose the DC offset corresponding to the test result signal that this makes use as the direct current offset
Correcting value.
With reference to Fig. 3, the depicted electronic device for Fig. 1 carries out another example of the test result signal of offset correction
The time frequency analysis schematic diagram of FFT spectrum.As shown in figure 3, X-axis represents the time, Y-axis represents the energy intensity numerical value of FFT spectrum
(the signal energy intensity for representing certain frequency), horizontal dotted line represent certain threshold TH, and 1201 system of offset determining means is first detectd
The maximum energy value for measuring the frequency spectrum of the test result signal is located at A ' points, using the corresponding energy value of A ' points and a ratio
Value (such as 1/4) and TH is calculated, that is, the numerical value of TH is that the corresponding energy value of A ' points removes 4, then compares the test knot
It is more than a local extremum number of the radio-frequency component of certain threshold TH in fruit signal, may detect that local extremum number is 8,
If the certain number is set as 4, the local extremum number of this test result signal is 8 and is more than certain number 4, because
This offset determining means 1201 will not choose the DC offset conduct corresponding to the test result signal that this makes use
The DC offset correction amount.
When a test result signal is converted from time-domain to frequency domain, the maximum energy value of the frequency spectrum of the test result signal
Frequency typically represent estimated physiological characteristic (i.e. palmic rate), no direct current offset interference ideal situation, this is most
Big energy value can be several times as much as the energy intensity value of other radio-frequency components, therefore use a ratio (such as 1/4) to distinguish this
The frequency of maximum energy value and the radio-frequency component of other energy intensities, and the setting of a certain number is actually also used, sentence
Direct current offset suffered by disconnected current test result signal is affected or smaller, such as the test result news shown in Fig. 2
Number, the local extremum number for being more than the radio-frequency component of certain threshold TH is 2, and offset determining means 1201 can determine that
The test result signal influenced by direct current offset it is smaller, if without other test result signals better than the test result interrogate
Number, then the corresponding DC offset of the test result signal can be selected as DC offset correction amount.Conversely, as shown in Figure 3
Test result signal, be more than the radio-frequency component of certain threshold TH local extremum number be 8, be more than the certain number
4, offset determining means 1201 can determine that the test result signal is affected by direct current offset, not select the test
As a result the corresponding DC offset of signal is as DC offset correction amount.Accordingly, offset determining means 1201 can determine
The quality of different test result signals come determine it is last used by DC offset correction amount.
It should be noted that when offset determining means 1201 carries out spectrum analysis, can only take in normal heartbeat frequency range
Energy intensity value it is for statistical analysis, such as Fig. 2 shown in frequency f1 to frequency f2 range.So this be not this case limit
System.
When offset determining means 1201 finally determines the correcting value of direct current offset, action noise corrects unit 1202
The direct current offset for then reducing or eliminating the PPG signals of one of corresponding wavelength of light come adjustment according to the correcting value, then gives birth to
Reason feature estimation unit 1203 is then according to neither estimating palmic rate, implementation with the PPG signals after the correction of wavelength of light
On, physiological characteristic estimates unit 1203 can be to neither turning with the PPG signals after the correction of wavelength of light from chrominance space is carried out
It changes, takes the correspondence chromatic value of two PPG signals, calculate a ratio of two groups of chrominance signals, then the ratio is converted to frequency domain
Judge palmic rate.For example, in an embodiment, if offset determining means 1201 determines red light wavelength and green wavelength respectively
Color channel physiological characteristic measure signal two offset correction amounts, then action noise correction unit 1202 use red light wavelength
The direct current offset of the PPG signals of red light wavelength, green wavelength is reduced or eliminated with the two offset corrections amount of green wavelength, and
Two PPG signals of red light wavelength, green wavelength after correction are then carried out chrominance space by physiological characteristic estimation unit 1203 respectively
Conversion takes the correspondence chromatic value of two PPG signals, calculates the ratio of two groups of chrominance signals, then the ratio is converted to frequency domain
Judge palmic rate.Aforesaid operations are also suitable for the PPG signals of blue light wavelength and green wavelength.
In order to make reader be apparent to the operation of the present invention, separately by the fortune of the processing circuit 120 of above-mentioned electronic device 100
It is described in Fig. 4 as process step, to contribute to reader to understand.If generally can reach identical as a result, not needing to be certain
It is carried out according to the step sequence in the flow shown in Fig. 4, and the step shown in Fig. 4 is not necessarily intended to be carried out continuously, that is, other are walked
Suddenly it also can be inserted into wherein：
Step 405：Start；
Step 410：Electronic device 100 is enabled to be in the environment of masking ambient light；
Step 415：The DC offset for measuring PPG signals under total darkness environment, as preliminary offset correction amount；
Step 420A：Gradually apply respectively and adjusts multiple and different direct current offset test volumes to neither with corresponding to wavelength of light
Color Channel, generate corresponding to different wavelength of light multiple and different PPG test result signals；
Step 425A：It is counter to estimate the DC offset correction for surveying different wavelength of light from different PPG test result signals
Amount；
Step 420B：According to the maximum energy value and a ratio of one of every PPG test result signals frequency spectrum, generate
One certain threshold；
Step 425B：Compare a local pole of the radio-frequency component in each test result signal more than the certain threshold
It is worth number, counter estimate measures a DC offset correction amount；
Step 430：Determine that, neither with the final offset correction amount of wavelength of light, correction is neither interrogated with the PPG of wavelength of light
Number；
Step 435：To neither with the PPG signals after the correction of wavelength of light from color space convert is carried out, taking two PPG
The correspondence chrominance signal value of signal；And
Step 440：A ratio of two groups of chrominance signals is calculated, then the ratio is converted to frequency domain to judge user's heartbeat
Frequency.
Furthermore in implementation, above-mentioned PPG signals include a matrix value, when calculating offset correction amount, can respectively by with
Single pixel cell size is that unit (pixel-based) is calculated the addition of sensor pixel value, by with single pixel column
Size is unit (column-based) by value progress phase Calais's calculating of multiple pixel columns or by with whole picture frame size
The value of entire picture frame is subjected to phase Calais's calculating for unit (frame-based), similar implementation variation such as meets this hair
Bright spirit.
The foregoing is merely the preferred embodiments of the invention, all equivalent changes done according to scope of the present invention patent with
Modification, should all belong to the covering scope of the present invention.
Claims (18)
1. a kind of electronic device measuring physiological characteristic, including：
An at least luminescence unit, to emit an at least light, which includes to correspond to one of first wave length the
One light composition and corresponding to the second light composition of an one of second wave length, which is different from the second wave length；
One image sensing circuit, can be in response in an at least light, and sensing generates an at least physiology characteristic quantity and surveys signal；And
One processing circuit is coupled to the image sensing circuit, and one is carried out partially to survey signal for an at least physiology characteristic quantity
Shift correction operates, and measures signal to generate the physiological characteristic after at least one correction, and special to the physiology after at least one correction
Sign measures signal and is calculated, and a physiology characteristic results are obtained with estimation.
2. the electronic device as described in claim the 1st, the wherein processing circuit estimate an at least light composition it
An at least DC offset correction amount, and according to an at least DC offset correction amount, signal is surveyed to an at least physiology characteristic quantity
Carry out a direct current offset correction operation.
3. it includes pair that the electronic device as described in claim the 2nd, wherein an at least physiology characteristic quantity, which survey signal,
Signal should be measured in the first physiological characteristic of one of the first light composition and correspond to the second physiology of one of the second light composition
Feature measures signal, which includes：
One determining means corresponds to estimate to correspond to the first DC offset correction amount of one of the first light composition and estimation
In the second DC offset correction amount of one of the second light composition；
One correction unit, to respectively according to the first DC offset correction amount and the second DC offset correction amount, to this
One physiology characteristic quantity surveys signal and second physiological characteristic measures signal and executes DC offset correction operation, after generating the correction
The first physiological characteristic measure the second physiological characteristic after signal and the correction and measure signal；And
One estimation unit, to measure the second physiological characteristic after signal and the correction according to the first physiological characteristic after the correction
Signal is measured, calculating estimates the physiological characteristic result.
4. the electronic device as described in claim the 3rd, wherein the estimation unit carry out a color space convert, respectively
The first physiological characteristic after the correction is measured to the second physiological characteristic measurement signal after signal and the correction and is converted to one first
Chrominance signal and one second chrominance signal, and one of first chrominance signal and second chrominance signal ratio are calculated, and will
The ratio is converted to a frequency domain, the physiological characteristic result is calculated.
5. the electronic device as described in claim the 2nd, the wherein processing circuit measure one when lying in masking environment light
Direct current offset one of uses preliminary DC offset correction amount as the offset correction operation.
6. the electronic device as described in claim the 5th, wherein the processing circuit system are multiple and different by gradually applying
An at least Color Channel corresponding to direct current offset test volume to this at least a light composition, to generate multiple and different test knots
Fruit signal, the test result signal different from those, counter estimate survey the DC offset correction amount, use the DC offset correction amount
The offset correction operation is executed to survey signal to an at least physiology characteristic quantity.
7. the electronic device as described in claim the 6th, the wherein processing circuit gradually apply multiple different direct currents
It is that those different direct current offset test volumes are generated by a fixation step size when offset test volume extremely should be to Color Channel,
Or a step size is adjusted to generate those different direct current offset test volumes by dynamic.
8. the electronic device as described in claim the 6th, wherein the processing circuit system is according to each test result signal
One of frequency spectrum a maximum energy value and a ratio, a certain threshold is generated, and according to each test result signal
In be more than the certain threshold radio-frequency component a local extremum number, corresponding to the test result signal different from those it
In those different DC offsets, counter estimate measures the DC offset correction amount, wherein when one of fc-specific test FC result signal
When local extremum number is less than a certain number, it is inclined which chooses one of the corresponding direct current of the fc-specific test FC result signal
Shifting amount is as the DC offset correction amount.
9. it includes a number that the electronic device as described in claim the 1st, wherein an at least physiology characteristic quantity, which survey signal,
Value matrix, it is unit which, which can use with pixel unit (pixel), take pixel column (column) as unit or with
Picture frame (frame) calculates the numerical value of the numerical matrix for unit, one of is used with calculating the offset correction operation
Offset correction amount.
10. a kind of method measuring physiological characteristic, including：
Using an at least luminescence unit, to emit an at least light, an at least light include correspond to a first wave length it
One first light composition and corresponding to the second light composition of an one of second wave length, which is different from the second wave length；
Using an image sensing circuit, in response in an at least light, sensing generates an at least physiology characteristic quantity and surveys signal；
Signal is surveyed for an at least physiology characteristic quantity and carries out an offset correction operation, it is special to generate the physiology after at least one correction
Sign measures signal；And
Signal is measured to the physiological characteristic after at least one correction to calculate, and a physiology characteristic results are obtained with estimation.
11. the method as described in claim the 10th, wherein the step for carrying out the offset correction operation includes：Estimation should
An at least DC offset correction amount for an at least light composition, and according to an at least DC offset correction amount, to this at least one
Physiological characteristic measures signal and carries out a direct current offset correction operation.
12. it includes correspondence that the method as described in claim the 11st, wherein an at least physiology characteristic quantity, which survey signal,
Signal is measured in the first physiological characteristic of one of the first light composition and corresponding to one of the second light composition the second physiology spy
Sign measures signal, and this method additionally comprises：
Estimation corresponds to the first DC offset correction amount of one of the first light composition and estimation correspond to second light at
The second DC offset correction amount of one of part；
Respectively according to the first DC offset correction amount and the second DC offset correction amount, which is measured and is interrogated
Number and second physiological characteristic measure signal execute the DC offset correction operation, generate the first physiological characteristic amount after the correction
The second physiological characteristic surveyed after signal and the correction measures signal；And
The second physiological characteristic after signal and the correction is measured according to the first physiological characteristic after the correction and measures signal, and calculating is estimated
Measure the physiological characteristic result.
13. the method as described in claim the 12nd, wherein calculating the step of estimating the physiological characteristic result and including：
A color space convert is carried out, the first physiological characteristic after the correction is measured to the second life after signal and the correction respectively
Reason feature measures signal and is converted to one first chrominance signal and one second chrominance signal；
Calculate one of first chrominance signal and second chrominance signal ratio；And
The ratio is converted to a frequency domain, the physiological characteristic result is calculated.
14. the method as described in claim the 11st, wherein estimating an at least direct current for an at least light composition
The step of offset correction amount includes：A direct current offset is measured when covering environment light one of to use as the offset correction operation
Preliminary DC offset correction amount.
15. the method as described in claim the 14th, wherein estimating an at least direct current for an at least light composition
The step of offset correction amount additionally comprises：
Led to by gradually applying at least color corresponding to multiple and different direct current offset test volumes to this at least a light composition
Road, to generate multiple and different test result signals；And
The test result signal different from those, counter estimate survey the DC offset correction amount.
16. the method as described in claim the 15th, wherein gradually applying the step of multiple and different direct current offset test volumes
Suddenly include：Those different direct current offset test volumes are generated by a fixation step size, or big by dynamic one step of adjustment
It is small to generate those different direct current offset test volumes.
17. the method as described in claim the 15th, wherein from those different test result signals it is counter estimate survey should
The step of DC offset correction amount includes；
According to the maximum energy value and a ratio of one of each test result signal frequency spectrum, a certain threshold is generated；With
And
According to a local extremum number of the radio-frequency component in each test result signal being more than the certain threshold, from those
In those corresponding different DC offsets of different test result signals, counter estimate measures the DC offset correction amount；
It is to choose the fc-specific test FC wherein when one of fc-specific test FC result signal local extremum number is less than a certain number
As a result one of the corresponding DC offset of signal is as the DC offset correction amount.
18. it includes a numerical value that the method as described in claim the 10th, wherein an at least physiology characteristic quantity, which survey signal,
Matrix, the processing circuit system can using as unit of pixel unit, as unit of pixel column or as unit of picture frame come to this
The numerical value of numerical matrix is calculated, and offset correction amount one of is used to calculate the offset correction operation.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110934324.2A CN113558596B (en) | 2017-04-27 | 2017-04-27 | Method for measuring physiological characteristics |
| CN201710290128.XA CN108784675B (en) | 2017-04-27 | 2017-04-27 | Electronic device and method capable of eliminating or avoiding offset interference and accurately measuring physiological characteristics |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710290128.XA CN108784675B (en) | 2017-04-27 | 2017-04-27 | Electronic device and method capable of eliminating or avoiding offset interference and accurately measuring physiological characteristics |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110934324.2A Division CN113558596B (en) | 2017-04-27 | 2017-04-27 | Method for measuring physiological characteristics |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108784675A true CN108784675A (en) | 2018-11-13 |
| CN108784675B CN108784675B (en) | 2021-09-03 |
Family
ID=64069807
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710290128.XA Active CN108784675B (en) | 2017-04-27 | 2017-04-27 | Electronic device and method capable of eliminating or avoiding offset interference and accurately measuring physiological characteristics |
| CN202110934324.2A Active CN113558596B (en) | 2017-04-27 | 2017-04-27 | Method for measuring physiological characteristics |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110934324.2A Active CN113558596B (en) | 2017-04-27 | 2017-04-27 | Method for measuring physiological characteristics |
Country Status (1)
| Country | Link |
|---|---|
| CN (2) | CN108784675B (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5355882A (en) * | 1992-07-09 | 1994-10-18 | Nihon Kohden Corporation | Pulse oximeter |
| WO2003068060A1 (en) * | 2002-02-15 | 2003-08-21 | Datex-Ohmeda, Inc. | Compensation of human variability in pulse oximetry |
| CN1692874A (en) * | 2004-05-08 | 2005-11-09 | 香港中文大学 | Finger ring type physiological information monitoring device |
| US20070027375A1 (en) * | 2002-06-20 | 2007-02-01 | Melker Richard J | Optimized gas supply using photoplethysmography |
| CN101730503A (en) * | 2007-03-15 | 2010-06-09 | 帝国创新有限公司 | Heart rate measurement |
| US20120265080A1 (en) * | 2011-04-15 | 2012-10-18 | Xiong Yu | Non-contact sensing of physiological signals |
| CN105249939A (en) * | 2014-07-17 | 2016-01-20 | 原相科技股份有限公司 | Physiology detection module capable of denoising and physiology detection method thereof |
| CN105491943A (en) * | 2014-06-30 | 2016-04-13 | 皇家飞利浦有限公司 | Photoplethysmography sensor apparatus and method |
| CN105832289A (en) * | 2015-01-30 | 2016-08-10 | 三星电子株式会社 | Method and equipment using Hilbert transform to estimate biophysiological rates |
| CN106456030A (en) * | 2014-06-13 | 2017-02-22 | 日东电工株式会社 | Device and method for removing artifacts in physiological measurements |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH10179541A (en) * | 1996-12-26 | 1998-07-07 | Toshiba Corp | Mr signal processing circuit and processing method for mr signal |
| JP3789024B2 (en) * | 1997-06-11 | 2006-06-21 | 株式会社東芝 | SQUID magnetometer |
| US7171251B2 (en) * | 2000-02-01 | 2007-01-30 | Spo Medical Equipment Ltd. | Physiological stress detector device and system |
| EP1611847A1 (en) * | 2004-06-28 | 2006-01-04 | Datex-Ohmeda, Inc. | Validating pulse oximetry signals in the potential presence of artifact |
| US8761853B2 (en) * | 2011-01-20 | 2014-06-24 | Nitto Denko Corporation | Devices and methods for non-invasive optical physiological measurements |
| JP6270287B2 (en) * | 2012-11-23 | 2018-01-31 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | Device and method for extracting physiological information |
| EP2896359B1 (en) * | 2014-01-07 | 2016-12-14 | Samsung Electronics Co., Ltd | Method and system for measuring heart rate in electronic device using photoplethysmography |
| JP6629242B2 (en) * | 2014-05-28 | 2020-01-15 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | Motion artifact reduction using multi-channel PPG signals |
| JP6615176B2 (en) * | 2014-07-22 | 2019-12-04 | コーニンクレッカ フィリップス エヌ ヴェ | Non-interfering skin tissue hydration measuring device and related method |
| TWI558375B (en) * | 2014-09-18 | 2016-11-21 | 義明科技股份有限公司 | Ppg signal processing device and method thereof |
| US20160174887A1 (en) * | 2014-12-18 | 2016-06-23 | Koninklijke Philips N.V. | System and method for extracting physiological information and medical instrument for use in the system |
| CN105796051B (en) * | 2014-12-31 | 2019-11-19 | 原相科技股份有限公司 | Three-dimensional physiology-detecting system and its operating method |
-
2017
- 2017-04-27 CN CN201710290128.XA patent/CN108784675B/en active Active
- 2017-04-27 CN CN202110934324.2A patent/CN113558596B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5355882A (en) * | 1992-07-09 | 1994-10-18 | Nihon Kohden Corporation | Pulse oximeter |
| WO2003068060A1 (en) * | 2002-02-15 | 2003-08-21 | Datex-Ohmeda, Inc. | Compensation of human variability in pulse oximetry |
| US20070027375A1 (en) * | 2002-06-20 | 2007-02-01 | Melker Richard J | Optimized gas supply using photoplethysmography |
| CN1692874A (en) * | 2004-05-08 | 2005-11-09 | 香港中文大学 | Finger ring type physiological information monitoring device |
| CN101730503A (en) * | 2007-03-15 | 2010-06-09 | 帝国创新有限公司 | Heart rate measurement |
| US20120265080A1 (en) * | 2011-04-15 | 2012-10-18 | Xiong Yu | Non-contact sensing of physiological signals |
| CN106456030A (en) * | 2014-06-13 | 2017-02-22 | 日东电工株式会社 | Device and method for removing artifacts in physiological measurements |
| CN105491943A (en) * | 2014-06-30 | 2016-04-13 | 皇家飞利浦有限公司 | Photoplethysmography sensor apparatus and method |
| CN105249939A (en) * | 2014-07-17 | 2016-01-20 | 原相科技股份有限公司 | Physiology detection module capable of denoising and physiology detection method thereof |
| CN105832289A (en) * | 2015-01-30 | 2016-08-10 | 三星电子株式会社 | Method and equipment using Hilbert transform to estimate biophysiological rates |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113558596B (en) | 2024-08-09 |
| CN113558596A (en) | 2021-10-29 |
| CN108784675B (en) | 2021-09-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20160338603A1 (en) | Signal processing device, signal processing method, and computer-readable recording medium | |
| US8498689B2 (en) | Real time motion information capture in an MRI environment | |
| CN113647921B (en) | Blood pressure measuring method, system, equipment and storage medium | |
| US20120283535A1 (en) | Method and system for pulse measurement | |
| Vinciguerra et al. | Progresses towards a processing pipeline in photoplethysmogram (PPG) based on SiPMs | |
| CN112806972A (en) | PPG test circuit and method, and wearable electronic device | |
| JP2014124454A (en) | Method and apparatus for measuring blood component, and medical equipment | |
| US4273442A (en) | Spectral photometer for medical applications | |
| CN114569101A (en) | Non-contact heart rate detection method and device and electronic equipment | |
| CN108784675A (en) | Can eliminate or avoid offset interference and can accurate measurement physiological characteristic electronic device and method | |
| CN103393425A (en) | Real-time signal quality evaluation algorithm for near-infrared brain functional imaging system | |
| US11025843B2 (en) | Device and method for the continuous and non-invasive determination of physiological parameters of a test subject | |
| US20230078228A1 (en) | Blood Oxygen Saturation Measurement Method And Apparatus | |
| JP2000155091A (en) | Subcutaneous fat measuring device, near infrared biological measuring device and near infrared biological measuring method | |
| Raposo et al. | Camera-based Photoplethysmography (cbPPG) using smartphone rear and frontal cameras: An experimental study | |
| CN105796051B (en) | Three-dimensional physiology-detecting system and its operating method | |
| JP2021023490A (en) | Biological information detection device | |
| US12070309B2 (en) | Method, computer programming product and electronic device for calculating blood oxygen saturation | |
| CN209951236U (en) | Optical image pulse measuring system | |
| Karlen et al. | Recognition of correct finger placement for photoplethysmographic imaging | |
| WO2001043624A2 (en) | Method for de-weighting motion-contaminated data in photoplethysmographic analyte measurements | |
| CN108814619B (en) | System and device for measuring high-precision blood oxygen saturation | |
| CN111166314A (en) | Heart rate measuring device and method based on combination of multi-wavelength light source and accelerometer | |
| JP4946265B2 (en) | Color measuring device, color measuring method and color error measuring method | |
| CN109770885A (en) | A kind of examing heartbeat fastly method based on preview frame |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |