Из вышесказанного следует, что задача анализа НАП на ЭКС и измерения их параметров является на сегодняшний день весьма актуальной задачей. Она может быть успешно решена при помощи вероятностных методов обработки электрокардиосигналов. Рассмотрим применение этих методов на примере анализа одного из видов НАП - поздних потенциалов желудочков сердца (ППЖ).
Поздние потенциалы желудочков [1,2], которые возникают в конце комплекса QRS и первой трети сегмента ST. Их наличие говорит о повышенной (до 50%) вероятности возникновения и развития у данного пациента аритмий сердца, в том числе и опасных для жизни. Поэтому своевременное обнаружение ППЖ чрезвычайно важная задача. Ее решение позволяет еще до наступления аритмии заранее принять меры к ее недопущению.
Однако, проблема регистрации ППЖ заключается в том, что на обычной электрокардиограмме они теряются в шумах, так как их максимальная амплитуда не превышает 30 мкВ. Для повышения отношения сигнал/шум широко применяется временное усреднение ЭКС. Опорной точкой при суммировании кардиоциклов служит вершина пика R. При регистрации ППЖ наиболее распространен метод Симсона [2].
ЭКС по методу Симсона снимается с поверхности тела человека по трем биполярным ортогональным отведениям X, Y, Z. Предварительная аналоговая обработка сигнала заключается в усилении сигнала в 2000 - 5000 раз и фильтрации его для каждого отведения в отдельности. Далее ЭКС преобразуется в цифровую форму с частотой дискретизации FД=1000 Гц и усредняется во времени по вершине пика R. Затем усредненные сигналы фильтруются цифровым фильтром высоких частот с частотой среза 25 Гц. Затем по ним вычисляется амплитуда вектора ЭДС сердца, которая носит название фильтрованного QRS-комплекса (FQRS).:
, (1)
где UX, UY, UZ - усредненные и фильтрованные ЭКС X, Y, Z соответственно.
Анализ ППЖ представляет собой вычисление трех амплитудно-временных параметров FQRS:
- Длительность фильтрованного QRS комплекса - DFQRS.
- Длительность низкоамплитудной (ниже 40 мкВ) части заднего фронта QRS комплекса - LAS.
- Среднеквадратичная амплитуда за последние 40 мс FQRS - RMS.
Критерием наличия ППЖ является превышение хотя бы двумя из трех вышеперечисленных параметров следующих порогов: DFQRS>120 мс; LAS>39 мс; RMS<25 мкВ [2]. В противном случае делается вывод об отсутствии ППЖ.
Значения параметров метода Симсона зависят от временного положения двух характерных точек FQRS - его начала ТН и конца ТК. Поэтому эффективность алгоритма анализа наличия ППЖ на ЭКС определяется величиной ошибки вычисления временного положения ТН и ТК - τош. Эта ошибка значительно уменьшается в результате применения вероятностных методов при определении точек.
Для нахождения точки ТК ищется участок DТС длиной 5 мс, среднее значение амплитуды которого mС. превышает порог:
. (2)
где, mШ и σШ - среднее значение и С.К.О. амплитуды участка шума на FQRS.
Найденный ΔТС содержит внутри себя точку ТК. Далее участок ΔТС делится пополам и точка половинного деления принимается за точку конца фильтрованного комплекса QRS - ТК. При этом средняя ошибка τош. составляет 1,2 мс [2], а максимальная ошибка - 2,5 мс. Этот способ (деление ΔТС пополам) хорошо работает при уровне шумов на FQRS не более 1,5 мкВ.
При большем уровне шумов из-за малой длины ΔТС он может быть найден неверно и истинная точка ТК окажется вне ΔТС; поэтому длину ΔТС надо увеличить. Но это приводит существенному возрастанию τош. В нашем исследовании (346 пациентов) при уровне шумов до 4 мкВ длина DТС составила 30 мс; только тогда ТК попадала внутрь ΔТС во всех случаях. При этом τош составила 7,2 мс в среднем, а максимальная ошибка - 15 мс. Это недопустимо т.к. примерно в половине случаев привело к неверному решению о наличии или отсутствии ППЖ на ЭКС.
Было предложено вместо деления ΔТС пополам для нахождения точки ТК применить статистический алгоритм проверки гипотез о временном положении точки ТК внутри DТС. Отсчеты FQRS следуют через 1 мс и ΔТС содержит 30 отсчетов. Поэтому были выдвинуты 30 гипотез о принадлежности точки ТК каждому из отсчетов. Для каждой из гипотез вычислялась целевая функция и по максимуму целевой функции принималось решение о временном положении точки конца комплекса QRS.
Построим целевую функцию на основе критерия максимального правдоподобия [3]:
, (3)
где Hi - гипотеза; w0 - плотность распределения вероятностей сигнала (заднего фронта FQRS); w1 - плотность распределения вероятностей шума (сегмента ST FQRS). Будем принимать как правильное решение ту гипотезу, для которой функция (3) окажется максимальна:
, (4)
В работе [4] были найдены плотности распределения вероятностей w0 и w1 с учетом различного характера изменений заднего фронта FQRS у разных пациентов (было выделено 3 класса) и сегмента ST (выделено 2 класса). Поэтому общее выражение (3) конкретизируется в нашем исследовании:
(5)
qjш и qjс - весовые коэффициенты выделенных классов для областей сегмента ST и заднего фронта FQRS соответственно.
При определении временного положения точки Tk с использованием алгоритма (5) были получены следующие результаты.
Ошибка определения Tk составила в среднем 1,1 мс, что вполне приемлемо с практической точки зрения для автоматизированного диагностического комплекса. Эта ошибка примерно равна ошибке стандартного метода Симсона, однако в нашем исследовании допускался гораздо больший уровень шумов (до 4 мкВ).
Всего с использованием статистического алгоритма (5) было исследовано 173 пациента. При этом в 150 случаях наблюдалось полное совпадение в оценке точки алгоритмом и врачами-экспертами. Ошибочных определений ППЖ наблюдалось всего 8 из 173; в том числе 5 случаев пропуска ППЖ и 3 случая ложного обнаружения ППЖ, т.е. всего было 4.6% неверных решений относительно ППЖ при применении статистического алгоритма.
Таким образом, в результате применения статистического алгоритма (3) существенно повысилась достоверность проводимого автоматического анализа ППЖ. Это, в свою очередь позволило улучшить качество и расширить возможности ранней диагностики различных аритмий сердца. Данный алгоритм ценен еще и тем, что он достаточно просто реализуется на компьютере и без существенных изменений может с успехом применяться для более точного анализа других видов НАП, например, таких как поздние потенциалы предсердий или вызванные потенциалы головного мозга.
ЛИТЕРАТУРА
- Иванов Г.Г. и др. Электрокардиография высокого разрешения: некоторые итоги 4-летних исследований Кардиология. 1994. № 5. С. 22 - 25.
- Simson M. B. Use of signal in the terminal QRS komplex to identify patients with ventricular tachycardia after myocardial infarction. Circulation. 1981. Vol. 64.P. 235 - 242.
- Левин Б. Р. Теоретические основы статистической радиотехники Т2. -М.: Советское радио, 1975. 390 с
- Седов С.С., Цибулькин Н.А., Щербакова Т.Ф. Плотности распределения вероятностей амплитуды вектора ЭДС сердца и их применение при диагностике поздних потенциалов желудочков сердца. Биомедицинская радиоэлектроника. 1998. №4. С. 34-39.