РАДЮЕЛЕКТРОН1КА ТА ТЕЛЕКОМУН1КАЦ11
УДК 004.93
Д. М. Пиза, А. А. Степаненко
СЕГМЕНТАЦИЯ УЛЬТРАЗВУКОВЫХ МЕДИЦИНСКИХ ИЗОБРАЖЕНИЙ НА ОСНОВЕ МЕТОДА НЕЧЕТКОЙ
КЛАСТЕРИЗАЦИИ
Рассмотрен новый подход к повышению эффективности процедуры визуального анализа ультразвуковых медицинских эхо-импульсных изображений. Представлены результаты экспериментальной проверки информационных возможностей метода.
ВВЕДЕНИЕ
Ультразвуковая интроскопия относится к одному из основных и распространенных методов медицинской диагностики [1], в рамках которой основное внимание уделяется повышению чувствительности и достовер-ности процедуры визуализации ультразвуковых эхо-импульсных изображений. В работе [2] рассматривались возможности повышения чувствительности выделения границ разделов низкоконтрастных структур и нейтрализации влияния многократных переотражений ультразвуковых импульсов, затрудняющих визуальный анализ эхо-импульсных изображений. Однако, с точки зрения определения однородности, шероховатости, регулярности, наибольший интерес представляет процедура сегментации непосредственно самого ультразвукового изображения. Вопросы сегментации различных типов изображений относятся к области задач компьютерного зрения [3], однако характерной особенностью ультразвуковой диагностики является наличие модуляционных спекл-шумов, поэтому применение любого метода сегментации ультразвуковых изображений должно учитывать наличие этого фактора и сопровождаться использованием процедуры фильтрации модуляционных шумов.
Еще одной характерной особенностью ультразвуковой диагностики является то обстоятельство, что на ультразвуковых изображениях границы анатомических областей являются «размытыми», поскольку акустические импедансы соседних биоструктур могут быть близкими либо плавно изменяющимися. Это обстоятельство приводит к необходимости использования нечетких методов кластеризации. Нечеткая кластеризация вводит понятие нечетких кластеров и функцию принадлежности объектов к ним, изменяющуюся в интервале [0,1], что позволяет оценить степень принадлежности объекта к тому либо иному классу. Основой нечеткой кластеризации служит метод FCM (Fuzzy C-Means), имеющий множество модификаций [4]. Метод
© Пиза Д. М., Степаненко А. А., 2008
24
ЕСМ базируется на использовании идей и математического аппарата нечеткой логики и широко используется применительно к задачам анализа медицинских радиологических изображений [5], однако информационные возможности метода применительно к анализу ультразвуковых изображений ранее не исследовались.
1 МЕТОД ИССЛЕДОВАНИЯ
В результате работы алгоритма ЕСМ каждому пикселю изображения ставится в соответствие вектор из функций принадлежности к каждому классу, на основе которого можно делать выводы о природе данного объекта, но лишь при нейтрализации влияния шумових факторов. В настоящее время отсутствует единый поход к решению этой задачи. В [6] предлагалось использовать винеровскую фильтрацию, а в [7] адаптивную анизотропную фильтрацию. Основной проблемой здесь является то обстоятельство, что для осуществления оптимальной фильтрации необходимо априорное знание спектральных характеристик как шума, так и самого незашумленного (т. е. идеального) изображения. В этой связи интерес представляет метод селективного сингулярного разложения автоморфного отображения (МССРАО) [8] зашумленных изображений, не требующий априорного знания спектральных характеристик, а лишь приблизительного знания интервала корреляции шума.
МССРАО относится к рамочным методам анализа изображений, причем размер квадратной рамки должен соответствовать интервалу корреляции. На основе МССРАО формируется многомерный массив, к которому применяется процедура декоррелирующего преобразования с формированием нового массива «собственных» изображений. В качестве отфильтрованного изображения принимается первое «собственное» изображение, соответствующее наибольшему сингулярному числу.
Приведенные в [8] экспериментальные результаты относились к анализу изображений геофизических полей и информационные возможности метода применительно к анализу ультразвуковых эхо-импульсных изображений остались невыясненными.
ISSN 1607-3274 «Радюелектрошка. 1нформатика. Управлшня» № 1, 2008
Д. М. Пиза, А. А. Степаненко: СЕГМЕНТАЦИЯ УЛЬТРАЗВУКОВЫХ МЕДИЦИНСКИХ ИЗОБРАЖЕНИЙ НА ОСНОВЕ МЕТОДА НЕЧЕТКОЙ КЛАСТЕРИЗАЦИИ
Поэтому для дальнейшего повышения эффективности медицинской диагностики необходимо исследование информационных возможностей МССРАО и нечеткой БСМ кластеризации применительно к задачам сегментации медицинских ультразвуковых эхо-импульсных изображений.
2 РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЙ
Нечеткая кластеризация является обобщенным вариантом традиционного алгоритма К-средних, преодолевающим ограничения принадлежности выборок только к одной группе. В нечеткой кластеризации используется нечеткое множество для задания степени принадлежности выборок только к одной группе. Численное значение степени принадлежности лежит в диапазоне [0-1], так что сумма всех значений принадлежности одного элемента (пикселя) ко всем группам равна единице.
Алгоритм БСМ кластеризации базируется на использовании процедуры итеративной минимизации целевой функции вида
N
J(и, V) = £ £
г = 1 п = 1
т |
и{п\Хп '
(1)
где V = {»1, ..., »с} - центры кластеров; и = [и1п] -матрица размером с х N, где щп есть г-я функция принадлежности п-го входа хп, т е [ 1,<») - параметр фаззификации, с - задаваемое число кластеров. Значения функций принадлежности должны удовлетворять следующим ограничениям
0 < игп < 1, г = 1, 2,..., с, п = 1, 2,..., N; (2)
£ игп = 1 п = 1 2,..„N;
г = 1
N
0 < £ иш < N, г = 1, 2,..., с.
п = 1
(3)
(4)
Алгоритм является итеративным и основан на последовательном вычислении следующих уравнений:
N
N
п = 1
£ игпхгп, г = 1 2,., с;
(5)
п = 1
1 /(т - 1)
£
] = 1
|хп »]\
1 /(т - 1)'
г = 1, 2, с, п = 1, 2, ..., N.
(6)
Для каждого пикселя ансамбля изображений, значение его функции принадлежности для каждого кластера зависит от расстояния до соответствующего центра кластера. Параметр фаззификации т уменьшает влияние небольших значений функции принадлежности. К достоинствам алгоритма нечетких С-средних можно отнести его гибкость, позволяющую работать со значениями функций принадлежности. К недостаткам можно отнести необходимость задания априорного числа кластеров и теоретическая неопределенность выбора начальных значений и параметра фаззификации, которая приводит к получению различных результатов кластеризации, а также чувствительность к шуму.
На рис. 1, а-г в координатах «время - пространство» представлено ультразвуковое эхо-импульсное изображение головного мозга ребенка (а) и первые три «собственных» изображения МССРАО (б-г). Такие эхо-импульсные изображения рассматриваються, как правило, в виде набора сигнальных трасс (столбцов изображения). Из рассмотрения исходного изображения на рис. 1, а видно, что оно является зашумлен-ным, тогда как на первом «собственном» изображении МССРАО (рис. 1, б) влияние модуляционной помехи фактически нейтрализовано. Рассматривая амплитудно-временные срезы (рис. 2, а и б) 50-й импульсной трассы (50-го столбца изображений рис. 1, а и б), заметно очищение сигнала (рис. 2, б) от излишних шумов (быстроизменяющейся амплитуды в пиках и спадах (рис. 2, а)).
На рис. 2, в, г показаны амплитудно-временные срезы второго и третьего «собственных» изображений МССРАО (рис. 1, в, г), соответствующих шумовым компонентам, и их отбрасывание фактически эквивалентно эффекту адаптивной фильтрации, поскольку никаких предположений о спектральных и статистических характеристиках модуляционных шумов не использовалось.
На рис. 3 представлены результаты применения метода нечеткой кластеризации для сегментации рис. 1, а и рис. 1, б на семь классов.
Из сопоставления рис. 1, а с рис. 3, а следует, что применение метода нечеткой кластеризации для сегментации зашумленного изображения фактически ничего не дает, тогда как его применение к отфильтрованному изображению (рис. 3, б) позволяет облегчить процедуру визуального анализа и интерпретации
и=
2
V:
V =
т
и
РАДЮЕЛЕКТРОН1КА ТА ТЕЛЕКОМУН1КАЦ11
а) б) в) г)
Рисунок 1 - Ультразвуковое изображение головного мозга:
а - оригинал; б-г - первые три «собственных» изображения, синтезированных на основе использования МССРАО (горизонтальная ось - пространственная; вертикальная ось - временная)
Рисунок 2 - Амплитудно-временные срезы 50-й сигнальной трассы изображений на рис. 1, а-г соответственно (горизонтальная ось - временная, вертикальная - относительная амплитуда)
26
1607-3274 «Радюелектрошка. 1нформатика. Управлшня» № 1, 2008
И. М. Прохорец: МАТЕМАТИЧЕСКИЕ МОДЕЛИ ДЛЯ АНАЛИЗА СТАТИЧЕСКИХ ХАРАКТЕРИСТИК МНОГОЭЛЕМЕНТНЫХ ПОЛУПРОВОДНИКОВЫХ ДЕТЕКТОРОВ
а) б)
Рисунок 3 - Результат применения метода нечеткой кластеризации: а - рис. 1, а; б - рис. 1, б
результатов диагностики посредством, например, варьирования числа классов.
ВЫВОДЫ
Нечеткая кластеризация является достаточно простым и удобным методом повышения эффективности процедуры визуального анализа ультразвуковых медицинских эхо-импульсных изображений. Необходимым условием успешного применения метода нечеткой кластеризации является нейтрализация влияния модуляци-
онных шумов. Метод селективного сингулярного разложения автоморфного отображения зашумленных изображений является эффективным способом фильтрации модуляционных шумов, поскольку не требует априорного знания их спектральных и статистических характеристик.
ПЕРЕЧЕНЬ ССЫЛОК
1. Физика визуализации изображений в медицине. Под ред. С. Уэбба. - М.: Мир. - 1991. - Т. 2. - 406 с.
2. Ахметшин А. М., Степаненко А. А. Повышение чувствительности ультразвуковой диагностики на основе метода параметрического спектрального анализа второго порядка. // Харьков: Клиническая информатика и телемедицина. - 2005. - Т. 2, № 1. - С. 98-100.
3. Форсайт Д., Понс Ж. Компьютерное зрение: современный подход. - М.: Вильямс. - 2004. - 926 с.
4. Леоненков А. Нечеткое моделирование в среде MAT-LAB и fuzzyTECH. - Санкт-Петербург: БХВ-Петербург. -2003. - 719 с.
5. Hall L. O., Bensaid A., Clarke L. P. A comparison of neural network and fuzzy clustering techniques in segmentation MRI of the brain // IEEE Trans. On Neural Networks. -1992. - Vol. 3. - P. 672-682.
6. Применение ультразвука в медицине. Под ред. К. Хил-ла. - М.: Мир. - 1989. - 567 с.
7. Karaman M., KutayA., Bozdagi G. An adaptive speckle suppression filter for medical ultrasonic imaging // IEEE Transaction on Medical Imaging. - 1995. - Vol. 14, No. 2. -P. 283-292.
8. Ахметшина Л. Г. Адаптивная фильтрация шумов в сигналах и изображениях: метод селективного сингулярного разложения автоморфного отображения // Искусственный интеллект. - 2005. - № 3. - С. 328-335.
Надшшла 24.03.2008
Розглянуто новий nidxid до тдвищення ефективност1 процедуры вiзуaльного aнaлiзу ультразвукових медичних ехоАмпульсних зображень. Представлено результати екс-периментальноЧ nеревiрки тформацшних можливостей методу.
A new approach to increasing efficiency of visual analysis ultrasonic medical pulse-echo image is considered. Results experimental investigation information possibilities of the new method are presented.
УДК 621.382.002
И. М. Прохорец
МАТЕМАТИЧЕСКИЕ МОДЕЛИ ДЛЯ АНАЛИЗА СТАТИЧЕСКИХ ХАРАКТЕРИСТИК МНОГОЭЛЕМЕНТНЫХ ПОЛУПРОВОДНИКОВЫХ
ДЕТЕКТОРОВ
На основании теории направленных графов рассмотрены математические модели стрип-детектора для измерения межстриповых сопротивлений на постоянном токе. Приведены примеры расчета совместных линейных уравнений для электрических цепей, используемых для измерения межстриповых сопротивлений.
© Прохорец И. М., 2008
ВВЕДЕНИЕ
Многоэлементные полупроводниковые детекторы, к которым относятся полосковые или стрип-детекторы, пиксельные детекторы, относятся к приборам, получив-