EP4008245A1 - Systeme permettant de determiner un indicateur representatif d'un parametre physiologique - Google Patents

Systeme permettant de determiner un indicateur representatif d'un parametre physiologique Download PDF

Info

Publication number
EP4008245A1
EP4008245A1 EP20211323.9A EP20211323A EP4008245A1 EP 4008245 A1 EP4008245 A1 EP 4008245A1 EP 20211323 A EP20211323 A EP 20211323A EP 4008245 A1 EP4008245 A1 EP 4008245A1
Authority
EP
European Patent Office
Prior art keywords
fitting
data values
fit
stage
determined
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
EP20211323.9A
Other languages
German (de)
English (en)
Inventor
Thomas POLESCHNER
Benjamin Stolze
Stephan Guido Maria REGH
Manuela Braun
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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 Koninklijke Philips NV filed Critical Koninklijke Philips NV
Priority to EP20211323.9A priority Critical patent/EP4008245A1/fr
Priority to CN202180081464.8A priority patent/CN116583219A/zh
Priority to JP2023533256A priority patent/JP2023551546A/ja
Priority to EP21819119.5A priority patent/EP4266992A1/fr
Priority to US18/039,580 priority patent/US20240090781A1/en
Priority to PCT/EP2021/083267 priority patent/WO2022117471A1/fr
Publication of EP4008245A1 publication Critical patent/EP4008245A1/fr
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/02028—Determining haemodynamic parameters not otherwise provided for, e.g. cardiac contractility or left ventricular ejection fraction
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/021—Measuring pressure in heart or blood vessels
    • A61B5/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
    • 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/7235—Details of waveform analysis
    • A61B5/7246—Details of waveform analysis using correlation, e.g. template matching or determination of similarity
    • 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/7235—Details of waveform analysis
    • A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/0205—Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/021—Measuring pressure in heart or blood vessels
    • A61B5/022—Measuring pressure in heart or blood vessels by applying pressure to close blood vessels, e.g. against the skin; Ophthalmodynamometers
    • A61B5/02225—Measuring pressure in heart or blood vessels by applying pressure to close blood vessels, e.g. against the skin; Ophthalmodynamometers using the oscillometric method

