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

Савельева, Татьяна Александровна

Комбинированный спектроскопический метод исследования сильнорассеивающих биологических сред тема автореферата и диссертации по физике, 01.04.21 ВАК РФ

Савельева, Татьяна Александровна АВТОР

кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Москва МЕСТО ЗАЩИТЫ

2013 ГОД ЗАЩИТЫ

01.04.21 КОД ВАК РФ

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

Автореферат диссертации на тему "Комбинированный спектроскопический метод исследования сильнорассеивающих биологических сред"

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

005537143

Савельева Татьяна Александровна

КОМБИНИРОВАННЫЙ СПЕКТРОСКОПИЧЕСКИЙ МЕТОД ИССЛЕДОВАНИЯ СИЛЬНОРАССЕИВАЮЩИХ БИОЛОГИЧЕСКИХ СРЕД

01.04.21 — Лазерная физика

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

7 НОЯ 2013

Москва-2013

005537143

Работа выполнена в Федеральном государственном бюджетном учреждении науки Институте общей физики им. A.M. Прохорова Российской академии наук

Научный руководитель:

Официальные оппоненты:

доктор физико-математических наук, профессор Лощенов Виктор Борисович Федеральное государственное бюджетное учреждение науки Институт общей физики им. АМ. Прохорова Российской академии наук, заведующий лабораторией

доктор физико-математических наук, профессор Кустов Евгений Федорович, Национальный исследовательский университет «МЭИ», профессор

кандидат физико-математических наук, Приезжев Александр Васильевич, Московский государственный университет имени М. В. Ломоносова, доцент

Ведущая организация:

Национальный исследовательский Саратовский государственный университет имени Н.Г. Чернышевского

Защита диссертации состоится 2 декабря 2013 г. на заседании диссертационного совета Д002.063.03 при ИОФ РАН по адресу г. Москва 119991, ул. Вавилова, 38.

С диссертацией можно ознакомиться в библиотеке ИОФ РАН. Автореферат разослан

(дата)

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

диссертационного совета_¿уу Сп?-^---Т.Б. Воляк

(подпись)

Общая характеристика работы

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

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

На настоящий момент при проведении флуоресцентного спектрального анализа ¡n vivo в некоторых случаях его сочетают с измерением спектров обратного диффузного отражения [1, 2]. Коррекция спектров флуоресценции с помощью результатов анализа спектров диффузного отражения обеспечивает определение состава флуорофоров в ткани. В ряде случаев спектр обратного диффузного отражения используется для независимого анализа поглощающих и рассеивающих свойств исследуемых сред [3]. Также следует отметить, что последовательное измерение спектров диффузного отражения и флуоресценции чревато возникновением артефактов, обусловленных взаимным смещением источника света, объекта и фотоприемника, и увеличением времени анализа, что может оказаться неприемлемым в- клинических условиях. Что касается систем спектрально-разрешенной визуализации, то они, как правило, используют наборы узкополосных фильтров или перестраиваемые фильтры для получения изображений объекта, либо основаны на методе сканирования всего поля изображения поточечно с использованием классической спектроскопической

системы с диспергирующим элементом [4]. У обоих подходов основным недостатком является низкая скорость анализа данных.

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

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

Для чего были решены следующие задачи:

1. Разработка методов регистрации комбинированных спектров диффузного отражения и флуоресценции с применением волоконно-оптических зондов и матричных фотоприемников.

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

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

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

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

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

Практическая значимость

Предложенный метод комбинированной спектроскопии был апробирован в НИИ нейрохирургии им. H.H. Бурденко при проведении интраоперационной навигации во время операций по удалению внутричерепных и интрамедуллярных опухолей. В качестве опухолевого маркера использовался эндогенный протопорфирин IX, индуцированный предварительным введением в организм пациента 5-аминолевулиновой кислоты. Метод показал повышение чувствительности и специфичности диагностики по сравнению как с ординарным спектроскопическим анализом, так и с анализом видеофлуоресцентных изображений с помощью операционного флуоресцентного микроскопа Karl Zeiss Opmi Pentero.

Система анализа мультиспектральных изображений, сопряжённая со стандартной щелевой лампой, была использована в ФГЪУ «НИИ глазных болезней» РАМН дня косвенной оценки уровня внутриглазного давления по степени оксигенации сосудов глазного дна и при проведении лабораторных исследований по эффективности фотодинамической терапии и приживлению тканей в ИОФ РАН.

Разработанные методы также нашли применение в научных исследованиях, выполняемых совместно с ФГБУ «ГНЦССП им. В.П. Сербского» Минздрава России. Основные положения, выносимые на защиту

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

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

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

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

