Математическое, алгоритмическое и программное обеспечение для системы внутрисердечной навигации

Размещено на http://allbest.ru

На правах рукописи

АВТОРЕФЕРАТ

диссертации на соискание ученой степени

кандидата технических наук

Математическое, алгоритмическое и программное обеспечение для системы внутрисердечной навигации

Специальность 05.13.18 - Математическое моделирование, численные методы и комплексы программ

Жарый Сергей Викторович

Томск - 2011

Работа выполнена в Томском государственном университете систем управления и радиоэлектроники (ТУСУР)

Научный руководитель - доктор технических наук профессор Шелупанов Александр Александрович

Официальные оппоненты: доктор физико-математических наук профессор Боровской Игорь Георгиевич (ТУСУР)

доктор технических наук профессор Кочегуров Владимир Александрович (Национальный исследовательский Томский политехнический университет)

Ведущая организация -Алтайский государственный университет (г. Барнаул)

Защита состоится «17» ноября 2011 г. в 15.15 на заседании диссертационного совета Д 212.268.02 при Томском государственном университете систем управления и радиоэлектроники по адресу: 634050, г. Томск, пр. Ленина, 40, ауд. 203

С диссертацией можно ознакомиться в библиотеке Томского государственного университета систем управления и радиоэлектроники по адресу:

634045, г. Томск, ул. Вершинина, 74.

Автореферат разослан «__» ________ 2011 г.

Ученый секретарь

диссертационного совета Мещеряков Р.В.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы. Системы навигации внутрисердечных электродов, предназначенные для отображения позиций эндокардиальных электродов при проведении малоинвазивных операций по устранению аритмий сердца, активно совершенствуются и обновляются. Используемый при лечении метод радиочастотной катетерной аблации требует выполнения электрофизиологического исследования с целью обнаружения точного положения источника аритмии.

Поиск источника проводится путём манипуляций под рентгеноскопическим контролем положением хирургического инструмента (катетера) внутри сердца, регистрации внутрисердечных электрокардиограмм и определении времени активации. Обнаруженный очаг аритмии инактивируется в процессе локального радиочастотного воздействия. Важной составляющей процедуры электрофизиологического исследования являются средства контроля и запоминания позиций инструмента (электродов) с анализируемыми параметрами.

В последнее время значительные усилия ученых (Shlomo Ben-Haim, Daniel Reisfeld, Wittkampf H.M. Frederik, Федотов Н.М., Андреев С.Ю.) и медицинских компаний (Biosense Webster, Saint Jude Medical, Boston Scientific, Электропульс) направлены на создание безопасных для пациента средств нерентгеноскопического контроля и навигации инструмента.

Несмотря на наличие средств и математических методов обработки данных в этой области исследований, существует ряд проблем, связанных с развитием теории и практической реализацией методов и средств.

Общим недостатком известных технических решений является невозможность комплексирования информации навигационных систем с регистрирующими и аблационными системами, и диагностическими электрокардиостимуляторами при формировании операционных электрофизиологических комплексов. Нет данных по распределению электрических полей, создаваемых в теле пациента для определения координат электродов, нет программных средств для исследования влияния размеров электродов и их положения на теле пациента на форму электропотенциального поля в сердце пациента.

Таким образом, моделирование электрических полей в теле человека для выбора оптимальных размеров генерирующих электродов, разработка численных методов построения объёмных моделей сердца и создание единого аппаратно-программного комплекса для выполнения операций для лечения нарушений ритма сердца с полной синхронизацией запоминаемой информации, для упрощения взаимодействия с электрофизиологическим комплексом и обеспечивающего доступ основной части населения РФ к высокотехнологичной медицинской помощи является актуальным.

Целью данной работы заключается в создании математического, алгоритмического и программного обеспечения электрофизиологического комплекса трёхмерной навигации эндокардиальных электродов для повышения качества работы хирургов при лечении нарушений ритма сердца и увеличения доступности операций катетерных аблаций нарушений ритма сердца.

Для достижения поставленной цели требовалось решить следующие задачи:

1. проанализировать существующие подходы к решению задачи моделирования электрических полей в теле человека;

2. провести моделирование электрических полей в теле человека для оценки искажений и выбора размера генерирующих электродов;

3. провести анализ существующих методов построения трёхмерных геометрических моделей;

4. выбрать или усовершенствовать один из методов построения трёхмерных геометрических моделей камер сердца;

5. создать алгоритмы обработки поступающей информации и отображения результатов работы;