Definitions

  • the invention relates to an apparatus, a method and a computer program for determining an indicator that is representative for a physiological parameter.
  • EP 2 759 257 B1 discloses a method for determining an indicator that is representative for a patient's volume responsiveness.
  • a sequence of pulse signals of a patient over a time period corresponding to a plurality of subsequent respiratory cycles of the patient is measured, wherein the sequence of measured pulse signals is detected by a non-invasive pulse measurement method using a pressure cuff, wherein the pressure applied in the pressure cuff is continuously increased or decreased in the course of a measurement period.
  • the measured pulse signals are represented as pulse signals oscillating around an average value thereof, wherein an envelope signal curve for the pulse signals is determined.
  • a fit envelope signal function is determined based on the previously determined envelope signal curve, wherein the fit envelope signal function is determined based on a predetermined first functional prototype and represents an idealized curve progression of the envelope signal curve over the plurality of respiratory cycles without comprising pulse variations caused by ventilation or respiration induced heart-lung interaction.
  • respiratory pulse variation signals are determined, which correspond to the pulse variations caused by ventilation or respiration induced heart-lung interaction over the plurality of respiratory cycles, wherein the determined respiratory pulse variation signals correspond to a difference between the envelope signal curve and the fit envelope signal function.
  • an envelope respiration curve is determined for the previously determined respiratory pulse variation signals and a fit envelope respiration function is determined based on the previously determined envelope respiration curve, wherein the fit envelope respiration function is determined based on a predetermined second functional prototype and represents an idealized curve progression of the envelope respiration curve over the plurality of respiratory cycles.
  • the indicator that is representative for the patient's volume responsiveness is determined based on the fit envelope signal function and the fit envelope respiration function.
  • the fit envelope signal function is determined by fitting a Cauchy-Lorentz function to the envelope signal function by using a Levenberg-Marquardt algorithm.
  • the fit envelope respiratory function is determined by fitting the Cauchy-Lorentz function to the envelope respiratory function also by using the Levenberg-Marquardt algorithm.
  • the Levenberg-Marquardt algorithm can be computationally very intensive and hence time consuming such that the time needed for determining the indicator that is representative for the patient's volume responsiveness is generally very long.
  • an apparatus for determining an indicator that is representative for a physiological parameter comprises:
  • the physiological parameter is preferentially a physiological parameter of a patient.
  • the physiological parameter is indicative of, or corresponds to, a patient's volume or fluid responsiveness.
  • the determination of the indicator that is representative for the physiological parameter is preferentially a determination of a fluid responsiveness parameter (FRP), especially of a pulse pressure variation (PPV), a systolic pressure variation (SPV) and/or a stroke volume variation (SVV).
  • FRP fluid responsiveness parameter
  • the determined indicator may be the respective FRP, i.e. a PPV, an SPV and/or an SVV.
  • the physiological parameter can also be indicative of, or correspond to, for instance, a blood pressure, such as an arterial blood pressure.
  • the determination of the indicator that is representative for the physiological parameter can also be a determination of one or more parameters related to the blood pressure.
  • the determined indicator can be the respective parameter related to the blood pressure.
  • the provision of the sequence of measured pulse signals can correspond to measuring the sequence of pulse signals of the patient, particularly to measuring signals corresponding to blood pulsations by continuously or semi-continuously recording a blood pressure or a pressure indicative of the blood pressure.
  • the sequence of measured pulse signals has been detected by a non-invasive pulse measurement method using a pressure cuff, wherein the pressure applied in the pressure cuff is continuously increased or decreased in the course of a measurement period.
  • the measured signals might then correspond to continuous or semi-continuous recordings of a pressure in the pressure cuff.
  • the provision of the sequence of measured pulse signals may also comprise a providing of a sequence of pulse signals which have been measured previously by, for instance, a physician.
  • the measured pulse signals may hereafter have been stored and this stored sequence might be provided.
  • Non-invasive pulse measurements methods include, for instance, oscillometric methods.
  • oscillometric pulse measurement methods usually a pressure cuff is inflated such that a pneumatic coupling between blood pulsations and a pressure sensor recording a pressure in the inflated volume is generated.
  • the sequence of pulse signals is measured using a pressure cuff that is adapted to provide a hydraulic coupling between blood pulsations and a pressure sensor of the pressure cuff that is adapted to record the pressure indicative of the blood pressure.
  • a pressure cuff of this kind is disclosed, for instance, in WO 2014/121945 A1 .
  • a pressure cuff which will subsequently be referred to also as a high-fidelity pressure cuff, allows for an improved quality of the recorded pressure signal as compared to oscillometric measurements relying on a pneumatic coupling.
  • the pressure cuff may be adapted to bring the pressure sensor into direct contact with the patient's tissue such that it is hydraulically coupled to blood pulsations via the tissue.
  • the pressure cuff may, for instance, be a high-fidelity upper arm pressure cuff.
  • the determination procedure that the indicator determination unit is adapted to carry out includes a determination of data values based on the provided sequence of measured pulse signals.
  • the pulse signals providing unit may be configured to provide the sequence of measured pulse signals as a continuous or a quasi-continuous digital signal, wherein the indicator determination unit may be adapted to determine the data values by discretizing the provided sequence of measured pulse signals and/or determining a start and an end point of the provided sequence of measure pulse signals.
  • the data values may then still be understood as a direct representation, or sampling, of the sequence of measured pulse signals.
  • Determining the data values based on the provided sequence of measured pulse signals may, however, additionally or alternatively also comprise other types of preprocessing of the measured pulse signals, such as, for instance, the determination of an envelope signal curve, of respiratory pulse variation signals and/or of a signal or curve derived from the respiratory pulse variation signals, such as an absolute respiratory pulse variation curve, like an envelope respiration curve or a another curve indicative of absolute values of the respiratory pulse variation signals, for the measured pulse signals.
  • the data values may then correspond to the envelope signal curve or the signal or curve derived from the respiratory pulse variation signals, respectively, i.e., for instance, to a discretized version of it.
  • the determination procedure further includes a fitting of a predetermined functional prototype to the determined data values.
  • the functional prototype can depend on one or more fit parameters to be modified during fitting, wherein an initial value for each of the one or more fit parameters can be determined.
  • the initial value for each or only some of the one or more fit parameters can be determined based on characteristics of the data values which are related to characteristics of the functional prototype and/or based on a value of the respective fit parameter known from a previous fitting.
  • the initial value might be predefined or determined based on characteristics of the measurement method, such as rate of change in pressure applied in a pressure cuff, for instance.
  • the functional prototype can be predetermined depending on an expected form of the sequence of measured pulse signals, wherein the form might be expected, for instance, based on the measurement method. Moreover, the functional prototype can also be predetermined depending on a chosen manner of determining the data values based on the provided sequence of measured pulse signals, wherein the chosen manner may correspond, for instance, to a defined preprocessing.
  • the pulse signals providing unit is configured to provide several sequences of measured pulse signals of a patient over time, wherein each of the sequences of measured pulse signals has been detected in a different time period and wherein the indicator determination unit is configured to determine a set of data values for each of the provided sequences of measured pulse signals.
  • the different time periods may, for instance, correspond to subsequent cycles of inflation and/or deflation of the pressure cuff.
  • any known fitting algorithm may be used.
  • a non-linear least squares algorithm may be used, particularly the Levenberg-Marquardt algorithm. It is understood, however, that the above indicated effect of reduced computational efforts can also be achieved for different fitting algorithms.
  • the indicator determination unit is configured to reduce the number of data values such that the reduced number of data values fulfills the Shannon-Nyquist criterion with reference to a bandwidth of the provided functional prototype and/or with reference to the data values themselves.
  • the number of data values may be reduced such that the reduced number of data values fulfills the Shannon-Nyquist criterion with reference to the provided sequence of pulse signals based on which the data values have been determined. This may particularly be the case if the data values are determined such that they are representative of the provided sequence of pulse signals.
  • the bandwidth of the provided functional prototype, particularly the frequencies contained therein, may, for instance, be estimated based on a time elapsed during the increase or decrease in pressure in the pressure cuff from a first pressure value to a second pressure value, wherein the first pressure value and the second pressure value may refer, for instance, to a diastolic and a systolic arterial pressure value, respectively. If the data values are determined such that they are representative of the respiratory cycles of the patient, the number of data values may be reduced such that the reduced number of data values fulfills the Shannon-Nyquist criterion with reference to the respiratory cycles of the patient.
  • a reduction of data values to be fitted can be done by selecting one data value per second.
  • the indicator determination unit is configured to reduce the number of data values by determining one data value at a predetermined frequency and/or by downsampling. In a further embodiment the indicator determination unit is configured to lowpass filter the data values before reducing the number of data values.
  • the functional prototype is a bell-shaped function. This is particularly preferred if the sequence of measured pulse signals has been detected using a pressure cuff, wherein the pressure applied in the pressure cuff is continuously increased or decreased in the course of a measurement period.
  • the bell-shaped function is preferentially a Cauchy-Lorentz function.
  • the bell-shaped function can have a first fit parameter being indicative of the height of the maximum of the bell-shaped function and a second fit parameter being indicative of the position of the maximum of the bell-shaped function, wherein the indicator determination unit can be configured to determine a height of a maximum of the data values and a position of the maximum of the data values as initial values for the first and second fit parameters.
  • Providing initial values for the fit parameters, which are close to the final values for the fit parameters, can reduce calculation time significantly.
  • This prediction can be done by analyzing the data values for characteristics that describe the functional prototype. Predetermination and preselection of characteristic patterns from the data values can make the data values fit better to the functional prototype and thereby reduce the total number of iterations during fitting. For example, by fitting, as explained above, to a bell-shaped functional prototype, in which at least two of the fit parameters are the height and position of the maximum, the height and position of the maximum of the data values can be used as initial values for the fit parameters.
  • a third fit parameter on which the bell-shaped functional prototype can depend might relate to the full width at half maximum.
  • An initial value for this third fit parameter can be determined without reference to the data values, namely a) based on the time elapsed during an increase in pressure in the pressure cuff from a first pressure value to a second pressure value or a decrease in pressure in the pressure cuff from the second pressure value to the first pressure value, wherein the first pressure value and the second pressure value may refer, for instance, to a diastolic and a systolic arterial pressure value, respectively, or b) directly as the width of a pressure range including the first and the second pressure value, wherein the pressure range may be, for instance, a range from 0.7 times the diastolic arterial pressure value to 1.3 times the systolic arterial pressure value or a range from 10 mmHg below the diastolic arterial pressure value to 20 mmHg above the systolic arterial pressure value.
  • the diastolic and systolic arterial pressure values that may be used for determining the third fit parameter may refer, for instance,
  • the indicator determination unit is configured such that the fitting is carried out in three stages, wherein in a first stage a first percentage of the data values, in a following second stage a higher, second percentage of the data values and in a following third stage 100 percent of the data values are used.
  • the first percentage can be 1/10 percent and the second percentage 1 percent.
  • the first percentage and second percentage can also be, for instance, percentages corresponding to 1/250 and 1 / 250 , i.e. 0.4 percent and approximately 6.32 percent, respectively.
  • the indicator determination unit is configured such that the fitting is carried out in several stages, wherein, after the fitting has been carried out in a stage, it is determined whether the fitted fit parameter is in conformance with an abort criterion, wherein, if this is the case, the fitting is terminated.
  • the indicator determination unit can be configured to use a deviation between a) a fitted value of the fit parameter, which resulted from a respective stage, and b) an expected value of the fit parameter as the abort criterion.
  • the fit algorithm is not yielding valid data, for instance, because the underlying data values are not forming a shape related to the functional prototype, analyzing the fit results of each stage before the final stage on validity can lead to the conclusion that not all stages have to be executed, thereby avoiding unnecessary calculation effort and therefore reducing runtime. This can be done, as explained above, when there is a certain range of expected values of the final fit parameters and the stages before the final stage yield fit parameters too far off these ranges.
  • the abort criterion indicates that the fit parameters resulting from a stage are already very well, for instance, by determining whether a deviation between fit parameters resulting from a current stage from fit parameters resulting from a previous stage is smaller than a predefined threshold, wherein in this case the fitting can be stopped, before the final stage has been carried out. Also this allows for a reduced calculation time.
  • fitting When fitting is carried out in several stages, i.e. in the case of multiple-staged fitting, and if the original number of data values is used for the fitting in the last stage, the final results, particularly the accuracy of the final determined fit parameters, can be the same as compared to a single-staged fitting to the original number of data values, even though the overall calculation time can be reduced. To allow for a further reduction in overall calculation time, fitting may also be carried out in several stages such that even for the fitting in the final stage less than the original number of data values is used.
  • the sequence of measured pulse signals may be detected over a configurable time period.
  • the time period is configured such that it covers a defined number of subsequent respiratory or ventilation cycles of the patient, wherein the predetermined number may be any number from one to ten.
  • the time period may be configured to have a length of between, for instance, 10 seconds and four minutes, preferably between 30 seconds and two minutes, more preferably about one minute.
  • the pulse rate of the patient depends on various factors, such as age, stress, et cetera.
  • the heart of an adult beats usually between 50 to 90 times per minute. Thus, a comparatively larger number of pulse variations caused by heartbeats can be detected in the detection time period.
  • the detection period should preferably not exceed three minutes so as to avoid adverse effects due to disturbed blood flow caused by the pressure of the pressure cuff.
  • the patient's respiration has a detectable influence on pulse variation. While the measured variation of the blood pressure is primarily derived from the function of the heart, i.e. from its cyclic contractions and relaxations, there is also another influential factor to be considered: Thus, two functions are superimposed: Variations of higher frequency caused by the function of the heart are superimposed by variations of lower frequency caused by the respiration or ventilation of the patient. Notably, such low frequency variations caused by the patient's respiration are not only detected with mechanically ventilated patients, but also with non-ventilated patients breathing spontaneously.
  • the apparatus can be equally applied to mechanical ventilated and spontaneous breathing non-ventilated patients, provided that the respiration - or ventilation - induced maneuver produces significant heart-lung-interaction.
  • the measured pulse signals may be represented as a function over time or, alternatively, as a function over a clamping pressure exerted by the pressure cuff, wherein the pressure applied in the pressure cuff is continuously increased or decreased in the course of the measurement period.
  • the envelope signal curve to which the first data values may correspond, is determined based on the sequence of measured pulse signals.
  • an envelope curve or function of a rapidly varying signal is considered to be a smooth curve outlining the extremes in amplitude of the rapidly varying signal.
  • the envelope curve or function may be determined by simply connecting the maxima or the minima of the rapidly varying signal.
  • the envelope curve of the blood pressure signals i.e. the envelope signal curve, is determined by continuously determining a distance dimension of the measured pulse signals from an average thereof and, thereafter, preferably by applying a low pass filter to the distance dimension.
  • the low pass filter has a cutoff frequency below the pulse rate of the patient.
  • the portion below the average value of the oscillating curve of measured pulse signals is preferably folded up to the upper portion.
  • the resulting curve is preferably flattened by using a low pass filter having a cutoff frequency below the pulse rate of the patient.
  • the resulting curve can be flattened in such a way that the area below the flattened curve remains unchanged as compared to the area under the non-flattened curve.
  • the flattened curve may additionally be multiplied by a predetermined value. If the flattened curve is multiplied, for example, by the square root of 2, the finally obtained envelope signal curve substantially is at the level of the upper extremes in amplitude of the measured pulse signals.
  • the distance dimension could also be defined as the squared signal or the extreme values within a certain region, for instance, of about one pulse width or any other metric function. More generally, the distance dimension could also be defined, for instance, as the absolute value of the difference between the signal and its average value to the power of n, wherein n may be 1 or any other number. In case of using the absolute value of the difference as a distance dimension, the result correlates well with the SVV as well as the PPV. Also with the absolute value of the difference to the power of 2 as a distance dimension, the result correlates well to SVV as well as PPV.
  • the absolute value of the difference to the power of 2 is used as the distance dimension if the SVV is to be determined, and the absolute value of the difference is used as the distance dimension if the PPV is to be determined.
  • the correlation may be further increased.
  • the maximum metric i.e. with infinite n as a distance dimension, the result correlates very well with PPV.
  • the maximum metric could be realized by searching the maximum value minus the minimum value of the difference within a moving window equal to the duration of a heartbeat.
  • the average value, which may be used as the basis of the distance dimension calculation, of the measured pulse signals can be determined as moving average over a period of one single pulse cycle of the patient.
  • the first element of the moving average is obtained by taking the average of the initial fixed subset of the number series.
  • the subset is modified by "shifting forward", that is excluding the first number of the series and including the next number following the original subset in the series. This creates a new subset of numbers, which is averaged. This process is repeated over the entire data series.
  • a window function is advantageous.
  • pulse signals measured by a continuous invasive blood pressure measurement method usually do not exhibit any bell-shaped form. If pulse signals are measured using a non-invasive blood pressure measurement method employing a high-fidelity pressure cuff, application of a window function may not be required.
  • measured pulse signals usually already exhibit a bell-shaped form, and, therefore, are - as such - well fittable with a functional prototype also exhibiting a bell-shaped form.
  • the window function is preferably a non-negative smooth bell-shaped curve, for example a Cauchy-Lorentz function.
  • the determination of the fit envelope signal function represents an idealized curve progression of the envelope signal curve with the object to exclude any pulse variation caused by ventilation or respiration induced heart-lung interaction.
  • the difference between the envelope signal curve and the fit envelope signal function reflects a modulation which is due to the respiration or ventilation of the patient.
  • respiratory pulse variation signals corresponding to the pulse variations caused by the respiration of the patient are determined.
  • the respiratory pulse variation signals are determined in such a way that the respiratory pulse variation signals oscillate around an average value thereof.
  • the area defined by the lower part of the curve of the respiratory pulse variation signals i.e. the area below the average value, substantially corresponds to the area defined by the upper part of the curve of the respiratory pulse variation signals, i.e. the area above the average value.
  • the absolute respiratory pulse variation curve may correspond to an envelope respiration curve, in which case the fit respiration function might also be regarded as a fit envelope respiration function.
  • the fit respiration function may then represent an idealized curve progression of the envelope respiration curve over the plurality of respiratory cycles.
  • the determination of the envelope respiration curve is based on the previously determined respiratory pulse variation signals.
  • the envelope respiration curve is calculated by continuously determining a distance dimension of the respiratory pulse variation signals from an average thereof and by thereafter preferably applying a low pass filter to the distance dimension.
  • the low pass filter has a cutoff frequency below the respiration frequency of the patient.
  • the portion below the average value of the oscillating curve of respiratory pulse variation signals is preferably folded up to the upper portion.
  • the resulting curve is preferably flattened by using a low pass filter having a cutoff frequency below the respiration frequency of the patient.
  • the curve can be flattened in such a way that the area below the flattened curve remains unchanged as compared to the area under the non-flattened curve.
  • the flattened curve may additionally be multiplied by a predetermined value. If the flattened curve is multiplied, for example, with the value of the square root of 2, the finally obtained envelope respiration curve substantially is at the level of the upper extremes in amplitude of the respiratory pulse variation signals.
  • the average value, which is used as the basis of the distance dimension calculation, of the respiratory pulse variation signals is determined as moving average over a period of one single respiration cycle of the patient.
  • the absolute respiratory pulse variation curve does not correspond to an envelope respiration curve, but to a different curve indicative of absolute values of the respiratory pulse variation signals. That is to say, it may be preferred that the fit respiration function is determined without first determining an envelope respiration curve for the previously determined respiratory pulse variation signals and then determining the fit respiration function based on the envelope respiration curve. Instead, it may be preferred to determine the fit respiration function by fitting the second functional prototype to data values, i.e. second data values, corresponding to a different kind of absolute respiratory pulse variation signals, i.e. to data values derived differently from the respiratory pulse variation signals. Nevertheless, such different kind of preprocessing might involve one or more of the steps described above with respective to the determination of the envelope respiration curve from the respiratory pulse variation signals.
  • the absolute respiratory pulse variation curve which may correspond to the second data values, is determined by considering absolute values of respective differences between the respiratory pulse variation signals and their average. These absolute values may then correspond to the second data values.
  • the fit respiration function may be determined by considering absolute values of respective differences between the respiratory pulse variation signals and their average, thereby folding the portion of the respiratory pulse variation signals below an average value thereof up to the upper portion, and by fitting the second functional prototype to the so obtained values.
  • the envelope respiration curve is preferably calculated analogously to the envelope signal curve, but is associated with the respiratory pulse variation signals instead of the measured pulse signals. More generally, the fit respiration function is preferably calculated based on the data values, i.e. the second data values, to which it is to be fitted, i.e.
  • the envelope respiration curve or the data values derived differently from the respiratory pulse variations signals for instance, analogously to the fit envelope signal function, with the exception that the data values to which it is to be fitted are different, namely are associated with the respiratory pulse variation signals instead of the measured pulse signals. That is, for example, if the Cauchy-Lorentz function is used as functional prototype to determine the fit envelope signal function, also the Cauchy-Lorentz function is preferably used as functional prototype to determine the fit respiration function.
  • how, for instance, before the fitting, the number of data values is reduced for the entire fitting may be determined based on characteristics of the data values, i.e. for instance, whether they correspond to an envelope curve or a curve derived from an oscillatory signal by considering its absolute values.
  • the indicator determination unit is configured to determine the absolute respiratory pulse variation curve by considering absolute values of respective differences between the respiratory pulse variation signals and their average, and to reduce the number of second data values by determining one second data value for each local maximum of the absolute respiratory pulse variation curve.
  • the indicator that is representative for the physiological parameter can be determined based on the fit envelope signal function and the fit respiration function.
  • the indicator that is representative for the physiological parameter is determined based on at least one parameter of the fit envelope signal function and at least one parameter of the fit respiration function.
  • an indicator that is representative for a patient's volume responsiveness can be determined based on a ratio between a maximum of the fit envelope signal function and a maximum of the fit respiration function. Such a ratio, or a function of this ratio, particularly its inverse, represents an appropriate indicator that is representative for the patient's volume responsiveness.
  • further data values are determined by determining a further absolute respiratory pulse variation curve for the respiratory pulse variation signals, wherein, in a further fitting, a further functional prototype, which depends on a further fit parameter to be modified during the fitting, is fitted to the further data values such that a further fit respiration function results, wherein the indicator that is representative for the physiological parameter is determined based further on the further fit respiration function.
  • the further fitting may be performed such that a), before the further fitting, the number of further data values is reduced for the entire further fitting, and/or b) an initial value is determined for the further fit parameter based on characteristics of the further data values which are related to characteristics of the further functional prototype and/or based on a value of the further fit parameter known from a previous further fitting, and/or c) the further fitting is carried out in several stages, wherein the number of further data values used for the further fitting is increased from stage to stage and wherein in a current stage a value of the further fit parameter determined in a previous stage is used as initial value for the further fit parameter in the current stage.
  • the further functional prototype is preferably the same as the second functional prototype.
  • the further absolute respiratory pulse variation curve is preferably different from the second absolute respiratory pulse variation curve. It is also possible that several such further fittings are performed, i.e. two, three or more further fittings, wherein the indicator that is representative for the physiological parameter can then be determined based further on the several further fit respiration functions.
  • a method for determining an indicator that is representative for a physiological parameter comprises:
  • a computer program for determining an indicator that is representative for a physiological parameter comprises program code means for causing an apparatus as defined by any of claims 1 to 11 to carry out the steps of the method as defined in claim 12.
  • Fig. 1 shows schematically and exemplarily an apparatus 100 for determining an indicator that is representative for a physiological parameter.
  • the apparatus 100 comprises a pulse signals providing unit 101 configured to provide a sequence of measured pulse signals of a patient over a time period corresponding to a plurality of subsequent respiratory cycles of the patient. The sequence of measured pulse signals has been detected by a non-invasive pulse measurement method.
  • the apparatus 100 further comprises an indicator determination unit 102 configured to carry out a determination procedure adapted to determine an indicator that is representative for a physiological parameter based on the provided sequence of pulse signals.
  • the determination procedure includes a determination of data values s0, S0, S0', S0" based on the provided sequence of pulse signals and the fitting of a provided functional prototype, which depends on a fit parameter to be modified during the fitting, to the determined data values s0, S0, S0', S0".
  • the indicator determination unit 102 is configured such that a), before the fitting, the number of data values s0, S0, S0', S0" is reduced for the entire fitting, and/or b) an initial value is determined for the fit parameter based on characteristics of the data values s0, S0, S0', S0" which are related to characteristics of the functional prototype and/or based on a value of the fit parameter known from a previous fitting, and/or c) the fitting is carried out in several stages, wherein the number of data values s0, S0, S0', S0" used for the fitting is increased from stage to stage and wherein in a current stage a value of a fit parameter determined in a previous stage is used as initial value for the fit parameter in the current stage.
  • the apparatus 100 may comprise a connection between the pulse signals providing unit 101 and the indicator determination unit 102 via which the pulse signals providing unit 101 can provide the sequence of measured pulse signals to the indicator determination unit 102.
  • Fig. 2 illustrates schematically and exemplarily how an indicator that is representative for a patient's volume responsiveness can be determined based on a provided sequence of measured pulse signals that has been detected using a pressure cuff, wherein the pressure applied in the pressure cuff is continuously increased in the course of the measurement period.
  • the horizontal axis of the graph shown in Fig. 2 indicates the pressure applied in the pressure cuff, and could therefore also be viewed as a time axis.
  • the vertical axis indicates a pressure measured by the pressure sensor of the pressure cuff.
  • the measured pulse signals are represented as pulse signals p0 oscillating around an average value thereof. Based on the pulse signals p0, data values s0 are determined by determining an envelope signal curve for the pulse signals p0.
  • the data values s0 may be regarded as being, or being values of, the envelope signal curve.
  • an envelope respiration curve is determined by determining an envelope curve for the previously determined respiratory pulse variation signals r0.
  • a second functional prototype g x g amp 1 + x ⁇ g max g bw 2 , which is in this case of the same bell-shaped Cauchy-Lorentz type as the first functional prototype, is provided.
  • x is again the pressure or time indicated on the horizontal axis
  • the parameter g amp is decisive for the amplitude of the bell-shaped curve of the functional prototype
  • the parameter g max is decisive for the location of the maximum on the time-axis or pressure-axis
  • the parameter g bw is decisive for the width at half maximum.
  • the fit envelope respiration function g0 might be understood more generally as a fit respiration function g0.
  • the magnitude of the pulse variations caused by the patient's respiration is represented by the parameter g 0 amp .
  • the parameter g 0 amp is "normalized" in equation (3), namely by division through f 0 amp , which indicates a magnitude of the pulse variations caused by heart beats. Therefore, the volume responsiveness indicator obtained according to equation (3) may also be viewed as being indicative for the conventional indicator obtained according to equation (4). Differences might arise from the circumstance that the procedure leading to equation (3) does not rely on single maximum/minimum values of measured pulse pressure variations corresponding to one single heartbeat, but instead takes all pulse signals measured in the measurement time into account.
  • the determination procedure described above with reference to Fig. 2 is an example of a determination procedure used for determining an indicator that is representative for a physiological parameter - in this case a patient's volume responsiveness - which involves a fitting algorithm and which can therefore be computationally very intensive and hence time consuming.
  • this problem can be particularly present, since two fittings are involved.
  • Figs. 3 to 10 illustrate how this problem can be alleviated using the apparatus 100.
  • Fig. 3 shows, like Figs. 4 to 11 , schematically and exemplarily a graph whose axes correspond to the axes of the graph shown in Fig. 2 .
  • the graph of Fig. 3 includes data values SO that have been determined based on a provided sequence of pulse signals, and a fit function F0 corresponding to a fitting of the provided functional prototype to the data values S0.
  • the data values SO are not the values, i.e., for instance, sampled values, of the measured pulse signals, but have been determined based on a function derived from the sequence of measured pulse signals, wherein the function might be viewed as corresponding to an envelope signal curve like the envelope signal curve s0 in Fig. 2 .
  • the form of the fit function F0 in Fig. 3 corresponds to that of a provided functional prototype, which is, in this case, again a Cauchy-Lorentz function.
  • the fit function F0 has been fitted to the data values S0, wherein the values of the fit parameters have been determined, i.e. modified starting from respective initial values, during the fitting.
  • all of the data values SO have been taken into account.
  • Fig. 4 corresponds to Fig. 3 , wherein, however, instead of the full set of data values SO and the corresponding fit function F0, a reduced number of data values SR and a corresponding fit function FR fitted to the reduced number of data values SR is plotted.
  • the data values SR were obtained by selecting every thousandth of the data values S0.
  • the number of data values SO used for the fitting has thereby been reduced to 0.1 percent.
  • the fitting was carried out using the same functional prototype as used to arrive at the graph of Fig. 3 , namely a Cauchy-Lorentz function.
  • Fig. 5 which relates to Figs. 3 and 4 , no data values are plotted, but only the fit functions F0 and FR resulting, respectively, from the fitting to the entire set of data values SO and the reduced set of data values SR. It can be seen from Fig. 5 that even though the number of data values used to arrive at the reduced fit function FR was reduced by 99.9 percent compared to the entire number of data values S0, the reduced fit function FR only deviates by a relatively small amount from the fit function F0. The time required for fitting the functional prototype to the reduced number of data values, however, was substantially reduced compared to the time required for fitting the functional prototype to the entire number of data values S0.
  • the indicator that is representative for the physiological parameter is determined from the final fit parameters resulting from the fitting, such as, like in above equation (3), the maximum of the fit function
  • the indicator can be determined substantially faster based on the fitting to the reduced number of data values resulting in the fit function FR.
  • the fit function FR deviates only by a relatively small amount from the fit function F0, the faster determination of the indicator does not necessarily negatively affect the accuracy of the determined indicator.
  • the values of the one or more fit parameters on which the provided functional prototype depends and whose values are to be determined by fitting can be estimated before the fitting, wherein the estimated values may be used as initial values of the fit parameters for the fitting. This is schematically and exemplarily illustrated by Figs. 6 and 7 .
  • Figs. 6 and 7 relate to an embodiment in which the functional prototype is, again, a Cauchy-Lorentz function.
  • the Cauchy-Lorentz function has a first fit parameter being indicative of the height of its maximum and a second fit parameter being indicative of the position of its maximum, as already explained in relation to above equations (1) and (2).
  • Fig. 6 illustrates an embodiment in which the indicator determination unit 102 is configured such that an initial value is determined for the fit parameter based on characteristics of the data values SO which are related to characteristics of the functional prototype.
  • the indicator determination unit 102 is configured to determine a height of a maximum of the data values SO and a position of the maximum of the data values S0, and to determine these values as initial values for the first and the second fit parameter, respectively.
  • the initially assumed fit function i.e. the fit function assumed at the beginning of the fitting, will in this way be well adapted to the data values SO such that the amount of fitting to be done can be decreased.
  • the amount of fitting needed can also be decreased by only considering those data values for the fitting that capture the characteristics of the entire set of data values that are relevant in respect of the functional prototype.
  • the initially assumed fit function be adapted to the data values to which it is to be fitted, but also the subset of data values to which the fit function is to be fitted can be adapted to the type of the fit function. This is schematically and exemplarily illustrated in Fig. 7 .
  • Fig. 7 shows a graph in which data values SO' are plotted that have been determined based on the sequence of measured pulse signals using a different kind of preprocessing as compared to the data values S0.
  • the data values S0' which may be referred to as second data values, might be viewed as being derived from signals corresponding to the respiratory pulse variation signals r0 shown in Fig. 2 , and hence as being derived from data values corresponding to the data values S0, which may be referred to as first data values.
  • the data values SO' have in this case been derived by considering absolute values of respective differences between signals corresponding to the respiratory pulse variation signals r0 as shown in Fig. 2 and an average value thereof, i.e.
  • the data values SO' correspond to an absolute respiratory pulse variation curve that is not an envelope respiration curve.
  • the data values SO' comprise a series of peaks which vary in amplitude, i.e. height. Since, in this embodiment, the characteristic of interest of the data values SO' is still the progression in amplitude of the peaks, i.e. as opposed to the shape of any single one of the peaks, the number of data values considered for fitting can be reduced to those data values which capture this characteristic of interest, i.e. which capture the progression in peak amplitudes, i.e. heights.
  • the data values considered for fitting have been reduced to the data values corresponding to the local maxima.
  • the resulting reduced set of data values is annotated by SR' in Fig. 7 .
  • SR' fitting to the entire set of data values SO' results in a fit function G0'
  • a fitting to the reduced set of data values SR' results in a fit function GR'.
  • the deviation of the fit function GR' from the fit function G0' is noticeable, due to the rather large variations in the original data values S0'. Nevertheless, GR' allows for a good prediction of the fit parameters corresponding to G0'.
  • the position of the maximum of GR' is close to the position of the maximum of G0', and also the full width at half maximum of GR' is a good estimate of the full width at half maximum of G0'. Still, it may be preferred to apply a predefined scaling factor on the data values SO' before fitting or on the fit function GR'. In this way, it can be taken into account that by considering only the local maxima SR' of SO' for fitting, the height of the fit maximum will always be larger than when fitting to all data values S0'.
  • a further additional or alternative way to reduce the time needed for fitting and therefore also for determining the indicator that is representative for the physiological parameter based on the fit function resulting from the fitting, is schematically and exemplarily illustrated in Figs. 8 to 11 .
  • Figs. 8 to 11 relate to an embodiment according to which the fitting is carried out in several stages, wherein the number of data values used for the fitting is increased from stage to stage and wherein in a current stage a value of a fit parameter determined in a previous stage is used as initial value for the fit parameter in the current stage.
  • the indicator determination unit 102 is configured such that the fitting is carried out in three stages, wherein in a first stage a first percentage of the data values, in a following second stage a higher, second percentage of the data values and in a following third stage 100 percent of the data values are used.
  • the data values SR1" corresponding to the first percentage of the entire set of data values S0" are plotted in the graph of Fig. 8 together with the fit function FR1" fitted to the data values SR1".
  • the original data values S0" correspond to the original data values SO from the embodiment described with respect to Figs. 3 to 5 .
  • the reduced set of data values SR1" corresponds to the reduced set of data values from the embodiment described with respect to Figs. 3 and 5 . That is to say, the first percentage is 0.1 percent in this case.
  • FIGs. 8 and 9 illustrate the first and, respectively, second stage of the entire fitting process
  • the third and final stage of the entire fitting process in this embodiment is illustrated in Fig. 10 .
  • the part of the fitting process carried out at the second stage comprised a fitting of the provided functional prototype to the data values SR2" that used, as initial values for the one or more fit parameters of the functional prototype, the final values resulting from the first stage, i.e. the fit parameters describing the fit function FR1".
  • the time needed for fitting at the second stage could be performed in less time as compared to the time that would have been needed if the functional prototype would have been fitted to the data values SR2" based on generic initial values for the fit parameters.
  • the fitting at the third stage required less time than what would have been needed for fitting the functional prototype to the data values SR3", which correspond to the entire set of data values S0", based on generic initial values for the fit parameters.
  • the indicator determination unit 102 may be configured such that the fitting is carried out in several stages, wherein, after the fitting has been carried out in a stage, it is determined whether the fitted fit parameter is in conformance with an abort criterion, wherein, if this is the case, the fitting is terminated.
  • the fitting may only be terminated if more than one or all of the fit parameters are in conformance with respective abort criteria.
  • abort criteria may be used to terminate fitting if the current fit function is already accurate enough, abort criteria may also be used to terminate fitting if it is concluded, for instance, after one of the stages of a multiple staged fitting process, that the fit results will likely not lead to any valid fit function at all.
  • the fit parameter corresponding to the position of the maximum of the function may be expected to lie in a certain range for a given set of data values.
  • the range may be defined by one or both of an upper limit and a lower limit, wherein the lower limit may correspond to the position of the first data value and the upper limit may correspond to a multiple of a signal length, the multiple being, for instance, any value from 1 to 5.
  • the fitting is carried out in three stages, wherein at the first stage 1/250 of the of the data values are used, at the second stage 1 / 250 of the data values are used and at the third stage all data values are to be used, the final stage does not need to be calculated if the position of the maximum of the fitted function resulting from the second stage is not within the previously defined range, because it can be expected that the final result, particularly the final position of the maximum resulting from the final stage, will not be a valid result.
  • Fig. 12 schematically and exemplarily illustrates a determination procedure 1100 according to a combination of previously described embodiments.
  • a first step 1101 data values s0, S0, S0', S0" are determined based on a sequence of measured pulse signals that has been provided.
  • initial values are determined for each of the one or more fit parameters of the functional prototype that is to be fitted to the data values s0, S0, S0', S0".
  • the initial values are determined based on characteristics of the data values s0, S0, S0', S0" that correspond to characteristics of the functional prototype for which the respective fit parameters are indicative.
  • a height of a maximum of the data values s0, S0, S0', S0" and a position of the maximum of the data values s0, S0, S0', S0" can be determined as initial values for the first and the second fit parameter.
  • multiple staged fitting is carried out using the initial values determined in step 1102 as initial values for the fit parameters in the first stage.
  • the multiple staged fitting carried out in step 1103 can, for instance, be carried out as described with reference to Figs. 7 to 10 .
  • a further set of data values needs to be considered for a further fitting, wherein the further set of data values would be determined based on a further sequence of measured pulse signals.
  • This decision may be made separately for the set of first data values determined for a first fitting and the set of second data values determined for a second fitting performed for a same sequence of measured pulse signals.
  • the decision whether more fitting is needed may be made in accordance with an abort criterion similar to the one possibly employed for determining, after any of the non-final stages of the multiple staged fitting, whether the multiple staged fitting is to be terminated.
  • step 1104 If it is decided in step 1104 that further fitting is needed, initial fit parameters for the further fitting are estimated in step 1105 based on the final fit parameters resulting from the previous fitting performed in step 1103.
  • Multiple staged fitting as performed in step 1103 is then repeated for the further set of data values, with the exception that, as initial values for the fit parameters in the first stage of the multiple staged fitting carried out in the repeated step 1105, no values describing characteristics of the further data values corresponding to characteristics of the functional prototype for which the fit parameters are indicative are used, but, instead, the final values of the fit parameters resulting from the multiple staged fitting to the previous set of data values are used.
  • the indicator determination unit 102 determines an initial value for a fit parameter based on a value of the fit parameter known from a previous fitting.
  • the cycle defined by steps 1103, 1104 and 1105 is repeated as long as it is decided in step 1104 or its repetitions that more fitting is needed, i.e., for instance, as long as an abort criterion is not met, and as long as not all provided sets of data values have been used for fitting. For instance, if four sequences of measured pulse signals have been provided by the pulse signals providing unit, or a single sequence of measured pulse signals corresponding to four inflations of the pressure cuff, and if, based on the provided sequences or, respectively, sequence, four sets of data values have been determined, each of which corresponds to one sequence or, respectively, inflation, then the cycle defined by steps 1103, 1104 and 1105 will be carried out a maximum number of four times.
  • step 1106 Whenever it is decided that no more fitting is needed, i.e., in step 1104 or its repetitions, a step 1106 follows.
  • the indicator that is representative for the physiological parameter is determined based on the final fit function resulting from the previous, i.e. last, fitting.
  • the apparatus 100 configured as described above with reference to Fig. 12 is used on a set of 540 representative records for calculating a PPV, an SVV or an SPV, wherein for each record 11 individual fits are calculated for 2 different kinds of input samples and wherein a bell-shaped Cauchy-Lorentz function is used as functional prototype, the following performance optimization can be observed.
  • the time required for the fitting of all records amounts to 64 seconds in total, the average time required for the fitting per record being 0.119 seconds. For each record, an average of 49 fitting iterations on all of the approximately 13000 input samples can be observed.
  • the final fit parameters of a previous fit were used as initial fit parameters of the fit to be done.
  • the time required for the fitting of all records amounts to 24 seconds, the average time required for the fitting per record being 0.044s, which is approximately 38 percent of the time required in the non-optimized case.
  • an average of 38 fitting iterations on approximately 50 reduced input samples, an average of 23 fitting iterations on approximately 800 reduced input samples and an average of 11 fitting iterations on all approximately 13000 input samples can be observed.
  • Fig. 13 illustrates schematically and exemplarily an embodiment of a method 200 for determining an indicator that is representative of a physiological parameter, such as a patient's volume responsiveness.
  • the method comprises, in a first part 201, providing a sequence of measured pulse signals of a patient over a time period corresponding to a plurality of subsequent respiratory cycles of the patient by a pulse signals providing unit 101, wherein the sequence of measured pulse signals has been detected by a non-invasive pulse measurement method.
  • the measurement might be performed using a pressure cuff, wherein the pressure applied in the pressure cuff is continuously increased or decreased in the course of a measurement period.
  • the method 200 comprises carrying out a determination procedure adapted to determine an indicator that is representative for the physiological parameter based on the provided sequence of pulse signals by an indicator determination unit 102, wherein the determination procedure includes a determination of data values s0, S0, S0', S0" based on the provided sequence of pulse signals and a fitting of a provided functional prototype, which depends on a fit parameter to be modified during the fitting, to the determined data values s0, S0, S0', S0".
  • the indicator determination unit 102 is configured such that a), before the fitting, the number of data values s0, S0, S0', S0" is reduced for the entire fitting, and/or b) an initial value is determined for the fit parameter based on characteristics of the data values s0, S0, S0', S0" which are related to characteristics of the functional prototype and/or based on a value of the fit parameter known from a previous fitting, and/or c) the fitting is carried out in several stages, wherein the number of data values used for fitting is increased from stage to stage and wherein in a current stage a value of a fit parameter determined in a previous stage is used as initial value for the fit parameter in the current stage.
  • the determination procedure carried out in the second part 202 of the method 200 can be a determination procedure of the kind described above, for example, with reference to Fig. 12 .
  • the apparatus 100 has mainly been described as being an apparatus for determining an indicator that is representative for a patient's volume responsiveness
  • the apparatus may also be an apparatus for determining an indicator that is representative for a patient's blood pressure, for instance.
  • the data values to which the functional prototype is fitted may correspond to a tissue pressure waveform (TPW) curve like the one disclosed in WO 2018/210931 A1 , wherein the functional prototype may be chosen such that it reflects a known basic form of such a TPW curve.
  • TPW tissue pressure waveform
  • the apparatus may be understood as being an apparatus for optimizing the runtime of fitting algorithms.
  • the pulse signals providing unit may be configured to provide any sequence of signals that has been detected, wherein the indicator determination unit 102 may be configured to determine data values based on the provided sequence of signals and to fit a provided functional prototype to the determined data values.
  • the functional prototype was a bell-shaped function, other functional prototypes may be used.
  • the functional prototype may be provided depending on a type of the detected sequence of signals based on which the data values are determined to which the functional prototype is to be fitted.
  • the runtime of a fitting algorithm can be reduced compared to simply applying the fitting algorithm to the entire set of data values by reducing, for instance, the total number of iterations in the fitting algorithm or the number of data values actually used for fitting.
  • the reduction is preferably done in a way that keeps the main characteristics of the data values in relation to the functional prototype.
  • the total number of iterations needed, and therefore the total time required for fitting can be reduced, for instance, by initially predicting parameters of the functional prototype that are preferably close to the final parameters to be determined by the fitting.
  • the final fit parameters of a previous fit can be used as initial fit parameters of the fit to be done.
  • the sets of data values may correspond to subsequent breathing or ventilation cycles of the patient, in the more general context of reducing the runtime of fitting algorithms, the different sets of data values may comprise similar characteristics due to other reasons than the cyclic nature of breathing or ventilation.
  • the runtime of a general fitting algorithm can also be reduced by iteratively, i.e. in stages, applying the fitting algorithm, starting with the reduced number of data values and increasing the number of data values in each subsequent stage, wherein, preferably, in the last stage the original, entire number of data values is used for fitting.
  • the preferred number of stages may be found depending on the intended application by, for instance, finding the minimum of the total runtime experimentally. In this way, the optimal combination of the number of data values used for the fitting in each stage and the number of stages can be found.
  • a single unit or device may fulfill the functions of several items recited in the claims.
  • the mere fact that certain measures are recited in mutually different dependent claims does not indicate that a combination of these measures cannot be used to advantage.
  • These procedures can be implemented as program code means of a computer program and/or as dedicated hardware, particularly as an embedded system.
  • a computer program may be stored/distributed on a suitable medium, such as an optical storage medium or a solid-state medium, supplied together with or as part of other hardware, but may also be distributed in other forms, such as via the Internet or other wired or wireless telecommunication systems.
  • a suitable medium such as an optical storage medium or a solid-state medium, supplied together with or as part of other hardware, but may also be distributed in other forms, such as via the Internet or other wired or wireless telecommunication systems.
  • the invention relates to an apparatus for determining an indicator that is representative of a physiological parameter like a fluid responsiveness parameter.
  • the indicator is determined based on a fitting of a functional prototype to data values determined from pulse signals measured over subsequent respiratory cycles.
  • the fitting process is accelerated by a) reducing the number of data values before fitting, b) determining an initial fit parameter value for the functional prototype based on characteristics of the data values and/or a fit parameter value known from a previous fitting, and/or c) carrying out the fitting in several stages, wherein the number of data values used is increased from stage to stage and a fit parameter value determined in a previous stage is used as initial fit parameter value in a current stage. This allows for a faster determination of, for instance, a fluid responsiveness parameter.
