CC BY
49
УПРАВЛЕНИЕ В МЕДИЦИНЕ И БИОЛОГИИ X
УДК 621.317.08; 681.2.001.5; 615,47
МЕТОДИКА ОБНАРУЖЕНИЯ МИКРОПОТЕНЦИАЛОВ ЭКГ
О. О. Жаринов,
ассистент А. П. Шепета,
д-р техн. наук, профессор
Санкт-Петербургский государственный университет аэрокосмического приборостроения (СПбГУАП)
Исследование тонкой структуры ЭКГ считается перспективным направлением повышения диагностических возможностей электрокардиографии. Экспериментально доказано [ 1 ], что ряд электрофи-зиологических нарушений в миокарде сопровождается появлением низкоамплитудных «зубцов», часто недетерминированной структуры и нестабильными во временихарактеристиками. Возможности регистрации микропотенциалов на обычной ЭКГ принципиально ограничены наличием множества помех, что предопределяет использование методов статистической теории обработки сигналов.
Постановка задачи. Использование системы из 3 ортогональных отведений позволяет рассматривать регистрируемый процесс как векторный
Л [А] = {г|<х>[^] Л<у>[/<] Л<г>[^]} ■ Как известно, ЭКГ является квазипериодическим сигналом, состоящим из отдельных повторяющихся кардиоциклов. Обобщенная математическая модель процесса, регистрируемого на протяжении некоторого £ -го кардиоцикла, может быть представлена в виде
ЛД/с] = 9^х(5[/с] + 5^[/с]) + ЙД/<],
1 = \1] /с = ГЙ, (1)
где информационным сигналом является низкоамплитудный компонент 5цД/с], по-разному (нерегулярно) проявляющийся в Ь кардиоциклах записи ЭКГ; все остальные составляющие: I[/с] — высокоамплитудный, повторяющийся, сигнал «обычной» ЭКГ (комплекс Р-ОЯ^-Т, регулярная активность); 9^ — матрица преобразования координат размера 3x3, описывающая механические перемещения источника сигналов; многокомпонентный случайный процесс п([к], в контексте решаемой задачи рассматриваются как помехи.
Выделение составляющих МП из смеси (1) с учетом высокой степени априорной неопределенности, широкого динамического диапазона, сложной структуры помехи и ее нестационарного характера и представляет собой сложнейшую задачу, решение которой требует многоэтапной обработки сигнала. Для приведения задачи анализа информационного сигнала к обычному виду смеси типа
г\^[к] = ^е[к] + п е[к],£ = \1, (2)
необходимо последовательно решать задачи исключения внеполосных составляющих помехи, об-
наружения кардиоциклов, оценивания параметров матрицы 91 £, оценивания сигнала регулярной активности (посредством усреднения L сигналов) и его компенсации.
Некоторые из перечисленных задач являются традиционными для автоматического анализа ЭКГ; другие, в частности необходимость оценивания матрицы 91^ [2], и собственно задача анализа выборки (2), ставятся редко, а известные методы их решения развиты недостаточно. Здесь приводится краткое описание предлагаемых подходов к реализации процедур обработки ЭКГ, имеющих ключевое значение.
Компенсация мультипликативных помех. Процессы механического перемещения источника сигналов преимущественно обусловлены дыхательной активностью, в результате чего сердце осуществляет повороты на углы вокруг ортогональных осей {фх, фу, ф2}. Оценить параметры матрицы преобразования (вращения координат), связанные нелинейной зависимостью и изменяющиеся во времени, довольно сложно. Некоторого упрощения удается достичь в случае спокойного неглубокого дыхания, когда на протяжении каждого кардиоцикла параметры можно считать постоянными и углы поворота малы, так что для любого из них sin (ф) ~ ф, Ф2 = 0, соэ(ф) = 1. В таком приближении матрица вращения представляется следующим образом:
(л \
1 -фг фу
Я= Фг 1 -Фх • (3)
,-фу Фх 1 ,
Таким образом, необходимо оценивать три параметра для каждой пары выборок, соот-
ветствующих различным кардиоциклам {т^ [/с]}^=1 и {ЧеЛк]}к=1- Поскольку информативны не абсо-
лютные значения углов, а относительные, то сигнал одной из выборок, например, ^-|-ой, принимается за начало отсчета, т. е. 9^ = I, Важно отметить, что умножение на матрицу (3) увеличивает евклидову норму вектора, ||91 х 5[/с]|| > Щк]\\. Это обстоятельство является, в общем, побочным эффектом приближения, но может оказаться полезным, поскольку открывает возможность «автоматического» учета изменения амплитуды сигнала вследствие изменения проводимости тканей при дыхании.
Оценки параметров матрицы по методу максимального правдоподобия [3] для случая широкополосного шума получаются в результате минимизации функционала:
N г
<Н = агдтіп X X
9Ї т=х,у,г к=1
Л/, [к]-
-
п=х,у,г
>Л^ [к]
(4)
Решение (4) может быть получено как прямыми методами, так и посредством любых итерационных процедур последовательного приближения. Применение последних может быть предпочтительнее, поскольку не возникает проблем с устойчивостью решения.
Обнаружение информационного сигнала. Проведение статистического синтеза метода обработки смеси (2) требует описания сигнала и помехи. Известно, что при нормальных условиях регистрации ЭКГ широкополосная составляющая аддитивной помехи (в основном, представленная электромиограммой) является стационарной и имеет гауссовское распределение. Наличие изо-электрической фазы ЭКГ, в которой отсутствует электрическая активность сердца, позволяет выделить обучающую выборку помехи. Адекватное описание каждой £ -й реализации информационного сигнала ^[к] предложить сложно, поэтому общий подход к решению задачи обнаружения состоит в использовании статистических тестов на однородность характеристик двух выборок. Класс возможных методов ограничивает взаимная коррели-рованность отсчетов помехи и небольшая длительность каждой выборки.
Обработку смеси (2) можно проводить по скалярному сигналу каждого отведения. Алгоритм обработки процесса в на фоне помехи с коррелированными отсчетами основан на декоррелирующем (выбеливающем) преобразовании. Эффективный вычислительный алгоритм реализации выбеливающего фильтра использует представление помехи параметрической авторегрессионной моделью, которая связывает некоторое количество дискретных отсчетов процесса во временной области линейным рекуррентным соотношением р
п[к] + У£агп[к-г] = Ь0х[к], (5)
г=1
где {х[к]} — так называемый порождающий процесс, дисперсия которого равна единице, а отсчеты
взаимно некоррелированны; |{аг}^,Ь0] — параметры модели; р — ее порядок. Оценки параметров модели (5) по обучающей выборке помехи получаются
в результате решения системы уравнений Юла-Уол-кера [4]. Тогда ЧПФ цифрового выбеливающего фильтра определяется следующим образом:
Л + ^агг~
г=1
которая однозначно определяет алгоритм обработки контрольной выборки. После применения выбеливающего преобразования к исходной реализации процесса (2), г\а[к] = )У;={г\[к]}, возникает задача обнаружения сигнала на фоне белого шума с единичной дисперсией. Для принятия решения
о наличии или отсутствии сигнала в выборке Шк]}?, необходимо вычислить статистику, позволяющую осуществить проверку гипотезы о величине дисперсии нормального БШ:
N
т='ЕЫк)у к=1
(6)
и сравнить ее значение с фиксированным пороговым уровнем. г л і
Известно, что оценки параметров \{аг}^,Ь0\ являются состоятельными и асимптотически эффективными. Показано [5], что в асимптотическом случае, при N»^, выражение для потенциально возможных характеристик обнаружения сигнала посредством рассмотренного алгоритма представляется в виде
1
^4лЛ/(1 + 2 д|
оо г (к-N(^ + q2))2
і ехР Ха 4Л/(1 + 2д?) ^ )
с/к,
(7)
где Ха является пороговым уровнем принятия решения для статистики (6) при уровне ложной тре-N / N
вогиРР=а; д2 = ^{э^к])2 / ^[п=[к])2 — отноше-к=1 / к=1 ние сигнал-шум для процесса {г1=[/с]}1 . Чтобы интерпретировать полученное выражение для отношения сигнал-шум в исходном процессе (2), N / N
<72 = 5>и)2/ ^{п[к])2, используется правило к=л / к=1 преобразования спектральных функций процессов и равенство Парсеваля:
я/Гд
я2 і |<%}(/ш)| |и/=(усо7д)| сУш
п/Т,
(8)
о
где С{5}(усо) — спектр сигнала {вц [/с]}™, 7Д —период дискретизации процессов.
Качественные показатели обнаружения зависят от соотношения спектров мощности сигнала и помехи. Наиболее неблагоприятным случаем для обнаружения является полное совпадение спектральных свойств помехи и сигнала.
Результаты эксперимента служат критерием объективной проверки истинности теоретических положений, изложенных в работе, и работоспособности синтезированных алгоритмов и процедур в реальных условиях помеховой обстановки. Практическая оценка качественных показателей алгоритмов затруднена априорной неопределенностью наличия нерегулярных микропотенциалов в реальной контрольной выборке ЭКГ. Обойти данное ограничение позволяет метод полунатурного моделирования: в исходный зарегистрированный процесс добавляется аддитивный компонент информационного сигнала, но никаких изменений в структуре и параметрах алгоритма не производится.
Конкретный пример практического воплощения принятой концепции эксперимента иллюстрируют реализации реально зарегистрированных и обработанных процессов, приведенные на рис. 1-4.
■ Рис. 1. Усредненные сигналы трех ортогональных отведений ЭКГ по Франку
■ Рис. 2. Сигналы в отведении X:
1 — реализация одного из кардиоциклов; 2 — усредненный сигнал {л<х> [/<•]}: регулярная активность; 3 — результат компенсации регулярных составляющих ЭКГ — контрольная выборка процесса, г|„ [/с] = = э11г[к] + пе[к]
■ Рис. 3. Обработка контрольной выборки процесса в отведении X, представленного на рис. 2:
1 — исходная контрольная выборка; 2 — выходной процесс выбеливающего фильтра {Лм 3 —
пороговый уровень С для Рг=0,01\ 4 —решающая статистика (6), вычисляемая на скользящем окне длительностью N = 50 отсчетов
■ Рис. 4. Полунатурное моделирование алгоритмов: в контрольную выборку процесса в отведении X, представленную на рис. 2 и рис. 3 введен искусственный аддитивный сигнал. Отношение сигнал-шум <7=7:
1 — сигнал: два периода синусоиды частотой 100 Гц;
2 — контрольная выборка; 3 — выходной процесс выбеливающего фильтра {г|ц4 — пороговый уровень С для Р)-= 0,01; 5 — решающая статистика (6), вычисляемая на скользящем окне длительностью N = 50 отсчетов
На рис. 1 в увеличенном масштабе отображения приведены усредненные сигналы в трех ортогональных отведениях по Франку. Усреднено
1 = 293 кардиоциклов. Из рисунка видны детали тонкой структуры ЭКГ: на протяжении зубца Р и в конечной части ОРЭ-комплекса. Указанные компоненты обычно и считаются микропотенциалами в традиционном понимании вопроса. Вообще они относятся к составляющим регулярных составляющих ЭКГ.
Из рис. 2 можно получить наглядное представление об эффективности компенсации регулярных составляющих на протяжении произвольно выбранного кардиоцикла. В результате вычитания усред-
ненного сигнала (эпюра 2) из исходной зашумленной реализаций кардиоцикла (эпюра 1) образуется контрольная выборка {г\^[к]}^ (процесс 3), в которой могут присутствовать составляющие нерегулярных микропотенциалов [/с]}^.
Рис. 3 иллюстрирует этапы обработки полученной контрольной выборки. Процесс на выходе выбеливающего фильтра представлен эпюрой 2. Визуальный анализ позволяет констатировать сходность его структуры с процессом типа белого шума. Вычисление решающей статистики (6) на скользящем окне длительностью N = 50 отсчетов (априорно предполагаемая длительность сигнала) соответствует процессу 4. Превышение порогового уровня, установленного исходя из Р[-= 0,01, не наблюдается ни в одной точке, т. е. составляющие нерегулярных микропотенциалов не обнаружены. Это подтверждает работоспособность предложенного подхода: не обнаружено никаких составляющих, которые могли бы быть побочными продуктами предыдущих этапов обработки реализации ЭКГ.
Проверку способности алгоритма к обнаружению сигнала иллюстрируют эпюры на рис. 4. В исходный процесс искусственно вводится информационный сигнал, представляющий собой 2 периода синусоидального сигнала частотой 100 Гц. Исходные данные выбраны, исходя из наихудших условий: длительность сигнала примерно равна периоду корреляции помехи, т. е. его структура мало отличается от одиночных выбросов помехи. Отношение сигнал-шум д составляет 1. Визуальный анализ полученной контрольной выборки, соответствующей процессу 2, не позволяет принять решение о наличии какого-либо информационного сигнала (для облегчения проведения сравнения участков процесса интервал наличия сигнала специально выделен). Но заметно, что дисперсия процесса на выходе выбеливающего фильтра возрастает (эпюра 3). На протяжении интервала наличия сигнала регистрируется устойчивое превышение ре-
шающей статистикой порогового уровня. Результаты моделирования позволяют констатировать отличное совпадение расчетных и экспериментально полученных рабочих характеристик алгоритма обнаружения.
Заключение. Предложенная методика обработки ЭКГ предназначена для проведения диагностики электрической нестабильности миокарда. Необходимая достоверность конечного результата достигается только при оптимальном построении каждого этапа алгоритма обработки. В случае обнаружения микропотенциалов следует решать задачу оценивания их параметров, энергетических и неэнергетических [3], что можно осуществить стандартными методами статистической теории, и соответствующей классификации по виду нерегулярности.
1. Waldo A. L., Kaiser G. A. A study of ventricular arrhythmias associated with acute myocardial infarction in the canine heart. // Circulation. 1973. — Vol. 47. — P. 1222-1228.
2. Лебедева H. А. Алгоритмы анализа многоканальной ЭКГ для классификации форм QRS-комплексов. / Ав-тореф. дисс ... к. т. н. — СПб: ЛЭТИ, 1999. — 16 с.
3. Куликов Е. И., Трифонов А. П. Оценка параметров сигналов на фоне помех. — М.: Сов. Радио, 1978. — 296 с.
4. Марпл С. Л. Цифровой спектральный анализ и его приложения: Пер. с англ. — М.: Мир, 1990. — 564 с.
5. Зайченко К. В., Жаринов О. О., Кулин А. Н., Жа-ринов И. О. Использование принципов адаптивной фильтрации для обнаружения низкоамплитудных биоэлектрических сигналов. // Вестник С.-З. регионального отделения АМТН. Под ред. Г. В. Анцева. / ООО «Агентство РДК-принт». — СПб, 2001. — Вып. 5. — С. 157-162.
ИЗДАТЕЛЬСТВО «ПОЛИТЕХНИКА» ПРЕДСТАВЛЯЕТ
Кухаре в Г. А.
Биометрические системы: методы и средства идентификации личности человека. — СПб.: Политехника, 2001. — 240 с., ил.
В книге дан обзор идей, методов и средств, используемых в системах идентификации и верификации личности человека на основе его биометрических характеристик. Проанализированы состояние и тенденции развития различных биометрических технологий. Подробно обсуждены биометрические системы, использующие информацию об отпечатках пальцев, геометрии руки, параметрах лиц и радужной оболочке глаз. Представлены элементы цифровой обработки изображений, наиболее часто используемые в различных биометрических технологиях, базовые методы обработки исходных данных в приложении к системам распознавания лиц, а также методики построения (синтез, анализ и реализация) систем распознавания людей по их лицам. Приведены словарь терминов, обширная библиография и адреса наиболее важных интернет-сайтов по биометрическим технологиям.
Книга предназначена для студентов и аспирантов соответствующих специальностей, а также практиков, интересующихся использованием представленных в ней биометрических систем.
