CN111444797B - Non-contact heart rate detection method - Google Patents

Non-contact heart rate detection method Download PDF

Info

Publication number
CN111444797B
CN111444797B CN202010179914.4A CN202010179914A CN111444797B CN 111444797 B CN111444797 B CN 111444797B CN 202010179914 A CN202010179914 A CN 202010179914A CN 111444797 B CN111444797 B CN 111444797B
Authority
CN
China
Prior art keywords
signal
skin color
heart rate
camera
signals
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010179914.4A
Other languages
Chinese (zh)
Other versions
CN111444797A (en
Inventor
岳洪伟
王洪涛
邓辅秦
许弢
李俊华
金迎迎
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuyi University
Original Assignee
Wuyi University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuyi University filed Critical Wuyi University
Priority to CN202010179914.4A priority Critical patent/CN111444797B/en
Publication of CN111444797A publication Critical patent/CN111444797A/en
Application granted granted Critical
Publication of CN111444797B publication Critical patent/CN111444797B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/22Source localisation; Inverse modelling
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/024Detecting, measuring or recording pulse rate or heart rate
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/21Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • G06F18/2134Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on separation criteria, e.g. independent component analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/26Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
    • G06V10/267Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion by performing operations on regions, e.g. growing, shrinking or watersheds
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/56Extraction of image or video features relating to colour
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/12Classification; Matching
    • G06F2218/14Classification; Matching by matching peak patterns

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Multimedia (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Cardiology (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Signal Processing (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Physiology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

According to the non-contact heart rate detection method, the chrominance characteristic S is firstly calculated for the signal acquired by the camera 1 (t) and S 2 (t) by the colorimetric characteristics S 1 (t) and S 2 (t) preprocessing and extracting each path of signal to obtain pulse signals h (t) and h 1 (t) and h 2 (t), and then the extracted pulse signal h 1 (t) and h 2 (t) adopting FastICA algorithm to further separate blind sources, and utilizing pulse signals h (t) and h 1 (t) and h 2 (t) determining the kurtosis value of the power spectrum as a blind source separation result, and in addition, when the kurtosis value of the power spectrum of the signal is close, selecting the signal by using the kurtosis value easily and wrongly, and using the signal-to-noise ratio as an index can be helpful for completing the signal selection; therefore, the invention can effectively process noise interference caused by illumination change and involuntary movement of a subject, obtain a signal with lower noise and improve the anti-interference performance and the measurement precision of detection.

Description

Non-contact heart rate detection method
Technical Field
The application relates to the technical field of signal processing, in particular to a non-contact heart rate detection method.
Background
The heart rate refers to the number of beats per minute of the heart, and modern medicine finds that the slower the heart rate is, the longer the life of the animal is; heart rate is an important index of health, currently, the most used clinical devices are patch type heart rate detectors, which extract heart rate from human body electrical signals, although the patch type heart rate detectors have high precision, the patch type heart rate detectors must be in contact with human body when in use, which causes great inconvenience to certain people, such as heart rate monitoring of newborn babies; therefore, the video-based non-contact heart rate detection has a wide application prospect in the field of medical health, however, the video-based non-contact heart rate detection is easily interfered by noise, such as illumination change, involuntary movement of a subject, and the like, which makes it difficult to effectively extract signals, and may cause the accuracy and stability of heart rate measurement to be reduced, and therefore, a heart rate detection method capable of resisting interference and improving the measurement accuracy is required.
Disclosure of Invention
The purpose of this application lies in solving one of the technical problem that exists among the prior art at least, for this reason, this application provides non-contact heart rate detection method, improves non-contact heart rate detection's interference immunity and measurement accuracy.
The technical scheme adopted by the invention for solving the technical problems is as follows:
according to an embodiment of the invention, there is provided a non-contact heart rate detection method, comprising the steps of:
acquiring a video signal, and selecting a skin color region ROI1 from a video frame image by adopting a skin color detection algorithm;
step two, dividing the skin color area ROI1 into at least two different skin color areas ROI2 and ROI3 to generate corresponding RGB observation signals;
thirdly, filtering and time normalization preprocessing are carried out on the RGB observation signals, and primary color separation is carried out on three color channels of the RGB observation signals to obtain three-channel data of green G (t), blue B (t) and red R (t);
step four, calculating the chromaticity characteristics S of the skin color area ROI1, the skin color area ROI2 and the skin color area ROI3 1 (t) and S 2 (t),S 1 (t)=G(t)-B(t),S 2 (t)=G(t)+B(t)-2R(t);
Step five, calculating the skin color region ROI1, the skin color region ROI2 and the skin color region ROI3 to be respectively alignedCorresponding pulse signals h (t), h 1 (t) and h 2 (t),h(t)=S 1 (t)+α·S 2 (t);
Step six, pulse signals h are corrected 1 (t) and pulse signal h 2 (t) further performing blind source separation to obtain pulse signal Y 1 (t) and Y 2 (t);
Step seven, calculating pulse signals h (t) and Y 1 (t) and Y 2 (t) kurtosis and signal-to-noise ratio of the power spectrum;
step eight, comparing the pulse signals h (t) and Y 1 (t) and Y 2 (t) kurtosis value and SNR of power spectrum, when pulse signals h (t) and Y 1 (t) or Y 2 (t) when the kurtosis value and the signal-to-noise ratio of a group of component signals are simultaneously maximum, the component is the acquired BVP signal; otherwise, the pulse signal h (t) is taken as the BVP signal.
The non-contact heart rate detection method provided by the embodiment of the invention at least has the following beneficial effects: in the invention, the signal collected by the camera is firstly calculated to obtain the chrominance characteristic S 1 (t) and S 2 (t) by the colorimetric characteristics S 1 (t) and S 2 (t) preprocessing and extracting each path of signal to obtain pulse signals h (t) and h 1 (t) and h 2 (t), and then the extracted pulse signal h 1 (t) and h 2 (t) adopting a FastICA algorithm to further carry out blind source separation, meanwhile, because the pulse signals have obvious quasi-periodicity, the power spectrum of the pulse signals can have very obvious spectral peaks in the heart rate frequency range, the corresponding kurtosis value is usually larger, the other signal component does not have obvious spectral peaks in the heart rate frequency range, and the corresponding kurtosis value of the power spectrum is relatively lower, so the pulse signals h (t) and h are utilized 1 (t) and h 2 (t) the kurtosis value of the power spectrum is used for judging the blind source separation result, in addition, when the kurtosis of the power spectrum of the signal is close, the signal is selected by mistake easily only by using the kurtosis value, so that the signal-to-noise ratio is also used as an index to be helpful for completing the signal selection; therefore, the invention can effectively process noise interference caused by illumination change and involuntary movement of a subject, and obtain a signal with lower noiseHigh anti-interference performance and high measurement precision.
According to some embodiments of the invention, the skin color region is segmented using approximately halving.
According to some embodiments of the invention, the kurtosis value is calculated by:
Figure BDA0002412144910000031
wherein,
Figure BDA0002412144910000032
x (n) is an input signal; and N is the number of sampling points.
According to some embodiments of the invention, the snr is calculated by:
Figure BDA0002412144910000033
wherein S (f) is the frequency spectrum of the extracted pulse signal, (U) t (f) The passband of the spectral peak.
According to some embodiments of the invention, the passband bandwidth of the spectral peak is 0.4Hz.
According to some embodiments of the invention, wherein
Figure BDA0002412144910000034
According to some embodiments of the invention, when the color type of the video frame collected by the camera is an RGB type, the three color channels are directly subjected to primary color separation; when the color type of the video frame collected by the camera is not the RGB type, the video frame is converted into the RGB color type, and then the primary colors are separated.
According to some embodiments of the invention, the sampling frequency of the camera is 30Hz/s.
According to some embodiments of the invention, the camera is located no more than 1 meter away from the subject.
According to some embodiments of the present invention, the photographing region of the camera is skin at any region of the subject.
Drawings
The above and/or additional aspects and advantages of the present invention will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which:
FIG. 1 is a schematic flow chart of a method according to an embodiment of the present invention;
FIG. 2 is a spectrum diagram of a pulse signal h (t) according to an embodiment of the present invention;
FIG. 3 shows a pulse signal h according to an embodiment of the present invention 1 (t) and h 2 (t) mixed spectrogram;
FIG. 4 shows a pulse signal Y according to an embodiment of the present invention 1 (t) a spectrogram;
FIG. 5 shows a pulse signal Y according to an embodiment of the present invention 2 (t) spectrogram.
Detailed Description
Reference will now be made in detail to the present preferred embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to like elements throughout.
In the description of the present invention, unless explicitly defined otherwise, terms such as calculating, obtaining, sampling, corresponding and the like should be understood in a broad sense, and those skilled in the art can reasonably determine the specific meanings of the terms in the present invention by combining the specific contents of the technical solutions.
As shown in fig. 1, a non-contact heart rate detection method provided in an embodiment of the present application includes the following steps:
a method of non-contact heart rate detection is provided, comprising the steps of:
acquiring a video signal, and selecting a skin color region ROI1 from a video frame image by adopting a skin color detection algorithm;
dividing a skin color region ROI1 into at least two different skin color regions ROI2 and ROI3 to generate corresponding RGB observation signals;
thirdly, filtering and time normalization preprocessing are carried out on the RGB observation signals, and primary color separation is carried out on three color channels of the RGB observation signals to obtain three-channel data of green G (t), blue B (t) and red R (t);
step four, calculating the chromaticity characteristics S of the skin color region ROI1, the skin color region ROI2 and the skin color region ROI3 1 (t) and S 2 (t),S 1 (t)=G(t)-B(t),S 2 (t)=G(t)+B(t)-2R(t);
Step five, calculating pulse signals h (t), h corresponding to the skin color region ROI1, the skin color region ROI2 and the skin color region ROI3 respectively 1 (t) and h 2 (t),h(t)=S 1 (t)+α·S 2 (t) in which
Figure BDA0002412144910000051
Step six, pulse signals h are corrected 1 (t) and pulse signal h 2 (t) further performing blind source separation to obtain pulse signal Y 1 (t) and Y 2 (t);
Step seven, calculating pulse signals h (t) and Y 1 (t) and Y 2 (t) kurtosis value and signal-to-noise ratio of the power spectrum, wherein the kurtosis value is calculated by the formula:
Figure BDA0002412144910000052
wherein,
Figure BDA0002412144910000053
x (n) is an input signal; n is the number of sampling points; the calculation formula of the signal-to-noise ratio is as follows:
Figure BDA0002412144910000054
wherein S (f) is the frequency spectrum of the extracted pulse signal, (U) t (f) A passband at the peak of the spectrum; the kurtosis factor is used for representing the smoothness degree of a waveform and describing the distribution of variables, the kurtosis of normal distribution is equal to 3, the distribution curve is relatively flat when the kurtosis is less than 3, and the distribution curve is relatively steep when the kurtosis is more than 3;
step eight, comparing the pulse signals h (t) and Y 1 (t) and Y 2 (t) kurtosis value and signal-to-noise ratio of power spectrum, when pulse signals h (t) and Y 1 (t) or Y 2 (t) one group ofWhen the kurtosis value and the signal-to-noise ratio of the component signal are simultaneously maximum, the component is the acquired BVP signal; otherwise, the pulse signal h (t) is taken as the BVP signal. In the invention, the signal collected by the camera is firstly calculated to obtain the chrominance characteristic S 1 (t) and S 2 (t) by the colorimetric characteristics S 1 (t) and S 2 (t) preprocessing and extracting each path of signal to obtain pulse signals h (t) and h 1 (t) and h 2 (t), and then the extracted pulse signal h 1 (t) and h 2 (t) adopting a FastICA algorithm to further carry out blind source separation, meanwhile, because the pulse signals have obvious quasi-periodicity, the power spectrum of the pulse signals can have very obvious spectral peaks in the heart rate frequency range, the corresponding kurtosis value is usually larger, the other signal component does not have obvious spectral peaks in the heart rate frequency range, and the corresponding kurtosis value of the power spectrum is relatively lower, so the pulse signals h (t) and h are utilized 1 (t) and h 2 (t) the kurtosis value of the power spectrum is used for judging the blind source separation result, in addition, when the kurtosis of the power spectrum of the signal is close, the signal is selected by mistake easily only by using the kurtosis value, so that the signal-to-noise ratio is also used as an index to be helpful for completing the signal selection; therefore, compared with the traditional bar code correction mode, the method firstly extracts the pulse signals and then carries out blind source separation; and the kurtosis value and the signal-to-noise ratio are used as indexes to select signals, so that noise interference caused by illumination change and involuntary movement of a subject can be effectively processed, signals with lower noise are obtained, and the anti-interference performance and the measurement accuracy of detection are improved.
In another embodiment of the present invention, the skin color area is divided into two regions by approximately bisecting the skin color area, and if the skin color area is divided into three regions, the skin color area is approximately halved.
As shown in fig. 2-5, in one embodiment of the present invention, the sampling frequency of the camera is set to 30Hz/s, and the distance from the camera to the tested object is 1 meter, but of course, other distances are also possible; shooting the skin of the tested object by a camera to obtain a video signal, and collecting the video signal from a video frame image by adopting a skin color detection algorithmAn interested skin color region ROI1, equally dividing the skin color region ROI1 into two different skin color regions ROI2 and ROI3 to generate corresponding RGB observation signals, carrying out filtering and time normalization pretreatment on the RGB observation signals, and carrying out primary color separation on three color channels of the RGB observation signals to obtain three-channel data of green G (t), blue B (t) and red R (t); calculating pulse signals h (t), h corresponding to the skin color region ROI1, the skin color region ROI2 and the skin color region ROI3 respectively 1 (t) and h 2 (t),h(t)=S 1 (t)+α·S 2 (t) in which
Figure BDA0002412144910000071
As shown in fig. 2, it is a spectrogram of a pulse signal h (t), and it can be seen from the spectrogram that the extracted pulse signal h (t) is severely interfered by noise, and a frequency value corresponding to a main peak of a frequency spectrum is not a pulse frequency measured this time; further, the pulse signals h are respectively corrected 1 (t) and h 2 (t) adopting a FastICA algorithm to further carry out blind source separation, and respectively obtaining independent component pulse signals Y of the two regions which are more accurate 1 (t) and Y 2 (t), as shown in FIG. 3, FIG. 3 shows the pulse signal h 1 (t) and h 2 (t) mixed spectrum, as shown in FIGS. 4-5, FIG. 4 is the independent component pulse signal Y 1 (t), FIG. 5 shows the pulse signal Y of the independent component 2 (t), the corresponding pulse wave signal Y can be seen from FIG. 4 1 (t) has a more pronounced pulse signal than the pulse wave signal h (t) of fig. 2. H (t) of FIG. 2 and Y of FIG. 4 1 (t) and Y of FIG. 5 2 The kurtosis of the power spectrum of (t) is 10.23,22.98 and 10.29, respectively, and the corresponding signal-to-noise ratios are 0.19,0.45 and 0.21, respectively. Thus, the independent component pulse signal Y 1 (t) has the highest kurtosis and signal-to-noise ratio, so that the pulse signal Y is an independent component 1 And (t) the corresponding time domain signal can be used as an effective pulse signal.
In another embodiment of the present invention, the video signal is captured by a camera, and specifically, a camera of a mobile phone can be used for shooting.
In another embodiment of the invention, when the color type of the video frame collected by the camera is an RGB type, the three color channels are directly subjected to primary color separation; when the color type of the video frame collected by the camera is not the RGB type, the video frame is converted into the RGB color type, and then the primary colors are separated.
In another embodiment of the present invention, the distance between the camera and the object to be photographed is not more than 1 meter, which is limited by the resolution and the lens focal length of a general camera, and the shooting distance is set to not more than 1 meter in this embodiment, however, when a camera with better resolution and lens focal length is selected, the distance between the camera and the object to be photographed can be more than 1 meter.
In another embodiment of the invention, the shooting area of the camera is the skin of any area of the object to be shot, when a skin area is selected for shooting, and when the limbs of the object to be tested do not move autonomously, the shooting area can leave the shooting visual field range, and the shooting area needs to be readjusted to recover the previous shooting area.
The above is only a preferred embodiment of the present invention, but the present invention is not limited to the above embodiments, and the technical effects of the present invention can be achieved by any means which are the same or similar, and all of which are within the scope of the present invention.

Claims (8)

1. The non-contact heart rate detection method is characterized by comprising the following steps of:
acquiring a video signal, and selecting a skin color region ROI1 from a video frame image by adopting a skin color detection algorithm;
step two, dividing the skin color area ROI1 into at least two different skin color areas ROI2 and ROI3 to generate corresponding RGB observation signals;
thirdly, filtering and time normalization preprocessing are carried out on the RGB observation signals, and primary color separation is carried out on three color channels of the RGB observation signals to obtain three-channel data of green G (t), blue B (t) and red R (t);
fourthly, calculating the chromaticity characteristics S of the skin color region ROI1, the skin color region ROI2 and the skin color region ROI3 1 (t) and S 2 (t),S 1 (t)=G(t)-B(t),S 2 (t)=G(t)+B(t)-2R(t);
Step five, calculating pulse signals h (t), h corresponding to the skin color region ROI1, the skin color region ROI2 and the skin color region ROI3 respectively 1 (t) and h 2 (t),h(t)=S 1 (t)+α·S 2 (t);
Step six, pulse signals h are corrected 1 (t) and pulse signal h 2 (t) further performing blind source separation to obtain pulse signal Y 1 (t) and Y 2 (t);
Step seven, calculating pulse signals h (t) and Y 1 (t) and Y 2 (t) kurtosis and signal-to-noise ratio of the power spectrum;
step eight, comparing the pulse signals h (t) and Y 1 (t) and Y 2 (t) kurtosis value and SNR of power spectrum, when pulse signals h (t) and Y 1 (t) or Y 2 (t) when the kurtosis value and the signal-to-noise ratio of a group of component signals are simultaneously maximum, the component is the acquired BVP signal; otherwise, taking the pulse signal h (t) as a BVP signal;
in step 7, the kurtosis value is calculated by the following formula:
Figure FDA0003998163160000011
wherein,
Figure FDA0003998163160000012
x (n) is an input signal; n is the number of sampling points;
in said step 5, wherein
Figure FDA0003998163160000021
2. The method for non-contact heart rate detection according to claim 1, wherein in the step 1, the skin color region is divided into regions by approximately bisection.
3. Contactless heart rate system of claim 1The detection method is characterized in that, in the step 7, the calculation formula of the signal-to-noise ratio is as follows:
Figure FDA0003998163160000022
wherein S (f) is the frequency spectrum of the extracted pulse signal, (U) t (f) The passband of the spectral peak.
4. The method of claim 3, wherein the passband bandwidth of the spectral peak is 0.4Hz.
5. The non-contact heart rate detection method according to claim 2, wherein when the color type of the video frame collected by the camera is an RGB type, the three color channels are directly subjected to primary color separation; when the color type of the video frame collected by the camera is not the RGB type, the video frame is converted into the RGB color type, and then the primary colors are separated.
6. The method of claim 5, wherein the sampling frequency of the camera is 30Hz/s.
7. The non-contact heart rate detection method according to claim 6, wherein the distance between the camera and the shot object is not more than 1 meter.
8. The non-contact heart rate detection method according to claim 7, wherein the photographing region of the camera is skin at any region of a subject to be photographed.
CN202010179914.4A 2020-03-16 2020-03-16 Non-contact heart rate detection method Active CN111444797B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010179914.4A CN111444797B (en) 2020-03-16 2020-03-16 Non-contact heart rate detection method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010179914.4A CN111444797B (en) 2020-03-16 2020-03-16 Non-contact heart rate detection method

Publications (2)

Publication Number Publication Date
CN111444797A CN111444797A (en) 2020-07-24
CN111444797B true CN111444797B (en) 2023-03-17

Family

ID=71654086

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010179914.4A Active CN111444797B (en) 2020-03-16 2020-03-16 Non-contact heart rate detection method

Country Status (1)

Country Link
CN (1) CN111444797B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112102949A (en) * 2020-08-31 2020-12-18 安徽大学 Facial blood flow pulse signal and facial organ motion feature extraction method
CN112863546A (en) * 2021-01-21 2021-05-28 安徽理工大学 Belt conveyor health analysis method based on audio characteristic decision
CN112869737B (en) * 2021-02-01 2022-10-28 浙江大学山东工业技术研究院 Non-contact human body blood oxygen saturation detection method
CN112965060A (en) * 2021-02-19 2021-06-15 加特兰微电子科技(上海)有限公司 Detection method and device for vital sign parameters and method for detecting physical sign points

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105962915A (en) * 2016-06-02 2016-09-28 安徽大学 Non-contact human body respiration rate and heart rate synchronous measurement method and system
CN109044322A (en) * 2018-08-29 2018-12-21 北京航空航天大学 A kind of contactless heart rate variability measurement method
CN109793506A (en) * 2019-01-18 2019-05-24 合肥工业大学 A kind of contactless radial artery Wave shape extracting method
CN110384491A (en) * 2019-08-21 2019-10-29 河南科技大学 A kind of heart rate detection method based on common camera

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105962915A (en) * 2016-06-02 2016-09-28 安徽大学 Non-contact human body respiration rate and heart rate synchronous measurement method and system
CN109044322A (en) * 2018-08-29 2018-12-21 北京航空航天大学 A kind of contactless heart rate variability measurement method
CN109793506A (en) * 2019-01-18 2019-05-24 合肥工业大学 A kind of contactless radial artery Wave shape extracting method
CN110384491A (en) * 2019-08-21 2019-10-29 河南科技大学 A kind of heart rate detection method based on common camera

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于单片机的智能SMS心率计设计;张建民 等;《现代电子技术》;20170331;第40卷(第6期);第1-2页 *

Also Published As

Publication number Publication date
CN111444797A (en) 2020-07-24

Similar Documents

Publication Publication Date Title
CN111444797B (en) Non-contact heart rate detection method
CN110269600B (en) Non-contact video heart rate detection method based on multivariate empirical mode decomposition and combined blind source separation
CN106778695B (en) Multi-person rapid heart rate detection method based on video
Subramaniam et al. Estimation of the Cardiac Pulse from Facial Video in Realistic Conditions.
CN109044322A (en) A kind of contactless heart rate variability measurement method
CN109820499B (en) High anti-interference heart rate detection method based on video, electronic equipment and storage medium
CN105147274A (en) Method for extracting heart rate from visible spectrum section face video signal
CN109350030B (en) System and method for processing human face video heart rate signal based on phase amplification
CN112233813A (en) Non-contact non-invasive heart rate and respiration measurement method and system based on PPG
CN111938622B (en) Heart rate detection method, device and system and readable storage medium
CN111127511B (en) Non-contact heart rate monitoring method
CN111281367A (en) Anti-interference non-contact heart rate detection method based on face video
CN113068015A (en) Endoscope image distortion correction system based on optical fiber probe
CN111839492A (en) Heart rate non-contact type measuring method based on face video sequence
CN115153473B (en) Non-contact heart rate detection method based on multivariate singular spectrum analysis
CN113591769B (en) Non-contact heart rate detection method based on photoplethysmography
KR102468654B1 (en) Method of heart rate estimation based on corrected image and apparatus thereof
CN112788200A (en) Method and device for determining frequency spectrum information, storage medium and electronic device
CN114569101A (en) Non-contact heart rate detection method and device and electronic equipment
WO2023165482A1 (en) Method and apparatus for heart rate detection
CN105046209A (en) Non-contact heart rate measurement method based on canonical correlation analysis
CN117064356A (en) Non-contact blood pressure monitoring device, method and system for enhancing robustness by utilizing spatial redundancy
Hu et al. Illumination robust heart-rate extraction from single-wavelength infrared camera using spatial-channel expansion
CN110321781B (en) Signal processing method and device for non-contact measurement
Wang et al. KLT algorithm for non-contact heart rate detection based on image photoplethysmography

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