EP20211323.9A 2020-12-02 2020-12-02 Systeme permettant de determiner un indicateur representatif d'un parametre physiologique Withdrawn EP4008245A1 (fr)

Priority Applications (6)

Application Number Priority Date Filing Date Title
EP20211323.9A EP4008245A1 (fr) 2020-12-02 2020-12-02 Systeme permettant de determiner un indicateur representatif d'un parametre physiologique
CN202180081464.8A CN116583219A (zh) 2020-12-02 2021-11-29 用于确定表示流体响应性参数的指标的装置
JP2023533256A JP2023551546A (ja) 2020-12-02 2021-11-29 流体反応性パラメータを表す指標を決定するための装置
EP21819119.5A EP4266992A1 (fr) 2020-12-02 2021-11-29 Appareil pour déterminer un indicateur représentatif d'un paramètre de réactivité de fluide
US18/039,580 US20240090781A1 (en) 2020-12-02 2021-11-29 Apparatus for determining an indicator representative for a fluid responsiveness parameter
PCT/EP2021/083267 WO2022117471A1 (fr) 2020-12-02 2021-11-29 Appareil pour déterminer un indicateur représentatif d'un paramètre de réactivité de fluide

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
EP20211323.9A EP4008245A1 (fr) 2020-12-02 2020-12-02 Systeme permettant de determiner un indicateur representatif d'un parametre physiologique

