Восстановление распределения вектора скорости кровотока в линейном и нелинейном акустических томографах тема автореферата и диссертации по физике, 01.04.06 ВАК РФ

Матвеев, Олег Владимирович АВТОР
кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Москва МЕСТО ЗАЩИТЫ
2014 ГОД ЗАЩИТЫ
   
01.04.06 КОД ВАК РФ
Диссертация по физике на тему «Восстановление распределения вектора скорости кровотока в линейном и нелинейном акустических томографах»
 
Автореферат диссертации на тему "Восстановление распределения вектора скорости кровотока в линейном и нелинейном акустических томографах"

На правах рукописи УДК 534.2 : 534.7 : 534.222

Матвеев Олег Владимирович

ВОССТАНОВЛЕНИЕ РАСПРЕДЕЛЕНИЯ ВЕКТОРА

СКОРОСТИ КРОВОТОКА В ЛИНЕЙНОМ И НЕЛИНЕЙНОМ АКУСТИЧЕСКИХ ТОМОГРАФАХ

Специальность: 01.04.06 - акустика

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

Москва-2014

005549802

005549802

Работа выполнена на кафедре акустики физического факультета Московского государственного университета им. М.В. Ломоносова.

Научные руководители: доктор физико-математических наук,

профессор В .А. БУРОВ, МГУ имени М.В.Ломоносова, физический факультет, г. Москва

кандидат физико-математических наук РУМЯНЦЕВА Ольга Дмитриевна, МГУ имени М.В.Ломоносова, физический факультет, г. Москва

Официальные оппоненты: доктор физико-математических наук,

профессор ПИРОГОВ Юрий Андреевич, МГУ имени М.В.Ломоносова, физический факультет, г. Москва

кандидат физико-математических наук, начальник отдела гидроакустических систем ПРУДНИКОВА Инесса Петровна, ОАО "Тетис Про", г. Москва

Ведущая организация: Нижегородский государственный

университет имени Н.И.Лобачевского, г. Нижний Новгород

Защита диссертации состоится *' 26 "июня 2014 г. в 11-00 часов на заседании Специализированного Совета Д.501.001.67 в МГУ им. М.В. Ломоносова по адресу: 119991, г.Москва, ГСП-1, Ленинские Горы, МГУ, физический факультет ___

С диссертацией можно ознакомиться в научной библиотеке физического факультета МГУ им М.В. Ломоносова.

Автореферат разослан " Иг- мая 2014 г.

Ученый секретарь Специализированного Совета Д.501.001.67; кандидат физико-математических на

А.Ф. КОРОЛЕВ

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы

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

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

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

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

Для достижения цели работы были поставлены и решены следующие задачи:

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

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

3. Оценка помехоустойчивости алгоритмов и границ их работоспособности.

4. Проведение модельных численных расчетов, имитирующих физический эксперимент.

5. Проведение физического эксперимента.

Научная новизна работы

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

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

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

Достоверность

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

Научная и практическая значимость работы

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

сочетании с оценками пространственных распределений линейных характеристик и нелинейных параметров биологической ткани.

Основные положения, выносимые на защиту

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

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

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

Апробация работы

Материалы докладывались на следующих конференциях:

• Международная конференция студентов, аспирантов и молодых ученых по фундаментальным наукам "Ломоносов-2003". Секция "Физика".

• Международная конференция студентов, аспирантов и молодых ученых по фундаментальным наукам "Ломоносов-2004". Секция "Физика".

• II Евразийский конгресс по медицинской физике и инженерии "Медицинская физика-2005" 21-24 июня 2005.

• X Всероссийская школа-семинар "Волновые явления в неоднородных средах". Звенигород, Московская область, 22-27 мая 2006.

• The 8th Pacific Symposium on Flow Visualization and Image Processing. Ed. I.A.Znamenskaya. Moscow: Moscow State University, 2011 (8-й мировой симпозиум по визуализации потоков и вычислению изображений им. И.А. Знаменской, МГУ, 2011).

Публикации

Основные результаты диссертации изложены в 12 работах [А1-А12], в том числе, четыре статьи опубликованы в журналах, включенных в рекомендованный список ВАК РФ [Al, А2, A3, А4]. Список публикаций приводится в конце автореферата.

Структура и объем диссертации

Диссертация состоит из введения, двух глав, заключения и списка литературы. Список литературы состоит из 102 наименований. Общий объем работы составляет 165 страниц, включая 32 рисунка.

Личный вклад автора.

Все результаты получены автором лично на основе модельных данных, полученных по методам, разработанным автором.

СОДЕРЖАНИЕ РАБОТЫ

В первой главе в разделе 1.1 приведено обоснование актуальности выбранной темы диссертационной работы (далее - работы), показана цель работы и поставлены задачи, решенные по ходу работы. Приведены формулировки научной новизны и значимости работы. Показан вклад автора, и дана формулировка, подтверждающая достоверность результатов работы. Раскрыты выносимые на защиту положения, перечислены конференции и публикации, где докладывались и изложены результаты работы. В разделе 1.2 представлен обзор научных публикаций по теме работы. Описано развитие современных методов УЗ диагностики от одноканальных систем до многоканальных систем, как линейных, так и нелинейных. Показаны связь и различие доплеровских средств диагностики и корреляционных методов.

