Способ и устройство анализа баллистокардиографических сигналов



Способ и устройство анализа баллистокардиографических сигналов
Способ и устройство анализа баллистокардиографических сигналов
Способ и устройство анализа баллистокардиографических сигналов
Способ и устройство анализа баллистокардиографических сигналов
Способ и устройство анализа баллистокардиографических сигналов
Способ и устройство анализа баллистокардиографических сигналов
Способ и устройство анализа баллистокардиографических сигналов
Способ и устройство анализа баллистокардиографических сигналов
Способ и устройство анализа баллистокардиографических сигналов
Способ и устройство анализа баллистокардиографических сигналов
Способ и устройство анализа баллистокардиографических сигналов
Способ и устройство анализа баллистокардиографических сигналов

 


Владельцы патента RU 2517583:

КОНИНКЛЕЙКЕ ФИЛИПС ЭЛЕКТРОНИКС Н.В. (NL)

Группа изобретений относится к медицине, а именно к кардиологии. Определяют оценочное начальное временное значение для первого сокращения сердца в баллистокардиографическом сигнале. Итерационно вычисляют оценочные значения для последующих сокращений сердца в баллистокардиографическом сигнале, используя оценочное начальное временное значение, при этом каждая итерация на этапе вычисления содержит оценку целевой функции, которая содержит взвешенную сумму множества оценочных функций. Причем каждый итерационный этап при вычислении оценочных значений для последовательных сердечных сокращений в баллистокардиографическом сигнале ограничен целевым интервалом после оценочного временного значения, найденного на предшествующем итерационном этапе вычисления. Для реализации способа используют устройство и машиночитаемый носитель, содержащий управляющую компьютерную программу. Группа изобретений позволяет проводить оценку баллистографического сигнала с помощью минимального количества датчиков, при этом не теряя точности диагностики, кроме того, обладает высокой чувствительностью для идентификации паттернов нерегулярности, например при аритмиях, и надежна при работе с артефактами, вызванными движением. 3 н. и 18 з.п. ф-лы, 10 ил.

 

Область техники, к которой относится изобретение

Изобретение относится к способу и устройству анализа баллистокардиографических сигналов, в частности к способу и устройству, предусматривающим обнаружение отдельных событий сокращения сердца в баллистокардиографических сигналах.

Уровень техники изобретения

Баллистокардиограф (BCG) измеряет перемещение тела человека, вызванное динамикой крови по мере ее прокачивания сердцем.

BCG имеет преимущества над электрокардиографом (ECG) в том, что измерение жизненно важных сигналов возможно без приклеивания электродов к телу или использования специальных датчиков в виде поясов, трикотажных изделий и т.п., которые требуется на себе носить. Он в особенности пригоден для измерения частоты пульса и данных по изменению частоты пульса для оценки качества сна, стрессового состояния и сердечной деятельности. В этих областях применения ненавязчивость BCG-мониторинга представляет собой первостепенную важность, поскольку датчики, находящиеся в непосредственном контакте с пациентом, неизбежно снижают качество сна.

В настоящее время алгоритмы анализа баллистокардиографических сигналов с целью определения частоты сокращений сердца используют спектральные методы либо методы определения во времени для обнаружения многократного проявления определенных паттернов, например, путем оценки автокорреляционной функции сигнала. При всех таких подходах следует рассматривать сегменты сигнала с продолжительностью в несколько секунд, так чтобы охватывать множество сокращений сердца. В результате доступна информация о среднем показателе сокращений сердца на некотором отрезке времени, однако точная информация в интервале между последовательными сокращениями сердца недоступна.

Представлены некоторые алгоритмы оценки последовательных сокращений сердца на основе баллистокардиографических сигналов, но такие подходы либо требуют большого и дорогостоящего массива датчиков для должной работы ("FFT averaging of multichannel BCG signals from bed mattress sensor to improve estimation of heart beat interval" Kortelainen J.M. и Virkkala J., Engineering in Medicine and Biology Society, 2007, EMBS 2007, 29-ая Ежегодная Международная конференция по IEEE, 22-26 Август 2007, страницы 6685-6688), либо вмешательства человека ("Automatic Ballistocardiogram (BCG) Beat Detection Using a Template Matching Approach" J.H. Shin, B. H. Choi, Y. G. Lim, D.U. Joeng и K.S. Park, Engineering in Medicine and Biology Society, 2008, EMBS 2008, 30-ая Ежегодная Международная конференция по IEEE, 21-24 Август 2008), либо использования различных схем датчиков при недостаточной точности измерений ("Estimation of Respiratory Waveform and Heart Rate Using an Accelerometer" D. H. Phan, S. Bonnet, R. Guillemaud, E. Castelli, N.Y. Pham Thi, Engineering in Medicine and Biology Society, 2008, EMBS 2008, 30-ая Ежегодная Международная конференция по IEEE, 21-24 Август 2008).

Остаются вопросы в отношении того, возможно ли использование этих алгоритмов на рынке, а также могут ли они работать в условиях высокой вариабильности баллистокардиографических сигналов, связанных с организмом конкретного пациента и различиями между пациентами.

Разумеется, можно получить оценочную информацию в интервале между последовательными сокращениями сердца стандартными способами с использованием электрокардиограммы (ECG), однако, как говорилось выше, такие способы требуют, чтобы пациент был "опутан проводами", а значит, они причиняют неудобства.

Как уже упоминалось, способы существующего уровня техники для расчета частоты сердечных сокращений на основе баллистокардиографических сигналов, полученных с одного датчика, могут дать оценку частоты сердечных сокращений в интервале времени в несколько секунд. Приложения в области регистрации стрессового состояния во время сна или апноэ сна, сердечной недостаточности и т.д. требуют информации об изменении частоты сердечных сокращений между последовательными сокращениями сердца, т.е. усредненной оценки в течение некоторого периода времени недостаточно. Кроме того, существующие ныне способы требуют регулярности сердечных сокращений в данном интервале для проведения точного анализа. В этой связи наличие определенной аритмии, такой как эктопические сокращения или перебои в работе сердца, либо искажают оценку частоты сердечных сокращений, либо просто остаются незамеченными. Таким образом, может быть обнаружена лишь низкочастотная составляющая изменения частоты пульса.

Сущность изобретения

