Разработка алгоритмов и устройств улучшения характеристик ультразвуковых эхотомоскопов для контроля и диагностики биологических объектов тема диссертации и автореферата по ВАК РФ 05.11.13, кандидат технических наук Нагулин, Сергей Николаевич

• Нагулин, Сергей Николаевич
• кандидат технических наук
• 2007

Оглавление диссертации кандидат технических наук Нагулин, Сергей Николаевич

Принятые сокращения и обозначения.

Введение.

Глава 1. Методы и средства построения ультразвуковых диагностических систем для медицинской диагностики.

1.1 Режимы зондирования и формирования ультразвукового изображения

1.2 Методы обработки сигналов в приемном тракте эхотомоскопа.

1.2.1 Формирование полутоновых ультразвуковых изображений.

1.2.2 Доплеровские измерения.

1.2.3 Цветовое доплеровское картирование.

1.3 Цифровое преобразование формата ультразвуковых изображений.

1.4 Программно-аппаратные средства ультразвуковых эхотомоскопов.

1.5 Выводы по главе 1.

Глава 2. Разработка и исследование методов режекции сигналов движущихся тканей.

2.1 Методы режекции сигналов движущихся тканей, используемые в ультразвуковых системах ЦДК.

2.2 Оценка параметров полезного сигнала и помехи в шумах методом максимальной энтропии.

2.3 Использование метода максимальной энтропии для режекции помехи в режиме ЦДК.

2.3.1 Применение спектральной оценки ММЭ для оценки параметров сигнала движущейся жидкой среды.

2.3.2 Оценка параметров сигнала движущейся жидкой среды на основе линейной экстраполяции вне интервала наблюдения.

2.3.3 Режекция помехи во временной области с помощью ММЭ.

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

Глава 3. Методы построения тракта обработки ультразвукового сигнала на основе итерационных алгоритмов.

3.1. Итерационный алгоритм поворота вектора в параметризованной системе координат.

3.1.1. Общий вид алгоритма ИПВ.

3.1.2 Сходимость алгоритма.

3.1.3. Вычисление элементарных функций с помощью алгоритма ИПВ.

3.2 Методическая ошибка алгоритма ИПВ.

3.3 Применение алгоритма ИПВ для построения цифрового приемного тракта ультразвукового сканера.

3.3.1 Этапы обработки ультразвукового эхо-сигнала с использованием операции поворота вектора.

3.3.2. Реализация операции поворота вектора.

3.3.3. Выполнение операции детектирования.

3.3.4. Сжатие динамического диапазона.

3.3.5. Преобразование форматов ультразвуковых изображений на основе алгоритма ИПВ.

3.4. Выводы по главе 3.

Глава 4. Реализация ультразвуковых комплексов с расширенной компьютерной обработкой.

4.1. Особенности построения ультразвуковых приборов с расширенной компьютерной обработкой.

4.2. Многофункциональный ультразвуковой эхотомоскоп «Сономед-500»

4.2.1. Фронтальный модуль.

4.2.2. Канал визуализации.

4.2.3. Канал спектральной обработки доплеровских сигналов.

4.2.4. Канал цветового доплеровского картирования.

4.3. Компенсация задержек системы фокусировки и сшивка фрагментов изображений в режиме динамического фокуса.

4.4. Фетальный монитор «Сономед-200».

4.4.1 Алгоритм вычисления частоты пульсаций.

4.4.2 Методы тестирования алгоритмов расчета ЧСС для ультразвукового фетального монитора.

4.5. Выводы по главе 4.

Введение диссертации (часть автореферата) на тему «Разработка алгоритмов и устройств улучшения характеристик ультразвуковых эхотомоскопов для контроля и диагностики биологических объектов»

В настоящее время для диагностики и контроля качества движущейся жидкости и биологических объектов в нефтехимии, технике и медицине, необходимо измерять линейную скорость движения жидкой среды, частоту пульсаций тканей, а также исследовать внутреннюю структуру биологических объектов. Для решения этих задач широко применяется ультразвуковой метод, на основе которого создано большое количество ультразвуковых эхотомоскопов таких известных фирм как Philips, General Electric, Siemens, Toshiba и др.

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

