ES2662198T3 - Procesamiento de señal para oximetría de pulso - Google Patents

Procesamiento de señal para oximetría de pulso Download PDF

Info

Publication number
ES2662198T3
ES2662198T3 ES06815720.5T ES06815720T ES2662198T3 ES 2662198 T3 ES2662198 T3 ES 2662198T3 ES 06815720 T ES06815720 T ES 06815720T ES 2662198 T3 ES2662198 T3 ES 2662198T3
Authority
ES
Spain
Prior art keywords
functions
input signal
pffj
function
subset
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.)
Active
Application number
ES06815720.5T
Other languages
English (en)
Inventor
Yoram Wasserman
Guy Russell Lowery
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.)
Conmed Corp
Original Assignee
Conmed Corp
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 Conmed Corp filed Critical Conmed Corp
Application granted granted Critical
Publication of ES2662198T3 publication Critical patent/ES2662198T3/es
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/024Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02416Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, 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/024Detecting, measuring or recording pulse rate or heart rate
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
    • A61B5/14551Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters for measuring blood gases
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
    • A61B5/1455Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters

Abstract

Un método de determinación del ritmo del pulso de un sujeto vivo a partir de una señal de entrada de oximetría de pulso, comprendiendo el método las etapas de: comparar la señal de entrada para cada uno de una pluralidad de diferentes desplazamientos de tiempo para cada elemento de un conjunto de funciones, (i) en donde el conjunto incluye funciones con una pluralidad de diferentes frecuencias entre 30 a 240 ciclos por minuto; (ii) en donde el conjunto contiene un primer subconjunto que incluye funciones en las que el tiempo de elevación es más corto que el tiempo de caída y un segundo subconjunto que incluye funciones en las que (a) el tiempo de elevación es aproximadamente el mismo que el tiempo de caída O (b) el tiempo de elevación es más largo que el tiempo de caída, (iii) en donde la forma de todas las funciones en el primer subconjunto es similar, y (iv) en donde la forma de todas las funciones en el segundo subconjunto es similar; seleccionar, basándose en los resultados obtenidos en la etapa de comparación, una función a partir del conjunto que mejor coincide con la señal de entrada; y usar la frecuencia de la función seleccionada como una estimación del ritmo del pulso del sujeto.

Description

