27 декабря 2012 г.
Image Echo частиц велосиметрии (EPIV) система, способная приобретения двумерного поля скорости в оптически непрозрачных жидкостей или через оптически непрозрачных геометрии описаны и проверки измерений в трубе не сообщается.
Цель данного эксперимента — получение мгновенных двумерных полей скорости в потоке Хагена-Пуазейля, также известном как ламинарный поток в трубе. С помощью метода эхо-particle image velocimetry (E-PIV) продемонстрированы измерения в рециркуляционном потоке в трубе, содержащей смесь воды и глицерина в соотношении 50:50. В используемой системе трубного потока для создания потока применяется постоянный напор, поддерживаемый аквариумным насосом.
Фазированный ультразвуковой датчик закрепляется на стенке трубы системы потока, после чего осуществляется потоковая передача ультразвуковых изображений в B-режиме. Жидкость насыщается полыми стеклянными сферами или трассирующими частицами, которые точно повторяют движения потока. Затем ультразвуковые изображения регистрируются, передаются на ПК и конвертируются в формат изображений, совместимый с коммерческим программным обеспечением для PIV-велосиметрии.
Алгоритмы кросс-корреляции применяются к последовательным изображениям ультразвукового B-режима для вычисления двумерных полей скорости. В конечном итоге векторные поля анализируются для расчета интересующих показателей потока, таких как среднегрупповое касательное напряжение и завихренность. Хотя этот метод позволяет получить представление об основах гидродинамики, он может использоваться и часто применяется в практических системах потоков, включая биомедицинские приложения, например, при изучении артериальных или межжелудочковых потоков.
Нас заинтересовала данная методика в связи с проводимыми нами экспериментами с использованием сжиженной биомассы и жидкостей. Для настройки системы EP IV сначала включите насосы. Это запустит рециркуляцию потока в трубках с постоянной скоростью.
Затем нанесите на ультразвуковой датчик топический гель на водной основе. Гель минимизирует потери при передаче ультразвукового луча между поверхностью датчика и трубой. После этого, используя специально сконструированный держатель датчика со сквозными трубными фитингами, закрепите ультразвуковой датчик на внешней стенке трубы и включите ультразвуковой аппарат.
После загрузки всех систем на экране ультразвукового сканера появится видеопоток изображений. Для проведения измерений EPIV с помощью линейного датчика по умолчанию используется 2D-режим. Сначала отвесьте соответствующий сухой вес полых стеклянных сфер номинальным диаметром 10 мкм, чтобы их конечная концентрация при добавлении в систему потока составила около 17 частей на миллион по весу.
Затем отберите определенный объем жидкости из резервуара и добавьте в нее частицы для приготовления концентрированного раствора частиц. При введении в систему потока эти частицы будут служить ультразвуковыми контрастными веществами или трассерными частицами. Добавьте концентрированный раствор частиц в систему рециркуляционного трубного потока, перемешав его в резервуарах с водой.
Затем стеклянные сферы можно наблюдать на мониторе ультразвукового аппарата. Через несколько минут стеклянные сферы равномерно распределятся по всей системе. Одной из наиболее сложных частей данной процедуры является получение четких ультразвуковых изображений с высоким разрешением. Чтобы максимально повысить качество этих изображений, мы настраиваем усиление, фокальные точки и динамический диапазон, основываясь на априорных оценках скорости потока.
Далее мы оптимизируем эти параметры в режиме реального времени, анализируя ультразвуковые изображения. Используйте ручку регулировки глубины на панели управления ультразвукового аппарата, чтобы установить глубину изображения на уровне трех сантиметров. Затем с помощью ручки усиления 2D отрегулируйте общее усиление изображения, чтобы увеличить его яркость так, чтобы частицы-затравки были четко видны на панели прибора.
Отрегулируйте ползунки временной компенсации усиления, чтобы ослабить рассеяние от стенок трубки и компенсировать затухание ультразвукового сигнала, связанное с глубиной. Это позволит устранить избыточные помехи в верхней и нижней частях стенок трубки в 2D-режиме. Ручки управления в верхней части панели управления (слева направо) отвечают за ширину, частоту фокусировки и частоту кадров.
Используйте эти регуляторы для дальнейшей настройки изображения, чтобы достичь максимально возможного физического разрешения, частоты и частоты кадров для проведения анализа. Затем установите рабочую частоту датчика на 10 мегагерц и частоту кадров на 49,5 кадра в секунду. Обратите внимание, что эти четыре параметра взаимозависимы.
Следовательно, при сканировании ультразвуковым изображением существует компромисс между пространственным и временным разрешением. Из-за ограниченного латерального разрешения стеклянные сферы будут размыты в латеральном направлении и будут выглядеть на изображении как эллипсоиды. После оптимизации параметров наступает время сбора данных на панели управления ультразвуковым прибором.
Нажмите кнопку «Новое исследование», чтобы начать новый эксперимент. В поле «Пациент» введите pipe flow в качестве фамилии и дату в качестве имени.
Введите номер теста в поле идентификатора пациента (patient ID). Ультразвуковое сканирование начнется автоматически после достижения предустановленного максимума от 1000 до 1 500 изображений. По мере продолжения сканирования начинается новый цикл сканирования.
Отрегулируйте параметры визуализации до тех пор, пока частица-трассер не окажется в четком фокусе, при этом в каждой области исследования должно находиться около 10 частиц. Чтобы перезапустить цикл записи сканирования, нажмите кнопку freeze на панели управления ультразвуковым сканером. После получения достаточного количества подходящих изображений нажмите кнопку freeze.
Затем нажмите кнопку C-петли на панели управления ультразвуковым аппаратом. Выберите все изображения, чтобы включить все ультразвуковые снимки в набор для анализа. После выбора необходимых изображений нажмите кнопку сохранения изображений, чтобы сохранить выбранный набор ультразвуковых снимков.
После сохранения изображений нажмите кнопку архивации на панели управления ультразвуковым сканером. Когда появится запрос, в маленьком окне выберите интересующий вас syn-цикл для сохранения на локальный жесткий диск. Затем с помощью курсора мыши выберите команду «завершить исследование».
Нажмите кнопку архивации и с помощью курсора мыши сначала выберите пункт «more», а затем «disc management». Это перенесет сохраненный цикл или циклы cyl на ПК с установленным программным обеспечением для измерения скорости частиц по изображениям (particle image velocity symmetry или PIV). После захвата и сохранения изображений ультразвуковой снимок необходимо преобразовать из файла стандарта DICOM (digital imaging communications in medicine) в файл изображения формата JPEG (joint photographic experts group) для последующего анализа.
Используйте скрипт MATLAB, запускающий файл DICOM to jpeg.m, для конвертации файлов DICOM в формат JPEG. Данный скрипт был разработан внутри организации; его можно получить в образовательных целях по указанному здесь веб-адресу. После конвертации файла откройте программное обеспечение Davi от Law Vision.
Дважды щелкните по иконке davi, выберите «new project» (новый проект), затем выберите PIV. На панели инструментов выберите «import images» (импорт изображений) и вариант «import via numbered files» (импорт через нумерованные файлы). Затем в выпадающем меню найдите папку, в которой хранятся ультразвуковые изображения в формате JPEG, и дважды щелкните по первому изображению из набора.
Это позволит импортировать все ультразвуковые изображения из данного пронумерованного набора для определения области интереса (ROI) для анализа, включающей только жидкость. Чтобы создать маску, примените инструмент маскирования, введя координаты прямоугольной области с помощью двух точек координат x и y на основе информации из DICOM-файла и данных о размере пикселей. Затем на главной панели управления в dvu перейдите на вкладку под текущим проектом, содержащую импортированные изображения.
Выберите таблицу с надписью batch processing. Это откроет окно обработки векторов в программе Davi для пакетной обработки из списка операций; используя дерево PIV plus PIV time series, выберите расчет параметров векторов и установите режим multipass с уменьшением размера области допроса с 64 pixels by 64 pixels до 12 pixels by 12 pixels с перекрытием 50%. Установите ограничение относительного диапазона векторов (relative vector range restriction) на all, а ограничение абсолютного диапазона векторов (absolute vector range restriction) — на пять pixels. Затем примените медианный фильтр для подавления шума и сглаживания векторных полей.
Затем для обработки векторов в меню параметров расчета векторов установите флажок, указывающий, что диапазон данных соответствует маскированной области (data range equals use masked area). Обратите внимание, что оптимальный подбор параметров расчета векторов зависит от геометрии потока, свойств потока, разрешения изображения, плотности частиц-трассеров и требуемого количественного анализа потока. После установки всех необходимых параметров в левой части экрана пакетной обработки выберите общее количество изображений для обработки.
Нажмите «начать обработку» (start processing). Это позволит вычислить поле смещения между последовательными ультразвуковыми изображениями с помощью алгоритмов кросс-корреляции. Для анализа обработанных данных экспортируйте векторные поля UCV из DAVO в виде txt-файлов.
Для этого на экране проекта в ветке JPEG-изображений выберите ветку векторного смещения. На панели инструментов перейдите на вкладку экспорта. Выберите тип файла.
Запросите текстовый файл. Выберите создание папки для экспорта и нажмите «экспорт». Затем откройте файл в MATLAB, запустив соответствующий скрипт MATLAB.
Экспортированные векторные поля именуются B-X-X-X-X-X.txt, где X представляет собой число, увеличивающееся от 1 до 99 999. Каждый файл содержит четыре столбца данных, которые можно просмотреть с помощью текстового редактора: один — координата X вектора на изображении, два — координата Y вектора на изображении, три — X-компонента смещения (потоковое смещение) и четыре — Y-компонента смещения, описывающая смещение по нормали к стенке. Для вычисления векторного поля скоростей U как функции от X и Y, где X и Y соответствуют пространственным координатам на ультразвуковом изображении, сначала преобразуют поле смещений D от X и Y, измеренное в пикселях, в поле смещений, измеренное в метрах, используя параметр масштабирования изображения M, указанный в единицах m/pixel.
Затем поле смещения делится на скорректированный с учетом развертки временной интервал между изображениями delta T, где delta T равен величине, обратной частоте кадров (указанной в кадрах в секунду), плюс поле смещения, деленное на время, за которое ультразвуковое изображение проходит всю ширину кадра. Таким образом, U от X и Y равно M умножить на D от X и Y, деленное на delta T. В формате DICOM изначально хранится структура файла, содержащая информацию, необходимую для вычисления параметра масштабирования изображения M и скорректированного с учетом развертки временного интервала Delta T. В настоящем исследовании M составляет 77 microns per pixel, FPS равен 49.5, а B равен 25, 047 pixels per second.
Наконец, вычислите усредненные по ансамблю векторные поля скорости, а также нормальные профили средней скорости и любые другие интересующие характеристики потока. Чтобы продемонстрировать эффективность EPIV и оценить неопределенность измерений, были получены двумерные мгновенные поля скорости ламинарного потока в трубе, как описано в данном видео; на этом мгновенном векторном графике векторы скорости показаны для каждой четвертой колонки, а фоновая цветовая контурная карта соответствует модулю скорости. Двумерное пространственное положение вектора скорости обозначено как D/D и X/D, где D — радиальное положение, измеренное от верхней стенки.
D — диаметр трубы, а X — положение по направлению потока, измеренное от входа в трубу. Видимая параболическая форма векторов скорости вдоль столбцов указывает на то, что измерения соответствуют ожидаемому профилю скорости для трубного потока. График усредненных по ансамблю векторов, рассчитанный путем усреднения 1000 мгновенных векторных графиков, представляет собой изображение среднего поля скоростей.
Это также позволит нивелировать случайные шумовые погрешности в мгновенных векторных полях. Векторы скорости направлены преимущественно вдоль потока. Максимальные значения скорости наблюдаются на осевой линии трубы.
Скорости снижались до нуля у стенок трубы, и поток был примерно симметричным; здесь представлен средний профиль продольной скорости вдоль радиуса трубы, полученный путем усреднения ансамблевого векторного графика по строкам в горизонтальном направлении. Также показан ожидаемый средний профиль скорости для ламинарного потока в трубе с учетом экспериментальных условий.
Согласованность между измерениями EPIV и ожидаемым профилем Хиггина-Позо наилучшая вблизи осевой линии трубы и наихудшая — у стенок трубы. Значительные расхождения у стенки, вероятно, обусловлены сильным отражением и преломлением ультразвуковых волн на изогнутой поверхности стенки трубы, что создает высокую интенсивность изображения в этих областях; высокая интенсивность у стенки маскирует интенсивность частиц, что приводит к ошибке измерения. Благодаря разработке этого метода исследователи, изучающие гидродинамику в технических или биологических потоковых системах, теперь могут получать пространственно-временные вариации поля скоростей в оптически непрозрачных жидкостях или через оптически непрозрачные геометрии.
После просмотра этого видео вы должны получить четкое представление о принципах работы EPIV, его ограничениях, а также о том, как собрать и эксплуатировать систему EPIV с использованием коммерческого аппарата ультразвукового исследования.
Просмотрите полный транскрипт и получите доступ к тысячам научных видео
В данной статье описывается система эхо-Particle Image Velocimetry (EPIV), предназначенная для регистрации двумерных полей скоростей в оптически непрозрачных жидкостях. Представлены результаты валидационных измерений в ламинарном потоке в трубе, демонстрирующие возможности данной системы.
Эхо-Particle Image Velocimetry (EPIV) позволяет проводить неинвазивные измерения полей скоростей жидкости с высоким разрешением в оптически непрозрачных системах, что способствует механистическому пониманию явлений переноса в сложных биологических и инженерных потоках. Эта возможность помогает снизить риски, связанные с допущениями в гидродинамике на ранних этапах валидации мишеней и разработки доклинических моделей, предоставляя количественные пространственно-разрешенные данные о переносе массы, импульса и энергии. Поля скоростей, полученные методом EPIV, повышают прогностическую достоверность в доклинических системах с ограниченным прямым оптическим доступом, таких как вязкие биоматериалы или фантомы, имитирующие ткани.
EPIV интегрируется в рабочий процесс поиска новых препаратов в качестве инструмента для оценки гидродинамики на основе выдвинутых гипотез и располагается между этапами раннего скрининга формуляций и доклинической функциональной валидации, особенно в тех случаях, когда оптическая непрозрачность ограничивает возможности традиционных методов визуализации.