Доплеровская оптическая когерентная томография (ДОКТ) – методика получения информации о потоке жидкости, которая сочетает принципы Доплеровской лазерной анемометрии и оптической когерентной томографии. ДОКТ является эффективным не инвазивным методом получения информации о потоке крови в сердечнососудистой системе человека. Данная методика основана на волновых свойствах света, который по мере прохождения в тканях отражается от частиц движущегося потока [4,8].
Для проведения ДОКТ исследований в простейшем случае используется интерферометр Михельсона, в котором электромагнитное излучение разбивается на две части, одна из которых идет в опорное плечо, где установлено зеркало, двигающееся с постоянной скоростью, а вторая направляется в объектное плечо, где располагается исследуемый объект. После отражения от зеркала в опорном плече и от различных оптических неоднородностей в плече образца световые пучки снова соединяются вместе, что приводит к возрастанию интерференции, в том случае, если разница длин оптических путей пучков меньше чем длина когерентности [1,4,5]. Если же оптическая неоднородность в плече образца также двигается с некоторой постоянной скоростью, то интерференционный сигнал претерпевает доплеровский сдвиг, 
, который пропорционален искомой скорости потока 
:
,
где 
– длина волны источника излучения, а 
– угол между направлением освещения и вектором скорости (Доплеровский угол) [7,9].
Классическая методика ДОКТ ограничена тем, что для определения значения скорости необходимо знать точное значения доплеровского угла, что довольно затруднительно при исследовании сильно рассеивающих сред, таких как кожные покровы или сетчатка глаза человека, изучение которых является одним из первостепенных приложений данной методики. Даже в случае простейших измерений скорости потока in vitro, неправильная оценка доплеровского угла может привести к 50 % ошибке полученных значений. Для решения данной проблемы ранее уже были предложены различные подходы, один из которых предполагает использование призмы Волластона в плече образца интерферометра для точного измерения доплеровского угла. Применение данного метода оказалось довольно затруднительным, так как равенство оптических путей в каждом канале на практике является зачастую недостижимым. Также для оценки поперечной скорости без предварительного определения доплеровского угла предлагалось использовать уширение доплеровского спектра, возникающее при пересечении частицами потока светового пучка, что, впрочем, не позволяет оценить действительную скорость потока, так как продольная составляющая скорости при этом не учитывается. Соответственно разработка методики получения скорости потока, не требующей знания доплеровского угла до проведения исследования, или позволяющей точно оценить его, является одним из важнейших шагов, необходимых для дальнейшего развития ДОКТ [10,11,13].
В представленной работе описаны методы, которые позволяют определить значение скорости потока жидкости, минуя указанные ограничения, накладываемые стандартными принципами ДОКТ. При этом для вычисления доплеровского угла используется среднеквадратическое отклонение 
уширения доплеровского спектра, что в совокупности с представленным способом фильтрации позволяет повысить точность исследования и преодолеть некоторые специфические ошибки, которые могут возникать при использовании описываемой методики.
Материалы и методы исследования
Двумерный вектор скорости может быть получен с помощью двух его компонентов, которые представляют собой проекции вектора полной скорости 
на направление двух пересекающихся осей (рис. 1). Измерение данных компонентов 
и 
получают путем двух последовательных измерений в положениях № 1 и № 2, отстоящих друг от друга на расстоянии 
. При известном расстоянии между положениями вектор скорости может быть представлен как:
Рис. 1. Схематичное изображение экспериментальной модели исследования потока
где 
– представляет разность углов 
и 
, 
а выражаются как:
Если расстояние недостаточно велико, ошибка измерения величины и направления потока может быть слишком большой, так как значение угла может быть близко к нулю и направление вектора скорости с помощью представленных выражений может быть рассчитано неверно. В то же время, это расстояние должно быть настолько мало, насколько это возможно, так как иначе данный подход может оказаться неприменим для исследования кровотока в сердечнососудистой системе человека [2,6,12].
Анализ эффективности применяемой методики проводился посредством соответствующего эксперимента, при котором производилось измерение скорости течения жидкости в прозрачном сосуде. В качестве источника излучения использовался суперлюминесцентный диод, длина волны излучения которого составляет λ=1300 нм, мощность 0.5 мВт, ширина оптического спектра Δλ=50 нм. Скорость сканирования в опорном плече составляет 2 мм/с, частота дискретизации АЦП – 50 кГц. В качестве потока использовалась жидкость с оптическими характеристиками, эквивалентными оптическим характеристикам крови, пропущенная через цилиндрический капилляр диаметром 1 мм. Скорость потока регулировалась насосом в пределах 15–40 мм/с, со значением доплеровского угла 84
. Диаметр фокусного пятна используемой фокусирующей линзы составляет 40 мкм, конфокальный параметр 1.5 мм. Для приближения экспериментальных условий к реальным, использован цифровой гауссовский шум, значение которого составляет 5–10 % от максимального значения полученной амплитуды.
Результаты и обсуждение
Спектрограмма полученного сигнала представлена на верхней части рис. 2а. Можно отметить появление артефактов и шумов различной природы, в частности возникновение отраженной составляющей, которая имеет такую же амплитуду и как и измеренный сигнал, а также аддитивного шума на частотах 2.5–4.5 кГц, что ведет к некорректному определению центроида (средневзвешенной скорости) доплеровских спектров и (рис. 2б), поэтому дополнительная компьютерная обработка результатов является необходимой.
а)
б)
Рис. 2. Спектрограммы полученного сигнала при доплеровском угле 84 (сверху) и при перпендикулярном сканировании (снизу) (а); центроиды доплеровских спектров (синий) и (красный) (б)
Одним из важных этапов компьютерной обработки является усреднение по соседним А-сканам. После усреднения значение отношения С/Ш возрастает в несколько раз пропорционально числу усредненных А-сканов. Данная процедура позволяет избавиться от дополнительно смоделированного гауссовского шума. При увеличении числа усредненных А-сканов соотношение С/Ш будет линейно возрастать до определенного предела, но при исследованиях in vivo усреднение целесообразно проводить только по небольшому числу 
соседних А-сканов, так как информация о скорости потока в небольших капиллярах при усреднении будет потеряна, если велико [3]. Также для устранения шумов использовался регрессионный анализ и частотная фильтрация (спектральное вычитание). Полученные после обработки центроиды и представлены на рис. 3.
Рис. 3. Центроиды доплеровских спектров и после компьютерной обработки
Программная обработка сигнала позволяет удалить большую часть аддитивного шума на частотах 2.5–4.5 кГц, не повреждая при этом информации об искомом сигнале, что позволяет более точно локализовать местонахождения начала и конца параболического профиля.
Определение доплеровского угла производится путем построения зависимости от частоты с помощью регрессионного анализа (рис. 4).
а)
б)
Рис. 4. Зависимости от частоты до (а) и после фильтрации помех (б).
Регрессионный анализ (красная линия) даёт хорошее соответствие, R2~0.8-0.9
Важной особенность представленной методики является то, что значение скорости потока может быть получено даже в том случае, если направление зондирующего пучка перпендикулярно направлению потока, что в классической ДОКТ технике сделать невозможно. В этом случае значение продольной скорости будет равно нулю, а уширение доплеровского спектра будет возникать по причине ограниченного времени пролёта и флуктуаций интереференционной картины, вызываемых частицами, движущимися через зондирующий пучок. Спектрограмма сигнала, получаемого при перпендикулярном сканировании, представлена в нижней части рисунка 2а, а доплеровские центроиды и на рис. 5.
Рис. 5. Центроиды доплеровских спектров и при перпендикулярном сканировании
Измеренные значения скорости и доплеровского угла соответствуют предустановленным с погрешностью 5–10 %.
Заключение
Результат данного исследования показывает, что представленная методика доплеровской оптической когерентной томографии, объединённая с цифровой обработкой сигнала спектральным вычитанием, усреднением и медианной фильтрацией позволяет с большой точностью определить значение скорости потока и доплеровского угла без использования призмы Волластона или иных методов вычисления скорости, требующих предварительного знания величины доплеровского угла. Следующим логичным шагом исследования является применение данной методики для исследования потоков in vivo, что позволит оценить возможность коммерческого применения представленных методов.
Рецензенты:
Фролов С.В., д.т.н., зав. кафедрой «Биомедицинская техника», ФГБОУ ВПО «Тамбовский государственный технический университет», г. Тамбов;
Туголуков Е.Н., д.т.н., профессор кафедры «Техника и технологии производства нанопродуктов», ФГБОУ ВПО «Тамбовский государственный технический университет», г. Тамбов.