Согласно первому аспекту изобретения предложен способ анализа баллистокардиографического сигнала для определения частоты сердечных сокращений, при этом способ содержит определение оценочного начального временного значения для первого сокращения сердца в баллистокардиографическом сигнале; итерационное вычисление оценочных значений для последующих сокращений сердца в баллистокардиографическом сигнале, используя оценочное начальное временное значение; при этом каждая итерация на этапе вычисления содержит оценку целевой функции, которая содержит взвешенную сумму множества оценочных функций; причем каждый итерационный этап при вычислении оценочных значений для последующих сердечных сокращений в баллистокардиографическом сигнале ограничен целевым интервалом после оценочного временного значения, найденного на предшествующем итерационном этапе вычисления.

Согласно второму аспекту изобретения предложен компьютерный программный продукт, содержащий управляющую компьютерную программу, которая при выполнении на компьютере или процессоре выполнена с возможностью реализации этапов вышеописанного способа.

Согласно третьему аспекту изобретения предложено устройство для использования совместно с блоком измерения баллистокардиографического сигнала у пациента, при этом устройство содержит средство получения баллистокардиографического сигнала от блока; и средство обработки информации для реализации вышеописанного способа применительно к полученному баллистокардиографическому сигналу.

В отличие от предшествующего уровня техники способ и устройство согласно изобретению рассчитывают продолжительность временного интервала между последовательными сокращениями сердца, как в случае ECG. Таким образом, можно удовлетворить требованиям, выдвигаемым в вышеупомянутых областях применения. Кроме того, для получения требуемых сигналов необходим только один датчик. Не требуются ни дорогостоящее оборудование, содержащее множество датчиков, ни контроль со стороны специалиста. Способ и система работают устойчиво с BCG-сигналом, обычно обладающим высокой степенью сложности, поскольку используется информация о характеристических событиях, долгосрочный прогноз и априорная информация о продолжительности сердечного цикла. Таким образом, способ и устройство достаточно устойчиво работают в ситуациях, при которых сигнал искажается незначительными артефактами движения, и обладают достаточно высокой чувствительностью для идентификации паттернов проявления нерегулярности, подобных аритмиям.

Краткое описание чертежей

Далее изобретение будет описано, приводимое лишь в качестве примера, со ссылкой на следующие чертежи, на которых:

на Фиг.1 показан типовой баллистокардиографический сигнал;

на Фиг.2 показана блок-схема алгоритма, иллюстрирующая способ согласно изобретению;

на Фиг.3 более подробно показан этап 105, представленный на Фиг.2;

на Фиг.4 показан график, иллюстрирующий пример оценочной функции, характеризующей высокочастотные составляющие;

на Фиг.5 показана спектрограмма баллистокардиографического сигнала;

на Фиг.6 показан график, иллюстрирующий пример оценочной функции для долгосрочного прогнозирования;

на Фиг.7 показан график, иллюстрирующий пример оценочной функции для функции распределения вероятностей применительно к априорной продолжительности сердечных сокращений;

на Фиг.8 показан график, иллюстрирующий комбинированную функцию и ее максимум;

на Фиг.9 показан график, иллюстрирующий результаты этапа 105, представленного на Фиг.2; а также

на Фиг.10 показан график, иллюстрирующий результаты этапа 105 и уточненную оценку на этапе 107, представленном на Фиг.2.

Подробное описание изобретения

На Фиг.1 показана типовая баллистокардиограмма (BCG), полученная с использованием фольгового датчика, помещенного ниже грудной клетки пациента, лежащего на кровати или на столе. Баллистокардиограмма регистрирует дыхательные движения и периодический паттерн сердечных сокращений, связанных с сердечным ритмом.

Можно видеть, что BCG имеет доминирующую низкочастотную составляющую (продолжительностью около четырех секунд), которая связана с дыхательными движениями пациента, а также меньшие по размеру флуктуации, обладающие более высокой частотой, которые вызваны механической активностью сердца.

Следует понимать, что пациент должен находиться в состоянии покоя, чтобы получить такой ясно выраженный BCG-сигнал. Более значительные перемещения приведут к преобладанию артефактов движения в BCG-сигнале, что существенно осложнит его анализ.

Первый этап в обработке BCG-сигнала заключается в разбивке его на подинтервалы, на которых прочие движения и возмущения препятствуют оценке частоты сердечных сокращений и частоты дыхания, и области, в которых оценка возможна. Такого рода разбивка достигается путем оценки энергетического уровня сигнала.

Кроме того, вклады от дыхательных движений и от механической активности сердца можно разделить путем использования фильтров. Например, полосовая фильтрация с помощью фильтра Баттерворта третьего порядка, имеющего нижнюю граничную частоту в пределах 0,04-0,08 Гц, а верхнюю граничную частоту - в пределах 0,50-0,70 Гц, позволяет получить дыхательную составляющую. Составляющая, определяемая сердечными сокращениями, может быть выделена путем фильтрации с помощью фильтра верхних частот (например, фильтра Баттерворта второго порядка с частотой отсечки в диапазоне 0,8-1,2 Гц). В способе, описанном ниже, используется составляющая, обусловленная сердечными сокращениями.

Способ согласно изобретению будет кратко описан со ссылкой на блок-схему алгоритма на Фиг.2.

На этапе 101 осуществляют сбор сигналов с помощью BCG.

На этапе 103 вычисляется начальное значение времени проявления первого сердечного сокращения от начала баллистокардиографического сигнала. Расчет этого оценочного показателя будет подробнее описан ниже, и это служит отправной точкой итерационной процедуры (показанной на этапе 105), в которой итеративно вычисляются оценочные показатели сердечных сокращений с пошаговым продвижением по времени при каждой оценке вплоть до достижения окончания баллистокардиографического сигнала. Процедура, выполняемая на этапе 105, будет подробнее описана ниже со ссылкой на Фиг.3.

Результатом этапа 105 является первая сегментация баллистокардиографического сигнала на интервалы сердечных сокращений.

На этапе 107 выполняется процедура уточнения, в которой используются данные прогнозирования, полученные на этапе 105, для расчета продолжительности заключительного интервала сердечных сокращений. Реализация этой процедуры будет подробнее описана ниже.

Сначала подробнее рассмотрим этап 103 со ссылкой на Фиг.3.