Также объектом контроля является внутренняя структура тканей биологических объектов. При выполнении сканирования структур объектов с применением конвексных ультразвуковых датчиков, осуществляется преобразование формата изображения из полярных координат в декартовы для отображения на мониторе. При этом наилучшее качество ультразвукового изображения обеспечивает метод R-0 интерполяции. Проблемы его реализации связаны с необходимостью вычисления нелинейных функций вида л[х, arctg(x), Их, которые реализуются только табличным методом с использованием элементов памяти большого объема. В связи с этим актуальна задача по разработке способов уменьшения используемых вычислительных ресурсов и сокращения объемов памяти при реализации алгоритма R-Q интерполяции на программируемых логических интегральных схемах (ПЛИС).

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

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

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

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

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

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

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

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

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

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

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

Научная новизна настоящей работы заключается в следующем:

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

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

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

- Разработан алгоритм преобразования формата ультразвукового изображения методом итерационного поворота вектора;

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

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

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

Практическая значимость и реализация результатов работы.

Впервые в нашей стране с участием автора создан многофункциональный эхотомоскоп с режимом цветового доплеровского картирования «Сономед-500» и фетальный монитор с автоматическим анализом кардиотокограмм «Сономед-200», о чем имеется акт внедрения. Данные модели ультразвуковых приборов успешно прошли приемочных технические испытания, клинические испытания, зарегистрированы в Министерстве здравоохранения и социального развития и серийно производятся на предприятии ЗАО «Спектромед» (г. Москва).

Апробация работы. Основные результаты работы докладывались и обсуждались на VII Международной научно-практической конференции «Фундаментальные и прикладные проблемы приборостроения, информатики, экономики и права» (Сочи, 2004 г.), II евразийском конгрессе "Медицинская физика-2005" (Москва, 2005 г.), 9-й Международной научно-технической конференции «Цифровая обработка сигналов и ее применение» - DSPA-2007 (Москва, 2007 г.).

Публикации. По теме диссертации опубликовано 7 печатных работ, из которых 4 - в изданиях, включенных в перечень ВАК.

Структура и объём работы. Диссертационная работа состоит из введения, четырёх глав, заключения, 5 приложений и списка литературы, содержащего 94 источника. Материал диссертации изложен на 136 стр. машинописного текста, содержит 6 таблиц, 48 рисунков.

Заключение диссертации по теме «Приборы и методы контроля природной среды, веществ, материалов и изделий», Нагулин, Сергей Николаевич

4.5. Выводы по главе 4

1. При мелкосерийном производстве оборудования использование микросхем - контроллеров шин сторонних производителей является более предпочтительным по сравнению с реализацией интерфейса в ПЛИС. Для ультразвуковых приборов класса «СОНОМЕД - 500» оптимальным выбором интерфейса взаимодействия с ПК является шина PCI, а для приборов класса

СОНОМЕД - 200» - USB 2.0. Это связано с тем, что данные шины обеспечивают необходимую скорость пропускания, а также тем, что существует большой выбор контроллеров этих шин и средств разработки.

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

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

4. Разработанная методика тестирования позволят сравнивать фетальные мониторы как с применением шаблонных (используя фетальный симулятор), так и с применением реальных записей и определять качество формирования кривой ЧСС без вмешательства в их электронику.

5. Прибор «СОНОМЕД - 200», построенный по технологии расширенной компьютерной обработки, обладает характеристиками, не уступающими приборам ведущих производителей, построенным без использования в своей основе ПК, а также дополнительным набором сервисных возможностей.

Заключение

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

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

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

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

Разработанный алгоритм расчета ЧСС для применения в фетальном мониторе основан на вычислении автокорреляционной функции. Выбор времени накопления сигнала для вычисления АКФ является серьезной практической проблемой. Для подавления шумов, присутствующих во входном сигнале, целесообразно увеличивать время накопления сигнала, но для отслеживания быстрых изменений периода сигнала, необходимо его сокращать. Разработанный алгоритм позволяет отслеживать кратковременные изменения частоты пульсаций и обеспечивает устойчивую работу при наличии в сигнале шумов и артефактов. Применение разработанных методов тестирования фетальных мониторов, позволило оптимизировать параметры алгоритма и обеспечить характеристики, не уступающие приборам ведущих мировых производителей.

Результаты диссертационной работы успешно внедрены и используются в приборах «Сономед-500» (многофункциональный эхотомоскоп с режимом цветового доплеровского картирования) и «Сономед-200» (фетальный монитор с автоматическим анализом кардиотокограмм). Данные модели ультразвуковых приборов успешно прошли приемочных технические испытания, клинические испытания, зарегистрированы в Министерстве здравоохранения и социального развития и серийно производятся на предприятии ЗАО «Спектромед» (г. Москва).

Список литературы диссертационного исследования кандидат технических наук Нагулин, Сергей Николаевич, 2007 год

1. Домаркас В. Й., Пилецкас Э.Л. Ультразвуковая эхоскопия. -Л.: Машиностроение. 1988

2. Осипов Л.В. Ультразвуковые диагностические приборы: Практическое руководство для пользователей. М.: Видар, 1999.

3. Лелюк В.Г., Лелюк С.Э. Ультразвуковая ангиология. -М.: Реальное время, 1999

4. Физика визуализации изображений в медицине: в 2-х томах. Т. 2: Пер. с англ. / Под ред. С. Уэба. М.: Мир, 1991.

5. Wells P.N.T., Ziskin Marvin С., "New Techniques and Instrumentation in Ultrasonography", Churchill Livingstone Inc., 1980.

6. Reneman Robert S., Hoeks Arnold P.G., "Doppler Ultrasound in the Diagnosis of Cerebrovascular Disease", John Wiley & Sons Ltd., 1982.

7. Evans D.H., "Doppler Ultrasound: Physics, Instrumentation, and Clinical Applications", John Wiley & Sons Ltd, 1989.

8. McDicken W.N., Diagnostic "Ultrasonics: Principles and Use of Instruments", Longman Group UK limited, 1991.

9. Babikian Viken L., Wechsler Lawrence R., "Transcranial Doppler Ultrasonography", Mosby-Year Book Inc., 1993.

10. Color duplex sonography: principles and clinical applications. Edited by Karl-Jurgen and Franz Fobbe. Thieme Medical Publishers, N. Y., 1995

11. Bjorn A. J. Angelsen, "Ultrasound Imaging Waves, Signals, and Signal Processing", Emantec, Norway, 2000

12. A. Boos, M. Schraag, "Method and apparatus for calculating the fetal heart", US Patent No 5,170,791,1992.

13. J. L. Muller, R. G. Hohlfeld, "System for measuring the period of a quasi-periodic signal", US Patent No 5,584,295, 1996.

14. S. L. Leon, B. S. Hopkins, "Fetal heart rate monitoring apparatusl", US Patent No 3,934,577, 1976.

15. R.J. Dunki-Jacobs and others, "Color flow imaging system utilizing a frequency domain adaptive wall filter", US Patent No 5,349,525, 1994.

16. P. J. Phillips, "Method and apparatus for medical ultrasound flash suppression", US Patent No 6,287,258, 2001.

17. C. Kasai, K. Namekawa, A. Koyano, R. Omoto "Real-Time Two-Dimensional Blood Flow Imaging Using an Autocorrelation Technique," IEEE Trans. Sonics Ultrasonics, vol. SU-32, pp. 458-464, 1985.

18. King-Yuen Wong, Zoran Banjanin, Dong-Chyuan Liu "Method and system for selectively smoothing color flow images in an ultrasound system", US Patent No 5,860,928, 1999.

19. S. F. Forestieri, R. S. Spratt, "Temporal filtering of color Doppler signal data", US Patent No 5,357,580, 1994.

20. J. A. Smith, Matthew Ellis, R. F. Saunders, A. L. Hall, "Color adaptive averaging", US Patent No 5,467,770, 1995.

21. I. M. Guracar, "Energy weighted parameter spatial/temporal filter", US Patent No 5,860,930, 1999.

22. Juin-Jet Hwang, L. S. Pflugrath "Ultrasonic diagnostic image flash suppression technique", US Patent No 5,782,769, 1998.

23. Moo Ho Bae, Ki Jong Lee, "Color imaging display method and apparatus for color flow imaging system", US Patent No 6,402,694, 2002.

24. R. N. Phelps, Z. B. Banjanin, Jin Kim, "Acoustic scan conversion method and apparatus for velocity flow", US Patent No 5,429,137, 1995.

25. Z. Z. Stroll et all, VLSI for image rotation, VLSI Signal Processing, IEEE Press, 1984.

26. Дж. Офир, Н.Ф. Маклад. Цифровые преобразователи УЗ- изображения для медицинской диагностики. ТИИЭР, т. 67, № 4, 1979, стр. 251-263.

27. D. Е. Robinson and Р. С. Knight, Interpolation scan conversion in pulse-echo ultrasound, Ultrasound Imaging, no. 4, pp. 297-310, 1982.

28. M. H. Lee, A study on performance improvements in the ultrasound B-scan system, Phys. D. Thesis, Seoul, Korea, 1985.

29. S. С. Leavitt et all, A scan conversion algorithm for displaying ultrasound images, HP Journal, vol. 34, no. 10, pp. 30-34,1983.

30. Нагулин H.E. Отечественные ультразвуковые медицинские приборы// Военно- медицинский журнал, т. CCCXXI, № 3, март, 2000, стр. 41-43.

31. О R. Andraka " A survey of CORDIC algorithms for FPGA based computers", in Proceedings of the 1998 ACM/SIGDA sixth international symposium on Field programmable gate arrays (FPGA '98), Feb. 22-24, 1998, Monterey, CA. pp. 191-200.

32. J. Vails, "Evaluation of CORDIC algorithms for FPGA design", Journal of VLSI signal processing, 32, pp. 207-222, 2002.

33. M. O'Donnell, W.E. Engeler, J.J. Bloomer, J.T. Pedicone, "Method and apparatus for digital phased array imaging", US Patent No 4,983,970, 1991.

34. M. O'Donnell, W.E. Engeler, J.J. Bloomer, J.T. Pedicone, "Dynamic phase focus for coherent imaging beam formation", US Patent No 5,111,695, 1992.

35. J.N. Wright, C.R. Cole, A. Gee, "Method and apparatus for receive beamformer system", US Patent No 5,685,308, 1997.

36. R. Schwann, G.Kappen, "CORDIC based postprocessing of ultrasound beamformer data", Proc. of the 9th International Student Conference on Electrical Engeineering "POSTER 2005", Prague, May, 26,2005.

37. Banerjee, A.S. Dhar and Swapna Banerjee, "FPGA realization of a CORDIC Based FFT processor for biomedical signal processing", Microprocessor & Microsystems, vol. 25/3, pp. 131-142, May, 2001.

38. Рабинер JI., Гоулд Б. Теория и применение цифровой обработки сигналов / Пер. с англ.; Под ред. Ю. И. Александрова. М.: Мир, 1978.

39. Оппенгейм А., Шафер Р. Цифровая обработка сигналов . М.: Связь, 1979.

40. Введение в цифровую фильтрацию. / Под ред. Р. Богнера и А. Константинидаса. -М.: Мир, 1976.

41. Хемминг Р. Цифровые фильтры. М.: Сов. Радио, 1980 г.

42. Гольденберг JI.M., Матюшкин Б.Д., Поляк М.Н. Цифровая обработка сигналов: Справочник. М.: Радио и связь, 1985 г.

43. Голд Б., Рэйдер Ч. Цифровая обработка сигналов: Пер. с англ. / Под ред. A.M. Трахтмана. -М.: Сов. радио, 1973.

44. Лэм Г., Аналоговые и цифровые фильтры: Расчет и реализация. -М.:Мир, 1982.

45. S. Bjserum, Н. Torp, and К. Kristoffersen, "Clutter filter design for ultrasound color flow imaging", IEEE Trans. Ultrason., Ferroelect., Freq. Contr., vol. 49, no. 2, pp. 204-216,2002.

46. Yoichi Ogasawara, Yasuhiko Abe, Ryoichi Kanda, "Ultrasound diagnostic apparatus", US Patent No 5,913,824, 1999.

47. C. Kargel, G. Hobenreich, B. Trummer, M. F. Insana, "Adaptive clutter rejection filtering in ultrasound strain-flow imaging", IEEE Trans. Ultrason., Ferroelect., Freq. Contr., vol. 50, no. 7, pp. 824 835, 2003.

48. E. S. Chornoboy, "Initialization for improved IIR filter performance", IEEE Trans. Signal Processing, vol. 40, no. 3, pp. 543 550, 1992.

49. A. P. Kadi and T. Loupas, "On the performance of regression and step-initialized IIR clutter filters for color Doppler systems in diagnostic medical ultrasound", IEEE Trans. Ultrason., Ferroelect., Freq. Contr., vol. 42, no. 5, pp. 927-937,1995.

50. S. M. Torres, D. S. Zrnic, "Ground Clutter Canceling with a Regression Filter", Journal of Atmospheric and Oceanic Technology, vol. 16, no. 10, pp. 1364-1372,1999.

51. D.E. Kruse and K.W. Ferrara, "A new high resolution color flow system using an eigendecomposition-based. adaptive filter for clutter rejection", IEEE Trans. Ultrason. Ferroelect. Freq. Contr., vol. 49, no. 5, pp. 1384 -1399, 2002.

52. H. Torp, "Clutter rejection filters in color flow imaging: a theoretical approach", IEEE Trans. Ultrason., Ferroelect., Freq. Contr., vol. 44, no. 2, pp. 417 -424, 1997.

53. A. Heimdal, and H. Torp, "Ultrasound Doppler measurements of low velocity blood flow: limitations due to clutter signals from vibrating muscles", IEEE Trans. Ultrason., Ferroelect., Freq. Contr., vol. 44, no. 4, pp. 873 881, 1997.

54. Ahmed Elnokrashy, Abou-Bakr M. Youssef and Yasser M. Kadah, "Nonparametric clutter rejection in Doppler ultrasound using principal component analysis," Proc. SPIE Medical Imaging, San Diego, February 2003.

55. Yasser M. Kadah, "Spatio-temporal analysis of color Doppler information using independent component analysis," Proc. SPIE Medical Imaging, San Diego, February 2002.

56. A. Hyvarinen and E. Oja, "Independent component analysis: algorithms and application" Neural Networks 13. pp. 411 -430,2000.

57. C.M. Gallippi, K.R. Nightingale, and G.E. Trahey, "Blind source separation-based adaptive filtering of physiological and ARFI-induced tissue, blood, and cyst fluid motion, in-vivo", Proc. IEEE Ultrason. Symp, 2003, vol.1, pp. 841- 846.

58. C.M. Gallippi, K.R. Nightingale, and G.E. Trahey, "BSS-based filtering of physiological and ARFI-induced tissue and blood motion", Ultrasound in Med. & Biol., vol. 29, no. 11, pp. 1583 1592,2003.

59. S. Bjserum, H. Torp, "Statistical evaluation of clutter filters in color flow imaging", Ultrasonics, no. 38, pp. 376 380,2000.

60. Lindsay I. Smith. A Tutorial on Principle Component Analysis. Personal Publication (Available at http://csnet.otago.ac.nz/cosc45362. /student tutorials/ principal components.pdf), February, 2002.

61. M.R. Sturgill, "Maximum entropy velocity estimator for ultrasonic flow imaging system", US Patent No 5,107,841.

62. R.B. Peterson, "Ultrasonic color flow imaging using autoregressive processing", US Patent No 5,226,420.

63. R.A. Hager, D.W. Clark, K.E. Thiele, J.F. Witt, "Ultrasound clutter filter", US Patent No 6,689,064.

64. J.P. Burg, "Maximum Entropy Spectral Analysis", 37th Annual Meeting, Soc. Explor. Geophys., Oklahoma, 1967.

65. Kay, Steven M., "Modern Spectral Estimation Theory and Application", Prentice Hall, 1988.

66. Z.B. Banjanin, J. Kim, C. Zobkiw, S.R. Marquis, "Ultrasonic Doppler imager having a spatially smoothed control signal for adaptive tissue rejection filter", US Patent No 5,494,037.

67. T.J, Ulrych, T.N. Bishop, "Maximum entropy spectral analysis and autoregressive decomposition", Rev. Geophys., vol. 13, No 1, pp. 183 200, 1975.

68. L.S. Griffits, "High resolution spectral estimates obtained using data extrapolation", Proc. IEEE Int. Conf. ASSP-8, pp. 584 587, 1980.

69. S.J. Johnsen, N. Andersen, "On power estimation in maximum entropy spectral analysis", Gephys., vol. 43, No 4, pp. 584 587, 1978.

70. N. Andersen, "On the calculation of filter coefficients for maximum entropy spectral analysis", Gephys., vol. 39, No 1, pp. 69 72, 1974.

71. Лутра A.K., Стайнберг Т. Изображение целей на основе метода максимальной энтропии, ТИИЭР, т. 70, JMb 1, стр. 115- 116, 1982.

72. W.D. White, "Angular spectra in radar applications", IEEE Trans., vol. AES-15, pp. 895-899, 1979.

73. T. Thorvaldsen, " Maximum entropy spectral analysis in antenna spatial filtering", IEEE Trans., vol. AP-28, pp. 556 560, 1980.

74. Гейбриел У.Ф. Спектральный анализ и методы сверхразрешения с использованием адаптивных решеток, ТИИЭР, т. 68, № 6, стр. 19-32, 1980.

75. А.К. Jain, S. Ranganath, "Extrapolation algorithms for discrete signals with application in spectral estimation", IEEE Trans. Acoust., Speech, Signal Proc., vol. ASSP-29, No 4, pp. 830-845, 1981.

76. H. Ocura, Y. Yosuda, "Spectral analysis and subtraction of noise in radar signals", IEEE Trans. Aerospace., Electron. Systems, vol. AES-17, No 1, pp. 62 -71,1981.

77. J.E. Shore, "Maximum cross-entropy spectral analysis", IEEE Trans. Acoust., Speech, Signal Proc., vol. ASSP-29, No 2, pp. 230-236, 1981.

78. E. K. Hung, R.W. Herring, "Simulation experiments to compare the signal detection properties of DFT and MEM spectra", IEEE Trans. Acoust., Speech, Signal Proc., vol. ASSP-29, No 5, pp. 1084- 1089, 1981.

79. T. Dyson, S.S. Rao, "Equal observation interval comparison of maximum entropy and weighted overlapped segment averaging spectrum estimation techniques", IEEE Trans. Acoust., Speech, Signal Proc., vol. ASSP-29, No 4, pp. 919-922,1981.

80. S.T. Alexander, E.H. Satorius, J.R. Zeidler, "Linear prediction and maximum entropy spectral analysis of finite bandwidth signals in noise", Proc. IEEE Int. Conf. ASSP-6, pp. 188 191,1978.

81. L. Marple, "High resolution autoregressive spectrum analysis using noise power cancellation", Proc. IEEE Int. Conf. ASSP-6, pp. 345 348, 1978.

82. S. Kay, "Improvement of autoregressive spectral estimates in presence of noise", Proc. IEEE Int. Conf. ASSP-6, pp. 357 360, 1978.

83. J.E. Voider, "The CORDIC trigonometric computing technique", IRE Trans. Electron. Comput., vol.EC- 3, pp. 330-334, Sept. 1959.

84. J.S. Walther "A unified algorithm for elementary functions", 1971 Proc, Joint Spring Comput. Conf., pp. 379-385, 1971.

85. Байков В.Д., Смолов В.Б. Специализированные процессоры; итерационные алгоритмы и структуры. М.: Радио и связь, 1985.

86. СБИС для распознавания образов и обработки изображений: Пер. с англ./Под. ред. К. Фу. -М: Мир. 1988.

87. М. Lohning, Т. Hentschel and G. Fettweis, "Digital down conversion in software radio terminals", Test and Control Technology, 40(5): 191 -194,2003.

88. Jarkko Vuori and Jorma Skytta, "Implementation of a digital phase-locked loop using CORDIC algorithm," in Proc. IEEE ISCAS. '96, pp. 164-167, 1996.

89. E. Antelo, J.D. Bruguera, T. Lang and E.L. Zapata, "Error analysis and reduction for angle calculation using the CORDIC algorithm", ", IEEE Trans. Comput., vol. 46, no. 11, pp. 1264- 1271, 1997

90. A.M. Despain, "Fourier transform computers using CORDIC iterations", IEEE Trans. Comput., vol. C-23, no. 10, pp. 993- 1001, 1974.

91. Нагулин H.E., Скосырев C.B., Шарапов A.A. Система формирования УЗ изображения в реальном масштабе времени// Доклады 3-ей Международной конференции «Радиоэлектроника в медицинской диагностике», 29 сентября -1 октября 1999 г., Москва, стр. 121 122.

92. Digital Receive Signal Processor AD6620 Data Sheet. Analog Devices, Inc., 1998.

93. M.B. Медведев, E.B. Юдина. Задержка внутриутробного развития плода. -М.: РАВУЗДПГ, 1998.