6. разработать программное обеспечение для электрофизиологического комплекса, использующее созданные алгоритмы.

Объектом исследования является комплекс эндокардиального картирования, предназначенный для лечения нарушений ритма сердца.

Предметом исследования является метод построения имитационных трёхмерных моделей камер сердца, определение зависимости неоднородности электрических полей от размеров генерирующих электродов.

Методы исследований математическое моделирование, численные итерационные методы решения систем дифференциальных уравнений и эмпирические методы исследований.

Достоверность результатов производимых расчётов подтверждается сравнением решения модельной задачи с решением, полученным другими авторами, и сравнение с аналитическим решением, а так же экспериментальными данными, полученными при внедрении разработанных изделий медицинского назначения в аккредитованных организациях.

Научная новизна работы состоит в следующем:

1. Модифицирован метода Шепарда применительно к случаю аппроксимации трёхмерной замкнутой поверхности по ограниченному набору точек.

2. Усовершенствован алгоритм построения имитационных моделей камер сердца. Предложенное усовершенствование алгоритма позволило уменьшить число итераций аппроксимации до 3 раза и понизить плотность точек аппроксимирующей поверхности на 8%, исключить влияние геометрического размера набора точек и изменения последовательности добавления точек в набор.

3. Создан новый алгоритм автоматического построения трёхмерных моделей камер сердца во время манипуляции эндокардиальным электродом.

4. Получены новые результаты при моделировании электрических полей в теле человека от электродов расположенных на грудной клетке. Установлено, что при существенном изменении площади генерирующих электродов (в 36 раз), нелинейность поля изменяется незначительно (от 0,1% до 0,7%).

Практическая значимость состоит в созданном комплексе программ обеспечении для электрофизиологической системы, позволяющем проводить малоинвазивнные хирургические вмешательства по устранению нарушений ритма сердца, методом катетерной аблации, использую разработанные алгоритмы для построения объёмных поверхностей.

Применение программного комплекса, реализующего предложенные в диссертационной работе подходы, позволило решить ряд практических задач, а именно:

· сократить время проведения «рентгеноскопии» во время проведения операций катетерных аблаций нарушений ритма сердца до 35%;

· упростить восприятие текущего положения эндокардиальных электродов, что повысило информированность оперирующего хирурга;

· уменьшить время необходимое для картирования и устранения эктопических очагов в предсердиях в среднем на 15%.

Основные теоретические результаты получены при выполнении исследований в рамках проектов финансируемых по гранту Минобрнауки: «Моделирование процесса радиочастотной аблации с охлаждаемым электродом» (2005г.).

По грантам РФФИ: «Создание методологии визуализации электродов-катетеров и исследование параметров электрической активности сердца» (2006-2008г.); «Методология поиска источников патологической активности миокарда в лечебно-диагностических медицинских системах» (2008-2010г.); «Создание анатомически точных моделей камер сердца методами ротационной рентгенографии и их интеграция в медицинские системы локации и управления внутрисердечными электродами» (2009-2011г.). По конкурсу молодых учёных Администрации Томской области: «Разработка роботизированного комплекса для диагностики и лечения сложных форм аритмий сердца» (2008-2009 г.).

Результаты работы использованы при выполнении прикладных хоздоговорных НИР, финансируемых фирмой ООО «Л.М.Э. «Биоток» и Администрацией Томской области: «Устройство контроля состояния сердечнососудистой системы человека» (2005-2006г.); «Разработка комплекса трехмерной локации электродов эндокардиальных катетеров» (2006-2007г.).

Положения, выносимые на защиту

1.Модифицированный метод Шепарда объёмной интерполяции поверхности камер сердца по данным измерений электропотенциальных полей в теле человека. Модификация за счет использования адаптивного степенного коэффициента, позволила уменьшить число итераций аппроксимации более чем в 3 раза, при повышении равномерности получаемой аппроксимирующей сетки за счёт нормализации входных данных, исключения влияния порядка добавления точек в набор, создания исходной сетки с равномерным распределением узлов.

2. Алгоритм автоматизированного набора восполняющих точек позволяет построить объемную модель поверхности камеры сердца, при проведении эндоскопических вмешательств. Построение трёхмерной модели сердца связано с добавлением восполняющих точек в изменяемый набор при любых перемещениях электрода хирургом внутри сердца. Алгоритм добавления основан на совместной работе алгоритмов объёмной интерполяции и методе построения выпуклой оболочки.