Когда происходит событие сердечного сокращения, известно, что следующее сердечное сокращение последует в пределах определенного интервала времени. Таким образом, в предпочтительном варианте осуществления данное обстоятельство используется путем ограничения интервала поиска в отношении следующего сердечного сокращения разумными с физиологической точки зрения величинами. Итак, для данного оценочного значения tn сердечного сокращения, момент времени tn+1 следующего сокращения сердца должен находиться в интервале [tn+tmin, tn+tmax]. Для tmin и tmax разумно выбрать значения около 0,5 секунды (что соответствует частоте сердечных сокращений, составляющей примерно 120 ударов в минуту) и до 2 секунд (что соответствует частоте сердечных сокращений, составляющей примерно 30 ударов в минуту). Специалист в данной области техники поймет, что можно выбрать альтернативные величины, приемлемые с физиологической точки зрения.

Предпочтительно формируются три оценочные функции λ(t), μ(t) и σ(t) (этапы 123, 125 и 127) для оценки того, насколько удачно проведена возможная оценка в отношении tn+1 по критериям, которые предпочтительно касаются соответственно наличия характеристических высокочастотных составляющих (оценочная функция λ(t)), сердечного ритма на длительном отрезке времени (μ(t)) и вероятностного распределения длин интервалов сердечных сокращений (σ(t)).

На этапе 129 наилучшее оценочное значение согласно этим оценочным функциям принимается за оценочное значение tn+1 и служит основой для расчета tn+2. Интервал поиска для tn+2 далее задается в виде [tn+1+tmin, tn+1+tmax], а оценочные функции λ(t), μ(t) и σ(t) корректируются.

Путем повторения данной итерационной процедуры способ обеспечивает обработку баллистокардиографического сигнала слева направо вдоль временной оси, начиная с исходного значения в начальной части баллистокардиографического сигнала и до момента окончания баллистокардиографического сигнала.

Далее описывается способ вычисления оценочных функций λ(t), μ(t) и σ(t) при заданном оценочном значении tn сердечных сокращений, и как они используются для расчета следующего итерационного этапа tn+1.

Характеристическая высокочастотная составляющая λ(t): Баллистокардиографический сигнал подвергается полосовой фильтрации для выделения высокочастотного контента. Возведение в квадрат и низкочастотная фильтрация результирующего сигнала позволяют получить новый сигнал с превалирующими максимальными уровнями там, где в исходном сигнале имеют место высокочастотные составляющие. Данная процедура подробнее описана в патентной заявке под названием "Method and apparatus for the analysis of ballistocardiogram signals", зарегистрированной в то же время, что и настоящая заявка, и тем же заявителем. Сигнал, полученный этим способом, представляет собой оценочную функцию для проявления характеристических высокочастотных составляющих в полном баллистокардиографическом сигнале и, будучи ограниченным интервалом [tn+tmin, tn+tmax], может служить в качестве λ(t) для этого конкретного этапа итерации. На Фиг.4 показан пример оценочной функции λ(t), когда полосовой фильтр выделяет полосу частот от 20 до 40 Гц, а низкочастотный фильтр имеет частоту отсечки 3,5 Гц.

Долгосрочное прогнозирование μ(t): Исходным моментом для вычисления функции μ(t) является частотно-временное распределение с хорошей частотной локализацией.

Используя возведенные в квадрат абсолютные значения спектрограммы, получаем

где BCG(s) представляет зарегистрированный баллистокардиографический сигнал, а h(s) - окно Хенинга h(s) продолжительностью 15 секунд

На спектрограмме измеряют частотное содержание для частоты f в момент времени t. Если заданы tn и возможное событие сердечного сокращения в момент времени t, оба момента времени определяют частоту сердечного ритма, соответствующую (t-tn)-1. Если это соответствует правильной частоте сердечного ритма, спектрограмма покажет повышенное значение для S(t, (t-tn)-1), в противном случае оно будет меньше. Таким образом, оценочная функция μ(t) выбирается в следующем виде:

μ(t)=S(t, (t-tn)-1) для t∈[tn+tmin, tn+tmax](3)

Часть спектрограммы баллистокардиографического сигнала показана на Фиг.5. Здесь боле светлые участки соответствуют более высоким значениям энергии на спектрограмме, а линия представляет функцию t→(t-tn)-1.

На Фиг.6 показана оценочная функция μ(t) для долгосрочного прогнозирования, построенная с использованием спектрограммы. Значения функции на линии, представленной на Фиг.5, используются для расчета оценочной функции на Фиг.6.

Альтернативный подход к определению μ(t) заключается в использовании способа определения периодичности, например, использования автокорреляционной функции (ACF). При таком подходе рассчитывают автокорреляцию

где BCG(s) - баллистокардиографический сигнал, а 2ε - длительность окна анализа, которая обычно изменяется в пределах от 5 до 20 секунд.

Вероятностная функция распределения σ(t): Длины интервалов между двумя последовательными сердечными сокращениями подчиняются распределению Гаусса. Таким образом, в предпочтительном варианте осуществления оценочная функция σ(t) может быть выбрана в следующем виде:

В предпочтительном варианте осуществления среднее значение m и стандартное отклонение s следующие: m=0,92 и s=0,4.

В дополнительных вариантах осуществления для выбора m и s может рассматриваться история ранее обнаруженных сердечных сокращений. В альтернативном варианте могут применяться стохастические модели, в которых информация о продолжительности интервалов сердечных сокращений в прошлом может использоваться для прогнозирования продолжительности последующего интервала сердечных сокращений, и которые известны специалистам в данной области техники.

Пример оценочной функции σ(t) вероятностного распределения показан на Фиг.7.

Для более совершенного контроля над присваиванием весовых коэффициентов между этими тремя функциями полезно нормализовать оценочные функции в отношении максимум-нормы или отобразить диапазон значений функций в интервале [0, 1] при помощи соответствующего аффинного отображения. Далее на этапе 129 находят максимальное значение многочлена αλ(t)+βμ(t)+χσ(t), где α, β и χ - скалярные величины, для t∈[tn+tmin, tn+tmax]. Превалирующее максимальное значение этого многочлена должно быть найдено, чтобы решить проблему многоцелевой оптимизации. Данное максимальное значение предпочтительно определяется посредством низкочастотной фильтрации результата сложения вслед за поиском максимума, и оно показано на Фиг.8. Определенное таким образом максимальное значение становится следующим оценочным значением сердечного ритма tn+1, и итерационный процесс продолжается от этого оценочного значения, как уже было описано, до момента достижения окончания сигнала.

