СПОСОБ ОПРЕДЕЛЕНИЯ ОПТИЧЕСКИХ И БИОФИЗИЧЕСКИХ ПАРАМЕТРОВ БИОТКАНИ
METHOD FOR DETERMINING OPTICAL AND BIOPHYSICAL TISSUE PARAMETERS

Изобретение относится к области медицинского приборостроения.
Известны диагностическое устройство [1] и диагностический комплекс [2] для измерения медико-биологических характеристик кожи и слизистых оболочек in vivo включающий посылку оптического излучения на тестируемую ткань и регистрацию отраженного тканью излучения, которые определяют соответственно степень насыщения кислородом гемоглобина артериальной крови, частоту пульса и состояние кровотока в системе микроциркуляции. При этом коэффициент рассеяния ткани не определяется, а значит результат измерения вышеуказанных параметров подвержен существенному влиянию рассеивающих свойств ткани.
Известны также способы определения оптических и биофизических параметров биологических тканей по отраженному излучению с использованием интегрирующей сферы [3, 4] и оптоволоконного устройства, содержащего два соосных волновода, расположенных на фиксированном расстоянии друг от друга, для посылки и регистрации излучения [5]. В обоих случаях обратная задача по интерпретации спектральных коэффициентов диффузного отражения (КДО) ткани осуществляется путем минимизации разницы между экспериментальными и рассчитанными в рамках модели биоткани спектрами КДО. В [3, 4] предполагается измерение таких параметров кожи, как толщина эпидермиса, концентрация меланина и степень оксигенации гемоглобина. В [5] определяются коэффициент поглощения и транспортный коэффициент рассеяния слизистых оболочек полости рта, пищевода, желудка и легких, концентрации в них кровеносных сосудов и насыщение крови кислородом. К недостаткам данных способов следует отнести узкую область корректности используемых методов расчета характеристик переноса излучения, неоднозначность решения обратной задачи, вследствие невозможности разделения вкладов рассеяния и поглощения ткани в измеряемый спектр, и необходимость использования большого объема априорной информации об исследуемой среде (в частности, для корректной оценки поглощающих свойств ткани необходимо располагать информацией об ее рассеивающих свойствах, и наоборот).
s
s
s
s
s
Известны способы и устройства для определения оптических параметров однородной [6, 7] биоткани с использованием измерений рассеянного тканью света на нескольких расстояниях L от точки освещения. Способы позволяют определять коэффициент поглощения k и транспортный коэффициент µ рассеяния ткани. В [6] это достигается путем сопоставления расчетных и экспериментальных пространственных профилей диффузного отражения. Для расчета пространственного профиля КДО - R(L) используется уравнение полиномиальной множественной регрессии между R(L) и оптическими параметрами среды - k и µ, установленное на основе измерений R(L) для множества калибровочных образцов с известными оптическими параметрами. Данный способ требует строгое совпадения показателей преломления и индикатрис рассеяния реальных тканей и калибровочных образцов, что редко достигается в действительности. Кроме того, способ может давать неадекватные или не соответствующие реальности результаты, если значения параметров k и µ выйдут за пределы области, охватываемой калибровочными образцами. В способе [7] для восстановления параметров k и µ из измерений R(L) используется специально обученная нейронная сеть. Обучающий массив данных получен на основе расчета R(L) методом Монте-Карло для множества комбинаций значений k и µ. При этом в расчетах не учитываются вариации показателя преломления среды и фактора анизотропии рассеяния, что приводит к неудовлетворительным результатам использования данного метода для расстояний L<2 мм. В способах [6, 7] не учитывается многослойное строение ткани, поэтому их использование для кожи и других тканей, состоящих из нескольких слоев с принципиально отличающимся биохимическим составом и оптическими свойствами, не является корректным. Кроме того, данные способы предполагают, что КДО измеряются в спектральной области малых поглощений хромофоров ткани (меланин, билирубин, гемоглобин, оксигемоглобин и т.п.) - λ=660-1000 нм. Это обуславливает высокую чувствительность коэффициентов поглощения, определяемых на основе данных способов, к погрешностям измерения диффузного рассеяния и делает необходимым использования больших расстояний между точками освещения и регистрации, что, в свою очередь, приводит к невысокому уровню сигнала на фоне шумов измерительного устройства и большой вероятности искажения результатов измерений неоднородностью объема ткани, с которым взаимодействует излучение.
Наиболее близким к заявляемому изобретению является способ определения концентрации гемоглобина, основанный на измерении диффузного отражения ткани на двух изобестических длинах волн окси- и деоксигемоглобина, включающий посылку излучения на образец ткани, измерение диффузного отражения от ткани, вычисление отношения интенсивностей диффузного отражения на двух изобестических длинах волн с использованием измеренного диффузного отражения и определение концентрации гемоглобина в ткани посредством подстановки вычисленного отношения в предопределенное аналитическое выражение, соответствующее используемым изобестическим длинам волн [8]. К недостаткам данного способа следует отнести влияние вариаций рассеивающих свойств ткани и ее многослойного строения на точность измерений гемоглобина, а также невозможность определения комплекса других важных биофизических и оптических параметров ткани. Кроме того, при использовании данного способа для кожи, результаты измерения гемоглобина будут существенно искажены влиянием пигментации кожи.
Предлагаемое изобретение направлено на решение задачи расширения функциональных возможностей за счет одновременного определения комплекса как оптических, так и биофизических параметров биологической ткани в режиме реального времени, повышение точности измерения данных параметров за счет учета их общей вариативности, а также за счет исключения калибровочных измерений для нормированного r(Z, λ) спектрально-пространственного профиля коэффициента диффузного отражения ткани и использования априорной информации.
0
0
0
0
Для решения данной задачи посылку излучения на ткань в одну или несколько точек осуществляют на длинах волн λ из диапазона 350-1600 нм, измеряют диффузное отражение P(L, λ) на длинах волн посылаемого излучения для каждой из точек освещения, где L - расстояние между точками освещения и регистрации диффузного отражения на поверхности ткани, определяют абсолютный R(L, λ) или нормированный r(L, λ) спектрально-пространственный профиль коэффициента диффузного отражения ткани, как R(L, λ)=P(L, λ)/(P (λ)A(L,λ)) или r(L, λ)=P(L, λ)/P(L, λ), где Р(λ) - мощность излучения, посылаемого на ткань; A(L, λ) - аппаратурная константа, определяемая экспериментально с использованием эталонного отражателя, или рассчитываемая по известным геометрическим параметрам и спектральным характеристикам измерительного устройства; L - одно из расстояний между точками освещения и регистрации диффузного отражения; а оптические и биофизические параметры (X) определяют на основе аналитических выражений, представляющих собой множественные регрессии между X и R(L, λ) или между X и r(L, λ), которые получают путем измерения или расчета методом Монте-Карло коэффициентов диффузного отражения для множества образцов биоткани или моделирующих ее фантомов с известными оптическими и биофизическими параметрами, накопления ансамбля реализации оптических и биофизических параметров биоткани и соответствующих им спектрально-пространственных профилей R(L, λ), r(L, λ) для возможных диапазонов вариаций оптических и биофизических параметров ткани.
В качестве биоткани в предлагаемом изобретении используют кожу, молочную железу, слизистую оболочку полости рта, пищевода, органов желудочно-кишечного тракта и легких как in vivo, так и in vitro. Определяют такие оптические параметры ткани как коэффициент поглощения, транспортный коэффициент рассеяния, фактор анизотропии рассеяния, и такие биофизические параметры как объемная концентрация меланина, интегральное содержание меланина в эпидермисе, концентрация общего гемоглобина, насыщение гемоглобина кислородом, концентрация билирубина, средний диаметр и объемная концентрация капилляров с кровью, концентрация и размер эффективных рассеивателей ткани, содержание воды в ткани.
Коэффициенты диффузного отражения определяют или на характеристических длинах волн спектров поглощения и рассеяния компонентов исследуемой ткани, или же во всем диапазоне посылаемого излучения с последующей интерполяцией или экстраполяцией полученных коэффициентов R(L, λ) и r(L, λ) к длинам волн, соответствующим аналитическим выражениям, на основе которых осуществляется переход от R(L, λ) и r(L, λ) к определяемым параметрам ткани.
Предопределенные аналитические выражения между профилями R(L, λ), r(L, λ) и определяемыми оптическими и биофизическими параметрами ткани в виде множественных регрессий устанавливают путем измерения или расчета методом Монте-Карло коэффициентов диффузного отражения для множества образцов биоткани или моделирующих ее фантомов с известными оптическими и биофизическими параметрами, накопления ансамбля реализации оптических и биофизических параметров биоткани и соответствующих им спектрально-пространственных профилей R (L, λ), r(L, λ) и использования регрессионного или других математических методов анализа полученного ансамбля реализации.
Оптические и биофизические параметры ткани можно также определять на основе множественных регрессий между X и линейно-независимыми компонентами профиля R(L, λ) или r(L, λ), определяемыми как проекции профиля R(L, λ) или r(L, λ) на пространство из собственных векторов его ковариационной матрицы.
Свойства, появляющиеся у заявляемого технического решения, следующие:
1) повышение точности определения оптических и биофизических параметры ткани за счет учета их общей вариативности, исключения использования априорной информации, устранения калибровочных измерений для нормированного r(L, λ) спектрально-пространственного профиля коэффициента диффузного отражения ткани
2) расширение функциональных возможностей за счет одновременного определения комплекса оптических и биофизических параметров ткани в режиме реального времени, что не достигается с использованием известных способов.
а
Сущность данного изобретения поясняется с помощью фиг.1-8. На фиг.1,  изображена блок-схема устройства, реализующая предлагаемый способ, включающая блок светового излучения (1), оптоволоконный зонд (5), оснащенный освещающими (3) и принимающими (4) волокнами, монохроматор (2) или спектрометр (7), устройство регистрации (8) и устройство обработки измерительной информации (9). В качестве альтернативы оптоволоконному зонду 5 для посылки света на образец биоткани 6 и сбора отраженного ей потока может использоваться светопередающая оптика. Устройство регистрации 8 может быть выполнено как прибор с зарядовой связью (ПЗС). Вместо ПЗС могут использоваться фотоумножитель (ФЭУ) или фотодиоды. Блоком светового излучения 1 может служить лампа белого света (ксеноновая, аргоновая или криптоновая) или светодиоды (СД). Альтернативой монохроматору 2 является набор сменных фильтров. Оптоволоконный зонд 5 может включать эластичную трубку, содержащую множество оптических волокон (т.е., по крайней мере, одно освещающее и не менее двух принимающих), позволяя тем самым регистрировать спектрально-пространственный профиль сигналов диффузного отражения относительно освещающих волокон. Зонд 5 может быть также оснащен твердым наконечником, в котором множество волокон расставляются в различные геометрические конфигурации «освещение-регистрация».
Свет от источника излучения 1, пропущенный через монохроматор 2, поступает в оптоволоконный зонд 5, в котором по одному или нескольким освещающим волокнам 3 попадает на биоткань 6. Принимающие волокна 4 в оптоволоконном зонде 5 собирают диффузно отраженный свет и передают его на спектрометр 7. Сигналы диффузного отражения от различных принимающих волокон пространственно разделяются при помощи ПЗС 8, позволяя, таким образом, одновременно измерять спектры диффузного отражения для различных расстояний между каналами освещения и приема излучения.
На фиг.1, б представлен один из возможных вариантов расположения каналов возбуждения и регистрации сигналов диффузного отражения. Оптическое излучение подается в каналы возбуждения, расположенные на различных расстояниях L от точек регистрации сигналов отражения (в данном варианте по двум каналам). Сравнение профилей рассеянного излучения от каналов возбуждения позволяет оценить степень неоднородности освещаемого объема и тем самым выбрать оптимальный для проведения измерений участок объекта. Диффузное отражение измеряют или во всем спектральном диапазоне посылаемого на ткань излучения с некоторым спектральным разрешением, или же, на характеристических длинах волн спектров поглощения и рассеяния компонентов исследуемой биоткани с последующей интерполяцией или экстраполяцией к длинам волн, соответствующим используемым аналитическим выражениям для перехода к определяемым параметрам биоткани. Для получения пространственных значений коэффициентов отражения диффузное отражение измеряют не менее чем на двух расстояниях от точек освещения. После измерения спектра отражения ткани, полученные данные поступают на блок обработки 9, позволяющий получать абсолютные R(L, λ) или нормированные r(L, λ) спектрально-пространственные профили КДО ткани, а также вычисляющий значения оптических и биофизических параметров ткани (X) на основе множественных регрессий между X и спектрально-пространственными профилями КДО ткани или их линейно-независимыми компонентами.
0
Регистрируемые сигналы P(L, λ) зависят от спектрально-пространственного профиля КДО ткани - R(L, λ), а также от аппаратурной константы A(L, λ) и спектральной мощности источника света P(λ):
std
где A(L, λ)=G(L)τ(λ)S(λ), L - расстояние между освещающим и принимающим волокнами, G(L) - числовая апертура оптических волокон, S(λ) - спектральная чувствительности приемника, τ(λ) - функция пропускания оптической системы. Профили R(L, λ) определяются путем сопоставления соответствующих им сигналов P(L, λ) с сигналами P(L, λ), соответствующими эталонному диффузному отражателю:
std
dark
где R - КДО эталонного диффузного отражателя; P - мощность сигнала в отсутствии освещения, соответствующая темновому току.
0
0
0
Для устранения необходимости проведения калибровочных измерений параметры ткани следует определять по нормированному профилю ее КДО R(L, λ)/R(L, λ), который можно получить в виде отношения сигналов диффузного отражения для пространственно разнесенных каналов регистрации, как r(L, λ)=P(L, λ)/P(L, λ), где L - расстояние между освещающим и одним из светопринимающих волокон. Профили r(L, λ) не зависят от спектральных характеристик компонент измерительного устройства, а определяются только разностью оптических путей световых потоков для двух пространственно разнесенных каналов регистрации.
Для получения множественных регрессий, на основе которых определяются параметры ткани X, используется ансамбль реализации X и профилей R(L, λ), получаемый путем моделирования процесса переноса излучения в исследуемой ткани при широкой вариации ее параметров и расчета R(L, λ), r(L, λ) методом Монте-Карло [9] с учетом конструктивных параметров измерительного устройства (геометрической конфигурации волокон и их числовых апертур). Таким образом, аналитические выражения, связывающие X с профилями КДО ткани, получают путем многократных решений прямой задачи методом Монте-Карло (расчета КДО ткани) при различных комбинациях параметров оптической модели ткани и последующего статистического анализа полученных данных.
Ниже, возможности и достоинства предлагаемого способа демонстрируются на примере кожной ткани и слизистых оболочек полых органов (полости рта, пищевода, желудочно-кишечного тракта и легких).
s
Оптическая модель слизистой оболочки. Известно, что основными параметрами, характеризующими распространение оптического излучения в рассеивающей среде, являются коэффициенты поглощения k, рассеяния β и индикатриса рассеяния или ее средний косинус g. Причем, в оптике биологических тканей, как очень мутных сред, для описания световых полей достаточно знать не сами величины β и g, а их комбинацию - транспортный коэффициент рассеяния µ=β(1-g).
Спектр коэффициента поглощения слизистой оболочки k(λ) моделируется в виде линейной комбинации спектров поглощения ее хромофоров. В видимой области основными поглощающими хромофорами биоткани являются оксигемоглобин и деоксигемоглобин. В связи с этим, выражение для коэффициента поглощения имеет вид:
t
0
0
t
tHb
Hb
Hb
-1
-1
где k(λ)=A·exp[-B[λ-λ)] - коэффициент поглощения обескровленной ткани, λ=632 нм, A (мм) и B (нм) - параметры, характеризующие спектральную зависимость k; F (г/л) - концентрация общего гемоглобина в биологической ткани; µ=64.5 мг/моль - молярная масса гемоглобина;  и ε - молярные коэффициенты поглощения соответственно оксигемоглобина и деоксигемоглобина [10]; S - насыщение гемоглобина кислородом; α - поправочный коэффициент, учитывающий эффект локализованного поглощения света кровеносными сосудами [11]. Суть введения коэффициента α состоит в учете того, что кровь не равномерно распределена в объеме ткани, а локализована в капиллярах. Для хаотически распределенных капилляров с диаметром D коэффициент α определяется на основе следующего выражения [11]:
Hb
где  - коэффициент поглощения крови, C=150 г/литр - концентрация гемоглобина в крови.
Транспортный коэффициент рассеяния биологических тканей в видимой области спектра с хорошей точностью аппроксимируется степенной зависимостью:
0
s
где λ=632 нм; C и ν - структурные параметры ткани, характеризующие объемное содержание и размер ее «эффективных» рассеивателей. Для описания процесса рассеяния излучения используется однопараметрическая функция Хеньи-Гринштейна [9] с фактором анизотропии рассеяния g.
tHb
t
Таким образом, оптическая модель слизистой оболочки определяется параметрами g, A, B, C, ν, F, D, S. Кроме того, для учета отражения падающего света от поверхности ткани, а также многократных переотражений вызванного им диффузного излучения между внутренними слоями и поверхностью ткани необходимо располагать информацией о ее показателе преломления n. В табл.1 приведены диапазоны вариаций модельных параметров, выбранные путем критического анализа результатов различных авторов [5, 12-22] для нормальных и опухолевых тканей организма (слизистых оболочек полости рта, пищевода, органов желудочно-кишечного тракта и легких).
-1
-1
-1
tis
Вариации модельных параметров будем осуществлять независимо друг от друга, но при этом будем контролировать, чтобы при каждой их комбинации значения оптических параметров ткани, получаемые по формулам (1) и (2), не выходили из диапазонов, наблюдаемых в эксперименте. Исходя их вышеуказанных литературных данных, выбраны следующие ограничения: 1) β'(700 нм)≥0.2 мм; 2) k(632 нм)=0.02-0.5 мм; 3) k(450 нм)≤0.35 мм; 4) Λ'(632 нм)=0.5-0.98, где Λ'=β'/(k+β') - транспортное альбедо однократного рассеяния ткани.
Табл.1
Диапазоны вариаций параметров оптической модели слизистой оболочки
Параметр
Диапазон
n<sub>t</sub>
1.35-1.45
A
0.01-0.1 мм<sup>-1</sup>
B
0.001-0.028 нм<sup>-1</sup>
C<sub>S</sub>
0.5-3.0 мм<sup>-1</sup>
x
1.0-2.5 мкм
g
0.5-0.95
F<sub>tHb</sub>
1.5-35 г/л
D
5-30 мкм
S
40-98%
epi
skin
Оптико-статистическая модель кожи. Верхний слой кожи - эпидермис с толщиной L, нижний - дерма, которая в оптическом плане считается бесконечно толстой. Показатель преломления слоев кожи n относительно воздуха считается одинаковым, поэтому френелевское отражение излучения имеет место только на границе раздела кожи с внешней средой.
t
-1
Полагаем, что основой кожи является слабопоглощающая, обескровленная ткань. Зависимость ее коэффициента поглощения k [см] от длины волны λ [нм] аппроксимируется следующим выражением:
В качестве индикатрисы рассеяния "эффективных" рассеивателей ткани (основными из которых являются упакованные в жгуты коллагеновые и эластиновые волокна дермы, а также кератиноциты и меланоциты эпидермиса) используем функцию Хеньи-Гринштейна, спектральная зависимость фактора анизотропии (среднего косинуса) которой может быть описана эмпирическим соотношением:
s
s
Спектр транспортного коэффициента рассеяния µ(λ) кожи в видимой и ближней ИК-областях спектра рассчитывается как суперпозиция спектров µ(λ), соответствующих рассеивателям Ми и Рэлея с размерами d≥λ и d<<λ соответственно:
0
Mie
где λ=400 нм; ρ - доля рассеяния Ми; x - параметр спектральной зависимости редуцированного коэффициента рассеяния Ми, зависящий от размера рассеивателей и их показателя преломления.
mel
epi
Коэффициенты поглощения слоев кожи определяются по известным правилам сложения оптических величин как сумма коэффициентов поглощения отдельных компонент с весами, равными их объемным концентрациям. Спектр коэффициента поглощения эпидермиса рассчитывается в соответствии с объемными концентрациями меланина f и воды w в его составе:
mel
w
где k и k - коэффициенты поглощения меланина и воды.
blood
derm
bil
Коэффициент поглощения дермы зависит от содержания в ней капилляров с кровью f, воды w и билирубина G:
bil
bil
blood
-1
где ε - молярный коэффициент поглощения билирубина в см/(моль/л); µ=585 г/моль - молярная масса билирубина; k - коэффициент поглощения крови; α - корректирующий фактор, учитывающий эффект локализованного поглощения света кровеносными сосудами.
blood
Коэффициент поглощения крови k в зависимости от длины волны излучения может быть представлен в виде суммы коэффициентов поглощения гемоглобина (окисленного и не окисленного) и билирубина:
Hb
Hb
Hb
dif
-1
где C [г/л] - концентрации общего гемоглобина в крови; µ=64500 г/моль - молярная масса гемоглобина;  и ε - молярные коэффициенты поглощения оки- и деоксигемоглобина в см моль/литр); S - насыщение крови кислородом (доля окисленного гемоглобина в общем гемоглобине); t=5 - отношение концентраций билирубина в крови и в окружающей ткани (полагается фиксированным).
skin
s
0
Mie
epi
mel
epi
blood
derm
bil
tHb
Таким образом, оптическая модель кожи определяется следующими параметрами: показатель преломления цельной кожи - n, параметры спектральной зависимости редуцированного коэффициента рассеяния - µ(λ), ρ и x; толщина эпидермиса L; объемные концентрации меланина и воды в эпидермисе -f и w; объемные концентрации капилляров с кровью, воды и билирубина в дерме - f, w и C; средний диаметр капилляров - D; концентрация общего гемоглобина в крови - C; насыщение крови кислородом - S. Диапазоны вариаций модельных параметров, выбранные путем критического анализа результатов различных авторов, приведены в табл.2.
Схема моделируемого эксперимента. Будем полагать, что излучение инжектируются в среду посредством передающего волновода с диаметром сердцевины 0.2 мм. Рассеянное тканью в заднее полупространство излучение попадает в принимающие волноводы, расположенные на расстояниях L=0.23, 0.46, 0.69, 0.92, 1.15 мм от центра передающего волновода (диаметр сердцевины всех волноводов 0.2 мм). Для кожи канал регистрации с L=1.15 мм не рассматривается из-за низкого уровня сигнала при λ≤500 нм на фоне шума за счет сильного поглощения меланином и билирубином.
Табл.2
Диапазоны вариаций параметров оптической модели кожи
Параметр
Диапазон
Параметр
Диапазон
n<sub>skin</sub>
1.4-1.5
f<sub>mel</sub>L<sub>epi</sub>, мкм
0.6-10
µ<sub>s</sub>(λ<sub>0</sub>), мм<sup>-1</sup>
3-11
f<sub>blood</sub>, %
0.2-7.0
x
0.5-1.0
C<sub>tHb</sub>, г/л
120-200
ρ<sub>Mie</sub>
0.1-0.6
S, %
40-98
µ<sub>s</sub>(λ<sub>0</sub>)/µ<sub>s</sub>(900)
3-10
C<sub>bil</sub> мг/л
0.0'1-50
L<sub>epi</sub>, мкм
50-150
D, мкм
5-30
s
s
6
5
3
Расчет R(L, λ) осуществляется по следующей схеме. Случайным образом выбираются значения модельных параметров из диапазонов, указанных в табл.1 и 2. Для каждой реализации параметров производится расчет k(λ) и µ(λ) по формулам (1)-(8) на 26 длинах волн (λ=450, 460, 470, 480, 495, 506, 514, 522, 530, 540, 548, 560, 564, 568, 574, 580, 586, 594, 600, 610, 620, 630, 640, 660, 676, 700 нм) для слизистой оболочки и на 30 длинах волн (λ=452, 460, 470, 480, 495, 506, 514, 522, 530, 540, 548, 560, 564, 568, 574, 580, 586, 594, 600, 610, 620, 630, 640, 660, 676, 700, 720, 740, 760, 780, 800 нм) для кожи. Указанные λ характеризуют основные особенности спектральных кривых поглощения хромофоров ткани в видимой области спектра. В соответствии с рассчитанными спектральными значениями оптических параметров, а также фактора анизотропии рассеяния и показателя преломления ткани методом Монте-Карло моделируется профиль R(L, λ). Фантом ткани представляет собой цилиндр с радиусом 10 мм и высотой 5 мм. Такие размеры выбираются с тем расчетом, чтобы, с одной стороны, устранить эффект границ, а с другой - уменьшить время счета. В расчете принимается, что излучение вводится в среду в направлении нормали к ее поверхности и равномерно распределено по площади сечения освещающего световода. Числовые апертуры всех волокон условно приняты равными 1.0 (разумеется, что при решении задачи интерпретации реальных оптоволоконных измерений расчет профилей R(L, λ) должен осуществляться в строгом соответствии со значениями апертур волокон и угловым распределением инжектируемого в среду излучения). Количество фотонов, вводимых в среду, зависит от транспортного альбедо однократного рассеяния среды - Λ=µ/k и варьируется от 6·10 при Λ≤0.5 до 2·10 при Λ≥10. Данная процедура повторяется для 2·10 реализации R(L, λ), что, как будет видно ниже, является достаточным для получения статистически-значимых результатов.
0
Следует отметить, что при использовании в качестве источников излучения светодиодов, расчет диффузного отражения в соответствующих им спектральных интервалах Δλ должен осуществляться с учетом реальной формы линии излучения P(λ) и спектральной чувствительности приемника S(λ). Для того, чтобы избежать трудоемких вычислений монохроматического КДО и его последующего интегрирования в интервале Δλ с учетом аппаратурных функций, можно моделировать длину волны каждого фотона как случайное число, равномерно распределенное на интервале интегрирования и по значениям счетчиков вылетающих из среды фотонов получать оценку значений искомого интеграла как:
λ
ph,i
i
i
где N - число подинтервалов на которые разбивается интервал интегрирования Δλ; N - общее число фотонов с длиной волны из i-го подинтервала, пропущенных через среду (начальный "вес" каждого фотона равен единице); r - суммарный "вес" фотонов с длиной волны из i-го подинтервала, вылетевших из среды; w - квадратурные коэффициенты, зависящие от используемого метода численного интегрирования.
λ
L
λ
L
λ
L
λ
L
Решение обратной задачи. На фиг.2 изображен пример смоделированных спектральных зависимостей КДО слизистой оболочки на пяти различных расстояниях от точки освещения, образующих в совокупности спектрально-пространственный профиль R(L, λ). На фиг.3 изображен аналогичный профиль для кожи. Обратная задача состоит в восстановлении модельных параметров из профилей R(L, λ), Совокупность значений R(L, λ), определяемых в моделируемом эксперименте, удобно представлять в виде вектора измерений r из NN компонент, где N и N - количество спектральных и пространственных каналов регистрации КДО ткани соответственно (для кожи N=30, N=4; для слизистой оболочки N=26, N=5). Для удобства работы с числами одного порядка будем использовать логарифмическое представление измерений, т.е. компонентами r будут lnR(L, λ).
С помощью метода главных компонент [23, 24] вектор r разлагается по системе ортогональных векторов, причем оптимальной аппроксимацией r является разложение по собственным векторам его ковариационной матрицы
k
mes
где  и σ - среднее значение и дисперсия r, рассчитанные по набору векторов r, K - число векторов в наборе (количество реализации параметров вышерассмотренной модели кожи); 1≤i, j≤N. Искомое разложение любой реализации r имеет вид:
mes
PC
n
PC
n
где V - матрица размером N×N, столбцами которой являются собственные векторы v матрицы S; N - число главных компонент; ξ - вектор из коэффициентов разложения ξ (линейно-независимых компонент), которые находятся по формуле
В связи с быстрой сходимостью рассматриваемого разложения на первые ортогональные функции (собственные векторы) приходится большая часть изменчивости компонент вектора r. Таким образом, метод главных компонент позволяет значительно уменьшить размерность исходных данных и выделить из них несколько линейно-независимых компонент, содержащих в себе столько же информации, сколько ее было в исходных данных. Очевидно, что при малом количестве измеряемой информации, когда суммарное количество спектральных и пространственных каналов регистрации не превышает количества линейно-независимых компонент в профилях R(L, λ), определять искомые параметры ткани можно напрямую из измеряемых значений R(L, λ) без извлечения из них линейно-независимых компонент, например при помощи полиномиальных регрессий
ij
i
j
ij
i
j
ij
i
j
ij
i
j
λ
L
ijm
a
где R=R(L, λ) или R=lnR(L, λ); r=r(L, λ) или r=lnr(L, λ); N и N - соответственно количество длин волн и расстояний от точки освещения, на которых осуществляются измерения диффузного отражения,  - коэффициенты регрессий, численные значения которых определяются методом наименьших квадратов; M - степень полинома; X - оптический или биофизический параметр (или его логарифм).
mes
При большом количестве оптических измерений N восстанавливать параметры ткани предпочтительно не из самих КДО, измеряемых непосредственно в эксперименте, а из коэффициентов их разложения по ортогональным функциям, например, по собственным векторам ковариационной матрицы (9). При этом для определения количественных значений оптических биофизических параметров биоткани удобно использовать полиномиальные регрессии
a
nm
n
PC
n
PC
L
λ
где  - коэффициенты регрессий; M≥1 - степень полинома; X - оптический или биофизический параметр (или его логарифм); ξ - линейно-независимые компоненты спектрально-пространственных профилей R(L, λ), r(L, λ) или их логарифмов, N - количество главных компонент (собственных векторов v), используемых для определения параметра X (N=1, …, NN). Использование в (14) логарифмического представления измерений и искомых параметров позволяет оперировать с числами одного порядка и устраняет возможность получения отрицательных значений определяемых параметров.
n
PC
PC
n
PC
*
Степень полиномов M в (14) выбирается, исходя из среднеквадратичной погрешности аппроксимации статистической связи между ξ и X. Как правило, для аппроксимации достаточно использовать М=3. Оптимальное количество главных компонент N в (14) для каждого из модельных параметров X определяется путем замкнутого численного эксперимента, заключающегося в следующем. Первоначально по формуле (11) определяются линейно-независимые компоненты всех реализации профилей R(L, λ) и устанавливаются регрессии между первыми N компонентами ξ и параметром X. Далее перебираются все реализации модельных параметров и для каждой реализации осуществляется расчет X с использованием (11), (14) при наложении на R(L, λ) случайных отклонений в пределах δR (моделирующих погрешности измерения КДО). Полученное в результате значение X сравнивается со значением X, соответствующем рассматриваемой реализации, и рассчитывается погрешность восстановления X. После перебора всех реализации вычисляется средняя погрешность восстановления X. Оптимальное значение N соответствует минимуму этой погрешности.
tHb
s
PC
tHb
s
PC
*
На фиг.4 приведены результаты вышеописанного численного эксперимента для следующих параметров слизистой оболочки: F - концентрация общего гемоглобина в ткани; S - насыщение гемоглобина кислородом; C и ν - параметры спектральной зависимости редуцированного коэффициента рассеяния ткани, характеризующие объемную концентрацию и размер «эффективных» рассеивателей ткани соответственно. Зависимости средних по ансамблю погрешностей восстановления модельных параметров от количества главных компонент получены при наложении на R(L, λ) случайных отклонений в пределах 1% (окружности), 5% (квадраты) и 10% (треугольники). Видно, что при малых погрешностях оптических измерений (δR=1%) погрешности восстановления модельных параметров - δX уменьшаются с увеличением N. При δR=5-10% для восстановления рассматриваемых параметров оптимально использовать от 9 до 12 главных компонент. На фиг.5 приведены результаты сопоставления значений параметров F S, C и ν, восстановленных из смоделированных профилей R(L, λ) при N=12, со значениями этих же параметров, соответствующих модельным реализациям R(L, λ). На рисунках также указаны средние погрешности восстановления параметров. Поскольку восстановление искомых параметров осуществлялось без наложения на профили R(L, λ) случайных отклонений, то разброс точек на диаграммах относительно прямой X=X характеризует чувствительность результатов восстановления параметра X к вариациям всех других модельных параметров.
tHb
blood
tHb
bil
epi
mel
epi
s
s
s
epi
tHb
bil
s
0
epi
0
blood
blood
bil
blood
blood
ν
На фиг.6 представлены результаты восстановления биофизических параметров кожи. Вышеописанный численный эксперимент проводился для наиболее важных диагностических параметров: концентрация общего гемоглобина в дерме F=fC (г/л); насыщение гемоглобина кислородом S; концентрация билирубина в крови C (мг/л); интегральное содержание меланина в эпидермисе OD=fL (мкм); параметры C и ν спектральной зависимости редуцированного коэффициента рассеяния кожи µ(λ)=C·(632/λ) в диапазоне 600-700 нм, характеризующие структурные свойства ткани (объемную концентрацию и размер «эффективных» рассеивателей). Параметр OD восстанавливался из абсолютных значений R(L, λ). Для восстановления параметров F, S, C, C и ν использовались нормированные профили КДО - r(L, λ), где L=0.23 мм. Невозможность восстановления параметра OD из профилей r(L, λ) связана с тем, что профили r(L, λ) практически не зависят от меланиновой пигментации кожи, поскольку оптические пути, проходимых светом в тонком верхнем слое кожи (эпидермисе), для рядом расположенных принимающих волокон примерно одинаковы и поэтому слабо влияют на отношение соответствующих им сигналов P(L, λ) и P(L, λ). Следует отметить, что использование для решения обратной задачи профилей r(L, λ) позволяет устранить влияние аппаратурных констант и мощности посылаемого излучения на точность восстановления искомых параметров, а значит устройство, реализующее измерения r(L, λ) не нуждается в калибровке. Погрешность восстановления параметра S зависит от концентрации капилляров с кровью f и резко возрастает с ее уменьшением. В связи с этим коэффициенты регрессии (3) для данного параметра получены на основе реализации r(L, λ), соответствующих значениям f более 1%. При получении регрессий для параметра C, напротив, использовались реализации r(L, λ) для которых f<3%, поскольку более высокие значения f не характерны для нормальной кожи.
tHb
bil
epi
epi
tHb
s
tHb
bil
Знание параметров F, S, C, OD, C и ν имеет большое значение для решения широкого круга задач диагностики и терапии различных заболеваний. Так, имеется ряд данных о различии параметров OD, F, C и ν в нормальных и опухолевых тканях, поэтому информация о них может использоваться для поиска и диагностики злокачественных новообразований (например, меланомы). Параметр F также может использоваться для обнаружения внутреннего кровоизлияния и для непрерывного мониторинга уровня гемоглобина в крови при хирургических операциях и переливаниях крови. Параметр C имеет огромное значение при мониторинге желтухи у новорожденных.
Определение концентрации гемоглобина. Для сравнения с прототипом рассмотрим также возможности предлагаемого способа по отношению к гемоглобину. Вышерассмотренная обратная задача решалась для большого количества спектральных и пространственных каналов регистрации КДО. Это с одной стороны позволяет повысить устойчивость решения обратной задачи, с другой - увеличивает время измерения, требует использования дорогостоящего спектрометрического оборудования и сложной оптики блока регистрации рассеянного излучения. Кроме того, для проведения абсолютных измерений R(L, λ) необходимо знание аппаратурных констант и мощности посылаемого на ткань излучения, т.е. устройство, реализующее данные измерения, требует калибровки. С точки зрения экономичности и простоты реализации наибольший интерес представляет метод определения параметров ткани, основанный на измерении ее КДО в небольшом количестве спектральных участков, соответствующих излучению нескольких светодиодов. Устранить необходимость калибровки таких измерений можно, если определять параметры ткани из разности логарифмов сигналов диффузного отражения P(L, λ) для двух пространственно разнесенных каналов регистрации, которая, при равенстве апертур принимающих волокон, определяется только разностью оптических путей соответствующих световых потоков.
tHb
2
1
В связи со сказанным выше рассмотрим возможность определения концентрации общего гемоглобина F по измерениям отношения сигналов r(λ)=P(L, λ)/P(L, λ) в небольшом количестве спектральных участков. Следует отметить, что именно локальное увеличение содержания гемоглобина в месте образования опухоли является наиболее существенным признаком патологических процессов [5, 21], поэтому разработка простого и эксплуатационно-надежного способа неинвазивного определения концентрации гемоглобина в биологических тканях является чрезвычайно актуальной задачей.
tHb
λ
tHb
λ
1
2
tHb
tHb
tHb
Минимальное количество спектральных участков, необходимое для восстановления F, можно определить на основе замкнутых численных экспериментов. Компонентами вектора измерений r являются спектральные значения r(λ), поэтому количество линейно-независимых компонент r соответствует количеству независимых длин волн оптического зондирования N. Анализ зависимости погрешности восстановления F от N проводился для отношения сигналов, соответствующего наибольшему пути между принимающими световодами в вышерассмотренной схеме оптических измерений, т.е. L=0.23 мм и L=1.15 мм. Результаты такого анализа показывают, что при погрешности измерения r(λ), равной 5-10%, для определения F достаточно четырех длин волн оптического зондирования. Таким образом, задача планирования измерений F сводится к выбору четырех оптимальных спектральных участков. Такой выбор осуществлялся путем компьютерного перебора всех возможных комбинаций из исходного набора λ (26 значений) и оценивания соответствующих им погрешностей восстановления F.
tHb
Очевидно, что интерпретация измерений r(λ) в четырех спектральных участках невозможна с использованием строгих математических методов решения обратных задач, поскольку количество измерений меньше количества модельных параметров. Однако предлагаемый способ позволяет решать обратную задачу и при малом количестве оптических измерений. В данном случае в качестве оператора решения обратной задачи (оператора преобразования r(λ) в F) использовались полиномиальные регрессии следующего вида:
n
tHb
n
tHb
tHb
n
n
tHb
n
tHb
n
1
2
3
4
tHb
Оптимальный набор λ в уравнении (15) определяется точностью восстановления F и устойчивостью данного уравнения к случайным отклонениям δr соответствующих спектральных значений r(λ). Влияние δr на результат восстановления F оценивается на модельном ансамбле реализации F и r(λ) путем проведения замкнутого численного эксперимента в котором для всех реализации r(λ) осуществляется расчет F с использованием (15) при наложении на r(λ) случайных отклонений в пределах δr=5%. После перебора всех реализации вычисляется средняя по ансамблю погрешность восстановления F, которая и служит критерием выбора оптимальной комбинации λ. Для рассматриваемого набора из 26 длин волн такой комбинацией является λ=480 нм; λ=574 нм; λ=586 нм; λ=630 нм. Соответствующая ей погрешность восстановления F составляет 9.2 и 10.8% при δr=1 и 5% соответственно. Коэффициенты регрессии (15), соответствующие выбранным длинам волн, приведены в табл.3.
Табл.3
Коэффициенты <i>a</i><sub>nm</sub> регрессии (15) при <i>a</i><sub>00</sub>=-4.1742
n
m=1
m=2
m=3
1
2.1791
0.3204
0.0183
2
-7.5811
-1.2871
-0.0677
3
4.6361
1.1887
0.0706
4
-3.1531
-1.5916
-0.1688
s
s
s
s
-1
-1
-1
Определение оптических параметров. Известно, что оптические параметры (ОП) биологических тканей (коэффициент поглощения k, транспортный коэффициент рассеяния µ и фактор анизотропии рассеяния g) несут обширную информацию об их структуре и биохимическом составе. Помимо этого, знание ОП тканей необходимо для оценки глубины проникновения в них излучения, определяющей выбор дозы облучения при фототерапии различных заболеваний. В связи с этим рассмотрим более детально способ определения ОП биологических тканей по измерениям их КДО на нескольких расстояниях L от точки освещения. Исходя из точности и оперативности определения ОП, несомненный практический интерес представляет возможность получения значений ОП на основе аналитического преобразования измеряемых пространственных профилей КДО - R(L). Оператор такого преобразования (оператор решения обратной задачи) может быть получен на основе смоделированного ансамбля реализаций k(λ), µ(λ), g(λ) и R(L, λ). Учитывая, что процедура восстановления ОП из измерений R(L) может выполняться безотносительно к длине волны зондирующего излучения, сформируем на основе смоделированных спектральных значений k(λ), µ(λ), g(λ) и R(L, λ) один общий ансамбль, соответствующий диапазонам вариаций ОП, характерным для большинства биологических тканей в видимой области спектра: 0.01≤k≤5 мм, 0.45≤µ≤6 мм, 0.5≤g≤0.95 мм. При этом из ансамбля исключены те реализации, для которых отношение КДО для наиболее близкого и наиболее удаленного каналов регистрации превышает 4 порядка, поскольку динамический диапазон соответствующих им сигналов выходит за пределы линейности применяемых приемников излучения.
i
s
Количество линейно-независимых компонент в измерениях R(L) (Z=0.23, 0.46, 0.69, 0.92, 1.15 мм) можно определить путем анализа собственных чисел матрицы ковариаций R(L). Каждое из собственных чисел l (i=1, …, 5) этой матрицы определяет относительный вклад соответствующего ему собственного вектора в вариации R(L). Доля, вносимая первыми n собственными векторами в суммарную дисперсию R(L) определяется как  [23, 24]. В нашем случае на три первых собственных вектора, соответствующих наибольшим собственным числам матрицы ковариаций R(L), приходится 99.99% суммарной дисперсии R(L). Этот факт является вполне понятным, если принять во внимание, что вариации профилей R(L) обусловлены влиянием трех основных ОП - k, µ и g. Таким образом, для определения ОП ткани в принципе можно ограничиться измерениями ее КДО на трех расстояниях от точки освещения. Однако мы все же будем использовать для восстановления ОП профили R(L), соответствующие пяти вышеуказанным значениям L, что позволяет повысить устойчивость решения обратной задачи и, тем самым, снизить требования к точности измерений R(L). При этом для решения обратной задачи выбраны операторы вида
s
1
2
3
n
n
mnk
n
a
где P - один из ОП (k, µ или g); ξ, ξ, ξ - линейно-независимые компоненты профилей R(L), для расчета которых используется формула (11), в которой r - вектор измерений с компонентами, равными R(L) (n=1, …, 5), a v - собственные векторы его ковариационной матрицы;  - коэффициенты, численные значения которых рассчитываются методом наименьших квадратов на основе ансамбля реализации ξ и P.
s
Результаты замкнутых численных экспериментов по восстановлению ОП с использованием формул (11), (16) позволяют сделать следующие выводы. Для восстановления параметра g погрешность измерения R(L) должна быть ~ 1-3%. Так, при δR=1% информативность профилей R(L) относительно данного параметра (как отношение априорной неопределенности g к апостериорной) составляет примерно 5.6, но уже при δR=7% апостериорная неопределенность параметра g практически равна априорной. Эти результаты не вызывают особого удивления, поскольку хорошо известно, что световые поля в сильно рассеивающих средах (к которым относятся биологические ткани) определяются не величинами их коэффициента рассеяния β и фактора анизотропии рассеяния g по отдельности, а произведением µ=β(1-g). Для более точного восстановления параметра g необходимы измерения R(L) для расстояний L от точки освещения, на которых преобладает однократное или малократное рассеяние. Следует отметить, что слабое проявление индикатрисы рассеяния в рассматриваемых профилях R(L), оправдывает правомерность используемого при моделировании предположения об отсутствии спектральной зависимости у параметра g.
s
s
s
s
n
mnk
a
Средняя по ансамблю реализации модельных параметров погрешность восстановления µ составляет 2.9 и 4.6% при δR=1 и 5% соответственно. Коэффициент поглощения k при тех же значениях δR восстанавливается с погрешностями 8.2 и 9.3%, причем наибольшие погрешности соответствуют спектральной области малых поглощений биологических тканей λ>600 нм. Последнее обстоятельство объясняется малыми оптическими путями между каналами возбуждения и регистрации (L≤1.15 мм), на которых изменение интенсивности зондирующего излучения за счет поглощения хромофорами ткани находится в пределах погрешности измерений и, кроме того, маскируется значительно более сильным влиянием рассеяния. Повысить точность восстановления параметра k можно путем использования дополнительных каналов регистрации, расположенных на расстояниях L≥1.15 мм от канала возбуждения. При этом максимальное значение L не должно превышать диаметр инструментального канала эндоскопа, типичные значения которого составляют 2-3 мм. В связи с этим рассмотрим, как изменится точность восстановления параметров k, µ, и g с использованием дополнительных каналов регистрации с L=1.38, 1.61, 1.84, 2.07 мм. Для этого разобьем ансамбль реализации параметров k, µ, g и R(L) на две части, соответствующие спектральным областям сильного (λ<600 нм) и слабого (λ≥600 нм) поглощения. Полагаем, что для восстановления k и µ в области слабого поглощения доступны все каналы регистрации, а для области сильного поглощения только каналы с L≤1.15 мм (поскольку регистрировать оптические сигналы с более удаленных участков ткани сложно из-за их низкого уровня на фоне «шума»). При этом для восстановления ОП в обоих случаях будем использовать формулы (11), (16). Собственные векторы v ковариационной матрицы вектора измерений и коэффициенты  для спектральных областей сильного и слабого поглощения ткани, входящие в формулы (11), (16) приведены в табл.4 и 5.
Табл.4
Собственные векторы ковариационной матрицы lnR(L)
L, мм
λ<600 нм
λ≥600 нм
v<sub>1</sub>
v<sub>2</sub>
v<sub>3</sub>
v<sub>1</sub>
v<sub>2</sub>
v<sub>3</sub>
0.23
0.4109
0.4605
0.4686
0.2338
0.2962
0.3423
0.46
-0.7142
-0.2815
0.0681
-0.5382
-0.4343
-0.3031
0.69
-0.5359
0.5001
0.4477
0.6809
0.0325
-0.2808
0.92
0.1809
-0.6251
0.4055
0.4009
-0.5226
-0.3453
1.15
0.0353
-0.2605
0.6411
-0.1652
0.6625
-0.6520
1.38
-
-
-
-0,0504
0.0682
0.4002
1.61
-
-
-
-0.0209
0.0348
0.0125
1.84
-
-
-
0.0298
-0.0683
-0.0022
2.07
-
-
-
0.0086
0.0031
-0.0862
Табл.5
Коэффициенты множественной регрессии (16)
<i>a</i><sub>mnk</sub>
λ<600 нм
λ≥600 нм
k
µ<sub>s</sub>
g
k
µ<sub>s</sub>
g
<i>a</i><sub>000</sub>
1.2200
1.2126
0.7254
0.0971
0.8206
0.7151
<i>a</i><sub>100</sub>
-0.1444
0.2759
-0.0542
-0.0347
0.1297
0.0852
<i>a</i><sub>200</sub>
-0.0143
0.0228
0.0085
-0.0029
-0.0235
-0.0056
<i>a</i><sub>300</sub>
0.0010
-0.0010
0.0008
0.0012
-0.0078
0.0036
<i>a</i><sub>010</sub>
-0.5224
-0.9009
0.5625
-0.0735
-0.4859
0.2186
<i>a</i><sub>020</sub>
-0.1589
0.1871
-0.0693
-0.0023
0.1128
-0.0099
<i>a</i><sub>030</sub>
0.0985
0.1511
-0.0555
0.0062
-0.0284
-0.0125
<i>a</i><sub>001</sub>
-0.3323
-0.0983
1.9220
-0.0277
0.2097
-1.3878
<i>a</i><sub>002</sub>
-3.7763
-4.9956
-3.4440
-0.3861
-1.5089
-2.1737
<i>a</i><sub>003</sub>
6.0127
9.7284
-0.4417
-0.3892
2.3523
-2.0063
<i>a</i><sub>110</sub>
0.1089
-0.2192
-0.0359
0.0265
-0.0527
-0.0337
<i>a</i><sub>210</sub>
-0.0009
-0.0089
-0.0036
-0.0011
0.0175
-0.0134
<i>a</i><sub>120</sub>
-0.0123
0.0026
0.0082
0.0007
-0.0450
0.0299
<i>a</i><sub>101</sub>
-0.0780
-0.3699
-0.0075
0.0780
0.4820
0.3514
a<sub>201</sub>
0.0044
-0.0870
0.0043
-0.0232
0.0868
-0.0530
<i>a</i><sub>102</sub>
0.1859
-0.6871
0.0531
0.1773
-0.4478
0.5923
a<sub>011</sub>
-0.9082
-1.0915
-1.1694
-0.0663
-0.2365
0.2550
<i>a</i><sub>021</sub>
0.8677
1.7832
-0.4849
-0.0174
0.5851
-0.0784
<i>a</i><sub>012</sub>
3.4995
7.7425
-1.3934
0.0123
-1.7248
0.1619
<i>a</i><sub>111</sub>
0.0000
0.0000
0.0000
0.0000
0.0000
0.0000
s
n
n
s
О точности восстановления параметров k, µ и g в условиях их общей вариативности можно судить по представленным на фиг.7 диаграммам, которые иллюстрируют результаты восстановления ОП слизистой оболочки для спектральных областей сильного (а, в, д) и слабого (б, г, е) поглощений. Данные диаграммы получены в предположении, что профили R(L) измеряются с погрешностью, равной 1%. В табл.6 даны средние погрешности восстановления ОП при наложении на значения R(L) (=1, …, 9) случайных отклонений δR в пределах 5 и 10%. Видно, что погрешности восстановления µ практически не отличаются от аналогичных погрешностей, соответствующих всей видимой области спектра. Погрешность же восстановления параметра k уменьшается более чем два раза. Результаты восстановления g по-прежнему являются сильно чувствительными к погрешностям измерений R(L), однако при δR=1-3% данный параметр можно восстанавливать с удовлетворительной точностью.
Табл.6
Погрешности восстановления ОП ткани с использованием регрессий (16)
λ
<600 нм
≥600 нм
δR, %
0
5
10
0
5
10
δk, %
2.4
3.3
5.0
4.7
5.5
7.4
δµ<sub>s</sub>, %
2.8
4.4
7.5
2.3
4.4
6.6
δg, %
2.5
7.0
12.8
2.1
3.5
6.2
Разработанный способ определения оптических и биофизических параметров биотканей и устройство для его реализации могут использоваться для поиска и диагностики злокачественных новообразований кожи (например, меланомы), а также в ходе проведения эндоспопических исследований слизистой оболочки полости рта, пищевода, органов желудочно-кишечного тракта и легких. Способ может также применяться для контроля газового состава крови в реанимации, токсикологии, при интенсивной терапии, для определения влияния на газовый состав гемоглобина факторов внешней среды (экологической обстановки, радиационного воздействия), для обнаружения внутреннего кровоизлияния, для непрерывного мониторинга уровня гемоглобина в крови при хирургических операциях и переливаниях крови, мониторинга желтухи у новорожденных, для определения состояния сосудистой системы и оценки жизнеспособности ткани при обморожениях и ожогах, а также в дерматологии и косметологии (диагностика и терапия таких изменений кожи как отек, красная волчанка, витилиго, солнечные лентинго и др.).
Разработанный способ позволяет в режиме реального времени осуществлять совместную интерпретацию спектральных и пространственных измерений диффузного отражения биотканей, что позволяет повысить точность определения как оптических, так и биофизических параметров ткани по сравнению со случаями раздельной интерпретации спектральной или пространственной составляющих информации. Повышается точность измерения данных параметров и за счет учета вариаций структурных и биохимических параметров биоткани, исключения использования априорной информации и устранения калибровочных измерений для нормированного r(L, λ) спектрально-пространственного профиля коэффициента диффузного отражения ткани
Литература.
1. Патент РФ №2234853, A61B 5/05, G01N 33/72; 27.08.2004.
2. Патент РФ №2337608, A61B 5/05, G01N 21/47; 10.11.2007.
3. Yudovsky, D. Rapid and accurate estimation of blood saturation, melanin content, and epidermis thickness from spectral diffuse reflectance / D. Yudovsky, L. Pilon // Appl. Opt. - 2010. - V.49, №10. - P.1707-1719.
4. Иванов, А.П. Спектры отражения света как средство диагностики структурных и биофизических параметров кожи / А.П. Иванов, В.В. Барун // Оптика и спектр. - 2008. - Т.104, №2. - С.344-351.
5. Bargo, P.R. In vivo determination of optical properties of normal and tumor tissue with white light reflectance and an empirical light transport model during endoscopy / P.R. Bargo [et. al.] // J. Biomed. Opt. - 2005. - Vol.10, №3. - P.034018-1-15.
6. Dam, J.S. Fiber-optic probe for noninvasive real-time determination of tissue optical properties at multiple wavelengths / J.S. Dam, C.B. Pedersen, T. Dalgaard, P.E. Fabricius, P. Aruna, S. Andersson-Engels // Appl. Opt. - 2001. - Vol.40 №7. P.1155-1164.
7. Kienle, A. Spatially resolved absolute diffuse reflectance measurements for noninvasive determination of the optical scattering and absorption coefficients of biological tissue / Kienle, L. Lilge, M.S. Patterson, R. Hibst, R. Steiner, B.C. Wilson // Appl. Opt. - 1996. - Vol.35, №15. P.2304-2314.
8. Междунар. заявка WO 2011068998/А2, G01N 33/49, G01N 33/72, G01N 21/47, A61B 5/1455; 09.06.2011.
9. Wang, L. MCML - Monte Carlo modeling of photon transport in multi-layered tissues / L. Wang, S.L. Jacques, L. Zheng // Computers Methods and Programs in Biomedicine. - 1995. - №47. - P.131-146.
10. Prahl, S.A. Optical absorption of hemoglobin [Электронный ресурс] / S.A. Prahl, - Режим доступа: http://omlc.ogi.edu/spectra/hemoglobin/index.html, свободный. - Название с экрана.
11. Барун, В.В. Оценка вклада локализованного поглощения света кровеносными сосудами в оптические свойства биологической ткани / В.В. Барун, А.П. Иванов // Оптика и спектр. - 2004. Vol.96, №6. P.1019-1024.
12. Saidi, I. Transcutaneous Optical Measurement of Hyperbilirubinemia in Neonates / I. Saidi // PhD thesis [Electronic resource]. - Rice University, - 1992. - Mode of access: http://scholarship.rice.edu/handle/1911/19082. - Date of access: 10.10.2011.
13. Jacques, S.L. Origins of tissue optical properties in the UVA, visible, and NIR regions / S.L. Jacques // Advances in Optical Imaging and Photon Migration. - 1996. - Vol.2. - P.364-369.
14. Jacques, S.L. Developing an optical fiber reflectance spectrometer to monitor bilirubinemia in neonates / S.L. Jacques, I. Saidi, A. Ladner, D. Oelberg // SPIE Proceedings. - 1997. - Vol.2975, №115. P.115-124.
15. Mourant, J.R. Mechanisms of light scattering from biological cells relevant to noninvasive optical-tissue diagnostics / J.R. Mourant [et. al.] // Appl. Opt. - 1998. - Vol.37, №16. - P.3586-3593.
16. Zonios, G. Diffuse Reflectance Spectroscopy of Human Adenomatous Colon Polyps In Vivo / G. Zonios [et. al.] // Appl. Opt. - 1999. - Vol.38, №31. - P.6628-6637.
17. Ghosh, N. Measurement of optical transport properties of normal and malignant human breast tissue / N. Ghosh [et. al.] // Appl. Opt. - 2001. - Vol.40, №1. - P.176-184.
18. Wei, H.J. Determination of optical properties of normal and adenomatous human colon tissues in vitro using integrating sphere techniques / H.J. Wei [et. al.] // World J. Gastroenterol. - 2005. - Vol.11, №16. - P.2413-2419.
24. Hidovic-Rowea, D. Modelling and va lidation of spectral reflectance for the colon / D. Hidovic-Rowea, E. Claridge // Phys. Med. Biol. - 2005. Vol.50. - P.1071-1093.
19. Salomatina E. Optical properties of normal and cancerous human skin in the visible and near-infrared spectral range / E. Salomatina [et. al.] // J. Biomed. Opt. - 2006. - Vol.11, №6. - P.064026-1-9.
20. Bashkatov, A.N. Optical properties of human stomach mucosa in the spectral range from 400 to 2000 nm: Prognosis for gastroenterology / A.N. Bashkatov [et. al.] // Medical Laser Application. - 2007. - Vol.22. - P.95-104.
21. Claridge, E. Quantifying mucosal blood volume fraction from multispectral images of the colon / E. Claridge, D. Hidovic-Rowe, P. Taniere, T. Ismail. // SPIE proceedings, 6511: Medical Imaging 2007: Physiology, Function, and Structure from Medical Images, (2007) 6511 ОС. DOI: 10.1117/12.709559.
22. Гираев, К.М. Спектры поглощения и рассеяния света тканями желудка при патологии / К.М. Гираев, Н.А. Ашурбеков, М.А. Лахина // Журн. прикл. спектр. - 2011. - Т.78, №1. - C.104-113.
23. Зуев, В.Е. Статистические модели температуры и газовых компонент земной атмосферы / В.Е. Зуев, B.C. Комаров. - Ленинград: Гидрометеоиздат., 1986. - С.47-50.
24. Veselovskii, I. Information content of multiwavelength lidar data with respect to microphysical particle properties derived from eigenvalue analysis /I. Veselovskii [et. al.] // Appl. Opt. - 2005. - Vol.44, №25. - P.5292-5303.