CN104188664B - Blood sugar test scaling method and system - Google Patents

Abstract

The invention discloses a kind of blood sugar test scaling method and system, this blood sugar test scaling method and system, the correlation coefficient of skin different depth tissue scatter's coefficient and temperature can be calculated, by finding the analysis of correlation coefficient, all responsive tissue regions is changed to temperature and blood sugar for human body.The high relevant range utilizing correlation analysis to obtain can as the demarcation region of noninvasive dynamics monitoring, for aspects such as follow-up blood sugar test.Utilize this blood sugar test scaling method, greatly simplify Woundless blood sugar and demarcate flow process, substantially reduce the nominal time, alleviate the misery of diabetics to a certain extent, reduce infection risk.

Description

Blood sugar test scaling method and system

Technical field

The present invention relates to blood sugar test field, particularly relate to blood sugar test scaling method and system.

Background technology

Diabetes are frequently-occurring diseases of mid-aged population, and along with the raising of people's living standard, the sickness rate of diabetes also rises day by day, and diabetes, tumor are classified as worldwide three disaster diseases by World Health Organization (WHO) together with cardiovascular and cerebrovascular disease.Find a kind of method that blood sugar concentration detects the prevention and therapy of diabetes is had very great significance.Current blood sugar concentration detection method is most widely used is have wound to measure.The method having wound to detect main application relies on electrochemical method to refer to that blood detects the blood glucose concentration value obtaining patient to patient.The method can realize detecting the blood glucose concentration value in a certain moment, but equally also there are some shortcomings.Such as extract and refer to that blood is more painful, measurement blood sugar concentration needs consumptive material etc.Also have Wicresoft's detection method in addition, Wicresoft's detection method mainly detects blood sugar concentration by detecting the tissue fluid extracted from skin, and the method can alleviate the misery of patient, but causes certain wound to patient equally.In sum, the device of Non-invasive detection blood sugar concentration and correlation method are highly significant.

Research shows, utilizes optical means can detect the corresponding relation of the intensity of the light of returning from skin reflex and the scattering coefficient of skin histology, and the inner blood glucose concentration value of the scattering coefficient of skin histology and biological tissue is closely related.Following single scattered light strength formula can the corresponding relation of the approximate description intensity of light of returning from skin reflex and the scattering coefficient of skin histology:

I _R＝I _oexp[-(μ _a+μ _s)L]

Wherein I _rfor the light intensity of returning from skin reflex, I _ofor projecting the light intensity of skin, μ _afor the absorptance of skin histology, μ _sfor the scattering coefficient of skin histology, L is the total optical path of light transdermal.As can be seen from the above equation, the intensity of the light of returning from skin reflex is the exponential damping with scattering coefficient and absorptance.

In skin histology, the refractive index of body fluid and the refractive index of organelle there are differences.This difference can cause skin histology to scattering of light phenomenon.Glucose is the main ingredient of of body fluid, and when blood sugar concentration changes, the refractive index of body fluid also can change thereupon, and this can cause tissue scatter's coefficient to change.At near infrared band, the scattering coefficient change that glucose causes wants the change of specific absorptivity much bigger, and therefore the change of blood sugar concentration mainly causes the change of skin histology scattering coefficient instead of the change of absorptance.As can be seen here, optical detection skin histology scattering coefficient can as of a noninvasive dynamics monitoring important means.

Skin histology is from organizational structure, comprise a lot of sweat glands, oils and fats gland and blood vessel, these organizational structuries are very strong for the absorption of infrared light, therefore at use OCT (OpticalCoherenceTomography, optical coherence tomography) when gathering relevant to skin histology data, can very weak or disappearance at these organizational structure regional signals.When making optically to carry out blood glucose concentration value detection, need to consider that these strong absorptive tissue parts are on the impact on blood sugar test precision of whole signal.Therefore, be necessary to find the rule that in skin histology, different tissues structural region changes by blood glucose concentration value, thus strong absorptive tissue structure is avoided when Non-invasive detection blood glucose, obtain the measurement data with the highly sensitive tissue of change of blood sugar, realize the method for being penetrated the change of skin detection blood sugar concentration by illumination.