Следует понимать, что каждый из вышеприведенных этапов 123, 125 и 127, когда они используются по отдельности, предрасположен к ошибкам. Например, отыскание высокочастотной составляющей на этапе 123 для отдельных сердечных сокращений не даст результатов в присутствии артефактов движения (т.е. если пациент чрезмерно подвижен). В других случаях расчет средней частоты сердечных сокращений на этапе 125 может вызывать сложности, поскольку аритмия не учитывается должным образом или полностью игнорируется. Наконец, вероятностный подход на этапе 127 может быть предпочтительно использован, только если имеется информация о предыстории.

Таким образом, выходные данные на каждом из этапов 123, 125 и 127 объединяются в единую функцию на этапе 129, и эта функция исследуется на максимум, чтобы обеспечить надежную оценку сердечного ритма, определяемого на основе последовательного чередования ударов пульса.

Вернемся к Фиг.2 и рассмотрим этап 103 более подробно.

Имея некоторую произвольную начальную точку, в способе, представленном на Фиг.2, обычно потребуется проведение всего нескольких итераций для достижения устойчивого режима и получения точных оценочных значений сердечного ритма.

В альтернативном варианте начальная точка может определяться как максимальное значение оценочной функции для высокочастотной составляющей в начальной части сигнала. Ограничительный интервал для оценочной функции может составлять [0 c, 1,5 с].

В дополнительном альтернативном варианте оценочные значения для начального момента времени t1 и следующего итерационного этапа t2 могут вычисляться одновременно, так чтобы t1 представлял возможную начальную точку сегмента сердечных сокращений, а t2 - возможную конечную точку. Для каждого t1 в диапазоне [0 c, 1,5 с] определяется Q(t1, t2)=αλ(t2)+βμ(t2)+χσ(t2) для t2 ∈ [t1+tmin, t1+tmax]. Другими словами, Q(t1, t2) определяет, насколько удачен выбор t1 и t2 в отношении одной и той же оценочной функции, используемой для итерационного метода. Два значения аргумента, при которых Q достигает максимума, дают наилучшие оценки для моментов времени t1 и t2 соответственно первых двух сердечных сокращений.

В еще одном альтернативном варианте один из двух вышеописанных способов используется для расчета предварительной начальной точки, а далее несколько раз выполняются обычные итерации до точки tm итерации (где, например, m=15). С этого момента, однако, алгоритм работает "в обратном направлении", а это означает, что в качестве поискового интервала выбирается интервал [tm-tmax, tm-tmin] вместо интервала [tm+tmin, tm+tmax]. Итерационное вычисление оценочных значений сердечных сокращений продолжается "в обратном направлении" до тех пор, пока не будет снова достигнуто начало сигнала. Последний результат оценки "в обратном направлении" далее служит начальной точкой t1. Здесь используется тот факт, что алгоритм дает верные оценочные значения за несколько итераций, даже если начальная точка была выбрана неудачно.

После выполнения этапов 103 и 105 становится возможной первая сегментация баллистокардиографических сигналов, при которой каждый сегмент по существу содержит характеристический паттерн только одного отдельного сердечного сокращения. В конечном итоге данная сегментация может быть использована для предоставления начального значения для способов детализации, которые находят точные моменты времени событий сердечных сокращений или продолжительность интервалов сердечных сокращений.

В одном альтернативном подходе уточненные оценочные значения могут вычисляться путем поиска характеристических признаков в баллистокардиографическом сигнале, близких к оценочным значениям, найденным на этапах 103 и 105. В баллистокардиографическом сигнале, который предпочтительно был подвергнут высокочастотной фильтрации, за точками, найденными с помощью этапов 103 и 105, следует волновой сигнал низкой частоты, но с относительно высокой амплитудой. Обнаружение локально наибольших максимумов в малой близлежащей области точек, найденных с помощью этапов 103 и 105, позволит провести более точную локализацию событий сердечных сокращений. После того как сердечные сокращения идентифицированы, можно рассчитать интервалы между последовательными сердечными сокращениями.

В еще одном альтернативном варианте имеется возможность использовать автокорреляционную функцию (ACF) или другие способы для определения периодичности во временном ряду для определения кардиоинтервала между двумя последовательными сердечными сокращениями. В этом случае этапы 103 и 105 уже предоставляют некоторую приблизительную информацию об интервале, в котором должна наблюдаться временная задержка Δt, для которой ACF достигает своего максимума (т.е. производится оценка не всех возможных значений Δt). Такое ограничение лишь имеющими значение интервалами для Δt требуется для обеспечения оценки сердечного ритма между последовательными сердечными сокращениями, выполняемой с помощью ACF. Это устраняет проблему, возникающую при рассмотрении всей теоретически возможной области поиска для Δt, поскольку ACF часто показывает побочные максимумы, не относящиеся к событиям сердечных сокращений.

Результаты процедуры сегментации можно видеть на Фиг.9. На Фиг.9 показан баллистокардиографический сигнал (после высокочастотной фильтрации) вместе с идентифицированными сегментами, а также события сердечных сокращений, обнаруженные в ECG-сигнале (помеченные крестами). При маркировке максимального значения сигнала в пределах 0,12 секунды после получения результата сегментации, характеристическое событие надежно определяется, как показано на Фиг.10.

Оценка продолжительности интервала между последовательными сокращениями сердца (без уточнений и с уточнениями соответственно) и сравнение с результатами ECG можно видеть в таблице 1, приведенной ниже.

Таблица 1
ECG 1,040 1,020 0,980 0,968 1,032 1,012 0,956 0,912 0,940
BCG
(не уточн.)
1,020 1,016 1,008 0,972 0,996 1,008 0,988 0,956 0,920
BCG (уточн.) 1,040 1,028 0,976 0,964 1,028 1,012 0,960 0,916 0,940

Специалисту в данной области техники станут очевидны многие конструктивные конфигурации и вариации способа. Трем различным оценочным функциям в целевой функции на этапе 129 могут быть присвоены весовые коэффициенты по-разному, чтобы адаптировать способ применительно к различным требованиям. Это можно осуществить путем соответствующего выбора скалярных величин α, β и χ. Например, сделав упор на долгосрочное прогнозирование по результатам выполнения этапа 125, можно повысить надежность оценки сердечных сокращений, когда ожидается регулярный сердечный ритм, а при необходимости обнаружить аритмию, будет лучше сосредоточить больше внимания на высокочастотных составляющих, определяемых на этапе 123. Длительность tmax - tmin целевого интервала для этапа 129 определения многокритериальной целевой функции, очевидно, влияет на оценочные значения и время вычислений и должна выбираться адекватно.

