Методы повышения информативности рентгеноскопических изображений при выполнении катетерных операций по лечению аритмий сердца

УДК 004.946: 616-073.759

Н.М. Федотов, С.В. Жарый, А.И. Оферкин

Методы повышения информативности рентгеноскопических изображений при выполнении катетерных операций по лечению аритмий сердца

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

Ключевые слова: рентгеноскопическое изображение, катетерная аблация, трёхмерная реконструкция.

Введение

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

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

Простым и логичным способом совмещения является наложение контуров контрастного изображения на текущее изображение. Выделение контуров предлагается выполнять полуавтоматически с использованием алгоритма Interactive Live-Wire Boundary Extraction [1]. Данный алгоритм позволяет пользователю приближённо указывать границы контрастного объекта, уточнение границ производится автоматически. В основе алгоритма Interactive Live-Wire Boundary Extraction лежит классический алгоритм поиска кратчайшего пути во взвешенном графе (алгоритм Дейкстры), где стоимость переходов определяется взвешенной суммой, получаемой из нескольких видов стандартных преобразований (Laplacian Zero-Crossing, фильтр Собеля и др.) - формула (1).

i ^ q)=ZfZ (q)+GfG (q)+DfD ^ q),

где p и q - два соседних пикселя изображения; fZ - значение функции Laplacian Zero-Crossing; fG - значение магнитуды градиента (фильтр Собеля); D - направление локального градиента; z ,g ,raD - веса коэффициентов, в оригинальной статье приводятся значения

roZ = 0,43, roG = 0,43, roD = 0,14.

Результаты работы алгоритма Дейкстры могут повторно использоваться на последующих шагах, при изменении целевой точки. Алгоритм предполагает ручное либо автоматическое (path cooling) закрепление сегментов пути.

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

а б

Рис. 1. Варианты наложения контура, полученного с помощью алгоритма на текущее изображение: а - изменение цвета области камеры; б - наложение контура камеры

Более удобным вариантом является получение трёхмерной модели из результатов компьютерной томографии с последующим наложением полученной картинки на текущее рентгеноизображение. Для сегментации изображений можно использовать алгоритмы, основанные на алгоритме Interactive Live-Wire Boundary Extraction или GraphCut [2].

Для трёхмерной реконструкции модели из сегментированных изображений можно использовать алгоритм марширующих кубов [3] (рис. 2). Алгоритм марширующих кубов применяется для отображения воксельных данных, является развитием ранее существовавшего алгоритма марширующих квадратов. Алгоритм марширующих кубов предполагает отображение набора треугольников для каждого вокселя в соответствии с имеющимися для данного вокселя «соседями», таким образом, может существовать всего 256 вариантов наборов треугольников, в зависимости от имеющихся или отсутствующих восьми соседей. Заранее вычисленный массив из 256 конфигураций полигонов можно получить поворотами и отражениями 15 различных конфигураций куба. Алгоритм известен тем, что он, несмотря на его очевидность, был запатентован [4]. Для совмещения изображений можно, как и в предыдущем варианте, использовать рентген контрастную метку, наложенную непосредственно перед томографическим исследованием. Предложенный метод имеет один серьёзный недостаток - томографические данные нельзя получить непосредственно во время операции на сердце.

Нашим коллективом ранее был предложен метод интраоперационной трёхмерной реконструкцией камер сердца по рентген изображениям [5]. Суть метода заключается в последовательном накоплении небольшого числа точек поверхности - порядка нескольких десятков точек. Для определения точек поверхности камеры сердца используется метод определения координат эндокардиальных электродов внутри сердца по двум рентгеновским проекциям (Рис. 3). После накопления некоторого набора точек можно применить метод трёхмерной реконструкции поверхности, который более подробно описан в статье [6] и является развитием метода степеней обратных расстояний. Недостатком метода рентгенореконструкции по набору точек является необходимость длительного накопления координат электродов и связанное с ним длительное время рентгеноизлучения.

120.0 kV 252.0 mA 1.3 mm/0.0:1 Tilt: 0.0 0.4 s

DFOV: 0.0 x 0.0cm W:860 L:-101 а б

Рис. 2. Пример трёхмерной реконструкции из томографических срезов: а - реконструкция левого предсердия; б - общая реконструкция грудной клетки

2 а

За

Ф"

1а

-f—-

!

26

3 Ъ

1Ъ

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

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

- разобьем двухмерные проекции на горизонтальные слои;

- для набора линейных проекций каждого слоя построим обратные проекции на двухмерном срезе;

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

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

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

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

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

Рис. 5. Пример трёхмерной реконструкции левого предсердия по 32 проекциям

Все предложенные методы были реализованы в программном обеспечении специализированной рентгенотелевизионной установки для малоинвазивной хирургии нарушений ритма сердца XR Биоток-15 [8].

Литература

1. Barrett W.A. Interactive Live-Wire Boundary Extraction / W.A. Barrett, E.N. Mor-tensen // Medical Image Analysis. - 1997. - Vol. 1, № 4. - P. 331-341.

2. Boykov Y. Interactive organ segmentation using graph cuts / Y. Boykov, M.-P. Jolly // In Medical Image Computing and Computer-Assisted Intervention. Proceedings of «MICCAI». - Pittsburgh, 2000. - P. 276-286.

3. Lorensen W.E. Marching Cubes: A high resolution 3D surface construction algorithm / W.E. Lorensen, H.E. Cline// ACM SIGGRAPH Computer Graphics. - 1987. -Vol. 21, № 4. - P. 163-169.

4. Cline H.E. System and method for the display of surface structures contained within the interior region of a solid body / H.E. Cline, W.E. Lorensen // US Patent 4,710,876. Dec. 1, 1987.

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

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

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

8. Жарый С.В. Программа для управления системой визуализации рентгенографических исследований и трёхмерной реконструкции сердца / С.В. Жарый, Н.М. Федотов, А.И. Оферкин. - Регистрационный номер 2010615380 от 20 августа 2010 г.

Федотов Николай Михайлович

Канд. техн. наук, зав. лаб. безопасных биомедицинских технологий ЦТБ ТУСУРа

Тел.: +7-913-829-27-53

Эл. почта: fmobile@rambler.ru

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

Аспирант каф. комплексной информационной безопасности

и электронно-вычислительных систем ТУСУРа

Тел.: (382-2) 56-09-23

Эл. почта: sergey.zhary@gmail.com

Оферкин Александр Иванович

Канд. мед. наук, зав. отделением сердечно-сосудистой хирургии СГМУ

Тел.: (382-2) 55-59-91

Эл. почта: oferkin@biotok.ru

Fedotov N.M., Zhary S.V., Oferkin A.I.

Methods for enhancement of information ability of fluoroscopic images while performing catheter operations for treatment for cardiac arrhythmias

The methods for creating three-dimensional objects from multi-plane x-ray survey,which is performed in high radiation mode in order to improve the information content of the fluoroscopic images while performing catheter operations for treatment for cardiac arrhythmias,are presented. Keywords: X-ray image, catheter ablation, three-dimensional reconstruction.