Способ визуализации областей объекта, содержащих микродвижения



Способ визуализации областей объекта, содержащих микродвижения
Способ визуализации областей объекта, содержащих микродвижения
Способ визуализации областей объекта, содержащих микродвижения
Способ визуализации областей объекта, содержащих микродвижения
Способ визуализации областей объекта, содержащих микродвижения
Способ визуализации областей объекта, содержащих микродвижения

Владельцы патента RU 2626310:

федеральное государственное бюджетное образовательное учреждение высшего образования "Нижегородская государственная медицинская академия" Министерства здравоохранения Российской Федерации (ФГБОУ ВО НижГМА Минздрава России) (RU)
Федеральное государственное бюджетное научное учреждение "Федеральный исследовательский центр Институт прикладной физики Российской академии наук (ИПФ РАН) (RU)

Изобретение относится к медицине, а именно к медицинской диагностике, и может быть использовано для получения и обработки изображений оптической интерферометрии. Осуществляют регистрацию набора интерферограмм, при этом сканирование зондирующим лучом осуществляют последовательно в главном и перпендикулярном к главному направлениях. Пересчитывают интерферограммы в наборе двумерных изображений. Вычисляют корреляционные ангиографические изображения. Проводят процедуры бинаризации с присвоением элементу изображения значения «1», если значение этого элемента превышает величину, установленную как уровень шума, и значения «0», если это условие не выполняется. Визуализируют полученный от перемножения бинарного массива данных и корреляционного ангиографического изображения результат в виде цветовой или яркостной амплитуды. Сканирование вдоль главного направления осуществляют с шагом dx, соотносящимся с диаметром зондирующего пучка D как dx<0.25 D. Вычисление корреляционного ангиографического изображении осуществляется следующим образом: двумерное структурное изображение ОКТ подвергают спектральному преобразованию вдоль координаты х, полученное распределение умножают на функцию окна, имеющую значение «0» в области низких частот и значение «1» в области частот, соответствующих скоростям движения, подлежащим визуализации. Полученное двумерное спектральное распределение подвергается процедуре обратного преобразования Фурье. Способ позволяет повысить точность и объективность визуализации областей, содержащих микродвижения, за счет корреляционного анализа элемента двумерного изображения в спектральной области с последующим восстановлением и цветовым кодированием информации изображения. 6 ил.

 

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

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

Известен способ выделения и визуализации информации о микродвижениях (патент US 8180134, В2 МПК7 G06K 9/00, опубл. 15.05.2012), при котором информация о движении выделяется за счет использования эффекта вырождения изображений оптической интерферометрии по знаку задержки, возникающего в условиях регистрации комплексной амплитуды зондирующей волны при наличии доплеровского сдвига частоты при рассеянии на движущихся рассеивателях.

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

Ближайшим аналогом разработанного способа является способ визуализации областей объекта, содержащих микродвижения, включающий регистрацию набора интерферограмм, при котором сканирование зондирующим лучом осуществляют последовательно в главном и перпендикулярном к главному направлениях, пересчет интерферограмм в набор двумерных изображений, вычисление корреляционного ангиографического изображения, проведение процедуры бинаризации с присвоением элементу изображения значения «1», если значение этого элемента превышает величину, установленную как уровень шума, и значения «0», если это условие не выполняется, и визуализацию полученного от перемножения бинарного массива данных корреляционного ангиографического изображения в виде цветовой или яркостной амплитуды (см. Joey Enfield, Enock Jonathan and Martin Leahy. «In vivo imaging of the microcirculation of the volar forearm using correlation mapping optical coherence tomography (cmOCT)» Joey Enfield, Enock Jonathan and Martin Leahy (Biomedical optics express, V. 2, I. 5, pp. 1184-1193)»).

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

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

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

Указанный технический результат достигается тем, что в известном способе визуализации областей объекта, содержащих микродвижения, включающем регистрацию набора интерферограмм, при этом сканирование зондирующим лучом осуществляют последовательно в главном и перпендикулярном к главному направлениях, пересчет интерферограмм в набор двумерных изображений, вычисление корреляционного ангиографического изображения, проведение процедуры бинаризации с присвоением элементу изображения значения «1», если значение этого элемента превышает величину, установленную как уровень шума, и значения «0», если это условие не выполняется; визуализацию полученного от перемножения бинарного массива данных корреляционного ангиографического изображения в виде цветовой или яркостной амплитуды, сканирование вдоль главного направления осуществляют с шагом dx, соотносящимся с диаметром зондирующего пучка D как dx<0.25 D, вычисление корреляционного ангиографического изображении осуществляется следующим образом: двумерное структурное изображение ОКТ подвергают спектральному преобразованию вдоль координаты х, полученное распределение умножают на функцию окна, имеющую значение «0» в области низких частот и значение «1» в области частот, соответствующих скоростям движения, подлежащим визуализации; полученное двумерное спектральное распределение подвергается процедуре обратного преобразования Фурье.