3.Результаты численного моделирования распределения электрического поля в теле человека от электродов, расположенных на грудной клетке показали, что преимущества от увеличения размеров генерирующих электродов нивелируются анатомическими особенностями расположения рабочей зоны.

4.Комплекс программ моделирования, интерполяции и съема сигналов с электродов позволяет повысить качество восприятия электрофизиологической информации, упростить управление электрофизиологическим комплексом, ускорить процесс воспроизведения 3D изображений камер сердца, что сократило время проведения процедуры и снизило дозу рентгеновского излучения, получаемого пациентом и медицинским персоналом. Впервые представлен аппаратно-программный комплекс, вся информация и управление которым сконцентрированы на одной рабочей станции. Это позволило снизить количество и разнородность управляющих элементов и концентрировать внимание хирурга на двух мониторах вместо 4-х традиционных.

Внедрение результатов. Результаты диссертационной работы внедрены в продукцию ООО «Л.М.Э. «Биоток», применяются в ГВОУ ВПО Сибирском государственном медицинском университете и Научно-исследовательском институте комплексных проблем сердечно-сосудистых заболеваний СО РАМН, в учебном процессе ФГБОУ ВПО Томского государственного университета систем управления и радиоэлектроники.

Апробация работы. Результаты работы докладывались и обсуждались на Втором Всероссийском съезде аритмологов (г. Москва, 2007 г., г. Москва) и IX-м Международном славянском Конгрессе по электростимуляции и клинической электрофизиологии сердца (г. Санкт-Петербург, 2010 г.), Международной научно-практической конференция «Электронные средства и системы управления» (г. Томск, 2004 г.), Научной сессии ТУСУР (г. Томск, 2008-2010 гг.). По результатам работы опубликовано 13 статей, из них 5 в журналах рекомендованных ВАК, получено два свидетельства о государственной регистрации программ для ЭВМ в Роспатенте (№ 2010610591 «Программа навигации внутрисердечных электродов и поиска источников аритмий в трёхмерной модели сердца» и №2010615380 «Программа управления компьютерной системой визуализации рентгенографических исследований и трёхмерной реконструкции камер сердца»), получен патент на изобретение №2422084.

Личный вклад автора:

1. Предложена идея нормировки входных данных и использования сетки с равномерным распределением узлов для метода построения моделей камер сердца.

2. Получены формулы численного решения уравнения Пуассона для электрических полей в среде с неоднородной проводимостью с использованием метода конечных разностей.

3. Проведено моделирование электрических полей в теле человека от электродов, расположенных на грудной клетке, сделаны выводы об отсутствии практического преимущества от использования генерирующих электродов с большой контактной поверхностью.

4. Предложен алгоритм, сокращающий время построения моделей камер сердца.

5. Разработано алгоритмическое обеспечение для системы внутрисердечной навигации.

6. Создано программное обеспечение электрофизиологического комплекса Биоток Space Vision.

Структура и объём работы. Диссертационная работа изложена на 137 страницах и состоит из введения, четырёх глав с выводами, заключения, списка литературы из 105 наименований, содержит 77 рисунков и 8 таблиц.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Введение содержит обоснование актуальности темы диссертационной работы, определение цели и задач исследования, описание научной новизны и практической ценности работы, перечень научных положений, выносимых на защиту.

В первой главе рассматриваются основные методы локации эндокардиальных электродов и производиться их сравнение, рассматривается подход к математическому моделированию электрических полей, проводиться обзор методов построения камер сердца и делается вывод о необходимости модификации одного из методов.

Для навигации электрода катетера ранее предлагались несколько принципов навигации: рентгеноскопический контроль, магнитные поля, электрические поля, ультразвуковые датчики. В данной работе был сделан выбор в пользу метода использования электрических полей. Электрические поля позволяют определять координаты любых электродов предназначенных для электрофизиологических исследований, однако, вследствие неоднородности электрической проводимости тканей тела человека, поля в сердце могут иметь искажения. Оценить эти искажения можно проведя математическое моделирование.

Математическим описанием электромагнитных полей являются уравнения Максвелла. Из которых следует, что при колебаниях электрического поля в проводниках с небольшой частотой, током смещения можно пренебречь и считать электромагнитное поле квази-стационарным. Подобное допущение является стандартным в работах связанных с расчётом электромагнитных полей. В таком случае распределение потенциалов электрического поля может быть рассчитано, используя закона Ома в дифференциальной форме - уравнение 1. Взяв дивергенцию от обеих частей уравнения и считая, что в моделируемой области нет источников тока, получим уравнение 2.

	(1)
	(2)

