EP3408772A1 - Estimation d'une concentration d'une substance au moyen de plusieurs capteurs non invasifs - Google Patents

Estimation d'une concentration d'une substance au moyen de plusieurs capteurs non invasifs

Info

Publication number
EP3408772A1
EP3408772A1 EP17702061.7A EP17702061A EP3408772A1 EP 3408772 A1 EP3408772 A1 EP 3408772A1 EP 17702061 A EP17702061 A EP 17702061A EP 3408772 A1 EP3408772 A1 EP 3408772A1
Authority
EP
European Patent Office
Prior art keywords
concentration
model
chemical constituent
spectrum
sensor
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.)
Withdrawn
Application number
EP17702061.7A
Other languages
German (de)
English (en)
Inventor
Terje VEGAR KARSTANG
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.)
Prediktor Medical AS
Original Assignee
Prediktor Medical AS
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 Prediktor Medical AS filed Critical Prediktor Medical AS
Publication of EP3408772A1 publication Critical patent/EP3408772A1/fr
Withdrawn legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0075—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by spectroscopy, i.e. measuring spectra, e.g. Raman spectroscopy, infrared absorption spectroscopy
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/053—Measuring electrical impedance or conductance of a portion of the body
    • A61B5/0531—Measuring skin impedance
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/145—Measuring 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/14532—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue for measuring glucose, e.g. by tissue impedance measurement
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/145—Measuring 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/1455—Measuring 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2505/00—Evaluating, monitoring or diagnosing in the context of a particular type of medical care
    • A61B2505/07—Home care
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2560/00—Constructional details of operational features of apparatus; Accessories for medical measuring apparatus
    • A61B2560/02—Operational features
    • A61B2560/0223—Operational features of calibration, e.g. protocols for calibrating sensors
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A61B5/7282—Event detection, e.g. detecting unique waveforms indicative of a medical condition

