This invention relates generally to the field of medical devices and, more particularly, to adjusting parameters used in pulse oximetry analysis.
This section is intended to introduce the reader to various aspects of art that may be related to various aspects of the present invention, which are described and/or claimed below. This discussion is believed to be helpful in providing the reader with background information to facilitate a better understanding of the various aspects of the present invention. Accordingly, it should be understood that these statements are to be read in this light, and not as admissions of prior art.
A pulse oximeter is a medical device that may be used to measure various blood characteristics, for example, the oxygen saturation of hemoglobin in pulsing blood and/or the pulse rate of a patient. To measure these characteristics, a non-invasive sensor may be used to pass light through a portion of blood perfused tissue and photo-electrically sense the absorption and scattering of light in the tissue. The amount of light absorbed and/or scattered is analyzed to estimate the amount of blood constituent in the tissue.
A detector signal resulting from measurement of the light describes the blood characteristics. As an example, pulses refer to the varying amount of arterial blood present in the tissue during a cardiac cycle. The varying amount of arterial blood yields cyclic attenuation of the light passing through the tissue. Accordingly, the detector signal from measurement of the light exhibits the familiar plethysmographic waveform.
Analysis of detector signals involves processes that use various parameters. As an example, the analysis may involve filtering estimates of hemoglobin saturation to improve the accuracy of the saturation estimates. As another example, the analysis may involve filtering of plethysmographic waveforms. The filtering may use parameters such as filter weights or coefficients to adjust the filtering process. As another example, the analysis may involve applying pulse qualification criteria to qualify or disqualify pulses. The pulse qualification criteria may include parameters used to adjust the pulse qualification.
It is desirable to provide a flexible and robust methodology for adjusting the parameters of oximetry analysis.
Certain aspects commensurate in scope with the originally claimed invention are set forth below. It should be understood that these aspects are presented merely to provide the reader with a brief summary of certain forms the invention might take and that these aspects are not intended to limit the scope of the invention. Indeed, the invention may encompass a variety of aspects that may not be set forth below.
In accordance with one aspect of the present invention, there is provided a method for adjusting a pulse qualification criterion. The method may include receiving a signal representing a plurality of pulses, where the signal is generated in response to detecting light scattered from blood perfused tissue. A filter parameter value of a filter parameter of a filter may be determined, where the filter may be operable to filter the signal. A pulse qualification criterion may be adjusted in accordance with the filter parameter value, where the pulse qualification criterion may be used for qualifying a pulse. The pulses may be evaluated according to the pulse qualification criterion.
In accordance with another aspect of the present invention, there is provided a method for adjusting a pulse qualification criterion. The method may include receiving a signal representing a plurality of pulses, where the signal may be generated in response to detecting light scattered from blood perfused tissue. Each pulse may have an amplitude and a period, and a subset of the pulses may have a plurality of amplitudes and a plurality of periods. An average amplitude may be determined from the plurality of amplitudes, and an average period may be determined from the plurality of periods. A pulse qualification criterion may be adjusted in accordance with the average amplitude and the average period. Pulses may be evaluated according to the pulse qualification criterion.
In accordance with another aspect of the present invention, there is provided a method for adjusting a filter weight of a saturation filtering process. The method may include receiving a signal representing a plurality of pulses, where the signal may be generated in response to detecting light scattered from blood perfused tissue. The light may comprise a red waveform and an infrared waveform. A ratio-of-ratios variability metric indicating the variation of a ratio-of-ratios may be determined. A ratio-of-ratios may represent the ratio of absorbances of the red waveform and the infrared waveform. A pulse quality metric indicating the quality of one or more pulses may be determined. A saturation filter weight may be adjusted in accordance with the ratio-of-ratios variability metric and the pulse quality metric, where the saturation filter weight may represent a weight used for a filtering process operable to filter a saturation estimate of the blood perfused tissue.
Certain exemplary embodiments are described in the following detailed description and in reference to the drawings in which:
The exemplary embodiments described below are best understood by referring to
Pulse oximeter 120 analyzes the detector signals. According to the illustrated embodiment, pulse oximeter 120 includes general processing and interface components such as a microprocessor 122, a ROM 146, a RAM 126, a display 128, and control inputs 154 coupled to an internal bus 124. Pulse oximeter 120 also includes components that operate to control the light that passes through tissue 112. According to the illustrated embodiment, pulse oximeter 120 includes a time processing unit (TPU) 130 and light drive circuitry 132. TPU 130 provides timing control signals to light drive circuitry 132. Light drive circuitry 132 controls when light source 110 is illuminated, and may control multiplexed timing if multiple light sources 110 are used.
Signals from detector 114 are received through amplifier 133. Pulse oximeter 120 includes components that operate to process the received signal. According to the illustrated embodiment, pulse oximeter 120 includes a switching circuit 134, an amplifier 136, a low pass filter 138, an analog-to-digital converter 140, and a queued serial module (QSM) 142. Switching circuit 134 controls the sampling of the signals in response to instructions from TPU 130. If multiple light sources are used, the sampling times may depend upon which of the light sources 110 are illuminated.
The signals from the switch 134 are passed through amplifier 136, low pass filter 138, and analog-to-digital converter 140. Digital data from the signals is then stored in a queued serial module (QSM) 142. The digital data may be downloaded to RAM 126 as QSM 142 is filled. In one embodiment, there may be multiple parallel paths of separate amplifier, filter, and analog-to-digital converters for different light wavelengths.
Microprocessor 122 calculates oxygen saturation based on the values of the received signals. ROM 146 may store coefficients used in the calculations. Detector/decoder 144 selects the appropriate coefficients according to signals received from encoder 116. Control inputs 154 receive input data and instructions, and may comprise, for instance, a switch on the pulse oximeter, a keyboard, or a port providing instructions from a remote host computer. Display 128 provides feedback and results of the analysis.
One or more components of pulse oximeter 120 may include appropriate input devices, output devices, mass storage media, processors, memory, or other components for receiving, processing, storing, and communicating information according to the operation of pulse oximeter 120. As an example, one or more components of pulse oximeter 120 may include logic, an interface, memory, or any suitable combination of the preceding. By way of example, “logic” may refer to hardware, software, firmware, or any suitable combination of the preceding. Certain logic may manage the operation of a device, and may comprise, for example, a processor. “Processor” may refer to any suitable device operable to execute instructions and manipulate data to perform operations.
“Interface” may refer to logic of a device operable to receive input for the device, send output from the device, perform suitable processing of the input or output or both, or any combination of the preceding, and may comprise one or more ports, conversion software, or both. “Memory” may refer to logic operable to store and facilitate retrieval of information, and may comprise Random Access Memory (RAM), Read Only Memory (ROM), a magnetic drive, a disk drive, a Compact Disk (CD) drive, a Digital Video Disk (DVD) drive, removable media storage, any other suitable data storage medium, or a combination of any of the preceding.
Modifications, additions, or omissions may be made to pulse oximeter 120 without departing from the scope of the invention. The components of pulse oximeter 120 may be integrated or separated according to particular needs. Moreover, the operations of pulse oximeter 120 may be performed by more, fewer, or other modules. Additionally, operations of pulse oximeter 120 may be performed using any suitable logic comprising software, hardware, other logic, or any suitable combination of the preceding.
The brief description of an exemplary pulse oximeter set forth above, serves as a basis for describing the exemplary methods for adjusting oximetry parameters. Any suitable oximeter, however, may be used.
Block 202 represents the operations of a Signal Conditioning subsystem. Block 202 may receive digitized Red and IR signals, and may output pre-processed Red and IR signals. The Signal Conditioning subsystem conditions signals to emphasize higher frequencies associated with the human plethysmograph and to attenuate lower frequencies associated with interference from noise sources. The derivative-filtered plethysmographs characteristically have a negative skewness. The signals may be conditioned by taking a first derivative to reduce or eliminate a baseline shift, and then low pass filtering with coefficients based on hardware characteristics. The Signal Conditioning subsystem may also divide the lowpass filtered values by the average of the respective IR or Red signals.
Block 204 represents the operations of a Pulse Identification and Qualification subsystem. Block 204 may receive pre-processed Red and IR signals, average pulse period, lowpass waveforms from the Signal Conditioning subsystem, and/or ensemble averaging filter weights. Block 204 may output pulse quality of individual pulses, degree of arrhythmia, pulse amplitude variations, and/or qualified pulse periods and age.
The Pulse Identification and Qualification subsystem identifies pulses and qualifies the pulses as likely arterial pulses. Pulses may be identified and qualified by applying a pre-trained neural net to the IR signals and/or Red signals. Signal metrics describing pulses may be compared with pulse qualification criteria in order to qualify the pulses. A signal metric represents a feature of a signal, and it may be used to classify a signal. For example, a signal metric may be used to qualify a signal. Example signal metrics may describe features of individual pulses, such as amplitude, period, shape, ratio-of-ratios, and/or rise time. Example signal metrics may describe features of a sequence of pulses such as the variability of the features of individual pulses (e.g., period variability).
A pulse qualification criterion may include parameters with which the signal quality metrics may be compared. The parameters may be adjusted in response to variables such as filter parameters or signal metrics. A filter parameter may refer to a parameter of a filter that may be modified to change the filtering. An example embodiment of a method that adjusts parameters in response to filter parameters is described with reference to
Block 206 represents operations that compute signal quality metrics. Block 206 may receive: raw digitized Red and IR signals; degree of arrhythmia, individual pulse quality, pulse amplitude variation; pre-processed Red and IR signals; and/or average pulse period. Block 206 may output lowpass and ensemble averaging filter weights, normalized pre-processed waveforms, and/or percent modulation. The signal quality metrics may be used to set parameters for other processes.
Block 208 represents operations that compute the average pulse period from the received pulses. Block 208 may receive qualified pulse periods and age, and it may output the average pulse period.
Block 210 represents the operations of a Lowpass Filter and Ensemble Averaging subsystem. Block 210 may receive normalized pre-processed Red and IR signals, average pulse period, and/or low pass filter weights and ensemble averaging filter weights. Block 210 may output filtered Red and IR signals and/or age. The Lowpass Filter and Ensemble Averaging subsystem filters and ensemble averages normalized and preprocessed signals processed by block 206. Ensemble averaging may involve attenuating frequencies that are not of interest. For example, ensemble averaging may involve attenuating frequencies that are not at the estimated pulse rate or harmonic. The Lowpass Filter and Ensemble Averaging subsystem may also track the age of the signal and/or filtering results.
Block 212 represents operations that estimate the ratio-of-ratios variance for the filtered waveforms. Block 212 may receive filtered Red and IR signals, age, calibration coefficients, and/or response mode. Block 212 may output a ratio-of-ratios variance. A ratio-of-ratios is the ratio of the absorbances of the red and infrared signals. A ratio-of-ratios variability metric indicates the variation of a ratio-of-ratios. According to one embodiment, a ratio-of-ratios variance may be adjusted according to a signal metric and a ratio-of-ratios variability metric. An example method is described with reference to
Block 218 represents the operations of a Low Pass Filter and Ensemble Averaging subsystem. Block 218 may operate in a substantially similar manner as block 210. Block 220 represents the operations of a Filtered Pulse Identification and Qualification subsystem. Block 220 may operate in a substantially similar manner as block 204. The Filtered Pulse Identification and Qualification subsystem calculates and qualifies the pulse periods from the filtered waveforms. The results from the subsystem may be used by block 222 if a pulse period is disqualified by block 204.
Block 222 represents the operations of an Average Pulse Periods and Calculate Pulse Rate subsystem. Block 222 may receive qualified pulse periods and age, and may output an average pulse period and/or a pulse rate.
Block 224 represents the operations that detect venous pulsation. Block 224 may receive the pre-processed Red and IR signals and age from block 202, and it may output pulse rate and an indication of venous pulsation. Block 226 represents the operations that detect sensor-off and loss of pulse amplitude.
Modifications, additions, or omissions may be made to system 200 without departing from the scope of the invention. The components of system 200 may be integrated or separated according to particular needs. Moreover, the operations of system 200 may be performed by more, fewer, or other modules. Additionally, operations of system 200 may be performed using any suitable logic comprising software, hardware, firmware, or any suitable combination of the preceding.
The method starts at step 300, where an ensemble averaging weight is established. An ensemble averaging weight may refer to a weight value that is used to calculate a weighted average of new samples and previously ensemble averaged samples from a previous pulse period. Any suitable ensemble averaging weight may be used. According to one embodiment, Ensemble_Averaging_Weight used by the Ensemble Averaging subsystem may be used. If the pulses are not ensemble averaged, the weight may be set to default value, for example, Ensemble_Averaging_Weight=1.0.
Pulses are received at step 302. Variables describing the pulses are updated at step 304. The variables may be updated in any suitable manner. For example, the variables may be updated for each sample, prior to the every-potential-pulse, according to the following operations:
1. Baseline represents the average of input samples, and may be updated according to the following equation:
Baselinet=Baselinet-1+c1Δt*(Curr_Sample−Baselinet-1)
where t represents a sample index, c1 represents a constant, Δt represents the sample interval given in seconds, Curr_Sample represents the current sample. Constant c1 may be any suitable value, such as c1=0.01 for a one-second response time, given sampling interval Δt=10 milliseconds.
2. Mean_Square represents the mean-square of input samples, and may be updated according to the following equation:
Mean_Squaret=Mean_Squaret-1+k*((Curr_Sample−Baselinet)2−Mean_Squaret-1)
where k represents a constant selected to yield a Mean_Square with a particular response time.
As an example, k may be determined according to the following equation to yield a Mean_Square with a response time of one second or one pulse, whichever is shorter:
where Avg_Period represents the average pulse period, and Potential_Pulse_Period represents the pulse period of a potential pulse. A short response time may allow Mean_Square to decline quickly when a motion artifact ends, which may minimize the likelihood of ignoring real pulses. During the first potential pulse, the above equation for k may give equal weight to each sample for up to one second, which may prevent initial underestimates of Mean_Squaret as it diverges from zero.
3. Gated_RMS represents the square root of Mean_Square, and may be updated according to the following equation:
At step 308, the pulse period may exceed a threshold that indicates that a pulse has just been missed or will occur shortly. The threshold may have any suitable value, such as two seconds. If the pulse period exceeds the threshold at step 308, the method proceeds to step 312, where mean square parameters are reduced. According to one embodiment, the mean square parameters may include Mean_Square and Gated_RMS, and may be reduced according to the following equations:
Gated_RMSt=(1−Δt)Gated_RMSt
Mean_Squaret=(1−Δt)Mean_Squaret
The method then proceeds to step 316. If the pulse period does not exceed the threshold at step 308, the method proceeds directly to step 316.
The minimum phase of a potential pulse is identified at step 316. The minimum phase may be identified when the value of Curr_Samplet has a specified relationship to previous or subsequent Curr_Sample values. For example, a minimum phase may be identified when Curr_Samplet exceeds samples that occurred in the previous 100 to 150 milliseconds. Steps 320 through 332 describe updating the noise gate parameter of a noise gate criterion. A noise gate criterion may refer to a pulse qualification criterion that qualifies a pulse according to a noise gate. A noise gate parameter may refer to a parameter that controls the threshold that identifies whether a waveform includes noise. For example, the noise gate parameter may be adjusted to reduce the noise gate level if the ensemble averaging weight indicates an increased degree of ensemble averaging.
According to one embodiment, the noise gate parameter Noise_Gate may be defined according to the following equation:
Noise_Gate=max(n*Gated_RMS,(c2+c3*E—A—W)*m*Noise_Floor/IR_DC) 4.
where n represents the noise gate multiplier, E_A_W represents Ensemble_Averaging_Weight, m represents a noise floor multiplier, IR_DC represents the current infrared (IR) A/D values, Noise_Floor represents a minimum noise level associated with the oximetry hardware, which determines a minimum value for the noise gate, and c2 and c3 represent any suitable constants, for example, c2=c3=0.5.
The noise floor multiplier of the noise gate parameter is set at step 320. The noise floor multiplier may be set in any suitable manner. As an example, noise floor multiplier m may initially be set to an unit value, for example, m=1.0. Subsequent values of noise floor multiplier m may be set according to the following equation:
where Pulse_Period represents the pulse period of the current pulse, l represents a parameter, and c4, c5, and c6 represent constants having any suitable values, for example, c4=0.5, c5=0.25, and c6=0.75.
Parameter l may be given by the following equation:
l=w((1−c5)*Avg_Period−Potential_Pulse_Period)+(1−w)*(c5*Avg_Period−min(|c4*Avg_Period−Potential_Pulse_Period|,c5*Avg_Period))
where:
w=(Frequency_Ratio−c7)/c8
w=max(0, min(1, w))
and where Frequency_Ratio represents the frequency content of plethysmograph relative to the pulse rate, and c7 and c8 represent constants having any suitable values, for example, c7=1.25 and c8=0.60.
According to one embodiment, if Avg_Period is zero, then noise floor multiplier m may be calculated using the previous Potential_Pulse_Period in place of Avg_Period, provided that the previous Potential_Pulse_Period is acceptable for modifying noise floor multiplier m. For example, the duration of the previous Potential_Pulse_Period may be required to satisfy a threshold indicating that the pulse periods are less likely to reflect noise. According to the example, if Avg_Period is zero, then Mean_Square is established according to the following equation:
Mean_Squaret=max(Mean_Squaret,c9*Potential_Pulse_Amp2)
where Potential_Pulse_Amp represents the amplitude of a pulse, and c9 represents a constant having any suitable value, for example, c9=0.12. Potential_Pulse_Amp may be given as Potential_Pulse_Amp=Potential_Pulse_Max−Potential_Pulse_Min, where Potential_Pulse_Max represents the maximum of the previous several samples at the end of the maximum phase preceding the minimum phase, and Potential_Pulse_Min represents the minimum of the previous several samples at the end of the minimum phase.
The noise gate multiplier n is established at step 328. Noise gate multiplier n may be established in any suitable manner. According to one embodiment, noise gate multiplier n may be established according to the following equation:
n=c10*m*bound(c11+Potential_Pulse_Skew,c12,c13), if Pulse_Period>Potential_Pulse_Period;
c10*m*bound(c14+Potential_Pulse_Skew,c12,c13), if Skew_Derivative Input Weight>c15; and
c10*m otherwise.
where Skew_Derivative_Input_Weight represents a weight used to combine the Curr_Sample waveform with its derivative to obtain a waveform having a more negative skew, Potential_Pulse_Skew represents the skewness of the samples in this combined waveform over the duration of the Potential_Pulse_Period, and c10, c11, c12, c13, c14, and c15 represent constants having any suitable values, for example, c10=0.85, c11=1.5, c12=0.4, c13=1.0, c14=2.0, and c15=8.0. The notation bound(a, b, c) is used to denote min(max(a, b), c).
The noise gate parameter is adjusted according to the noise floor multiplier and the noise gate multiplier at step 332. After updating the noise gate parameter, the method ends.
Modifications, additions, or omissions may be made to the method without departing from the scope of the invention. The method may include more, fewer, or other steps. Additionally, steps may be performed in any suitable order without departing from the scope of the invention.
The method starts at step 350, where an ensemble averaging weight is established. Pulses are received at step 354. Steps 358 through 370 describe evaluating the received pulses in accordance with pulse qualification criteria. The pulse qualification criteria may be applied in any suitable order. According to one embodiment, the pulse qualification criteria are applied in the order of the steps, and if a pulse fails one criterion, subsequent criteria are not applied.
Pulse qualification criteria comprising one or more sample criteria are applied at step 358. Any suitable sample criteria SCi may be used, for example:
SC1: Pulse_Min<Pulse_Avg<Pulse_Max
SC2: Avg_Min<Pulse_Avg<Pulse_Max
where Pulse_Min represents the minimum of the previous several samples at the end of the minimum phase, Pulse_Avg represents the average of the samples of a period, and Pulse_Max represents the maximum of the previous several samples at the end of the maximum phase. Pulse_Avg may be expressed as
where Pulse_Sum represents the sum of the samples of a pulse, and Pulse_Period represents the pulse period of the current pulse. Pulse_Sum may be updated at every sample and reset to zero after calculating Pulse_Avg.
Pulse qualification criteria comprising one or more output criteria are applied at step 362. Any suitable output criteria OCi may be used, for example:
OC1: Pulse_Qual_NN_Outputt>NN_Thresh
OC2: Pulse_Qual_NN_Output_Integralt≧c8, if no pulses have been qualified yet
where Pulse_Qual_NN_Output represents the output of the pre-trained pulse qualification neural net, NN_Thresh represents a neural net threshold that may be used to qualify a pulse, Pulse_Qual_Output_Integral represents the integral of Pulse_Qual_NN_Outputt, and c8 represents a constant. The output of the neural net may range from 0 to 1.
In one example, NN_Thresh may be given by the following equation:
NN_Thresh=c1+c2*(Period_Ratio−c3); and
c4, if there are no qualified pulses, for example, if Avg_Period=0 where Period_Ratio represents the ratio of the current period and the average period, and c1, c2, c3, and c4 represent constants having any suitable values, for example, c1=0.4, c2=0.25, c3=1.0, and c4=0.5. Period_Ratio may be given as
In one example, Pulse_Qual_Output_Integral may be given by the following equation:
Pulse_Qual_Output_Integralt=min(Pulse_Qual_Output_Integralt-1+Pulse_Qual—NN_Outputt−c5,c6); and
0, if Pulse_Qual_NN_Outputt≦c7
where c5, c6, and c7 represent constants having any suitable values, for example, c5=0.5, c6=100, and c7=0.5, and the subscript t−1 denotes the previous pulse. In one example, constant c8 may be set such that the initial estimate for Avg_Period is determined from either a very good pulse (for example, >95% probability of being an arterial pulse) or multiple acceptable successive pulses.
Pulse qualification criteria comprising one or more initialization criteria are applied at step 366. Any suitable initialization criterion ICi may be used, for example:
IC1: |Pulse_Period−Last_Pulse_Period|≦c9*Initial_Pulse_Count*Pulse_Period
where Pulse_Period represents the pulse period of the current pulse, Last_Pulse_Period represents the Pulse_Period for the previous non-ignored pulse, Initial_Pulse_Count represents the number of non-ignored pulses found since re-initialization, and c9 represents a constant having any suitable value, for example, c9=0.1. Initial_Pulse_Count may be incremented after a pulse is evaluated according to the pulse qualification criteria. Criterion IC1 requires that consecutive pulse periods agree closely, but gradually relaxes its threshold with each successive pulse, which may reduce or prevent inaccurate initial pulse estimates. According to one embodiment, the first pulse may be rejected by requiring Initial_Pulse_Count>1.
IC2: Path_Ratio<c10, if no pulses have been qualified yet where Path_Ratio represents the ratio of the path length and the pulse amplitude, and c10 represents a constant having any suitable value. Path_Ratio may be given by
where Path_Length represents the sum of sample-to-sample differences for the samples of a pulse. Constant c10 may be selected to reduce the likelihood of initializing the pulse-rate estimate to values that are too low under noisy conditions. For example, c10=4.0.
A pulse period criterion is adjusted according to the ensemble averaging weight at step 368. A pulse period criterion may refer to a pulse qualification criterion that qualifies a pulse according to a pulse period. Any suitable pulse period criterion PPCi may be used. According to one embodiment, a pulse with a shorter-than-average period may be disqualified if the ensemble averaging weight indicates that qualified pulses are less likely to have shorter-than-average period. As an example, pulse period criterion PPC1 may be used and adjusted according to the following:
PPC1: (Avg_Period−Pulse_Period)≦2*Avg_Period*min(Ensemble_Averaging_Weightt,Ensemble_Averaging_Weightt-1), if Avg_Period>0
where Avg_Period represents an estimate of the average pulse period, Pulse_Period represents the pulse period of the current pulse. Pulse qualification criteria comprising the pulse period criterion are applied at step 370.
The pulses are qualified or disqualified according to the evaluations at step 374. The pulses may be qualified or disqualified according to the evaluations in any suitable manner. According to one embodiment, a pulse that satisfies the above pulse qualification criteria is designated as qualified, and a pulse that fails to satisfy any of the above pulse qualification criteria is designated as non-qualified. After qualifying the pulses, the method ends.
Modifications, additions, or omissions may be made to the method without departing from the scope of the invention. The method may include more, fewer, or other steps. Additionally, steps may be performed in any suitable order without departing from the scope of the invention.
The method starts at step 400, where pulses are received. Steps 412 through 424 describe evaluating the received pulses in accordance with pulse qualification criteria. The pulse qualification criteria may be applied in any suitable order. According to one embodiment, the pulse qualification criteria are applied in the order of the steps, and if a pulse fails one criterion, subsequent criteria are not applied.
Pulse qualification criteria comprising one or more sample criteria are applied at step 412. Any suitable sample criterion SCi may be used, for example, SC1 and SC2. Pulse qualification criteria comprising one or more output criteria are applied at step 416. Any suitable output criterion OCi may be used, for example, OC1 and OC2. Pulse qualification criteria comprising one or more initialization criteria are applied at step 420. Any suitable initialization criterion ICi may be used, for example, IC1 and IC2.
Pulse qualification criteria comprising one or more pulse amplitude criteria are applied at step 424. Any suitable pulse amplitude criterion PSCi may be used. According to one embodiment, a pulse amplitude criterion may disqualify a pulse with an amplitude that is much smaller than (e.g., more than three times smaller) an average amplitude and a period that is much greater that an average period. The pulse amplitude criterion may reduce the likelihood of qualifying instrument noise when the pulse amplitude suddenly disappears.
An example pulse amplitude criterion PSC1 may be expressed as:
PSC1: Sum_Amp_Diff≧−c1; OR
Period_Ratio≦c2; OR
Period_Ratio≦(c1+Sum_Amp_Diff/c3), provided that Pulse_Period>Avg_Period, pulses have been previously qualified, and a significant number of pulses have been previously evaluated;
where Sum_Amp_Diff represents the sum of amplitudes, Period_Ratio represents the ratio of the current period and the average period, and c1, c2, and c3 represent constants having any suitable values, for example, c1=2.3, c2=1.3, and c3=6.
According to one embodiment, Sum_Amp_Diff may be defined according to the following equation:
Sum_Amp_Diff=Amp_Diff2t+Amp_Diff2t-1
where:
Amp_Diff2=ln(Pulse_Ampt)/Avg_Ln_Pulse_Ampt
where Avg_Ln_Pulse_Amp represents a variable-weight filtered version of ln(Pulse_Amp), with a low initial value, and Pulse_Amp represents the amplitude of a pulse in the Curr_Sample waveform.
The pulses are qualified or disqualified according to the evaluations at step 428. The pulses may be qualified or disqualified according to the evaluations in any suitable manner. According to one embodiment, a pulse that satisfies the above pulse qualification criteria is designated as qualified, and a pulse that fails to satisfy any of the above pulse qualification criteria is designated as non-qualified. After the pulses are qualified, the method ends.
Modifications, additions, or omissions may be made to the method without departing from the scope of the invention. The method may include more, fewer, or other steps. Additionally, steps may be performed in any suitable order without departing from the scope of the invention.
The method starts at step 450, where a saturation weight is established. The method may be applied to any suitable filter that filters saturation estimates. According to one embodiment, the filter may comprise a Kalman filter. For a Kalman filter, the saturation s may be updated for each sample from an estimation error P and Kalman gain K according to the following equation:
st=st-1+Kt(vt−utst-1)
where:
Kt=PtutRt-1−1
Pt=((Pt-1++Q)−1+ut2Rt-1−1)−1
and u and v represent inputs, Q represents the variance of saturation s, and R represents the variance of noise. Inputs u and v may represent inputs calculated from the IR and Red plethysmographs utilizing coefficients defined by encoder 116. The ratio v/u may converge to the oxygen saturation estimate in the absence of artifacts.
According to one embodiment, noise variance R may be given by the following equation:
where N represents a suitable number of samples, such as one pulse period or one second.
Saturation weight wt may be given by the following equation:
Saturation weight wt depends on Kalman gain K, which in turn depends on noise variance R. Accordingly, a change in the noise variance may yield a change in the saturation weight. Noise variance R may be dependent, at least in part, on the filtered saturation output of the Kalman filter.
A pulse quality metric is determined at step 462. The pulse quality metric Pulse_Qual may reflect the quality of the most recently evaluated pulse. According to one embodiment, the pulse quality metrics range from 0 to 1.
The saturation filtering may be in one of multiple response modes based on the response time of the filtering at step 466. As an example, the response modes may include a fast response mode, normal response mode, or slow response mode, where the fast response mode has a time of less than 4 seconds, normal response mode has a time of greater than 4 and less than 6, or slow response mode has a time of more than 6 seconds. The response mode may be user-selected or automatically selected.
If the saturation filtering is in a fast or a normal response mode at step 466, the method proceeds to step 468, where the noise variance is adjusted. According to one embodiment, the noise variance in the fast response mode may be adjusted according to the following equation:
where b is a constant selected in accordance with the response mode. According to the embodiment, the noise variance in the normal response mode may be adjusted according to the following equation:
where ρ represents a second noise variance estimate. The second noise-variance estimate is designed to be less sensitive to rapid physiological saturation changes.
If the saturation filtering is in the slow response mode at step 466, the method proceeds to step 470, where the noise variance is adjusted in accordance with pulse quality. According to one embodiment, the noise variance in the fast response mode may be adjusted according to the following equation:
where:
temp=max(min(c1*Pulse_Qual,c2),c3)
and c1, c2, and c3 represent any suitable constants, for example, c1=2, c2=1, and c3=0.01. According to the above equation, a change in pulse quality may yield a change in noise variance. As discussed above, a change in the noise variance may yield a change in the saturation weight. Accordingly, saturation weight may be adjusted in accordance with the pulse quality.
The saturation weight is adjusted in accordance with the noise variance c at step 474. After adjusting the saturation weight, the method ends.
Modifications, additions, or omissions may be made to the method without departing from the scope of the invention. The method may include more, fewer, or other steps. Additionally, steps may be performed in any suitable order without departing from the scope of the invention.
The method starts at step 550, where a saturation weight is established. A ratio-of-ratios variability metric is determined at step 558. Any suitable ratio-of-ratios variability metric may be used. According to one embodiment, the ratio-of-ratios
where Red represents red input values, IR represents infrared input values, and N represents the time period divided by the sample time.
According to the embodiment, a ratio-of-ratios variability metric RoR_Variance_Per may be given by the following equation:
A pulse quality metric is determined at step 562. The pulse quality metric Pulse_Qual may reflect the quality of the most recently evaluated pulse.
The saturation filtering may be in one of multiple response modes based on the response time of the filtering at step 566. As an example, the response modes may include a fast response mode, normal response mode, or slow response mode. If the saturation filtering is in a fast or a normal response mode at step 566, the method proceeds to step 568, where the saturation weight is calculated from the ratio-of-ratios variance metric. According to one embodiment, the saturation weight may be adjusted by first calculating an intermediate value Sat_Noise_Per according to the following equation:
and where dSat_dRoR represents the derivative of the saturation with respect to the ratio-of-ratios at a ratio-of-ratios value equal to RoR. This derivative may be calculated using the calibration coefficients selected by the oximeter based on the signals from the encoder 116. According to one embodiment, if Sat_Noise_Per is greater then a maximum threshold, Sat_Noise_Per may be set to a maximum value.
The saturation weight may be adjusted using the intermediate value Sat_Noise_Per. If Sat_Noise_Per is greater than a threshold c4, such as c4=0.09, then saturation weight may be adjusted according to the following equation:
wt=wt*c4/Sat_Noise_Per
According to the embodiment, before the first second of filtering has elapsed, the minimum filter weight wm may be set at wm=1.0. Otherwise, minimum filter weight wm may be set at wm=0.5 divided by the number of seconds.
If the saturation filtering is in the slow response mode at step 566, the method proceeds to step 570, where the saturation weight is calculated from the ratio-of-ratios variance metric and then adjusted in accordance with the pulse quality metric. The saturation weight wt may be calculated in a manner similar to that of step 568, except that wt may be adjusted according to the pulse quality metric. For example, the equation for RoR_Noise_Per may be modified to:
RoR_Noise_Per RoR_Variance_Per*RoR/temp
where temp=max(min(c1*Pulse_Qual, c2), c3), and c1, c2, and C3 represent any suitable constants, for example, c1=2, c2=1, and c3=0.1.
Saturation is calculated at step 574. Saturation may be calculated according to the following equation:
Saturationt=wt*Raw_Saturationt+(1−wt)*Saturationt-1
where Raw_Saturationt is calculated from
Modifications, additions, or omissions may be made to the method without departing from the scope of the invention. The method may include more, fewer, or other steps. Additionally, steps may be performed in any suitable order without departing from the scope of the invention.
As an example, the above-described saturation weight might be used to filter the ratio-of-ratios prior to saturation calculation, rather than calculating saturation first and then filtering the ratio-of-ratios. According to one embodiment, a ratio-of-ratios weight may be established in accordance with the ratio-of-ratios variability metric and the pulse quality metric. The ratio-of-ratios weight may be established in a manner similar to that of step 570. The ratio-of-ratios may be filtered in accordance with the ratio-of-ratios weight, and the saturation calculated according to the filtered ratio-of ratios.
Certain embodiments of the invention may provide one or more technical advantages. A technical advantage of one embodiment may be that a noise gate parameter of a noise gate criterion may be adjusted in accordance with an ensemble averaging weight. Ensemble averaging filtering may reduce the noise level, which may allow for a lowered noise gate. Accordingly, the ensemble averaging weight, which indicates the degree of ensemble averaging, may be used to adjust the noise gate parameter.
Another technical advantage of one embodiment may be that a pulse period criterion may be adjusted in accordance with an ensemble averaging weight. According to the embodiment, an ensemble averaging weight may indicate that a qualified pulse is not likely to have a pulse period that is substantially shorter than an average pulse period. Accordingly, a pulse period criterion may be adjusted to disqualify a pulse having a pulse period that is substantially shorter than an average pulse period.
Another technical advantage of one embodiment may be that a pulse amplitude criterion may be adjusted in accordance with an average pulse amplitude and an average pulse period. According to one embodiment, a pulse amplitude criterion may disqualify a pulse with an amplitude that is much smaller than (e.g., more than three times smaller) an average amplitude and a period that is much greater that an average period. The pulse amplitude criterion may reduce the likelihood of qualifying instrument noise when the pulse amplitude suddenly disappears.
Another technical advantage of one embodiment may be that a saturation weight may be adjusted in accordance with a ratio-of-ratios variability metric and a pulse quality metric. The metrics may be used to adjust a saturation weight, which may enhance accuracy under conditions such as motion artifact, low perfusion, or electronic interference.
While the invention may be susceptible to various modifications and alternative forms, specific embodiments have been shown by way of example in the drawings and have been described in detail herein. However, it should be understood that the invention is not intended to be limited to the particular forms disclosed. Rather, the invention is to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the invention as defined by the following appended claims.
This application is a continuation of U.S. application Ser. No. 11/261,012 filed on Oct. 28, 2005, now abandoned, which is incorporated herein by reference in its entirety.
Number | Name | Date | Kind |
---|---|---|---|
3638640 | Shaw | Feb 1972 | A |
4714341 | Hamaguri et al. | Dec 1987 | A |
4805623 | Jöbsis | Feb 1989 | A |
4807631 | Hersh et al. | Feb 1989 | A |
4911167 | Corenman et al. | Mar 1990 | A |
4913150 | Cheung et al. | Apr 1990 | A |
4936679 | Mersch | Jun 1990 | A |
4938218 | Goodman et al. | Jul 1990 | A |
4971062 | Hasebe et al. | Nov 1990 | A |
4972331 | Chance | Nov 1990 | A |
4974591 | Awazu et al. | Dec 1990 | A |
5028787 | Rosenthal et al. | Jul 1991 | A |
5065749 | Hasebe et al. | Nov 1991 | A |
5084327 | Stengel | Jan 1992 | A |
5119815 | Chance | Jun 1992 | A |
5122974 | Chance | Jun 1992 | A |
5167230 | Chance | Dec 1992 | A |
5190038 | Polson et al. | Mar 1993 | A |
5246003 | DeLonzor | Sep 1993 | A |
5247931 | Norwood | Sep 1993 | A |
5263244 | Centa et al. | Nov 1993 | A |
5275159 | Griebel | Jan 1994 | A |
5279295 | Martens et al. | Jan 1994 | A |
5297548 | Pologe | Mar 1994 | A |
5355880 | Thomas et al. | Oct 1994 | A |
5372136 | Steuer et al. | Dec 1994 | A |
5385143 | Aoyagi | Jan 1995 | A |
5390670 | Centa et al. | Feb 1995 | A |
5413099 | Schmidt et al. | May 1995 | A |
5469845 | DeLonzor et al. | Nov 1995 | A |
5482036 | Diab et al. | Jan 1996 | A |
5483646 | Uchikoga | Jan 1996 | A |
5485847 | Baker, Jr. | Jan 1996 | A |
5553614 | Chance | Sep 1996 | A |
5564417 | Chance | Oct 1996 | A |
5575285 | Takanashi et al. | Nov 1996 | A |
5611337 | Bukta | Mar 1997 | A |
5630413 | Thomas et al. | May 1997 | A |
5645059 | Fein et al. | Jul 1997 | A |
5645060 | Yorkey | Jul 1997 | A |
5680857 | Pelikan et al. | Oct 1997 | A |
5692503 | Kuenstner | Dec 1997 | A |
5730124 | Yamauchi | Mar 1998 | A |
5743263 | Baker, Jr. | Apr 1998 | A |
5758644 | Diab et al. | Jun 1998 | A |
5779631 | Chance | Jul 1998 | A |
5782757 | Diab et al. | Jul 1998 | A |
5786592 | Hök | Jul 1998 | A |
5830136 | DeLonzor et al. | Nov 1998 | A |
5830139 | Abreu | Nov 1998 | A |
5831598 | Kauffert et al. | Nov 1998 | A |
5842981 | Larsen et al. | Dec 1998 | A |
5846190 | Woehrle | Dec 1998 | A |
5853364 | Baker, Jr. et al. | Dec 1998 | A |
5865736 | Baker, Jr. et al. | Feb 1999 | A |
5871442 | Madarasz et al. | Feb 1999 | A |
5873821 | Chance et al. | Feb 1999 | A |
5920263 | Huttenhoff et al. | Jul 1999 | A |
5995855 | Kiani et al. | Nov 1999 | A |
5995856 | Mannheimer et al. | Nov 1999 | A |
5995859 | Takahashi | Nov 1999 | A |
6011986 | Diab et al. | Jan 2000 | A |
6035223 | Baker, Jr. | Mar 2000 | A |
6064898 | Aldrich | May 2000 | A |
6081742 | Amano et al. | Jun 2000 | A |
6083172 | Baker, Jr. et al. | Jul 2000 | A |
6088607 | Diab et al. | Jul 2000 | A |
6120460 | Abreu | Sep 2000 | A |
6134460 | Chance | Oct 2000 | A |
6150951 | Olejniczak | Nov 2000 | A |
6154667 | Miura et al. | Nov 2000 | A |
6163715 | Larsen et al. | Dec 2000 | A |
6181958 | Steuer et al. | Jan 2001 | B1 |
6181959 | Schöllermann et al. | Jan 2001 | B1 |
6230035 | Aoyagi et al. | May 2001 | B1 |
6266546 | Steuer et al. | Jul 2001 | B1 |
6285895 | Ristolainen et al. | Sep 2001 | B1 |
6312393 | Abreu | Nov 2001 | B1 |
6339715 | Bahr et al. | Jan 2002 | B1 |
6353750 | Kimura et al. | Mar 2002 | B1 |
6393311 | Edgar, Jr. et al. | May 2002 | B1 |
6397091 | Diab et al. | May 2002 | B2 |
6411833 | Baker, Jr. et al. | Jun 2002 | B1 |
6415236 | Kobayashi et al. | Jul 2002 | B2 |
6419671 | Lemberg | Jul 2002 | B1 |
6438399 | Kurth | Aug 2002 | B1 |
6461305 | Schnall | Oct 2002 | B1 |
6463311 | Diab | Oct 2002 | B1 |
6466809 | Riley | Oct 2002 | B1 |
6487439 | Skladnev et al. | Nov 2002 | B1 |
6501974 | Huiku | Dec 2002 | B2 |
6501975 | Diab et al. | Dec 2002 | B2 |
6519486 | Edgar, Jr. et al. | Feb 2003 | B1 |
6526301 | Larsen et al. | Feb 2003 | B2 |
6544193 | Abreu | Apr 2003 | B2 |
6546267 | Sugiura et al. | Apr 2003 | B1 |
6549795 | Chance | Apr 2003 | B1 |
6564077 | Mortara | May 2003 | B2 |
6580086 | Schulz et al. | Jun 2003 | B1 |
6591122 | Schmitt | Jul 2003 | B2 |
6594513 | Jobsis et al. | Jul 2003 | B1 |
6606509 | Schmitt | Aug 2003 | B2 |
6606511 | Ali et al. | Aug 2003 | B1 |
6615064 | Aldrich | Sep 2003 | B1 |
6618042 | Powell | Sep 2003 | B1 |
6622095 | Kobayashi et al. | Sep 2003 | B2 |
6647280 | Bahr et al. | Nov 2003 | B2 |
6654621 | Palatnik et al. | Nov 2003 | B2 |
6654624 | Diab et al. | Nov 2003 | B2 |
6658276 | Kianl et al. | Dec 2003 | B2 |
6658277 | Wasserman | Dec 2003 | B2 |
6662030 | Khalil et al. | Dec 2003 | B2 |
6668183 | Hicks et al. | Dec 2003 | B2 |
6671526 | Aoyagi et al. | Dec 2003 | B1 |
6671528 | Steuer et al. | Dec 2003 | B2 |
6675031 | Porges et al. | Jan 2004 | B1 |
6678543 | Diab et al. | Jan 2004 | B2 |
6684090 | Ali et al. | Jan 2004 | B2 |
6690958 | Walker et al. | Feb 2004 | B1 |
6697658 | Al-Ali | Feb 2004 | B2 |
6701170 | Stetson | Mar 2004 | B2 |
6708048 | Chance | Mar 2004 | B1 |
6711424 | Fine et al. | Mar 2004 | B1 |
6711425 | Reuss | Mar 2004 | B1 |
6714245 | Ono | Mar 2004 | B1 |
6721584 | Baker, Jr. et al. | Apr 2004 | B2 |
6725074 | Kastle | Apr 2004 | B1 |
6731274 | Powell | May 2004 | B2 |
6785568 | Chance | Aug 2004 | B2 |
6793654 | Lemberg | Sep 2004 | B2 |
6801797 | Mannheimer et al. | Oct 2004 | B2 |
6801798 | Geddes et al. | Oct 2004 | B2 |
6801799 | Mendelson | Oct 2004 | B2 |
6810277 | Edgar, Jr. et al. | Oct 2004 | B2 |
6829496 | Nagai et al. | Dec 2004 | B2 |
6836679 | Baker, Jr. et al. | Dec 2004 | B2 |
6850053 | Daalmans et al. | Feb 2005 | B2 |
6863652 | Huang et al. | Mar 2005 | B2 |
6873865 | Steuer et al. | Mar 2005 | B2 |
6889153 | Dietiker | May 2005 | B2 |
6898451 | Wuori | May 2005 | B2 |
6898452 | Al-Ali et al. | May 2005 | B2 |
6939307 | Dunlop | Sep 2005 | B1 |
6947780 | Scharf | Sep 2005 | B2 |
6949081 | Chance | Sep 2005 | B1 |
6961598 | Diab | Nov 2005 | B2 |
6983178 | Fine et al. | Jan 2006 | B2 |
6987994 | Mortz | Jan 2006 | B1 |
6993371 | Kiani et al. | Jan 2006 | B2 |
6996427 | Ali et al. | Feb 2006 | B2 |
7003339 | Diab et al. | Feb 2006 | B2 |
7006856 | Baker, Jr. et al. | Feb 2006 | B2 |
7016715 | Stetson | Mar 2006 | B2 |
7024235 | Melker et al. | Apr 2006 | B2 |
7027849 | Al-Ali | Apr 2006 | B2 |
7030749 | Al-Ali | Apr 2006 | B2 |
7035697 | Brown | Apr 2006 | B1 |
7039538 | Baker, Jr. | May 2006 | B2 |
7047056 | Hannula et al. | May 2006 | B2 |
7072702 | Edgar, Jr. et al. | Jul 2006 | B2 |
7079880 | Stetson | Jul 2006 | B2 |
7127278 | Melker et al. | Oct 2006 | B2 |
7130671 | Baker, Jr. et al. | Oct 2006 | B2 |
7162306 | Caby et al. | Jan 2007 | B2 |
7194293 | Baker, Jr. | Mar 2007 | B2 |
7209775 | Bae et al. | Apr 2007 | B2 |
7215984 | Diab et al. | May 2007 | B2 |
7215986 | Diab et al. | May 2007 | B2 |
7221971 | Diab et al. | May 2007 | B2 |
7236811 | Schmitt | Jun 2007 | B2 |
7254433 | Diab et al. | Aug 2007 | B2 |
7263395 | Chan et al. | Aug 2007 | B2 |
7272426 | Schmid | Sep 2007 | B2 |
7302284 | Baker, Jr. et al. | Nov 2007 | B2 |
7315753 | Baker, Jr. et al. | Jan 2008 | B2 |
7328053 | Diab et al. | Feb 2008 | B1 |
7336983 | Baker, Jr. et al. | Feb 2008 | B2 |
7343187 | Stetson | Mar 2008 | B2 |
7373193 | Al-Ali et al. | May 2008 | B2 |
7376453 | Diab et al. | May 2008 | B1 |
7383070 | Diab et al. | Jun 2008 | B2 |
7392075 | Baker, Jr. | Jun 2008 | B2 |
7454240 | Diab et al. | Nov 2008 | B2 |
7534212 | Baker, Jr. | May 2009 | B2 |
20010005773 | Larsen et al. | Jun 2001 | A1 |
20010020122 | Steuer et al. | Sep 2001 | A1 |
20010039376 | Steuer et al. | Nov 2001 | A1 |
20010044700 | Kobayashi et al. | Nov 2001 | A1 |
20020026106 | Khalil et al. | Feb 2002 | A1 |
20020035318 | Mannheimer et al. | Mar 2002 | A1 |
20020038079 | Steuer et al. | Mar 2002 | A1 |
20020042558 | Mendelson | Apr 2002 | A1 |
20020049389 | Abreu | Apr 2002 | A1 |
20020062071 | Diab et al. | May 2002 | A1 |
20020111748 | Kobayashi et al. | Aug 2002 | A1 |
20020133068 | Huiku | Sep 2002 | A1 |
20020156354 | Larson | Oct 2002 | A1 |
20020161287 | Schmitt | Oct 2002 | A1 |
20020161290 | Chance | Oct 2002 | A1 |
20020165439 | Schmitt | Nov 2002 | A1 |
20020198443 | Ting | Dec 2002 | A1 |
20030023140 | Chance | Jan 2003 | A1 |
20030055324 | Wasserman | Mar 2003 | A1 |
20030060693 | Monfre et al. | Mar 2003 | A1 |
20030139687 | Abreu | Jul 2003 | A1 |
20030144584 | Mendelson | Jul 2003 | A1 |
20030220548 | Schmitt | Nov 2003 | A1 |
20030220576 | Diab | Nov 2003 | A1 |
20040010188 | Wasserman | Jan 2004 | A1 |
20040054270 | Pewzner et al. | Mar 2004 | A1 |
20040087846 | Wasserman | May 2004 | A1 |
20040107065 | Al-Ali | Jun 2004 | A1 |
20040127779 | Steuer et al. | Jul 2004 | A1 |
20040171920 | Mannheimer et al. | Sep 2004 | A1 |
20050085735 | Baker, Jr. et al. | Apr 2005 | A1 |
20050113656 | Chance | May 2005 | A1 |
20050143634 | Baker, Jr. et al. | Jun 2005 | A1 |
20050197549 | Baker, Jr. | Sep 2005 | A1 |
20050197551 | Al-Ali et al. | Sep 2005 | A1 |
20050197552 | Baker, Jr. | Sep 2005 | A1 |
20050209517 | Diab et al. | Sep 2005 | A1 |
20060030766 | Stetson | Feb 2006 | A1 |
20060258923 | Al-Ali et al. | Nov 2006 | A1 |
20060258924 | Al-Ali et al. | Nov 2006 | A1 |
20060258925 | Al-Ali et al. | Nov 2006 | A1 |
20060258927 | Edgar, Jr. et al. | Nov 2006 | A1 |
20060270920 | Al-Ali et al. | Nov 2006 | A1 |
20060281983 | Al-Ali et al. | Dec 2006 | A1 |
20070225581 | Diab et al. | Sep 2007 | A1 |
20070225582 | Diab et al. | Sep 2007 | A1 |
20070249918 | Diab et al. | Oct 2007 | A1 |
20070291832 | Diab et al. | Dec 2007 | A1 |
20080004514 | Diab et al. | Jan 2008 | A1 |
20080033265 | Diab et al. | Feb 2008 | A1 |
20080033266 | Diab et al. | Feb 2008 | A1 |
20080045823 | Diab et al. | Feb 2008 | A1 |
Number | Date | Country |
---|---|---|
19640807 | Sep 1997 | DE |
1491135 | Dec 2004 | EP |
3170866 | Jul 1991 | JP |
3238813 | Oct 1991 | JP |
4191642 | Jul 1992 | JP |
4332536 | Nov 1992 | JP |
7124138 | May 1995 | JP |
7136150 | May 1995 | JP |
10216115 | Aug 1998 | JP |
2003194714 | Jul 2003 | JP |
2003210438 | Jul 2003 | JP |
2003275192 | Sep 2003 | JP |
2003339678 | Dec 2003 | JP |
2004008572 | Jan 2004 | JP |
2004113353 | Apr 2004 | JP |
2004135854 | May 2004 | JP |
2004194908 | Jul 2004 | JP |
2004202190 | Jul 2004 | JP |
2004248819 | Sep 2004 | JP |
2004290545 | Oct 2004 | JP |
WO9101678 | Feb 1991 | WO |
WO9221281 | Dec 1992 | WO |
WO9309711 | May 1993 | WO |
WO9316629 | Sep 1993 | WO |
WO9403102 | Feb 1994 | WO |
WO9512349 | May 1995 | WO |
WO9749330 | Dec 1997 | WO |
WO9842249 | Oct 1998 | WO |
WO9842251 | Oct 1998 | WO |
WO9843071 | Oct 1998 | WO |
WO9932030 | Jul 1999 | WO |
WO0021438 | Apr 2000 | WO |
WO2005009221 | Feb 2005 | WO |
Number | Date | Country | |
---|---|---|---|
20090253971 A1 | Oct 2009 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 11261012 | Oct 2005 | US |
Child | 12483810 | US |