Уравнение (2) относится к классу дифференциальных уравнений в частных производных эллиптического типа. Для полного решения задачи уравнения (2) необходимо дополнить граничными условиями.

Приближённое решение задачи Коши может быть получено такими методами как: методы стрельбы, сеточные методы и методы разложения. Основным недостатком методов стрельбы является то, что если исходное уравнение содержит быстрорастущие решения, то решение уравнения будет получено с высокой погрешностью.

Один из наиболее универсальных методов, получивший в настоящее время широкое распространение для приближённого решения уравнений математической физики, является метод конечных разностей, заключающийся в переходе от дифференциальных уравнений к линейным алгебраическим уравнениям путём замены производных конечно-разностными функциями. Суть метода разложения или метода проекций состоит в разложении неизвестной функции в функциональный ряд по выбранной из тех или иных соображений системе базисных функций. Выбор базисных функций для метода проекций в случае с неоднородной средой может быть затруднителен. Наиболее удобным в данном случае является использование метода сеток.

Для решения подобных задач целесообразно использовать итерационные методы, например метод релаксации. Непосредственный расчёт электрических полей и анализ результатов приводиться во второй главе. электрофизиологический аритмия сердце электрод

Задачи создания трёхмерных моделей по набору точек в настоящее время особенно актуальны, в связи с развитием методов получения наборов точек поверхностей реальных объектов. Для использования в системе навигации, к методу построения трёхмерной поверхности предъявляются следующие требования: алгоритм должен работать как при последовательном добавлении точек, так и с уже готовым набором точек; необходимо, чтобы поверхность проходила через все исходные точки; необходимо сглаживание модели; поверхность должна быть замкнутой; поверхность может быть вогнутой.

По ряду причин для удовлетворения предъявленным требованиям проще всего использовать модификацию метода Шепарда. Метод Шепарда имеет некоторые преимущества перед другими методами: позволяет строить вогнутые поверхности; предоставляет возможность использования произвольной исходной сетки.

Метод, используемый в системе Carto (США), имеет следующие недостатки: возможность накопления ошибки для той части точек, которая напрямую не меняется на этапе «окончательного выравнивания»; один и тот же набор точек может быть получен в результате разной последовательности действий, что приведёт к различию в координатах тех точек, которые напрямую не меняются на этапе «окончательного выравнивания»; используемая модификация метода Шепарда аппроксимирует поверхность, «плотно» стягивая точки сетки к точкам исходного набора, что ухудшает «качество» поверхности между исходными точками; алгоритм инициализации сетки создаёт не равномерную сетку, что может приводить к различной плотности узлов сетки на выходе алгоритма.

Этап окончательного выравнивания, добавленный в метод системы Carto, гарантирует прохождения полученной модели через все исходные точки. На этапе основного приближения сетки применяется вариация метода Шепарда, для расчёта вектора смещения узлов сетки используется формула (3).

	(3)

где - вектор смещения j-й точки сетки,

- коэффициент стягивания,

- степенной коэффициент расстояния между точками,

- малая величина для предотвращения деления на ноль,

- вектор между j-й точкой сетки и i-й точкой исходного набора,

- расстояние между j-й точкой сетки и i-й точкой исходного набора,

-величина, характеризующая плотность точек исходного набора, для i-й точки исходного набора, рассчитывается по формуле (4).

	(4)

где -номер точки исходного набора , для которой вычисляется плотность,

и - координаты i-й и k-й точек исходного набора точек.

В диссертационной работе предлагается модифицировать метод Шепарда и создать алгоритм со структурой, сходной с предложенной Daniel Reisfeld в патенте, который будет обладать следующими свойствами, представленными ранее: алгоритм инициализации исходной сетки должен создавать равномерно распределённый набор точек, расположенный на поверхности сферы; получаемая модель не должна зависеть от последовательности получения набора точек; уменьшение плотности точек аппроксимированной поверхности; полученная трёхмерная модель должна проходить через все точки исходного набора.

Во второй главе описывается численный расчёт распределения потенциалов электрического поля в неоднородной среде, приводиться описание информации об электрической проводимости тела человека и результатов моделирования.

Для перехода к конечно разностной схеме воспользуемся формулой преобразования производной произведения для уравнения 2:

	(5)

Воспользовавшись формулами аппроксимации первой (6) и второй производных (7), разрешив это уравнение относительно и применив метод релаксации для решения системы уравнений, было получено уравнение для трёхмерного случая.

	(6)
	(7)