OCT is continue ultra sonic imaging, X ray CT (ComputedTomography, computed tomography), MRI (MagneticResonanceImaging, nuclear magnetic resonance) after biomedical imaging technology of new generation, low coherence interference technology, the product that combines of confocal microscope principle and superhet Detection Techniques, non-intruding can be realized, high sensitivity, high-resolution carry out imaging to tissue.The ultimate principle that it utilizes low-coherent light to interfere, accurately can be measured amplitude and the relative phase of reflected light to the backscatter signals of the low-coherent light of incidence by detection biological tissue different depth aspect, obtain the microstructure features on organization internal depth direction, pass through horizontal scanning again, data and the image of biological tissue's two dimension or three dimensional structure can be obtained.By scanning the OCT two dimension that obtains or three-dimensional data can use single scattered light strength formula to calculate scattering coefficient with change in depth.

In organism optical, generally skin histology is regarded as the model of a stratiform, simply skin histology can be divided into three layers: skin surface, prickle cell layer and skin corium.Think that the composition structure of skin in layer is identical simultaneously, therefore can be described the propagation in skin histology by light propagation in the depth direction, namely can such as, by calculating the optical parametric of skin histology in the degree of depth, scattering coefficient.There are some researches prove, when the blood glucose concentration value in human body changes, the scattering coefficient of skin histology can change.Therefore by using the minor variations of OCT detection of skin organization internal scattering coefficient, the functional relationship between skin scattering-in coefficient and blood glucose concentration value can be set up, and utilizes this functional relationship to realize the detection of blood glucose concentration value.But because skin internal structure is very complicated, the sensitivity that skin histology changes blood glucose concentration value on different depth is different.Therefore this factor is considered, the functional relationship set up between skin scattering-in coefficient and blood glucose concentration value needs to find and the scattering coefficient on the maximally related skin depth position of blood glucose concentration value, and find the process with the scattering coefficient on the maximally related skin depth position of blood glucose concentration value, be commonly referred to as blood sugar test and demarcate.

At present, the method that blood sugar test is demarcated refers to that the blood glucose value of blood or venous blood and synchronization tissue scatter coefficient set up corresponding relation by gathering several times, the relevant range obtaining skin different depth by correlation calculations carries out follow-up blood sugar test but this kind of scaling method required time is long, generally need 2-3 hour, and need repeatedly there is wound blood sampling, not only increase the misery of testee, also increase infection risk.

Summary of the invention

Based on this, be necessary for the problems referred to above, a kind of blood sugar test scaling method and system be provided, this blood sugar test scaling method and the system calibrating time short, and blood sampling number of times few, reduce misery and the infection risk of testee.

A kind of blood sugar test scaling method, comprises the steps: to determine surveyed area: carry out temperature modulation to described surveyed area, and obtains temperature value and the scattering coefficient of described surveyed area; The correlation coefficient of temperature value described in analytical calculation and scattering coefficient, obtains demarcating region; At least two not blood glucose values are in the same time gathered in temperature modulation process; According to demarcation region, calculate the functional relationship of scattering coefficient and corresponding blood glucose value.

Wherein in an embodiment, described temperature modulation is carried out to described surveyed area, and the step of the temperature value and scattering coefficient that obtain described surveyed area comprises: default temperature modulation scope is equally divided into several probe temperatures; The temperature of modulating described surveyed area successively equals described probe temperature, when the temperature of surveyed area equals described probe temperature, obtains scattering coefficient corresponding to the temperature value of surveyed area and skin different depth region thereof.

Wherein in an embodiment, described temperature modulation is carried out to described surveyed area, and obtain the step obtaining described scattering coefficient in the temperature value of described surveyed area and the step of scattering coefficient and be: certain 1 p in setting skin _i1for starting point, the thickness of unit area is Interval; Calculate and obtain with p _i1for starting point, thickness be the integral multiple of Interval depth areas in one group of scattering coefficient be μ _i12, μ _i13..., μ _i1n.