Publications (1)

Publication Number Publication Date
EP4008245A1 true EP4008245A1 (fr) 2022-06-08

Family

ID=73694901

Family Applications (2)

Application Number Title Priority Date Filing Date
EP20211323.9A Withdrawn EP4008245A1 (fr) 2020-12-02 2020-12-02 Systeme permettant de determiner un indicateur representatif d'un parametre physiologique
EP21819119.5A Pending EP4266992A1 (fr) 2020-12-02 2021-11-29 Appareil pour déterminer un indicateur représentatif d'un paramètre de réactivité de fluide

Family Applications After (1)

Application Number Title Priority Date Filing Date
EP21819119.5A Pending EP4266992A1 (fr) 2020-12-02 2021-11-29 Appareil pour déterminer un indicateur représentatif d'un paramètre de réactivité de fluide

Country Status (5)

Country Link
US (1) US20240090781A1 (fr)
EP (2) EP4008245A1 (fr)
JP (1) JP2023551546A (fr)
CN (1) CN116583219A (fr)
WO (1) WO2022117471A1 (fr)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5577508A (en) * 1995-01-13 1996-11-26 Johnson & Johnson Medical, Inc. Determination of oscillometric blood pressure by linear approximation
US5704362A (en) * 1993-08-13 1998-01-06 Johnson & Johnson Medical, Inc. Method for oscillometric blood pressure determination employing curve fitting
WO2014121945A1 (fr) 2013-02-08 2014-08-14 Up-Med Gmbh Système de mesure de la pression sanguine comprenant une partie rigide empêchant la formation de plis
CN105342590A (zh) * 2015-11-02 2016-02-24 浙江漫思网络科技有限公司 一种基于示波法的血压包络的估计方法
EP2759257B1 (fr) 2013-01-25 2016-09-14 UP-MED GmbH Procédé, unité logique et système permettant de déterminer un paramètre représentatif de la réactivité de volume du patient
WO2018210931A1 (fr) 2017-05-17 2018-11-22 Up-Med Gmbh Procédé pour la détermination non invasive d'au moins une valeur de pression artérielle, dispositif de mesure et système pour la détermination non invasive de la pression artérielle

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5704362A (en) * 1993-08-13 1998-01-06 Johnson & Johnson Medical, Inc. Method for oscillometric blood pressure determination employing curve fitting
US5577508A (en) * 1995-01-13 1996-11-26 Johnson & Johnson Medical, Inc. Determination of oscillometric blood pressure by linear approximation
EP2759257B1 (fr) 2013-01-25 2016-09-14 UP-MED GmbH Procédé, unité logique et système permettant de déterminer un paramètre représentatif de la réactivité de volume du patient
WO2014121945A1 (fr) 2013-02-08 2014-08-14 Up-Med Gmbh Système de mesure de la pression sanguine comprenant une partie rigide empêchant la formation de plis
CN105342590A (zh) * 2015-11-02 2016-02-24 浙江漫思网络科技有限公司 一种基于示波法的血压包络的估计方法
WO2018210931A1 (fr) 2017-05-17 2018-11-22 Up-Med Gmbh Procédé pour la détermination non invasive d'au moins une valeur de pression artérielle, dispositif de mesure et système pour la détermination non invasive de la pression artérielle

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
POOI LIM ET AL: "Improved Measurement of Blood Pressure by Extraction of Characteristic Features from the Cuff Oscillometric Waveform", SENSORS, vol. 15, no. 6, 16 June 2015 (2015-06-16), pages 14142 - 14161, XP055472873, DOI: 10.3390/s150614142 *

