УСТРОЙСТВО ДЛЯ ОПРЕДЕЛЕНИЯ КОНЦЕНТРАЦИИ ГЕМОГЛОБИНА И СТЕПЕНИ ОКСИГЕНАЦИИ КРОВИ В СЛИЗИСТЫХ ОБОЛОЧКАХ
DEVICE FOR MEASURING HAEMOGLOBIN CONCENTRATION AND DEGREE OF BLOOD OXYGENATION IN MUCOUS MEMBRANES

Изобретение относится к области медицинского приборостроения.
Диагностика опухолевых заболеваний является одной из актуальнейших задач современной онкологии и хирургии. Правильный диагноз на ранних этапах патологического процесса в сочетании с эффективными методами лечения позволяет повысить вероятность благоприятного исхода лечения, увеличить выживаемость и улучшить качество жизни онкологических больных.
Из всех существующих методов ранней диагностики опухолей, поражающих слизистую оболочку полых органов, включая дыхательные пути и пищеварительный тракт, ведущим остается эндоскопическое исследование. Современные эндоскопические системы формируют высококачественное изображение слизистой оболочки в реальном времени непосредственно на дисплее компьютера. Однако изображения слизистой оболочки в естественном цвете не позволяют достигнуть высокой точности диагностики, в связи с различием характеристик изображения для различных анатомических участков, а также влиянием на них спектральных характеристик эндоскопической системы и неравномерной подсветки слизистой оболочки. Кроме того, обычное эндоскопическое исследование основывается на субъективной оценке изображения слизистой оболочки врачом, исходя из его опыта, квалификации и психоэмоционального состояния (при утомлении - цветовосприятие врача снижается, при возбуждении - усиливается) [1].
Во многих опухолях плотность кровеносных сосудов становится более высокой по сравнению с нормальной тканью, поэтому знание концентрации гемоглобина в ткани может помочь отличить доброкачественные новообразования от злокачественных. Известно также, что злокачественные новообразования характеризуются более интенсивными обменными процессами и более интенсивным потреблением кислорода. Поэтому концентрация кислорода в венулах, отводящих кровь из патологических регионов, ниже нормы. В связи с этим, информация о распределении концентрации гемоглобина в слизистой оболочке и степени оксигенации крови (доля окисленного гемоглобина в общем гемоглобине), выводимая на дисплее компьютера в ходе проведения эндоскопии, могла бы значительно повысить современные возможности ранней диагностики онкологических заболеваний.
Известно устройство для определения концентрации гемоглобина и степени оксигенации крови в слизистых оболочках [2], включающее широкополосный источник света, освещающее оптическое волокно, принимающее оптическое волокно, спектрометр и блок анализа спектра диффузного отражения ткани.
В данном устройстве освещающее и принимающее оптические волокна подводятся к исследуемому участку ткани через инструментальный канал эндоскопа. В связи с этим данное устройство не обладает оперативностью и его использование приводит к увеличению времени проведения болезненной для пациента процедуры эндоскопии. Контакт оптических волокон с исследуемой поверхностью приводит к искажению результатов измерений за счет изменений объемного коэффициента рассеяния ткани и ее микроциркуляции, возникающих даже при слабом нажиме на ткань. К недостаткам данного устройства относится также и сложность его калибровки.
Наиболее близким к заявляемому изобретению является устройство для определения содержания гемоглобина в слизистых оболочках [3], включающее эластичный зонд; широкополосный источник излучения; освещающее оптическое волокно, подводящее излучение к исследуемой ткани; блок регистрации обратного излучения ткани, включающий устройство захвата изображения и спектрометр; блок обработки, связанный с вышеупомянутой системой регистрации, выполняющий измерение и обработку спектра обратного излучения ткани.
Главным недостатком данного устройства является локальность спектральных измерений обратного рассеяния ткани, выполняемых на его основе. Данные измерения позволяют судить о состоянии ткани только в одной ее точке, соответствующей нескольким пикселям изображения ткани. Между тем для более точного диагноза требуется выявлять распределение гемоглобина и оксигенации крови в пораженных участках ткани и в их окрестностях.
Предлагаемое изобретение направлено на решение задачи расширения функциональных возможностей устройства за счет получения количественных значений гемоглобина и оксигенации крови в пораженных участках ткани и в их окрестностях, повышения точности диагностики онкологических заболеваний слизистых оболочек.
Для решения данной задачи устройство для определения концентрации гемоглобина и степени оксигенации крови в слизистых оболочках, включающее источник излучения, освещающее оптическое волокно, эластичный зонд, блок регистрации изображения, блок обработки изображения, дополнительно содержит блоки управления излучателями и распределения каналов посылки излучения, два взаимно ортогональных поляризационных фильтра, цифровой кабель, причем источник излучения выполнен из набора излучателей на разных длинах волн и связан с микропроцессорным блоком управления излучателями и блоком распределения каналов посылки излучения, выход которого соединен со входом освещающего оптического волокна, расположенного в эластичном зонде, содержащем наконечник, в котором расположены два взаимно ортогональных поляризационных фильтра, один из которых связан с выходом освещающего волокна, второй - с блоком регистрации изображения, представляющим собой монохромную ПЗС-матрицу с установленной перед ней собирающей линзой, выход которого соединен цифровым кабелем, расположенным в зонде, - с блоком обработки изображения, определяющим количественные значения концентрации гемоглобина и степени оксигенации крови во всех точках изображения слизистой оболочки, получаемого на ПЗС-матрице.
Источник излучения может быть выполнен на основе широкополосного излучателя, а блок регистрации изображения - на основе цветной ПЗС-матрицы, позволяющей получать изображения в трех и более спектральных участках.
Сущность данного изобретения поясняется с помощью фиг.1-8.
Фиг.1 - Функциональная схема предлагаемого устройства.
Фиг.2 - Нормированные спектры диффузного отражения слизистой оболочки легких в норме, рассчитанные в рамках предложенной модели (кривые), и сравнение с экспериментальными данными (окружности).
Фиг 3 - Нормированные спектры диффузного отражения слизистой оболочки легких в норме при онкологическом заболевании, рассчитанные в рамках предложенной модели (кривые), и сравнение с экспериментальными данными (окружности).
k
Фиг.4 - Сопоставление заданных (ось абсцисс) и восстановленных из сигналов r(λ) (ось ординат) значений концентрации гемоглобина в слизистой оболочке.
k
Фиг.5 - Сопоставление заданных (ось абсцисс) и восстановленных из сигналов r(λ) (ось ординат) значений степени оксигенации крови в слизистой оболочке.
k
Фиг.6 - Используемые при моделировании профили относительной спектральной чувствительности цветной ПЗС-матрицы S(λ) (кривые 1, 2, 3 для k=В, G, R соответственно) и мощности широкополосного излучателя P(λ) (кривая 4).
BR
GR
Фиг.7 - Смоделированная совокупность нормированных RGB-сигналов (r, r) и соответствующих им концентраций гемоглобина.
BR
GR
Фиг.8 - Сопоставление заданных (ось абсцисс) и восстановленных из сигналов r и r (ось ординат) значений концентрации гемоглобина в слизистой оболочке.
Функциональная схема предлагаемого устройства изображена на фиг.1, где 1 - микропроцессорный блок управления излучателями, 2 - источник излучения, 3 - блок распределения каналов посылки излучения, 4 - освещающее оптическое волокно, 5 - эластичный зонд, 6 - наконечник зонда, 7 - исследуемая ткань, 8 и 9 - взаимно ортогональные поляризационные фильтры, 10 - собирающая линза, 11 - ПЗС-матрица, 72 - цифровой кабель, 13 - блок обработки изображения ткани.
Данное устройство функционирует следующим образом. Микропроцессорный блок управления излучателями 1 организует чередующуюся посылку источником 2 излучения с различной длиной волны. Источник излучения 2 состоит из набора излучателей на различных длинах волн, например лазерных источников или светодиодов. Блок распределения каналов посылки излучения 3 подводит излучение от каждого излучателя к входу освещающего оптического волокна 4. В простейшем случае в качестве блока 3 может использоваться разветвленное оптическое волокно, в котором несколько волокон, каждое из которых соединено с одним из излучателей, сводятся в одно общее волокно 4. Потоки излучения на различных длинах волн направляются на исследуемую ткань 7 по оптическому волокну 4, находящемуся внутри эластичного зонда 5 для эндоскопии внутри полости тела. Отраженное от ткани излучение собирается в блоке регистрации изображения, выполненном на основе ПЗС-матрицы 11 с расположенной перед ней линзой 10. При смене спектрального состава освещающего излучения на выходе блока регистрации последовательно формируются спектральные изображения ткани в отраженном от нее свете с различными длинами волн. Для исключения попадания на светочувствительные элементы матрицы зеркально отраженного света и устранения бликов в изображении ткани, в наконечнике 6 зонда 5 установлены два взаимно ортогональных поляризационных фильтра 8 и 9, один - на выходе освещающего оптоволокна 4, второй - перед линзой 10. Естественное излучение, прошедшее через поляризационный фильтр 8, становится поляризованным. Поскольку излучение, отражаемое поверхностью ткани, сохраняет исходную поляризацию, использование скрещенных поляризаторов позволяет блокировать эту паразитную составляющую. Часть излучения, проникшая вглубь ткани, практически полностью утрачивает свою первоначальную поляризацию вследствие многократного рассеяния. При этом рассеянное обратно излучение частично проходит через поляризационный фильтр 9 и попадает на блок регистрации изображения. Сигналы от светочувствительных элементов матрицы по цифровому кабелю 12, расположенному в зонде, передаются на блок обработки изображения 13, выполняющий анализ последовательности спектральных изображений исследуемой ткани и определяющий количественные значения концентрации гемоглобина и степени оксигенации крови во всех точках изображения ткани в соответствии с нижеприведенным алгоритмом. Блок обработки изображения 13 может быть компьютером или микропроцессором, способным выдавать свои результаты любому внешнему устройству для их визуализации.
В наиболее простом варианте реализации вышеописанного устройства источник излучения может быть выполнен на основе широкополосного излучателя, такого как белый светодиод или лампа. При этом спектральную селекцию излучения, отраженного от исследуемой ткани, можно осуществлять при помощи ПЗС-матрицы, позволяющей получать изображения в нескольких спектральных участках, например в синем, зеленом и красном. С точки зрения спектральных измерений это равносильно использованию в устройстве (например, на выходе широкополосного излучателя) нескольких фильтров, выделяющих из белого света спектральные составляющие.
Алгоритм количественной интерпретации спектральных изображений слизистой оболочки. Спектральные сигналы на выходе ПЗС-матрицы (спектральные изображения ткани) определяются в соответствии с формулой
<mrow><mi>V</mi><mo>=</mo><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><msub><mi>Λ</mi><mi>k</mi></msub></mrow><mo>)</mo></mrow><mo>=</mo><mi>G</mi><mrow><mo>(</mo><mrow><mi>θ</mi><mo>,</mo><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo></mrow><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mrow><msub><mi>λ</mi><mi>l</mi></msub></mrow><mrow><msub><mi>λ</mi><mi>r</mi></msub></mrow></munderover><mrow><msub><mi>P</mi><mi>k</mi></msub></mrow></mrow></mstyle><mrow><mo>(</mo><mi>λ</mi><mo>)</mo></mrow><mi>τ</mi><mrow><mo>(</mo><mi>λ</mi><mo>)</mo></mrow><mi>S</mi><mrow><mo>(</mo><mi>λ</mi><mo>)</mo></mrow><mi>R</mi><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>λ</mi></mrow><mo>)</mo></mrow><mi>d</mi><mi>λ</mi><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo stretchy="false">(</mo><mn>1</mn><mo stretchy="false">)</mo></mrow>
k
k
l
r
где Λ - спектральный участок, соответствующий полосе излучения k-го излучателя; P(λ) - спектральная мощность k-го излучателя; х, y - координаты элемента поверхности ткани; G(θ, х, y) - геометрический фактор, зависящий от пространственного распределения освещенности ткани и угла 6 между нормалью к поверхности ткани и осью собирающей линзы; R(x, y, λ) -спектральный коэффициент диффузного отражения ткани, т.е. отношение потока излучения с длиной волны λ, отраженного элементарной площадкой с центром в точке (х, y), к падающему на нее потоку; τ(λ) - спектральное пропускание оптического волокна и собирающей линзы; S(λ) - спектральная чувствительность матрицы; λ и λ - левая и правая границы диапазона спектра, используемого блоком регистрации.
k
k
k
0
k
0
k
Для устранения влияния на регистрируемые спектральные изображения V(x, y, Λ) неравномерной освещенности ткани и ее положения относительно блока регистрации будем оперировать с изображениями r(x, y, Λ)=V(x, y, Λ)/V(x, y, Λ), получаемыми путем нормировки изображений V(x, y, Λ) на опорный спектральный слой V(x, y, Λ). Карты двумерных распределений параметров слизистой оболочки получаются в результате анализа спектральных сигналов для каждого пикселя изображения r(x, y, Λ).
Учитывая высокое пространственное разрешение современных ПЗС-матриц, использование традиционных методов решения обратных задач, основанных на многократном решении уравнения переноса излучения в исследуемой среде при различных значениях ее параметров, не позволяет получить необходимую для практики оперативность построения искомых параметрических карт. Кроме того, для использования данных методов необходимо, чтобы количество спектральных измерений значительно превышало количество неизвестных параметров среды, что является трудно достижимым в условиях проведения эндоскопии.
k
i
i
k
Для оперативного решения обратных задач оптики светорассеивающих сред при малом объемом измеряемой информации наиболее подходящим является метод множественных регрессий. Данный метод основан на построении явного аналитического выражения F, связывающего измеряемые оптические характеристики объекта s=(S) с его искомыми параметрами p=(р), как p=F(s,a), где а - вектор параметров аналитического выражения. Для построения оператора F используется физико-математическая модель исследуемого объекта, позволяющая по заданным значениям вектора параметров p рассчитывать вектор результатов измерений s. На основе этой модели формируется «обучающий» ансамбль реализации p и s и выполняется его регрессионный анализ для нахождения параметров а выбранного аналитического выражения. Простейшим аналитическим выражением, позволяющим предсказывать нужный параметр p по измеренным сигналам S=(k=1, …, K), является множественная линейная регрессия:
<mrow><msub><mi>p</mi><mi>i</mi></msub><mo>=</mo><msub><mi>a</mi><mi>i</mi></msub><mo>+</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><mrow><msub><mi>a</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>⋅</mo><msub><mi>S</mi><mi>k</mi></msub><mo>.</mo></mrow></mstyle><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo stretchy="false">(</mo><mn>2</mn><mo stretchy="false">)</mo></mrow>
j
k
В действительности параметры объекта и измеряемые сигналы редко связаны линейной зависимостью, поэтому регрессия (2), как правило, не позволяет получить решение обратной задачи с приемлемой точностью. В тоже время в большинстве случаев регрессионная связь между величинами p и S хорошо описывается полиномами с заданной степенью N:
<mrow><msub><mi>p</mi><mi>i</mi></msub><mo>=</mo><msub><mi>a</mi><mi>i</mi></msub><mo>+</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><msub><mi>n</mi><mn>1</mn></msub></mrow></munder><mrow><mo>…</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><msub><mi>n</mi><mi>K</mi></msub></mrow></munder><mrow><msub><mi>a</mi><mrow><msub><mi>n</mi><mn>1</mn></msub><mo>…</mo><msub><mi>n</mi><mi>K</mi></msub><mo>,</mo><mi>i</mi></mrow></msub><mo>⋅</mo><msubsup><mi>S</mi><mn>1</mn><mrow><msub><mi>n</mi><mi>l</mi></msub></mrow></msubsup><mo>…</mo><msubsup><mi>S</mi><mi>K</mi><mrow><msub><mi>n</mi><mi>K</mi></msub></mrow></msubsup></mrow></mstyle></mrow></mstyle><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo stretchy="false">(</mo><mn>3</mn><mo stretchy="false">)</mo></mrow>
<mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>Л</mi></munderover><mrow><msub><mi>n</mi><mi>i</mi></msub><mo>≤</mo><mi>N</mi><mo>.</mo></mrow></mstyle></mrow>
где
tHb
k
1
tHb
2
k
k
k
j
tHb
k
k
tHb
k
k
<mrow><msub><mi>a</mi><mrow><msub><mi>k</mi><mn>1</mn></msub><mo>…</mo><msub><mi>k</mi><mi>N</mi></msub><mo>,</mo><mi>j</mi></mrow></msub></mrow>
Для решения задачи оперативной реконструкции двумерных распределений концентрации гемоглобина F(x,y) и степени оксигенации крови S(x,y) в слизистой оболочке по ее по нормированным спектральным изображениям r(x, y, Λ) также можно использовать полиномиальные множественные регрессии вида (3), в которых p=F(x, y), p=S(x, y) и S=r(x, y, Λ). Поскольку изображения r(х, y, Λ) свободны от влияния пространственного распределения освещенности ткани и ее положения относительно блока регистрации, то коэффициенты a и
не зависят от геометрических координат x и y. В связи с этим далее зависимость величин F, S и r от координат x, y для краткости будем опускать, подчеркивая тем самым единообразие алгоритма количественной интерпретации спектральных сигналов r(Λ) для всех точек зарегистрированного изображения слизистой оболочки. Численные значения коэффициентов регрессий будем определять с использованием «обучающего» ансамбля реализации r(Λ), F и S, соответствующего аппаратурным функциям конкретной эндоскопической системы. Данный ансамбль рассчитывается на основе модели переноса света в биоткани при широкой вариации ее параметров. Таким образом, априорная информация о вариациях структурных и биохимических параметров ткани в предлагаемом методе количественной интерпретации ее мультиспектральных изображений вносится на этапе численного моделирования сигналов r(Λ) и отражается в значениях коэффициентов множественных регрессий между r(Λ) и искомыми параметрами ткани.
s
Оптическая модель слизистой оболочки. Известно, что основными параметрами, характеризующими распространение оптического излучения в рассеивающей среде, являются коэффициенты поглощения k, рассеяния β и индикатриса рассеяния или ее средний косинус g. Причем, в оптике биологических тканей, как очень мутных сред, для описания световых полей достаточно знать не сами величины β и g, а их комбинацию - транспортный коэффициент рассеяния µ=β(1-g).
Спектр объемного коэффициента поглощения слизистой оболочки k(λ) моделируется в виде линейной комбинации спектров поглощения ее хромофоров. В видимой области основными поглощающими хромофорами биоткани являются оксигемоглобин и деоксигемоглобин, поэтому выражение для объемного коэффициента поглощения слизистой оболочки имеет вид:
<mrow><mi>k</mi><mrow><mo>(</mo><mi>λ</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>k</mi><mi>t</mi></msub><mrow><mo>(</mo><mi>λ</mi><mo>)</mo></mrow><mo>+</mo><mi>α</mi><mrow><mo>(</mo><mrow><msub><mi>D</mi><mi>V</mi></msub><mo>,</mo><mi>λ</mi></mrow><mo>)</mo></mrow><mfrac><mrow><msub><mi>F</mi><mrow><mi>t</mi><mi>H</mi><mi>b</mi></mrow></msub></mrow><mrow><msub><mi>M</mi><mrow><mi>H</mi><mi>b</mi></mrow></msub></mrow></mfrac><mn>1</mn><mi>n</mi><mrow><mo>(</mo><mrow><mn>10</mn></mrow><mo>)</mo></mrow><mrow><mo>[</mo><mrow><mi>S</mi><msub><mi>ε</mi><mrow><mi>H</mi><mi>b</mi><msub><mi>O</mi><mn>2</mn></msub></mrow></msub><mrow><mo>(</mo><mi>λ</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><mi>S</mi></mrow><mo>)</mo></mrow><msub><mi>ε</mi><mrow><mi>H</mi><mi>b</mi></mrow></msub><mrow><mo>(</mo><mi>λ</mi><mo>)</mo></mrow></mrow><mo>]</mo></mrow><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo stretchy="false">(</mo><mn>4</mn><mo stretchy="false">)</mo></mrow>
t
0
0
t
tHb
Hb
Hb
V
V
-1
-1
<mrow><msub><mi>ε</mi><mrow><mi>H</mi><mi>b</mi><msub><mi>O</mi><mn>2</mn></msub></mrow></msub></mrow>
где k(λ)=A·exp[-B(λ-λ)] - объемный коэффициент поглощения обескровленной ткани, λ=600 нм, А (мм) и В (нм) - параметры, характеризующие спектральную зависимость k; F (г/л) - концентрация общего гемоглобина в биологической ткани; M=64500 г/моль - молярная масса гемоглобина;
и ε - молярные коэффициенты поглощения оксигемоглобина и деоксигемоглобина соответственно [4]; S - степень оксигенации крови; α(D,λ) - поправочный коэффициент, учитывающий эффект локализованного поглощения света кровеносными сосудами [5]; D - средний диаметр кровеносных сосудов.
Согласно данным [2, 6-8] транспортный коэффициент рассеяния биологических тканей в видимой области спектра с хорошей точностью аппроксимируется степенной зависимостью
<mrow><msub><mi>μ</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>λ</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>C</mi><mi>s</mi></msub><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><msub><mi>λ</mi><mn>0</mn></msub></mrow><mi>λ</mi></mfrac></mrow><mo>)</mo></mrow></mrow><mi>ν</mi></msup><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo stretchy="false">(</mo><mn>5</mn><mo stretchy="false">)</mo></mrow>
0
s
где λ=632 нм; C и ν - структурные параметры ткани, характеризующие объемное содержание и размер ее «эффективных» рассеивателей. Для описания процесса рассеяния излучения используется однопараметрическая функция Хеньи-Гринштейна [9] с фактором анизотропии рассеяния g.
t
s
tHb
s
s
-1
С учетом показателя преломления ткани и; оптическая модель слизистой оболочки определяется 9 параметрами: n, g, А, В, C, ν, F, D, S. В табл.1 приведены диапазоны их вариаций, выбранные путем критического анализа результатов различных авторов для нормальных и опухолевых тканей слизистых оболочек полости рта, пищевода, органов желудочно-кишечного тракта и легких. Модельные параметры могут варьироваться независимо друг от друга, но при каждой их комбинации спектральные значения оптических параметров ткани k(λ), µ(λ) не должны выходить из диапазонов, наблюдаемых в эксперименте. Исходя из вышеуказанных литературных данных, выбраны следующие ограничения: 1) 0.04≤k(632 нм)≤0.5 мм; 2) 1.0≤Λ(632 нм)≤50, где Λ=µ/k.
3
s
t
Путем независимых вариаций модельных параметров с учетом выполнения вышеуказанных ограничений, сформирован ансамбль из 10 реализаций k(λ), µ(λ), g и n. Для каждой реализации методом Монте-Карло [9] рассчитан коэффициент диффузного отражения R(λ) на 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 нм), характеризующих основные особенности спектральных кривых поглощения хромофоров биологической ткани в видимой области спектра. Расчет R(λ) выполнялся для однородного фантома ткани толщиной 4 мм с использованием 500000 фотонов.
Табл.1. Диапазоны вариаций параметров оптической модели слизистой оболочки
Параметр
Диапазон
n<sub>t</sub>
1.35-1.45
g
0.5-0.95
А
0.01-0.10 мм<sup>-1</sup>
В
0.001-0.028 нм<sup>-1</sup>
C<sub>s</sub>
1.5-2.5 мм<sup>-1</sup>
ν
1.0-2.5
F<sub>tHb</sub>
0.5-40 г/литр
D<sub>V</sub>
5-30 мкм
S
40-99%
s
Следует отметить, что использование метода Монте-Карло для расчета коэффициентов R(λ) требует больших затрат машинного времени, что является особенно критичным при большой выборке реализации R(λ). В связи с этим нами разработан простой аналитический метод расчета R(λ), позволяющий в явном виде рассчитать коэффициент диффузного отражения среды с известными значениями показателя преломления n, коэффициента поглощения k и транспортного коэффициента рассеяния µ. Для этого используется выражение
<mrow><mi>R</mi><mo>=</mo><mrow><mo>(</mo><mrow><mi>n</mi><mo>,</mo><mi>k</mi><mo>,</mo><msub><mi>μ</mi><mi>s</mi></msub></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><msub><mi>r</mi><mi>c</mi></msub></mrow><mo>)</mo></mrow><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><msub><mi>r</mi><mi>d</mi></msub></mrow><mo>)</mo></mrow><msub><mi>R</mi><mi>c</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mo>,</mo><msub><mi>μ</mi><mi>s</mi></msub></mrow><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>−</mo><msub><mi>r</mi><mi>d</mi></msub><msub><mi>R</mi><mi>d</mi></msub><mrow><mo>(</mo><mrow><mi>k</mi><mo>,</mo><msub><mi>μ</mi><mi>s</mi></msub></mrow><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo stretchy="false">(</mo><mn>6</mn><mo stretchy="false">)</mo></mrow>
c
d
c
d
c
d
s
где r - коэффициент отражения нормально падающего излучения от поверхности среды; r - коэффициент отражения поверхности среды, освещенной диффузным потоком изнутри [10]; R и R - коэффициенты отражения среды без внешней границы при направленном и диффузном освещении соответственно. Коэффициенты R и R рассчитываются при помощи аналитических выражений, связывающих их с оптическими параметрами среды k и µ:
<mrow><msub><mi>R</mi><mrow><mi>c</mi><mo>,</mo><mi>d</mi></mrow></msub><mrow><mo>(</mo><mrow><mi>k</mi><mo>,</mo><msub><mi>μ</mi><mi>s</mi></msub></mrow><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mrow><mo>−</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>m</mi><mo>−</mo><mn>0</mn></mrow><mn>3</mn></munderover><mrow><msup><mi>δ</mi><mi>m</mi></msup></mrow></mstyle><mrow><mo>(</mo><mrow><msub><mi>A</mi><mi>m</mi></msub><msup><mi>k</mi><mi>m</mi></msup><mo>+</mo><msub><mi>B</mi><mi>m</mi></msub><msubsup><mi>μ</mi><mi>s</mi><mi>m</mi></msubsup></mrow><mo>)</mo></mrow><mo>−</mo><mi>k</mi><mi>δ</mi><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>m</mi><mo>−</mo><mn>0</mn></mrow><mn>3</mn></munderover><mrow><msub><mi>C</mi><mi>m</mi></msub></mrow></mstyle><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><msub><mi>μ</mi><mi>s</mi></msub></mrow><mi>k</mi></mfrac></mrow><mo>)</mo></mrow></mrow><mi>m</mi></msup></mrow><mo>)</mo></mrow><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo stretchy="false">(</mo><mn>7</mn><mo stretchy="false">)</mo></mrow>
<mrow><mi>δ</mi><mo>=</mo><mrow><mn>1</mn><mo>/</mo><mrow><msqrt><mrow><mn>3</mn><mi>k</mi><mrow><mo>(</mo><mrow><mi>k</mi><mo>+</mo><msub><mi>μ</mi><mi>s</mi></msub></mrow><mo>)</mo></mrow></mrow></msqrt></mrow></mrow></mrow>
m
m
m
s
s
s
s
-1
-1
-1
-1
где
. Значения коэффициентов A, B и C, получены путем аппроксимации расчетов R(n, k, µ) методом Монте-Карло и приведены в табл.2. Диапазоны значений параметров n, k, µ для которых правомерна аппроксимация (7) характерны для большинства биологических тканей в видимой области спектра: 0.01≤k≤8 мм, 1.0≤µ≤5 мм, 0.3≤µ/k≤110 мм, 1.35≤n≤1.45 мм. Средняя погрешность расчета коэффициента диффузного отражения по формулам (6), (7) относительно методом Монте-Карло составляет 1%, максимальная не превышает 6%.
Табл.2. Коэффициенты A<sub>m</sub>, B<sub>m</sub> и C<sub>m</sub> выражения (7)
m
коэффиценты для R<sub>c</sub>
коэффциенты для R<sub>d</sub>
A<sub>m</sub>
B<sub>m</sub>
C<sub>m</sub>
A<sub>m</sub>
B<sub>m</sub>
C<sub>m</sub>
0
4.0190
-2.9329
0.0000
1.8691
-1.1374
0.0000
1
3.9648
0.1549
0.9697
1.8515
0.3382
0.3684
2
1.5931
-0.6950
-0.6948
0.9910
-0.8309
-0.8309
3
0.3371
0.3371
1.5931
0.9704
0.9704
0.9911
k
0
0
0
t
s
tHb
V
0
-2
-1
-2
-1
-1
Для подтверждения корректности используемой оптической модели слизистой оболочки и разработанного метода расчета коэффициента диффузного отражения проведено сравнение результатов моделирования спектра R(λ) с экспериментальными данными. Для этого использовались результаты in-vivo измерений спектра диффузного отражения нормальной и раковой ткани легких, проводимых в работе [3] на базе широкополосного источника света, собирающего оптического волокна и спектрометра. Поскольку в условиях проведения эндоскопии практически невозможно соблюсти идентичные условия посылки и регистрации излучения для исследуемой ткани и калибровочного образца, то абсолютные измерения R(λ) неизбежно отягощены влиянием мультипликативных геометрических параметров и использование таких измерений для валидации модельных расчетов R(λ) не является корректным. Однако для решаемой задачи абсолютные значения R(λ) не представляют особого интереса, поскольку нормированные изображения r(x, y, Λ) зависят лишь от относительной спектрального хода R(λ). В связи с этим для валидации модельных расчетов использовались нормированные экспериментальные зависимости R(λ)/R(λ), где λ=700 нм. Данные зависимости для двух анатомических участков слизистой оболочки легких изображены на фиг.2, 3. Там же приведены результаты расчета R(λ)/R(λ) вышеописанным методом при следующих значениях модельных параметров (числитель - нормальная ткань, знаменатель - опухоль): n=1.35/1.35; А=3.2/5.2(10 мм), B=2.1/1.9(10 нм), C=1.72/1.63 (мм), ν=1.0/1.0, F=4.05/10.8 (г/л), D=5.0/5.0 (мкм), S=0.78/0.98. Как видно, используемый аналитический метод расчета R(λ) позволяет достаточно хорошо (в пределах погрешностей измерения и оцифровки экспериментальных данных) воспроизвести экспериментальные зависимости R(λ)/R(λ). К сожалению, точные значения параметров ткани, соответствующие экспериментальным данным, не известны, поэтому сделать вывод о точности их количественных оценок не представляется возможным. Тем не менее, полученные значения параметров характеризуют структуру и биохимический состав ткани в количественном аспекте и, вне всякого сомнения, могут использоваться для количественной интерпретации мультиспектральных эндоскопических изображений. Так, в частности, найденные значения модельных параметров демонстрируют существенные различия содержания гемоглобина и степени оксигенации крови в нормальной и опухолевой тканях, возможность выявления которых является ключевым моментом при ранней диагностике злокачественных новообразований.
k
k
k
k
k
k
0
k
k
0
k
k
Выбор оптимальных спектральных участков. На основе вышеописанной модели сформирован набор «обучающий» данных по параметрам слизистой оболочки и спектральным коэффициентам ее диффузного отражения R(λ) (26 спектральных значений). Расчет соответствующих им сигналов r(Λ) требует знания аппаратурной функций A(λ), учитывающей мощности излучателей, чувствительности матрицы и пропускания оптики. При известной функции A(λ) коэффициенты аналитических выражений, связывающих сигналы r(Λ) с искомыми параметрами слизистой оболочки, могут быть легко определены методом наименьших квадратов. Однако перед этим необходимо определить оптимальное количество спектральных измерений r(Λ), позволяющих определять искомые параметры ткани практически с такой же точностью, как и на основании анализа всего спектра ее диффузного отражения. Для этого проведем анализ линейно-независимых компонент, содержащихся в сигналах r(λ)=A(λ)R(λ)/R(λ) соответствующих узким спектральным участкам Λ≡λ где k=1, …, 25; λ=700 нм. Значения λ в порядке возрастания индекса k указаны выше. Функцию A(λ) без ограничения общности можно положить равной единице.
i
i
k
k
0
k
2
<mrow><mrow><mo>(</mo><mrow><msqrt><mrow><msub><mi>l</mi><mn>5</mn></msub></mrow></msqrt><mo>≈</mo><mn>2.6</mn><mi>%</mi></mrow><mo>)</mo></mrow></mrow>
Известно, что системой ортонормированных векторов, позволяющих оптимальным образом аппроксимировать любую реализацию случайного вектора x, является система собственных векторов его ковариационной матрицы [11]. При этом отношение суммы соответствующих собственных чисел l к следу матрицы (суммарной дисперсии компонент x) характеризует точность этой аппроксимации и определяет меру значимости процессов, описываемых этими собственными векторами в общей изменчивости вектора х. Число линейно-независимых компонент x определяется номером собственного числа, для которого l>δx, где δx - погрешность измерений x [12]. Результаты анализа собственных чисел ковариационной матрицы r(λ) показывают, что при δr>2% сигналы r(λ) содержат не более пяти линейно-независимых компонент
. Следовательно, использование шести спектральных участков (с учетом нормировочной длины волны λ) позволяет извлечь большую часть информации, содержащейся в r(λ).
k
tHb
k
k
Оптимальные значения λ выбирались путем компьютерного перебора различных комбинаций из вышеуказанного набора длин волн и оценки соответствующих им погрешностей восстановления параметров F и S. Для каждой комбинации λ методом наименьших квадратов определялись коэффициенты регрессии между рассматриваемыми параметрами и r(λ). Использовалась следующее уравнения регрессии:
<mrow><mi>y</mi><mo>=</mo><msub><mi>a</mi><mrow><mn>00</mn></mrow></msub><mo>+</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover><mrow><msub><mi>a</mi><mrow><mi>k</mi><mi>m</mi></mrow></msub><msup><mrow><mrow><mo>[</mo><mrow><mn>1</mn><mi>n</mi><mi>r</mi><mrow><mo>(</mo><mrow><msub><mi>λ</mi><mi>k</mi></msub></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow></mrow><mi>m</mi></msup></mrow></mstyle></mrow></mstyle><mo>.</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo stretchy="false">(</mo><mn>8</mn><mo stretchy="false">)</mo></mrow>
tHb
k
tHb
k
tHb
k
k
где К - количество спектральных слоев в мультиспектральном изображении ткани, y - параметр ткани (F или S). Далее перебирались все смоделированные реализации r(λ) и для каждой из них вычислялись параметры F и S с использованием регрессий (8) в которых на значения r(λ) накладывались случайные отклонения в пределах 1% (для оценки чувствительности регрессий к погрешностям оптических измерения). Полученные значения параметров F и S сравнивались с их истинными значениями и вычислялись погрешности восстановления искомых параметров. Оптимальной комбинацией λ принята та, которой соответствует наибольшая информативность сигналов r(λ) относительно искомых параметров. В качестве меры информативности бралась величина
<mrow><mi>Ф</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mi>δ</mi><msubsup><mi>F</mi><mrow><mi>t</mi><mi>H</mi><mi>b</mi></mrow><mi>n</mi></msubsup></mrow></mstyle></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mfrac><mrow><mrow><mo>|</mo><mrow><msubsup><mi>F</mi><mrow><mi>t</mi><mi>H</mi><mi>b</mi></mrow><mi>n</mi></msubsup><mo>−</mo><msubsup><mover accent="true"><mi>F</mi><mo>¯</mo></mover><mrow><mi>t</mi><mi>H</mi><mi>b</mi></mrow><mi>n</mi></msubsup></mrow><mo>|</mo></mrow></mrow><mrow><msubsup><mi>F</mi><mrow><mi>t</mi><mi>H</mi><mi>b</mi></mrow><mi>n</mi></msubsup></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mi>δ</mi><msup><mi>S</mi><mi>n</mi></msup></mrow></mstyle></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mfrac><mrow><mrow><mo>|</mo><mrow><msup><mi>S</mi><mi>n</mi></msup><mo>−</mo><mover accent="true"><mi>S</mi><mo>¯</mo></mover></mrow><mo>|</mo></mrow></mrow><mrow><msup><mi>S</mi><mi>n</mi></msup></mrow></mfrac></mrow></mstyle></mrow></mstyle><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo stretchy="false">(</mo><mn>9</mn><mo stretchy="false">)</mo></mrow>
<mrow><mi>δ</mi><msubsup><mi>F</mi><mrow><mi>t</mi><mi>H</mi><mi>b</mi></mrow><mi>n</mi></msubsup></mrow>
<mrow><msub><mover accent="true"><mi>F</mi><mo>¯</mo></mover><mrow><mi>t</mi><mi>H</mi><mi>b</mi></mrow></msub></mrow>
<mover accent="true"><mi>S</mi><mo>¯</mo></mover>
<mrow><msubsup><mi>F</mi><mrow><mi>t</mi><mi>H</mi><mi>b</mi></mrow><mo>*</mo></msubsup></mrow>
n
tHb
tHb
k
tHb
k
tHb
где индекс n обозначает номер реализации модельных параметров; N - полное количество реализации;
и δS - относительные погрешности восстановления параметров F и S на основе регрессий (8);
и
- средние значения параметров F и S. Каждое слагаемое в выражении (9) представляет собой отношение априорной и апостериорной неопределенностей искомого параметра. Для рассматриваемого набора длин волн максимум Ф достигается при следующих значениях λ: 450, 480, 580, 594, 610 нм. На фиг.4, 5 представлены результаты сопоставления известных значений параметров F и S со значениями
и S*, восстановленными из сигналов r(λ) по формуле (8) в условиях общей вариативности всех других модельных параметров. Средние погрешности восстановления параметров F и S, соответствующие данной комбинации длин волн, составляют 12.3 и 1.9% соответственно. Следует отметить, что возможны и другие комбинации длин волн, обеспечивающие примерно такую же точность определения рассматриваемых параметров.
Таким образом, обработка спектральных изображений ткани по формулам регрессий (3) или (8) включает лишь использование простейших арифметических операций и может быть легко запрограммирована на микропроцессоре. При этом не требуется решение уравнения переноса излучения и использование сложных математических алгоритмов решения некорректных обратных задач. Для восстановления параметров ткани не требуется также и использование априорной информации или предположений о фиксированном значении одного или нескольких параметров ткани, поскольку ансамбль реализации модельных параметров полностью охватывает область реальных значений параметров ткани.
Количественная интерпретация цветных изображений слизистой оболочки. Наиболее простой вариант реализации заявляемого устройства может быть основан на использовании широкополосного источника света и ПЗС-матрицы, позволяющей получать изображения в нескольких спектральных участках, например в синем, зеленом и красном (цветной RGB матрицы). Сигналы на выходе RGB матрицы определяются в соответствии с формулой
<mrow><mi>V</mi><mo>=</mo><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><msub><mi>Λ</mi><mi>k</mi></msub></mrow><mo>)</mo></mrow><mo>=</mo><mi>G</mi><mrow><mo>(</mo><mrow><mi>θ</mi><mo>,</mo><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo></mrow><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mrow><mn>400</mn></mrow><mrow><mn>700</mn></mrow></munderover><mrow><mi>P</mi><mrow><mo>(</mo><mi>λ</mi><mo>)</mo></mrow><mi>R</mi><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>λ</mi></mrow><mo>)</mo></mrow><mi>τ</mi><mrow><mo>(</mo><mi>λ</mi><mo>)</mo></mrow><msub><mi>S</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>λ</mi><mo>)</mo></mrow><mi>d</mi><mi>λ</mi><mo>,</mo></mrow></mrow></mstyle><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo stretchy="false">(</mo><mn>10</mn><mo stretchy="false">)</mo></mrow>
k
где P(λ) - спектральная мощность источника света; x, y - координаты элемента поверхности ткани; G(θ, x, y) - геометрический фактор, зависящий от пространственного распределения освещенности ткани и угла θ между нормалью к поверхности ткани и осью собирающей линзы; R(x, y, λ) - спектральный коэффициент диффузного отражения ткани; τ(λ) - спектральное пропускание оптического волокна и собирающей линзы; S(λ) - спектральная чувствительность матрицы в красном (R), зеленом (G) и синем (B) спектральном участках; k=R, G, В.
k
Рассмотрим возможность использования сигналов (10) для определения концентрации гемоглобина в слизистой оболочке. Рассмотрение проведем на примере устройства со спектральными характеристиками Р(λ) и S(λ), приведенными на фиг.6. Функцию пропускания τ(λ) в диапазоне λ=400-700 нм будем полагать равной единице.
<mrow><msub><mi>r</mi><mrow><mi>B</mi><mi>R</mi></mrow></msub><mo>=</mo><mrow><mrow><mi>V</mi><mrow><mo>(</mo><mrow><msub><mi>Λ</mi><mi>B</mi></msub></mrow><mo>)</mo></mrow></mrow><mo>/</mo><mrow><mi>V</mi><mrow><mo>(</mo><mrow><msub><mi>Λ</mi><mi>R</mi></msub></mrow><mo>)</mo></mrow></mrow></mrow></mrow>
<mrow><msub><mi>r</mi><mrow><mi>G</mi><mi>R</mi></mrow></msub><mo>=</mo><mrow><mrow><mi>V</mi><mrow><mo>(</mo><mrow><msub><mi>Λ</mi><mi>G</mi></msub></mrow><mo>)</mo></mrow></mrow><mo>/</mo><mrow><mi>V</mi><mrow><mo>(</mo><mrow><msub><mi>Λ</mi><mi>R</mi></msub></mrow><mo>)</mo></mrow></mrow></mrow></mrow>
BR
GR
tHb
BR
GR
k
k
BR
GR
Из сигналов (10), можно получить только два нормированных сигнала, независящих от пространственного распределения освещенности ткани и положения ткани относительно блока регистрации изображения, -
и
. Для получения аналитического выражения, связывающего сигналы r и r с концентрацией гемоглобина в ткани F проведен расчет r и r в рамках вышеописанной модели слизистой оболочки при различных комбинациях модельных параметров. При каждой комбинации параметров по формулам (4)-(7) рассчитывались спектральные профили R(λ) и выполнялось их интегрирование с учетом аппаратурных функций P(λ) и S(λ) для получения RGB сигналов V(Λ) и нормированных сигналов r и r.
tHb
BR
GR
tHb
BR
GR
tHb
На фиг.7 представлен ансамбль рассчитанных сигналов и соответствующих им значений F. Видно, что представленные точки (r, r, F) расположены вблизи некоторой поверхности в трехмерном пространстве. Разброс точек относительной этой поверхности обусловлен общей вариативностью рассеивающих и поглощающих свойств соединительной ткани, диаметра капилляров с кровью и других параметров слизистой оболочки. Для аппроксимации поверхности, образованной точками (r, r, F), использовался полином вида
<mrow><msub><mi>F</mi><mrow><mi>t</mi><mi>H</mi><mi>b</mi></mrow></msub><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mn>0</mn><mo>≤</mo><mi>m</mi><mo>+</mo><mi>p</mi><mo>≤</mo><mn>3</mn></mrow></munder><mrow><msub><mi>a</mi><mrow><mi>m</mi><mi>p</mi></mrow></msub><msup><mrow><mrow><mo>[</mo><mrow><mn>1</mn><mi>n</mi><msub><mi>r</mi><mrow><mi>B</mi><mi>R</mi></mrow></msub></mrow><mo>]</mo></mrow></mrow><mi>m</mi></msup><msup><mrow><mrow><mo>[</mo><mrow><mn>1</mn><mi>n</mi><msub><mi>r</mi><mrow><mi>G</mi><mi>R</mi></mrow></msub></mrow><mo>]</mo></mrow></mrow><mi>p</mi></msup><mo>,</mo><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mo stretchy="false">(</mo><mn>11</mn><mo stretchy="false">)</mo></mrow></mstyle></mrow>
mp
tHb
BR
GR
tHb
k
Коэффициенты a определялись по методу наименьших квадратов. На фиг.8 сопоставлены точные и рассчитанные с использованием аппроксимирующей функции (11) значения концентрации гемоглобина. Средняя погрешность оценки значения F с использованием полинома (11) составляет 12.9%. Таким образом, рассматриваемые измерения сигналов r и r позволяют определять значение F с точностью, сравнимой с точностью определения данного параметра на основе вышерассмотренных оптимальных измерений r(λ) - Аналитическое выражение (11) может быть легко включено в программу компьютерной обработки RGB-сигналов (10), что позволит конвертировать цветное изображение исследуемого участка ткани в высококонтрастную картину распределения в нем гемоглобина.
Таким образом, заявляемое устройство позволяет в реальном масштабе времени получать двумерные распределения концентрации гемоглобина и степени оксигенации крови в слизистых оболочках по их спектральным изображениям. Данное устройство может с успехом использоваться в ходе проведения эндоскопических исследований слизистых оболочек полости рта, пищевода, органов желудочно-кишечного тракта и легких. Получаемые на его основе данные позволят от субъективного (визуального) исследования слизистой оболочки перейти к объективной (количественной) оценке ее состояния и в значительной мере повысить точность ранней диагностики онкологических заболеваний.
Литература.
1. Tajri, H. Proposal for consensus terminology in endoscopy: how should different endoscopic imaging techniques be grouped and defined? / H.Tajri, H.Niwa // Endoscopy.2008. - Vol.40, №9. - P.775-778.
2. 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.
3. Patents US 2009/0270702 A1, A61B 5/1455, A61B 6/00, Oct. 29, 2009.
4. Prahl, S.A. Optical absorption of hemoglobin [Электронный ресурс] / S.A.Prahl, - Режим доступа: http://omlc.ogi.edu/spectra/hemoglobin/index.html, свободный. - Название с экрана.
5. Барун, В.В. Оценка вклада локализованного поглощения света кровеносными сосудами в оптические свойства биологической ткани / В.В.Барун, А.П.Иванов // Оптика и спектр. - 2004. - Vol.96, №6. - P.1019-1024.
6. Bashkatov, A.N. Optical properties of human skin, subcutaneous and mucous tissues in the wavelength range from 400 to 2000 nm / A.N.Bashkatov, E.A.Genina, V.I.Kochubey, V.V.Tuchin // J. Phys. D: Appl. Phys. - 2005. - Vol.38. P.2543-2555.
7. Giraev, K.M. Optical absorption and scattering spectra of pathological stomach tissues / Giraev K.M., Ashurbekov N.A., Lakhina M.A. // J. Appl. Spectr. - 2011. - Vol.78, №1. P.95-102.
8. 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.
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. Egan W.G., Hilgerman T.W., Reichman J. Determination of absorption and scattering coefficients for nonhomogeneous media // Appl. Opt., 1973. Vol.12, №8. P.1816-1823.
11. Эсбенсен, К. Анализ многомерных данных, пер. с англ. под ред. О.Радионовой. - Черноголовка: ИПХФ РАН, 2005. С.50-98.
12. 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.