Wherein in an embodiment, the correlation coefficient of temperature value and scattering coefficient described in described analytical calculation, the step obtaining demarcating region is: adopt Pearson's product moment correlation method to calculate the coefficient R of temperature value and scattering coefficient,

$R=\frac{n\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i{T}_i-\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{T}_i}{\sqrt{n\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i^2-{\left(\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i\right)}^2}\sqrt{n\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{T_i}^2-{\left(\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{T}_i\right)}^2}}$

Wherein, variable μ is scattering coefficient, and T is temperature value.

Wherein in an embodiment, the correlation coefficient of temperature value and scattering coefficient described in described analytical calculation, the step obtaining demarcating region is: the absolute value of described correlation coefficient is more close to 1, the linear relationship of scattering coefficient and temperature value is stronger, and the skin depth region corresponding to scattering coefficient is demarcates region.

A kind of blood sugar test calibration system, comprising: surveyed area determination module, for determining surveyed area; Temperature modulation module, for carrying out temperature modulation to described surveyed area, and obtains temperature value and the scattering coefficient of described surveyed area; Correlating module, for the correlation coefficient of temperature value described in analytical calculation and scattering coefficient, obtains demarcating region; Blood glucose acquisition module, for carrying out the blood glucose value gathering at least two moment in temperature modulation process; Functional relationship computing module, for according to demarcation region, calculates the functional relationship of scattering coefficient and corresponding blood glucose value.

Wherein in an embodiment, described temperature modulation module comprises: division unit, for default temperature modulation scope is equally divided into several probe temperatures; Modulating unit, the temperature for modulating described surveyed area successively equals described probe temperature, when the temperature of surveyed area equals described probe temperature, obtains scattering coefficient corresponding to the temperature value of surveyed area and skin different depth region thereof.

Wherein in an embodiment, described modulating unit comprises: setup unit, for setting certain 1 p in skin _i1for starting point, the thickness of unit area is Interval; Computing unit, obtains for calculating with p _i1for starting point, thickness be the integral multiple of Interval depth areas in one group of scattering coefficient be μ _i12, μ _i13..., μ _i1n.

Wherein in an embodiment, described correlating module adopts Pearson's product moment correlation method to calculate the coefficient R of temperature value and scattering coefficient,

$R=\frac{n\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i{T}_i-\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{T}_i}{\sqrt{n\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i^2-{\left(\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i\right)}^2}\sqrt{n\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{T_i}^2-{\left(\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{T}_i\right)}^2}}$

Wherein, variable μ is scattering coefficient, and T is temperature value.

Wherein in an embodiment, in described correlating module, the absolute value of correlation coefficient is more close to 1, and the linear relationship of scattering coefficient and temperature value is stronger, and the skin depth region corresponding to scattering coefficient is demarcates region.

Above-mentioned blood sugar test scaling method and system, can calculate the correlation coefficient of skin different depth tissue scatter's coefficient and temperature, changes all responsive tissue regions by finding the analysis of correlation coefficient to temperature and blood sugar for human body.The high relevant range utilizing correlation analysis to obtain can as the demarcation region of noninvasive dynamics monitoring, for aspects such as follow-up blood sugar test.Utilize this blood sugar test scaling method, greatly simplify Woundless blood sugar and demarcate flow process, substantially reduce the nominal time, alleviate the misery of diabetics to a certain extent, reduce infection risk.

Accompanying drawing explanation

Fig. 1 is the blood sugar test scaling method flow chart in an embodiment;

Fig. 2 carries out temperature modulation to described surveyed area in an embodiment, and obtains temperature value, the distance degree of depth of skin surface and the flow chart of scattering coefficient step of described surveyed area simultaneously;

Fig. 3 is the skin normalization one dimension surface of intensity distribution;

Fig. 4 is blood sugar test calibration system Organization Chart in an embodiment;

Fig. 5 is the Organization Chart of temperature modulation module in an embodiment;

Fig. 6 is the dependency demarcate analysis figure of hands skin of forearm different depth regions scatter coefficient and temperature value.

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.

