CN105550502A - Data processing method based on health monitoring - Google Patents

Data processing method based on health monitoring Download PDF

Info

Publication number
CN105550502A
CN105550502A CN201510897273.5A CN201510897273A CN105550502A CN 105550502 A CN105550502 A CN 105550502A CN 201510897273 A CN201510897273 A CN 201510897273A CN 105550502 A CN105550502 A CN 105550502A
Authority
CN
China
Prior art keywords
data
value
fusion
physiological data
time
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
Application number
CN201510897273.5A
Other languages
Chinese (zh)
Other versions
CN105550502B (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.)
Beijing Qi Huang clinical research center
Original Assignee
Nanjing Post and Telecommunication 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 Nanjing Post and Telecommunication University filed Critical Nanjing Post and Telecommunication University
Priority to CN201510897273.5A priority Critical patent/CN105550502B/en
Publication of CN105550502A publication Critical patent/CN105550502A/en
Application granted granted Critical
Publication of CN105550502B publication Critical patent/CN105550502B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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/0205Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
    • A61B5/02055Simultaneously evaluating both cardiovascular condition and temperature
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment

Abstract

The invention discloses a data processing method based on health monitoring. According to the method, firstly, measurement accuracy is ensured through a multi-touch multi-sensing data collection method; secondly, multiple groups of physiological data are screened according to the Grubbs criterion so as to remove collected gross error data; and finally collected data are fused through combining a batch estimation algorithm and a least squares algorithm so as to reduce influences of different measurement precisions of different sensors on results. According to the method of the invention, the physiological information of a user can be obtained accurately so as to provide valid basis for real time detection of a health monitoring system.

Description