Also Published As

Publication number Publication date
US20240090781A1 (en) 2024-03-21
EP4266992A1 (fr) 2023-11-01
CN116583219A (zh) 2023-08-11
JP2023551546A (ja) 2023-12-08
WO2022117471A1 (fr) 2022-06-09

Similar Documents

Publication Publication Date Title
EP3654837B1 (fr) Appareil et procédé permettant de déterminer la pression sanguine d'un sujet
US11298033B2 (en) Non-invasive blood pressure monitors and methods of operating the same
CN107920758B (zh) 用于确定血压值的方法、装置和计算机程序
EP2759257B1 (fr) Procédé, unité logique et système permettant de déterminer un paramètre représentatif de la réactivité de volume du patient
US11406305B2 (en) Signal quality metric for cardiovascular time series
US20120157791A1 (en) Adaptive time domain filtering for improved blood pressure estimation
US20070129638A1 (en) Method and system for cuff pressure reversions
WO2017129495A1 (fr) Module et procédé de mesure de fréquence du pouls
CN110573067A (zh) 无创肱动脉血压测量
EP4008245A1 (fr) Systeme permettant de determiner un indicateur representatif d'un parametre physiologique
US20220202298A1 (en) Unsupervised real-time classification for arterial blood pressure signals
WO2023072730A1 (fr) Dispositif, système et procédé d'étalonnage d'un substitut de tension artérielle pour une utilisation dans la surveillance de la tension artérielle d'un sujet
US20090221924A1 (en) Method and system for non-invasive blood pressure determination
EP4252633A1 (fr) Appareil permettant de déterminer un indicateur représentatif d'un paramétre physiologique
GB2461164A (en) Oscillatory blood pressure estimation with artefact analysis
EP4173556A1 (fr) Dispositif, système et procédé d'étalonnage d'un substitut de pression sanguine pour une utilisation dans la surveillance de la pression sanguine d'un sujet
EP4014837A1 (fr) Procédé, appareil et produit programme informatique pour l'analyse d'un signal d'onde d'impulsion
US20220386964A1 (en) Processing physiological signals to determine health-related information
JP2014180514A (ja) 呼吸情報推定装置及び方法並びにプログラム
CN117015338A (zh) 动脉压推定装置、动脉压推定系统及动脉压推定方法

Legal Events

Date Code Title Description
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: THE APPLICATION HAS BEEN PUBLISHED

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

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: 20221209