В предпочтительном варианте осуществления скалярные величины имеют следующие значения: α=1, β=0,6 и χ=0,4, которые позволяют получить достаточно точную универсальную оценку сердечных сокращений. Целевой интервал (т.е. tmax-tmin) должен составлять около 1,2 секунды, чтобы сократить время вычислений.

В альтернативном варианте осуществления, если способ следует использовать для пациентов с аритмией, например, скалярные величины могут принимать следующие значения: α=1, β=0 (так что долгосрочный прогноз практически отсутствует или игнорируется) и χ=0,2. Целевой интервал может быть задан чуть более длинным, чем в общем случае (т.е. tmax-tmin=2), чтобы верно идентифицировать аритмии.

Расширение многокритериальной целевой функции αλ(t)+βμ(t)+χσ(t) путем добавления к этой функции дополнительных оценочных функций, умноженных на соответствующие скалярные величины, - удобный способ использовать дополнительные источники информации.

Кроме того, способ может быть использован для оценки сердечных сокращений практически в реальном масштабе времени, лишь с задержкой на интервал времени, необходимый для выполнения долгосрочного прогноза и вычислений.

Вышеописанный способ позволяет выделить из баллистокардиографических сигналов интервалы между чередующимися сердечными сокращениями. Обеспечение оценочного значения сердечного ритма с учетом поочередных сердечных сокращений с использованием описанного способа может заменить стандартные ECG-устройства в различных областях применения, например при лечении сердечной недостаточности, выявлении аритмии, диагностики и лечении мерцания предсердий. Кроме того, анализ с учетом чередования ударов пульса позволяет точно вычислить вариабильность сердечного ритма, что необходимо для анализа состояния сна и стрессовых нагрузок. Потенциальные возможности выявления аритмии приведут к более широкому использованию BCG-подхода в медицинских учреждениях для мониторинга неострых состояний. Это представляет особый интерес, поскольку число коек в блоках интенсивной терапии ограничено, и желательно найти простое и экономичное решение по мониторингу для общей больничной палаты.

Наконец, можно прогнозировать обморочные состояния, связанные с деятельностью сердца, а значит, их избежать, используя представленный алгоритм, возможно в сочетании с другими схемами датчиков.

Изобретение может использоваться в домах престарелых с медицинским обслуживанием, больницах и в условиях домашнего наблюдения. Во всех случаях общее преимущество BCG над ECG заключается в ненавязчивом мониторинге пациентов без необходимости крепления к пациенту электродов и подобных элементов.

Хотя изобретение было описано в отношении способа или алгоритма, следует понимать, что изобретение может быть реализовано в виде BCG-системы (т.е. компьютерного устройства в сочетании с устройством измерения BCG-сигналов) или в виде автономной компьютерной системы или программы. Следует понимать, что BCG-система может предоставить баллистокардиографический сигнал в аналоговой или цифровой форме в устройство согласно настоящему изобретению, а устройство согласно настоящему изобретению соответственно может быть выполнено с возможностью получения этого сигнала. Например, BCG-система может подавать баллистокардиографический сигнал на устройство в аналоговой форме, а устройство может содержать фильтр для устранения ступенчатых искажений и аналого-цифровой преобразователь для передачи цифрового представления баллистокардиографического сигнала в программируемый процессор цифровых сигналов в этом устройстве. В альтернативном варианте BCG-система может быть оснащена аналого-цифровым преобразователем, так что баллистокардиографический сигнал поступает в устройство (конкретно в процессор цифровых сигналов в этом устройстве) в цифровой форме. Устройство может принимать баллистокардиографический сигнал с использованием любого соответствующего средства, например посредством проводного или беспроводного соединения с BCG-системой.

Хотя изобретение было подробно проиллюстрировано и описано на чертежах и в предшествующем описании, такое представление и описание следует рассматривать в качестве иллюстративного или служащего в качестве примера, но не ограничивающего; изобретение не ограничивается раскрытыми вариантами осуществления.

Специалисты в данной области техники смогут увидеть и реализовать раскрытые варианты осуществления в измененном виде, воплощая на практике заявленное изобретение, изучив чертежи, описание и прилагаемую формулу изобретения. В формуле изобретения термин "содержащий" не исключает присутствия других элементов или этапов, а единственное число не исключает множественного числа.

Единственный процессор или другой блок может выполнять функции нескольких элементов, приведенных в формуле изобретения. Тот факт, что определенные признаки упоминаются во взаимно отличных зависимых пунктах формулы изобретения, не означает, что сочетание этих признаков нельзя использовать с выгодой. Ни один из номеров ссылки в формуле изобретения не следует рассматривать как ограничивающий объем изобретения. Программное обеспечение может храниться/обеспечиваться на соответствующем носителе информации, таком как оптический носитель информации или твердотельный носитель информации, поставляемом совместно с другими аппаратными средствами или в качестве их части, но может также обеспечиваться в иных формах, например через Интернет или посредством других проводных или беспроводных телекоммуникационных систем.

1. Способ анализа баллистокардиографического сигнала для определения частоты сердечных сокращений, содержащий этапы, на которых:
определяют оценочное начальное временное значение для первого сокращения сердца в баллистокардиографическом сигнале;
итерационно вычисляют оценочные значения для последующих сокращений сердца в баллистокардиографическом сигнале, используя оценочное начальное временное значение;
при этом каждая итерация на этапе вычисления содержит оценку целевой функции, которая содержит взвешенную сумму множества оценочных функций;
причем каждый итерационный этап при вычислении оценочных значений для последовательных сердечных сокращений в баллистокардиографическом сигнале ограничен целевым интервалом после оценочного временного значения, найденного на предшествующем итерационном этапе вычисления.

2. Способ по п.1, в котором оценка целевой функции содержит идентификацию максимального значения целевой функции, при этом максимальное значение целевой функции соответствует сердечному сокращению.

3. Способ по п.1, в котором одна из оценочных функций отражает наличие высокочастотных составляющих в баллистокардиографическом сигнале.