A kind of data processing method based on health monitoring
Technical field
The invention belongs to portable medical technical field, particularly relate to a kind of data processing method based on health monitoring.
Background technology
Wearable technology is the brand new ideas grown up recent years, and it is widely used in the aspect such as clinical monitoring, family health care, emergent rescue, special population monitoring, Psychological Assessment, athletic training.And wearable technology is exactly the clothing of daily physiologic information detection technique and the daily wearing of people, ornaments are merged mutually in the application of medical domain, thus is not affecting user's daily life and under ensureing comfortable prerequisite, the characteristic parameter of user detected.But, in the physiological parameter detected, due to various environmental factor, many rough error parameters can be mixed.And in image data process, the angle of measurement is single, can not ensure the accuracy of data.
In addition, in some wearable devices, adopt single data fusion method, the physiologic information gathered is processed, when sensor experiences failure, the data that gather of Timeliness coverage can not there is error, finally cause the physiologic information that accurately can not reflect user.
Summary of the invention
In order to overcome above-mentioned the deficiencies in the prior art, the present invention proposes a kind of data processing method based on health monitoring.First the method utilizes the method for the many sensings of multiconductor to ensure that the accuracy of measurement, then screens the physiological data gathered, and merges, and final acquisition one group accurately can reflect the data of user's physiologic information.On the one hand, the multiconductor many sensing measurements mode in the method, learning process has good theoretical performance and ensures, accuracy of measurement is high; On the other hand, it is the reference scheme that wearable watch-dog provides a kind of data processing, has good application prospect.
In order to solve the problems of the technologies described above, the technical solution used in the present invention is:
Based on a data processing method for health monitoring, comprise the steps:
Step 1: the M class physiological data being gathered person to be measured by sensor, uses L mindividual m+1 class physiological data sensor, gather the m+1 class physiological data of person to be measured, its detailed process is:
Step 1-1, gather for the 1st time of physiological data, each sensor is not measured in the same time respectively the N number of of same contact, obtains N group physiological data, often organizes data and is expressed as:
A 1 k m [ Q ] = [ a 1 k m 0 , a 1 k m 1 , ... a 1 k m i , ... , a 1 k m ( Q - 1 ) ] ;
Step 1-2, gather for the 2nd time of physiological data, each sensor is respectively at L mthe N number of of individual contact does not measure in the same time, obtains N group physiological data, often organizes data and is expressed as:
A 2 k m [ Q ] = [ a 2 k m 0 , a 2 k m 1 , ... , a 2 k m i , ... , a 2 k m ( Q - 1 ) ] ;
Wherein, M is natural number, according to actual conditions value; M is the label of physiology data type, 0≤m≤M-1; L mrepresent the quantity of m+1 class physiological data sensor; N is the group number of each image data, and N is L mintegral multiple; K is the label often organizing data, 0≤k≤N-1; Q is the total number of measured value often organizing data; J represents the label of data acquisition order, j ∈ { 1,2}; I is the label of middle element, 0≤i≤Q-1, for in the i-th+1 element;
Step 2, to m+1 class physiological data element, utilize Grobus criterion carry out first time screening, concrete steps are as follows:
Step 2-1, calculates middle element arithmetic mean a j k m ‾ = 1 Q Σ i = 0 Q - 1 a j k m i , And then calculate its standard deviation
σ j k m = Σ i = 0 Q - 1 ( a j k m i - a j k m ‾ ) 2 Q - 1 ;
Step 2-2, by Ge Luobusi tables of critical values, chooses a critical value d (Q, θ), and the size of d (Q, θ) is determined by the value of Q and significance θ; Wherein significance θ gets 0.05;
Step 2-3, the g (i) that distributed by Ge Luobusi carries out comparing for Q time with d (Q, θ), wherein, if g (i)>=d (Q, θ), be in gross error, reject; Result after rejecting is designated as wherein, Q' is after middle element screens, total number of remaining element, Q'≤Q;
Step 3, to the data after first time screening element, carry out first time data fusion, adopt the method in batches estimated, calculate the fusion value often organizing data concrete formula is as follows:
R j k m = σ e 2 f ‾ + σ f 2 e ‾ σ e 2 + σ f 2
Wherein, for middle element superscript the mean value of element set for middle element superscript the mean value of element set, for middle element superscript the variance of element set, for middle element superscript the variance of element set; expression is not more than maximum integer;
Step 4, for different j, k, each group data fusion value difference, represents all fusion values set ψ: wherein, v is the label of element in described set ψ, 0≤v≤2N-1; Again utilize Grobus criterion, carry out programmed screening to element in set ψ, reject wherein gross error, the result after screening is designated as set for described set total number of middle element, 0≤S m≤ 2N; for described set middle element numerals, described set in individual element with represent;
Step 5, adopts least square method of weighting, to described set element carry out second time data fusion, be specially:
By described set element according to formula:
R ~ m = Σ v ‾ = 0 S m - 1 w v ‾ r v ‾
Merge, wherein, for last fusion value, for weights coefficient, computing formula be:
w v ‾ = 1 σ v ‾ Σ v ‾ = 0 S m - 1 1 σ v ‾
And weights coefficient meets: wherein, for the standard deviation of element in corresponding array; Last fusion value as the end value after the process of m+1 class physiological data.
Beneficial effect: a kind of data processing method based on health monitoring that the present invention proposes, first the method utilizes the collecting method of the many sensings of multiconductor, ensure that the accuracy of measurement.Secondly, to the many groups physiological data collected, screen according to Grobus criterion, to remove the rough error data collected.Finally, by being combined algorithm for estimating and least-squares algorithm in batches, the data collected being merged, reducing different sensors because of measuring accuracy difference, on the impact that result causes.The method is conducive to the physiologic information obtaining user accurately, provides effective foundation for health monitoring system detects in real time.
Accompanying drawing explanation
Fig. 1 is that a kind of data processing method based on health monitoring performs schematic flow sheet.
Fig. 2 is the partial crit value table of Grobus criterion.
Fig. 3 is LS weighted fusion algorithm algorithm model figure.
Embodiment
In order to a kind of data processing method based on health monitoring that more detailed description the present invention proposes, by reference to the accompanying drawings, illustrate as follows:
Fig. 1 shows a kind of data processing method based on health monitoring and performs schematic flow sheet.This method mainly comprises three contents: one is the collecting method utilizing the many sensings of multiconductor, gathers user's physiologic information; Two is the many groups physiological datas to collecting, and carries out twice screening in front and back; Three is be combined algorithm for estimating and least-squares algorithm in batches to merge the data after screening.
1. gathered the M class physiological data of person to be measured by sensor, M gets 2, and physiological data can be specially temperature data and blood pressure data.Use infrared sensor and thermistor (temperature) sensor to gather the temperature data of person to be measured, use resistive pressure sensor and capacitance pressure transducer, to gather the blood pressure data of person to be measured, its detailed process is:
1st step: first time single contact collection, each sensor is not measured in the same time respectively four of same contact, obtains 4 groups of physiological datas, often organizes data and be expressed as:
2nd step: second time multi-contact data collection, each sensor is not measured in the same time at four of four contacts respectively, obtains 4 groups of physiological datas, often organizes data and is expressed as:
M is the label of physiology data type, because only take temperature and blood pressure parameter, so 0≤m≤1; Namely the first kind is temperature data, and Equations of The Second Kind is blood pressure data.K is the label often organizing data, 0≤k≤3; J represents the label of data acquisition order, j ∈ { 1,2}; I is the label of middle element, 0≤i≤3, in the i-th+1 element.
It should be noted that, this enforcement gathers at every turn and only gathers 4 groups of data, often organize data and contain 4 measured values, also more multi-group data can be gathered according to actual conditions, gather from the physiological data of more angle to user, be that the accuracy of data will be higher because the data gathered are more, the health monitoring of user could accurately.
2. pair many groups collected physiological data, carries out first time screening:
1st step: differentiation and the rejecting of the data separate collected " Grobus criterion " being carried out to gross error; Concrete grammar is:
Calculate middle element arithmetic mean and then calculating standard deviation:
σ j k m = Σ i = 0 3 ( a j k m i - a j k m ‾ ) 2 3 - - - ( 1 )
2nd step: by the tables of critical values of Fig. 2 Grobus criterion, choose a critical value d (Q, α), its size is determined by the value of Q and α; Wherein significance α gets 0.05; Table look-up known d (4,0.05)=1.46, according to Ge Luobusi distribution g (i):
g ( i ) = | a j k m i - a j k m ‾ | σ j k m - - - ( 2 )
G (i) and d (4,0.05) are carried out 4 times compare, if g (i)>=d (4,0.05), be in gross error, reject; Result after screening is designated as wherein, Q' is after middle element screens, total number of remaining element.
3. adopt algorithm for estimating in batches to carry out first time fusion to each group of data, detailed process is as follows:
Adopt the method in batches estimated, to often organizing data respectively according to formula:
R j k m = σ e 2 f ‾ + σ f 2 e ‾ σ e 2 + σ f 2 - - - ( 3 )
Calculate the fusion value of data, wherein, for fusion value, for middle element superscript the mean value of element set, for middle element superscript the mean value of element set, for middle element superscript the variance of element set, for middle element superscript the variance of element set;
4: for different j, k, each group data fusion value difference, represents all fusion values set ψ: wherein, v is the label of element in described set ψ, 0≤v≤7; Again utilize Grobus criterion, carry out programmed screening to element in set ψ, reject wherein gross error, the result after screening is designated as set s mfor described set total number of middle element, 0≤S m≤ 8; for described set middle element numerals, described set in individual element with represent;
It should be noted that, why carry out differentiation and the rejecting of twice gross error, because primary process is just carried out for single-sensor, if one or several sensor degradation wherein, the data obtained is to sensor it doesn't matter itself gross error, and can fusion results be drawn, secondary data processing just can weed out having the sensing data of damage.
5: on the basis of previous step, adopt least square method of weighting, to set middle element carries out second time data fusion, and detailed process is:
By described set element according to formula:
R ~ m = Σ v ‾ = 0 S m - 1 w v ‾ r v ‾ - - - ( 4 )
Merge, wherein, for last fusion value, for weights coefficient, computing formula be::
w v ‾ = 1 σ v ‾ Σ v ‾ = 0 S m - 1 1 σ v ‾ - - - ( 5 )
And weights coefficient meets: for the standard deviation of element in corresponding array; Last fusion value as the end value after the process of m+1 class physiological data, this example is the last processing costs of body temperature, blood pressure.
It should be noted that, LS weighted fusion algorithm algorithm is adopted to carry out secondary fusion to data, because merge between the sensor measurement data estimating to be applicable to equally accurate in batches, owing to using many sensing measurements data, the precision of sensor is not exclusively the same, in order to make the result of fusion more excellent, so use minimum weight blending algorithm to carry out secondary fusion to data.Finally, the data of representative of consumer physiological parameter are drawn.
Comprehensive foregoing teachings, the execution flow process that can obtain this method is:
The first step, utilizes the collecting method of the many sensings of multiconductor, gathers user's physiologic information.
Second step, the many groups collected physiological data, carries out first time data screening according to Grobus criterion.
3rd step, to the often group data after screening, utilizes algorithm for estimating in batches to carry out first time data fusion.
4th step, to the result after fusion, according to the grouping of measurement parameter classification, carries out secondary data screening in group.
5th step, to the data after postsearch screening, adopts least square method to carry out the fusion of data secondary, the final measured value of last result representative of consumer physiological parameter.