Before carrying out optics noninvasive dynamics monitoring, all need to demarcate separately, and nominal data is only applicable to individual and can not be general, this is because the individual difference of people and the randomness of skin interior tissue species distribution determine.

Please refer to Fig. 1, in one embodiment, a kind of blood sugar test scaling method is provided, comprises the steps:

S110: determine surveyed area.

This step determines surveyed area on the person to be detected, generally can on the medial surface of hands forearm, select one piece of region to be fixed as surveyed area.

S120: temperature modulation is carried out to described surveyed area, and the temperature value and the scattering coefficient that obtain described surveyed area.Please refer to Fig. 2, also comprise the steps: in this step

S121: default temperature modulation scope is equally divided into several probe temperatures.

When carrying out temperature modulation to surveyed area, first can determine a default temperature modulation scope, and this temperature modulation scope preset mainly is equally divided into several probe temperatures by this step.In the present embodiment, the temperature modulation scope preset is 30 DEG C-40 DEG C, and this temperature modulation scope is equally divided into 11 probe temperatures, is followed successively by 30 DEG C, 31 DEG C, 32 DEG C, 33 DEG C, 34 DEG C, 35 DEG C, 36 DEG C, 37 DEG C, 38 DEG C, 39 DEG C, 40 DEG C.

S122: the temperature of modulating described surveyed area successively equals described probe temperature, when the temperature of surveyed area equals described probe temperature, obtains scattering coefficient corresponding to the temperature value of surveyed area and skin different depth region thereof.

It is first 30 DEG C by the temperature modulation of surveyed area.When the temperature of surveyed area reaches 30 DEG C, obtain the scattering coefficient of skin different depth within the scope of the temperature value of surveyed area now, surveyed area.Afterwards being 31 DEG C by the temperature modulation of surveyed area, when the temperature of surveyed area reaches 31 DEG C, obtain the scattering coefficient of skin different depth within the scope of the temperature value of surveyed area now, surveyed area.Circulation like this, till being 40 DEG C by the temperature modulation of surveyed area.And the probe temperature circulated successively within the scope of complete temperature modulation only needs the time of about 15 minutes, and need compared with 2-3 hour with traditional scaling method, detection time substantially reduces.

S130: the correlation coefficient of temperature value described in analytical calculation and scattering coefficient, obtains demarcating region.

This step mainly adopts Pearson's product moment correlation method (Personproduct-momentcorrelation) to calculate the coefficient R of scattering coefficient and temperature.

$R=\frac{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i{T}_i-n\stackrel{\‾}{\mathrm{\μ}}\stackrel{\‾}{T}}{(n-1)s_{\mathrm{\μ}}{s}_T}=\frac{n\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i{T}_i-\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{T}_i}{\sqrt{n\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i^2-{\left(\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i\right)}^2}\sqrt{n\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{T_i}^2-{\left(\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{T}_i\right)}^2}}$

Wherein, variable μ is scattering coefficient, and T is temperature value.The excursion of coefficient R is-1 to 1, and when the linear relationship of Two Variables strengthens, coefficient R is tending towards 1 or-1.

Please refer to Fig. 3, skin surface is maximum 1, A-B section is epidermal area region, and B-C section is the juncture area of epidermal area and skin corium, and C point is following is skin corium region.Concrete analytical calculation process is as follows:

Certain 1 p in (a) setting skin _i1for starting point, the thickness of unit area is Interval, along skin depth direction (abscissa) respectively calculated thickness be Interval × 1, Interval × 2 ..., Interval × n region inner tissue scattering coefficient, namely calculate p _i1p _i2region (i.e. p _i1to p _i2between region) in scattering coefficient μ _i12, p _i1p _i3region (i.e. p _i1to p _i3between region) in scattering coefficient μ _i13..., p _i1p _inregion (i.e. p _i1to p _inbetween region) in scattering coefficient μ _i1n.Finally can obtain with p _i1for starting point, thickness is one group of scattering coefficient μ in the region of the integral multiple of Interval _i12, μ _i13..., μ _i1n.

