CN104545808B - Method and system for calculating delay time of blood sugar - Google Patents

Method and system for calculating delay time of blood sugar Download PDF

Info

Publication number
CN104545808B
CN104545808B CN201410723343.0A CN201410723343A CN104545808B CN 104545808 B CN104545808 B CN 104545808B CN 201410723343 A CN201410723343 A CN 201410723343A CN 104545808 B CN104545808 B CN 104545808B
Authority
CN
China
Prior art keywords
blood glucose
calculating
scattering coefficient
time
time delay
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
CN201410723343.0A
Other languages
Chinese (zh)
Other versions
CN104545808A (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 CHINA LIGHT TECHNOLOGY CO.,LTD.
Original Assignee
SUZHOU OPTORING TECHNOLOGY Co Ltd
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 SUZHOU OPTORING TECHNOLOGY Co Ltd filed Critical SUZHOU OPTORING TECHNOLOGY Co Ltd
Priority to CN201410723343.0A priority Critical patent/CN104545808B/en
Publication of CN104545808A publication Critical patent/CN104545808A/en
Application granted granted Critical
Publication of CN104545808B publication Critical patent/CN104545808B/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/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • A61B5/14551Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biomedical Technology (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Optics & Photonics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Abstract

The invention discloses a method and a system for calculating the delay time of blood sugar. According to the method and the system for calculating the delay time of the blood sugar, the delay time of the blood sugar in different depth areas of the skin can be accurately calculated, and a depth area suitable for serving as a blood sugar detection area can be accurately selected according to the delay time, so that the detection accuracy of the blood sugar can be improved.

Description

Calculate the method and system of blood glucose time delay
Technical field
The present invention relates to blood sugar test technical field, more particularly to a kind of calculate the method for blood glucose time delay and be System.
Background technology
Diabetes are a kind of systemic diseases being basic biochemical character with persistent high blood sugar, and it is based on carbohydrate metabolism disturbance Show it is possible to cause multiple complications.If diabetes do not obtain enough control, some can be caused acute concurrent Disease, such as hypoglycemia, ketoacidosiss, non-ketone hyperosmolar coma.Blood sugar test is one of important Testing index of diabetes, is Avoid the generation of complication, client need real-time self-monitoring of blood glucose.Optics noninvasive dynamics monitoring technology mainly includes at present Near infrared spectrum, in (remote) infrared spectrum, optical coherence tomography, Raman spectrum, optoacoustic spectroscopy and polarized light optically-active etc..Its advantage It is not use biological pests, such as gathers blood, body fluid etc. without stabbing human body, long-term continuous detecting can be carried out.Cause This, have very big research and practical value, is the trend of blood glucose measurement research and development from now on.
Oct (optical coherence tomography, optical coherence tomography) is a kind of based on weak relevant former The noncontact of reason, Non-invasive technology, resolution can reach micron dimension.By measuring inside skin histology because of glucose Light attenuation coefficient change caused by value difference, can calculate its concentration of glucose.Human body skin is mainly by epidermal area, skin corium Constitute with subcutaneous tissue.Because the component of organization of each layer is different, its optical property parameter is with the dependency of change of blood sugar, variable quantity Also there is difference.The internal main scattering object of skin is cell membrane, nucleus and protein aggregate, and these organizational structuries are in skin Distribution be random, the ratio that it accounts in different skin regional location and exist position be different, therefore in different skins Can there is certain difference in the one-dimensional signal of skin regional location.It is thus desirable to according to skin histology depth direction structure distribution Make a concrete analysis of.And the advantage of oct is, essence can be carried out to the optical property parameter change in skin histology different depth region Really measure, thus finding and change of blood sugar maximally related tissue depth region, carrying out demarcating, predicting.
The blood glucose value being used for demarcating in noninvasive dynamics monitoring generally refers to blood or venous blood.When blood sugar for human body becomes quick During change, the change of blood sugar in its dermis of skin layer tissue can lag behind the change of blood sugar in peripheral blood or venous blood.And noinvasive blood Sugar detection technology changes to calculate blood glucose value generally by optical property parameter in detection skin corium, that is, measure is skin corium The change of blood sugar situation of tissue.Inspection for hyperglycemia or hypoglycemic, caused by this delay due to physiology It is breakneck for surveying result inaccurate.
Content of the invention
Based on this it is necessary to be directed to the problems referred to above, provide a kind of method and system calculating blood glucose time delay, this calculating The method and system of blood glucose time delay can calculate blood glucose time delay so that blood glucose test results are more accurate.
A kind of method calculating blood glucose time delay, comprise the steps: to obtain in the same depth of sample not in the same time Blood glucose value and accordingly obtain the time;Calculate first scattering coefficient of the correspondence each acquisition time in described same depth;Meter Calculate the second scattering coefficient in described same depth, described second scattering coefficient is from first acquisition time, to Not calculating in the same time in the time period that acquisition time terminates afterwards;Calculate the second scattering coefficient related to blood glucose value Property;Dependency described in com-parison and analysis, draws blood glucose time delay.
Wherein in an embodiment, described calculating the second scattering coefficient with the step of the dependency of blood glucose value is: according to Pearson correlation calculates the dependency of the second scattering coefficient and blood glucose value.
Wherein in an embodiment, dependency described in described com-parison and analysis, show that the step of blood glucose time delay is: blood It is equal to corresponding time during correlation coefficient maximum sugared time delay.
Wherein in an embodiment, the step of the second scattering coefficient in the described same depth of described calculating is: every two The interval time in the individual adjacent moment of described second scattering coefficient of calculating is identical.
Wherein in an embodiment, before the described not blood glucose value step in the same time obtaining in the same depth of sample also Comprise the steps: the dependency of blood glucose value and scattering coefficient on analytical calculation different depth, and select scattering coefficient and blood glucose The depth areas of value correlation maximum.
A kind of system calculating blood glucose time delay, including blood glucose value acquisition module, for obtaining the same depth of sample Not blood glucose value in the same time, and obtain the described acquisition time;Scattering coefficient computing module, for calculating described same depth Correspondence each acquisition time the first scattering coefficient;Scattering coefficient comparing module, for calculating second in described same depth Scattering coefficient, described second scattering coefficient is from first acquisition time, to last acquisition time terminate when Between not calculating in the same time in section;Correlation coefficient computing module is related to blood glucose value for calculating the second scattering coefficient Property;Com-parison and analysis module, for dependency described in com-parison and analysis, draws blood glucose time delay.
Wherein in an embodiment, described correlation coefficient computing module calculates the second scattering according to Pearson correlation Coefficient and the dependency of blood glucose value.
Wherein in an embodiment, in described com-parison and analysis module, be equal to blood glucose time delay correlation coefficient maximum when The corresponding time.
Wherein in an embodiment, in described scattering coefficient comparing module, the adjacent calculating of each two described second dissipates The interval time penetrating the moment of coefficient is identical.
Wherein in an embodiment, the described system calculating blood glucose time delay also includes: depth selecting module, is used for The dependency of blood glucose value and scattering coefficient on analytical calculation different depth, and select scattering coefficient and blood glucose value correlation maximum Depth areas.
The method and system of above-mentioned calculating blood glucose time delay can calculate blood glucose on skin different depth region exactly Time delay, according to can accurately select the depth areas being suitable as blood sugar test region this time delay such that it is able to Improve the accuracy of blood sugar test.
Brief description
Fig. 1 is the flow chart of the method calculating blood glucose time delay;
Fig. 2 is oct 3-D view;
Fig. 3 is oct one-dimensional signal figure;
Fig. 4 is skin light attenuation coefficient with blood glucose delay variation trendgram;
Fig. 5 is the system framework figure calculating blood glucose time delay.
Specific embodiment
Understandable for enabling the above objects, features and advantages of the present invention to become apparent from, below in conjunction with the accompanying drawings to the present invention Specific embodiment be described in detail.Elaborate a lot of details in order to fully understand this in the following description Bright.But the present invention can be much to implement different from alternate manner described here, and those skilled in the art can be not Similar improvement is done, therefore the present invention is not limited by following public specific embodiment in the case of running counter to intension of the present invention.
Refer to Fig. 1, provide a kind of method calculating blood glucose time delay in one embodiment, specifically include as Lower step:
S110: the dependency of blood glucose value and scattering coefficient on analytical calculation different depth, and select scattering coefficient and blood glucose The depth areas of value correlation maximum.
Refer to Fig. 2 and Fig. 3 it is necessary first to all a-scan in the 3-D view of oct machine scanning acquisition will be had to be put down All, such that it is able to obtain oct one-dimensional depth signal figure (i.e. Fig. 3).(in other embodiments, differ from the beginning of skin surface Determine to start it is also possible to from skin surface to from the beginning of bottom offset any distance, be not specifically limited herein from skin surface), successively Selected depth be 1*w, 2*w ..., n*w (n=1,2,3 ...) detection depth areas.In each depth areas, acquire not T in the same timeiBlood glucose value gli, and it is calculated each moment tiScattering coefficient μti(i is the integer more than 0).Further according to such as Lower formula (1), you can obtain the correlation coefficient of scattering coefficient and blood glucose value on different depth region.
r = n σ i = 1 n μ t i gl i - σ i = 1 n μ t i σ i = 1 n gl i n σ i = 1 n μ t i 2 - ( σ i = 1 n μ t i ) 2 n σ i = 1 n gl i 2 - ( σ i = 1 n gl i ) 2 - - - ( 1 )
Wherein, μtiFor scattering coefficient, gliFor blood glucose value.Will be according on above-mentioned formula calculated different depth region Scattering coefficient is compared with the correlation coefficient of blood glucose value, thus selecting detection depth areas.Specifically, the standard of selection can With self-defined, for example, correlation coefficient can be selected maximum, that is, the depth areas corresponding when correlation coefficient is maximum are as follow-up Calculate the detection depth areas of blood glucose time delay;As long as it is corresponding in certain numerical intervals to define correlation coefficient r Depth areas are as the detection depth areas of subsequent calculations blood glucose time delay.
This step is to calculate the previous step of blood glucose time delay, and not necessarily step.If executing this step, can There is provided facility for subsequent calculations blood glucose time delay.Certainly, do not execute this step, random when subsequent calculations blood glucose time delay Choose the calculating that detection depth areas can also realize blood glucose time delay.
S120: obtain the not blood glucose value in the same time in the same depth of sample and obtain the time accordingly.Described in this step Same depth be depth areas selected by previous step s110, and subsequently several step is all enterprising in this depth areas Row detection and analysis calculates.In t1、t2、…、tiIn the moment, the blood glucose value gathering the corresponding moment respectively is designated as gl1、gl2、…、gli, And write down the acquisition time, wherein, i is the integer more than 0.t1For obtaining blood glucose value gliFirst acquisition the time.
S130: calculate first scattering coefficient of the correspondence each acquisition time in described same depth.Corresponding to previous step S120 obtains each acquisition time of blood glucose value, calculates the first scattering coefficient μ of each acquisition time correspondingti.I.e. in t1、 t2、…、tiIn the moment, corresponding first scattering coefficient is μ respectivelyt1、μt2、…、μti.
S140: calculate the second scattering coefficient in described same depth.The second scattering coefficient μ in this steptlagBe Obtain time (t from first1) start, time period (the i.e. t terminating to last acquisition time1To tiBetween) in difference Moment calculates.Each two adjacent calculating the second scattering coefficient μtlagMoment interval time identical.With t1During for initiateing Carve (identical with the initial time of previous step), δ t is interval time, so, corresponding calculating the second scattering coefficient μtlagRight The moment answered is followed successively by t1、t1+δt、t1+2*δt、t1+ 3 δ t ..., corresponding second scattering coefficient μtlag, it is designated as successively μtlag1、μtlag2、μtlag3、μtlag4….
S150: calculate the dependency of the second scattering coefficient and blood glucose value.In this step, calculated according to Pearson correlation Second scattering coefficient μtlagDependency with blood glucose value.Pearson correlation computing formula equation below (2):
r = n σ i n ( μ t i - μ t l a g ) gl i - σ i n ( μ t i - μ t l a g ) σ i n gl i n σ i n ( μ t i - μ t l a g ) 2 - ( σ i n ( μ t i - μ t l a g ) ) 2 n σ i n gl i 2 - ( σ i n gl i ) 2 - - - ( 2 )
Wherein, μtiExpression sampling instant is tiWhen corresponding the first scattering coefficient, gliFor tiThe blood glucose value in moment, μtlag For the second scattering coefficient calculating in step s140.Work as tlagFor t1When, calculate correlation coefficient now, be designated as r1; Work as tlagFor t1During+δ t, calculate correlation coefficient now, be designated as r2;tlagFor t1During+2* δ t, calculate phase relation now Number, is designated as r3;..., the like, obtain one group of correlation coefficient.
S160: dependency described in com-parison and analysis, draws blood glucose time delay.Each t that previous step is drawnlagCorresponding Correlation coefficient (one group of correlation coefficient) be analyzed comparing.On this skin depth region, blood glucose time delay is one group of phase Close the maximum t corresponding to correlation coefficient in coefficientlag.
After the skin different depth blood glucose time delay being obtained by the method for above-mentioned calculating blood glucose time delay, Ke Yixuan Select the depth areas that time delay is short and dependency is high as the depth areas of blood glucose real-time detection, to reach raising blood sugar test The purpose of accuracy.Or scattering coefficient can be directly selected and carry out follow-up blood glucose with blood glucose dependency highest depth areas Prediction, by compensating calculating time delay in the calculated blood glucose in this region, to obtain accurately overall blood glucose prediction Trendgram.
The randomness being distributed due to individual difference and the skin interior tissue material of people, by detecting calculated skin The blood glucose of skin different depth is only applicable to individual and can not be general time delay.I.e. different people carries out optics noninvasive dynamics monitoring Before, it is required for first carrying out the calculating of blood glucose time delay.
Refer to Fig. 4, in ogtt experiment, certain testee's skin scattering coefficient is with the delay trendgram of change of blood sugar. Skin depth region is 200um~325um, and this region is located at the papillary layer of the skin corium first half, and is in negative with change of blood sugar Guan Xing.Figure 4, it is seen that scattering coefficient and blood sugar for human body have certain retardance, and time delay is 17 minutes.
The method of above-mentioned calculating blood glucose time delay can accurately calculate the delay of blood glucose on skin different depth region Time, after learning this blood glucose time delay, then carry out blood sugar test, you can reach the purpose improving blood sugar test accuracy.
Refer to Fig. 5, additionally provide a kind of system 100 calculating blood glucose time delay in one embodiment.This meter The system 100 calculating blood glucose time delay includes depth selecting module 110, blood glucose value acquisition module 120, scattering coefficient computing module 130th, scattering coefficient comparing module 140, correlation coefficient computing module 150 and com-parison and analysis module 160.Depth selecting module 110 For the dependency of blood glucose value on analytical calculation different depth and scattering coefficient, and select scattering coefficient with blood glucose value dependency Big depth areas;Blood glucose value acquisition module 120 is used for obtaining not blood glucose value in the same time in the same depth of sample and corresponding The acquisition time;What scattering coefficient computing module 130 was used for calculating correspondence each acquisition time in described same depth first dissipates Penetrate coefficient;Scattering coefficient comparing module 140 is used for calculating the second scattering coefficient in described same depth, described second scattering system Number μtlagIt is in not calculating in the same time within first acquisition time, the time period terminating to last acquisition time Go out;Correlation coefficient computing module 150 is used for calculating the second scattering coefficient μtlagDependency with blood glucose value;Com-parison and analysis module 160 are used for dependency described in com-parison and analysis, draw blood glucose time delay.
Depth selecting module 110 is used for the dependency of blood glucose value and the second scattering coefficient on analytical calculation different depth, and Select the depth areas of scattering coefficient and blood glucose value correlation maximum.Firstly the need of the 3-D view that the scanning of oct machine will be had to obtain In all a-scan carry out averagely, such that it is able to obtain oct one-dimensional depth signal figure (i.e. Fig. 3).From the beginning of skin surface, successively Selected depth be 1*w, 2*w ..., n*w (n=1,2,3 ...) detection depth areas.In each depth areas, collect each Individual moment tiBlood glucose value gli, and it is calculated each moment tiScattering coefficient μti(i is the integer more than 0).Further according to such as Lower formula (1), you can obtain the correlation coefficient of scattering coefficient and blood glucose value on different depth region.
r = n σ i = 1 n μ t i gl i - σ i = 1 n μ t i σ i = 1 n gl i n σ i = 1 n μ t i 2 - ( σ i = 1 n μ t i ) 2 n σ i = 1 n gl i 2 - ( σ i = 1 n gl i ) 2 - - - ( 1 )
Wherein, μtiFor scattering coefficient, gliFor blood glucose value.Will be according on above-mentioned formula calculated different depth region Scattering coefficient is compared with the correlation coefficient of blood glucose value, thus selecting detection depth areas.Specifically, the standard of selection can With self-defined, for example, correlation coefficient can be selected maximum, that is, the depth areas corresponding when correlation coefficient is maximum are as follow-up Calculate the detection depth areas of blood glucose time delay;As long as it is corresponding in certain numerical intervals to define correlation coefficient r Depth areas are as the detection depth areas of subsequent calculations blood glucose time delay.
Blood glucose value acquisition module 120 is used for obtaining the not blood glucose value in the same time in the same depth of sample and corresponding acquisition Time.In t1、t2、…、tiIn the moment, the blood glucose value gathering the corresponding moment respectively is designated as gl1、gl2、…、gli, and when writing down acquisition Between, wherein, i is the integer more than 0.t1For obtaining blood glucose value gliFirst acquisition the time.
Scattering coefficient computing module 130 is used for calculating the first scattering system of the correspondence each acquisition time in described same depth Number μti.Each corresponding to collection blood glucose value obtains the time, calculates corresponding each acquisition time corresponding first scattering coefficient μti.I.e. in t1、t2、…、tiIn the moment, corresponding first scattering coefficient is μ respectivelyt1、μt2、…、μti.
Scattering coefficient comparing module 140 is used for calculating the second scattering coefficient in described same depth, described second scattering Coefficient is in not calculating in the same time within first acquisition time, the time period terminating to last acquisition time 's.Second scattering coefficient μtlagIt is from first acquisition time (t1) start, the time period terminating to last acquisition time (i.e. t1To tiBetween) in do not calculate in the same time.Each two adjacent calculating the second scattering coefficient μtlagMoment between Identical every the time.With t1It is interval time for initial time, δ t, so, corresponding calculating the second scattering coefficient μtlagCorresponding Moment is followed successively by t1、t1+δt、t1+2*δt、t1+ 3 δ t ..., corresponding second scattering coefficient μtlag, it is designated as μ successivelytlag1、 μtlag2、μtlag3、μtlag4….
Correlation coefficient computing module 150 is used for calculating the dependency of the second scattering coefficient and blood glucose value.According to Pearson's phase Closing property calculates the second scattering coefficient μtlagDependency with blood glucose value.Pearson correlation computing formula equation below (2):
r = n σ i n ( μ t i - μ t l a g ) gl i - σ i n ( μ t i - μ t l a g ) σ i n gl i n σ i n ( μ t i - μ t l a g ) 2 - ( σ i n ( μ t i - μ t l a g ) ) 2 n σ i n gl i 2 - ( σ i n gl i ) 2 - - - ( 2 )
Wherein, μtiExpression sampling instant is tiWhen corresponding the first scattering coefficient, gliFor tiThe blood glucose value in moment, μtlag The second scattering coefficient calculating for scattering coefficient comparing module 140.Work as tlagFor t1When, calculate phase relation now Number, is designated as r1;Work as tlagFor t1During+δ t, calculate correlation coefficient now, be designated as r2;tlagFor t1During+2* δ t, calculate this When correlation coefficient, be designated as r3;..., the like, obtain one group of correlation coefficient.
Com-parison and analysis module 160 is used for dependency described in com-parison and analysis, draws blood glucose time delay.Correlation coefficient is calculated Module 150 show that each correlation coefficient is analyzed comparing.On this skin depth region, blood glucose is related for one group for time delay The maximum t corresponding to correlation coefficient in coefficientlag.
The system 100 of above-mentioned calculating blood glucose time delay can calculate prolonging of blood glucose on skin different depth region exactly The late time, according to can accurately select the depth areas being suitable as blood sugar test region this time delay such that it is able to improve The accuracy of blood sugar test.
Embodiment described above only have expressed the several embodiments of the present invention, and its description is more concrete and detailed, but simultaneously Therefore the restriction to the scope of the claims of the present invention can not be interpreted as.It should be pointed out that for those of ordinary skill in the art For, without departing from the inventive concept of the premise, some deformation can also be made and improve, these broadly fall into the guarantor of the present invention Shield scope.Therefore, the protection domain of patent of the present invention should be defined by claims.

Claims (10)

1. a kind of method calculating blood glucose time delay is it is characterised in that comprise the steps:
Obtain the not blood glucose value in the same time in the same depth of sample and obtain the time accordingly;
Calculate first scattering coefficient of the correspondence each acquisition time in described same depth;
Calculate the second scattering coefficient in described same depth, described second scattering coefficient is to open from first acquisition time Begin, not calculating in the same time in the time period terminating to last acquisition time;
Calculate the dependency of the second scattering coefficient and blood glucose value;
Dependency described in com-parison and analysis, draws blood glucose time delay.
2. the method calculating blood glucose time delay according to claim 1 is it is characterised in that described calculating second scattering system Number with the step of the dependency of blood glucose value is:
Calculate the dependency of the second scattering coefficient and blood glucose value according to Pearson correlation.
3. the method calculating blood glucose time delay according to claim 1 is it is characterised in that phase described in described com-parison and analysis Guan Xing, show that the step of blood glucose time delay is:
Blood glucose is equal to the corresponding time during correlation coefficient maximum time delay.
4. the method calculating blood glucose time delay according to claim 1 is it is characterised in that the described same depth of described calculating The step of the second scattering coefficient on degree is:
The interval time in each two adjacent moment of described second scattering coefficient of calculating is identical.
5. the method calculating blood glucose time delay according to claim 1 is it is characterised in that the same depth of described acquisition sample Also comprise the steps: before not blood glucose value step in the same time on degree
The dependency of blood glucose value and scattering coefficient on analytical calculation different depth, and select scattering coefficient with blood glucose value dependency Big depth areas.
6. a kind of system calculating blood glucose time delay is it is characterised in that include
Blood glucose value acquisition module, for obtaining the not blood glucose value in the same time in the same depth of sample, and when obtaining described acquisition Between;
Scattering coefficient computing module, for calculating first scattering coefficient of the correspondence each acquisition time in described same depth;
Scattering coefficient comparing module, for calculating the second scattering coefficient in described same depth, described second scattering coefficient is In not calculating in the same time within first acquisition time, the time period terminating to last acquisition time;
Correlation coefficient computing module, for calculating the dependency of the second scattering coefficient and blood glucose value;
Com-parison and analysis module, for dependency described in com-parison and analysis, draws blood glucose time delay.
7. the system calculating blood glucose time delay according to claim 6 is it is characterised in that described correlation coefficient calculates mould Tuber calculates the dependency of the second scattering coefficient and blood glucose value according to Pearson correlation.
8. the system calculating blood glucose time delay according to claim 6 is it is characterised in that described com-parison and analysis module In, blood glucose is equal to the corresponding time during correlation coefficient maximum time delay.
9. the system calculating blood glucose time delay according to claim 6 is it is characterised in that described scattering coefficient compares mould In block, the interval time in each two adjacent moment of described second scattering coefficient of calculating is identical.
10. the system calculating blood glucose time delay according to claim 6 is it is characterised in that described calculating blood glucose postpones The system of time also includes:
Depth selecting module, for the dependency of blood glucose value on analytical calculation different depth and scattering coefficient, and selects scattering system Number and the depth areas of blood glucose value correlation maximum.
CN201410723343.0A 2014-12-02 2014-12-02 Method and system for calculating delay time of blood sugar Active CN104545808B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410723343.0A CN104545808B (en) 2014-12-02 2014-12-02 Method and system for calculating delay time of blood sugar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410723343.0A CN104545808B (en) 2014-12-02 2014-12-02 Method and system for calculating delay time of blood sugar

Publications (2)

Publication Number Publication Date
CN104545808A CN104545808A (en) 2015-04-29
CN104545808B true CN104545808B (en) 2017-01-18

Family

ID=53063676

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410723343.0A Active CN104545808B (en) 2014-12-02 2014-12-02 Method and system for calculating delay time of blood sugar

Country Status (1)

Country Link
CN (1) CN104545808B (en)

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6931327B2 (en) * 2003-08-01 2005-08-16 Dexcom, Inc. System and methods for processing analyte sensor data
WO2014052136A1 (en) * 2012-09-26 2014-04-03 Abbott Diabetes Care Inc. Method and apparatus for improving lag correction during in vivo measurement of analyte concentration with analyte concentration variability and range data

Also Published As

Publication number Publication date
CN104545808A (en) 2015-04-29

Similar Documents

Publication Publication Date Title
WO2018107915A1 (en) General noninvasive blood glucose prediction method based on timing analysis
CN104958075B (en) Use Raman spectrum non-invasive measurement skin thickness and blood sugar concentration and its calibration method
CN103149177B (en) Device and method for detecting biological tissue of pressure modulation near infrared spectrum
Blank et al. Clinical results from a noninvasive blood glucose monitor
CN104224196B (en) The method of non-invasive measurement blood component concentration
CN103948393B (en) A kind of near infrared no-wound detection method of blood constituent content and device
KR20170035675A (en) Method and apparatus for predicting analyte concentration
Li et al. Classification of diabetes and measurement of blood glucose concentration noninvasively using near infrared spectroscopy
Zhang et al. “Guide star” assisted noninvasive photoacoustic measurement of glucose
DE10318393A1 (en) Pulse rate measurement method and device for recording pulse pressure and blood flow
CN111466921A (en) Noninvasive blood glucose detector based on multi-source information sensing and fusion and detection method
CN104188664B (en) Blood sugar test scaling method and system
CN113662538A (en) Non-invasive blood glucose detection method based on time-frequency domain comprehensive analysis
RU2405418C1 (en) Method for diagnostics of endothelial dysfunction in diabetes mellitus patients
US20050151976A1 (en) Method for monitoring of analytes in biological samples using low coherence interferometry
JP2016010717A (en) Concentration quantification apparatus
JP5521199B2 (en) Concentration determination apparatus, concentration determination method, and program
JP5626879B2 (en) Concentration determination apparatus, concentration determination method, and program
Pleitez et al. Infrared reflectometry of skin: Analysis of backscattered light from different skin layers
CN104545808B (en) Method and system for calculating delay time of blood sugar
Chen et al. Non-invasive blood glucose measurement based on mid-infrared spectroscopy
US20120116236A1 (en) Hydration and blood flow adjusted glucose measurement
JP5626880B2 (en) Concentration determination apparatus, concentration determination method, and program
CN106963392B (en) Optical noninvasive blood glucose detection two-dimensional correlation calibration method and system
JP2004321325A (en) Method of quantitating blood glucose level

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20211104

Address after: Room 33036, No. 3, Shuangqing Road, Haidian District, Beijing 100085

Patentee after: BEIJING CHINA LIGHT TECHNOLOGY CO.,LTD.

Address before: 215021 unit 1-b202, creative industrial park, No. 328, Xinghu street, Suzhou Industrial Park, Jiangsu Province

Patentee before: SUZHOU OPTORING TECHNOLOGY Co.,Ltd.

TR01 Transfer of patent right