CN104545808A - 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
CN104545808A
CN104545808A CN201410723343.0A CN201410723343A CN104545808A CN 104545808 A CN104545808 A CN 104545808A CN 201410723343 A CN201410723343 A CN 201410723343A CN 104545808 A CN104545808 A CN 104545808A
Authority
CN
China
Prior art keywords
blood glucose
time
calculating
scattering coefficient
depth
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
CN201410723343.0A
Other languages
Chinese (zh)
Other versions
CN104545808B (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, particularly relate to a kind of method and system calculating blood glucose time delay.
Background technology
The systemic disease of diabetes to be a kind of with persistent high blood sugar be basic biochemical character, it take carbohydrate metabolism disturbance as main manifestations, and can cause multiple complications.If diabetes do not obtain enough control, some acute complicationses can be caused, as hypoglycemia, ketoacidosis, non-ketone hyperosmolar coma.Blood sugar test is one of important Testing index of diabetes, in order to avoid the generation of complication, and the real-time self-monitoring of blood glucose of client need.Current optics noninvasive dynamics monitoring technology mainly comprise near infrared spectrum, in (far) infrared spectrum, optical coherence tomography, Raman spectrum, optoacoustic spectroscopy and polarized light optically-active etc.Its advantage does not use biological pests, need not stab human body and such as gather blood, body fluid etc., can carry out long-term continuous detecting.Therefore, having very large 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 noncontact based on weak relevant principle, Non-invasive technology, and resolution can reach micron dimension.Inner because of the light attenuation coefficient change caused by dextrose equivalent difference by measuring skin histology, its concentration of glucose can be calculated.Human body skin is formed primarily of epidermal area, skin corium and subcutaneous tissue.Because the component of organization of each layer is different, its optical property parameter also has difference with the dependency of change of blood sugar, variable quantity.Skin internal main wants scattering object to be cell membrane, nucleus and protein aggregate, the distribution in skin of these organizational structuries is random, be different in the position of its ratio accounted for of different skin regional location and existence, therefore can there is certain difference at the one-dimensional signal of different skin regional location.Therefore need to make a concrete analysis of according to the structure distribution of skin histology at depth direction.And the advantage of OCT is, accurately can measures the optical property parameter change in skin histology different depth region, thus find the maximally related tissue depth region with change of blood sugar, carry out demarcating, predicting.
Blood glucose value for demarcating in noninvasive dynamics monitoring is generally and refers to blood or venous blood.When blood sugar for human body is at Rapid Variable Design, 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 noninvasive dynamics monitoring technology is generally calculate blood glucose value by detecting optical property parameter change in skin corium, what namely measure is the change of blood sugar situation of skin corium tissue.For hyperglycemia or hypoglycemic, this testing result postponing to cause due to physiology is inaccurate is breakneck.
Summary of the invention
Based on this, be necessary for the problems referred to above, provide a kind of method and system calculating blood glucose time delay, the method and system of this calculating blood glucose time delay can calculate blood glucose time delay, make blood glucose test results more accurate.
Calculate the method for blood glucose time delay, comprise the steps: to obtain the not blood glucose value in the same time in the sample same degree of depth and corresponding acquisition time; Calculate the first scattering coefficient of each acquisition time of correspondence in the described same degree of depth; Calculate the second scattering coefficient in the described same degree of depth, described scattering coefficient is from first acquisition time, not calculating in the same time within the time period that last acquisition time terminates; Calculate the dependency of the second scattering coefficient and blood glucose value; Dependency described in com-parison and analysis, draws blood glucose time delay.
Wherein in an embodiment, the step of the dependency of described calculating second scattering coefficient and blood glucose value is: the dependency calculating the second scattering coefficient and blood glucose value according to Pearson correlation.
Wherein in an embodiment, dependency described in described com-parison and analysis, show that the step of blood glucose time delay is: blood glucose to equal time delay correlation coefficient maximum time corresponding time.
Wherein in an embodiment, the step of the second scattering coefficient in the described same degree of depth of described calculating is: the interval time in the moment of every two adjacent described second scattering coefficients of calculating is identical.
Wherein in an embodiment, also comprise the steps: the dependency of blood glucose value and scattering coefficient on analytical calculation different depth before not blood glucose value step in the same time in the same degree of depth of described acquisition sample, and select scattering coefficient and the maximum depth areas of blood glucose value dependency.
Calculating the system of blood glucose time delay, comprise blood glucose value acquisition module, for obtaining the not blood glucose value in the same time in the same degree of depth of sample, and obtaining described acquisition time; Scattering coefficient computing module, for calculating the first scattering coefficient of each acquisition time of correspondence in the described same degree of depth; Scattering coefficient comparing module, for calculating the second scattering coefficient in the described same degree of depth, described second scattering coefficient is from first acquisition time, not calculating in the same time within the time period that last acquisition time terminates; Calculation of correlation factor 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.
Wherein in an embodiment, described Calculation of correlation factor module calculates the dependency of the second scattering coefficient and blood glucose value according to Pearson correlation.
Wherein in an embodiment, in described com-parison and analysis module, blood glucose to equal time delay correlation coefficient maximum time corresponding time.
Wherein in an embodiment, in described scattering coefficient comparing module, the interval time in the moment of every two adjacent described second scattering coefficients of calculating is identical.
Wherein in an embodiment, the system of described calculating blood glucose time delay also comprises: the degree of depth selects module, the dependency of blood glucose value and scattering coefficient on analytical calculation different depth, and selects scattering coefficient and the maximum depth areas of blood glucose value dependency.
The method and system of above-mentioned calculating blood glucose time delay can calculate the time delay of blood glucose on skin different depth region exactly, can accurately select the depth areas being suitable as blood sugar test region according to this time delay, thus the accuracy of blood sugar test can be improved.
Accompanying drawing explanation
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 that skin light attenuation coefficient is with blood glucose delay variation trendgram;
Fig. 5 is the system framework figure calculating blood glucose time delay.
Detailed description of the invention
For enabling above-mentioned purpose of the present invention, feature and advantage become apparent more, are described in detail the specific embodiment of the present invention below in conjunction with accompanying drawing.Set forth a lot of detail in the following description so that fully understand the present invention.But the present invention can be much different from alternate manner described here to implement, those skilled in the art can when without prejudice to doing similar improvement when intension of the present invention, therefore the present invention is by the restriction of following public specific embodiment.
Please refer to Fig. 1, provide a kind of method calculating blood glucose time delay in one embodiment, specifically comprise the steps:
S110: the dependency of blood glucose value and scattering coefficient on analytical calculation different depth, and select scattering coefficient and the maximum depth areas of blood glucose value dependency.
Please refer to Fig. 2 and Fig. 3, first need to be averaged there being OCT machine to scan all A-scan in the 3-D view obtained, thus OCT one dimension depth signal figure (i.e. Fig. 3) can be obtained.From skin surface (in other embodiments, not necessarily from skin surface, also can from skin surface to bottom offset any distance, be not specifically limited) herein, successively selected depth be 1*w, 2*w ..., n*w (n=1,2,3 ...) detect depth areas.In each depth areas, acquire not t in the same time iblood glucose value gl i, and calculate each moment t iscattering coefficient μ ti(i be greater than 0 integer).Again according to following formula (1), the correlation coefficient of scattering coefficient and blood glucose value on different depth region can be obtained.
R = n Σ i = 1 n μ ti gl i - Σ i = 1 n μ ti Σ i = 1 n gl i n Σ i = 1 n μ ti 2 - ( Σ i = 1 n μ ti ) 2 n Σ i = 1 n gl i 2 - ( Σ i = 1 n gl i ) 2 - - - ( 1 )
Wherein, μ tifor scattering coefficient, gl ifor blood glucose value.The correlation coefficient of scattering coefficient and blood glucose value on the different depth region obtained according to above-mentioned formulae discovery is compared, thus selects detection depth areas.Concrete, the standard of selection can be self-defined, such as, correlation coefficient can be selected maximum, and namely corresponding when correlation coefficient is maximum depth areas is as the detection depth areas of subsequent calculations blood glucose time delay; As long as the detection depth areas of coefficient R depth areas corresponding in certain numerical intervals as subsequent calculations blood glucose time delay also can be defined.
This step is the previous step calculating blood glucose time delay, and and optional step.If perform this step, then can be subsequent calculations blood glucose and provide convenient time delay.Certainly, do not perform this step, when subsequent calculations blood glucose time delay, random selecting detects the calculating that depth areas also can realize blood glucose time delay.
S120: obtain the not blood glucose value in the same time in the sample same degree of depth and corresponding acquisition time.The said same degree of depth of this step is the depth areas selected by previous step S110, and follow-up several step all carries out detection analytical calculation in this depth areas.At t 1, t 2..., t iin the moment, the blood glucose value gathering the corresponding moment is respectively designated as gl 1, gl 2..., gl i, and write down acquisition time, wherein, i be greater than 0 integer.T 1for obtaining blood glucose value gl ifirst acquisition time.
S130: the first scattering coefficient calculating each acquisition time of correspondence in the described same degree of depth.Obtain each acquisition time of blood glucose value corresponding to previous step S120, calculate the first scattering coefficient μ of each acquisition time corresponding ti.Namely at t 1, t 2..., t iin the moment, the first corresponding respectively scattering coefficient is μ t1, μ t2..., μ ti.
S140: calculate the second scattering coefficient in the described same degree of depth.The second scattering coefficient μ in this step tlagfrom first acquisition time (t 1) start, to the time period that last acquisition time terminates (i.e. t 1to t ibetween) in not calculate in the same time.Every two adjacent calculating second scattering coefficient μ tlaginterval time in moment identical.With t 1for initial time (identical with the initial time of previous step), Δ t are interval time, so, corresponding calculating second scattering coefficient μ tlagthe corresponding moment is followed successively by t 1, t 1+ Δ t, t 1+ 2* Δ t, t 1+ 3 Δ t ..., the second corresponding scattering coefficient μ tlag, be designated as μ successively tlag1, μ tlag2, μ tlag3, μ tlag4
S150: the dependency calculating the second scattering coefficient and blood glucose value.In this step, calculate the second scattering coefficient μ according to Pearson correlation tlagwith the dependency of blood glucose value.The following formula of Pearson correlation computing formula (2):
R = n Σ i = 1 n ( μ ti - μ tlag ) gl i - Σ i = 1 n ( μ ti - μ tlag ) Σ i = 1 n gl i n Σ i = 1 n ( μ ti - μ tlag ) 2 - ( Σ i = 1 n ( μ ti - μ tlag ) ) 2 n Σ i = 1 n gl i 2 - ( Σ i = 1 n gl i ) 2 - - - ( 2 )
Wherein, μ tiexpression sampling instant is t itime corresponding the first scattering coefficient, gl ifor t ithe blood glucose value in moment, μ tlagfor the second scattering coefficient calculated in step S140.Namely t is worked as lagfor t 1time, calculate correlation coefficient now, be designated as R 1; Work as t lagfor t 1during+Δ t, calculate correlation coefficient now, be designated as R 2; t lagfor t 1during+2* Δ t, calculate correlation coefficient now, be designated as R 3; , the like, obtain one group of correlation coefficient.
S160: dependency described in com-parison and analysis, draws blood glucose time delay.To each t that previous step draws lagcorresponding correlation coefficient (one group of correlation coefficient) carries out com-parison and analysis.On this skin depth region, the blood glucose t of time delay corresponding to correlation coefficient maximum in one group of correlation coefficient lag.
After the skin different depth blood glucose time delay obtained by the method for above-mentioned calculating blood glucose time delay, time delay can be selected short and the depth areas that detects in real time as blood glucose of the depth areas that dependency is high, to reach the object improving blood sugar test accuracy.Or directly selective scattering coefficient and the highest depth areas of blood glucose dependency can carry out follow-up blood glucose prediction, compensate calculating time delay by the blood glucose calculated in this region, to obtain overall blood glucose prediction trendgram accurately.
Due to the individual difference of people and the randomness of skin interior tissue species distribution, the blood glucose of the skin different depth obtained by detection computations is only applicable to individual and can not be general time delay.Namely, before different people carries out optics noninvasive dynamics monitoring, the calculating carrying out blood glucose time delay is first all needed.
Please 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 positioned at the papillary layer of the skin corium first half, and is negative correlation with change of blood sugar.As can be seen from Figure 4, there is certain retardance in scattering coefficient and blood sugar for human body, and time delay is 17 minutes.
The method of above-mentioned calculating blood glucose time delay can calculate the time delay of blood glucose on skin different depth region accurately, after learning this blood glucose time delay, then carries out blood sugar test, can reach the object improving blood sugar test accuracy.
Please refer to Fig. 5, additionally provide a kind of system 100 calculating blood glucose time delay in one embodiment.The system 100 of this calculating blood glucose time delay comprises the degree of depth and selects module 110, blood glucose value acquisition module 120, scattering coefficient computing module 130, scattering coefficient comparing module 140, Calculation of correlation factor module 150 and com-parison and analysis module 160.The degree of depth selects the dependency of module 110 blood glucose value and scattering coefficient on analytical calculation different depth, and selects scattering coefficient and the maximum depth areas of blood glucose value dependency; Blood glucose value acquisition module 120 is for obtaining not blood glucose value in the same time in the same degree of depth of sample and corresponding acquisition time; Scattering coefficient computing module 130 is for calculating the first scattering coefficient of each acquisition time of correspondence in the described same degree of depth; Scattering coefficient comparing module 140 for calculating the second scattering coefficient in the described same degree of depth, described second scattering coefficient μ tlagfrom first acquisition time, not calculating in the same time within the time period that last acquisition time terminates; Calculation of correlation factor module 150 is for calculating the second scattering coefficient μ tlagwith the dependency of blood glucose value; Com-parison and analysis module 160, for dependency described in com-parison and analysis, draws blood glucose time delay.
The degree of depth selects the dependency of module 110 blood glucose value and second scattering coefficient on analytical calculation different depth, and selects scattering coefficient and the maximum depth areas of blood glucose value dependency.First need to be averaged there being OCT machine to scan all A-scan in the 3-D view obtained, thus OCT one dimension depth signal figure (i.e. Fig. 3) can be obtained.From skin surface, successively selected depth be 1*w, 2*w ..., n*w (n=1,2,3 ...) detect depth areas.In each depth areas, collect each moment t iblood glucose value gl i, and calculate each moment t iscattering coefficient μ ti(i be greater than 0 integer).Again according to following formula (1), the correlation coefficient of scattering coefficient and blood glucose value on different depth region can be obtained.
R = n Σ i = 1 n μ ti gl i - Σ i = 1 n μ ti Σ i = 1 n gl i n Σ i = 1 n μ ti 2 - ( Σ i = 1 n μ ti ) 2 n Σ i = 1 n gl i 2 - ( Σ i = 1 n gl i ) 2 - - - ( 1 )
Wherein, μ tifor scattering coefficient, gl ifor blood glucose value.The correlation coefficient of scattering coefficient and blood glucose value on the different depth region obtained according to above-mentioned formulae discovery is compared, thus selects detection depth areas.Concrete, the standard of selection can be self-defined, such as, correlation coefficient can be selected maximum, and namely corresponding when correlation coefficient is maximum depth areas is as the detection depth areas of subsequent calculations blood glucose time delay; As long as the detection depth areas of coefficient R depth areas corresponding in certain numerical intervals as subsequent calculations blood glucose time delay also can be defined.
Blood glucose value acquisition module 120 is for obtaining not blood glucose value in the same time in the same degree of depth of sample and corresponding acquisition time.At t 1, t 2..., t iin the moment, the blood glucose value gathering the corresponding moment is respectively designated as gl 1, gl 2..., gl i, and write down acquisition time, wherein, i be greater than 0 integer.T 1for obtaining blood glucose value gl ifirst acquisition time.
Scattering coefficient computing module 130 is for calculating the first scattering coefficient μ of each acquisition time of correspondence in the described same degree of depth ti.Correspond to each acquisition time gathering blood glucose value, calculate the first scattering coefficient μ that each acquisition time corresponding is corresponding ti.Namely at t 1, t 2..., t iin the moment, the first corresponding respectively scattering coefficient is μ t1, μ t2..., μ ti.
Scattering coefficient comparing module 140 is for calculating the second scattering coefficient in the described same degree of depth, described second scattering coefficient is from first acquisition time, not calculating in the same time within the time period that last acquisition time terminates.Second scattering coefficient μ tlagfrom first acquisition time (t 1) start, to the time period that last acquisition time terminates (i.e. t 1to t ibetween) in not calculate in the same time.Every two adjacent calculating second scattering coefficient μ tlaginterval time in moment identical.With t 1for initial time, Δt is interval time, so, corresponding calculating second scattering coefficient μ tlagthe corresponding moment is followed successively by t 1, t 1+ Δt, t 1+ 2* Δ t, t 1+ 3 Δ t ..., the second corresponding scattering coefficient μ tlag, be designated as μ successively tlag1, μ tlag2, μ tlag3, μ tlag4
Calculation of correlation factor module 150 is for calculating the dependency of the second scattering coefficient and blood glucose value.The second scattering coefficient μ is calculated according to Pearson correlation tlagwith the dependency of blood glucose value.The following formula of Pearson correlation computing formula (2):
R = n Σ i n ( μ ti - μ tlag ) gl i - Σ i n ( μ ti - μ tlag ) Σ i n gl i n Σ i n ( μ ti - μ tlag ) 2 - ( Σ i n ( μ ti - μ tlag ) ) 2 n Σ i n gl i 2 - ( Σ i n gl i ) 2 - - - ( 2 )
Wherein, μ tiexpression sampling instant is t itime corresponding the first scattering coefficient, gl ifor t ithe blood glucose value in moment, μ tlagfor the second scattering coefficient that scattering coefficient comparing module 140 calculates.Namely t is worked as lagfor t 1time, calculate correlation coefficient now, be designated as R 1; Work as t lagfor t 1during+Δ t, calculate correlation coefficient now, be designated as R 2; t lagfor t 1during+2* Δ t, calculate correlation coefficient now, be designated as R 3; , the like, obtain one group of correlation coefficient.
Com-parison and analysis module 160, for dependency described in com-parison and analysis, draws blood glucose time delay.Calculation of correlation factor module 150 is shown that each correlation coefficient carries out com-parison and analysis.On this skin depth region, the blood glucose t of time delay corresponding to correlation coefficient maximum in one group of correlation coefficient lag.
The system 100 of above-mentioned calculating blood glucose time delay can calculate the time delay of blood glucose on skin different depth region exactly, can accurately select the depth areas being suitable as blood sugar test region according to this time delay, thus the accuracy of blood sugar test can be improved.
The above embodiment only have expressed several embodiment of the present invention, and it describes comparatively concrete and detailed, but therefore can not be interpreted as the restriction to the scope of the claims of the present invention.It should be pointed out that for the person of ordinary skill of the art, without departing from the inventive concept of the premise, can also make some distortion and improvement, these all belong to protection scope of the present invention.Therefore, the protection domain of patent of the present invention should be as the criterion with claims.

Claims (10)

1. calculate the method for blood glucose time delay, it is characterized in that, comprise the steps:
Obtain the not blood glucose value in the same time in the sample same degree of depth and corresponding acquisition time;
Calculate the first scattering coefficient of each acquisition time of correspondence in the described same degree of depth;
Calculate the second scattering coefficient in the described same degree of depth, described scattering coefficient is from first acquisition time, not calculating in the same time within the time period that last acquisition time terminates;
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 for calculating blood glucose according to claim 1 time delay, it is characterized in that, the step of the dependency of described calculating second scattering coefficient and blood glucose value is:
The dependency of the second scattering coefficient and blood glucose value is calculated according to Pearson correlation.
3. the method for calculating blood glucose according to claim 1 time delay, is characterized in that, dependency described in described com-parison and analysis, show that the step of blood glucose time delay is:
Blood glucose to equal time delay correlation coefficient maximum time corresponding time.
4. the method for calculating blood glucose according to claim 1 time delay, it is characterized in that, the step of the second scattering coefficient in the described same degree of depth of described calculating is:
The interval time in the moment of every two adjacent described second scattering coefficients of calculating is identical.
5. the method for calculating blood glucose according to claim 1 time delay, is characterized in that, also comprises the steps: before the not blood glucose value step in the same time in the same degree of depth of described acquisition sample
The dependency of blood glucose value and scattering coefficient on analytical calculation different depth, and select scattering coefficient and the maximum depth areas of blood glucose value dependency.
6. calculate the system of blood glucose time delay, it is characterized in that, comprise
Blood glucose value acquisition module, for obtaining the not blood glucose value in the same time in the same degree of depth of sample, and obtains described acquisition time;
Scattering coefficient computing module, for calculating the first scattering coefficient of each acquisition time of correspondence in the described same degree of depth;
Scattering coefficient comparing module, for calculating the second scattering coefficient in the described same degree of depth, described second scattering coefficient is from first acquisition time, not calculating in the same time within the time period that last acquisition time terminates;
Calculation of correlation factor 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 of calculating blood glucose according to claim 6 time delay, it is characterized in that, described Calculation of correlation factor module calculates the dependency of the second scattering coefficient and blood glucose value according to Pearson correlation.
8. the system of calculating blood glucose according to claim 6 time delay, is characterized in that, in described com-parison and analysis module, blood glucose to equal time delay correlation coefficient maximum time corresponding time.
9. the system of calculating blood glucose according to claim 6 time delay, is characterized in that, in described scattering coefficient comparing module, the interval time in the moment of every two adjacent described second scattering coefficients of calculating is identical.
10. the system of calculating blood glucose according to claim 6 time delay, is characterized in that, the system of described calculating blood glucose time delay also comprises:
The degree of depth selects module, the dependency of blood glucose value and scattering coefficient on analytical calculation different depth, and selects scattering coefficient and the maximum depth areas of blood glucose value dependency.
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 true CN104545808A (en) 2015-04-29
CN104545808B 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)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100161269A1 (en) * 2003-08-01 2010-06-24 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

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100161269A1 (en) * 2003-08-01 2010-06-24 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

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
ANANDA BASU ET.AL.: "Time Lag of Glucose From Intravascular to Interstitial Compartment in Humans", 《DIABETES》 *
苏亚等: "光学相干层析无创血糖检测中相关性分析及标定", 《中国激光》 *