Для получения исходной информации о проводимости среды были использованы данные компьютерной томографии.

Ткани были разделены на группы, для которых была установлена табличная проводимость. Использовались граничные условия второго рода, на границе моделируемой области производная потенциала считалась равной 0.

Важными результатами являлись оценка линейности изменения потенциалов, формы эквипотенциальных поверхностей и виртуальная форма простого геометрического объекта, размещённого в рабочей области.

При разработке навигационной системы предполагалось, что для генерации навигационного поля на грудную клетку могут быть использованы стандартные одноразовые электроды для снятия ЭКГ, диаметр такого электрода 5 мм.

В виду их небольшого размера создаваемые такими электродами поля могут иметь выраженное искривление в зоне интереса, поэтому было проведено моделирование генерирующих электродов с большими контактными поверхностями.

Сравнения полей от электродов разной величины выполнялось с помощью оценки коэффициента корреляции между распределением потенциала вдоль оси и его линейной аппроксимацией (таблица 1).

Таблица 1 - Изменения коэффициента корреляции между изменением потенциала в рабочей области и линейной аппроксимацией

размер электрода, мм2	коэффициент корреляция для оси x	коэффициент корреляция для оси y
49,92Ч51	0,9992	0,9913
99,84Ч99	0,9993	0,9963
149,76Ч150	0,9998	0,9985
199,68Ч198	0,9999	0,999
299,52Ч300	0,9999	0,9992

Примеры эквипотенциальных поверхностей для электродов размером 5Ч5см² и 30Ч30см² приведены на рисунке 1. Как видно из таблицы 1 увеличение коэффициента корреляции при увеличении размера электродов незначительно. Помимо сравнения линейности распределения потенциалов проводилось сравнение виртуальных реконструкций прямоугольника, расположенного в левом предсердии, для электродов разной величины.

Все виртуальные реконструкции переводились к новому базису неортогональных векторов для исключения влияние линейных изменений.

Сравнение взаимного пересечения площадей виртуальных реконструкций не выявило преимуществ от использования генерирующих электродов какого-либо размера.

Одной из проблем при использовании навигационной системы является изменение сопротивления контактной поверхности, в связи с высыханием токопроводящего проводящего геля.

а)
б)

Рис. 1. Форма эквипотенциальных поверхностей электрических полей, генерируемых электродами разной величины:

а) электроды 49,92Ч51 мм²; б)299,52Ч300 мм²

Моделирование изменения сопротивления контактной поверхности в 5,10,25 и 50 раз показало, что максимальное изменение потенциалов для стандартных электродов в направлении слева-направо составило 0,098, для направления грудь-спина 0,071, для электрода наибольшей контактной поверхностью в направлении слева-направо 0,035, в направлении грудь-спина 0,044 (для нормированных значений). Таким образом, можно говорить о том, что использование электродов с большой контактной поверхностью позволит уменьшить влияние изменения сопротивления контактной поверхности. Исключить влияние изменения сопротивления возможно с применением токовых генераторов для навигационного поля.

а)	б)
в)	г)

Рис. 2. Эквипотенциальные поверхности на фантомах от электродов. а) фантом с лёгкими, электроды размером 49,92Ч51 мм²; б) фантом без лёгких, электроды размером 49,92Ч51 мм²; в) фантом с лёгкими, электроды размером 149,76Ч150 мм²; г) фантом без лёгких, электроды размером 149,76Ч150 мм²

В работе было проведено математическое моделирование на фантомах. Результаты моделирования приведены на рисунке 2. Анализ приведённых рисунков позволяет сделать вывод о том, что «выпрямление» эквипотенциальных поверхностей происходит в области сердца. Таким образом, «выпрямление» эквипотенциальных поверхностей в рабочей зоне слабо зависит от размеров генерирующих электродов и обусловлено анатомическим строением грудной клетки человека.

Таким образом, результаты моделирования электрических полей от электродов с различной площадью показали, что увеличение контактной поверхности не значительно влияют на эксплуатационные качества электрического поля, при этом, использование электродов с большой контактной поверхностью позволяет уменьшить влияние изменения сопротивления контактной поверхности электродов. Рекомендуется применять электроды предназначенные, для длительного использования на пациенте. Использование стабилизации по току позволит уменьшить влияние изменения сопротивления контактных поверхностей электродов на генерируемые поля.