B () sets the temperature value (namely corresponding multiple probe temperatures) that records as T ₁, T ₂..., T _m, the 3-D view in the correspondence blood sampling moment obtained by OCT measurement is s ₁, s ₂..., s _m.Then correspond to p in every width image _i1p _i2the correlation coefficient that region can calculate temperature and scattering coefficient is R _i12, corresponding p _i1p _i3the coefficient R in region _i13..., p _i1p _inthe coefficient R in region _i1n.Namely depth areas p is corresponded to _i1p _i2, p _i1p _i3..., p _i1p _in, finally can obtain one group of corresponding coefficient R _i12, R _i13..., R _i1n.

C (), according to the displacement Offset of setting, starting point is by p _i1move to p _{(i+1) 1}.According to the computational methods of above-mentioned steps (a) (b), can obtain corresponding to depth areas is p _{(i+1) 1}p _{(i+1) 2}, p _{(i+1) 1}p _{(i+1) 3}..., p _{(i+1) 1}p _{(i+1) n}one group of coefficient R _{(i+1) 12}, R _{(i+1) 13}..., R _{(i+1) 1n}.

Be then interval by displacement Offset, with p _i1, p _{(i+1) 1}..., p _{(i+m) 1}for calculating initial position point, the correlation coefficient of m group in skin different depth region can be obtained by above-mentioned computational methods.By each group of coefficient R on different depth region _{(i+1) 12}, R _{(i+1) 13}..., R _{(i+1) 1n}the correlation matrix be constructed as follows:

$\left[\begin{array}{cccc}{R}_{i12}& R_{i13}& \·\·\·& R_{i1n}\\ R_{(i+1)12}& R_{(i+1)13}& \·\·\·& R_{(i+1)1n}\\ \·& \·& & \·\\ \·& \·& \·\·\·& \·\\ \·& \·& & \·\\ R_{(i+m)12}& R_{(i+m)13}& \·\·\·& R_{(i+m)1n}\end{array}\right]$

The excursion of correlation coefficient is-1 to 1, and the absolute value of correlation coefficient is more large more close to 1, then illustrate that the linear relationship of scattering coefficient and temperature value is stronger.When determining to demarcate region, the corresponding skin depth region of the maximum absolute value of correlation coefficient (being scattering coefficient the strongest with the linear relationship of temperature) can be selected for demarcate region.Can also preset a certain reduced value (such as setting reduced value is 0.8), the absolute value choosing all correlation coefficienies is greater than the skin depth region of the correspondence of this reduced value as demarcation region.

S140: gather at least two not blood glucose values in the same time in temperature modulation process.

Although temperature is different with the mechanism of action of blood glucose to tissue scattering coefficient, closely similar in the relevant range distribution in skin depth direction.Its reason is the distribution situation of the inner different tissues material of skin due to the essence reflection of: relevant range.Scattering coefficient, for the suck tissue of variations in temperature sensitivity, also has hypersensitivity for change of blood sugar equally.Therefore, can be demarcated for blood glucose by the demarcation region of the scattering coefficient that calculates and temperature and detect.

It is the process temperature of surveyed area being set to each probe temperature when temperature modulation is carried out to surveyed area, in the process, two moment (when the temperature by surveyed area is arranged to certain two probe temperature) at least will be selected to detect the blood glucose value of detected personnel.Generally can select to detect the blood glucose value when temperature of surveyed area being set to first probe temperature and last probe temperature.Concrete, conventional blood sugar detection method can be adopted to obtain blood glucose value, such as, gather the detection method referring to blood or venous blood.Compared with traditional scaling method, blood sampling number of times obviously reduces, thus reduces the misery of tested personnel and the risk of infection.

S150: according to demarcation region, calculate the functional relationship of scattering coefficient and corresponding blood glucose value.

Above-mentioned steps has been determined to demarcate region, namely has extremely sensitive depth areas to the change of temperature.Equally, this demarcation region change to blood glucose value also has hypersensitivity.By linear fit method, calculate the functional relationship of demarcating the blood glucose value that region records in corresponding scattering coefficient and step S140, this functional relationship obtained can be used for follow-up noninvasive dynamics monitoring.Concrete, linear fit method is as method of least square or partial least square method method.