5
10
15
20
25
30
35
40
45
50
55
60
65
DESCRIPCION
Procesamiento de señal para oximetría de pulso Referencia cruzada a solicitudes relacionadas
La presente solicitud reivindica el beneficio de la Solicitud de Patente Provisional de Estados Unidos N.° 60/722.257, presentada el 30 de septiembre de 2005.
Antecedentes
La presente solicitud se refiere a la determinación del ritmo del pulso de una señal biológica pulsátil en presencia de ruido. Las técnicas descritas en el presente documento son particularmente útiles para el procesamiento de señales desde sensores de oximetría de pulso.
La oximetría de pulso es un procedimiento diagnóstico no invasivo para la medición del nivel de saturación de oxígeno en la sangre arterial de un paciente. La oximetría de pulso se basa en el principio de pasar la energía luminosa de al menos dos longitudes de onda en un medio fisiológico absorbente de la luz, adquirir la luz emitida reflejada (o transmitida), en respuesta a la absorción de luz, y calcular el nivel de saturación de oxígeno a partir de las señales adquiridas. Los oxímetros de pulso típicos tienen dos componentes principales: un sensor fijado a la piel del paciente para la adquisición de señales, y una unidad de procesamiento para procesar las señales adquiridas con el fin de determinar la saturación de oxígeno de la sangre arterial y el ritmo del pulso. Desafortunadamente, los sistemas de oximetría de pulso convencionales pueden ser sensibles al ruido, lo que puede dar como resultado lecturas inestables, mediciones imprecisas del ritmo del pulso y de la saturación de oxígeno y falsas alarmas. El ruido puede ser particularmente difícil de gestionar con sistemas de oximetría de pulso por reflexión, en los que los niveles de señal son mucho más pequeños que en los sistemas de oximetría de pulso por transmisión.
Los documentos US 2003/055325, US 6.002.952 y US 6.684.090 muestran métodos para la determinación del ritmo del pulso con un sistema de oximetría de pulso correspondiente. Sin embargo, ninguna de las técnicas anteriores enseña una solución simple para tener una estimación precisa del ritmo del pulso y la localización del pulso.
Sumario de la invención
Cada una de las funciones contenidas dentro de un conjunto preseleccionado de funciones se compara con la señal de entrada con muchos desplazamientos de tiempo diferentes, y se selecciona la combinación de función/ desplazamiento de tiempo que mejor coincide con la señal de entrada. La frecuencia de la función de mejor coincidencia se usa entonces como la mejor estimación de la frecuencia de la señal de entrada.
Breve descripción de los dibujos
La FIG. 1 es un diagrama de bloques que representa un método para determinar el ritmo del pulso de una señal pulsátil.
La FIG. 2A es un primer ejemplo de un conjunto de funciones de frecuencia predefinidas (S-PFF) para su uso con el método de la FIG. 1.
La FIG. 2B es un segundo ejemplo de un S-PFF para su uso con el método de la FIG. 1.
La FIG. 2C es un tercer ejemplo de un S-PFF para su uso con el método de la FIG. 1.
La FIG. 3 es un cuarto ejemplo de un S-PFF para su uso con el método de la FIG. 1.
La FIG. 4 es una señal de entrada y una función coincidente seleccionada a partir del S-PFF.
Descripción detallada de las realizaciones preferidas
Aunque las realizaciones preferidas divulgadas en el presente documento se describen en el contexto de oximetría de pulso, la presente invención puede usarse asimismo en otros contextos, incluyendo, pero sin limitarse a, procesamiento de señales para señales biológicas pulsátiles obtenidas a partir de fuentes distintas de las sondas de oximetría.
Para calcular el nivel de saturación de oxígeno (SaÜ2), a partir de las señales adquiridas en un sistema de oximetría de pulso, es frecuentemente preferible determinar primero el ritmo del pulso en las señales y localizar la localización del pulso. La obtención de una estimación precisa del ritmo del pulso y la localización del pulso es importante para calcular el nivel de saturación de oxígeno correcto en señales ruidosas.
5
10
15
20
25
30
35
40
45
50
55
60
65
La figura 1 es un diagrama de bloques de un planteamiento preferido para determinar la ritmo del pulso de una señal de entrada pulsátil 20 mediante la resolución de un problema de optimización que implica un Conjunto de Funciones de Frecuencia Predefinidas (S-PFF, del inglés “Set of Pre-Defined Frequency Functions”). La señal de entrada pulsátil 20 y un S-PFF 22 entran en un operador matemático 24 para procesamiento. El operador matemático 24 se emplea para transformar un problema de estimación de frecuencia en un problema de optimización. La salida 26 del sistema es el ritmo del pulso determinado para la señal de entrada.
La señal de entrada 20 puede ser una señal en bruto obtenida a partir de la fuente (por ejemplo, sonda de oximetría), o puede ser una preprocesada. En una realización preferida, los componentes de c.c. y de muy baja frecuencia de la señal de entrada se eliminan mediante un filtro de preprocesamiento apropiado antes de entrar en la etapa de procesamiento. La eliminación del componente de c.c. asegura que la señal procesada pulsa alrededor de cero y tiene picos positivos y negativos.
El S-PFF es un conjunto de N funciones periódicas que se usa en la determinación de la frecuencia de la señal de entrada. Preferentemente, las funciones del S-PFF se seleccionan para asemejarse, al menos aproximadamente, a las señales (sin ningún ruido añadido) que se esperan en la aplicación deseada. El intervalo de frecuencias de las funciones del S-PFF, el número N de funciones en un S-PFF dado, y el tamaño de los saltos entre frecuencias sucesivas se seleccionan preferentemente basándose en las señales de entrada esperadas y basándose también en la precisión requerida del resultado obtenido. La selección adelantada de funciones contenidas dentro de un S-PFF dado puede llevarse a cabo usando un proceso de entrenamiento adecuado basándose en un conjunto de datos que representa las señales que se encontrarán durante el uso posterior.
La FIG. 2A es un ejemplo de un primer S-PFF 30 que incluye un conjunto de N formas de onda periódicas 31-33, cada una teniendo una forma simétrica similar pero con diferentes frecuencias f - fN. La FIG. 2B es un ejemplo de un segundo S-PFF 40 que incluye un conjunto de N formas de onda periódicas 41-43, cada una teniendo una forma asimétrica similar pero con diferentes frecuencias f1 - fN. En realizaciones alternativas, puede usarse un S-PFF digital 50, tal como se muestra en la FIG. 2C. Este tercer S-PFF 50 incluye también un conjunto de N formas de onda periódicas f - fN 51-53, cada una teniendo un patrón similar pero con diferentes frecuencias. En estas realizaciones, las funciones 51-53 pueden tomar solamente 3 valores: {-1, 0, +1}. Una ventaja del uso de este tipo de S-PFF es que hace más simple y más rápida la realización de operaciones matemáticas, en comparación con un S-PFF analógico. Debe enfatizarse que los S-PFF 30, 40, 50 representados en las FIGS. 2A-C son meramente ejemplos de unos S-PFF que pueden usarse, y no han de interpretarse como limitativos de la presente invención.
En el contexto de la oximetría de pulso, la señal de entrada que se está procesando es una señal fisiológica pulsátil que se origina a partir de un sensor óptico, más grandes cantidades de ruido. El ritmo de la pulsación corresponde al ritmo cardiaco o ritmo del pulso del paciente. Dado que el ritmo del pulso esperado para casi todos los pacientes estará entre 30 y 240 latidos por minuto (lpm), un intervalo apropiado de frecuencias para un S-PFF usado para oximetría de pulso variaría desde ^ Hz (correspondiente a 30 lpm) a 4 Hz (correspondiente a 240 lpm). Dado que muchas aplicaciones médicas requieren que el ritmo cardiaco se determine con una precisión y resolución de 1 lpm, un S-PFF apropiado para oximetría de pulso sería un conjunto de 211 frecuencias diferentes (N=211) variando desde 30 lpm (^ Hz) a 240 lpm (4 Hz) en saltos de 1 lpm. Si se desea una mayor resolución o precisión, N puede incrementarse más allá de 211. Opcionalmente, los saltos entre las diversas frecuencias pueden disponerse no linealmente con saltos de lpm más pequeños a frecuencias más bajas y saltos mayores en las frecuencias más altas (por ejemplo, en una progresión logarítmica para conseguir una precisión en porcentaje uniforme).
En oximetría, la forma de las señales fisiológicas pulsátiles en pacientes saludables tendrá normalmente un tiempo de elevación más rápido y un tiempo de caída más lento. Como resultado, un S-PFF adecuado para oximetría de pulso sería el 40 mostrado en la FIG. 2B, que también tiene un tiempo de elevación más rápido y un tiempo de caída más lento y por lo tanto coincide aproximadamente con la señal esperada. Como se ha explicado anteriormente, 211 frecuencias (es decir, N=211) variando desde ^ Hz a 4 Hz en saltos de 1 lpm es un ejemplo de una distribución de frecuencias apropiada para aplicaciones de oximetría de pulso. En algunos casos, sin embargo, la forma de la forma de onda pulsátil fisiológica diferirá (por ejemplo, basándose en ciertas condiciones fisiológicas), lo que puede dar lugar a un rendimiento degradado cuando la forma de onda de entrada no coincide claramente con las formas de onda en el S-PFF. Es por lo tanto ventajoso incorporar diversas formas diferentes de onda dentro del S-PFF para tener en cuenta dichas diferencias.
La FIG. 3 es un ejemplo de un S-PFF 60 preferido para su uso en aplicaciones de oximetría de pulso que se utiliza para manejar señales de entrada con diferentes formas de onda. Este S-PFF 60 incluye un subconjunto de funciones 61-63 en 211 diferentes frecuencias que varían desde 30 lpm a 240 lpm con tiempos de elevación más rápidos y tiempos de caída más lentos, para coincidir con las señales fisiológicas pulsátiles esperadas de la mayor parte de pacientes saludables. Este S-PFF 60 también incluye un segundo subconjunto de funciones 64-66 en cada una de las 211 diferentes frecuencias con tiempos de elevación más lentos y tiempos de caída más rápidos, para coincidir con señales pulsátiles fisiológicas con formas correspondientes de pacientes cuyas formas de onda se desvían de la forma de onda de elevación rápida más común. Este S-PFF 60 incluye también un tercer subconjunto de funciones 67-69 en cada una de las 211 diferentes frecuencias con tiempos de elevación y tiempos de caída coincidentes para que coincidan con señales fisiológicas pulsátiles con formas correspondientes a pacientes cuyas
5
10
15
20
25
30
35
40
45
50
55
60
formas de onda se desvían de la norma. De ese modo, el S-PFF 60 preferido contiene un total de 633 funciones (3 * 211). Obsérvese que una duración adecuada para cada una de las funciones 61-69 dentro del S-PFF 60 para su uso en oximetría de pulso es de 3-5 segundos de duración.
Otro S-PFF preferido para su uso en aplicaciones de oximetría de pulso incluye un subconjunto de funciones 61-63 a 211 frecuencias diferentes que varían desde 30 lpm a 240 lpm con tiempos de elevación más rápidos y tiempos de caída más lentos (para coincidir con las señales fisiológicas pulsátiles esperadas desde la mayor parte de pacientes saludables), y un segundo subconjunto de funciones 67-69 a cada una de las 211 frecuencias diferentes con tiempos de elevación y tiempos de caída coincidentes. Este S-PFF preferido contiene por lo tanto un total de 422 funciones (2 * 211). 3-5 segundos es también una duración adecuada para cada una de las funciones 61-63, 67-69 dentro de este S-PFF.
Opcionalmente, pueden añadirse subconjuntos adicionales de funciones (no mostrado) para coincidir con condiciones fisiológicas específicas. Por ejemplo, un subconjunto adicional de 211 funciones que coinciden con las señales fisiológicas pulsátiles esperadas de pacientes con una válvula mitral deficiente puede añadirse a uno de los S-PFF explicados anteriormente (la forma de onda esperada de dichos pacientes difiere de la de pacientes saludables debido a que el flujo sanguíneo va inicialmente desde el ventrículo izquierdo a la aurícula izquierda. Solo después de que se acumule la presión suficiente se abre la válvula aórtica.)
Opcionalmente, pueden usarse diferentes S-PFF para diferentes clases de pacientes tales como adultos, niños, bebés y fetos. Por ejemplo, el intervalo de frecuencia del S-PFF puede desplazarse a frecuencias más altas para bebés (en comparación con el intervalo anteriormente descrito de 30-240 lpm).
Volviendo ahora a la FIG. 1, la señal de entrada 20 contiene la señal fisiológica pulsátil subyacente más ruido. La señal de entrada 20 y el S-PFF 22 entran en una etapa de procesamiento, que incluye un operador matemático F 24 que se diseña para estimar la frecuencia de la señal de entrada. El operador matemático F convierte el problema de estimación de frecuencia en un problema de optimización, y su implementación específica depende de la aplicación disponible. Tal como es el caso con el S-PFF, el operador matemático F puede optimizarse y ensayarse en una sesión de entrenamiento empleada sobre un conjunto de datos típicos.
En una realización preferida de la invención, el operador matemático F se define como el siguiente operador no lineal:
(1) F = MAX j {a,*/i(S-PFFj (t), I(t-i)) + a2*/2(S-PFFj(t), I(t-i))
+ a3*/3(S-PFFj(t),I(t-i)) ... + an*/n(S-PFFj(t),I(t-i))}
En la que:
ft-fn es un conjunto de n funciones lineales o no lineales;
S-PFF es un conjunto de funciones de frecuencia predefinidas, tal como se ha detallado anteriormente; j es un índice que indica una función S-PFF seleccionada de entre las N funciones S-PFF definidas (nótese que n es independiente de N);
ai-an son constantes de ponderación seleccionadas;
* es el operador de multiplicación;
I es la señal de entrada;
t es un índice que indica una muestra específica de S-PFFj e I (es decir estas funciones en el tiempo t); y i es un valor de desplazamiento seleccionado, que indica el desplazamiento actual entre la función S-PFFj y la señal de entrada I.
El operador matemático F representa un problema de optimización, en el que se busca un máximo, tal como se describe por la ecuación (1). En esta forma, cada una de las funciones contenidas dentro del S-PFF se compara con la señal de entrada en muchos desplazamientos de tiempo diferentes, y se selecciona la combinación función/ desplazamiento de tiempo que mejor coincide con la señal de entrada (es decir, mediante la selección de la combinación en la que ocurre el máximo). Una granularidad adecuada para el desplazamiento de tiempo, entre las diversas comparaciones es un desplazamiento que corresponde a una muestra. La frecuencia de la función de mejor coincidencia se usa entonces como la mejor estimación de la frecuencia de la señal de entrada. Por ejemplo, en la resolución de F el máximo se halló a (j=12, i=5), esto indicaría que el S-PFF12 es la función que proporciona el resultado máximo. Así la frecuencia de la función S-PFF12 se seleccionaría como la frecuencia de la señal de entrada I para producir la salida. En el contexto de la oximetría de pulso, la frecuencia será el ritmo del pulso (es decir, el ritmo cardiaco).
Además de optimizar el S-PFF para las señales esperadas, como se ha descrito anteriormente, las funciones fi-fn y constantes de ponderación ai-an, que son parte de F pueden definirse y optimizarse en una etapa de entrenamiento para encajar mejor en una aplicación específica y las características de las señales esperadas en esa aplicación. La
5
10
15
20
25
30
35
40
45
50
optimización puede realizarse usando un método de optimización estándar adecuado tal como MLS (mínimos cuadrados) o máxima probabilidad.
En una realización preferida que es adecuada para determinar el ritmo del pulso en sistemas de oximetría de pulso usando las 633 funciones S-PFF 60 preferidas (mostradas en la figura 3), n = 3 de modo que el operador matemático F incluye solamente tres funciones fi, f2 y Z3. En esta realización, la función fi(i, j) es la suma de los productos, obtenida mediante multiplicación de las funciones S-PFFj y la señal de entrada I, mientras se desplazan en i muestras, y la constante de ponderación ai es +1, dando como resultado la siguiente ecuación: fi(i, j) = £tS-PFFj(t)*I(t-i). La función f2(i, j) es la suma de diferencias entre los valores I(tk) consecutivos de la señal de entrada I y la constante de ponderación a2 es -1. Los valores {tk} representan los índices de los máximos de la función S-PFFj(t). La función f¡(i, j) es la suma de las diferencias entre los valores consecutivos I(tm), de la señal de entrada y la constante de ponderación a-¡ es también -1. Los valores {tm} representan los índices de los mínimos de la función S-PFFj(t). Con este conjunto particular de tres funciones, fi es principalmente responsable de la estimación de la frecuencia o del ritmo del pulso de la señal de entrada. Las funciones f2 y f3 son responsables principalmente de impedir que las frecuencias armónicas de la señal de entrada interfieran con la estimación de la frecuencia.
Una vez se selecciona la mejor función dentro del S-PFF, los índices {tk} representan la localización del máximo de los pulsos mientras los índices {tm} representan la localización de los mínimos de los pulsos. Mediante la resta los valores de las señales adquiridas en los índices {tm} de los valores en los índices {tk}, puede obtenerse la amplitud de los pulsos. Estos valores pueden usarse para el cálculo del nivel de saturación de oxígeno usando métodos de oximetría estándar.
Opcionalmente, el ritmo del pulso y/o las localizaciones de los máximos y mínimos de los pulsos pueden usarse para recoger un segmento particular de la onda pulsátil para procesamiento adicional. Por ejemplo, en un sistema de oximetría de pulso, una vez se identifican los máximos, mínimos y ritmo del pulso, pueden localizarse las diversas fases del ciclo cardíaco (dado que la parte de elevación corresponde al inicio del ciclo cardiaco, y la parte de descenso corresponde al final del ciclo cardiaco). Se deduce que los resultados mejorados se obtienen cuando las muestras que corresponden a la parte de elevación se usan para los cálculos de SaO2.
Pasando ahora a la FIG. 4, este planteamiento puede extenderse adicionalmente tomando un gran número de muestras durante cada ciclo cardiaco. En el ejemplo ilustrado la señal incidente en la entrada 70 tiene una frecuencia de 1 Hz, que corresponde a 60 lpm, y tiene una fase de elevación relativamente rápida y una fase de caída más lenta. Mediante la aplicación de las técnicas de procesamiento de señal descritas anteriormente usando el S-PFF 60 (mostrado en la FIG. 3), debería seleccionarse la coincidencia de la función 75 de 60 lpm (1 Hz) como se ha descrito anteriormente. Si el muestreo es suficientemente rápido para proporcionar 12 mediciones de SaO2 por segundo, la superposición de la muestra 71 de la señal de entrada 70 sobre la función 75 seleccionada a partir de los S-PFF permite la diferenciación entre las diversas fases del ciclo cardiaco. Por ejemplo, pueden identificarse muestras en la fase de elevación de la señal de entrada 70 mediante la selección de las muestras que corresponden a la parte de elevación 76 de la función 75 seleccionada. En el ejemplo ilustrado, la muestra 71' se identificaría como correspondiente a la parte de elevación de la señal de entrada 70. Esa muestra 71' se usaría posteriormente para realización de los cálculos de SaO2.
Las realizaciones anteriormente descritas proporcionan un planteamiento preciso y eficiente para determinar la frecuencia de una señal pulsátil de entrada. Este planteamiento puede suceder con éxito en el cálculo del ritmo del pulso incluso para valores de SNR (señal a ruido) muy bajos, lo que tiene lugar cuando la amplitud de la señal reflejada adquirida es baja y el nivel de ruido es alto. Aunque el planteamiento se ha descrito anteriormente en el contexto de la determinación del ritmo cardiaco y de las localizaciones del pulso a partir de señales derivadas de oximetría de pulso, debería reconocerse que estos planteamientos pueden usarse también para procesar otras señales fisiológicas pulsantes. Más aún, las funciones, frecuencia, resolución, duración, e intervalos de frecuencia expuestos anteriormente son meramente ejemplos, y como tales no han de ser interpretados como limitativos de la presente invención.