В третьей главе был рассмотрен комплекс программ электрофизиологического комплекса Биоток Space Vision. В состав комплекса входят следующие системы: электрофизиологическая система Биоток-1000; навигационная система Биоток-3D; система электрокардиостимуляции Биоток ЭКСД-01Л; система эндокардиальной аблации Биоток-50ЭД. Для работы комплекса были разработаны следующие программы: программа контроля радиочастотного (РЧ) деструктора; программа следящего монитора системы для электрофизиологических исследований (ЭФИ); программа блока контроля и управления ЭФИ системы и навигационной системы.

Хотя программное обеспечение ЭФИ системы и навигации представляет собой единое целое - большую часть алгоритмов можно рассматривать отдельно для ЭФИ системы и навигационной системы.

Отдельно рассмотрен метод построения трёхмерных моделей камер сердца. Для построения виртуальной геометрической поверхности камеры сердца, предлагается изменить метод, использованный в системе Carto XP (США), следующим образом: изменить метод генерации сфероида для равномерного распределения точек сфероида, применить предварительное нормирование входных координат, убрать этап искажения поверхности, изменить этап приближения сетки и отказаться от использования данных с предыдущих шагов накопления данных.

Для генерации равномерно распределённых точек сфероида может быть использован алгоритм построения геодезической сферы, в котором первым приближением к сфероиду считается икосаэдр, а последующее увеличение числа точек производиться делением сторон треугольников пополам.

Известно, что метод Шепарда чувствителен к размерам исходного набора точек. Так если размеры объекта меньше 1, смещение узлов сетки может быть и вовсе незаметным. Для предотвращения этого явления в алгоритме построения трёхмерной поверхности используется нормирование входных данных.

Основное изменение метода построения трёхмерных поверхностей заключается в модификации метода Шепарда для аппроксимации сетки к исходному набору точек. Для этого параметр k в формуле (3) заменяется адаптивным коэффициентом:

,	(9)

где - коэффициент плотности i-й точки исходного набора.

На рисунке 3 изображена работа модификации метода Шепарда, применяемой в Carto XP с коэффициентами стягивания 0,8 и 0,5, и предлагаемой модификации метода Шепарда с коэффициентом стягивания 0,5.

а)	б)
в)

Рис. 3. Сравнение результатов работы оригинального алгоритма и модифицированного алгоритмов:

а - две итерации оригинального алгоритма с коэффициентом Cf=0,8; б - три итерации оригинального алгоритма с коэффициентом Cf=0,5; в - две итерации модифицированного алгоритма с коэффициентом Cf=0,5

Как видно из рисунка плотность точек аппроксимирующей сетки в предлагаемом методе меньше. Это позволяет более качественно сглаживать поверхность камер сердца, либо при сравнимом качестве сглаживания существенно (более чем в 3 раза) уменьшать число итераций стягивания.

Обычно процесс построения модели камеры сердца предполагает последовательное перемещение электрода, установка его в нужные позиции и запоминание точек. Этот процесс можно упростить - сделать его автоматизированным.

Для определения возможности добавления точки в модель камеры сердца предлагается использовать алгоритм построения выпуклой оболочки. Если точка входит в построенную выпуклую оболочку, то эта точка добавляется в модель камеры сердца, если же точка не входит в выпуклую оболочку, то эта точка удаляется из камеры сердца.

Предложенный алгоритм предполагает, что точки, добавленные вручную не должны удаляться, независимо от вхождения их в выпуклую оболочку.

Для отображения полученных данных используется библиотека OpenGL и язык описания шейдеров GLSL. Для системы навигации был разработан ряд алгоритмов отображения модели. Также в состав комплекса входят программа контроля радиочастотного деструктора и программа следящего монитора ЭФИ.

На рисунке 4 представлен скриншот программы навигационной системы.

Таким образом, был разработан комплекс программ и алгоритмы построения объёмных моделей камер сердца, обеспечивающие функциональность электрофизиологического комплекса, создан уникальный наборы инструментов для управления отображением объёмных моделей камер сердца.

Новый алгоритм построения трёхмерных поверхностей позволяет уменьшить число итераций до 3 раз, по сравнению с аналогом.



Рис 4. Скриншот программы навигационной системы. В окне навигационной системы изображена модель камеры правого предсердия и коронарного синуса, а также плоскость среза модели

В четвёртой главе рассматриваются вопросы достоверности результатов математического моделирования электрических полей в теле человека и клинической апробации электрофизиологического комплекса.