Claims (1)

1. based on a data processing method for health monitoring, it is characterized in that, comprise the steps:
Step 1: the M class physiological data being gathered person to be measured by sensor, uses L mindividual m+1 class physiological data sensor, gather the m+1 class physiological data of person to be measured, its detailed process is:
Step 1-1, gather for the 1st time of physiological data, each sensor is not measured in the same time respectively the N number of of same contact, obtains N group physiological data, often organizes data and is expressed as:
Step 1-2, gather for the 2nd time of physiological data, each sensor is respectively at L mthe N number of of individual contact does not measure in the same time, obtains N group physiological data, often organizes data and is expressed as:
Wherein, M is natural number, according to actual conditions value; M is the label of physiology data type, 0≤m≤M-1; L mrepresent the quantity of m+1 class physiological data sensor; N is the group number of each image data, and N is L mintegral multiple; K is the label often organizing data, 0≤k≤N-1; Q is the total number of measured value often organizing data; J represents the label of data acquisition order, j ∈ { 1,2}; I is the label of middle element, 0≤i≤Q-1, for in the i-th+1 element;
Step 2, to m+1 class physiological data element, utilize Grobus criterion carry out first time screening, concrete steps are as follows:
Step 2-1, calculates middle element arithmetic mean and then calculate its standard deviation
Step 2-2, by Ge Luobusi tables of critical values, chooses a critical value d (Q, θ), and the size of d (Q, θ) is determined by the value of Q and significance θ; Wherein significance θ gets 0.05;
Step 2-3, the g (i) that distributed by Ge Luobusi carries out comparing for Q time with d (Q, θ), wherein, if g (i)>=d (Q, θ), be in gross error, reject; Result after rejecting is designated as wherein, Q' is after middle element screens, total number of remaining element, Q'≤Q;
Step 3, to the data after first time screening element, carry out first time data fusion, adopt the method in batches estimated, calculate the fusion value often organizing data concrete formula is as follows:
Wherein, for middle element superscript the mean value of element set, for middle element superscript the mean value of element set, for middle element superscript the variance of element set, for middle element superscript the variance of element set; expression is not more than maximum integer;
Step 4, for different j, k, each group data fusion value difference, represents all fusion values set ψ: wherein, v is the label of element in described set ψ, 0≤v≤2N-1; Again utilize Grobus criterion, carry out programmed screening to element in set ψ, reject wherein gross error, the result after screening is designated as set s mfor described set total number of middle element, 0≤S m≤ 2N; for described set middle element numerals, ; Described set in individual element with represent;
Step 5, adopts least square method of weighting, to described set element carry out second time data fusion, be specially:
By described set element according to formula:
Merge, wherein, for last fusion value, for weights coefficient, computing formula be:
And weights coefficient meets: wherein, for the standard deviation of element in corresponding array; Last fusion value as the end value after the process of m+1 class physiological data.
CN201510897273.5A 2015-12-08 2015-12-08 A kind of data processing method based on health monitoring Active CN105550502B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510897273.5A CN105550502B (en) 2015-12-08 2015-12-08 A kind of data processing method based on health monitoring

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510897273.5A CN105550502B (en) 2015-12-08 2015-12-08 A kind of data processing method based on health monitoring