Please refer to Fig. 6, calibrated and calculated scope is that skin surface is within the scope of subcutaneous 1 millimeter.Vertical coordinate represents the start position demarcating region, displacement Offset=0.1 millimeter, abscissa is expressed as the depth bounds of calculated relevant range, here the initial gap Interval=0.1 millimeter set, namely from start position, along skin depth direction, calculated thickness is 0.1 millimeter (Interval × 1), 0.2 millimeter (Interval × 2), 0.3 millimeter (Interval × 3) respectively ... scattering coefficient in region and the dependency demarcating blood glucose value, terminate until be calculated to when the degree of depth is the final position of 1 millimeter.Utilize the Matlab program of writing to carry out correlation analysis computing, in the dependency demarcate analysis figure calculated, a-quadrant represents positive correlation region (namely demarcating region).

Above-mentioned blood sugar test scaling method can calculate the correlation coefficient of skin different depth tissue scatter's coefficient and temperature, changes all responsive tissue regions by finding the analysis of correlation coefficient to temperature and blood sugar for human body.The high relevant range utilizing correlation analysis to obtain can as the demarcation region of noninvasive dynamics monitoring, for aspects such as follow-up blood sugar test.Utilize this blood sugar test scaling method, greatly simplify Woundless blood sugar and demarcate flow process, substantially reduce the nominal time, alleviate the misery of diabetics to a certain extent.

Please refer to Fig. 4, also provide a kind of blood sugar test calibration system 100 in one embodiment, this blood sugar test calibration system comprises surveyed area determination module 110, temperature modulation module 120, correlating module 130, blood glucose acquisition module 140 and functional relationship computing module 150.Surveyed area determination module 110 is for determining surveyed area, temperature modulation module 120 is for carrying out temperature modulation to surveyed area, and obtain temperature value and the scattering coefficient of surveyed area, correlating module 130 is for the correlation coefficient of temperature value described in analytical calculation and scattering coefficient, obtain demarcating region, blood glucose acquisition module 140 is for carrying out the blood glucose value gathering at least two moment in temperature modulation process, functional relationship computing module 150, for according to demarcation region, calculates the functional relationship of scattering coefficient and corresponding blood glucose value.

Please refer to Fig. 5, temperature modulation module 120 comprises division unit 121 and modulating unit 122, and division unit 121 is for being equally divided into several probe temperatures by default temperature modulation scope.When carrying out temperature modulation to surveyed area, first can determine a default temperature modulation scope, be generally that this temperature modulation scope preset is equally divided into several probe temperatures.In the present embodiment, the temperature modulation scope preset is 30 DEG C-40 DEG C, and this temperature modulation scope is equally divided into 11 probe temperatures, is followed successively by 30 DEG C, 31 DEG C, 32 DEG C, 33 DEG C, 34 DEG C, 35 DEG C, 36 DEG C, 37 DEG C, 38 DEG C, 39 DEG C, 40 DEG C.

Modulating unit 122 equals probe temperature for the temperature in modulation detection region successively, when the temperature of surveyed area equals probe temperature, obtains temperature value and the scattering coefficient of surveyed area.It is first 30 DEG C by the temperature modulation of surveyed area.When the temperature of surveyed area reaches 30 DEG C, obtain temperature value and the scattering coefficient of surveyed area.Be 31 DEG C by the temperature modulation of surveyed area more afterwards, when the temperature of surveyed area reaches 31 DEG C, obtain temperature value and the scattering coefficient of surveyed area.Circulation like this, till being 40 DEG C by the temperature modulation of surveyed area.And the probe temperature circulated successively within the scope of complete temperature modulation only needs the time of about 15 minutes, need compared with 2-3 hour with traditional scaling method, detection time substantially reduces.

Correlating module 130, for the correlation coefficient of temperature value described in analytical calculation and scattering coefficient, obtains demarcating region.The main coefficient R adopting Pearson's product moment correlation method (Personproduct-momentcorrelation) to calculate scattering coefficient and temperature in this module.