Для подтверждения правильности программной реализации расчётов в среде Delphi было проведено сравнение c результатами, полученными Поршневым С. В. в среде MathCAD. Сравнение показало полное совпадение полученных результатов. Также было проведено сравнение с аналитическим расчётом изменения потенциалов в линейном проводнике в соответствии с законом Ома. Погрешность численного решения составила менее 0,1%. Полученные результаты позволили сделать вывод о правильности работы алгоритмов, созданных для численного моделирования электрических полей в среде с неоднородной проводимостью.

Для подтверждения работоспособности электрофизиологического комплекса проводилась клиническая апробация в отделении сердечнососудистой хирургии клиник Сибирского государственного медицинского университета. Было проведено 120 операций катетерных аблаций нарушений ритма сердца, подробные анонимизированные протоколы 5 операций приведены в приложении. Результатом всех проведённых вмешательств было устранение нарушения ритма сердца. Трёхмерная реконструкция, проведённая во время одного из вмешательств, приведена на рисунке 5.

Рис. 5. Аблационные точки в области наиболее ранней активации, располагающейся на задней стенке левого предсердия: 1 - аблационный электрод; 2 - многополюсный электрод в коронарном синусе; 3 - область наиболее ранней активации (красный цвет на модели); 4 - красные пятна - области аблационных воздействии

Работоспособность системы также подтверждается актами внедрения из лечебных учреждений (ГВОУ ВПО Сибирский государственный медицинский университет и Научно-исследовательский институт комплексных проблем сердечно-сосудистых заболеваний СО РАМН).Используемые стандартные подходы для математического моделировании электрических полей в теле человека, результаты решения модельной задачи совпадают с ранее опубликованными Поршневым С.В. Сравнение с аналитическим решением задачи распределения потенциалов в одномерном проводнике подтвердило достоверность использованных при моделировании подходов. Разработанный комплекс программ внедрён в компании ООО «Л.М.Э. «Биоток» и выпускается серийно. Точность возвращения в ранее запомненную позицию электрода составляет менее 1 мм.

В заключении излагаются основные результаты теоретических исследований и практических разработок, представленных в диссертационной работе.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

1. Исследованы существующие подходы к решению задачи математического моделирования электрических полей в теле человека. Описаны условия численного эксперимента, указаны граничные условия.

2. Проведено моделирование электрических полей в теле человека для оценки искажений и выбора размера генерирующих электродов. Выявлено незначительное улучшение «линейности» распределения полей в рабочей зоне при увеличении размера электродов. Выявлено отсутствие изменения виртуальных реконструкций прямоугольника, расположенного в рабочей области, при переходе для каждого поля к новому базису не ортогональных векторов.

3. Сформулированы требования к методу построения трёхмерных моделей камер сердца. Проведён анализ существующих методов построения трёхмерных поверхностей. Исходя из сформулированных требований из первой главы и описанных преимуществ одного из методов, решено модифицировать алгоритм на основе метода Шепарда.

4. Представлена модификация метода аппроксимации трёхмерных поверхностей и новые алгоритмы (построения трёхмерных моделей камер сердца и автоматического выбора точек для построения трёхмерных моделей камер сердца).

5. Созданы алгоритмы обработки поступающей информации и отображения результатов работы.

6. Разработано программное обеспечение для электрофизиологического комплекса, использующее созданные алгоритмы.

7. Проведена клиническая апробация электрофизиологического комплекса. Экспериментальные исследования показали возможность эффективного использования представленного решения в клинической практике. Получены акты внедрения из ГВОУ ВПО Сибирского государственного медицинского университета и Научно-исследовательского института комплексных проблем сердечно-сосудистых заболеваний СО РАМН.

8. Результаты диссертационного исследования используются в учебном процессе, что подтверждается актом о внедрении в учебный процесс ФГБОУ ВПО Томского государственного университета систем управления и радиоэлектроники.

ОСНОВНЫЕ ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ

В журналах, рекомендованных ВАК:

1. Жарый, С.В. Методы построения и визуализации трехмерных моделей по конечному множеству точек для систем внутрисердечной навигации / Н.М. Федотов, С.В. Жарый, А.И. Оферкин, А.А. Шелупанов // Вестник Сибирского государственного аэрокосмического университета. - 2006 г. - С. 24-28.

2. Жарый, С.В. Безопасные режимы радиочастотной аблации источников аритмий сердца / Н.М. Федотов, А.А. Шелупанов, А.И. Оферкин, А.С. Коблош, С.В. Жарый // Вестник Сибирского государственного аэрокосмического университета. - 2006 г. - С. 28-35.