Publications (2)

Publication Number Publication Date
CN105550502A true CN105550502A (en) 2016-05-04
CN105550502B CN105550502B (en) 2018-04-03

Family

ID=55829690

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510897273.5A Active CN105550502B (en) 2015-12-08 2015-12-08 A kind of data processing method based on health monitoring

Country Status (1)

Country Link
CN (1) CN105550502B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106202860A (en) * 2016-06-23 2016-12-07 南京邮电大学 A kind of mood regulation service push method and wearable collaborative supplying system
CN107436997A (en) * 2017-07-03 2017-12-05 上海百纬健康科技有限公司 The analysis system and method for a kind of physiological data
CN113806934A (en) * 2021-09-16 2021-12-17 浙江衡玖医疗器械有限责任公司 Hemispherical temperature sensor prediction model construction method and device and electronic equipment

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101477143A (en) * 2009-01-13 2009-07-08 国网电力科学研究院 Electronic current mutual inductor data processing method based on multi-sensor data amalgamation technology
CN101923789A (en) * 2010-03-24 2010-12-22 北京航空航天大学 Safe airplane approach method based on multisensor information fusion
CN102567640A (en) * 2011-12-29 2012-07-11 上海电机学院 Method for monitoring mine gas
WO2013134845A1 (en) * 2012-03-13 2013-09-19 Hongyue Luo Wearable miniature health monitoring system and method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101477143A (en) * 2009-01-13 2009-07-08 国网电力科学研究院 Electronic current mutual inductor data processing method based on multi-sensor data amalgamation technology
CN101923789A (en) * 2010-03-24 2010-12-22 北京航空航天大学 Safe airplane approach method based on multisensor information fusion
CN102567640A (en) * 2011-12-29 2012-07-11 上海电机学院 Method for monitoring mine gas
WO2013134845A1 (en) * 2012-03-13 2013-09-19 Hongyue Luo Wearable miniature health monitoring system and method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吴杰: "《基于虚拟仪器的温室监测系统的研究》", 《中国优秀硕士学位论文全文数据库信息科技辑》 *
范满红 等;: "《基于多传感器数据融合的温湿度监测系统》", 《压电与声光》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106202860A (en) * 2016-06-23 2016-12-07 南京邮电大学 A kind of mood regulation service push method and wearable collaborative supplying system
CN106202860B (en) * 2016-06-23 2018-08-14 南京邮电大学 A kind of mood regulation service push method
CN107436997A (en) * 2017-07-03 2017-12-05 上海百纬健康科技有限公司 The analysis system and method for a kind of physiological data
CN113806934A (en) * 2021-09-16 2021-12-17 浙江衡玖医疗器械有限责任公司 Hemispherical temperature sensor prediction model construction method and device and electronic equipment

