text: Patent US20050027205 - Combining measurements from breathing rate sensors - Google PatentsSearch Images Maps Play YouTube News Gmail Drive More »Sign inAdvanced Patent SearchPatentsA method for combining breathing rate measurements from two or more independent measurement channels. Independent measurements of breathing rate, for instance by impedance pneumography and photoplethysmography, can be combined to derive an improved measurement eliminating artefacts on one channel. A...http://www.google.com/patents/US20050027205?utm_source=gb-gplus-sharePatent US20050027205 - Combining measurements from breathing rate sensorsAdvanced Patent SearchPublication numberUS20050027205 A1Publication typeApplicationApplication numberUS 10/498,673PCT numberPCT/GB2002/005684Publication dateFeb 3, 2005Filing dateDec 13, 2002Priority dateDec 14, 2001Also published asDE60235917D1, EP1463447A1, EP1463447B1, US7318808, WO2003051198A1Publication number10498673, 498673, PCT/2002/5684, PCT/GB/2/005684, PCT/GB/2/05684, PCT/GB/2002/005684, PCT/GB/2002/05684, PCT/GB2/005684, PCT/GB2/05684, PCT/GB2002/005684, PCT/GB2002/05684, PCT/GB2002005684, PCT/GB200205684, PCT/GB2005684, PCT/GB205684, US 2005/0027205 A1, US 2005/027205 A1, US 20050027205 A1, US 20050027205A1, US 2005027205 A1, US 2005027205A1, US-A1-20050027205, US-A1-2005027205, US2005/0027205A1, US2005/027205A1, US20050027205 A1, US20050027205A1, US2005027205 A1, US2005027205A1InventorsLionel Tarassenko, Catherine MasonOriginal AssigneeLionel Tarassenko, Mason Catherine LauraExport CitationBiBTeX, EndNote, RefManPatent Citations (13), Referenced by (41), Classifications (4), Legal Events (5) External Links: USPTO, USPTO Assignment, EspacenetCombining measurements from breathing rate sensors
  The present invention relates to a method and apparatus for measuring the breathing rate of a subject, such as a human or animal, and in particular to a way of combining measurements from two or more breathing rate sensors in order to provide an improved measurement of breathing rate. There is no clinically acceptable method for the non-invasive measurement of breathing rate. A nasal thermistor provides a simple and inexpensive means of tracking respiration but it is obtrusive and is not deemed to be acceptable for patients on the general ward or even the Coronary Care Unit (CCU). Electrical impedance plethysmography (also known as impedance pneumography) can be used to provide an indirect measure of respiration by measuring the changes in electrical impedance across the chest with breathing. In this method a small-amplitude, high-frequency current is injected into the body through a pair of surface electrodes and the resulting voltage is demodulated to obtain impedance measurements. The electrical impedance increases as high-resistivity air enters the lungs during inspiration but part of the change is also due to the movement of the electrodes on the chest wall. When monitored in a clinical environment, the impedance plethysmography (IP) signal is often very noisy and is seriously disrupted by patient movement or change in posture. As a consequence, it has not been considered reliable enough to provide respiration information for regular use on the ward. Respiratory information is found in other signals recorded from patients with non-invasive sensors. For example, both the electrocardiogram (EGG) and photoplethysmogram (PPG) waveforms are modulated by the patient's breathing. The PPG signal represents the variation in light absorption across a finger or earlobe with every heart beat. This signal is measured at two wavelengths (usually in the red and near infra-red parts of the spectrum) in a standard pulse oximeter. FIG. 2 shows 30-second sections of IP (FIG. 2A) and PPG (FIG. 2B) signals recorded from a patient in CCU. Referring first to the IP signal in FIG. 2A, the main peaks I, each corresponding to a breath, can be identified. However, changes in electrical impedance with the heart beat (as opposed to respiration) are apparent between t=1607 and t=1 608 (marked as "A") and probably conceal a peak caused by a breath. Referring to FIG. 2B, the modulation of the PPG caused by respiration (marked "I") can be seen, although there are small movement artefacts B, C at t=1617 and t=1629. FIG. 3 shows the IP and PPG signals from the same patient, for a five-minute period. The effect, marked D, E, of movement artefact on the finger probe from which the PPG signal is recorded becomes very obvious at t=1770 and t=1920. Clearly the artefacts caused by movement and heartbeat would affect the measurement of breathing rate from these signals. One might consider removing these artefacts by some type of filtering and thresholding. FIG. 4B shows the effect of applying a Finite Impulse Response (FIR) low-pass filter with a pass-band cut-off at 0.33 Hz (corresponding to 20 breaths per minute) and a stop-band cut-off of 1 Hz (−50 dB) to the IP waveform of FIG. 4A. The cardiac-synchronous changes are filtered out (see the output of the filter between t=1605 and t=1610) and the respiratory cycle is clear. There are as many peaks (i.e. breaths) in FIG. 4B as there are in FIG. 4D, which shows the result of tracking the peaks of the PPG signal shown in FIG. 4C and interpolating (with straight lines) between each peak. The modulation envelope picked out by the peak tracking produces another respiratory waveform. The same information is again displayed on a five-minute timescale in FIG. 5. The breathing rate can be estimated by calculating the interval between two successive peaks of these waveforms, inverting the result and multiplying it by a factor of 60 in order to obtain an estimate of breathing rate in breaths per minute. The results of these computations over the five-minute period are shown in FIGS. 5E and 5F respectively. The breathing rate is approximately 18 breaths/minute throughout the period, but the occurrence of peaks caused by cardiac as opposed to respiratory changes in the IP signal of FIGS. 5A and 5B gives rise to erroneously high breathing rates at t=1720, 1730, 1780, 1795, 1830 and 1860. Unfortunately no amount of optimisation of the FIR low-pass filter characteristics will ever completely remove the cardiac-synchronous information which is occasionally very prominent. For instance, in a hyperventilating patient, the breathing rate may be as high as 40 breaths per minute, which is similar to a slow heart rate. So separating signal from noise on a fixed basis is impossible. Similarly, the two major instances D, E of movement artefact in the PPG signal of FIGS. 5C and 5D at t=1770 and t=1920 cause erroneously low estimates of breathing rates at D′ and E′ because there is a significant delay between the last valid peak and the first peak after the movement artefact. The estimate of breathing rate from the PPG signal is also more variable because occurrences of even slight movement artefact affect the tracking of the peaks of the PPG signal. The present invention provides a method and apparatus for improving measurement of breathing rate by combining two measurements of it in a way which allows valid changes in the breathing rate to be distinguished from artefacts. In accordance with the invention two measurements of breathing rate made in different ways are combined with weights based on the amount of "confidence" in the measurement, to give an improved measurement or estimate of the actual breathing rate. The two measurements may, for example, be obtained using impedance pneumography and photoplethysmography, though other signals which include respiratory information (for instance ECG) can also be used instead of either of these signals, or in addition to them. In more detail the invention provides a method of measuring breathing rate of a subject comprising the steps of: predicting the value of each of two independent measurements of the breathing rate, making two independent measurements of the breathing rate to produce two measured values, calculating the respective differences between the predicted values and the measured values, and combining the two measured values with weights determined by said differences. The steps of prediction, measurement, calculation and combination may be repeated continuously. The predicted value for each of the independent measurements may be based on the preceding predicted value and the difference between the preceding predicted value and the preceding measurement. The predicted value for each of the measurements may be calculated using a linear or non-linear predictive model, and the model may be adaptive, adapting in dependence upon the amount of process noise in the measurements. In the combining of the two measured values, the weight of each value may vary inversely with the square of the difference between the predicted value and the measurement. In one example the two measured values may be combined according to the formula:— BR = BR 1 σ 2 2 σ 1 2 + σ 2 2 + BR 2 σ 1 2 σ 1 2 + σ 2 2 where BR1 and BR2 are the two measured values of breathing rate and σ1 and σ2 are the differences between the two measured values and their respective predicted values. The predicted values for the respective measurements may be based on respective models of the system, and the models may include estimates for process noise and sensor noise. The respective models may be mutually independent and may include the same estimates for process noise and sensor noise. The models may be Kalman filters. Measurements for which the differences between both measured values and their predicted values exceed a predetermined threshold may be discarded. Further, artefacts, for instance caused by movement or heartbeat in the measurements may be identified based on the values of the differences between the measured values and their predicted values. This identification may be used to discard sections of the signal. Thus with the present invention a prediction is made for each breathing rate measurement and the actual measurement is compared with its prediction. The difference is computed, which is termed the "innovation", and this innovation is used to calculate a weight which will be given to that measurement when it is combined with the other measurement, also weighted according to its innovation. The weights are calculated so that if the innovation on one measurement channel is high, whereas the innovation on the other measurement channel is low, the measurement from the low innovation channel is more heavily weighted. This is because a high level of innovation from one channel coinciding with a low innovation on the other charnel is regarded as indicative of an artefact on the higher innovation channel. It will be appreciated that the invention can be embodied using computer software, and thus the invention extends to a computer program for controlling and executing the method, or parts of it, and to a computer readable storage medium carrying the program. The invention also extends to corresponding apparatus for carrying out the method.
  An embodiment of the invention will now be described in which the invention is applied in the medical field for the measurement of breathing rate using the two signals described above, namely the impedance pneumography and photoplethysmography signals. The results of applying this to the data of FIGS. 2 to 5 will be illustrated. FIG. 1 schematically illustrates a breathing rate measurement apparatus in accordance with an embodiment of the invention. A PPG signal is obtained, as illustrated from a finger of a patient, though the signal can also be obtained from the earlobe, using a conventional PPG sensor 1 which is driven by and supplies its output to PPG apparatus 3. An IP signal is obtained by conventional IP equipment 7 using two electrodes 5 (though four electrodes may be used, with a pair for current injection and a pair for voltage measurement). The IP apparatus 7 and PPG apparatus 3 supply their signals to a processor 10 which processes the signals in the manner to be described below, and outputs a display of breathing rate on display 12. The display may be in graphical or numeric form, and may be combined with a display of others of the patient's vital signs. It will be appreciated that in this embodiment the IP apparatus 7 and PPG apparatus 3 produce an output to the processor 10 of raw signals of the types illustrated in FIGS. 6A and C. The processor 10 processes these signals to produce the signals of FIGS. 6B and D to G. However, the functions of the parts of the apparatus shown in FIG. 1 may be combined in different ways. The processor 10 in essence runs a Kalman filter on each breathing rate signal, and then fuses the filtered signals to produce a fused estimate of the breathing rate. The Kalman filter is a generic framework for analysing linear dynamical systems. Using the process model, the next state xt+1 is computed from the current state xt using the transition matrix A, assuming first-order (Markov) dynamics. The observation model relates the observations y1 to the state x1 of the system via the observation matrix C. The process and observation noise are assumed to be independent and to have zero-mean, Gaussian probability distributions. In the normal physiological state, the breathing rate can be modelled with a scalar Kalman filter which assumes that the time to the next breath is the same as the interval since the previous breath plus some process noise characterising normal or breath-to-breath variability. Thus: x t+l =Ax t +w(Process model) y t =Cx t+ν(Observation model) In this embodiment it is assumed that C=i, i.e. the breathing rate is both the measurement and the state describing the process, with ν as the sensor noise model. It is also assumed that the next breathing rate only varies by a small amount with respect to the current one, i.e. A=1 and w is the variance associated with breath-by-breath variability. A scalar Kalman filter is run separately on each breathing rate signal, using the same process and measurement noise models for both channels. For each channel, the next state xpred is predicted using the Kalman filter equation (in which the next state is equal to the previous state plus the Kalman gain times the innovation). In this embodiment, the next measurement ypred is the same as the next state xpred (since C=1) and, after the next measurement of breathing rate, yt+1, the innovation εt+1 is computed, where: y t+1 =y pred+εt+1 εt+1 the difference between the actual value and the predicted value, should normally be a zero-mean white noise sequence. The square of the innovation, or variance, σt+1 2=εt+1 2, is the inverse of the "confidence" which is associated with the prediction. A robust estimate of the breathing rate is now obtained by mixing the two breathing rate estimates in inverse proportion to the variance associated with each one: BR = σ 2 2 σ 1 2 + σ 2 2 BR1 + σ 1 2 σ 1 2 + σ 2 2 BR2 where BR1=yt+1 and σ1 2=σt+1 2 for channel 1 (filtered IP), with equivalent expressions for BR2 and channel 2 (PPG peaks). The innovation variance and the Kalman gain are also calculated to update the state estimate in readiness for the next cycle. Ai example of an implementation of this model in MATLAB is given in Appendix 1. That example is general and will work for vector quantities though in this embodiment the quantities are scalar. It can be seen from Appendix 1 that the predicted breathing rate for each new measurement cycle (xnew) is equal to the previously predicted value (xpred) plus the Kalman gain K times the innovation e. The Kalman gain K is derived from the predicted variance Vpred and the measurement variance R. The predicted variance is derived from the previous predicted variance and the process noise Q. To start the process off it is initialised using an initial value of the breathing rate as 15 and an initial value of the state variance of 40. The process noise in the Q in this embodiment is set to 10 and the measurement noise variance R is set to 100. It will be clear from the implementation that, as normal with a Kalman filter, the variance and Kalman gain are not dependent on the measurement values. The measurement values are only used in the new prediction of breathing rate via the innovation e. Thus it will be noted that for the constant values of Q and R used in this example the Kalman gain K tends to about 0.2 and the state co-variance V tends to about 27. However, K can be made adaptive by modifying the values for the variance constants. Q and R, preferably the process variance Q, according to the type of process being encountered. The advantages of this method can be appreciated from four possible contexts: 1. normal breathing: low innovations on both channels; both measurements BR1 and BR2 are weighted equally. 2. valid change seen on both channels: the subject begins to breathe more quickly or more slowly due to a physiological change; although the innovation is high, it is high for both channels, and so both (valid) measurements are again weighted equally. 3. artefact on one channel: high innovation on one channel only; the information from that channel is ignored because it is given a low weighting (high variance). 4. artefact on both channels: the information is corrupted on both channels. Prolonged movement artefact, however, is characterised by high values of innovation on both channels for a sustained period of time and this can be the basis for discarding sections of data corrupted by movement artefact. Note that a fused estimate is computed every time a new measurement of BR1 or BR2 becomes available, i.e. twice during the respiration cycle (once for each sensor). It should also be noted that in the event of a problem with an electrode such as it becoming detached or falling off, or of signal loss through some other cause, this would result in a high innovation and thus a low weighting for that channel. FIG. 6G illustrates the result of combining the breathing rate signals of FIGS. 6E and F in this way. FIG. 6G illustrates the result of applying the invention to the breathing rate signals of FIGS. 6E and F obtained respectively from the impedence pneumography and PPG signals of FIGS. 6A to D. Thus the original IP signal of FIG. 6A is filtered using a finite impulse response low-pass filter to produce the processed signal of FIG. 6B. The PPG signal of FIG. 6C is subjected to peak tracking as described above to produce a respiratory waveform of FIG. 6D. Each of these is used to estimate a breathing rate, by calculating the interval between two successive peaks, inverting the result and multiplying it by a factor of 60. These two breathing rate estimates are shown in FIGS. 6E and F. The two breathing rate estimates BR1 and BR2 are subjected to the Kalman filtering and fusing the process described above and this gives the fused breathing rate estimate of FIG. 6G. It can be seen that apparent changes in breathing rate which appear only on one channel (for instance marked F on FIG. 6E or the changes D′ and E′ on FIG. 6F) are removed from the fused breathing rate estimate. However, a change in breathing rate on both channels, such as marked G in FIGS. 6E and 6F, appears as G′ on the fused estimate of FIG. 6G. Appendix 1 load -ascii ip_resprate; data_file = ip_resprate; time = data_file(:,1); br = data_file(:,2); start = 1 stop = size(data_file),1); fprintf('number of breaths detected = %d \n', stop); br_limit = 40; % X(t+1) = A X(t)+noise(Q) - process model with Q as variance of noise w % Y(t) = C X(t)+noise(R) - measurement model with R as variance of noise ν ss = 1; % state size - sets to one dimensional, ie scalar though routine works for vectors os = 1; % observation size - sets to one dimensional, ie scalar A = [1]; % assume x(t+1)=x(t) C = [1]; % assume y=x Q = 10.0*eye(ss); % process noise - eye is the identity matrix in MATLAB -here just unity R = 100.0*eye(os); % measurement noise variance initx = [15]; % initial state value(BR of 15 bpm) initV = 40*eye(ss); % initial state variance xnew = initx; % - initialisation Vnew = initV; no_update = 0; n = 1; for i = start:stop % - start of cycle x = xnew; % update from previous cycle V = Vnew; % update from previous cycle xpred = A*x; % prediction of state Vpred = A*V*A′ + Q; % prediction of state covariance, A′ is transpose of A ypred(i) = C*xpred; % prediction of measurement y(i) = br(i); % "make measurement" t(i) = time(i); % if y(i) == br_limit no_update = 1; %artefact % else no_update = 0; % end e = y(i) − ypred(i); % calculate innovation innov(i) = e; % for plotting sigma2(i) = e*e; % variance for saving S = C*Vpred*C′ + R; % innovation covariance Sinv = inv(S); % invert S to compute Kalman gain K = Vpred*C′ *Sinv; % Kalman gain matrix if no_update == 0 % only update if no artefact xnew = xpred + K*e; % update state by the innovation controlled by the Kalman gain Vnew = (eye(ss) − K*C)*Vpred; % update state covariance end if i − (n*10) == 0 fprintf( {grave over ( )}\n %d0 breaths processed {grave over ( )}, n); n = n + 1; end end Patent CitationsCited PatentFiling datePublication dateApplicantTitleUS4791573 *Sep 2, 1986Dec 13, 1988Sanders Associates, Inc.State-deviation-estimation circuit employing a phase-locked-loop phase referenceUS4949710 *Oct 6, 1988Aug 21, 1990Protocol Systems, Inc.Method of artifact rejection for noninvasive blood-pressure measurement by prediction and adjustment of blood-pressure dataUS5355893 *Sep 13, 1993Oct 18, 1994Mick Peter RVital signs monitorUS5396893 *Feb 14, 1991Mar 14, 1995Oberg; Ake P.Method and apparatus for analyzing heart and respiratory frequencies photoplethysmographicallyUS5626140 *Nov 1, 1995May 6, 1997Spacelabs Medical, Inc.System and method of multi-sensor fusion of physiological measurementsUS5853364 *Jun 7, 1996Dec 29, 1998Nellcor Puritan Bennett, Inc.Method and apparatus for estimating physiological parameters using model-based adaptive filteringUS5919141 *Sep 13, 1996Jul 6, 1999Life Sensing Instrument Company, Inc.Vital sign remote monitoring deviceUS6299582 *Mar 18, 1997Oct 9, 2001Data Sciences International, Inc.Respiration monitoring system based on sensed blood pressure variationsUS6839659 *Jun 8, 2001Jan 4, 2005Isis Innovation LimitedSystem and method for acquiring dataUS6928370 *Jul 5, 2001Aug 9, 2005Oxford Biosignals LimitedHealth monitoringUS20030187337 *Jun 8, 2001Oct 2, 2003Lionel TarassenkoCombining measurements from different sensorsUS20040148140 *May 29, 2002Jul 29, 2004Lionel TarassenkoPatient condition displayUS20050027205 *Dec 13, 2002Feb 3, 2005Lionel TarassenkoCombining measurements from breathing rate sensors* Cited by examinerReferenced byCiting PatentFiling datePublication dateApplicantTitleUS7318808Dec 13, 2002Jan 15, 2008Isis Innovation LimitedCombining measurements from breathing rate sensorsUS7647185Jun 8, 2001Jan 12, 2010Oxford Biosignals LimitedCombining measurements from different sensorsUS7978064Jul 12, 2011Proteus Biomedical, Inc.Communication system with partial power sourceUS8473306Apr 1, 2010Jun 25, 2013Ottawa Hospital Research InstituteMethod and apparatus for monitoring physiological parameter variability over time for one or more organsUS8515513Oct 30, 2009Aug 20, 2013Covidien LpSystem and method for facilitating observation of monitored physiologic dataUS8545417 *Sep 14, 2009Oct 1, 2013Sotera Wireless, Inc.Body-worn monitor for measuring respiration rateUS8616203Apr 20, 2012Dec 31, 2013Ric Investments, LlcNasal and oral patient interfacesUS8622922 *Sep 14, 2009Jan 7, 2014Sotera Wireless, Inc.Body-worn monitor for measuring respiration rateUS8674825Mar 13, 2009Mar 18, 2014Proteus Digital Health, Inc.Pharma-informatics systemUS8696588 *Oct 24, 2007Apr 15, 2014Dr�ger Medical GmbHDevice and method for determining a respiration rateUS8718193Nov 19, 2007May 6, 2014Proteus Digital Health, Inc.Active signal processing personal health signal receiversUS8730031Jul 11, 2011May 20, 2014Proteus Digital Health, Inc.Communication system using an implantable deviceUS8740802Dec 30, 2009Jun 3, 2014Sotera Wireless, Inc.Body-worn system for measuring continuous non-invasive blood pressure (cNIBP)US8740807Sep 14, 2009Jun 3, 2014Sotera Wireless, Inc.Body-worn monitor for measuring respiration rateUS8755871Nov 30, 2011Jun 17, 2014Covidien LpSystems and methods for detecting arrhythmia from a physiological signalUS8795170Nov 29, 2006Aug 5, 2014Venture Gain LLCResidual based monitoring of human healthUS8808188Dec 30, 2009Aug 19, 2014Sotera Wireless, Inc.Body-worn system for measuring continuous non-invasive blood pressure (cNIBP)US8816847Jun 3, 2011Aug 26, 2014Proteus Digital Health, Inc.Communication system with partial power sourceUS8868453Nov 4, 2010Oct 21, 2014Proteus Digital Health, Inc.System for supply chain managementUS8888700 *Apr 19, 2010Nov 18, 2014Sotera Wireless, Inc.Body-worn monitor for measuring respiratory rateUS8912908Jul 11, 2011Dec 16, 2014Proteus Digital Health, Inc.Communication system with remote activationUS8924235Jun 11, 2013Dec 30, 2014Ottawa Hospital Research InstituteMethod and apparatus for monitoring physiological parameter variability over time for one or more organsUS8945005Oct 25, 2007Feb 3, 2015Proteus Digital Health, Inc.Controlled activation ingestible identifierUS8956287May 2, 2007Feb 17, 2015Proteus Digital Health, Inc.Patient customized therapeutic regimensUS9014779Jan 28, 2011Apr 21, 2015Proteus Digital Health, Inc.Data gathering systemUS9060708Jul 25, 2014Jun 23, 2015Proteus Digital Health, Inc.Multi-mode communication ingestible event markers and systems, and methods of using the sameUS9060746May 8, 2014Jun 23, 2015Covidien LpSystems and methods for detecting arrhythmia from a physiological signalUS9083589Mar 6, 2014Jul 14, 2015Proteus Digital Health, Inc.Active signal processing personal health signal receiversUS20050027205 *Dec 13, 2002Feb 3, 2005Lionel TarassenkoCombining measurements from breathing rate sensorsUS20080139955 *Oct 24, 2007Jun 12, 2008Drager Medical Ag & Co. KgDevice and method for determining a respiration rateUS20110066037 *Mar 17, 2011Matt BanetBody-worn monitor for measuring respiration rateUS20110066038 *Sep 14, 2009Mar 17, 2011Matt BanetBody-worn monitor for measuring respiration rateUS20110257489 *Apr 19, 2010Oct 20, 2011Sotera Wireless, Inc.Body-worn monitor for measuring respiratory rateUS20130109996 *May 2, 2013Case Western Reserve UniversityExpert system to facilitate source localization of brain electrical activityUSRE41236Jul 3, 2001Apr 20, 2010Seely Andrew J EMethod and apparatus for multiple patient parameter variability analysis and displayDE102006057709B4 *Dec 7, 2006Apr 2, 2015Dr�ger Medical GmbHVorrichtung und Verfahren zum Bestimmen einer AtemfrequenzDE102006060819A1Dec 21, 2006Jul 3, 2008Fresenius Medical Care Deutschland GmbhPatient's respiration rate determining method, involves determining momentary respiration rates, where weights of individual respiration rates depend on difference between respective respiration rates and estimated valueEP2359903A1 *Nov 19, 2009Aug 24, 2011Tohoku UniversitySignal processing device, signal processing method, signal processing program, computer-readable recording medium on which is recorded a signal processing program, and radiation therapy deviceWO2008019294A2Aug 2, 2007Feb 14, 2008Ric Investments LlcNasal and oral patient interfaceWO2013044073A1 *Sep 21, 2012Mar 28, 2013Nellcor Puritan Bennett IrelandSystems and methods for determining respiration information from a photoplethysmographWO2014149061A1 *Mar 22, 2013Sep 25, 2014Ipgdx, Llc (Michigan Company)Kalman algorithm electrical impedance device, and methods of making and using same* Cited by examinerClassifications U.S. Classification600/529International ClassificationA61B5/08Cooperative ClassificationA61B5/0816European ClassificationA61B5/08RLegal EventsDateCodeEventDescriptionJul 19, 2004ASAssignmentOwner name: ISIS INNOVATION LIMITED, UNITED STATESFree format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:TARASSENKO, LIONEL;MASON, CATHERINE LAURA;REEL/FRAME:016450/0345;SIGNING DATES FROM 20040602 TO 20040709May 20, 2008CCCertificate of correctionDec 1, 2009ASAssignmentOwner name: OXFORD BIOSIGNALS LIMITED, UNITED KINGDOMFree format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:ISIS INNOVATION LIMITED;REEL/FRAME:023594/0705Effective date: 20091126Jun 15, 2011FPAYFee paymentYear of fee payment: 4Jul 1, 2015FPAYFee paymentYear of fee payment: 8RotateOriginal ImageGoogle Home - Sitemap - USPTO Bulk Downloads - Privacy Policy - Terms of Service - About Google Patents - Send FeedbackData provided by IFI CLAIMS Patent Services