Большинство появившихся в последние десятилетия систем медицинской акустоскопии - ультразвуковые многоканальные системы активной локации (УЗИ) и разрабатывавшиеся в последние годы многоракурсные томографические системы [1]. Эти диагностические комплексы представляют собой, как правило, линейные когерентные системы активной локации, основанные на аддитивном (когерентном) сложении сигналов, принимаемых элементами антенной решетки. Важным достоинством этих систем является линейная зависимость выходного отношения сигнал/помеха от такого же входного отношения. В когерентных системах, применяемых для диагностики состояния кровеносной системы исследуемого органа, используется доплеровский эффект. Развитие таких систем шло от одноканальных систем [2] к многоканальным [1,3] системам, на которых, в частности, основаны классические УЗИ-системы. Когерентно-аддитивные методы требуют точного априорного знания обрабатываемого сигнала и для получения результата используют операцию согласованной фильтрации. Некогерентные методы не требуют точного знания вида сигналов и основаны лишь на пространственно-временных свойствах их корреляционных функций.

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

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

Более детальное описание процедуры восстановления карты кровеносных сосудов и вектора скорости кровотока заключается в следующем. Используется круговая антенная решетка, где поочередно происходит излучение зондирующих импульсов одним из преобразователей. Каждый излучатель 5 с радиус-вектором г5. последовательно генерирует три или более одинаковых импульса, которые можно объединить в две пары, обозначаемые далее индексами I и II, соответственно. Временной интервал Ггер между первым и

вторым импульсами фиксированной пары одинаков для обеих пар. Рассеянные сигналы м5С, приходящие на приемник Я с радиус-вектором гк и соответствующие фиксированной паре посылок, деструктивно объединяются, формируя два разностных сигнала: Дм^ (г5 ,гд,/) , Дм" (г5 ,гд ,0 , где I -время. Это простейший вариант селекции, подавляющий рассеяние от малоподвижных объектов. В сигналах Дг/'с и Дн" остается информация только от движущихся рассеиветелей. Допустимое полное время измерений при фиксированном излучателе ограничено временем разрушения фиксированного взаимного расположения кластеров крови.

На основе значений Д и51.=Дг/'с и Дм5С = Дм", предварительно сфазированных на текущую точку г, строятся, соответственно, два изображения 2 = 21{ г5,г) и 2 = 2п(г5,г) для каждого фиксированного излучателя (индекс Б). Перебор всех излучателей дает две серии изображений. Например, при аддитивном методе построения этих изображений разностные сигналы объединяются взвешенным суммированием по приемникам:

2 (г5 ,г) = X АИ5С(Гх ,Гд ,/5Л(г)) 3(Гх ,Гя ,Г). (1)

я

Здесь значение /57г(г) равно сумме времени распространения сигнала от излучателя 5 до текущей точки г и времени распространения от точки г до приемника Я. При этом учитывается дополнительная задержка, которая обеспечивает фазирование на максимум огибающей комплексного сигнала Аих. Множители 3(гх,гл,г) компенсируют как геометрическую расходимость волн, так и их поглощение в фоновой среде.

Рис. 1. Мультипликативно-корреляционная схема томографирования системы кровоснабжения.

Область восстановления движущихся рассеивателей

Приемник Л

Приемоизлучающие преобразователи

елекция движущихся рассеивателей

Коррелятор

Усреднение по парам приемников при фиксированном излучателе

Пространственная корреляция спекл-структур

Усреднение по излучателям

Уточнение картины кровеносных сосудов и оценка вектора скорости кровотока

Усреднение по излучателям |

Восстановление картины кровеносных сосудов

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

предварительно, например, методом измерения всех времен пролета с последующим уточнением на основе полных томографических данных [4]. Тем самым, восстановлению вектора скорости кровотока должно предшествовать восстановление полной сглаженной томографической картины распределений скорости звука и поглощения, что осуществляется в рамках той же самой томографической схемы.

Поскольку в построении изображений типа (1) участвуют непосредственно сами сигналы, а не их модули, то эти изображения обладают спекл-структурами, перемещающимися вдоль кровеносных сосудов при движении по ним крови. Такое перемещение определяется пространственно-корреляционным сравнением небольших фрагментов двух изображений Z1 и Z11 при одном и том же излучателе. Коррелируемые фрагменты Z'(rs,r + р,), p¡ и Zn(rs,r + р2), р2 е Р2, имеют общий фиксированный центр г, а

переменный вектор р, или р2 характеризует координаты точек в этих фрагментах относительно центра. Результатом является взаимная корреляционная функция Г(г5, г, Аг), построенная для каждой точки изображения г, где г и Аг - координаты центра и относительного смещения коррелируемых фрагментов:

Г(гя,г,Дг) = J [zI(rs,r + p)-ZI(rs,r)][zi;i0lm.(r.,r + p + Ar)],í/p. (2) pe'í

Здесь Z1 (rs,г) = — | Z'(rs,r + p)í/p-среднее значение меньшего фрагмента

Si ре/]

с площадью S¡. При этом в (2) фигурируют значения большего фрагмента, нормированные на свое среднеквадратичное значение:

vii / i Zn(r?,r+p+Ar)

Z„prml(rj..r+p+Ar)= Смысл этой нормировки