Definitions

  • the present invention relates to methods and devices for use in estimating the concentration of a chemical constituent in a material, preferably without the need for invasive or direct measurement.
  • a great deal of control systems operate on the basis of knowing the levels of various substances within a material, often where the material is made of multiple different substances. This is typically measured by taking a sample of the material. This sample is then used in various chemical and physical tests to determine its composition. Such direct measurements may be termed "invasive" methods.
  • This process has the advantage of providing a high accuracy reading of the blood glucose level.
  • invasive means In particular, they are inconvenient and tend to require significant training to safely use. More problematically, even when used properly, taking blood samples may cause tissue damage, and can result in infection.
  • optical sensors may provide an estimate of the level of the substance (such as glucose) in the material (such as blood).
  • substance such as glucose
  • non-invasive methods are more convenient, they tend to be of low accuracy. This may be because, for example, a sensor measures not only the desired substance, but also other undesired substances in the material.
  • a method for determining a concentration of a chemical constituent in a material comprising:
  • a method for determining a concentration of a chemical constituent in a material comprising: obtaining a first measured spectrum from a first noninvasive sensor, the first measured spectrum relating to a first characteristic of the material; obtaining a second measured spectrum from a second noninvasive sensor, the second measured spectrum relating to a second characteristic of the material; calculating a first intermediate estimate of the concentration of a chemical constituent in the material based on the first measured spectrum; calculating a second intermediate estimate of the concentration of a chemical constituent in the material based on the second measured spectrum; determining a first and second residual spectrum, the residual spectra defined as the deviation between the measured spectrum and an estimated spectrum, the estimated spectrum being the component of the measured spectrum corresponding to the intermediate estimate of the concentration of the chemical constituent; and calculating an improved estimated concentration of the chemical constituent in the material based on the first intermediate estimate, the second intermediate estimate and the residual spectra r- ⁇ , r 2 .
  • Calculating the intermediate estimates of the level of the substance in the material can be performed in any suitable manner.
  • calculating a first intermediate estimate y x of the concentration of a chemical constituent in the material comprises: providing the first measured spectrum u)i as input to a predetermined first calibration model describing the relationship between the corresponding characteristic and concentration of the chemical constituent; and receiving an output from the first calibration model as the first intermediate estimate y x and wherein calculating a second intermediate estimate y 2 of the concentration of the chemical constituent in the material comprises: providing the second measured spectrum ⁇ 2 as input to a predetermined second calibration model describing the relationship between the corresponding characteristic and concentration of the chemical constituent; and receiving an output from the calibration model as the second intermediate estimate y 2 level of the substance in the material.
  • the loading matrix P may be estimated from a separate calibration data set using regression.
  • the calibration model may then be utilised, following a sensor reading from a material, to fit the estimated spectrum (corresponding to the estimated concentration of the chemical constituent) to a measured spectrum ⁇ and accordingly output an intermediate estimate of the concentration according to that sensor. In this way, an initial reading from the sensor can be used to provide an intermediate estimate of the concentration using a calibration model specific to that sensor in a straightforward manner.
  • the residual spectrum may then provide a means to further adjust the estimate using a function of the residual as a correction term.
  • the calibration model may be configured to model the change in the level of the substance in the material over time.
  • a model is easily adaptable and updatable, so that the calculation of the estimated level can be kept highly accurate.
  • the first correction term can be formulated as and the second correction term can be formulated as
  • S t is a first coefficient
  • ⁇ 2 is a second coefficient
  • r 1 is a first residual spectrum
  • r 2 is a second residual spectrum
  • rsd is the residual standard deviation of a residual spectrum.
  • y may be the estimate from a physical model, such as a metabolism model, which may be combined with the first and second intermediate estimates to further improve the accuracy of the final estimate.
  • one or more of the first scaling factor, the second scaling factor, the third scaling factor and the actual measurement of the level of the substance in the material are calculated by means of measurement of the level of the substance in the material.
  • one or more of the first correction term, the second correction term and the bias term are calculated by means of a recursive least squares regression or a Kalman filter.
  • the approach can be used for any substance and any material. However, it has been found particularly useful where the substance is glucose and the material is the blood of a user. This avoids the need for conventional invasive glucose measurement.
  • the first non-invasive sensor comprises an optical sensor configured to measure absorbance of light and the second sensor noninvasive comprises a bio-impedance sensor configured to measure the impedance of skin of the user. Such an arrangement provides measurements which lead to an accurate estimate.
  • the method further comprises: receiving a third reading representing a quantity of food ingested by the user; and receiving a fourth reading representing a quantity of exercise performed by the user; wherein calculating an estimated level of the substance in the material further comprises: inputting the third reading and the fourth reading into the model; wherein the output of the model is further based on the third reading and the fourth reading.
  • the model can also be adapted to model the metabolism of the user. This allows the estimate to have an improved predictive power.
  • a device configured to perform the method of the first aspect.
  • Figure 1 shows a method for non-invasively obtaining an estimated concentration of a chemical constituent in a material
  • Figure 2 shows a method for making use of a sensor fusion model in obtaining the estimated level of the substance in the material
  • Figure 3 shows the measured, estimated and ideal residual spectra corresponding to one sensor
  • Figure 4 shows a method for calibrating the sensor fusion model using actual measurements
  • Figure 5 shows a modified method for calibrating the sensor fusion model which based on a physical model.
  • FIG. 1 An overview of a method for non-invasively obtaining an estimated concentration of a chemical constituent in a material is shown in Figure 1.
  • the method is described in terms of two sensors (and consequently two readings). However, any number of sensors may be present. As a rule, a greater number of sensors will tend to provide a more accurate estimate.
  • Each of the sensors should be based on a different measurement principle or measuring position, and preferably measure different characteristics.
  • a first reading is obtained from a first sensor.
  • the first reading relates to a first characteristic of the material.
  • a second reading is obtained from a second sensor.
  • the second reading relates to a second characteristic of the material.
  • the second characteristic is different from the first characteristic.
  • the reading takes the form of a spectrum comprising a series of intensities across a plurality of values of the sensor parameter.
  • the sensor may be an optical sensor, where the characteristic being read is the absorbance of near infra-red light, the material and the input parameter is the wavelength of light output by the sensor.
  • the sensor parameter is the frequency used in the impedance measurement and the sensed characteristic is the impedance of the material.
  • the reading is a spectrum comprising the bioimpedance intensities across a range of measured frequencies.
  • the measured spectrum may consist of 10 intensity readings at 10 distinct values of the sensor parameter.
  • the form of the measured spectrum ⁇ received by the sensor is dependent on the concentration of the chemical constituents within the measured material. Some part of this measured spectrum is due to the particular chemical constituent of interest. The remaining part of the measured spectrum may be due to the response of other constituents of the measured material.
  • the sensors are chosen such that the desired substance in the material is particularly responsive to the measured characteristic. That is the sensors are chosen such that the measured spectrum changes significantly with the concentration of the substance of interest in the material and a significant portion of the measured spectrum is due to the chemical constituent to be measured.
  • the sensor may be configured to maximise any such effect. For example, if the desired substance is glucose, it is known that the absorbance of light in blood differs based on the level of glucose in the blood. Because of this, an optical sensor may be selected to measure the characteristic of absorbance of light in the blood. Alternatively, it is known that the impedance measured at a user's skin differs based on the glucose level of their blood. Because of this, a bio-impedance sensor may be selected to measure the characteristic of impedance.
  • the first and second sensors are selected such that their respective characteristics are not correlated to each other over the same materials and substances. That is, there should be no known simple equation (such as a linear equation) that would map the variation (with concentration of the substance) in the measured spectrum from the first sensor to the variation in the measured spectrum of the second sensor.
  • steps 110 and 120 are described sequentially, they may be performed simultaneously.
  • a first intermediate estimate of the level of the substance in the material is calculated. This is based on the first reading obtained at step 110.
  • a second intermediate estimate of the level of the substance in the material is calculated. This is based on the second reading obtained at step 120.
  • steps 130 and 140 are described sequentially, they may be performed simultaneously. Steps 130 and 140 may comprise using a formula, algorithm, look-up table, model or the like in order to calculate an intermediate estimate of the concentration based on the reading from the sensor. Alternatively, a machine learning approach (such as a neural network) may be utilised with a suitable training set to determine a relation between the measured spectrum and the estimated concentration of the chemical constituent within the material. The selection of the appropriate approach will be made easily using the common knowledge in the art of the relationships between the measured characteristic, the substance and the material.
  • steps 130 and 140 each use a model, calibrated to the spectral response of the corresponding sensor, which takes as an input the measured spectrum and outputs an estimated value of the concentration of the chemical constituent in the material - the "intermediate estimate".
  • This model may be referred to as a calibration model since it is calibrated using a calibration data set to relate the spectral response of the sensor to a concentration of the particular chemical constituent for which it is calibrated.
  • the calibration model is a partial least squares regression model, a principal component regression model or other similar modelling approaches based on latent variables.
  • Such approaches may proceed by fitting an estimated spectrum to the measured spectrum using one of the approaches listed above.
  • the form of the model may be developed using a separate calibration data set.
  • the calibration model may be used to obtain an estimate of the chemical constituent in the material (based on the fitted estimated spectrum) and a residual spectrum.
  • the residual spectrum is the difference between the measured and estimated spectrum, which is due to other substances being present in the material which have not been calibrated for - the calibration model being calibrated to model the spectral response of the chemical constituent of interest.
  • an estimated level of the substance in the material is calculated. This is based on the first intermediate estimated concentration and the second intermediate estimated concentration.
  • the calculation may preferably also use a residual spectrum where the residual spectrum is the component of the measured spectrum which corresponds to the deviation between the measured spectrum and estimated spectrum.
  • the final estimated concentration of the chemical constituent within the material is based on a sensor fusion model wherein the model uses the estimates from each sensor and the corresponding residual spectra to output an improved, final estimated concentration value.
  • a first step in obtaining the estimates of a chemical constituent is to identify a calibration model.
  • the identification of this model is done on a separate calibration data set.
  • this model can be a principal component regression model, a partial least squares regression model or other modeling approaches based on latent variables.
  • the calibration model is then used for obtaining the estimates of the chemical constituents and the residual spectra.
  • W TP + R
  • T the scores matrix
  • P the loading matrix.
  • the residuals are contained in R.
  • the loading matrix P may be estimated from a the calibration data set using regression.
  • the calibration model may then be utilised, following a sensor reading from a material, to fit an estimated spectrum (corresponding to the estimated concentration of the chemical constituent) to a measured spectrum ⁇ and accordingly output an intermediate estimate of the concentration according to that sensor.
  • the calibration model is applied to the measured spectrum ⁇ from the material using the form of the matrix P determined from the separate calibration data set.
  • the model is used to fit an estimated spectrum to the measured spectrum and output an intermediate estimate of the concentration associated with the calculated estimated spectrum.
  • the output of each sensor can be seen to be a measure of the corresponding characteristic of the overall material.
  • the estimated spectrum may be expressed as tP, where P is the loading matrix which has orthogonal rows and is estimated from the calibration data set using partial least squares regression or other methods based on latent variables.
  • the scores, t is estimated as raw responses weighted by the factor loadings P and is defined as the projection coordinates in the multidimensional space.
  • the measured spectrum ⁇ is shown for a sensor configured according to a parameter. For example this may be the measured absorbance intensities across each of 10 near infra-red wavelengths output by an optical sensor or the measured bioimpedance intensities across a plurality of frequencies. Ideally, all the variation in the spectra should be described by the calibration model. However, because the material itself may respond to the sensor and the material may contain uncalibrated substances (as shown by residual spectrum portion r), the overall measured spectrum ⁇ , combined with the calibration models may not give an accurate estimate of the chemical constituent. In this case, the residual portion of the spectrum is that which corresponds to the unmodelled constituents in the material and the material itself.
  • the intermediate estimates from each sensor are therefore provided to a sensor fusion model which allows for a more accurate final estimate of the concentration to be found based on the combined output of the sensors.
  • the individual intermediate estimates from each sensor may not be particularly accurate due to the other components present in the material. Taken alone, it is therefore difficult to estimate the concentration of the chemical constituent with sufficiency accuracy.
  • the principle here is that since each measured charateritic from the plurality of sensors has a differing relationship with the components present in the material, readings from multiple sensors may be used to more accurately identify the response due to the chemical constituent of interest.
  • a calibration model for the relationship between the spectral response and concentration for each sensor is used to calculate a value of the estimated concentration for each sensor. Therefore an estimated value of the concentration y 1 is calculated from the measured spectrum ⁇ from the first sensor and an estimated value of the concentration y 2 is calculated from the measured spectrum ⁇ 2 from the second sensor.
  • a sensor fusion model may then be used which relates an actual measured value of the chemical concentration y to the estimated values from each sensor.
  • step 150 comprises making use of one or more sensor fusion models, each of which takes one or more of the intermediate estimates as inputs.
  • These models may be in series, so that the output of a first model can be used as an input to a second model and so on, such that the final model provides the estimated level of the substance in the material as an output.
  • the models may be in parallel, so that the outputs of two or more models are combined, for example by taking an average (mean) value. In some cases, a combination of models in series and parallel may be used.
  • the residual spectra r - the difference between the measured and estimated spectra for each sensor - may also be input into the model to provide the final estimate of the chemical concentration in the material.
  • a correction term which is a function of the residual specttrum, ⁇ ( ), may be introduced into the calculation which corresponds to the portion of the intermediate estimate that is caused by other substances in the materials, as described below.
  • step 150 using a model is shown in Figure 2.
  • one or more of the intermediate estimates and the residual spectra are input into the sensor fusion model.
  • the model processes the inputs. These may be combined with an internal state of the model.
  • the internal state may be a pre-defined starting state for the model, if the model has not previously received an input. As a result of the processing, the model produces an output which is based on the intermediate estimates.
  • step 220 may comprise using a function which models the expected change of the level of the substance in the material over time. This is particularly applicable where the level of the substance changes at a reasonably consistent rate (such as linearly or logarithmically with respect to time). For example, a half-life of a radioactive substance in a substrate or the metabolism of a biological substance in a living being would benefit from such a model.
  • the output from the model is received. This output can be used as the estimated level of the substance in the material.
  • the model may be a general model which is configured for a wide range of situations and uses. However, in preferred embodiments, the model may be configured to be specific for the particular material or application, in order to give a more accurate result. For example, a general model may map the metabolism of a drug in a nominal average person, whereas a specific model may map the metabolism of a drug in the particular user.
  • a sensor fusion model relates an actual measured value of the chemical concentration y to the estimated values from each sensor.
  • ⁇ 3 ⁇ 4 and a 2 are scaling factors-these will typically be scalar values, but may in some cases be a function respectively; gi and g 2 are error terms which can relate to a measurement inaccuracy; and
  • /(ii) and /(r 2 ) are correction terms which may be a function of the residual spectra r.
  • An exemplary calibration method for the sensor fusion model is shown in Figure 4. This calibration step may be performed before or after each estimate calculation (as in step 150), or may only occasionally be performed.
  • one or more actual measurements y of the concentration of a chemical constituent in the material are obtained. For example, in the case of blood glucose, a series of invasive measurements of the user's blood glucose level are taken.
  • the scaling factors are calculated based on the one or more actual measurements. In some cases, this is done by taking sample readings from each of the sensors in combination with a conventional invasive measurement.
  • the scaling factors can be derived using a recursive least squares regression or a Kalman filter.
  • step 320 need not be performed before each calibration.
  • the scaling factors, ⁇ 3 ⁇ 4 and 2 may be precalculated for a number of subsequent uses. They may alternatively be predetermined, allowing step 320 to be omitted.
  • the correction terms, ⁇ ( ⁇ ) and f 2 (r 2 ), which are a function of the residual spectra, are calculated.
  • the correction terms, ⁇ ) and f 2 (r 2 ) can be derived from the residual spectra, using a recursive least squares regression or a Kalman filter as described in the above section.
  • the correction terms may be used as part of step 150, such that the estimated concentration is calculated based, in part, on the correction terms.
  • the correction terms are proportional to the residual standard deviation of the spectra ⁇ . These measurements will typically be derived from a sum of each of the components of the material. However, the contribution of each material to the overall value will differ, depending on the type of measurement used.
  • the residual standard deviation provides a good approximation of the error due to non-calibration for response of other substances in the material creating the residual spectrum.
  • the prediction error is proportional to the residual standard deviation (rsd) of the fitted spectrum containing an un-modelled chemical component.
  • a similar calibration technique may also be used to provide feedback to each model.
  • the output of the model may not be precisely the same as an actual measurement. It is possible to calibrate the model to improve its accuracy by feeding back the correction terms to the model.
  • ⁇ ⁇ is a scaling factor, which is preferably 1
  • ⁇ ⁇ is the estimated level of the substance in the material
  • ⁇ is a bias represents an inaccuracy in the output, ⁇ , of the model in comparison with the actual measurement, y.
  • ⁇ ⁇ may be an estimated concentration of the chemical constiuent in the material, as derived from a physical model such as a metabolism model
  • step number is identical to a step in Figure 4, the same description should be used. This is omitted here for brevity.
  • the scaling factor, ⁇ ⁇ is calculated. This may be calculated in the same manner as the scaling factors ⁇ ; .
  • the bias term ⁇ is calculated. This may be calculated in the same manner as the other correction terms ⁇ ( ).
  • one or more of the correction terms ⁇ ( ) or the bias ⁇ are fed back to the model.
  • the internal state of the model may be updated based on the calculated correction terms, in order to improve the accuracy of the model in further use.
  • the sensor fusion model may be used during measurement to output an improved estimate by using the determined parameters to adjust the intermediate estimates, taking into account the correction terms based on the residual spectra.
  • the estimated level being outside of a pre-determined acceptable range is considered an alarm condition. As a consequence of this, an alarm or the like may be raised.
  • the estimated level is used for automated control of a system, such that the operation of the system is updated periodically based on the estimated level.
  • users with diabetes may be provided with an artificial pancreas configured to administer insulin to the user.
  • the estimated level may be provided to a control system of the artificial pancreas in order to provide for automated (or at least partly automated) insulin delivery.
  • a device may be provided which is configured to perform the method described above.
  • the device comprises at least two sensors, a processor and a memory.
  • the sensors will be located so that they can monitor a characteristic of the material.
  • the output of the sensors will be passed to the processor which performs the method noted above.
  • the processor may also implement the model, using the memory to store the internal state of the model.
  • One particularly beneficial application of the present invention is in the field of blood monitoring.
  • Various health-related characteristics of a person can be determined based on the levels of certain substances in their blood.
  • blood glucose levels are particularly important to people who have diabetes. If a person with diabetes has overly high blood sugar, they may need to urgently have insulin administered to avoid dangerous health consequences.
  • alanine ascorbate
  • lactate lactate
  • triacetin urea
  • a user is provided with two or more non-invasive sensors, generally on or near the user's skin. Each sensor is configured to measure a different characteristic of the user.
  • the sensors are in communication with a monitoring device.
  • one sensor is an optical sensor which is configured to measure the absorbance of near-infrared and/or visible light.
  • This can be arranged near a light emitting unit which is arranged to emit near-infrared and/or visible light towards the user. This may be positioned on the user's wrist or finger, or another part on the user's body.
  • one sensor is a bio-impedance sensor, which is configured to measure the impedance at a portion of the user's skin.
  • Readings are periodically or continuously obtained from the two or more sensors by the monitoring device. Based on these readings, intermediate estimates of the blood glucose level in the user's blood are calculated in with the method of Figure 2 using well-known mappings between blood glucose levels and sensor readings.
  • model can be configured to model the changes in blood glucose over time, for example due to ordinary metabolism in the user.
  • the monitoring device may be re-calibrated based on one or more actual measurements of a user's blood glucose. This allows the calibration method of Figures 4 and 5 to be performed. Preferably, this would be done once a day. A longer period would still allow the approach to function, though there is an increasing risk of inaccurate estimates being calculated.
  • the monitoring device may be configured to receive information from one or more other sources. This received information may be provided as further inputs to the model. The model can then be updated to take account of the information.
  • the monitoring device may be configured to receive exercise information related to exercise performed by the user. This may be received from an accelerometer associated with the monitoring device. The exercise information can then be passed to the model.
  • the model may be configured to reflect that a high level of exercise may, for example, cause the blood glucose level to reduce more quickly with respect to time.
  • the monitoring device may additionally or alternatively be configured to receive food information related to food consumed by the user. This may involve a lookup service, where a textual description of food is entered by a user and is converted to a number of nutrient parameters such as carbohydrates, protein, fats and the like.
  • the model may be configured to reflect that certain nutrients (such as carbohydrates) tend to cause an increase in blood glucose level shortly after consumption.
  • the calculated estimated level may be used to control a medical device, for example an insulin delivery system. Based on the estimated level, the insulin delivery system may calculate and/or deliver an appropriate quantity of insulin to the user.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Public Health (AREA)
  • Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Surgery (AREA)
  • Molecular Biology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Epidemiology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Optics & Photonics (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Signal Processing (AREA)
  • Primary Health Care (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Emergency Medicine (AREA)
  • Physiology (AREA)
  • Dermatology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

L'invention porte sur un procédé pour déterminer un taux d'une substance dans un matériel. Le procédé comprend les étapes consistant : à obtenir un premier spectre mesuré à partir d'un premier capteur non invasif, le premier spectre mesuré se rapportant à une première caractéristique du matériel ; à obtenir un second spectre mesuré à partir d'un second capteur non invasif, le second spectre mesuré se rapportant à une seconde caractéristique du matériel ; à calculer une première estimation intermédiaire ŷ x de la concentration d'un constituant chimique dans le matériel sur la base du premier spectre mesuré ; à calculer une seconde estimation intermédiaire ŷ 2 de la concentration d'un constituant chimique dans le matériel sur la base du second spectre mesuré ; à déterminer un premier et un second spectre résiduel r1, r2, les spectres résiduels étant définis comme l'écart entre le spectre mesuré et un spectre estimé, le spectre estimé étant la composante du spectre mesuré correspondant à l'estimation intermédiaire de la concentration du constituant chimique ; et à calculer une concentration estimée améliorée du constituant chimique dans le matériel sur la base de la première estimation intermédiaire ŷ 1, la seconde estimation intermédiaire ŷ 2 et les spectres résiduels r1, r2.
EP17702061.7A 2016-01-25 2017-01-25 Estimation d'une concentration d'une substance au moyen de plusieurs capteurs non invasifs Withdrawn EP3408772A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP16152594 2016-01-25
PCT/EP2017/051573 WO2017129634A1 (fr) 2016-01-25 2017-01-25 Estimation d'une concentration d'une substance au moyen de plusieurs capteurs non invasifs

Publications (1)

Publication Number Publication Date
EP3408772A1 true EP3408772A1 (fr) 2018-12-05

Family

ID=55405106

Family Applications (1)

Application Number Title Priority Date Filing Date
EP17702061.7A Withdrawn EP3408772A1 (fr) 2016-01-25 2017-01-25 Estimation d'une concentration d'une substance au moyen de plusieurs capteurs non invasifs

Country Status (3)

Country Link
EP (1) EP3408772A1 (fr)
CN (1) CN107924715A (fr)
WO (1) WO2017129634A1 (fr)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3636141A1 (fr) 2018-10-10 2020-04-15 Prediktor Medical AS Capteur de glucose sanguin vestimentaire
KR20200119501A (ko) 2019-04-10 2020-10-20 삼성전자주식회사 생체정보 추정 장치 및 방법

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6441388B1 (en) * 1998-10-13 2002-08-27 Rio Grande Medical Technologies, Inc. Methods and apparatus for spectroscopic calibration model transfer
JP4636762B2 (ja) * 1999-08-31 2011-02-23 シーエムイー テレメトリックス インク． 分光装置を較正する方法
US6575905B2 (en) * 2000-09-22 2003-06-10 Knobbe, Lim & Buckingham Method and apparatus for real-time estimation of physiological parameters
JPWO2003079900A1 (ja) * 2002-03-25 2005-07-21 山越　憲一 非観血血液成分値測定装置及び方法
US6968221B2 (en) * 2003-03-14 2005-11-22 Futrex, Inc. Low-cost method and apparatus for non-invasively measuring blood glucose levels
US20050192492A1 (en) * 2004-02-27 2005-09-01 Ok-Kyung Cho Blood sugar level measuring apparatus
CN101696931B (zh) * 2009-10-15 2011-08-31 清华大学深圳研究生院 物质浓度无损检测方法及其装置
CN103239239B (zh) * 2013-04-23 2014-11-26 天津大学 一种定幅值的动态光谱数据提取方法

Also Published As

Publication number Publication date
CN107924715A (zh) 2018-04-17
WO2017129634A1 (fr) 2017-08-03

Similar Documents

Publication Publication Date Title
US11437141B2 (en) Performance reports associated with continuous sensor data from multiple analysis time periods
US11096611B2 (en) Method for providing a signal quality degree associated with an analyte value measured in a continuous monitoring system
US20190261905A1 (en) Advanced calibration for analyte sensors
US5068536A (en) Method for providing custom calibration for near infrared instruments for measurement of blood glucose
EP2051620B1 (fr) Procédé et dispositif pour surveiller un paramètre physiologique
CA2853236C (fr) Procede et/ou systeme de surveillance d'un analyte present dans plusieurs compartiments
Freckmann et al. Standardization process of continuous glucose monitoring: traceability and performance
US20090240440A1 (en) Non-Invasive Glucose Monitoring
US20140024907A1 (en) Method and system to indicate hyperglycemia or hypoglycemia for people with diabetes
JP5726286B2 (ja) 自己完結型の体外診断装置
RU2012117828A (ru) Способ тестирования аналита и устройство для контроля для пациентов, страдающих сахарным диабетом
WO1999056616A1 (fr) Mesure non invasive du glucose dans le sang
EP3236851B1 (fr) Procédé et système pour déterminer un niveau de glucose dans le sang d'un patient
EP4098194A1 (fr) Procédé de calcul de la sensibilité d'étalonnage d'un capteur destiné à être introduit dans un corps
Barcelo-Rico et al. Adaptive calibration algorithm for plasma glucose estimation in continuous glucose monitoring
Haller et al. Adverse impact of temperature and humidity on blood glucose monitoring reliability: a pilot study
Vettoretti et al. Modeling the error of factory-calibrated continuous glucose monitoring sensors: application to Dexcom G6 sensor data
EP3408772A1 (fr) Estimation d'une concentration d'une substance au moyen de plusieurs capteurs non invasifs
CA2952759C (fr) Procede et dispositif de determination d'un taux de glucose de fluide corporel d'un patient et produit programme d'ordinateur
US8734347B2 (en) Analytical method and investigation system
JP6521307B2 (ja) 近赤外分光法を用いた生体特性値測定方法及び測定装置
US20180172618A1 (en) Verifying Operation of a Meter
EP3636141A1 (fr) Capteur de glucose sanguin vestimentaire
Bailey Blood Glucose Monitoring: Glycated Hemoglobin, Fructosamine, Meters, and Sensors

Legal Events

Date Code Title Description
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: UNKNOWN

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE INTERNATIONAL PUBLICATION HAS BEEN MADE

PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: REQUEST FOR EXAMINATION WAS MADE

17P Request for examination filed

Effective date: 20180731

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

DAV Request for validation of the european patent (deleted)
DAX Request for extension of the european patent (deleted)
STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

17Q First examination report despatched

Effective date: 20201125

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN

18D Application deemed to be withdrawn

Effective date: 20210407