Claims (19)

  1. 5
    10
    15
    20
    25
    30
    35
    40
    45
    50
    55
    60
    REIVINDICACIONES
    1. Un método de determinación del ritmo del pulso de un sujeto vivo a partir de una señal de entrada de oximetría de pulso, comprendiendo el método las etapas de:
    comparar la señal de entrada para cada uno de una pluralidad de diferentes desplazamientos de tiempo para cada elemento de un conjunto de funciones,
    (i) en donde el conjunto incluye funciones con una pluralidad de diferentes frecuencias entre 30 a 240 ciclos por minuto;
    (ii) en donde el conjunto contiene un primer subconjunto que incluye funciones en las que el tiempo de elevación es más corto que el tiempo de caída y un segundo subconjunto que incluye funciones en las que (a) el tiempo de elevación es aproximadamente el mismo que el tiempo de caída O (b) el tiempo de elevación es más largo que el tiempo de caída,
    (iii) en donde la forma de todas las funciones en el primer subconjunto es similar, y
    (iv) en donde la forma de todas las funciones en el segundo subconjunto es similar;
    seleccionar, basándose en los resultados obtenidos en la etapa de comparación, una función a partir del conjunto que mejor coincide con la señal de entrada; y
    usar la frecuencia de la función seleccionada como una estimación del ritmo del pulso del sujeto.
  2. 2. El método de la reivindicación 1, que comprende adicionalmente la etapa de eliminar los componentes de c.c. y de muy baja frecuencia de la señal de entrada, en donde la etapa de eliminación se implementa antes de la etapa de comparación.
  3. 3. El método de la reivindicación 1, en el que cada una de las funciones del conjunto tiene una duración de entre aproximadamente 3 y aproximadamente 5 segundos.
  4. 4. El método de la reivindicación 3, en el que el conjunto contiene al menos 211 funciones.
  5. 5. El método de la reivindicación 1, en el que el conjunto contiene un primer subconjunto de funciones en las que el tiempo de elevación es más corto que el tiempo de caída y la forma de todas los funciones en el primer subconjunto es similar.
  6. 6. El método de la reivindicación 1, en el que las frecuencias de las funciones en el primer subconjunto y en el segundo subconjunto se distribuyen en frecuencia para proporcionar resolución de al menos un ciclo por minuto.
  7. 7. El método de la reivindicación 6, en el que el conjunto contiene adicionalmente un tercer subconjunto de funciones en las que el tiempo de elevación es más largo que el tiempo de caída, y la forma de todas las funciones en el tercer subconjunto es similar.
  8. 8. El método de la reivindicación 1, en el que las etapas de comparación y de selección se implementan basándose en la fórmula
    imagen1
    en la que:
    fi-fn es un conjunto de n funciones lineales o no lineales;
    S-PFF es un conjunto de funciones de frecuencia predefinidas;
    j es un índice que indica una función S-PFF seleccionada de entre las N funciones S-PFF definidas; ai-an son constantes de ponderación seleccionadas;
    * es el operador de multiplicación;
    I es la señal de entrada;
    t es un índice que indica una muestra específica de S-PFFj e I (es decir estas funciones en el tiempo t); y i es un valor de desplazamiento seleccionado, que indica el desplazamiento actual entre la función S-PFFj y la señal de entrada I.
  9. 9. El método de la reivindicación 1, en el que las etapas de comparación y de selección se implementan basándose en la fórmula
    F = MAX y {/XS-PFFj(t), I(t-i)) -/2(S-PFFj(t), I(t-i)) -/3(S-PFFj(t), I(t-i))}
    5
    10
    15
    20
    25
    30
    35
    40
    45
    50
    55
    60
    en la que:
    la función f1(i, j) es la suma de los productos, obtenida mediante multiplicación de las funciones S-PFFj y la señal de entrada I, mientras se desplazan en i muestras, dando como resultado la ecuación: fi(i, j) = £tS-PFFj(t)*I(t-i). la función f2(i, j) es la suma de diferencias entre los valores I(tk) consecutivos de la señal de entrada I, donde los valores {tk} representan los índices de los máximos de la función S-PFFj(t);
    la función /3(/, j) es la suma de las diferencias entre los valores consecutivos I(tm), de la señal de entrada, donde los valores {tm} representan los índices de los mínimos de la función S-PFFj(t);
    S-PFF es un conjunto de funciones de frecuencia predefinidas;
    j es un índice que indica una función S-PFF seleccionada de entre las N funciones S-PFF definidas;
    I es la señal de entrada;
    t es un índice que indica una muestra específica de S-PFFj e I (es decir estas funciones en el tiempo t); y i es un valor de desplazamiento seleccionado, que indica el desplazamiento actual entre la función S-PFFj y la señal de entrada I.
  10. 10. Un método de procesamiento de una señal de entrada de oximetría de pulso, comprendiendo el método las etapas de:
    comparar la señal de entrada para cada uno de una pluralidad de diferentes desplazamientos de tiempo para cada elemento de un conjunto de funciones,
    (i) en donde el conjunto incluye funciones con una pluralidad de diferentes frecuencias entre 30 a 240 ciclos por minuto;
    (ii) en donde el conjunto contiene un primer subconjunto que incluye funciones en las que el tiempo de elevación es más corto que el tiempo de caída y un segundo subconjunto que incluye funciones en las que (a) el tiempo de elevación es aproximadamente el mismo que el tiempo de caída O (b) el tiempo de elevación es más largo que el tiempo de caída,
    (iii) en donde la forma de todas las funciones en el primer subconjunto es similar, y
    (iv) en donde la forma de todas las funciones en el segundo subconjunto es similar;
    seleccionar, basándose en los resultados obtenidos en la etapa de comparación, una función a partir del conjunto que mejor coincide con la señal de entrada;
    identificar la parte de elevación de la señal de entrada basándose en una alineación en el tiempo con la parte de elevación de la función seleccionada; y
    calcular un nivel de saturación de oxígeno basándose en una o más muestras tomadas a partir de la parte de la señal de entrada que se identificó en la etapa de identificación.
  11. 11. El método de la reivindicación 10, en el que la señal de entrada tiene una tasa de muestreo de al menos 12 muestras por segundo.
  12. 12. El método de la reivindicación 10, que comprende adicionalmente la etapa de eliminar los componentes de c.c. y de muy baja frecuencia de la señal de entrada, en donde la etapa de eliminación se implementa antes de la etapa de comparación.
  13. 13. El método de la reivindicación 10, en el que cada una de las funciones del conjunto tiene una duración de entre aproximadamente 3 y aproximadamente 5 segundos.
  14. 14. El método de la reivindicación 10, en el que el conjunto contiene al menos 211 funciones.
  15. 15. El método de la reivindicación 10, en el que al menos una función del primer subconjunto y al menos una función del segundo subconjunto son de una frecuencia similar.
  16. 16. El método de la reivindicación 10, en el que las frecuencias de las funciones en el primer subconjunto y en el segundo subconjunto se distribuyen en frecuencia para proporcionar resolución de al menos un ciclo por minuto.
  17. 17. El método de la reivindicación 16, en el que el conjunto contiene adicionalmente un tercer subconjunto de funciones en las que el tiempo de elevación es más largo que el tiempo de caída y la forma de todas las funciones en el tercer subconjunto es similar.
  18. 18. El método de la reivindicación 10, en el que las etapas de comparación y de selección se implementan basándose en la fórmula
    F = MAX y ía,*/,(S-PFFj(t), I(t-i)) + a2*/2(S-PFFj(t), I(t-i))
    + a3*/3(S-PFFj(t), I(t-i)) ... + an*/n(S-PFFj(t), I(t-i))},
    5
    10
    15
    20
    25
    30
    en la que:
    fi-fn es un conjunto de n funciones lineales o no lineales;
    S-PFF es un conjunto de funciones de frecuencia predefinidas;
    j es un índice que indica una función S-PFF seleccionada de entre las N funciones S-PFF definidas; ai-an son constantes de ponderación seleccionadas;
    * es el operador de multiplicación;
    I es la señal de entrada;
    t es un índice que indica una muestra específica de S-PFFj e I (es decir estas funciones en el tiempo t); y i es un valor de desplazamiento seleccionado, que indica el desplazamiento actual entre la función S-PFFj y la señal de entrada I.
  19. 19. El método de la reivindicación 10, en el que las etapas de comparación y de selección se implementan basándose en la fórmula
    F = MAX u {/l(S-PFFj(t), I(t-i)) -/2(S-PFFj(t), l(t-i)) -/3(S-PFFj(t), I(t-i))},
    en la que:
    la función f1(i, j) es la suma de los productos, obtenida mediante multiplicación de las funciones S-PFFj y la señal de entrada I, mientras se desplazan en i muestras, dando como resultado la ecuación: fi(i, j) = £tS-PFFj(t)*I(t-i). la función f2(i, j) es la suma de diferencias entre los valores I(tk) consecutivos de la señal de entrada I, donde los valores {tk} representan los índices de los máximos de la función S-PFFj(t);
    la función Z3(/, j) es la suma de las diferencias entre los valores consecutivos I(tm), de la señal de entrada, donde los valores {tm} representan los índices de los mínimos de la función S-PFFj(t);
    S-PFF es un conjunto de funciones de frecuencia predefinidas;
    j es un índice que indica una función S-PFF seleccionada de entre las N funciones S-PFF definidas;
    I es la señal de entrada;
    t es un índice que indica una muestra específica de S-PFFj e I (es decir estas funciones en el tiempo t); y i es un valor de desplazamiento seleccionado, que indica el desplazamiento actual entre la función S-PFFj y la señal de entrada I.
