Способ выявления аномалий физической плотности при акустической визуализации

Авторы патента:

Громов Александр Игоревич (RU)

Кульберг Николай Сергеевич (RU)

Фин Виктор Александрович (RU)

Леонов Денис Владимирович (RU)

A61B8/08 - обнаружение органических изменений в организме, например различных новообразований, кист, опухолей (A61B 8/02-A61B 8/06 имеют преимущество)

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

Государственное бюджетное учреждение здравоохранения города Москвы "Научно-практический центр медицинской радиологии Департамента здравоохранения города Москвы" (ГБУЗ "НПЦМР ДЗМ") (RU)

Изобретение относится к обнаружению аномалий физической плотности. Техническим результатом является повышение достоверности акустической визуализации. Способ заключается в том, что проводят измерения сигнала в режиме цветового доплеровского картирования, причем на основе полученных значений измеренного сигнала вычисляют модуль коэффициента корреляции действительной и мнимой частей отсчетов сигнала; выбирают порог по модулю коэффициента корреляции действительной и мнимой частей отсчетов сигнала; путем сравнения с порогом строят маску, которая принимает ненулевые значения на участках, где модуль коэффициента корреляции действительной и мнимой частей отсчетов сигнала превышает некоторый порог; при отображении аномалий физической плотности на экране прибора цветом кодируют энергию доплеровских сигналов, для которых маска принимает ненулевые значения, то есть сигналов, соответствующих упругим колебаниям, исходящим от аномалий физической плотности. 2 н. и 5 з.п. ф-лы, 4 ил.

Область техники, к которой относится изобретение

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

Уровень техники

В ультразвуковой медицинской диагностике для обнаружения аномалий физической плотности часто используется В-режим. При этом искомые объекты представляют собой яркие области на экране ультразвукового прибора, за которыми может быть видна акустическая тень. Но точность такого метода обнаружения невелика [1]. Для увеличения точности известно использование режима цветового доплеровского картирования кровотока. При этом на аномалиях физической плотности наблюдается мерцающий артефакт - явление, проявляющееся в быстрой смене окрашенных пикселей в окрестности гиперэхогенного объекта. Согласно статистическим данным применение мерцающего артефакта может повысить вероятность истинного обнаружения кальцинатов на 20-37% [1].

Возможно и дополнительное повышение вероятности обнаружения объекта, если использовать добавочный источник низкочастотных акустических колебаний. Эта методика называется акустической резонансной визуализацией [2, 3] и заключается в возбуждении резонанса объекта под действием сигнала на частоте, близкой к частоте их собственных колебаний. В патенте [2] предлагается дооснастить ультразвуковой медицинский диагностический прибор генератором гармонических колебаний для возбуждения микрокальцинатов в режиме цветового доплеровского картирования. Это приведет к увеличению интенсивности мерцающего артефакта.

Однако режим цветового доплеровского картирования в медицинских устройствах предназначен для картирования кровотока. Мерцающий артефакт является скорее «случайной находкой». Часто он мешает визуализации кровотока. Многие производители ультразвуковых медицинских диагностических приборов стремятся от него избавиться. Поэтому разумным выходом является создание режима, предназначенного для обнаружения аномалий физической плотности [4, 5].

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

Раскрытие изобретения

Задачей изобретения является повышение информативности акустической визуализации при поиске аномалий физической плотности в исследуемом объекте.

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

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

Описание чертежей

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

Осуществление изобретения

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

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

где - это ансамбль, состоящий из N комплексных входных отсчетов, называемый также пачкой. Для получения ансамбля исследуемый объект облучается N раз с некоторым интервалом T_PRF. k - номер импульса в пачке.

A_k - комплексная амплитуда сигналов от быстродвижущихся объектов (например, кровотока), - доплеровский сдвиг частоты, обусловленный движением объекта со скоростью ν_A. На фиг. 1 дано представление сигнала, содержащего только компоненту A_k, на комплексной плоскости;

B_k - комплексная амплитуда сигналов от малоподвижных объектов, которая, как правило, на 20-30 дБ превосходит амплитуду сигналов от кровотока.

- доплеровский сдвиг частоты, обусловленный движением тканей и колебаниями датчика. Предполагается, что скорость этого движения ν_B ≪ ν_A ;