Предлагаемый способ поясняется графическим материалом.

На фиг. 1 представлено типичное изображение ОКТ для модельной среды, содержащей область, рассеиватели которой находятся в движении.

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

На фиг. 3 представлено схематичное изображение пространственной фильтрации изображения ОКТ по фиг. 1 в главном направлении. 1 - усредненный нормированный спектр стационарной области объекта, 2 - гауссовская аппроксимация нормированного спектра стационарной области объекта, 3 - уровень отсечки, 4 - форма режекционного фильтра.

На фиг. 4 представлен результат пространственной фильтрации изображения ОКТ по фиг. 1 в главном направлении.

На фиг. 5 представлен пример визуализации двумерного изображения ОКТ по фиг. 1 с окрашиванием области наличия движения по фиг. 4.

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

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

Предлагаемый способ осуществляют следующим образом.

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

Сканирование вдоль главного направления осуществляют с шагом dx, соотносящимся с диаметром зондирующего пучка D как dx<0.25 D, т.е. сканирование осуществляют с шагом, много меньшим диаметра зондирующего пучка, что обеспечивает неизотропность спекл-структуры, составляющей изображение (фиг. 1) - элементы становятся вытянутыми вдоль главного направления. При этом вследствие того что отдельные строки регистрируются в различные моменты времени, наличие динамики в отдельных областях объекта (выделено пунктиром на фиг. 1) приводит к наблюдаемому изменению поперечного масштаба отдельных элементов изображения.

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

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

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

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

Пример конкретной реализации способа визуализации областей объекта, содержащих микродвижения

Для конкретной реализации способа визуализации областей объекта, содержащих микродвижения, был выбран объект - ухо мыши линии BalbC.