Also Published As

Publication number Publication date
CN104545808B (en) 2017-01-18

Similar Documents

Publication Publication Date Title
Blank et al. Clinical results from a noninvasive blood glucose monitor
von Lilienfeld-Toal et al. A novel approach to non-invasive glucose measurement by mid-infrared spectroscopy: The combination of quantum cascade lasers (QCL) and photoacoustic detection
CN101557752B (en) Method for the glucose concentration in pulsational blood
US6456870B1 (en) Non-invasive method of determining skin thickness and characterizing layers of skin tissue in vivo
CN104224196B (en) The method of non-invasive measurement blood component concentration
CN103149177B (en) Device and method for detecting biological tissue of pressure modulation near infrared spectrum
CN104958075A (en) Non-invasive measurement method of skin thickness and blood sugar concentration using Raman spectrum, and calibration method
US20150342508A1 (en) Non-invasive optical measurement of blood analyte
CN103948393B (en) A kind of near infrared no-wound detection method of blood constituent content and device
JP5674093B2 (en) Concentration determination apparatus, concentration determination method, and program
US20100160754A1 (en) Method and Apparatus for Quantification of Optical Properties of Superficial Volumes Using Small Source-to-Detector Separations
Li et al. Classification of diabetes and measurement of blood glucose concentration noninvasively using near infrared spectroscopy
CN104188664B (en) Blood sugar test scaling method and system
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
US20010041829A1 (en) Non-invasive method of determining skin thickness and characterizing layers of skin tissue in vivo
US20120116236A1 (en) Hydration and blood flow adjusted glucose measurement
CN104545808B (en) Method and system for calculating delay time of blood sugar
Colombo et al. Coherent fluctuations in time-domain diffuse optics
CN106963392B (en) Optical noninvasive blood glucose detection two-dimensional correlation calibration method and system
JP5626880B2 (en) Concentration determination apparatus, concentration determination method, and program
Arabi et al. Early Diagnosis of Diabetic Vasculopathy by Thermo-Regulatory Vascular Impairment
CN103884556B (en) One freezes Tilapia Fillet quality determination method

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
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.