4. Способ по п.2, в котором одна из оценочных функций отражает наличие высокочастотных составляющих в баллистокардиографическом сигнале.

5. Способ по п.3, в котором упомянутую оценочную функцию, отражающую наличие высокочастотных составляющих в баллистокардиографическом сигнале, определяют с помощью этапов, на которых:
фильтруют баллистокардиографический сигнал для выделения высокочастотных составляющих;
возводят в квадрат отфильтрованный баллистокардиографический сигнал;
фильтруют возведенный в квадрат и отфильтрованный баллистокардиографический сигнал для получения конечного сигнала.

6. Способ по п.4, в котором упомянутую оценочную функцию, отражающую наличие высокочастотных составляющих в баллистокардиографическом сигнале, определяют с помощью этапов, на которых:
фильтруют баллистокардиографический сигнал для выделения высокочастотных составляющих;
возводят в квадрат отфильтрованный баллистокардиографический сигнал;
фильтруют возведенный в квадрат и отфильтрованный баллистокардиографический сигнал для получения конечного сигнала.

7. Способ по п.5 или 6, в котором этап определения оценочного начального временного значения для первого сердечного сокращения в баллистокардиографическом сигнале содержит этапы, на которых:
выбирают точку максимума упомянутой оценочной функции, которая отражает наличие высокочастотных составляющих в баллистокардиографическом сигнале, при этом точка максимума выбрана из интервала в начальной части баллистокардиографического сигнала.

8. Способ по любому из пп.1-6, в котором этап определения оценочного начального временного значения для первого сердечного сокращения в баллистокардиографическом сигнале содержит этапы, на которых:
оценивают целевую функцию для оценочного начального значения и оценочного возможного временного значения второго сердечного сокращения в баллистокардиографическом сигнале; и
выбирают оценочное начальное временное значение и оценочное возможное временное значение в качестве значений, при которых целевая функция достигает максимума для первых двух сердечных сокращений.

9. Способ по любому из пп.1-6, в котором этап определения оценочного начального временного значения для первого сердечного сокращения в баллистокардиографическом сигнале содержит этап, на котором:
выбирают значение времени в начальной части баллистокардиографического сигнала в качестве оценочного начального временного значения.

10. Способ по любому из пп.1-6, в котором после выполнения множества итераций этапа вычисления способ дополнительно содержит уточнение оценочного начального временного значения путем:
выполнения дополнительных этапов вычисления, при этом каждый целевой интервал для проведения уточнения ограничен целевым интервалом в малой близлежащей области оценочного временного значения, найденного на предшествующем итерационном этапе вычисления.

11. Способ по любому из пп.1-6, в котором одна из оценочных функций отражает долгосрочное прогнозирование частоты сердечных сокращений.

12. Способ по п.11, в котором упомянутую оценочную функцию, отражающую долгосрочное прогнозирование частоты сердечных сокращений, определяют с помощью этапа, на котором:
оценивают спектрограмму баллистокардиографического сигнала во множестве временных точек, при этом временные точки зависят от оценочного временного значения, найденного на предшествующем итерационном этапе вычисления.

13. Способ по п.11, в котором упомянутую оценочную функцию, отражающую долгосрочное прогнозирование частоты сердечных сокращений, определяют с помощью этапа, на котором:
оценивают автокорреляционную функцию баллистокардиографического сигнала.

14. Способ по любому из пп.1-6, в котором одна из оценочных функций отражает информацию о стохастических характеристиках интервалов сердечных сокращений.

15. Способ по п.14, в котором упомянутую оценочную функцию, отражающую информацию о стохастических характеристиках интервалов сердечных сокращений, определяют с помощью этапа, на котором:
рассчитывают вероятностное распределение, отражающее вероятностное распределение длительности интервалов сердечных сокращений.

16. Способ по п.14, в котором упомянутую оценочную функцию, отражающую информацию о стохастических характеристиках интервалов сердечных сокращений, определяют с помощью этапа, на котором:
оценивают синхронизацию обнаруженных ранее сердечных сокращений.

17. Способ по любому из пп.1-6, дополнительно содержащий этап, на котором:
используют оценочные временные значения последующих сердечных сокращений и баллистокардиографический сигнал для уточнения идентификации сердечных сокращений.

18. Способ по п.17, в котором этап использования оценочных временных значений содержит этапы, на которых:
применяют фильтр верхних частот к баллистокардиографическому сигналу;
идентифицируют максимальные значения в отфильтрованном баллистокардиографическом сигнале, связанные с оценочными временными значениями; и
идентифицируют сердечные сокращения в виде волнового сигнала низкой частоты, но с высокой амплитудой, вслед за идентифицированными максимальными значениями баллистокардиографического сигнала.

19. Способ по п.17, в котором этап использования оценочных временных значений содержит этап, на котором:
используют способ обнаружения периодичности во временном ряду для идентификации кардиоинтервала между двумя последовательными сердечными сокращениями.

20. Машиночитаемый носитель, содержащий управляющую компьютерную программу, которая при выполнении на компьютере или процессоре выполнена с возможностью реализации этапов способа по пп.1-19.

21. Устройство для использования совместно с блоком измерения баллистокардиографического сигнала у пациента, при этом устройство содержит:
средство получения баллистокардиографического сигнала от блока; и
средство обработки информации для реализации способа, охарактеризованного в каждом из пп.1-19, применительно к полученному баллистокардиографическому сигналу.



 

Похожие патенты:

Группа изобретений относится к медицине. Способ получения биологической информации осуществляют с помощью системы получения биологической информации.

Группа изобретений относится к медицине. Респираторный монитор содержит первый датчик, выполненный с возможностью генерирования сигнала мониторинга респираторно-связанных движений, указывающего на респираторно-связанные движения; второй датчик, выполненный с возможностью генерирования сигнала мониторинга шумов, указывающего на респираторно-связанные шумы; а также синтезатор сигналов, выполненный с возможностью синтеза сигнала респираторного монитора на основе сигнала мониторинга респираторно-связанных движений и сигнала мониторинга респираторно-связанных шумов.

Изобретение относится к медицине, в частности к функциональной диагностике и диагностике в онкологии, и может быть использовано для скрининг-диагностики рака простаты.

Изобретение относится к медицинской технике, а именно к устройствам и способам определения равновесия при движении тела. Устройство содержит блок получения сигнала, связанного с движением тела, блок определения равновесия при движении тела и блок вывода результата.