ES06815720.5T 2005-09-30 2006-09-29 Procesamiento de señal para oximetría de pulso Active ES2662198T3 (es)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US72225705P 2005-09-30 2005-09-30
US722257P 2005-09-30
PCT/US2006/037936 WO2007041252A2 (en) 2005-09-30 2006-09-29 Signal processing for pulse oximetry

Publications (1)

Publication Number Publication Date
ES2662198T3 true ES2662198T3 (es) 2018-04-05

Family

ID=37906721

Family Applications (1)

Application Number Title Priority Date Filing Date
ES06815720.5T Active ES2662198T3 (es) 2005-09-30 2006-09-29 Procesamiento de señal para oximetría de pulso

Country Status (10)

Country Link
US (2) US8005523B2 (es)
EP (1) EP1928309B1 (es)
JP (1) JP5203210B2 (es)
CN (1) CN101282685B (es)
AU (1) AU2006297202B2 (es)
CA (1) CA2623270C (es)
ES (1) ES2662198T3 (es)
HK (1) HK1125018A1 (es)
IL (2) IL189943A (es)
WO (1) WO2007041252A2 (es)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009064979A2 (en) * 2007-11-14 2009-05-22 Conmed Corporation Method and apparatus for processing a pulsatile biometric signal
EP2303121B1 (en) * 2008-07-11 2012-06-20 Koninklijke Philips Electronics N.V. Automatic transmission power adjustment for doppler radar
US8216136B2 (en) 2009-03-05 2012-07-10 Nellcor Puritan Bennett Llc Systems and methods for monitoring heart rate and blood pressure correlation
US8636667B2 (en) 2009-07-06 2014-01-28 Nellcor Puritan Bennett Ireland Systems and methods for processing physiological signals in wavelet space
US8494786B2 (en) 2009-07-30 2013-07-23 Covidien Lp Exponential sampling of red and infrared signals
JP6393984B2 (ja) * 2013-12-19 2018-09-26 富士通株式会社 脈拍計測装置、脈拍計測方法及び脈拍計測プログラム
JP6488722B2 (ja) * 2015-01-26 2019-03-27 富士通株式会社 脈波検出装置、脈波検出方法及び脈波検出プログラム
JP6309659B1 (ja) * 2017-01-27 2018-04-11 フクダ電子株式会社 生体信号解析装置およびその制御方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4193393A (en) * 1977-08-25 1980-03-18 International Bio-Medical Industries Diagnostic apparatus
US4960126A (en) * 1988-01-15 1990-10-02 Criticare Systems, Inc. ECG synchronized pulse oximeter
US6002952A (en) 1997-04-14 1999-12-14 Masimo Corporation Signal processing apparatus and method
US6463311B1 (en) * 1998-12-30 2002-10-08 Masimo Corporation Plethysmograph pulse recognition processor
US6684090B2 (en) 1999-01-07 2004-01-27 Masimo Corporation Pulse oximetry data confidence indicator
IL130371A (en) * 1999-06-08 2004-06-01 Oridion Medical Ltd Capnography waveform interpreter
US6850787B2 (en) * 2001-06-29 2005-02-01 Masimo Laboratories, Inc. Signal component processor
IL145445A (en) * 2001-09-13 2006-12-31 Conmed Corp A method for signal processing and a device for improving signal for noise
DE60143115D1 (de) 2001-09-28 2010-11-04 Csem Ct Suisse Electronique Verfahren und Vorrichtung zur Pulsmessung
US6650918B2 (en) 2002-02-22 2003-11-18 Datex-Ohmeda, Inc. Cepstral domain pulse oximetry
US6763256B2 (en) 2002-08-16 2004-07-13 Optical Sensors, Inc. Pulse oximeter
US7016715B2 (en) 2003-01-13 2006-03-21 Nellcorpuritan Bennett Incorporated Selection of preset filter parameters based on signal quality
US7341560B2 (en) * 2004-10-05 2008-03-11 Rader, Fishman & Grauer Pllc Apparatuses and methods for non-invasively monitoring blood parameters