Also Published As

Publication number Publication date
CN105550502B (en) 2018-04-03

Similar Documents

Publication Publication Date Title
CN206026334U (en) Motion amount detection device and intelligent wearable equipment comprising same
CN101847069A (en) Multi-point touch detection method of touch screen
CN105550502A (en) Data processing method based on health monitoring
CN102697225A (en) Electronic measurement clothes and measuring method
CN105243285A (en) Big data health forecast system
CN102670190A (en) Heart rate variability nonlinear characteristic-based automatic diagnosis method for congestive heart failure
CN103838963A (en) Bra pressure comfort evaluation method
López-Nava et al. Estimation of temporal gait parameters using Bayesian models on acceleration signals
CN107966161A (en) Walking detection method based on FFT
CN105771224B (en) A kind of locomotion evaluation system based on multisensor
Flutur et al. Smart chair system for posture correction
CN107411702A (en) A kind of method and system for testing Wrist wearable type terminal heart rate detection precision
CN205120130U (en) FFT calculates equipment of step number and human motion consumption calorie based on improve
CN204085569U (en) A kind of hot comfort instrument for measuring index based on Internet of Things and neural network
CN207236771U (en) A kind of infrared body temperature detector and infrared body temperature detecting system
CN205333175U (en) Healthy electronic scale based on thing networking
CN105528857B (en) A kind of intelligent remote sign data harvester
CN106618570A (en) Skin biochemical index detection method and system based on biological dielectric spectrum
CN106971361A (en) A kind of professional ability remote test system and professional ability remote test method
Pan et al. Heterogeneous sensor data fusion for human falling detection
CN113825253B (en) Intelligent electronic scale monitoring system and method based on big data
CN110338798A (en) Human body respiration amount assessment system and method based on RFID label tag
CN109284936A (en) A kind of electricity quality evaluation method based on the output of cloud graph visualization
CN115752800A (en) Medical wireless body temperature monitoring system based on internet
CN105138835A (en) Human body composition prediction method based on physiological information entropy

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
TR01 Transfer of patent right

Effective date of registration: 20181102

Address after: 100071 6 floor 612, 1 building, 9 Guang'an Road, Fengtai District, Beijing.

Patentee after: Beijing Qi Huang clinical research center

Address before: 210003 new model road 66, Gulou District, Nanjing, Jiangsu

Patentee before: Nanjing Post & Telecommunication Univ.

TR01 Transfer of patent right