Материалы диссертации были доложены на ряде научных конференций: II Троицкая конференция «Медицинская физика и инновации в медицине»; Отечественные противоопухолевые препараты 2007, 2008, 2009; Advanced Laser Technologies 2007, 2008, 2009; International symposium on Laser medical application (LMA'2010), Москва; XIV International School for Young Scientists and Students on Optics, Laser Physics & Biophotonics Saratov Fall Meeting SFM'2010, 2011; научно-практическая конференция «Фотодинамическая терапия и флуоресцентная диагностика» 2011, Санкт-Петербург; LPHYS'll, Sarajevo; Laser Florence 2011; SPIE Photonics West, San Francisco, 2012.

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

Содержание работы

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

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

пространственным разрешением, что приобретает особое значение в микрохирургии и нейрооикологии.

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

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

Спектроскопические измерения проводились с помощью лазерного спектроанализатора ЛЭСА-01-БИОСПЕК [5], модифицированного системой фильтров и использованием волоконно-оптического зонда с двумя осветительными каналами. В диапазоне спектра 500-600 нм реализована регистрация спектров обратного диффузного отражения, позволяющая производить количественную оценку концентрации гемоглобина в оксигенированной и редуцированной форме. Интенсивность диффузно отраженного назад лазерного излучения, используемого для возбуждения флуоресценции опухолевого маркера, служит для анализа рассеивающих свойств исследуемых тканей. В красном и ближнем инфракрасном диапазонах регистрируется сигнал флуоресценции протопорфирина IX и продуктов фотораспада.

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

флуоресценцию содержащихся в ней флуорофоров, что приводит к регистрации флуоресцентного излучения.

г - расстояние от источника

Рис. 1. Схема проведения спектроскопических измерений с использованием оптоволоконной системы.

Взаимодействие излучения с биологическими средами на микроуровне описывается в терминах комплексного показателя преломления (ñ = n + i ■ к) и относительного размера (х = 2лпа/А) основных структурных элементов ткани (с характерным размером а). Эти параметры позволяют с помощью уравнений Максвелла вычислить такие параметры как сечение поглощения (аа) и рассеяния (ст5) света отдельной частицей, а также определить угловую зависимость рассеяния (Р(0)).

При переходе от отдельных частиц к их ансамблям используют уже усредненные по некоторому элементарному объёму оптические характеристики, такие как коэффициент поглощения (ца = с • аа), коэффициент рассеяния = с • as) и фактор анизотропии (g = (cos(0)>). Для сред, характеризующихся анизотропным рассеянием, вводят также так называемый редуцированный коэффициент рассеяния (|i's = ц5(1 — g)).

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

(ТПИ) в среде позволяет, последовательно пройдя по всем введенным выше масштабным описаниям, восстановить фундаментальные параметры компонентов среды. Одной из важных характеристик распространения излучения в сильнорассеивающей среде является глубина проникновения света в вещество, определяемая в диффузионном приближении ТПИ как б = 1/цеГГ = 1/

л/3Иа(Иа + С1

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

Воздух _ _ Регистрация 0Д0 в

дольней зоне

ПЗС

Равномерно освещенная поверхность

Ч^Г

У/УЛ

Регистрация 0Д0 в ближней зоне

Область наиболее вероятного присутствия детектируемых (ротонов

Рис. 2. Схема моделирования регистрации спектрально-разрешенных изображений -формирование сигнала одного пикселя изображения. Регистрации мультиспекгральных изображений также соответствует особая геометрия измерений. Для реализации метода было разработано устройство смены

фильтров, рассчитанное на комплект из четырех фильтров с учетом количества исследуемых параметров (любые четыре из набора 542 нм, 560 нм, 576 нм, 586 нм, 596 нм и фильтр для регистрации флуоресцентного излучения с максимумом пропускания на 710 нм). В качестве приемника изображений использовались камера Видеоскан-415-USB и фотоаппарат Canon EOS 5D Mark II. Рассмотрены варианты установки фильтрового колеса в осветительном и приемном каналах. На рис. 2 приведена схема формирования сигнала одного пикселя изображения как в дальней, так и в ближней зоне наблюдения. В работе рассмотрено влияние числовой апертуры приемника в различных положениях относительно поверхности ткани на значения детектируемого сигнала диффузного отражения.

В силу высокой плотности рассеивающих элементов биологические ткани принято рассматривать как мутные среды, поэтому в работе подробно рассмотрены различные приближения теории переноса излучения и влияния на результат измерений геометрии эксперимента. В случае сложной геометрии оптической системы и/или исследуемой среды зачастую используют имитационное моделирование распространения света в сильнорассеивающих средах. Авторский программный код, реализующий имитационное моделирование методом Монте-Карло, приведен в Приложении к диссертационной работе. Особенностью разработанного алгоритма является учет длины пути света в ткани от осветительного до приемного волокна и глубины зондирования. Гистограммы распределения этих величин подвергались аппроксимации логнормальной зависимостью, мода которой соответствовала использованным при поиске описательной аналитической зависимости средним значениям длины пути и глубины зондирования.

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

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

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

Алгоритм анализа комбинированных спектров включает выделение из сигнала компонент, обусловленных рассеянием и поглощением, их интерпретацию и учет их влияния на спектр флуоресценции. Спектральная зависимость коэффициента поглощения позволяет восстановить информацию о концентрации исследуемых хромофоров. Для этого производится деконволюция спектра поглощения на известные компоненты, такие как гемоглобин в оксигенированной и редуцированной форме и Пп IX. Редуцированный коэффициент рассеяния используется в качестве критерия оценки плотности тканей, а также для количественного анализа содержания флуоресцирующего маркера в рассеивающей среде по спектру флуоресценции. Спектр флуоресценции, регистрируемый в длинноволновой области видимого спектра, нормируется на интенсивность диффузно-отраженного лазерного излучения с учетом поправки на отличия рассеивающих свойств среды в диапазоне флуоресцентного и возбуждающего излучения. Предложенный алгоритм не ограничивает свою функциональность 5-АЛК индуцированным протопорфирином IX и может быть использован для анализа ряда флуоресцентных маркеров и фотосенсибилизаторов, в том числе тетрапиррольной природы, обладающих способностью к возбуждению и флуоресценции в соответствующих диапазонах спектра.

Алгоритм восстановления данных основан на аналитической аппроксимации, полученной с помощью численного моделирования зависимости коэффициента диффузного отражения от оптических свойств исследуемой среды. Ниже приведены результаты моделирования для значений числовой апертуры источника и приемника излучения ЫА==0.37 для ряда расстояний между ними 0.25, 0.5, 1, 1.5 мм (рис.3). Показатель поглощения изменяется при этом в диапазоне от 0.5 до 8 см"1, что соответствует физиологическим значениям нормального белого вещества и астроцитарных опухолей различной степени озлокачествления [6].

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

Глубина зондирования, мм

^ ♦......0.25 мм

-Ч8Ь~ 0.5 мм —А— 1 мм —л— 1.5 мм

Длина пути, мм

—*—0.25 мм -НШ—0.5 мм

^ — ¿г~ 1 ММ

1.5 мм

0 5

Коэффициент поглощения, см-1

О 5

Коэффициент поглощения, см-1

Рис. 3. Влияние коэффициента поглощения на коэффициент диффузного отражения (Яс1), глубину зондирования и длину пути света в ткани для различных расстояний между источником и приемником.

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

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

Яй(ца,с1) = ВД-е"1^ (1)

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

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

0.05

0.04

0.03

0.02

0.01

Глубина зондирования, мм

—♦— 0.25 мм

0.5 мм —к— 1 мм 1.5 мм

Длина пути, мм

♦ 0.25 мм «Ж—0.5 мм ~!г- 1 мм 1.5 мм

Лг- *г

Коэффициент рассеяния, см-1

0 500

Коэффициент рассеяния, см-1

Рис. 4. Влияние коэффициента рассеяния на коэффициент диффузного отражения (%/), глубину зондирования и длину пути света в ткани для различных расстояний между источником и приемником.

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

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

экспоненциальная зависимость, содержащая в качестве сомножителя коэффициент рассеяния, что обеспечило рост сигнала в проксимальной к источнику области: Яа = А • ц'5 ■ (2)

где А, В| -коэффициенты, зависящие от расстояния между источником и приемником излучения.

Модель показала хорошее соответствие с данными численного эксперимента (коэффициент корреляции >0.997) (рис.5).

4.Е+04 З.Е+04 З.Е+04 2.Е+04 2.Е+04 1_Е+04

5.Е+03 0.Е+00

число фотонов

5.Е+04

число фотонов

5.Е+04 4X404 «+04 З.Е+04 '"'•.. З.Е+04 2.Е-04

1.Е+04 5.Е+03

2*35_- .^Г-г.-ут .=■ о.Е+сю -

0 2 4 6 8 0

Коэффициент поглощения, см-1

♦ 0.25,чм-Монте-Карло

я 0.5мм-Монте-Карло

* 1мм- Монте-Карло

X 1.5мм-Монте-Карло ......0.25мм-------^ экспонента модиф.

---0.5мм- экспонента

моднф.

--1 мм - экспонента

~ - - А модиф. ~ ' ш '* — ■ -1.5мм-экспонента 100 модиф.

Редуцированный коэффициент рассеяния, см-1

Рис. 5. Аппроксимация зависимостей сигнала диффузного отражения от рассеяния и поглощения, полученных в результате численного моделирования.

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

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

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

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

в Монте-Карло 1 Ш in vivo

Край опухоли Центр Край опухоли Центр (B ill) опухоли Ш- (IV) опухоли (IV)

III)

Рис. 6. Распределение относительного коэффициента диффузного отражения в зависимости от морфологического статуса ткани по результатам исследования внутриопухолевых локальных изменений (II, II, IV - степень озлокачествления по ВОЗ).

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

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

Апробация комбинированного спектроскопического метода в клинических условиях осуществлялась в НИИ нейрохирургии им. Н.Н. Бурденко в рамках клинических исследований эффективности флуоресцентной диагностики с использованием 5-AJIK индуцированного Пп IX. В данной работе приведено рассмотрение распределения исследуемых параметров в операционном поле при удалении опухоли для пациентов с глиобластомой. На рис. 7 приведены спектральные зависимости, получаемые интраоперационно с помощью комбинированного метода.

1 Оксигенаиия

Гемоглобин * Рассеяние ■ Флуоресценция

500 550 600 650 700 Длина волны, ем

Рис. 7. Пример спектральных зависимостей, регистрируемых интраоперационно, и

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

нормальный мозг, 2 - край опухоли, 3 — центр опухоли.

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

анализа, проводимого in vivo в зоне удаления (резекции) опухоли, и гистологической

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

соответствии исследуемых спектроскопических характеристик типу ткани.

Спектроскопические измерения, производимые во время операции, показали

плавные изменения исследуемых характеристик в направлении от здоровой ткани к

центру опухоли, подчиняющиеся определенным закономерностям. В тканях с преобладанием некроза было показано отсутствие накопления Пп IX, высокий уровень рассеяния и кровенаполнения, низкий уровень оксигенации. Для этих участков уместнее оперировать термином «концентрация гемоглобина» вместо термина «кровенаполнение», поскольку в некротических тканях высокое значение этого параметра, скорее всего, достигается за счет геморрагий. Для образцов ткани из центра опухоли с микроскопическими включениями некроза был показан флуоресцентный контраст по сравнению со здоровой корой равный 7±1.7, уменьшение плотности ткани более чем вдвое по сравнению с нормой, кровенаполнение, превышающие в среднем в 6 раз и уровень оксигенации составляющий 60% от значений соответствующих параметров интактного мозга. Ткань из центра опухоли без признаков некроза показала сходные характеристики за исключением кровенаполнения, которое оказалось ниже. Зоне инфильтрации соответствовал максимальный уровень накопления Пп IX (в 10±3 раз выше, чем в норме), наименьшая плотность ткани, кровенаполнение вдвое превышающее норму и самый высокий уровень оксигенации, близкий к таковому интактного мозга. Избирательное накопление Пп IX в опухолевых тканях по сравнению с нормальными обусловлено рядом ферментативных нарушений в быстро пролиферирующих клетках опухоли. Также было показано [7], что накопление Пп IX коррелирует с индексом пролиферации, что объясняет высокий уровень в активно растущей части опухоли в зоне инфильтрации и несколько заниженные его значения в центре опухоли. В целом можно заключить, что более малигнизованные участки опухоли характеризуются более высоким уровнем кровенаполнения и пониженным уровнем оксигенации по сравнению с зоной инфильтрации опухолевых клеток в здоровый мозг.

Использование комбинированного спектроскопического метода привело к повышению чувствительности интраоперационной демаркации границ опухоли по сравнению с методом видеофлуоресцентной визуализации с использованием микроскопа Zeiss OPMI Pentero с 78 % до 96 %, специфичности - с 60 % до 82 %.

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

системы фильтров. Представлены различные варианты реализации системы. Один из вариантов повторяет функциональность рассмотренного в третьей главе комбинированного спектроскопического метода, объединяя анализ гемоглобина в диапазоне 500-600 нм по сигналу диффузного отражения и анализ Пп IX в окрестности 710 нм по сигналу флуоресценции. Второй вариант реализации системы подразумевает одновременное использование двух источников возбуждения флуоресценции Пп1Х с длинами волн 405 нм и 630 нм, что соответствует полосе Соре и одной из линий в О-полосе спектра поглощения Пп IX, с целью повышения эффективности возбуждения флуоресценции и преодоления экранирующего действия крови. Также в четвертой главе рассмотрены вопросы спектрального и пространственного разрешения системы визуализации, проведено имитационное моделирование взаимодействия излучения с рассеивающей средой с учетом геометрии используемой оптической системы. Показаны отличия зависимости сигнала диффузного отражения от оптических параметров для системы визуализации по сравнению с оптоволоконной системой.

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

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

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

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

А = 1п ( 1/К ) = А0 + Аг ■ Ца/ТТГ (3)

где А — коэффициент экстинкции; Л - регистрируемый сигнал диффузного отражения; А0 и А1 - подстроенные коэффициенты, зависящие от геометрии измерений.

5 6 ё 5

Монте-Карло * £ 4

А = 0.1364ца+0.37 а

=г 3

0 2 4 6

Коэффициент поглощения, см~1

-Монте-Карло -А=3/5С№Т(ц'з)

200 400 600 800

Коэффициент рассеяния, см-1

Рис. 8. Аппроксимация зависимости коэффициента экстинкции от оптических параметров

среды.

Оценка точности вычисления с помощью этой модели степени оксигенации

при аппроксимации экспериментальных данных предложенной теоретической

зависимостью по всему диапазону доступных длин волн при измерении спектра с

помощью спектроанализатора и по 4-м характерным длинам волн с учётом полос

пропускания интерференционных фильтров показала соответствие в пределах 5%.

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

Рис. 9. Визуализация степени оксигенации сосудов кожи пальцев руки (степень оксигенации представлена стандартной цветной шкалой, соответствующей равномерному разбиению диапазона степени оксигенации от 0 до 100%).

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

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

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

Разработаны методы комбинированного спектроскопического анализа с применением волоконно-оптических зондов и матричных фотоприемников, алгоритмы анализа комбинированных спектров и спектрально-разрешенных изображений. Создано программное и аппаратное обеспечение, показавшее по результатам апробации в клинике повышение чувствительности с 78% до 88%, специфичности - с 60% до 82% по сравнению с методом видеофлуоресцентной визуализации с использованием микроскопа.

Основные публикации по теме диссертации с участием автора:

1. А. А. Стратонников, Г. А. Меерович, А.В.Рябова, Т.А.Савельева, В.Б. Лощенов, "Использование спектроскопии обратного диффузного отражения света для мониторинга состояния тканей при фотодинамической терапии", Квант, электроника, 2006, 36(12), 1103-1110.

2. Савельева Т.А., Стратонников А.А., Меерович Г.А., Лощенов В.Б., «Разработка метода и оборудования для видеофлуоресцентной диагностики рака желудка», Альманах клинической медицины. 2006. № 12. С. 31-31.

3. Стратонников А.А., Кузнецов А.В., Меерович Г.А., Рябова А.В., Савельева Т.А., Лощенов В.Б., «Использование спектроскопии обратного диффузного отражения света ин виво для диагностики, контроля за терапией и в спортивной медицине», Альманах клинической медицины. 2006. № 12. С. 38-38.

4. Tatiana A. Savelieva; Aleksander A. Stratonmkov; Victor В. Loschenov; Multi-spectral imaging of oxygen saturation. Proc. SPIE 7022, Advanced Laser Technologies 2007, 702205 (June 05,2008); dol: 10.1117/12.803905.

5. Nickolay N. Brysin; Kirill G. Linkov; Alexander A. Stratonnikov; Tatiana A. Savelieva; Victor B. Loschenov; The video fluorescent device for diagnostics of cancer of human reproductive system. Proc. SPIE 7022, Advanced Laser Technologies 2007, 702206 (June 05,2008); doi:10.1117/12.803906.

6. Брысин H.H., Линьков К.Г., Стратонников A.A., Савельева Т.А., Рябова А.В., Беляева Л.А., Лощенов В.Б. Видеофлуоресцентный комплекс для диагностики репродуктивной системы человека Лазерная медицина, Т. 13, 2009 г., стр 49-53.

7. N. Kalyagina, V. Loschenov, D. Woif, С. Daul, W. Blondel, T. Savelieva, "Experimental and Monte Carlo investigation of visible diffuse reflectance imaging sensitivity to diffusing particle size changes in an optical model of a bladder wall", Applied Physics B, vol. 105(3), pp. 631-639,2011.

8. T.A. Савельева, A.B. Рябова, И.В. Андреева, H.A. Калягина, В.И. Конов, В.Б. Лощенов, Комбинированный спектроскопический метод определения концентрации флуорофоров в сильнорассеивающих средах, Краткие сообщения по физике, ноябрь 2011.

9. Т. Savelieva, M.V. Loschenov, K.G. Linkov, The multispectral visualization of glial brain tumors containing PP IX in diffuse and laser induced fluorescent light, Lasers in Medical Science, Volume 26, Supplement 1, 2011.

10. T.A. Savelieva, M.V. Loschenov, K.G. Linkov, The multispectral visualization of glial brain tumors containing PP IX in diffuse and laser induced fluorescent light, Lasers in Medical Science, Volume 26, Supplement 1,2011.

11. А.А. Потапов, А.Г. Гаврилов, В.А. Охлопков, Л.В. Шишкина, С.А. Горяйнов, В.А. Шурхай, П.В. Зеленков, В.Б. Лощенов, Т.А. Савельева, П.В. Грачев, М.Н. Холодцова, С.Г. Кузьмин, Д.А. Гольбин, Использование флюоресцентной диагностики и лазерного спектрального анализа в хирургии внутричерепных менингиом, Вопросы нейрохирургии имени Н.Н. Бурденко, т. 76, No. 4, 2012.

12. А.А. Потапов, А.Г. Гаврилов, С.А. Горяйнов, Д.А. Гольбин, П.В. Зеленков, Г.Л. Кобяков, В.А. Охлопков, В.Ю. Жуков, Л.В. Шишкина, В.А. Шурхай, В.Б. Лощенов, Т.А. Савельева, П.В. Грачев, М.Н. Холодцова, С.Г. Кузьмин, Г.Н. Ворожцов, Интраоперационная флуоресцентная диагностика и лазерная спектроскопия в хирургии глиальных опухолей головного мозга, Вопросы нейрохирургии имени Н.Н. Бурденко, т. 76, No. 5, 2012.

13. С.А. Горяйнов, А.А. Потапов, Д.А. Гольбин, П.В. Зеленков, Г.Л. Кобяков, А.Г. Гаврилов, В.А. Охлопков, В.А. Шурхай, Е.В. Шелеско, В.Ю. Жуков, В.Б. Лощенов, Т.А. Савельева, С.Г. Кузьмин, Флуоресцентная диагностика и лазерная биоспектросокпия как один из методов мультимодальной нейронавигации в нейрохирургии, Вопросы нейрохирургии имени Н.Н. Бурденко, т. 76, No. 6, 2012.

14. Потапов А.А., Горяйнов А.А., Лощенов В.Б., Савельева Т.А., Гаврилов А.Г., Охлопков В.А., Жуков В.Ю., Зеленков П.В., Гольбин Д.А., Шурхай В.А., Шишкина JI.B., Грачев П.В., Холодцова М.Н., Кузьмин С.Г., Ворожцов Г.Н., Чумакова А.П. Интраоперационная комбинированная спектроскопия (оптическая биопсия) глиом головного мозга, опубликована в журнале Вопросы нейрохирургии т. 77, № 1, 2013.

15. А.Н. Коновалов, А.А. Потапов, А.Г. Гаврилов, В.А. Шурхай, С.А. Горяйнов, В.А. Охлопков, Г.Л. Калинин, Д.В. Фомичев, М.А. Кутин, В.Ю. Жуков, Д.А. Гольбин, П.В. Зеленков, М.А. Хить, В.Б. Лощенов, Т.А. Савельева, М.Н. Холодцова, Г.Н. Ворожцов, С.Г. Кузьмин, В.П. Чехонин, Г.М Юсубалиева, В.П. Баклаушев. Глава 2. Современные технологии в нейрохирургии. // В издании «Современные технологии и клинические исследования в нейрохирургии» под ред. А.Н. Коновалова. Стр. 55-112.

16. T. Savelieva, N. Kalyagina, M. Kholodtsova, A. Potapov, S. Goryainov, V.B. Loschenov. Numerical modelling and in vivo analysis of fluorescent and laser light backscattered from glial brain tumours. // SPIE Proceedings, Vol. 8230, 82300L-1-L9, 2012.

17. Sergey S. Model; Tatiana A. Savelieva; Kirill G. Linkov; System for determining the concentration and visualization of the spatial distribution of photosensitizers based on tetrapyrrole compounds in the tissues of the human ocular fundus . Proc. SPIE 8699, Saratov Fall Meeting 2012: Optical Technologies in Biophysics and Medicine XIV; and Laser Physics and Photonics XIV, 86990E (February 26,2013); doi:10.1117/12.2016675.

Цитируемая литература

1. Robert S. Bradley and Maureen S. Thorniley, A review of attenuation correction techniques for tissue fluorescence, J. R. Soc. Interface (2006) 3, 1-13.

2. G. Palte et al. FLUORESCENCE AND REFLECTANCE SPECTROSCOPY FOR PROTOPORPHYRIN IX QUANTIFICATION IN TISSUE-LIKE MEDIA, Photonic Therapeutics and Diagnostics VII, edited by N. Kollias, et al., Proc. of SPIE Vol. 7883, 788340, 2011.

3. Pablo A.Valdes et al., Combined fluorescence and reflectance spectroscopy for in vivo quantification of cancer biomarkers in low- and high-grade glioma surgery, Journal of Biomedical Optics 16(11), 116007, November 2011.

4. Elizabeth M. C. Hillman, Optical brain imaging in vivo: techniques and applications from animal to man, J Biomed Opt. 2007; 12(5): 051402.

5. Loschenov V.B., Konov V.I., Prokhorov A.M. Photodynamic Therapy and Fluorescence Diagnostics. // Laser Physics. 2000. V. 10. No.6. P. 1188-1207.

6. Biomedical Photonics Handbook. Ed. Vo-Dinh T. Boca Raton, FL: CRC Press, 2003.

7. Valdes P.A. et al. Neuro-Oncology, 2011, 13, 8, 846-856.

Формат 60x90/16. Заказ 1727. Тираж 100 экз.

Печать офсетная. Бумага для множительных аппаратов.

Отпечатано в ООО "ФЭД+", Москва, Ленинский пр. 42, тел. (495)774-26-96

Текст научной работы диссертации и автореферата по физике, кандидата физико-математических наук, Савельева, Татьяна Александровна, Москва

Работа выполнена в лаборатории Лазерной биоспектроскопии Центра естественно-научных исследований Института общей физики им. A.M. Прохорова Российской академии наук

КОМБИНИРОВАННЫЙ СПЕКТРОСКОПИЧЕСКИЙ МЕТОД ИССЛЕДОВАНИЯ СИЛЬНОРАССЕИВАЮЩИХ БИОЛОГИЧЕСКИХ СРЕД

Диссертация на соискание ученой степени кандидата физико-математических наук

Специальность: 01.04.21 - лазерная физика

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

Савельева Татьяна Александровна

Научный руководитель: зав. лаб. лазерной биоспектроскопии, д. ф.-м. н., проф. В.Б. Лощенов

Москва-2013

Оглавление

ВВЕДЕНИЕ............................................................................................................".................................4

Актуальность работы.........................................................................................................................4

Цель и задачи......................................................................................................................................5

Практическая значимость..................................................................................................................6

Основные положения, выносимые на защиту.................................................................................6

Апробация работы..............................................................................................................................7

Содержание работы............................................................................................................................7

ГЛАВА 1. ОСОБЕННОСТИ СПЕКТРОСКОПИЧЕСКОГО АНАЛИЗА БИОЛОГИЧЕСКИХ СРЕД ШШО.........................................................................................................................................9

Введение..............................................................................................................................................9

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

Влияние биохимических особенностей биологических тканей на их оптические свойства.. 9

Влияние многократного рассеяния и геометрии измерений...................................................11

Влияние структурных особенностей биологических тканей на их оптические свойства ..14

Рассмотрение существующих спектроскопических методов анализа биологических сильнорассеивающих сред..............................................................................................................16

Сопоставление разработанных методов с существующими подходами к решению задачи спектрального анализа биологических сильнорассеивающих сред...........................................20

Заключение к первой главе.............................................................................................................23

ГЛАВА 2. МАТЕРИАЛЫ И МЕТОДЫ ИССЛЕДОВАНИЯ...........................................................24

Экспериментальная установка........................................................................................................24

Метод одновременной регистрации спектров диффузного отражения и флуоресценции с использованием оптоволоконного зонда....................................................................................24

Метод регистрации спектралъно-разрешенных изображений..............................................26

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

Взаимодействие электромагнитного излучения с независимыми частш(ами сферической и цилиндрической формы................................................................................................................28

Распространение электромагнитного излучения в сильнорассеивающих средах................32

Определение входных параметров математического моделирования взаимодействия света с нервными тканями.........................................................................................................43

Методы декомпозиции спектров....................................................................................................49

Объекты исследования....................................................................................................................50

Оптические фантомы биологических тканей, имитирующие их рассеивающие, поглощающие свойства и состав флуорофоров.......................................................................50

Клинические исследования...........................................................................................................51

Заключение ко второй главе............................................................................................................52

ГЛАВА 3. СПЕКТРОСКОПИЧЕСКИЙ АНАЛИЗ НЕРВНЫХ ТКАНЕЙ С ИСПОЛЬЗОВАНИЕМ ОПТОВОЛОКОННОГО ЗОНДА........................................................................................................53

Введение............................................................................................................................................53

Рассмотрение влияния геометрии измерений и оптических свойств среды на регистрацию диффузно-отраженного и флуоресцентного излучения...............................................................54

Разработка алгоритма разделения поглощающих и рассеивающих свойств.......................54

Апробация алгоритма определения концентрации гемоглобина по восстановленному спектру поглощения.....................................................................................................................67

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

Сопоставление результатов математического моделирования и измерения рассеивающих свойств глиом in vivo...................................................................................................................78

Учет рассеивающих свойств при определении концентрации протопорфирина IX..........82

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

Заключение к третьей главе............................................................................................................87

ГЛАВА 4. АНАЛИЗ СПЕКТРАЛЬНО-РАЗРЕШЕННЫХ ИЗОБРАЖЕНИЙ БИОЛОГИЧЕСКИХ ТКАНЕЙ................................................................................................................................................88

Введение............................................................................................................................................88

Пространственное разрешение системы визуализации................................................................89

Спектральное разрешение системы визуализации.......................................................................94

Разработка алгоритма разделения поглощающих и рассеивающих свойств при анализе спектрально-разрешенных изображений.......................................................................................96

Апробация метода анализа спектрально-разрешенных изображений в диффузно-отраженном и флуоресцентном свете на оптических фантомах.....................................................................100

Вычисление концентрации гемоглобина на основе анализа спектрально-разрешенных изображений в диффузно-отраженном свете.......................................................................100

Регистрация и анализ видеофлуоресцентного сигнала..........................................................101

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

Заключение к четвертой главе......................................................................................................106

ЗАКЛЮЧЕНИЕ..................................................................................................................................107

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

ЛИТЕРАТУРА....................................................................................................................................111

ВВЕДЕНИЕ

Актуальность работы

На настоящий момент при проведении флуоресцентного спектрального анализа in vivo в некоторых случаях его сочетают с измерением спектров диффузного отражения (ДО) [1, 2]. Коррекция спектров флуоресценции с помощью результатов анализа спектров диффузного отражения обеспечивает определение состава флуорофоров в ткани. В ряде случаев спектр диффузного отражения используется для независимого анализа поглощающих и рассеивающих свойств исследуемых сред [3]. Также следует отметить, что последовательное измерение спектров диффузного отражения и флуоресценции чревато возникновением артефактов, обусловленных взаимным смещением источника света, объекта и фотоприемника, и увеличением времени анализа, что может оказаться неприемлемым в клинических условиях. Что касается систем спектрально-разрешенной визуализации, то они, как правило, используют наборы узкополосных фильтров или перестраиваемые фильтры для получения изображений объекта, либо основаны на методе сканирования всего поля изображения поточечно с использованием классической спектроскопической системы с диспергирующим элементом [4]. У обоих подходов основным недостатком также является низкая скорость анализа данных.

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

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

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

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

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

Цель и задачи

Для чего были решены следующие задачи:

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

спектроскопический анализ в живом организме в реальном времени.

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

Практическая значимость

Система анализа мультиспектральных изображений, сопряжённая со стандартной щелевой лампой, была использована в ФГБУ «НИИ глазных болезней» РАМН для косвенной оценки уровня внутриглазного давления по степени оксигенации сосудов глазного дна и при проведении лабораторных исследований по эффективности фотодинамической терапии и приживлению тканей в ИОФ РАН.

Разработанные методы также нашли применение в научных исследованиях, выполняемых совместно с ФГБУ «ГНЦССП им. В.П. Сербского» Минздрава России.

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

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

2. Для восстановления информации о содержании в тканях гемоглобина в оксигенированной и дезоксигенированной форме по спектрально-разрешенным изображениям

может быть использовано линейное приближение зависимости коэффициента экстинкции от концентрации поглотителей.

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

Апробация работы

Содержание работы

Работа состоит из введения, четырех глав, заключения и списка литературы. Общий объем диссертации - 118 страниц, включая 54 рисунка и список литературы из 151 наименования.

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