Способ обработки электрокардиосигнала для персональных носимых кардиомониторов

Изобретение относится к области медицины, в частности к методам автоматического анализа биомедицинской информации, и предназначено для предварительной обработки электрокардиосигнала (ЭКС), поиска R-пиков, определения R-Rинтервалов и других информативных параметров электрокардиосигнала, например показателей вариабельности сердечного ритма в персональных носимых кардиомониторах.

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

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

Известны способы предварительной обработки ЭКС, наиболее совершенным из которых считается метод Пана - Томпкинса, которые основаны на предварительной фильтрации ЭКС, взятии первой или второй производной от ЭКС, ее возведении в квадрат или куб с последующей фильтровой обработкой [Рангайян P.M. Анализ биомедицинских сигналов, пер. с англ. А.Н. Калиниченко под ред. Немирко. М. Физматлит. 2007; Компьютерные методы автоматического анализа ЭКГ в системах кардиологического наблюдения. Калиниченко А.Н., докт. диссерт. 12Д5М 2003, СПб, 2008, заявка RU №2014127196, МКИ А61В 5/0452, опубл. 10.02.2016]. Недостаток этих способов заключается в том, что результирующий сигнал не имеет никакого сходства с исходным ЭКС, т.к. относительно медленные зубцы Р и Т подавляются, а выходной сигнал приобретает дополнительный шумовой характер, для устранения которого требуется дополнительная операция сглаживания, что требует привлечения дополнительных вычислительных ресурсов, увеличения длительности обработки и, как следствие, приводит к увеличению энергопотребления. Кроме того, известные способы характеризуются относительно плохими результатами обнаружения проблемных сигналов, в частности, содержащих дрейф изолинии.

Известны также способы предварительной обработки ЭКС (пат. RU 2345709, МПК А61В 5/0432, опубл. 10.02.2009, патент RU 2373945, МПК А61В 5/0402, опубл. 27.11.2009, заявка RU 2009 105 130, МПК А61В 5/0402, опубл. 27.08.2010), использующие встроенные средства фильтрации сигнала ЭКГ, расположенные между датчиками ЭКГ и блоком обработки сигналов, что, однако, искажает исходную форму ЭКС, а следовательно, снижает достоверность диагностики ЭКС. Кроме того, это создает дополнительную групповую задержку, которая приводит к паразитному временному сдвигу в отсчетах, для устранения которого приходится использовать дополнительную операцию сглаживания с использованием адаптивных фильтров, что также приводит к необходимости увеличения длительности обработки и вычислительной мощности. При этом выходной QRS-комплекс также не имеет сходства с исходным сигналом, что приводит к снижению достоверности выявления информативных параметров ЭКС.

Известны способы, использующие вейлет-преобразования, нейронные сети или спектральную обработку функций Уолша [Сравнительный анализ алгоритмов обнаружения R-зубца ЭКС, А.В. Рослякова, П.Г. Чупраков, Экспериментальная медицина и клиническая диагностика, Вятский медицинский вестник, №2, 2012, 29-35; Патент РФ 2242164, МПК А61В 5/0402 опубл. 20.12.2004]. Однако для качественной обработки ЭКС известные способы требуют очень высоких вычислительных мощностей и, следовательно, высокого энергопотребления. Поэтому в персональных носимых мониторах для длительного мониторинга они не могут быть использованы.

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

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

Сущность изобретения заключается в том, что с датчиков ЭКГ снимают электрокардиосигналы, усиливают их, проводят первичную аналоговую обработку, затем аналоговый сигнал дискретизируют по амплитуде и по времени для перевода в цифровой сигнал, корректируют дрейф изолинии путем определения скользящего среднего на прямоугольном временном окне длительностью N1 в 30 дискретных отсчетов с последующим его вычитанием из входного сигнала, затем проводят адаптивную нормализацию сигналов R-пиков, для чего на каждом шаге дискретизации выбирают скользящий максимум Amax(i) на прямоугольном окне N₂ в 300 дискретных отсчетов и умножают значения каждого текущего отсчета окна N₂ на коэффициент Gi=1/Amax(i), далее устраняют шумовые и высокочастотные наводки путем оценки для каждого текущего отсчета i разности ΔS_i между площадью сигналов S_N3 в рамках скользящего прямоугольного окна N₃ в 15 дискретных отсчетов симметрично по семь отсчетов вправо и влево относительно текущего отсчета и ее средним значением S_{N1 ср.} на скользящем окне N₁, и если разность ΔSi составляет величину более двух - оставляют исходную форму сигнала, а менее двух - заменяют на среднее значение сигнала в скользящем окне N₁, затем осуществляют локализацию положения R-пиков пороговым детектированием по уровню 0,5, выбором пиков и их расположения на временном окне из каждых пяти рядом расположенных отсчетов, а также удалением нетипичных кардиосигналов, далее проводят извлечение информации, относящейся к сердечным сокращениям пользователя.

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

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

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

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

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

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

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