$R=\frac{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i{T}_i-n\stackrel{\‾}{\mathrm{\μ}}\stackrel{\‾}{T}}{(n-1)s_{\mathrm{\μ}}{s}_T}=\frac{n\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i{T}_i-\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{T}_i}{\sqrt{n\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i^2-{\left(\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{\μ}}_i\right)}^2}\sqrt{n\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{T_i}^2-{\left(\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{T}_i\right)}^2}}$

Wherein, variable μ is scattering coefficient, and T is temperature value.The excursion of coefficient R is-1 to 1, and when the linear relationship of Two Variables strengthens, coefficient R is tending towards 1 or-1.

Although temperature is different with the mechanism of action of blood glucose to tissue scattering coefficient, closely similar in the relevant range distribution in skin depth direction.Its reason is the distribution situation of the inner different tissues material of skin due to the essence reflection of, relevant range.Scattering coefficient, for the suck tissue of variations in temperature sensitivity, also has hypersensitivity for change of blood sugar equally.Therefore, can be demarcated for blood glucose by the high relevant range of the scattering coefficient that calculates and temperature and detect.

Modulating unit 122 comprises setup unit and computing unit, and setup unit is for setting certain 1 p in skin _i1for starting point, the thickness of unit area is Interval, and computing unit obtains for calculating with p _i1for starting point, thickness be the integral multiple of Interval depth areas in one group of scattering coefficient be μ _i12, μ _i13..., μ _i1n.Each group of coefficient R on the acquisition skin different depth region of its correspondence _{(i+1) 12}, R _{(i+1) 13}..., R _{(i+1) 1n}the correlation matrix formed and determine to demarcate the concrete steps method in region and identical in above-mentioned blood sugar test scaling method, repeats no more herein.

Blood glucose acquisition module 140 for gathering at least two not blood glucose values in the same time in temperature modulation process.It is the process temperature of surveyed area being set to each probe temperature when temperature modulation is carried out to surveyed area, in the process, two moment (when the temperature by surveyed area is arranged to certain two probe temperature) at least will be selected to detect the blood glucose value of detected personnel.Generally can select to detect the blood glucose value when temperature of surveyed area being set to first probe temperature and last probe temperature.Concrete, conventional blood sugar detection method can be adopted to obtain blood glucose value, such as, gather the detection method referring to blood or venous blood.Compared with traditional scaling method, blood sampling number of times obviously reduces, thus reduces the misery of tested personnel and the risk of infection.

Functional relationship computing module 150, for according to demarcation region, calculates the functional relationship of scattering coefficient and corresponding blood glucose value., according to the demarcation region that correlating module has been determined, namely to the change of temperature, there is extremely sensitive depth areas.Equally, this demarcation region change to blood glucose value also has hypersensitivity.By linear fit method, calculate the functional relationship of demarcating the blood glucose value that region records in corresponding scattering coefficient and blood glucose acquisition module, this functional relationship obtained can be used for follow-up noninvasive dynamics monitoring.Concrete, linear fit method is as method of least square or partial least square method method.

Above-mentioned blood sugar test calibration system can calculate the correlation coefficient of skin different depth tissue scatter's coefficient and temperature, changes all responsive tissue regions by finding the analysis of correlation coefficient to temperature and blood sugar for human body.The high relevant range utilizing correlation analysis to obtain can as the demarcation region of noninvasive dynamics monitoring, for aspects such as follow-up blood sugar test.Utilize this blood sugar test calibration system, greatly simplify Woundless blood sugar and demarcate flow process, substantially reduce the nominal time, alleviate the misery of diabetics to a certain extent.

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 (8)

1. a blood sugar test scaling method, is characterized in that, comprises the steps:

Determine surveyed area;

Temperature modulation is carried out to described surveyed area, and obtains temperature value and the scattering coefficient of described surveyed area;

The correlation coefficient of temperature value described in analytical calculation and scattering coefficient, obtains demarcating region;