Произвели регистрацию изображения оптической интерферометрии, для чего направили зондирующее излучение на объект исследования, осуществили последовательную регистрацию строковых элементов изображения оптической интерферометрии при изменении положения зондирующего луча на поверхности объекта исследования вдоль главного и перпендикулярного ему. Сканирование вдоль главного направления осуществляли с шагом dx, много меньшим диаметра зондирующего пучка D (dx<0.15 D). Полученный набор интерферограмм подвергли преобразованию Фурье в направлении оси зондирующего пучка, в результате чего получили двумерные изображения исследуемого объекта. Полученные изображения подвергли преобразованию Фурье в главном направлении. Полученные изображения были подвергнуты процедуре фильтрации посредством перемножения полученного спектрального распределения и оконной функции режекционного фильтра. Полученные двумерные спектральные распределения были подвергнуты процедуре обратного преобразования Фурье в главном направлении, в результате чего на полученном ангиографическом изображении сохранились только области, в которых было отмечено наличие микродвижений (фиг. 6.в). Полученные области соответствуют положению сосудов в исследуемом объеме. После бинаризации полученного ангиографического изображения и цветового кодирования элементов структурного изображения было произведено сопоставление взаимного расположения сосудов и структуры исследуемого объекта с его прямым фотографическим изображением (фиг. 6.а) и изображением, полученным с использованием флуоресцентной микроскопии (стериомикроскоп Axio Zoom.V16 (Zeiss, Germany) с флуорофором FITC, конъюгированном с декстраном 150 kD (Sigma, USA). Сравнение показало, что структурные элементы, представленные на изображении, полученном с использованием способа визуализации областей объекта, содержащих микродвижения, соответствуют сетке кровеносных сосудов в объекте исследования.

Способ визуализации областей объекта, содержащих микродвижения, включающий регистрацию набора интерферограмм, при этом сканирование зондирующим лучом осуществляют последовательно в главном и перпендикулярном к главному направлениях, пересчет интерферограмм в набор двумерных изображений, вычисление корреляционного ангиографического изображения, проведение процедуры бинаризации с присвоением элементу изображения значения «1», если значение этого элемента превышает величину, установленную как уровень шума, и значения «0», если это условие не выполняется, визуализацию полученного от перемножения бинарного массива данных корреляционного ангиографического изображения в виде цветовой или яркостной амплитуды, отличающийся тем, что сканирование вдоль главного направления осуществляют с шагом dx, соотносящимся с диаметром зондирующего пучка D как dx<0.25D, вычисление корреляционного ангиографического изображении осуществляется следующим образом: двумерное структурное изображение ОКТ подвергают спектральному преобразованию вдоль координаты х, полученное распределение умножают на функцию окна, имеющую значение «0» в области низких частот и значение «1» в области частот, соответствующих скоростям движения, подлежащим визуализации, полученное двумерное спектральное распределение подвергается процедуре обратного преобразования Фурье.



 

Похожие патенты:

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

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

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

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

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

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

Изобретение относится к области обнаружения движущегося объекта. Технический результат – повышение точности обнаружения движущегося объекта или его возможного кандидата.

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

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

Изобретение относится к технологиям обнаружения движущихся объектов в последовательности видеоизображений. Техническим результатом является повышение эффективности обнаружения движущихся объектов в последовательности видеоизображения, за счет анализа качества изменения изображения.
Изобретение относится к офтальмологии. Проводят оптическую когерентную томографию (ОКТ) с измерением толщины перипапиллярного слоя нервных волокон сетчатки (СНВС) каждого глаза по четырем сегментам - височному, верхнему, носовому и нижнему.
Изобретение относится к медицине, рентгенологии, методам лечения больных туберкулезом легких без бактериовыделения. По данным рентгеновской компьютерной томографии органов грудной клетки (РКТ ОГК) определяют наиболее крупный фокус туберкулезного инфильтрата – основной фокус.
Изобретение относится к медицине, офтальмологии, может быть использовано для оценки эффективности аппаратного лечения миопии высокой степени. При помощи оптической когерентной томографии (ОКТ)-ангиографии определяют микроциркуляцию макулярной области.
Изобретение относится к медицине, онкологии и хирургии и может найти применение при лечении резектабельных злокачественных нейроэндокринных опухолей (НЭО) головки поджелудочной железы (ПЖ).
Изобретение относится к медицине, онкологии и хирургии и может найти применение при лечении злокачественных нейроэндокринных опухолей (НЭО) головки поджелудочной железы (ПЖ).
Изобретение относится к медицине, онкологии и хирургии и может найти применение при лечении резектабельных первичных местнораспространенных злокачественных нейроэндокринных опухолей (НЭО) головки поджелудочной железы (ПЖ).

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

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

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

Изобретение относится к медицине клинической кардиологии и диагностике, может быть использовано для количественной оценки нарушений перфузии миокарда правого желудочка (ПЖ) по данным однофотонно-эмиссионной компьютерной томографии (ОФЭКТ). Проводят однофотонно-эмиссионную компьютерную томографию (ОФЭКТ) и измеряют соотношения интенсивности накопления радиофармпрепарата (РФП) правым и левым желудочками (ЛЖ). При этом производят обведение ПЖ полуавтоматическим способом с помощью режимов "ограничения" и "маскирования" с получением количественных полярных карт перфузии ПЖ в дополнение к аналогичным картам ЛЖ. Определяют количественные параметры нарушений перфузии ПЖ аналогично ЛЖ, включающие данные об абсолютной статистике счета в каждом сегменте ПЖ. Далее вычисляют коэффициент RV/LV(s), используя данные количественной полярной карты:RV/LV(s)=CmaxRV/CmaxLV,где Cmax RV, LV - максимальное значение числа импульсов в сегментах, относящихся к ПЖ и ЛЖ соответственно. При значениях RV/LV(s)<0.3 состояние ПЖ трактуют как нормальное, значения 0.3<RV/LV(s)<0.5 трактуют как пограничные, значения RV/LV(s)>0.5 трактуют как признаки перегрузки и гипертрофии ПЖ. Способ обеспечивает повышение чувствительности, точности количественной оценки и диагностики нарушений перфузии миокарда ПЖ, характеризуется легкой воспроизводимостью и операторонезависимостью. 4 ил., 3 пр.

Изобретение относится к медицине, а именно к медицинской диагностике, и может быть использовано для получения и обработки изображений оптической интерферометрии. Осуществляют регистрацию набора интерферограмм, при этом сканирование зондирующим лучом осуществляют последовательно в главном и перпендикулярном к главному направлениях. Пересчитывают интерферограммы в наборе двумерных изображений. Вычисляют корреляционные ангиографические изображения. Проводят процедуры бинаризации с присвоением элементу изображения значения «1», если значение этого элемента превышает величину, установленную как уровень шума, и значения «0», если это условие не выполняется. Визуализируют полученный от перемножения бинарного массива данных и корреляционного ангиографического изображения результат в виде цветовой или яркостной амплитуды. Сканирование вдоль главного направления осуществляют с шагом dx, соотносящимся с диаметром зондирующего пучка D как dx<0.25 D. Вычисление корреляционного ангиографического изображении осуществляется следующим образом: двумерное структурное изображение ОКТ подвергают спектральному преобразованию вдоль координаты х, полученное распределение умножают на функцию окна, имеющую значение «0» в области низких частот и значение «1» в области частот, соответствующих скоростям движения, подлежащим визуализации. Полученное двумерное спектральное распределение подвергается процедуре обратного преобразования Фурье. Способ позволяет повысить точность и объективность визуализации областей, содержащих микродвижения, за счет корреляционного анализа элемента двумерного изображения в спектральной области с последующим восстановлением и цветовым кодированием информации изображения. 6 ил.

Наверх