Сущность изобретения поясняется чертежами, где на фиг. 1 дана блок-схема для первичной обработки электрокардиосигнала (ЭКС); на фиг. 2 показан дрейф изолинии; на фиг. 3 - электрокардиосигнал после стабилизации изолинии, на фиг. 4 - электрокардиосигнал до нормализации сигналов R-пиков, на фиг. 5 - ЭКС после адаптивной нормализации сигналов R-пиков; на фиг. 6 - ЭКС до сглаживания по площади, на фиг. 7 - ЭКС после сглаживания по площади.

Блок-схема устройства для реализации способа содержит два датчика ЭКГ 1, 2, подключенных к усилителю 3, выход которого связан с последовательно соединенными блоком обработки аналоговых сигналов 4, аналого-цифровым преобразователем 5, блоком цифровой обработки ЭКС 6.

Способ осуществляется следующим образом.

Размещают датчики ЭКГ 1, 2 в области сердца, например во 2-м или 5-м модифицированных отведениях. С датчиков ЭКС 1, 2 снимают электрокардиосигналы, усиливают их в усилителе 3, проводят первичную аналоговую обработку, в частности попарную дифференциальную обработку в блоке обработки аналоговых сигналов 4, затем аналоговый сигнал в аналого-цифровом преобразователе 5 дискретизируют по амплитуде и по времени для перевода в цифровой сигнал.

Цифровой сигнал подают в блок цифровой обработки ЭКС 6. Последующую цифровую обработку ЭКС проводят в прямоугольных окнах. При этом длительность окон выбирают в зависимости от этапа обработки ЭКС.

Предварительно корректируют дрейф изолинии (фиг. 2) оцифрованного электрокардиосигнала. Для этого определяют скользящее среднее на прямоугольном окне N₁ длительностью 30 дискретных отсчетов и его вычитают из входного сигнала, устраняя тем самым пьедестал дрейфа.

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

Далее проводят адаптивную нормализацию сигналов R-пиков (фиг. 4), т.е. приводят их к уровню нормированной единицы (фиг. 5). Для этого определяют скользящий максимум A_max(i) на прямоугольном окне N₂ в 300 дискретных отсчетов (размер окна должен быть менее длительности следования R-пиков), и умножают после этого значения текущих отсчетов на коэффициенты усиления G_i=1/A_max(i). Нормированные сигналы пприведены на фиг. 5.

Далее проводят сглаживание по площади нормализованного сигнала для исключения шумоподобных составляющих электрокардиосигнала, обусловленных наводками и мимикрией мышц (фиг. 6). Для этого определяют площадь S_N3 сигналов для каждого текущего отсчета i в рамках окна N₃ по семь отсчетов вправо и влево относительно текущего отсчета и среднюю площадь S_{N1 ср.} этого сигнала в рамках окна N₁, т.е. в 30 дискретных отсчетов. Далее вычисляют абсолютную разность этих площадей ΔS_i. Если ΔS_i составляет более двух - оставляют исходную форму сигнала, а менее двух - этот сигнал заменяют на среднее значение сигнала в скользящем окне N₁. Таким образом, большие выбросы сигнала, такие как R-пики, на участке плюс минус семь точек остаются, а остальные шумоподобные помеховые сигналы приводятся к среднему по 30 точкам, т.е. фактически устраняются. Сигнал после сглаживания приведен на фиг. 7.

Далее осуществляют локализацию положения R-пиков. Для этого проводят пороговое детектирование по уровню 0,5, затем определяют максимальное значение из каждых пяти рядом расположенных отсчетов и фиксируют его расположение на временном окне. Найденные таким образом максимумы сигнала интерпретируют как R-пики. После этого удаляют нетипичные R-циклы ЭКС и проводят извлечение информации, относящейся к сердечным сокращениям пользователя.

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