Also Published As

Publication number Publication date
IL189943A (en) 2012-04-30
AU2006297202A1 (en) 2007-04-12
CA2623270A1 (en) 2007-04-12
HK1125018A1 (en) 2009-07-31
EP1928309B1 (en) 2018-01-17
EP1928309A4 (en) 2012-10-10
US8005523B2 (en) 2011-08-23
CA2623270C (en) 2014-07-08
JP5203210B2 (ja) 2013-06-05
AU2006297202B2 (en) 2011-12-22
CN101282685A (zh) 2008-10-08
CN101282685B (zh) 2012-04-25
IL189943A0 (en) 2008-08-07
US20110270059A1 (en) 2011-11-03
EP1928309A2 (en) 2008-06-11
IL213211A0 (en) 2011-07-31
JP2009509668A (ja) 2009-03-12
IL213211A (en) 2012-06-28
WO2007041252A2 (en) 2007-04-12
WO2007041252A3 (en) 2007-11-15
US20070203417A1 (en) 2007-08-30

Similar Documents

Publication Publication Date Title
ES2662198T3 (es) Procesamiento de señal para oximetría de pulso
US9814418B2 (en) Sine saturation transform
RU2498768C2 (ru) Устройство и способ для определения физиологического параметра
US8507281B2 (en) Method of removing abnormal data and blood component spectroscopy analysis system employing the same
EP2552303B1 (en) Sensor device and method for measuring and determining a pulse arrival (pat) time
RU2010146145A (ru) Неинвазивный способ и устройство для определения стадий поверхностного и глубокого сна
US20150105666A1 (en) Narrow band feature extraction from cardiac signals
ES2276594B1 (es) Metodo para procesar señales fotopletismograficas obtenidas de una persona o animal, y oximetro que utiliza dicho metodo.
RU2268639C2 (ru) Способ пульсометрической оценки функционального состояния и характера вегетативной регуляции сердечно-сосудистой системы человека
Aarthi et al. Fingertip based estimation of heart rate using photoplethysmography
US20190290174A1 (en) Biological information measurement device and non-transitory computer readable medium
JP2018068556A (ja) 胸腔内圧推定装置
JP2016140373A (ja) 脈波計測装置、および脈波計測方法
US20190223731A1 (en) Transesophageal or transtracheal cardiac monitoring by optical spectroscopy
Mathieu et al. Using smart wearables to monitor cardiac ejection
JP7170459B2 (ja) 血行検出装置、方法、及びプログラム
Figini Development of a cuff-less Blood monitoring device
Huotari et al. Arterial pulse waves measured with EMFi and PPG sensors and comparison of the pulse waveform spectral and decomposition analysis in healthy young and elderly subjects
Chaurasia et al. Development of a Low Cost Heart Rate Monitoring and Transmission System using PPG Signal Processing for Wearable Devices
Ackshaiyaa et al. Continuous monitoring of blood pressure from photophlethysmograph signal
Andruschenko et al. A novel algorithm for “ac/dc” ratio calculation of pulse signals using energy criteria of the pulse
RU2575308C2 (ru) Автоматический способ измерения и обработки кровяного давления
Huotari et al. Photoplethysmographic Measurements of Finger/Toe Arterial Pulse Waveforms and Their Compound Time Domain Analysis