CC BY
24
ПРИБОРЫ, СИСТЕМЫ И ИЗДЕЛИЯ МЕДИЦИНСКОГО НАЗНАЧЕНИЯ
УДК 615.47:616-072.7 DOI 10.21685/2307-5538-2020-2-8
М. Н. Крамм
ЭКВИВАЛЕНТНЫЙ ЭЛЕКТРИЧЕСКИЙ ГЕНЕРАТОР СЕРДЦА ДЛЯ НЕИНВАЗИВНОЙ ЭЛЕКТРОКАРДИОДИАГНОСТИКИ
M. N. Kramm
EQUIVALENT ELECTRIC HEART GENERATOR FOR NON-INVASIVE ELECTROCARDIAC DIAGNOSTICS
Аннотация. Актуальность и цели. Рассматривается способ повышения эффективности неинвазивной кардиодиагностики, основанный на переходе от многоканальных электрокардиосигналов и связанных с ними карт распределения электрического потенциала на поверхности торса к пространственно-временным характеристикам эквивалентного электрического генератора сердца (ЭЭГС). Материалы и методы. Проводится реконструкция распределения эквивалентных электрических источников в сердце в течение кардиоцикла. При обработке используются как записанные электрокардиосигналы, так и координаты электродов на поверхности торса человека. Результаты. Представлены карты распределений электрического потенциала на поверхности сердца в различные моменты времени для ЭЭГС поверхностного типа и пространственные треки временной динамики электрической активности сердца в объеме сердца для ЭЭГС дипольного типа. Выводы. Результаты реконструкции ЭЭГС показывают возможность наблюдения карт пространственно-временной динамики электрической активности сердца и предназначены для повышения информативности электрокардиографических обследований.
Abstract. Background. We consider a method for increasing the effectiveness of noninvasive cardiodiagnostics, based on the transition from multichannel electrocardiosignals and related electrical potential distribution maps on the torso surface to the spatial-temporal characteristics of the equivalent electric heart generator (EEHG). Materials and methods. We carry out reconstruction of the distribution of the equivalent electrical heart sources of the cardiac cycle. When processing we use both recorded electrocardiosignals and the coordinates of the electrodes, placed on the human torso surface. Results. We obtain maps of electric potential distributions on the heart surface at various time points for surface type EEHG and spatial tracks of the temporal dynamics of the heart electrical activity in the volume of the heart for dipole type EEHG. Conclusions. The results of the EEHG reconstruction show the possibility of observing the spatio-temporal dynamics maps of the electrical heart activity and are intended to increase the informativeness of electrocardiographic examinations.
Ключевые слова: электрокардиографические отведения, электрический потенциал, торс, сердце, реконструкция, эквивалентный электрический генератор сердца, _карты распределений._
© Крамм М. Н., 2020
Измерение. Мониторинг. Управление. Контроль
Keyword s: electrocardiographic leads, electric potential, torso, heart, reconstruction,
equivalent electric heart generator, maps of distributions.
Введение
Современная неинвазивная электрокардиология широко использует информационные технологии (ИТ), обеспечивающие регистрацию, сбор, хранение и обработку кардиографической информации. В настоящее время этот своеобразный технологический «конвейер» неин-вазивной электрокардиологии помогает врачу при формировании заключения о состоянии сердца пациента [1, 2].
В качестве важного этапа этого конвейера следует рассматривать решение обратной задачи электрокардиографии путем реконструкции эквивалентного электрического генератора сердца (ЭЭГС), которое позволяет неинвазивным путем визуализировать характеристики электрической активности в области сердца [3-5]. При этом повышение достоверности элек-трокардиодиагностики в отношении получения новой диагностической информации связано с одновременным представлением кардиологам пространственно-временных характеристик как ЭЭГС поверхностного типа — для диагностики нарушения процессов проводимости, так и ЭЭГС дипольного типа — для диагностики ишемии. Предлагается переход от распределений потенциала электрического поля на поверхности торса (карт потенциала на поверхности торса) - к картам распределения источников электрической активности на поверхности эпикарда и пространственно временным характеристикам электрической активности ЭЭГС дипольного типа в объеме сердца [6]. Существенно также, что такой переход ориентирован на наблюдение более мелких деталей в динамике электрической активности сердца.
Постановка задачи
Многоканальные электрокардиосигналы (ЭКС) регистрируются с помощью системы электродов, распределенных на поверхности грудной клетки в количестве 40-60. После цифровой фильтрации мы располагаем для каждого отсчетного момента времени кардиоцикла массивом значений электрического потенциала в точках расположения электродов. Поскольку электрический потенциал, зарегистрированный с некоторого электрода, является интегральной характеристикой источников электрической активности в сердце, то для получения более детальной кардиографической информации мы решаем обратную задачу - проводим реконструкцию эквивалентного электрического генератора сердца. При этом существенной информацией являются координаты электродов на поверхности торса и опорных точек на поверхности эпикарда. Координаты электродов определяются путем измерения антропометрических параметров на торсе. Координаты точек на поверхности эпикарда определяются путем реконструкции этой поверхности на основе известных координат электродов и электрических потенциалов на этих электродах [7].
Реконструкция эквивалентного электрического генератора сердца поверхностного типа
Электрические потенциалы регистрировались с помощью многоэлектродной системы отведений, показанной на рис. 1 [3]. Электроды были размещены на эластичных поясах на торсе в количестве около 40 в 4-5 рядов (рис. 1,а). После антропометрических измерений определялись координаты электродов х, у и z (рис. 1,а и б).
Пусть r = [xt,yi,zi], i = 1..Ne — векторы с координатами электродов на поверхности торса. Далее осуществляется переход к достаточно мелкой сетке граничных элементов на поверхности торса с координатами центров этих элементов [хк,ук,zk], к = 1..Nb, причем Nb » Ne. На следующем этапе необходимо путем интерполяции перейти от электрических потенциалов фь, i = 1. .Ne, измеренных в точках расположения электродов ri, к потенциалам фк, к = 1..Nb на мелкой сетке граничных элементов с центрами гк =[хк, ук, zk ] . Переход к мелкой сетке граничных элементов путем выполнения интерполяции важен для расширения частотного диапазона пространственного спектра распределения потенциала на эпикарде [8].
66
а)
б)
Рис. 1. Многоэлектродная система отведений: а - расположение электродов на торсе; б - координаты электродов
Наиболее распространенным способом отображения электрической активности сердца на базе ЭЭГС поверхностного типа (ЭЭГС ПТ) является визуализация распределения электрического потенциала на поверхности эпикарда [7]. В этом случае говорят о распределении поверхностной плотности дипольного момента электрических источников, которое отражает распределение эквивалентного трансмембранного потенциала на поверхности сердца.
В работе [9], основываясь на теореме Грина [10], показано, что
2п ■ и(М) = | g(Р)
1
Щ(М, Р)
■I Ф( Р) ■ I;
(
1
\
Щ(М, Р)
I«
(
1
\
Щ(М, Р)
ёБ,
где ф - электрический потенциал; М - точка, которая может находиться на поверхности торса Бв или на поверхности эпикарда Бн ; Р - текущая точка интегрирования; Щ(М,Р) -расстояние между точками М и Р; Э/Эп - производная по направлению нормали, внешней
Эф
по отношению к области между поверхностью торса и поверхностью эпикарда; g = уэп ■
При этом в случае нахождения точки М на поверхности интегрирования эта точка должна быть «выколота», т.е. соответствующий интеграл следует рассматривать как несобственный интеграл второго рода.
В рамках метода граничных элементов (МГЭ) [11] в результате дискретизации поверхностей Бв и Бн получаем систему линейных алгебраических уравнений, которая может быть решена итерационным методом Зейделя [12] с последовательным использованием обратных матриц.
В качестве примера рассмотрим результаты реконструкции ЭЭГС ПТ по реальным ЭКС. ЭКС зарегистрированы с помощью 36 электродов, размещенных в четыре ряда на поверхности грудной клетки. Поверхность торса представлялась с помощью сетки, содержащей 30 элементов по вертикали и 60 элементов в горизонтальной плоскости. Использовалась реалистическая триангуляционная модель эпикарда, содержащая 1584 элемента. Реконструкция проводилась в характерные моменты времени, соответствующие вершинам Р, Щ и Т зубцов.
Результаты реконструкции ЭЭГС ПТ для данных моментов времени представлены на рис. 2 в виде карт распределения электрического потенциала на поверхности эпикарда.
Рис. 2. Потенциал на поверхности эпикарда для вершины Р-зубца (а); Р-зубца (б); Т-зубца (в)
Для большей наглядности полутоновая палитра на представленных рисунках выбрана с одинаковым контрастом, который соответствует различным диапазонам изменения электрического потенциала в мВ. Как и следовало ожидать, наибольшие значения потенциала наблюдаются для вершины Р-зубца, наименьшие - для вершины Р-зубца. Наиболее заметна диполь-ная структура электрического источника для Р-зубца. Располагая значениями потенциалов в различные моменты времени кардиоцикла, можно наблюдать ЭКС непосредственно на поверхности эпикарда.
Реконструкция эквивалентного электрического генератора сердца дипольного типа
С целью наблюдения временной динамики электрической активности сердца ЭЭГС дипольного типа (ЭЭГС ДТ) реконструируется для каждого отсчетного момента времени кардиоцикла. Исходными данными для реконструкции ЭЭГС для момента tk является вектор измеренных ЭКС Uk = (U1(tk),...,Un(tk),..UN (tk)), где Un(tk) - ЭКС, снимаемый с n-го
электрода, N¡ - количество используемых электродов. Также должны быть измерены координаты электродов (xln, yln, zln), n e (1, N¡) в системе координат торса. ЭЭГС дипольного типа
для всех временных отсчетов кардиоцикла tk характеризуется вектором параметров
sk = (Xsk, ysk, Zsk, Mxk, Myk, Mzk ) где (Xsk, ysk, Zsk ) - ко°рдинаты ЭЭЩ (Mxk , Myk, Mzk ) -
проекции вектора дипольного момента ЭЭГС [6].
Поскольку задача реконструкции относится к классу математически некорректных задач [13], то для ее решения необходимо решить оптимизационную задачу - найти для всех временных отсчетов кардиоцикла 4 вектор параметров ЭЭГС, при котором достигается минимум функционала
"а=| К - и (нк )||2 + I2,
где и(нк) = (и1(,нк),...,ип(нк),...им (нк)) - вектор ЭКС, рассчитанных по дипольной модели ЭЭГС с вектором параметров ,нк ; ип (нк) - модельный ЭКС для п-го электрода; ак - коэффици-
х к - х у . - у г . - г М . М . М .
н к се У ак ^ се ак се хк у к г к
ент регуляризации; нк =
(х -хе у к -усе г„,, -М„,, Мук М_,, \ Ян ' Ян ' Ян 'М1ШМ1ШМти
нормированный век-
тор параметров модели ЭЭГС; (хсе,усек,гсе) - координаты центра модели эпикарда пациента;
Лн - усредненный радиус эпикарда; Мш - модуль вектора дипольного момента ЭЭГС диполь-ного типа для временного отсчета максимума Л-зубца кардиоцикла. Вопрос выбора коэффициента регуляризации рассмотрен нами в работе [14].
Временная динамика поведения ЭЭГС ДТ в ходе кардиоцикла может наглядно отображаться пространственными кривыми, характеризующими движение источника в пространстве и движение конца вектора дипольного момента. В качестве примера на рис. 3 показаны треки -линии движения источника для Р-, Л- и Т-зубцов во фронтальной плоскости.
Выводы и заключение
Представленные алгоритмы ориентированы на анализ динамики электрической активности сердца путем пространственно-временного картирования эквивалентных электрических источников сердца. Автором предложено синэнэргетическое объединение двух подходов к описанию электрической активности сердца: 1) эквивалентный электрический генератор поверхностного типа (ЭЭГС ПТ), характеризующий электрические процессы на поверхности сердца и описываемый пространственно-временными картами электрического потенциала на поверхности эпикарда; 2) эквивалентный электрический генератор дипольного типа (ЭЭГС ДТ), интегрально характеризующий электрические процессы в объеме сердца и описываемый пространственно-временными характеристиками - треками движения ЭЭГС и динамикой изменения модуля и направления вектора дипольного момента. Если подход с использованием
Измерение. Мониторинг. Управление. Контроль
ЭЭГС ПТ ориентирован на диагностику нарушения процессов проводимости, то подход с ЭЭГС дипольного типа ориентирован на диагностику ишемии.
70
Библиографический список
1. Inverse problem of Electrocardiography / P. W. Macfarlane, A. van Oosterom, O. Pahlm, P. Kligfield, M. Janse, J. Camm // Comprehensive Electrocardiology. - 2th ed. - Springer-Verlag London Limited, 2011. - P. 1130-1280.
2. Мурашко, В. В. Электрокардиография / В. В. Мурашко, А. В. Струтынский. - Москва : МедПресс, 2019. - 360 с.
3. Титомир, Л. И. Неинвазивная электрокардиотопография / Л. И. Титомир, В. Г. Трунов, Э. А. И. Айду. - Москва : Наука, 2003. - 198 с.
4. Potyagaylo, D. Binary optimization for source localization in the inverse problem of ECG / D. Potyagay-lo, E. G.Cortes, W. H. W. Schulze, O.Dossel // Medical & Biological Engineering & Computing. -2014. - Vol. 52. - P. 717-728.
5. Verification of a new technique for non-invasive electrophysiological examination of the heart, based on solving the inverse problem of electrocardiography / A. Sh. Revishvili, V. V. Kalinin, O. S. Lyadzhina, et al // Herald of arrhythmology. - 2008. - Vol. 51. - P. 1-13.
6. Патент РФ 2651068. Способ неинвазивного определения электрофизиологических характе-ристик сердца / Бодин О. Н., Бодин А. Ю., Жихарева Г. В., Крамм М. Н., Палютина Ю. А., Стрелков Н. И., Черников А. И. - 2018.
7. Reconstruction of equivalent electrical sources on heart surface / G. V. Zhikhareva, M. N. Kramm, O. N. Bodin, R. Seepold, A. I. Chernikov, Y. A. Kupriyanova, N. A. Zhuravleva // 6th International work-conference (IWBBIO 2018, Granada, Spain, April 25-27) : proceedings. - Granada, Spain, 2018. - Part I.
8. Крамм, М. Н. Анализ влияния выбора количества электродов на карты распределения электрического потенциала на поверхностях торса и квазиэпикарда / М. Н. Крамм // Измерение. Мониторинг. Управление. Контроль. - 2019. - № 3 (33). - С. 61-68.
9. Титомир, Л. И. Математическое моделирование биоэлектрического генератора сердца / Л. И. Ти-томир, П. Кнеппо. - Москва : Наука ; Физматлит, 1999. - 447 с.
10. Баскаков, С. И. Электродинамика и распространение радиоволн / С. И. Баскаков. - Москва : ЛИБРОКОМ, 2012. - 416 с.
11. Бреббия, К. Методы граничных элементов / К. Бреббия, Ж. Теллес, Л. Вроубел. - Москва : Мир, 1987. - 524 с.
12. Ланцош, К. Практические методы прикладного анализа / К. Ланцош. - Москва : ГИФМЛ, 1961. -524 с.
13. Тихонов, А. Н. Методы решения некорректных задач / А. Н. Тихонов, В. Я. Арсенин. - Москва : Наука, 1979. - 285 с.
14. Крамм, М. Н. Регуляризация способа реконструк-ции эквивалентного электрического генератора сердца / М. Н. Крамм // Измерение. Мониторинг. Управление. Контроль. - 2018. - № 3 (25). -С. 86-92.
References
1. Macfarlane P. W., van Oosterom A., Pahlm O., Kligfield P., Janse M., Camm J. Comprehensive Electro-cardiology. 2th ed. Springer-Verlag London Limited, 2011, pp. 1130-1280.
2. Murashko V. V., Strutynskiy A. V. Elektrokardiografiya [Electrocardiography]. Moscow: MedPress, 2019, 360 p. [In Russian]
3. Titomir L. I., Trunov V. G., Aydu E. A. I. Neinvazivnaya elektrokardiotopografiya [Non-invasive electrocardiography]. Moscow: Nauka, 2003, 198 p. [In Russian]
4. Potyagaylo D., Cortes E. G., Schulze W. H. W., Dössel O. Medical & Biological Engineering & Computing. 2014, vol. 52, pp. 717-728.
5. Revishvili A. Sh., Kalinin V. V., Lyadzhina O. S. et al. Herald of arrhythmology. 2008, vol. 51, pp. 1-13.
6. Patent RF 2651068. Sposob neinvazivnogo opredeleniya elektrofiziologicheskikh kharakteristik serdtsa [Patent RF 2651068. Method for noninvasive determination of electrophysiological characteristics of the heart]. Bodin O. N., Bodin A. Yu., Zhikhareva G. V., Kramm M. N., Palyutina Yu. A., Strelkov N. I., Chernikov A. I. 2018. [In Russian]
7. Zhikhareva G. V., Kramm M. N., Bodin O. N., Seepold R., Chernikov A. I., Kupriyanova Y. A., Zhuravleva N. A. 6th International work-conference (IWBBIO 2018, Granada, Spain, April 25-27): proceedings. Granada, Spain, 2018, part I.
8. Kramm M. N. Izmerenie. Monitoring. Upravlenie. Kontrol' [Measurement. Monitoring. Management. Control]. 2019, no. 3 (33), pp. 61-68. [In Russian]
9. Titomir L. I., Kneppo P. Matematicheskoe modelirovanie bioelektricheskogo generatora serdtsa [Mathematical modeling of a bioelectric heart generator]. Moscow: Nauka; Fizmatlit, 1999, 447 p. [In Russian]
10. Baskakov S. I. Elektrodinamika i rasprostranenie radiovoln [Electrodynamics and radio wave propagation]. Moscow: LIBROKOM, 2012, 416 p. [In Russian]
11. Brebbiya K., Telles Zh., Vroubel L. Metody granichnykh elementov [Boundary element methods]. Moscow: Mir, 1987, 524 p. [In Russian]
12. Lantsosh K. Prakticheskie metody prikladnogo analiza [Practical methods of applied analysis]. Moscow: GIFML, 1961, 524 p. [In Russian]
13. Tikhonov A. N., Arsenin V. Ya. Metody resheniya nekorrektnykh zadach [Methods for solving incorrect problems]. Moscow: Nauka, 1979, 285 p. [In Russian]
14. Kramm M. N. Izmerenie. Monitoring. Upravlenie. Kontrol' [Measurement. Monitoring. Management. Control]. 2018, no. 3 (25), pp. 86-92. [In Russian]
Крамм Михаил Николаевич
кандидат технических наук, доцент, кафедра основ радиотехники, Национальный исследовательский университет «МЭИ»
(Россия, г. Москва, ул. Красноказарменная, 14) E-mail: KrammMN @mail.ru
Kramm Mikhail Nikolayevich
candidate of technical sciences, associate professor, sub-department of radio engineering fundamentals, National Research University «MPEI» (14 Krasnokazarmennaya street, Moscow, Russia)
Образец цитирования:
Крамм, М. Н. Эквивалентный электрический генератор сердца для неинвазивной электрокарди-одиагностики / М. Н. Крамм // Измерение. Мониторинг. Управление. Контроль. - 2020. - № 2 (32). -С. 65-71. - БО! 10.21685/2307-5538-2020-2-8.