I j |zn(rs, r+p'+Ar)|2£/p' V p'^i

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

Далее функция Г усредняется по всем излучателям при фиксированных г и Аг:

Г(г,Дг)=ХГ(г5,г,Дг) .

5

При каждом фиксированном г ищется вектор сдвига Ar =Armax, соответствующий положительному главному максимуму функции ReT(r, Аг), определяемому как

Гтах (r) ^ Re Г(г,Дгшах) = max Re Г(г, Ar) > 0,

Ar

и Ref(r,Armax)»| Imf(r,Armax)|.

Единственность главного максимума функции Rer(r, Дг) при фиксированном г обеспечивается накоплением большого количества слагаемых. В итоге вектор скорости v(r) в каждой точке изображения оценивается из очевидного соотношения v(r)=Armax/rrep .

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

F(r) = ^Z(rs,r), полученной суммированием промежуточных изображений

s

Z(rs,r).

Условия на время повторения Ггср имеют вид:

Т,ср>^ Vmh/rep^O.

где т5С - длительность рассеянного сигнала; vmin - минимальное значение |v(r)|, которое желательно оценить; Дг0 - минимальное смещение спекл-структур за время повторения зондирующих импульсов, которое можно надежно оценить в используемой томографической системе. При L = (0.14-0.2)м, с0 = 1500 м/с, скорости в достаточно крупных сосудах vmin = (0.05-ь0.2)м/с, Дг0 =Я.0/4 и несущей частоте =1.6 МГц (Я.0 - длина волны на несущей частоте) имеем Ггер s (14)мс, что вполне приемлемо в рассматриваемых томографических устройствах.

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

реального перемещения крови за суммарное время измерений при всех излучателях. В-третьих, процесс измерений при посылках от разных излучателей может быть синхронизован с ритмикой сердца. Тогда можно выделить, например, значение максимальной (за период биений сердца) скорости и направление кровотока в каждой точке изображения. Если для определения скорости кровотока используется излучение одного или небольшого количества излучателей, то весь процесс съема данных приемниками в количестве 200+300 штук занимает несколько миллисекунд. Если же используется весь набор излучателей (200 ч-300 штук), то весь процесс непосредственного съема данных занимает уже до 0.5сек, т.е. сравним с интервалом сердцебиений. Поэтому его следует разбить на несколько групп (по 10-н50 излучателей), привязав моменты посылок излучений в каждой последующей группе к одному и тому же временному интервалу соответствующего нового кардиоцикла. Именно некогерентный характер корреляционной обработки изображений позволяет провести такое разбиение.

Результаты модельных исследований по оценке векторной картины кровотока и первых проведенных экспериментов с движущейся модельной средой подтверждают перспективность описываемого метода в акустической томографии. Так, в разделе 2.2 приводятся результаты численного моделирования процедуры оценки вектора скорости в случае постоянного и переменного распределений вектора скорости кровотока в поперечных сечениях кровеносных сосудов, а также результаты первых экспериментов с движущейся модельной средой. Исходные параметры при численном моделировании выбирались близкими к реальным характеристикам томографа и имеющимся данным по акустическим характеристикам крови. Предполагалось, что в антенной решетке томографа используется равномерное круговое расположение приемоизлучающих преобразователей, хотя, как уже отмечалось, большое количество (200 + 300) равномерно расположенных преобразователей на практике может быть заменено гораздо меньшим количеством преобразователей (20+30), расположенных неравномерно (по специально подобранной схеме), в сочетании с дискретизованным вращением антенной решетки [5].

На рис. 2 приведен результат модельного восстановления аддитивно-корреляционным методом распределения вектора скорости кровотока в случае скорости, постоянной по абсолютной величине |у|, но изменяющейся по направлению (рис. 2а). Предлагаемый способ восстановления позволяет восстановить каждую из компонент вектора скорости (рис. 26): при этом в отсутствие помех хорошо воспроизводится как | V | (рис. 2е), так и собственно векторное поле V (рис. 2г).

Рис. 2. Восстановление вектора скорости кровотока в кровеносном сосуде петлеобразной формы (толщина сосуда скорость |у| =20 см/с; время повторения зондирующих импульсов Ггер = 1 мс; используются 64 преобразователя):

(а)-модель кровеносного сосуда; восстановленные оценки у-компоненты (б) и модуля (в) скорости кровотока; векторное поле скорости кровотока, восстановленное в отсутствие погрешностей в данных рассеяния (г), а также в присутствии помех в виде случайного движения фоновой ткани со скоростью =1А см/с (д) и в виде случайного аппаратурного шума, приводящего к относительным среднеквадратичным

20 см/с

20 см/с

.У^О

Возможность восстановления оценки неравномерного распределения скорости в поперечном сечении кровеносного проиллюстрирована на модели разветвленного кровеносного сосуда (рис. 3). Толщина основной части сосуда составляет ЗА,0, толщина двух ответвлений 2Х0 и /.0 (рис. За). Видно, что пространственное распределение максимумов функции Гтах(г) четко характеризует область локализации движущихся рассеивателей (рис. 36), и качество воспроизведенного вектора скорости кровотока оказывается достаточно хорошим (рис. Зв).

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

Относительную среднеквадратичную ошибку, привносимую по тем или иным причинам (движениями фона или аппаратурными шумами) в исходные рассеянные сигналы и^ (г5 ,гд,?) при фиксированной посылке, можно оценить

как

1

£ Я_

III

5 Я

где и5С=и5(с9)(г5,гл,/=Сх) - рассеянные

сигналы в отсутствие помех и шумов; — момент времени, соответствующий максимуму рассеянного сигнала для данной пары "источник 5-приемник Л";

г/"с'"ж - рассеянный сигнал в присутствии помех или шумов в тот же момент . Соответствующая относительная ошибка в разностных сигналах

составляет аАи =

5 Я

где Аи,с, Д«Г

разностные

рассеянные сигналы в отсутствие и в присутствии помех, рассматриваемые в момент времени, соответствующий максимуму |Ли5с| Ддя данной пары"источник З1-приемник Я".

Рис. 3. Восстановление вектора скорости кровотока, параболически распределенного в поперечных сечениях разветвленного кровеносного сосуда (скорость у изменяется в поперечных сечениях от б см/с до 30 см/с; Ггер = Змс; 64 преобразователя):

(а) - модель кровеносного сосуда;

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

ЗД„

30 см/с

Гтах(г) М

у/Х0

Присутствие помех от смещения малоподвижного фона в кровотоке моделировалось в виде потока рассеивателей, координаты которых флуктуируют в процессе движения. При этом положение каждого фонового рассеивателя изменялось случайным образом (как по величине, так и по направлению) в пределах участка площади % х х, где у_ — заданный линейный размер, составляющий сотые или десятые доли от Х0. Ни при одном из методов восстановления смещение %= ^о /100, приводящее при количестве приемоизлучающих преобразователей К = 32 к относительным среднеквадратичным ошибкам в данных аи ~ 4.4 и аАи ~ 0.44 почти не искажает результата восстановления, по сравнению с отсутствием помех (рис. 2(3). Оказалось, что сильные искажения появляются только в случае X = А.0/30 (аи-5.5, аАи «2.8 при К = 64), что соответствует уже достаточно заметной скорости движения фона уЬв = 2.4см/с, составляющей около 12% от

восстанавливаемой скорости движения кровотока.

Другой тип помех - аппаратурный шум - моделировался в виде нормально распределенного случайного процесса, который привносился независимо к каждой реализации рассеянного поля и3('=1), Здесь, в

отличие от помех, порождаемых малоподвижным фоном, для среднеквадратичных ошибок всегда имеет место неравенство стДц » с„, так как шум в разностных сигналах Аивс не компенсируется, и его относительный уровень становится существенно выше, чем в м8С. При предположениях, близких к реальным, когда отношение сигнал/шум выше 40дБ (ст„=10~2, сДц = 2), результат восстановления остается удовлетворительным (рис. 2е). Таким образом, устойчивость предлагаемого метода к помехам обоих типов приемлема для реализации томографической схемы на практике.

В разделе 2.3 представлены результаты физических экспериментов, подтверждающих работоспособность предлагаемого метода определения скорости кровотока в составе медицинского акустического томографа. Антенная решетка рассматриваемого томографа имеет разреженную структуру и состоит из 26 приемоизлучающих преобразователей, расположенных на кольце неравномерно, однако в строго определенных, заранее рассчитанных угловых позициях. Таким образом, реально имеется 26 излучателей и 26 приемников. Кольцевая решетка совершает дискретизованное круговое движение; угловой шаг дискретизации составляет 2 л/256. В каждом фиксированном угловом положении антенной решетки один преобразователь излучает, и все преобразователи принимают приходящий сигнал; потом излучает следующий преобразователь, и т.д. После перебора всех 26 излучателей происходит поворот в новую угловую позицию антенной решетки. В результате, общий объем экспериментальной информации будет эквивалентен информации в случае неподвижной решетки, содержащей 256 приемоизлучающих преобразователей.

В первом физическом эксперименте в область озвучивания помещалась трубочка из тонкого полиэтилена, имеющая диаметр 10 мм. Трубочка располагалась в плоскости томографирования. По трубочке протекала имитационная жидкость, представляющая собой взвесь зубного порошка в киселе при концентрации около 2 грамм порошка на 200 грамм киселя. Результаты, полученные с помощью обработки аддитивно-корреляционным методом данных, измеренных только при одном фиксированном положении антенной решетки (вместо 256 возможных положений), представлены на рисунке 4. Местоположение трубочки (рис. 46), направление движения и абсолютное значение скорости (рис. 4в, г) оцениваются при восстановлении правильно, что свидетельствует о работоспособности метода.

В следующем эксперименте использовалась медицинская силиконовая трубочка с внутренним диаметром 4 мм и толщиной стенки 1 мм. Условия проведения эксперимента оставались прежними. Результаты обработки данных, полученных для 250 различных угловых положений антенной решетки, представлены на рис. 5. Распределение максимумов функции корреляции (рис. 5а) стало более равномерным в области локализации трубочки, по сравнению с рис. 46. В карте восстановленной скорости (рис. 56) отсутствуют "изъяны", по сравнению с рис. 4 в, г). Таким образом, увеличение количества измеренных данных улучшает результат восстановления, подтверждая оценки, полученные на модельных данных.

Рис. 4. Эксперимент с движением жидкости в трубочке, ориентированной по диагонали (а-в) и вдоль оси ОХ (г); экспериментальные данные получены при единственном угловом положении антенной решетки: фрагмент одного из коррелируемых изображений при фиксированном излучателе (а); распределение максимумов функции корреляции Гтах(г) (б); карта скорости у(г) (в, г).

ReZ'(i"j,«,r)

Гтах(г)

Рис. 5. Эксперимент с движением жидкости в трубочке, ориентированной вдоль оси ОХ ; экспериментальные данные получены для 250 угловых положений антенной решетки: распределение максимумов функции корреляции (а) и карта скорости (б).

— Й(г) — 2.2 см/с

Гщах (Г)

В разделе 2.4 даны выводы ко второй главе.

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

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

Рп

параметр второго порядка е(г) определяется как s = 1 +

2 с20

д2Р Эр2.

> где р0

Р=Ро

невозмущенная плотность среды, с0 - невозмущенная скорость звука, Р -мгновенное давление среды, р - мгновенная плотность среды. В случае рассматриваемой схемы нелинейной томографии второго порядка на окружности располагаются два плоских излучателя и плоский приемник; все преобразователи ориентированы на центральную область этой окружности. Угол между акустическими осями преобразователей достаточно произволен (рис.6). Исследуемый объект располагается внутри области пересечения прожекторных зон преобразователей. Излучаются два кодированных широкополосных неколлинеарных сигнала. Когда нелинейная неоднородность с нелинейным параметром второго порядка е(г) присутствует в области пересечения первичных волн, возникают нелинейные эффекты взаимодействия этих волн, и в результате появляются комбинационные сигналы, которые регистрируются на приемнике. Поскольку вид кодирования первичных сигналов известен априори, то и кодировка сигналов на комбинационных частотах также известна.

Рис. 6. Томографическая схема, предназначающаяся для восстановления пространственных распределений нелинейного параметра и вектора скорости кровотока:

к,, к2-волновые Излучатель 1 "" ^

векторы на частотах />< /2/

регистрируются нелинейно рассеянные комбинационные сигналы с волновым вектором к ± на частотах /, + /2.

Излучатель 2

Приемник

Метод восстановления распределения скорости кровотока v(r) основан на определении вектора перемещения спекл-структуры, возникающей при оценке картины распределения нелинейного параметра движущейся крови, путем пространственно-корреляционного сравнения спекл-структур. Для этого на первом этапе, путем согласованной фильтрации комбинационных сигналов по методике [6], строятся изображения s(r), каждое из которых является оценкой распределения значений нелинейного параметра в среднем за короткий интервал времени и имеет вид спекл-структуры. Для практической реализации весьма важен вопрос о компенсации сигнала, рассеянного от неподвижной окружающей биоткани. Трудность заключается в том, что компенсация должна осуществляться для рассеянного сигнала, порождаемого однократным излучением первичных сигналов большой длительности (свыше 0.1 с). Поэтому предлагается использовать следующую процедуру согласованной фильтрации в сочетании с селекцией движущихся рассеивателей. Регистрируемый на плоском приемнике с радиус-вектором центра у комбинационный сигнал р(у, t) с

длительностью Т разбивается на короткие интервалы {р'''(у, /)}, равной длительности Г5иЬ; i - номер текущего интервала. По тому же самому принципу эталонный отклик р6(у|г;/) разбивается на {pg (у | г;/)}, . Эталонный отклик />5(у|г;г) для точки пространства г - это комбинационный сигнал на приемнике, порождаемый неподвижным пробным точечным рассеивателем, находящимся в точке г. Согласованная фильтрация сигнала pU){yj) осуществляется для каждого фиксированного интервала i. Используется фильтр, осуществляющий одновременно согласованную фильтрацию и

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

U(y|r;/), t е[0, ТьиЪ/2] Рй\(У I г; /) = <! .В результате отражение от

^sub]

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

нелинейного параметра второго порядка: е '(г) = —-=-^— .

J/$(y|r;/)|^(y|r;/)J dt

Следующим этапом обработки является пространственная корреляция фрагментов пар изображений {sj'fr), ¿^'(г)} за фиксированный период времени между начальным i -м и j -м интервалом. Процедура корреляции

аналогична используемой в линейных томографических системах, - см. (2). Результатом корреляции фрагментов при их взаимном сдвиге Ar является двумерная взаимная корреляционная функция Г^(г,Дг), построенная для окрестности каждой точки изображения г. Далее функция Г^'-,'(г,Аг) усредняется по множеству коррелируемых пар, соответствующих одному и тому же Ггер. Для каждого фиксированного г ищется вектор сдвига Ar = Armax,

соответствующий главному максимуму Гтах(г) этой усредненной функции. Область пространственной локализации функции Гтах(г) определяет положение кровеносного сосуда и его форму. В итоге, внутри сосуда искомый вектор скорости оценивается как v(r) = Агшах/Ггср . В предлагаемом подходе

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

В разделе 3.2 приводятся результаты численного моделирования. Одновременно показано, что аддитивное объединение спекл-структур, возникающих при различных углах приема или при различных положениях приемоизлучающей системы как целой не требует априорного знания карты сосудов для обеспечения определенной взаимной ориентации излучателей, приемника и сосуда. Первичные плоские волны моделировались в виде сигналов, фазоманипулированных по псевдослучайному закону; /¡° = 2 МГц,

/2° = 1.8 МГц - средние частоты, и /, е (1 ч-3) МГц, /2 е (0.9 н-2.7) МГц -излучаемые частотные полосы. Акустические оси двух плоских излучателей составляют, соответственно, угол 120° и 240° по отношению к оси X. Ось X совпадает с биссектрисой угла между направлениями излучения первичных волн и является линией одинакового кода регистрируемых рассеянных сигналов на суммарных комбинационных частотах. Кровеносный сосуд был ориентирован вдоль оси X и имел ширину 4Х°+ (где - длина волны на средней суммарной комбинационной частоте. Принимаемый комбинационный сигнал обрабатывается в полосе (0.5 1.5)/+° с центральной суммарной частотой /+° =3.8 МГц.

Нелинейное рассеяние моделировалось для движущихся рассеивателей (кластеры крови) и, одновременно, для неподвижных фоновых рассеивателей (биологическая ткань). Важно, что неподвижные рассеиватели помещались не только вне кровеносного сосуда, но и внутри него, поскольку озвученный при томографировании слой имеет конечную толщину. Диапазон случайного распределения величины нелинейного параметра задавался одинаковым (от 3.5 до 4.5) для элементарных рассеивателей внутри сосуда и для фоновой ткани. Имитировалось равномерное движение эритроцитов со скоростью | v |= 20 см/с. Принятый комбинационный сигнал длительности Г = 40мс разбивался на 80 интервалов длительности Г5иЬ = 0.5 мс, для каждого из которых строились

промежуточные изображения с(,)(г). Особенности ориентационной

изменчивости изображений иллюстрируются на рис.7а, 16. Анизотропия высокочастотного пространственного заполнения спекл-структур порождена анизотропией разрешающей способности [6]. Для последующего модельного восстановления вектора скорости использовался угол приема 190° (рис.7б), при котором характерные размеры элементов спекл-структуры достаточно малы. Итоговая оценка векторного поля скорости у(г) (рис.7в) хорошо воспроизводит истинное распределение. В процедуре компенсации влияния фона важны одновременно оба фактора - селектирующий фильтр (без селекции оценка у(г) получается неудовлетворительной — см. рис.7г) и накопление результатов пространственной корреляции по парам изображений.

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

В разделе 3.3 представлены выводы к третьей главе.

Рис. 7. Результат восстановления прямолинейном сосуде толщиной томографии:

спекл-структура, возникающая после селективной согласованной фильтрации комбинационного сигнала для углов излучения 120° и 240° при угле приема 100° (а), 190° (б); восстановленное векторное поле скорости с учетом селекции движущихся рассеивателей (в) и без селекции (г); (д) — аддитивное объединение спекл-структур, показанных на рис. а и б.

б

у(г) 20 см/с »

- —— - - ■

в

у/1!

у(г) 50 см/с

2.5

-2.5

I У У •

— ^ ✓ 4

. — — X -

-2.5

2.5 5

г

вектора скорости кровотока в 4А? на основе данных нелинейной

Приемник

-6 -4 -2 0 2 4 6 8

Излучатель 2 х / А,+

Приемник /излучатель I

у'А

Приемник'

\'8 -4 О

Излучатель 2

д

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

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

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

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

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

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

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

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

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

Список цитируемой литературы

1. DuricN., LittrupP., Chandiwala-Mody P., Li С., SchmidtS., Мус L., Rama О., Bey-Knight L„ Lupinacci J., Ranger В., Szczepanski A., West E. In-vivo imaging results with ultrasound tomography: Report on an ongoing study at the Karmanos Cancer Institute // Medical Imaging 2010: Ultrasonic Imaging, Tomography, and Therapy. Eds. D'hooge J., McAleavey S.A., 2010. V. 7629. P. 76290M.

2. FerraraK.W., ZagarB. High resolution 3D color flow mapping of breast vasculature // IEEE Ultrasonic symposium. 1995. P.1467-1470.

3. Svensson W.E., Humphries K., Stanley P., Holt P., Cosgrove D.O., Forouhi P. Combined/split three-dimensional greyscale/colour Doppler imaging shows detail of breast lesions vascular morphology // Ultrasound in Medicine & Biology. 2003. V.29. P.105-106.

4. Burov V.A., Zotov D.I., Karavay M.F., Rumyantseva O.D. Ultrasound tomography of soft biological tissues containing strong inhomogeneities // Physics of Wave Phenomena. 2013. V. 21. N 1. P. 74-80.

5. Пархоменко П.П., Каравай М.Ф. и др. Ультразвуковой томограф и кольцевая антенная решетка для ультразвукового томографа // Патент на изобретение № 2145797. Приоритет от 23.06.1999. Москва, 2000.

6. Буров В.А., Евтухов С.Н., Ткачева A.M., Румянцева О.Д. Акустическая томография нелинейного параметра с помощью малого числа преобразователей // Акустич. журн. 2006. Т.52. №6. С.760-776.

Список работ, опубликованных по теме диссертации

А1. Буров В.А., Евтухов С.Н., Матвеев О.В., Румянцева ОД. Методы и возможности некогерентной корреляционной акустической томографии // Биомедицинские технологии и радиоэлектроника. 2005. №4-5. С.55-63. А2. Буров В.А., Матвеев О.В., Нестерова Е.В., Румянцева О.Д. Восстановление картины распределения вектора скорости кровотока в процессе акустического томографирования //Акуст. журн. 2006. Т.52. №5. С.607-623.

A3. Буров В.А., Евтухов С.Н., Матвеев О.В., Румянцева О.Д. Восстановление вектора скорости кровотока в процессе томографирования акустического нелинейного параметра // Известия Российской Академии Наук. Серия Физическая. 2006. Т.70. №12. С.1694-1699.

А4. Буров В.А., Матвеев О.В., Румянцева О.Д. Пространственно-корреляционный метод восстановления вектора скорости кровотока в нелинейной томографии //Акустич. журн. 2010. Т.56. №2. С.268-276.

А5. Буров В. А., Матвеев О.В., Нестерова Е.В. Корреляционно-томографическое восстановление вектора скорости кровотока // Сборник трудов XV сессии Российского Акустического Общества. Т.З. М: ГЕОС, 2004. С. 83-87.

А6. Бурое В.А., Матвеев О.В., Нестерова Е.В. Восстановление распределения вектора скорости кровотока в акустическом томографе // II Евразийский конгресс по медицинской физике и инженерии "Медицинская физика-2005" 21-24 июня 2005. Сборник материалов. М.: Ассоциация медицинских физиков России. 2005. С.215-216.

А7. Матвеев О.В. Восстановление картины кровотока в процессе активного корреляционного томографирования // Международная конференция студентов, аспирантов и молодых ученых по фундаментальным наукам "Ломоносов-2003". Секция "Физика". Сборник тезисов, М.: Физический факультет МГУ, 2003. С. 117-119.

А8. Матвеев О.В., Нестерова Е.В. Корреляционно-томографическое восстановление распределения вектора скорости кровотока // Международная конференция студентов, аспирантов и молодых ученых по фундаментальным наукам "Ломоносов-2004". Секция "Физика". Сборник тезисов, М.: Физический факультет МГУ, 2004. С. 45-48.

А9. Burov V.A., Matveev О. V., Evtukhov S.N., Rumyantseva O.D. Reconstruction of spatial distribution of bloodstream vector in ultrasound linear and nonlinear tomographs // PSFVIP-8: The 8th Pacific Symposium on Flow Visualization and Image Processing. Ed. I.A.Znamenskaya. Moscow: Moscow State University, 2011. PSFVIP8-006.

A10. Буров В.А., Евтухов C.H., Матвеев О.В., Ткачева A.M. Томографическое восстановление распределения нелинейного параметра в биологических средах // Сборник трудов XVI сессии Российского Акустического Общества. Т.З. М: ГЕОС, 2005. С.И5-118.

All. Буров В.А., Евтухов С.Н., Матвеев О.В. Восстановление вектора скорости кровотока в процессе томографирования нелинейного параметра // Труды X Всероссийской школы-семинара "Волновые явления в неоднородных средах". Звенигород, Московская область, 22-27 мая 2006. Секция 1 "Распространение акустических и гидродинамических волн". М.: физич. ф-т МГУ, 2006. С.36-38.

А12. Burov V.A., Matveev O.V., Evtukhov S.N., Rumyantseva O.D. Reconstruction of blood velocity vector in nonlinear acoustical tomography // Acoustical Imaging. Eds. Nowicki A., Litniewski J., Kujawska T. Dordrecht, Heidelberg, London, New York: Springer. 2012. V. 31. P. 211-221.

Матвеев Олег Владимирович Восстановление распределения вектора скорости кровотока в линейном и нелинейном акустических томографах Формат 60x90/16 Тираж 100 экз. Усл.п.л. 1.75 Подписано в печать 23.05.2013 Заказ № 182 Типография ООО «Генезис» 8 (495) 434-83-55 119571, г. Москва, пр-т Вернадского, 86

 
Текст научной работы диссертации и автореферата по физике, кандидата физико-математических наук, Матвеев, Олег Владимирович, Москва

МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

им. М.В. Ломоносова

Физический факультет, кафедра акустики

На правах рукописи УДК 534.2 : 534.7 : 534.222

04201458369

МАТВЕЕВ Олег Владимирович

ВОССТАНОВЛЕНИЕ РАСПРЕДЕЛЕНИЯ ВЕКТОРА

СКОРОСТИ КРОВОТОКА В ЛИНЕЙНОМ И НЕЛИНЕЙНОМ АКУСТИЧЕСКИХ ТОМОГРАФАХ

Специальность: 01.04.06 - акустика

Диссертация на соискание ученой степени кандидата физико-математических наук

Научные руководители:

доктор физ.- мат. наук, проф. Буров Валентин Андреевич

кандидат физ.- мат. наук Румянцева Ольга Дмитриевна

Москва, 2014

СОДЕРЖАНИЕ

Глава 1. ВВЕДЕНИЕ...........................................................................................

1.1. Актуальность темы и цели исследования. Структура диссертационной работы...................................................................................................................

1.2. Общее состояние проблемы (по материалам научных публикаций)......

Глава 2. ВОССТАНОВЛЕНИЕ РАСПРЕДЕЛЕНИЯ ВЕКТОРА СКОРОСТИ КРОВОТОКА В ПРОЦЕССЕ АКТИВНОГО КОРРЕЛЯЦИОННОГО ТОМОГРАФИРОВАНИЯ..........................................................................................25

2.1. Аддитивно-корреляционный и мультипликативно-корреляционный алгоритмы восстановления общей картины распределения кровотока и вектора его скорости.................................................................................................................25

2.2. Численное моделирование корреляционно-томографических методов восстановления вектора скорости кровотока...........................................................45

2.2.1. Прямая задача рассеяния от малоподвижной фоновой среды и движущихся объектов...............................................................................................45

2.2.2. Сравнительный анализ результатов восстановления аддитивно-корреляционным и мультипликативно-корреляционным методами....................52

2.2.3. Помехоустойчивость решения при влиянии медленных движений фоновой среды и аппаратурных шумов.................................................................................118

2.2.4. Сравнение методов восстановления скорости кровотока.......................123

2.3. Экспериментальное обнаружение временного сдвига сигналов, рассеянных на кластеризованной движущейся жидкости, и восстановление карты кровотока .....................................................................................................................................124

2.4. Выводы к главе 2...............................................................................................133

Глава 3. ПРОСТРАНСТВЕННО-КОРРЕЛЯЦИОННЫЙ МЕТОД ВОССТАНОВЛЕНИЯ ВЕКТОРА СКОРОСТИ КРОВОТОКА В СХЕМЕ

НЕЛИНЕЙНОЙ ТОМОГРАФИИ...........................................................................134

3.1. Пространственно-корреляционная процедура восстановления скорости рассеивателей............................................................................................................134

3.2. Численное моделирование пространственно-корреляционного метода......142

3.3. Выводы к главе 3...............................................................................................153

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ............................................................154

СПИСОК ЛИТЕРАТУРЫ........................................................................................156

Глава 1. ВВЕДЕНИЕ

Работа посвящена получению численной оценки полного вектора скорости

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

1.1. Актуальность темы и цели исследования. Структура диссертационной

работы.

Актуальность темы

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

Большинство появившихся в последние десятилетия систем медицинской акустоскопии - ультразвуковые многоканальные системы активной локации (УЗИ) и разрабатывавшиеся в последние годы многоракурсные томографические системы [1—4; 62, 67, 82]. Эти диагностические комплексы представляют собой, как правило, линейные когерентные системы активной локации, основанные на аддитивном (когерентном) сложении сигналов, принимаемых элементами антенной решетки. Важным достоинством этих систем является линейная зависимость выходного отношения сигнал/помеха от такого же входного отношения. В когерентных системах, применяемых для диагностики состояния кровеносной системы исследуемого органа, используется доплеровский эффект. Развитие таких систем шло от одноканальных систем [5-8] к многоканальным [9-13,49, 60-72, 82] системам, на которых, в частности, основаны классические УЗИ-системы. Когерентно-аддитивные методы требуют точного априорного знания обрабатываемого сигнала и для получения результата

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

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

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

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

При достижении этой цели в работе были поставлены и решены следующие задачи:

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

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

3. Оценка помехоустойчивости алгоритмов и границ их работоспособности.

4. Проведение модельных численных расчетов, имитирующих физический эксперимент.

5. Проведение физического эксперимента.

Научная новизна.

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

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

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

Значимость

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

Вклад автора

Все результаты получены автором лично на основе модельных данных, полученных по методам, разработанным автором.

Достоверность.

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

Положения, выносимые на защиту

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

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

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

Апробация

Материалы докладывались на следующих конференциях:

• Международная конференция студентов, аспирантов и молодых ученых по фундаментальным наукам "Ломоносов-2003". Секция "Физика".

• Международная конференция студентов, аспирантов и молодых ученых по фундаментальным наукам "Ломоносов-2004". Секция "Физика".

• II Евразийский конгресс по медицинской физике и инженерии "Медицинская физика-2005" 21-24 июня 2005.

• X Всероссийская школа-семинар "Волновые явления в неоднородных средах". Звенигород, Московская область, 22-27 мая 2006.

• The 8th Pacific Symposium on Flow Visualization and Image Processing. Ed. I.A.Znamenskaya. Moscow: Moscow State University, 2011 (8-й мировой симпозиум по визуализации потоков и вычислению изображений им. И.А. Знаменской, МГУ, 2011).

Публикации

Результаты диссертации опубликованы в двенадцати работах [ 16, 37, 53, 73, 74, 75, 76, 77, 78, 79, 88, 95], в том числе, четыре статьи опубликованы в журналах, включенных в рекомендованный список ВАК РФ [ 16, 53, 78, 79].

Структура диссертационной работы

Работа состоит из введения, двух глав, заключения и списка литературы. Список литературы состоит из 102 наименований. Общий объем работы составляет 165 страниц, включая 32 рисунка.

В первой главе в разделе 1.1 приведено обоснование актуальности выбранной темы диссертационной работы (далее - работы), показана цель работы и поставлены задачи, решенные по ходу работы. Приведены формулировки научной новизны и значимости работы. Показан вклад автора, и дана формулировка, подтверждающая достоверность результатов работы. Раскрыты выносимые на защиту положения, перечислены конференции и публикации, где докладывались и изложены результаты работы. В разделе 1.2 представлен обзор научных публикаций по теме работы. Описано развитие современных методов УЗ диагностики от одноканальных систем до многоканальных систем, как линейных, так и нелинейных. Показаны связь и различие доплеровских средств диагностики и корреляционных методов.

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

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

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

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

1.2. Общее состояние проблемы (по материалам научных публикаций).

В схемах традиционных звуковых систем при одностороннем озвучивании лишь отчасти удается получить оценку полного вектора скорости за счет пространственной модуляции поперечной структуры диаграмм направленности [17, 18, 19]. При этом измерение продольной компоненты скорости (за счет эффекта Доплера) и поперечной компоне