Изобретение относится к медицине. Портативное устройство для бесконтактной выборочной проверки жизненных показателей пациента содержит: датчик расстояния для последовательного обнаружения изменений расстояния во времени относительно грудной клетки пациента, калькулятор частоты дыхания для определения дыхательной активности на основе обнаруженных изменений расстояния во времени.

Изобретение относится к медицине, судебной медицине и предназначено для идентификации личности неопознанных трупов и их фрагментов. Изобретение также может быть использовано при необходимости прижизненной идентификации человека в случае изменения внешности.

Группа изобретений относится к областям биометрии, электроники, медицины и психологии и может быть использовано для получения психофизиологической информации о живых объектах и измерения, контроля и коррекции психофизиологического состояния человека.

Изобретение относится к способу и системе, обеспечивающим определение возраста пользователя в сети по данным большого объема. Техническим результатом является обеспечение возможности точной фильтрации пользователей сети по возрасту.

Группа изобретений относится к области медицины, а именно к травматологии и ортопедии, педиатрии, нейрохирургии, а также невропатологии, и предназначена для осуществления коррекции и профилактики функциональных нарушений осанки путем информирования пользователя о недопустимом изменении положения его позвоночника в сагиттальной или фронтальной плоскости с целью выработать у пользователя навык правильной осанки, тем самым способствуют лечению ортопедических заболеваний и мобильных деформаций позвоночника.

Изобретение относится к медицинской технике, а именно к устройствам для исследования движения тела человека. В первом варианте устройство выполнено с возможностью установки на голове пользователя в области его височной и/или жевательной мускулатуры и включает датчик Холла, по меньшей мере, один постоянный магнит, установленные с возможностью взаимного смещения в упруго деформируемом корпусе, и блок управления и обработки информации.

Изобретение относится к средствам определения движения тела. Устройство содержит средство определения ускорения и вычислительное средство для вычисления движения тела на основании данных ускорения, участок закрепления/раскрепления для закрепления рабочей части на основном блоке устройства или раскрепления от него, причем вычислительное средство выполнено с возможностью выполнения процедуры определения закрепления/раскрепления на основании изменения ускорения, при закреплении закрепляемой рабочей части на участке закрепления/раскрепления или при раскреплении от него, и выполнения вычисления движения тела на основании определенного закрепления/раскрепления, при переключении в режим, соответствующий состоянию после закрепления/раскрепления. Во втором варианте выполнения устройства дополнительно имеется запоминающее средство для хранения данных и дисплейное средство для отображения результата вычисления, а также множество участков закрепления/раскрепления в зависимости от типа закрепляемой рабочей части, причем вычислительное средство выполнено также с возможностью определения типа закрепляемой рабочей части на основании изменения ускорения, и выполнения вычисления движения тела на основании определенного закрепления/раскрепления и типа, при переключении в режим, соответствующий типу закрепленной закрепляемой рабочей части. В третьем варианте участок закрепления/раскрепления содержит направляющую для сдвига или поворота закрепляемой рабочей части, когда закрепляемая рабочая часть закреплена на основном блоке устройства или раскреплена от него, и ударный участок, с которым одна часть закрепляемой рабочей части сталкивается при закреплении и раскреплении закрепляемой рабочей части вдоль направляющей. В четвертом варианте устройство имеет корпус для вмещения средства определения ускорения и вычислительного средства, при этом корпус содержит ударную рабочую часть, подлежащую удару при столкновении, а вычислительное средство сконфигурировано с возможностью определения данных ускорения при ударе, нанесенном по ударной рабочей части, и определения информации о нанесенном ударе из данных ускорения. В способе определения движения тела определяют изменение ускорения, которое возникает при закреплении рабочей части на основном блоке устройства или при раскреплении от него, на основании данных ускорения, и движение тела на основании данных ускорения после закрепления или раскрепления. Использование изобретения позволяет повысить точность измерения. 5 н. и 4 з.п. ф-лы, 16 ил.

Изобретение относится к области эргономики, психологии труда, медицине и может быть использовано для диагностики функционального состояния человека, а именно к исследованию и оценки усталости глаз пользователя компьютера, и искривления позвоночника пользователя ПК. Предварительно с помощью программного средства фиксируют нормальное вертикально ориентированное положение позвоночника пользователя, отклонение положения позвоночника от нормального его положения на недопустимую величину. Производят коррекцию положения позвоночника путем его возврата в нормальное вертикально ориентированное положение. Нормальное вертикально ориентированное положение позвоночника пользователя фиксируют встроенной фронтальной видеокамерой с последующей передачей видеоизображения на программный сервис компьютера. С помощью программы и видеокамеры определяют проекцию расстояния между центрами глаз, зафиксированное на видеоизображении при комфортном расстоянии между пользователем и экраном компьютера, которое составляет не менее 2-х диагоналей экрана. Путем съемки видеокамерой регистрируют динамику изменения проекции расстояния между центрами глаз в кадре во времени. Если расстояние между пользователем и экраном компьютера составляет менее 2-х диагоналей экрана, выводят предупреждение на экран. В качестве предупреждения выводят наложение на изображение экрана визуальных эффектов и/или звуковых и/или вибросигналов. При этом чем меньше расстояние между пользователем и экраном, тем сильнее воздействие визуальных эффектов в виде размытия экрана и/или звуковых и/или вибросигналов. Для самоконтроля пользователя при фиксации изменения расстояния между центрами глаз в проекции кадра в качестве визуального эффекта используют затенение экрана и/или изображение анимационных персонажей. После 30-минутной работы программы выводят на экран сообщение о необходимости перерыва и отдыха для глаз и выводят игровые расслабляющие упражнения для глаз или рекомендацию о перерыве на 5 минут, после чего визуальный эффект экрана и дополнительные сигналы выключают. Для экономии программно-аппаратного ресурса используют алгоритм разложения видеоряда на серию кадров с последующим анализом лишь каждого пятого кадра. Способ позволяет предотвратить усталость глаз и одновременного искривления позвоночника пользователя за счет поддержания безопасного расстояния от экрана до глаз человека. 3 з.п. ф-лы, 1 ил., 1 пр.

