УДК 621.391:004.93:616.12-073.97(075.8)621.38:61(0.75.8)
ВЕЙВЛЕТ-АНАЛИЗ ЭЛЕКТРОКАРДИОСИГНАЛОВ ДЛЯ ВЫЯВЛЕНИЯ ПРИЗНАКОВ ПОСТТРАВМАТИЧЕСКОЙ МИОКАРДИОДИСТРОФИИ1
Иванушкина Н.Г.1, к.т.н., доц.; Иванько Е.О.1, к.т.н.; Чеснокова О.В.1;
7 -?
Чайковский И.А. ', к.м.н., д.м.н.(ФРГ)
1 Национальный технический университет Украины «Киевский политехнический институт», Киев, Украина, [email protected] Институт кибернетики НАН Украины, Киев, Украина 3Национальный военно-медицинский клинический центр «Главный военный клинический госпиталь» МО Украины, Киев, Украина
WAVELET ANALYSIS OF THE ELECTROCARDIOSIGNALS FOR DETECTION OF THE POSTTRAUMATIC MYOCARDIAL DYSTROPHY FEATURES
Ivanushkina N., Ivanko K., Chesnokova O., Chaikovskiy I.,
National Technical University of Ukraine "Kyiv Polytechnic Institute", Kyiv, Ukraine Institute for Cybernetics of NAS of Ukraine, Main Military Medical Clinical Center «MMCH»
of Ukraine, Kyiv, Ukraine
Введение
В настоящее время наблюдается рост травматизма при дорожно-транспортных происшествиях, катастрофах, и, особенно, в зонах боевых действий. Травмы приводят не только к нарушению функций поражённого сегмента, но также и к общей реакции организма - травматической болезни. Как известно, травмы любой локализации вызывают заболевания сердечно-сосудистой системы разной степени тяжести [1,2].
Наиболее часто травматическая болезнь со стороны сердечнососудистой системы сопровождается миокардиодистрофией. Понятие «ми-окардиодистрофия» включает обратимые на ранних стадиях нарушения обмена веществ и энергообеспечения в миокарде, ведущие к недостаточности функций сердца при отсутствии структурных изменений [3].
Исследования на клеточном уровне показали, что при разнообразии факторов, воздействующих на сердце, в основе повреждения миокарда лежат нарушения энергообеспечения и транспорта ионов, с которыми связаны биоэнергетические процессы в кардиомиоците. Установлено, что патологические процессы, ведущие к нарушению транспорта ионов Ca2+ и к возникновению дефицита богатых энергией фосфорных соединений, играют основную роль в развитии миокардиодистрофии. При разной этиологии миокардиодистрофии увеличение проницаемости мембран миофибрилл
1 http://radap.kpi.ua/radiotechnique/article/view/1221
для Са2+ или снижение активности кальциевого насоса приводит к повышению концентрации Са2+ в саркоплазме, которое вызывает нарушение расслабления миофибрилл.
Процессы, связанные с кальцием, могут приводить либо к обратимому нарушению функции и структуры миокарда либо к возникновению в нем очагов некробиоза с дальнейшим рубцеванием и становиться основой не-коронарогенного кардиосклероза. Дистрофия миокарда также проявляется недостаточностью сократительной функции сердца; реже наблюдаются нарушения сердечного ритма и совсем редко — нарушения проводимости.
С помощью электрокардиографии можно рано установить изменения, связанные с миокардиодистрофией, так как на ЭКГ отражаются многие нарушения обмена веществ. Признаками миокардиодистрофии на ЭКГ в ранний посттравматический период могут быть: удлинение интервала Q — Т, укорочение длительности и уменьшение амплитуды зубца Т в сравнении с должными значениями, общее уменьшение амплитуды зубцов и комплексов ЭКГ. Однако, в поздний травматический период изменения амплитудно-временных параметров ЭКГ становятся незначительными.
Для выявления низкоамплитудных компонент электрокардиосигналов, которые могут отражать дисфункцию миокарда на клеточном уровне, предлагается применять неинвазивный метод электрокардиографии высокого разрешения (ЭКГ ВР). Целью работы является выявление признаков мио-кардиодистрофии на основе анализа электрокардиосигналов системами ЭКГ ВР.
Особенности посттравматической миокардиодистрофии
Работа основана на клинических исследованиях 2-х групп пациентов: раненных военнослужащих (группа «раненные), которые принимали участие в проведении АТО на востоке Украины, и спортсменов (группа «здоровые»). Военнослужащие (в возрасте в среднем - 33±8 лет) находились на стационарном лечении в хирургических отделениях Национального военно-медицинского клинического центра «Главный военный клинический госпиталь» МО Украины з наличием травматических повреждений средней или тяжёлой степени из-за действия факторов взрыва.
При травматической болезни после минно-взрывной травмы развивается гипоксемия(уменьшение концентрации кислорода в крови) и эндоток-сикоз (поступление в кровоток продуктов распада тканей и сгустков крови), которые сопровождаются нарушениями работы сердечно-сосудистой системы. Как известно, от 6,4 % до 59,6 % случаев травм разной локализации характеризуются развитием посттравматической миокардиодистро-фии [4].
Проявление миокардиодистрофии после взрывных ранений зависит от многих факторов: области поражения (сердца или других органов и тканей), наличия ожогов или химических отравлений, от степени стрессового
состояния организма, и может иметь локальный или диффузный характер.
Согласно [5] уделяется значительное внимание разработке методов диагностики и лечения болезни сердца для раннего травматического периода. Однако, для нормального исхода заболевания и полного выздоровления необходимо развитие методов и средств для оценивания состояния сердца в поздний травматический период на основе обнаружения тонких проявлений активности сердца.
Метод электрокардиографии высокого разрешения
ЭКГ ВР представляет неинвазивный метод (рис. 1), который включает в себя регистрацию электрокардиограммы с частотой дискретизации не менее 1 кГц и цифровую обработку сигнала.
Благодаря высокой частоте дискретизации и высокой разрядности (более 12-ти разрядов) аналого-цифрового преобразователя становится возможным диагностировать низкоамплитудные (до 20 мкВ) [7] высокочастотные составляющие ЭКГ сигнала, которые появляются при нарушениях в работе сердца и не видны на стандартной электрокардиограмме.
Рис. 1. Этапы обработки ЭКГ ВР сигнала: а) исходный зарегистрированный ЭКГ сигнал
первого стандартного отведения, Ь) ансамбль ЭКГ сигнала из кардиоциклов с синхронизацией по К-зубцу, с) выборочный кардиоцикл без усреднения, d) усреднённый
кардиоцикл
Цифровая обработка электрокардиосигнала включает в себя разделение исходного электрокардиосигнала (рис. 1 а) на кардиоциклы с синхронизацией по ^-зубцу (рис. 1 Ь), а также фильтрацию и временное усреднение ансамбля кардиоциклов, представленных на рис. 1 Ь. Фильтрация выполняется для устранения шумов и помех: сетевой наводки с частотой 50 Гц, мышечного тремора с частотой 35 Гц и дрейфа изолинии с частотой до 1 Гц.
Для сравнения, выборочный кардиоцикл (без усреднения) представлен на рис. 1 с. Временное усреднение позволяет значительно снизить соотношение сигнал/шум, так как полезный сигнал периодически повторяется от кардиоцикла к кардиоциклу, а шумовая составляющая во времени хаотична. Следовательно, при усреднении кардиоциклов полезный сигнал сохраняется, а шумовые составляющие частично компенсируют друг друга (рис. 1. d).
Временной анализ
При временном анализе в системах ЭКГ ВР определяются значения низкоамплитудных компонент (поздних потенциалов предсердий или желудочков) на определённом участке усреднённого кардиоцикла (терминальная часть зубца P или зубца 5), которые являются маркерами развития тахиари-тмий сердца.
Однако, такой анализ может быть информативным, если точно знать, на каком участке ЭКГ отражается то или иное нарушение работы сердца. При миокардиодистрофии характер нарушений в сердце может быть как локальным, так и диффузным, поэтому изменения электрокардиограммы неспецифичны и трудно сказать, в какой области сердца и, соответственно, на каком временном интервале электрокардиограммы присутствуют отклонения. В связи с этим, признаки миокардиодистрофии необходимо искать в преобразованных базисах.
Вейвлет-анализ
При клинических исследованиях выполнен дискретный вейвлет-анализ усреднённых кардиоциклов ЭКГ сигналов пациентов 2-х групп: «раненные» и «здоровые».
Вейвлет-преобразование обеспечивает двумерную развёртку, при этом масштаб и координата рассматриваются как независимые переменные, что позволяет анализировать ЭКГ сигнал сразу в двух пространствах — масштабном и временном [7, 8].
При исследовании сигналов с помощью дискретного вейвлет-анализа выполняется их представление в виде совокупности последовательных приближений грубой (аппроксимирующей) ат(р) и уточненной (детализирующей) йт(р) составляющих
с последующим их уточнением итерационным методом. Каждый шаг уточнения соответствует определенному масштабу ат (т.е. уровню т) анализа (декомпозиции) и синтеза (реконструкции) сигнала [7, 8].
Численные эксперименты
В работе использованы электрокардиосигналы второго стандартного отведения, зарегистрированные системой ЭКГ ВР с частотой дисретизации 1 кГц. Поскольку во всех отведениях записаны одни и те же сигналы электрической активности сердца, только снятые с разных точек на поверхности организма человека, анализ нескольких отведений лишь поможет уточнить данные, полученные при анализе одного отведения.
Согласно теории вейвлетов, наиболее информативный результат достигается при использовании материнской функции, морфология которой напоминает исследуемую функцию. Для вейвлет-анализа выбрана материнская функция Мейера, которая при разных временных масштабах очень похожа на зубцы и комплексы ЭКГ сигнала, что даёт возможность исследовать изменения отдельных участков ЭКГ на разных уровнях детализации. Разложение проведено до 9-го уровня.
Для численного анализа результатов вейвлет-разложения были найдены интегральные параметры всех девяти аппроксимирующих ат и детализирующих ёт составляющих согласно формулам
I — V™ \а .|
1ат ¿-1]=1\и-т] \ >
1 йт — =1 \ ^ т у \ ■
где т = 1...9 — уровень разложения, п — количество отсчетов уровня разложения, 1ат, — интегралы под кривыми ат и ёт соответственно.
Поскольку при миокардиодистрофии нарушаются процессы, связанные с транспортировкой ионов Са2+, изменения электрокардиосигналов при вейв-лет-разложении будут видны на высших уровнях разложения [8, 9]. Наиболее показательные результаты интегральных параметров аппроксимирующих составляющих ат и детализирующих составляющих ёт для пациентов из группы «здоровые» и группы «раненные» приведены в табл. 1. и показаны на рис. 2.
Таблица I
Минимальное, максимальное, среднее значения и диапазон изменений значений интеграла (показателей для 8-го и 9-го уровней детализации и 9-го уровня
аппроксимации)
«Здоровые» «Раненные»
Макс. Среднее Мин. Диапаюн Макс. Среднее Мин. Диапазон
170,37 83,39 32,94 137,43 30,82 21,02 9,69 21,12
45,33 18,95 6,65 38,68 6,54 3,19 1,52 5,02
15,37 9,06 1,79 13,58 3,17 2,16 0,87 2,29
а) б) в)
Рис. 2 Значения интегральных параметров аппроксимирующей составляющей а) а9 и детализирующих составляющих б) d8, в) d9 для пациентов из группы «здоровые» и группы
«раненные»
Примеры вейвлет-разложения усреднённых кардиоциклов пациентов из обеих групп представлены на рис. 3.
а)
б)
Рис. 3. Вейвлет-разложение усреднённых кардиоциклов пациентов из групп а) «здоровые» и б) «раненные»
Обсуждение результатов
Обычно считают, что миокардиодистрофия на электрокардиограмме влияет на длительность и амплитуду S-, T-зубцов ST-интервала [3]. Метод ЭКГ ВР позволяет регистрировать сигналы электрической активности сердца с высокой разрешающей способностью по времени и амплитуде, а также рассчитать таблицу амплитудно-временных параметров (амплитуды и длительности зубцов и интервалов) для усреднённой ЭКГ каждого отведения. По полученным параметрам определено, что у 38% пациентов снижена амплитуда Р-зубца и у 85% - амплитуда S-зубца, у 39% - амплитуда Т-зубца, у 54% пациентов значительно увеличена длительность Т-зубца. В то же время, амплитуда усреднённого кардиоцикла у «здоровых» пациентов лежит в диапазоне 1,2...2,5 мкВ, а у «раненных» — 0,3...0,8 мкВ. Данные результаты подтверждают тот факт, что при миокардиодистрофии может быть поражена любая область сердца, а электрокардиосигнал, отражающий патологические процессы, не имеет специфических признаков. Таким образом, для диагностики миокардиодистрофии требуется анализ всего кардиоцикла.
В полученных результатах вейвлет-анализа показано различие численных показателей практически на всех уровнях, но наиболее информативные параметры приведены для высоких уровней разложения (табл. 1), которые отражают низкочастотные составляющие электрокардиосигнала. С помощью вейвлет-анализа можно не только детектировать наличие патологии «миокардиодистрофия», но и разработать алгоритм классификации электро-кардиосигналов для группы «раненные» по степени тяжести.
Выводы
В работе выполнен вейвлет-анализ усреднённых кардиоциклов для пациентов из группы «здоровые» и группы «раненные». На основе дискретного вейвлет-разложения предложена методика выявления признаков посттравматической миокардиодистрофии. В качестве признаков принимаются числовые параметры, которые определяются как интегральные характеристики для самых высоких уровней вейвлет-разложения усреднённых кардиоциклов.
Перечень ссылок
1. Руководство по кардиологии / под ред. В. Н. Коваленко. - К. : МОРИОН, 2008. -1424 с.
2. ilhan E. Electrocardiographic features of patients with earthquake related posttraumatic stress disorder / E. ilhan, A. Kaplan, T. S. Güven9, M. Biteker, E. Karabulut, S. I§ikh // World Journal of Cardiology. - 2013. - Vol. 5, No 3. - pp. 60-64.
3. Исаева А. Д. Миокардиодистрофия. Причины, симптомы, признаки, диагностика и лечение патологии [Электронный ресурс] / А. Д. Исаева. - Режим доступа: http://www.polismed.com/articles-miokardiodistrofija-prichiny-simptomy-priznaki.html.
4. Шок: Теория, клиника, организация противошоковой помощи / под общ. ред. Г. С. Мазуркевича, С. Ф. Багненко. - СПб. : Политехника, 2004, - 539 с.
5. Военно-полевая хирургия [Электронный ресурс] / под ред. Е. К. Гуманенко. -М.: ГЭОТАР-Медиа, 2008. - Режим доступа : http://www.studmedlib.ru/book/ ISBN9785970409114.html
6. 1ванушкша Н. Г. Цифрова обробка низькоампштудних компонент електрокардь осигналiв / Н. Г. 1ванушкша, К. О. 1ванько. - Миколшв : ФОП Швець В.Д., 2014. - 184с.
7. Смоленцев Н. К. Основы теории вейвлетов. Вейвлеты в MATLAB / Н. К. Смоле-нцев. - М. : ДМК Пресс, 2005. - 304 с.
8. Haque A. K. Detection of small variations of ECG features using wavelet / A. K. Haque, M. H. Ali, M. A. Kiber, M. T. Hasan // ARPN Journal of Engineering and Applied Sciences. - 2009. - Vol. 4, No 6. - pp. 27-30.
9. Ishikawa Y. Wavelet Theory-Based Analysis of High-Frequency, High-Resolution Electrocardiograms: A New Concept for Clinical Uses / Y. Ishikawa // Progress in Biomedical Research. - 2002. - Vol. 7, No 3. - pp. 179-184.
References
1. Kovalenko V. N. ed. (2008) Rukovodstvo po kardiologii [Guidelines for Cardiology]. Kyiv, MORION, 1424 p.
2. ilhan E., Kaplan A., Guven9 T.S., Biteker M., Karabulut E., and I§ikli S. (2013) Electrocardiographic features of patients with earthquake related posttraumatic stress disorder. World Journal of Cardiology, Vol. 5, No 3, pp. 60-64.
3. Isaeva A. D. Myocardiodystrophy. Causes, symptoms, signs, diagnosis and treatment of disease. Available at: http://www.polismed.com/articles-miokardiodistrofija-prichiny-simptomy-priznaki.html.
4. Mazurkevicha G. S. eds, Barnenko S. F. (2004) Shok: Teoriya, klinika, orranizatsiya protivoshokovoi pomoshchi [Shock: Theory, clinic, organization antishock aid]. - SPb, Politekhnika, 539 p.
5. Gumanenko E. K. eds. (2008) Voenno-polevaya khirurgiya [Military Field Surgery]. Moskow, GEOTAR-Media
6. Ivanushkina N. H., and Ivanko K. O. (2014) Tsyfrova obrobka nyzkoamplitudnykh komponent elektrokardiosyhnaliv [Digital processing of low amplitude components of electrocardiosignal]. Mykolaiv, FOP Shvets V.D., 184 p.
7. Smolentsev N. K. (2005) Osnovy teorii veivletov. Veivlety v MATLAB [Fundamentals of wavelet theory. Wavelets in MATLAB]. Moskow, DMK Press, 304 p.
8. Haque A. K., Ali M. H., Kiber M. A., and Hasan M. T. (2009) Detection of small variations of ECG features using wavelet. ARPN Journal of Engineering and Applied Sciences, Vol. 4, No 6, pp. 27-30.
9. Ishikawa Y. Wavelet Theory-Based Analysis of High-Frequency, High-Resolution Electrocardiograms: A New Concept for Clinical Uses. Progress in Biomedical Research, 2002, vol. 7, № 3, pp. 179-184.
1ванушюна Н. Г., 1ванько С. О., Чеснокова О. В., Чайковський I. А. Вейвлет-анал1з електрокард1осигнал1в для виявлення ознак посттравматичноИ м1окард1одистро-
фИ. Розглянуто особливост1 впливу травматичног хвороби на серцево-судинну систему. Виконано числовг експерименти з досл1дження електрокардгосигналгв вгйськово-службовщв тсля м1нно-вибухового поранення. У процес1 клнчних досл1джень реестращя та обробка електрокард1осигнал1в виконаш за допомогою системи елек-трокардюграфи високого розр1знення. На основ1 проведеного вейвлет-аналгзу електро-кард1осигнал1в запропоновано 1нтегральн1 параметри для виявлення ознак посттрав-матичног мюкардюдистрофи.
Ключов1 слова: посттравматична мiокардiодистрофiя, система електро-кардюграфи високого розрiзнення, вейвлет-аналiз електрокардiосигналiв, низькоам-плтудт компоненти ЕКГ.
Иванушкина Н. Г., Иванько Е. О., Чеснокова О. В., Чайковский И. А. Вейвлет-анализ электрокардиосигналов для выявления признаков посттравматической миокардиодистрофии. Рассмотрены особенности влияния травматической болезни на сердечно-сосудистую систему. Выполнены численные эксперименты по исследованию электрокардиосигналов военнослужащих после минно-взрывных ранений. В процессе клинических исследований регистрация и обработка электрокардиосигналов выполнены с помощью системы электрокардиографии высокого разрешения. На основе проведенного вейвлет-анализа электрокардиосигналов предложены интегральные параметры для выявления признаков посттравматической миокардиодистрофии.
Ключевые слова: посттравматическая миокардиодистрофия, система электрокардиографии высокого разрешения, вейвлет-анализ электрокардиосигналов, низкоамплитудные компоненты ЭКГ.
Ivanushkina N., Ivanko K., Chesnokova O., Chaikovskiy I., Wavelet analysis of the elec-trocardioignals for detection of the posttraumatic myocardial dystrophy features.
Introduction. The features of the traumatic disease's impact on the cardiovascular system are reviewed. Numerical experiments on electrocardiosignals of the soldiers injuried by mine explosion were held. The nature of the posttraumatic myocardial dystrophy is described.
HR ECG signals' analysis. During clinical studies the registration and processing of electrocardiosignals were performed by the high resolution electrocardiography system. The time and amplitude analysis of HR ECG signals was held, and its low efficiency was approved. The integral parameters for detection of the posttraumatic myocardial dystrophy features based on the wavelet analysis of the electrocardiosignals are proposed.
Conclusion. Since the myocardial dystrophy is caused by the violation of the processes related to the transport of ions Ca2+, the changes of electrocardiosignals are visible at the highest levels of wavelet decomposition. Thus, the 8th and 9th levels of the 9-level wavelet decomposition of averaged cardiocycles showed the best result for diagnostic.
Keywords: posttraumatic myocardial dystrophy, high resolution electrocardiography system, wavelet analysis of the electrocardiosignals, low-amplitude components of ECG.