At least two not blood glucose values are in the same time gathered in temperature modulation process;

According to demarcation region, calculate the functional relationship of scattering coefficient and corresponding blood glucose value.

2. blood sugar test scaling method according to claim 1, is characterized in that, describedly carries out temperature modulation to described surveyed area, and the step of the temperature value and scattering coefficient that obtain described surveyed area comprises:

Default temperature modulation scope is equally divided into several probe temperatures;

The temperature of modulating described surveyed area successively equals described probe temperature, when the temperature of surveyed area equals described probe temperature, obtains scattering coefficient corresponding to the temperature value of surveyed area and skin different depth region thereof.

3. blood sugar test scaling method according to claim 1, is characterized in that, describedly carries out temperature modulation to described surveyed area, and obtains the step obtaining described scattering coefficient in the temperature value of described surveyed area and the step of scattering coefficient and be:

Certain 1 p in setting skin _i1for starting point, the thickness of unit area is Interval;

Calculate and obtain with p _i1for starting point, thickness be the integral multiple of Interval depth areas in one group of scattering coefficient be μ _i12, μ _i13..., μ _i1n.

4. blood sugar test scaling method according to claim 3, is characterized in that, the correlation coefficient of temperature value and scattering coefficient described in described analytical calculation, and the step obtaining demarcating region is:

Pearson's product moment correlation method is adopted to calculate the coefficient R of temperature value and scattering coefficient,

$R=\frac{n\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\μ}}_i{T}_i-\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\μ}}_i\underset{i=1}{\overset{n}{\Σ}}{T}_i}{\sqrt{n\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\μ}}_i^2-{\left(\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\μ}}_i\right)}^2}\sqrt{n\underset{i=1}{\overset{n}{\Σ}}{T_i}^2-{\left(\underset{i=1}{\overset{n}{\Σ}}{T}_i\right)}^2}}$

Wherein, variable μ is scattering coefficient, and T is temperature value.

5. a blood sugar test calibration system, is characterized in that, comprising:

Surveyed area determination module, for determining surveyed area;

Temperature modulation module, for carrying out temperature modulation to described surveyed area, and obtains temperature value and the scattering coefficient of described surveyed area;

Correlating module, for the correlation coefficient of temperature value described in analytical calculation and scattering coefficient, obtains demarcating region;

Blood glucose acquisition module, for carrying out the blood glucose value gathering at least two moment in temperature modulation process;

Functional relationship computing module, for according to demarcation region, calculates the functional relationship of scattering coefficient and corresponding blood glucose value.

6. blood sugar test calibration system according to claim 5, is characterized in that, described temperature modulation module comprises:

Division unit, for being equally divided into several probe temperatures by default temperature modulation scope;

Modulating unit, the temperature for modulating described surveyed area successively equals described probe temperature, when the temperature of surveyed area equals described probe temperature, obtains scattering coefficient corresponding to the temperature value of surveyed area and skin different depth region thereof.

7. blood sugar test calibration system according to claim 6, is characterized in that, described modulating unit comprises:

Setup unit, for setting certain 1 p in skin _i1for starting point, the thickness of unit area is Interval;

Computing unit, obtains for calculating with p _i1for starting point, thickness be the integral multiple of Interval depth areas in one group of scattering coefficient be μ _i12, μ _i13..., μ _i1n.

8. blood sugar test calibration system according to claim 5, is characterized in that, described correlating module adopts Pearson's product moment correlation method to calculate the coefficient R of temperature value and scattering coefficient,

$R=\frac{n\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\μ}}_i{T}_i-\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\μ}}_i\underset{i=1}{\overset{n}{\Σ}}{T}_i}{\sqrt{n\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\μ}}_i^2-{\left(\underset{i=1}{\overset{n}{\Σ}}{\mathrm{\μ}}_i\right)}^2}\sqrt{n\underset{i=1}{\overset{n}{\Σ}}{T_i}^2-{\left(\underset{i=1}{\overset{n}{\Σ}}{T}_i\right)}^2}}$

Wherein, variable μ is scattering coefficient, and T is temperature value.