Изобретение относится к медицинской технике. Устройство для определения функционального состояния опорно-двигательного аппарата содержит регистратор параметров опорно-двигательного аппарата. Регистратор включает датчики веса и поддерживаемую опорными элементами опорную пластину для стоп с установленным под пластиной датчиком изображения отпечатка подошвенной поверхности стоп, подключенным к компьютеру. Опорная пластина выполнена из оптически прозрачного материала, а опорные элементы выполнены в виде стоек, снабженных датчиками веса. Информационные выходы датчиков веса связаны с компьютером, выполненным с возможностью регистрации и одновременного отображения в одной системе координат изображения отпечатка подошвенной поверхности стоп и данных о положении центра давления на каждой из стоп и общего центра давления стоп. Применение изобретения позволит расширить функциональные возможности устройства, сократить временные затраты при проведении исследования, повысить точность определения положения центров давления по отношению к положению стоп за счет обеспечения возможности оперативного комплексирования результатов компьютерной плантографии и стабилометрии. 4 з.п. ф-лы, 5 ил.
Изобретение относится к медицине, педиатрии, перинатологии, детской неврологии и может быть использован для ранней диагностики нарушения движений у детей с перинатальным поражением нервной системы. В качестве признака раннего проявления нарушений движений диагностируют гиперкинезы мышц глазных яблок путем сочетанной оценки наличия симптома псевдо-Грефе и повышения уровня непрямого билирубина более 275 ммоль/л в течение первого месяца жизни ребенка. Способ обеспечивает раннюю диагностику нарушений движения у грудных детей с перинатальной патологией нервной системы, раннее распознавание симптомов поражения подкорковых ядер головного мозга за счет выявления наиболее патогномоничных признаков раннего начала болезни. 1 пр.

Изобретение относится к области медицины, а именно к судебной медицине, и может быть использовано для определения давности повреждений на трупе. Определяют коэффициент теплопроводности мягких тканей области повреждений, концентрации алкоголя в крови трупа, расчет. Расчет давности повреждений выполняется по формулам, различающимся для случаев наличия либо отсутствия алкоголя в крови трупа. Способ позволяет повысить точность определения давности повреждений на трупе, в том числе повреждений, расположенных на слизистой оболочке рта. 2 пр.

Изобретение относится к области медицины, а именно к судебной медицине, и может быть использовано для определения давности повреждений у живых лиц. Измеряют температуру области повреждения и неповрежденного участка и температуру окружающей среды. Давность повреждений у живых лиц определяют по соответствующей формуле. Способ позволяет повысить точность определения давности повреждений у живых лиц, в том числе повреждений, расположенных на слизистой оболочке рта. 2 пр.

Изобретение относится к медицине, а именно к судебной медицине, и может быть использовано для определения давности локального повреждения мягких тканей по температуре области мягких тканей. Способ заключается в проведении дистанционной термографии, с помощью которой определяют максимальную и минимальную температуры области локального повреждения мягких тканей и интактного участка. Кроме того, определяют площадь локального повреждения мягких тканей, однородную по температуре, и коэффициент излучения. На основании полученных данных рассчитывают давность локального повреждения мягких тканей. Способ позволяет повысить точность и расширить функциональные возможности определения давности локального повреждения мягких тканей (кровоподтека, внутримышечной инъекции), за счет максимально полного учета теплового состояния тела пациента в области локального повреждения мягких тканей. 2 пр.

Изобретение относится к медицине, а именно к психофизиологии, биомедицинским и психологическим исследованиям. Способ исследования двигательных и когнитивных функций человека реализуют в виде стабилометрического исследования, где испытуемый, управляя позой, выполняет инструкцию по удержанию собственного центра давления на стабилометрическую платформу в заданной зоне. Исследование включает наличие двух последовательных равных по длительности фаз, в первой из которых отсутствует биологическая обратная связь (БОС) по опорной реакции с помощью искусственного контура БОС, а во второй - присутствует. По итогам теста проводят сопоставление полученных в двух фазах количественных параметров, связанных с управлением позой, и количественных параметров, отражающих степень выполнения инструкции. По соотношению стабилометрических показателей судят о функциональном состоянии испытуемого. Способ обеспечивает объективную оценку состояния человека за счет количественного исследования способности выполнять заданное инструкцией управление позой тела и вниманием - сочетание двигательной и когнитивной задач в едином тесте, позволяет объективно оценить влияние когнитивной составляющей на управление позой при включении БОС по опорной реакции. 3 ил.

Изобретение относится средствам для бесконтактного мониторинга дыхания пациента. Способ обнаружения изменения от выдоха до вдоха пациента или наоборот включающий этапы излучения электромагнитного сигнала в сторону пациента и приема отраженного от пациента сигнала, преобразования отраженного сигнала с получением первого сигнала, сдвига по фазе отраженного электромагнитного сигнала и преобразования его с получением второго сигнала, обнаружение с помощью вычислительного блока одновременных первых переходов через ноль во временной производной первого сигнала и во временной производной второго сигнала, одновременных вторых переходов через ноль во временной производной первого сигнала и во временной производной второго сигнала, и одновременных третьих переходов через ноль во временной производной первого сигнала и во временной производной второго сигнала, определения первого и второго векторов и вычисления их скалярного произведения в качестве индикаторного значения для изменения от выдоха до вдоха пациента или наоборот, сравнения индикаторного значения с предварительно определенным пороговым значением и указания изменения от выдоха до вдоха пациента или наоборот, если индикаторное значение является меньшим, чем пороговое значение. Устройство для осуществления способа включает двухканальный доплеровский радиолокационный датчик и вычислительный блок. Использование изобретения позволяет повысить надежность измерения и повысить точность обнаружения частоты дыхания. 2 н. и 4 з.п. ф-лы, 6 ил.
Изобретение относится к медицине, травматологии и ортопедии, реабилитации и неврологии. Проводят оценку функционального укорочения нижней конечности, для чего при ходьбе пациента измеряют мощность шага левой и правой нижних конечностей и при асимметрии показателей до 1% относительно средней мощности шага укорочение считают компенсированным. При асимметрии от 2% до 8% - субкомпенсированным, а при асимметрии более 9% - декомпенсированным. Способ позволяет повысить точность оценки функционального укорочения нижней конечности, подобрать коррекцию, исключающую перегрузку больной конечности и дискомфорт, замедлить прогрессирование дегенеративных процессов суставов, улучшить походку, обеспечить профилактику нестабильности компонентов эндопротеза после эндопротезирования крупных суставов нижних конечностей. 2 пр.
Наверх