C_k - амплитуда сигнала от аномалии физической плотности, мощность которого на 0-10 дБ превосходит амплитуду компоненты B_k. Предполагается, что аномалия физической плотности может совершать движения вместе с окружающими ее малоподвижными объектами. Следовательно, ей также соответствует доплеровский сдвиг частоты, равный ω_B. Помимо этого, аномалия физической плотности может совершать вынужденные колебания под действием давления сканирующих импульсов, это приводит к дополнительному изменению фазы сигнала, обозначенному здесь ϕ_k. Для оценки величины этого сдвига мы используем следующую формулу:

где R ≪ - амплитуда вынужденных колебаний аномалии физической плотности вдоль луча; ω_C - частота этих колебаний, обычно она сравнима с частотой повторения импульсов пачки. На фиг. 2 дано представление сигнала, содержащего только компоненту C_k, на комплексной плоскости;

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

E_k - комплексная компонента, характеризующая тепловой шум и шумы квантования и дискретизации в тракте приема и первичной обработки сигнала. Ее амплитуда на 10-15 дБ меньше амплитуды сигнала A_k;

ϕ₀ - начальная фаза последовательности.

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

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

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

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

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

Сигналы от аномалий физической плотности C_k и D_k содержат компоненты, которые также не подавляются ФВЧ, предназначенным для исключения малоподвижных объектов. При этом они имеют достаточную мощность, чтобы пройти маску (3), изначально рассчитанную на обнаружение только компоненты A_k. В результате, наряду с картированием движущихся участков, в доплеровском режиме на изображении появляются дополнительные окрашенные области. Именно это явление получило название мерцающего артефакта.

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

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

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

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

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

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

D_k и повысит вероятность обнаружения объекта.

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

Источники информации

1. Громов А.И., Кубова С. Ю. Ультразвуковые артефакты. М.: Видар, 2007.

2. Seghal С. Apparatus for imaging an element within a tissue and method therefor. United States Patent №5,997,477. 1999.

3. Weinstein S.P., Seghal C, Conant E.F., Patron J.A. Microcalcifications in Breast Tissue Phantoms Visualized with Acoustic Resonance Coupled with Power Doppler US: Initial Observations. Radiology. July 2002: 265-269.

4. Lu W. Ultrasonic Detection and Expulsion of Kidney Stones [dissertation]. Seattle: Department of Bioengineering, University of Washington, 2012.

5. Леонов Д.В., Кульберг H.C., Громов А.И., Морозов С.П., Ким С.Ю. Исследование причин возникновения мерцающего артефакта в доплеровских режимах ультразвукового медицинского диагностического устройства// Акустический журнал. 2018. №1. С. 100-111.

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

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

- выбирают порог по модулю коэффициента корреляции действительной и мнимой частей отсчетов сигнала;

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

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

3. Способ по п. 1 или 2, отличающийся тем, что при картировании используют лишь один оттенок, окрашивающий часть кадра, в которой маска принимает ненулевые значения.

4. Способ по п. 1 или 2, отличающийся тем, что картируемым параметром является среднее значение частоты колебаний аномалии физической плотности.

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

- на основе полученных значений измеренного сигнала вычисляют среднеквадратичное отклонение отсчетов сигнала;

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

- выбирают порог по модулю коэффициента корреляции действительной и мнимой частей отсчетов сигнала;

- выбирают порог по среднеквадратичному отклонению отсчетов сигнала;

- выбирают порог по модулю коэффициента попарной корреляции отсчетов сигнала;

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

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

6. Способ по п. 5, отличающийся тем, что при картировании используется лишь один оттенок, окрашивающий часть кадра, в которой маска принимает ненулевые значения.

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

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

Сегментация крупных объектов из нескольких трехмерных видов // 2663649

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

Ультразвуковая рама для головы для экстренной медицинской помощи // 2663646

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

Ультразвуковая рама для головы для экстренной медицинской помощи // 2663646

Способ компрессионной ультразвуковой эндоректальной эластографии новообразований прямой кишки // 2662901

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

Способ проведения пункционно-аспирационной и трепан-биопсии опухолей мягких тканей под ультразвуковым контролем // 2662647

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

Способ диагностики гемартроза коленного сустава // 2661864

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

Способ многоимпульсной эластографии // 2659627

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

Система ультразвуковой визуализации и способ ультразвуковой визуализации // 2659021

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

Система трехмерной ультразвуковой визуализации // 2657855

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