УДК 61.007+004.932.72 DOI: 10.20998/2411-0558.2017.21.13
Л.С. ФАЙНЗИЛЬБЕРГ, д-р техн. наук, доц., Международный научно-учебный центр информационных технологий и систем НАН Украины и МОН Украины, Киев,
К.Б. ОРИХОВСКАЯ, асп., Международный научно-учебный центр информационных технологий и систем НАН Украины и МОН Украины, Киев
НОВЫЙ ПОДХОД К ОБНАРУЖЕНИЮ ЭФФЕКТА ЭЛЕКТРИЧЕСКОЙ АЛЬТЕРНАЦИИ СЕРДЦА ПО ОДНОКАНАЛЬНОЙ ЭЛЕКТРОКАРДИОГРАММЕ
Рассматривается метод обнаружения эффекта электрической альтернации сердца, основанный на анализе перестановочной энтропии, которая вычисляется по массиву показателей одноканальной электрокардиограммы в скользящих окнах. Представлены результаты исследования метода на модельных и реальных данных. Ил.: 3. Табл.: 1. Библиогр.: 10 назв.
Ключевые слова: электрическая альтернация сердца; перестановочная энтропия; электрокардиограмма.
Постановка проблемы. Один из предикторов внезапной сердечной смерти, который в последнее время получил широкую известность в клинических исследованиях, основан на анализе, так называемой, электрической альтернации сердца [1 - 3]. В то же время, по мнению специалистов, существующие компьютерные алгоритмы не обеспечивают требуемую надежность обнаружения этого эффекта в реальных клинических условиях.
Сложность решения задачи обусловлена тем, что реальные сигналы c наличием и отсутствием электрической альтернации внешне практически неразличимы (рис. 1) [4].
Рис. 1. Электрокардиограмма со случайным искажением (вверху) и альтернацией
(внизу) зубца Т
Согласно [5], диагностически значимый уровень альтернации зубца
© Л.С. Файнзильберг, К.Б. Ориховская, 2017 144
Т составляет всего 5 - 10 мкВ. Надежное обнаружение таких тонких изменений сигнала на фоне возможных искажений, вызванных случайными возмущениями, требует применения эффективных компьютерных алгоритмов.
Анализ литературы. Один из известных подходов к построению компьютерных алгоритмов анализа альтернации зубца основан на алгоритме быстрого преобразования Фурье и обнаружении характерного пика на частоте, равной половине частоты сердечных сокращений [6]. В то же время такой алгоритм может давать приемлемые результаты лишь в том случае, когда амплитуда зубца Т изменяется по синусоидальному закону.
Для повышения достоверности результатов в работе [7] предложено оценивать степень альтернации на основе спектральных составляющих разложения по функциям Уолша. Такой алгоритм позволяет обнаружить эффект альтернации, который проявляется не только в чередовании амплитуд, но и симметрий зубца Т .
В клинической практике наблюдаются также другие типы альтернации, в частности, чередование зубцов Т с разной продолжительностью [8]. Возможны также сочетания различных типов альтернации на одной электрокардиограмме (ЭКГ).
Практика показывает, что в реальных ситуациях фаза чередования зубцов может изменяться в случайные моменты времени. В таких ситуациях анализ альтернации на основе спектральных разложений приводит к ошибочным результатам.
Кроме того, спектральные методы обладают существенным недостатком: для получения корректных оценок требуется с помощью дополнительных стимулирующих катетеров в течение продолжительного времени (более 5 минут) поддерживать стабильные значения частоты скращения сердца (ЧСС) на высоких уровнях (не менее 100 уд/мин), что не всегда безопасно для пациента [6].
Поэтому получили известность методы анализа альтернации зубца Т, основанные на обработке ЭКГ во временной области, в частности, на оценке усредненных значений амплитуд четных и нечетных зубцов Т [5]. Однако такой алгоритм эффективен лишь при четком чередовании зубцов различных типов.
Для расширения возможностей компьютерного обнаружения эффекта электрической альтернации сердца в [9] предложен более совершенный алгоритм, основанный на анализе формы гистограммы, характеризующей распределение параметров зубцов. Признаком альтернации является двумодальность гистограммы и выполнение дополнительных условий, налагаемых на подмножества индексов элементов ЭКГ.
Этот алгоритм основан на предположении, что при электрической альтернации амплитуды зубцов T порождаются двумя случайными величинами Yj и Yg, которые с ограниченными дисперсиями варьируют относительно "своих" средних значений M{Yj} Ф M{Yg }, причем
Ya е [Yj™ YAmax], Yg е [Ygmin, Ygmax] (1)
и
Ymax < Ymin или Vmax < Vmin (2)
YA < YB или YB < ya • (2)
Однако из-за случайных возмущений в реальных условиях ограничение (2) может не выполняться.
Цель статьи - разработать надежный алгоритм обнаружения эффекта альтернации, основанный на анализе хаотичности последовательностей амплитудно-временных параметров элементов ЭКГ, и оценить эффективность алгоритма на модельных и реальных данных.
Математические основы предлагаемого метода. Пусть
аь a2,•••, ам (3)
- массив значений параметров, наблюдаемых на последовательности сердечных циклов, например, массив значений амплитуд зубцов T •
Разобьем массив (3) на L последовательных окон, в каждом таком l -м окне оценим хаотичность Hl значений исследуемого показателя и вычислим отношение Hl к хаотичности Hi элементов в первом окне:
hl = -100%, l =1,...,L , (4)
H1
полагая, что Hi Ф 0.
Для оценки Hl будем использовать алгоритм вычисления перестановочной энтропии
5
Hl ="Z p(p jl )lo§2 p(p jl h (5)
j=1
в которой p(p jl) - частота появления паттерна j -го класса в l -м окне.
Перестановочная энтропия (5) характеризует разнообразие формы пяти классов паттернов (рис. 2), образованных тройкой последовательных значений ат_1, am, am+1 массива (3).
000Н0
Р1 %2 %3 п4 п5
Рис. 2. Классы паттернов перестановочной энтропии
Для обеспечения плавности изменения энтропии процедура (4) реализуется при сдвиге / +1 -го окна по отношению / -му на одну точку (режим скользящего окна).
Применение специальных вычислительных процедур позволяет
отобразить последовательность Н/ на плоскости И(/), Л(/), где /5(/) -оценка первой производной к(/) в / -й точке, и определить площадь £ выпуклой оболочки, построенного таким образом фазового портрета энтропии (ФПЭ) (рис. 3).
А(/),% 101 ■
&9
94
■г*'"
X
г,-
аэа
0.73 0.47
21 -А 13 -¡>.05 о.оа
/г(/),%/с
-□.05-□
£
хч
\
/
В [У
06 0.17 О^Ж 0.61 а,83 1,06
Рис. 3. ФПЭ (слева) и его выпуклая оболочка (справа)
Модельные эксперименты. Для оценки эффективности предложенного подхода использовалась генеративная модель порождения искусственных ЭКГ реалистической формы [10], с помощью которой генерировались тестовые сигналы с различными значения уровня альтернации амплитуд зубца Т на фоне 15 % случайных искажений.
Несмотря на то, что модельные сигналы были практически неразличимы, площадь £ выпуклой оболочки фазового портрета энтропии (ФПЭ) монотонно уменьшалась по мере нарастания альтернации зубца Т , т.е. по мере увеличения доли регулярной составляющей в изменчивости амплитуд.
При уровне альтернации 60 мкВ площадь £ уменьшилась более чем на 30 % (от 0,89 ед. до 0,62 ед.) по сравнению с сигналом без альтернации (табл. 1).
Таблица 1
Площади S выпуклых оболочек ФПЭ
Адекватные результаты получены также и при обработке реальных ЭКГ. Площадь выпуклой оболочки ФПЭ, вычисленная по ЭКГ пациента с альтернацией продолжительностей RR -интервалов была почти в 4 раза ниже, чем у здорового пациента с внешне похожей ЭКГ.
Выводы. Предложен новый метод обнаружения электрической альтернации сердца, который, в отличие от ранее известных, основан на оценке изменения значения перестановочной энтропии по ходу наблюдения сигнала. Применение специальных вычислительных процедур позволяет оценить первую производную перестановочной энтропии и вычислить площадь выпуклой оболочки ее фазового портрета
(ФП).
Проведенные исследования подтвердили, что площадь выпуклой оболочки ФП позволяет надежно обнаружить электрическую альтернацию сердца, в том числе RR -интервалов и симметрии зубца T , и использовать этот показатель в качестве предиктора внезапной сердечной смерти.
Список литературы: 1. Файнзильберг Л. С. Основы фазаграфии / Л. С. Файнзильберг. -Киев: Освита Украины, 2017. - 264 с. 2. Sakabe K. Comparison of T-wave alternans and QT interval dispersion to predict ventricular tachyarrhythmia in patients with dilated cardiomyopathy and without antiarrhythmic drugs: a prospective study / K. Sakabe, T. Ikeda, T. Sakata // Japanese Heart Journal. - 2001. - Vol. 42. - №. 4. - P. 451-457.
3. Bloomfield D.M. Microvolt T-Wave Alternans Distinguishes Between Patients Likely and Patients Not Likely to Benefit From Implanted Cardiac Defibrillator Therapy / D.M. Bloomfield, R.S. Steinman // Circulation. - 2004. - №. 110. - P. 1885-1889.
4. Predescu D. Microvolt T wave alternans (MTWA) - a new non-invasive predictor of sudden cardiac death / D. Predescu, P. Mitrut, A. Giuca // Romanian Journal of Internal Medicine. - 2004. - Vol. 42. - №. 3. - P. 647-656. 5. Nearing B.D. Modified moving average analysis of T-wave alternans to predict vetricular fibrillation with high accuracy / B.D. Nearing, R.L. Verrier // Journal of Applied Physiology. - 2002. - Vol. 92. - P. 541-549. 6. Rosenbaum D.S. Electrical alternans and vulnerability to ventricular arrhythmias
148
/ D.S. Rosenbaum, L.E. Jackson, J.M. Smith // New England Journal of Medicine. - 1994. -№ 330. - P. 235-241. 7. Мельник О.В. Метод выявления альтернации T-зубца электрокардиосигнала / О.В. Мельник // Биомедицинская радиоэлектроника. - 2008. -№ 7. - С. 56-59. 8. Verrier R.L. Median beat analysis of T-wave alternans to predict arrhythmic death after myocardial infarction: results from the Autonomic Tone and Reflexes after Myocardial Infarction Study / R.L. Verrier, B.D. Nearing, M.T. LaRovere, G.Pinna, M.A. Mittleman, J.T. Bigger, P.J. Schwartz // Circulation. - 2000. - № 102. - P. 711-713. 9. Файнзшьберг Л.С. Споаб ощнювання електрично1 альтернацп серця / Л.С. Файнзшьберг, В.П. Сакун, В.В. Соловйов // Патент Украши на винахвд № 101918. -2013. - Бюл. № 9. 10. Fainzilberg L.S. T-Wave Alternats Modeling on Artificial Electrocardiogram with Internal and External Perturbation / L.S. Fainzilberg, T.Yu. Bekler // Journal of Automation and Information Sciences. - 2012. - Vol. 44. - Issue 7. - P. 1-14. References:
1. Fainzilberg, L.S. (2017), Fasegraphy basics, Osvita Ukrainy, Kyiv, 264 p.
2. Sakabe, K., Ikeda, T. and Sakata, T. (2001), "Comparison of T-wave alternans and QT interval dispersion to predict ventricular tachyarrhythmia in patients with dilated cardiomyopathy and without antiarrhythmic drugs: a prospective study", Japanese Heart Journal, Vol. 42, No. 4, pp. 451-457.
3. Bloomfield, D.M. and Steinman, R.S. (2004), "Microvolt T-wave alternans distinguishes between patients likely and patients not likely to benefit from implanted cardiac defibrillator therapy", Circulation, Vol. 110, pp. 1885-1889.
4. Predescu, D., Mitrut, P. and Giuca, A. (2004), "Microvolt T wave alternans (MTWA) - a new non-invasive predictor of sudden cardiac death", Romanian Journal of Internal Medicine, Vol. 42, No. 3, pp. 647-56.
5. Nearing, B.D. and Verrier, R.L. (2002), "Modified moving average analysis of T-wave alternans to predict vetricular fibrillation with high accuracy", Journal of Applied Physiology, Vol. 92, pp. 541-9.
6. Rosenbaum, D.S., Jackson, L.E. and Smith, J.M. (1994), "Electrical alternans and vulnerability to ventricular arrhythmias", New England Journal of Medicine, Vol. 330, pp. 235-41.
7. Melnik, O.V. (2008), "The Method for Detection of Electrical Cardiosignal T-Wave Alternation", Biomedical Radioelectronics, Vol. 7, pp. 56-9.
8. Verrier, R.L., Nearing, B.D., LaRovere, M.T., Pinna, G., Mittleman, M.A., Bigger, J.T. and Schwartz, P.J. (2000), "Median beat analysis of T-wave alternans to predict arrhythmic death after myocardial infarction: results from the autonomic tone and reflexes after myocardial infarction study", Circulation, Vol. 102, pp. 711-713.
9. Fainzilberg L.S., Sakun V.P., Soloviov V.V. (2013), Method for evaluating the electrical heart alternation, Patent UA, No. 101918.
10. Fainzilberg, L.S. and Bekler, T.Yu. (2012), "T-Wave alternats modeling on artificial electrocardiogram with internal and external perturbation", Journal of Automation and Information Sciences, Vol. 44, No. 7, pp. 1-14.
Статью представил д-р техн. наук, проф. каф. ВТП НТУ "ХПИ" Поворознюк А. И.
Поступила (received) 14.04.2017
Fainzilberg Leonid, Doctor of Technics.
International Research and Training Center for Information Technology
and Systems, NAS and MES of Ukraine
Str. Heroyiv Dnipra, 36-17, Kyiv, Ukraine, 04214
BiCHUK Нацiонапbного технiнного yHiBepcumemy "Xni", 2017, № 21 (1243)
Tel.: (067) 715-15-83, e-mail: fainzilberg@voliacable.com ORCID ID: 0000-0002-3092-0794
Orikhovska Kseniia, PhD stud.
International Research and Training Center for Information Technology
and Systems, NAS and MES of Ukraine
Str. Borshagivska, 139/141-20, Kyiv, Ukraine, 03056
Tel.: (095) 619-90-03, e-mail: kseniaor@gmail.com
УДК 61.007+004.932.72
Новий пщхщ до виявлення ефекту електричноТ альтернацп серця по одноканальнш електрокардiограмi / Файнзшьберг Л.С., Орiховська К.Б.
// Вiсник НТУ "ХШ". Серiя: 1нформатика та моделювання. - Харюв: НТУ "ХП1". - 2017. - № 21 (1243). - С. 144 - 151.
Розглядасться метод виявлення ефекту електрично! альтернацп серця, заснований на аналiзi перестановно! ентропп, яка обчислюеться по масиву показниюв одноканально! електрокардiограми в ковзних вшнах. Представленi результати дослiдження методу на модельних i реальних даних. 1л.: 3. Табл.: 1. Бiблiогр.: 10 назв.
Ключовi слова: електрична альтернащя серця; електрокардiограма; перестановна ентротя.
УДК 61.007+004.932.72
Новый подход к обнаружению эффекта электрической альтернации сердца по одноканальной электрокардиограмме / Файнзильберг Л.С., Ориховская К.Б. // Вестник НТУ "ХПИ". Серия: Информатика и моделирование. - Харьков: НТУ "ХПИ". - 2017. - № 21 (1243). - С. 144 - 151.
Рассматривается метод обнаружения эффекта электрической альтернации сердца, основанный на анализе перестановочной энтропии, которая вычисляется по массиву показателей одноканальной электрокардиограммы в скользящих окнах. Представлены результаты исследования метода на модельных и реальных данных. Ил.: 3. Табл.: 1. Библиогр.: 10 назв.
Ключевые слова: электрическая альтернация сердца; электрокардиограмма; перестановочная энтропия.
УДК 61.007+004.932.72
A new approach to detecting the electrical alternation effect of the heart on a single-channel electrocardiogram / Fainzilberg L.S., Orikhovska K.B.
// Herald of the National Technical University "KhPI". Subject issue: Information Science and Modelling. - Kharkov: NTU "KhPI". - 2017. - №. 21 (1243). - P. 144 -151.
The method of detecting the electrical alternation effect of the heart, based on the permutation entropy analysis, which is calculated from the array of single-channel electrocardiogram parameters in sliding windows, is considered. The results of the research method on model and real data are presented. Figs.: 3. Tabl.: 1. Refs.: 10 titles.
Keywords: electric alternation of the heart; electrocardiogram; permutation entropy.