3. Жарый, С.В. Методы повышения достоверности определения координат в системах локации и визуализации эндокардиальных электродов / Н.М. Федотов, А.А. Шелупанов, С.В. Жарый, А.И. Оферкин // Медицинская техника. - 2007. - №4. - С. 3-7.

4. Жарый, С.В. Программно-технический комплекс для исследования характеристик проводящей системы сердца / Н.М. Федотов, А.И. Оферкин, С.В. Жарый // Медицинская техника. - 2008. - №1. - С. 27-33

5. Жарый, С.В. Имитационное моделирование и визуализация электрической активности миокарда / Н.М. Федотов, С.В. Жарый, А.А. Шелупанов, А.И. Петш, А.С. Коблош // Известия Томского политехнического университета. - 2009. - Т. 315. - № 5. - C. 98-104.

Другие публикации:

6. Жарый, С.В. Новые функции генераторов для радиочастотной катетерной аблации / Н.М. Федотов, С.В. Жарый, А.С. Коблош, В.А. Бейнарович // Материалы Международной научно-практической конференции «Электронные средства и системы управления».- Ч. 3. - Томск: Изд. Института оптики атмосферы СО РАН. - 2004. - С. 66-69.

7. Жарый, С.В. Расширение функциональных возможностей радиочастотного деструктора / Н.М. Федотов, С.В. Жарый, В.А. Бейнарович // «Интеллектуальные системы в управлении, проектировании и образовании». -Выпуск 4/Под ред. Шелупанова. - Томск: Изд. Института оптики атмосферы СО РАН. - 2005 г. - С. 44-50.

8. Жарый, С.В. Коррекция искажений при построении трёхмерной модели сердца / Н.М. Федотов, С.В. Жарый // «Интеллектуальные системы в управлении, проектировании и образовании». - Выпуск 5/Под ред. Шелупанова. Томск: Изд. Института оптики атмосферы СО РАН. - 2006 г. - С. 28-34.

9. Zharyi, S.V. Methods for increasing the reliability of coordinate determination by the location and imaging systems of endocardial electrodes / N.M. Fedotov, A.A. Shelupanov, S.V. Zharyi and A.I. Oferkin // Biomedical Engineering. - Springer US. - 2007. - vol. 41. - no. 4. - P. 145-149

10. Жарый, С.В. Использование частотного картирования при аблации фибрилляции предсердий / А.И. Офёркин, Н.М. Федотов, С.В. Жарый и др. // Анналы аритмологии. - 2007. - № 3., - С. 57-57.

11. Zharyi S.V. Hardware-software system for monitoring the characteristics of the cardiac conduction system / N.M. Fedotov, A.I. Oferkin and S.V. Zharyi. // Biomedical Engineering. - Springer US. - 2008. - vol. 42. - no. 1. - P. 26-31

12. Жарый, С.В. Трехмерная ротационная атрио-вентрикулография при катетерной аблации аритмий / А.И. Оферкин, С.В. Жарый, Н.М. Федотов, А.И. Петш, И.В. Гущин, А.И. Буллер // Анналы аритмологии. - 2009. - № 2., - С. 137-137.

13. Жарый, С.В. Ротационная рентгенография для трёхмерной реконструкции камер сердца и сосудов. Материалы конгресса с международным участием «Кардиостим-2010» / С.В. Жарый, Н.М. Федотов, А.И. Оферкин // Вестник аримологии. - 2010. - С.186-186.

14. Жарый С.В. Свидетельство о государственной регистрации программы для ЭВМ «Программа навигации внутрисердечных электродов и поиска источников аритмий в трёхмерной модели сердца» / С.В. Жарый, Н.М. Федотов, А.И. Оферкин. - № 2010610591; заявл. 18.11.2009; опубл. 15.01.2010.

15. Жарый С.В. Свидетельство о государственной регистрации программы для ЭВМ «Программа управления компьютерной системой визуализации рентгенографических исследований и трёхмерной реконструкции сердца» / С.В. Жарый, Н.М. Федотов, А.И. Оферкин. - № 2010615380; заявл. 20.05.2010; опубл. 20.08.2010.

16. Жарый С.В. Патент на изобретение «Устройство слежения за электродами внутри тела пациента и способ его реализации» / Н.М. Федотов, С.В. Жарый, А.И. Оферкин. - № 2422084; заявл. 18.11.2009; зарег. 27.07.2011.

Размещено на Allbest.ru