главнаядиссертациябиблиотекассылкипоиск
Белявцев Алексей Анатольевич
E-mail :
allweb@mail.ru
Тема магистерской работы :
«Разработка специализированной компьютерной системы диагностики клеток на основе анализа изображений»
Руководитель : проф. Скобцов Ю. А.

Автореферат магистерской работы

Актуальность задачи выявления опухолевых клеток определяется возрастанием процента раковых заболеваний, которые являются второй по частоте причиной смертности в развитых странах. Многочисленные медико-социологические исследованиями показывают, что примерно 33% людей либо имеют онкологические заболевания, либо с определенной вероятностью могут ими заболеть. По статистике 40% больных раком при различном лечении выживают на протяжении не менее 5 лет. Для диагностирования и прогнозирования прогрессирования опухолей обязательным является гистологическое исследование атипичных клеточных структур и определение гистологической принадлежности опухоли. Использование математических методов и алгоритмов обработки гистологических изображений дает возможность, как ускорить процесс анализа препаратов, так и повысить точность результатов исследования. Сложность получения качественных гистологических препаратов существенно снижает эффективность компьютерного анализа.

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

  • Фон препарата, наличие и характер межуточного вещества;
  • Количество и расположение клеток, образование комплексов или структур, сохранность клеточных границ;
  • Размеры и формы клеток;
  • Ядро - форма и размеры, расположение и окрашиваемость;
  • Ядерно-цитоплазматическое соотношение;
  • Характер строения хроматина;
  • Ядрышки - количество, форма, размеры, четкость границ;
  • Цитоплазма - объем, окраска, четкость границ, секреция включения, вакуолизация;
  • Наличие многоядерных клеток, фигур деления (атипичные митозы).

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

Клетка

  • Размер клеток превосходит размеры клеток той ткани, которая явилась источником опухолевого роста: клетки могут быть гигантских размеров. Изменение размера не абсолютный показатель - размеры клеток могут быть не изменены;
  • Изменение формы клеток - не полностью или мало соответствует той форме, которая свойственная клеткам нормальной ткани. Форма может быть самая причудливая;
  • Может определяться отчетливый клеточный полиморфизм, т. е. клетки разных размеров и формы;
  • Изменяется ядерно-цитоплазматическое соотношение и как правило в пользу ядра (за счет увеличения последнего);
  • Появление многоядерных клеток, с отчетливым ядерным полиморфизмом;
  • Атипичное расположение ядра;
  • Диссоциация в созревании ядра и цитоплазма (молодое ядро в зрелой клетке, зрелое ядро в незрелой цитоплазме).

Ядро

  • Увеличение размера ядра;
  • Изменение формы ядра, иногда ядра приобретают причудливую форму, ядерный полиморфизм;
  • Контур ядра неправильный, неравномерно извилистый, иногда с глубокими вырезками, нередко грубо очерчен;
  • Ядерная мембрана неравномерно утолщена, разрывы.

Строение хроматина

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

Ядрышки

  • чаще всего определяются;
  • размеры увеличены;
  • форма неправильная, нередко полиморфизм;
  • число увеличено.

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

  • разработать алгоритмы сегментации и выделения клеток;
  • классифицировать изображения гистологических объектов;
  • программно реализовать алгоритмы сегментации гистологических объектов и проверить их на конкретных примерах.

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

Методы цифровой обработки изображения

Предварительный анализ изображений позволяет сделать вывод о том, что:

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

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

Классификация основных методов обработки и распознавания изображений.

Рис.1 - Основные процедуры и методы распознавания изображений

Изменение контраста

Слабый контраст обычно вызван малым динамическим диапазоном изменения яркости, либо сильной нелинейностью в передаче уровней яркости. Простейшим методом контрастирования является функциональное отображение градации яркости. На практике очень часто используют линейные функциональные отображения. Если в результате неравномерности освещения при фотографировании или изготовлении фотографий, возникает ситуация, когда различные участки изображения обладают разным контрастом. В таком случае для изменения контраста используют адаптивные алгоритмы контрастирования. Примером может служить алгоритм локального усиления контраста. Высокая эффективность работы алгоритма достигается в том случае, если на снимке присутствуют области с явно завышенным или заниженным контрастами. Суть алгоритма состоит в том, что снимок рассматривается как набор некоторого числа локальных областей, и эти области обрабатываются с учетом их характеристик. Обработка ведется в следующей последовательности: вычисляется коэффициент усиления срезов плотности Р отдельно для каждого локального участка изображения. И осуществляется обработка каждого пикселя изображения. Если Р=1, то над локальным участком изображения никакого действия не производится (если р отлично от единицы, то осуществляется повышение контраста локальной области). Первоначально вычисляется контраст в анализируемой точке относительно ближайшей окрестности. Затем значение относительного контраста складывается с единицей, и полученное значение принимается в алгоритме как коэффициент усиления Р, а далее производится вычисление по формуле:

формула

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

Сглаживание шумов

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

Подчеркивание границ

Методы сглаживания изображений могут устранять шум очень эффективно. Существенным недостатком алгоритмов сглаживания является смаз изображения (т.е. снижение четкости контурных элементов), при этом величина смаза пропорциональна размеру маски, используемой для сглаживания. Для однозначного анализа изображений, особенно при вычислении геометрических характеристик структурных элементов, очень важно убрать смаз с контуров объектов в изображении, то есть усилить разницу между градациями яркости контурных элементов объекта и соседних элементов фона. В этом случае при обработке изображений используются методы подчеркивания контуров. Обычно подчеркивание границ осуществляется методом высокочастотной пространственной фильтрации. Характеристики фильтров задаются в виде маски, в которой среднее значение должно быть равно нулю. Еще одним методом подчеркивания границ является так называемое статическое дифференцирование. В этом методе значение яркости каждого элемента делится на статистическую оценку среднеквадратического.

Медианная фильтрация

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

Сегментация изображений

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

В общем случае сегментация представляет собой операцию разбиения конечного множества плоскости, на которой определена функция исходного изображения f(x,y) на k непустых связанных подмножеств si (i=[1,k]) в соответствии с некоторым предикатом P, определяемом на множестве S={s1,s2,…,sk} и принимающий истинные значения, когда любая пара точек из каждого подмножества si удовлетворяет некоторому критерию однородности (например, критерий однородности, основанный на оценке максимальной разности яркости отдельного пикселя и среднего значения яркости, вычисленного по соответствующей области).

Методы сегментации можно разделить на следующие основные группы:

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

Таблица 1 - Классификация гистологических объектов и методов их сегментации

Вид объекта

Хар-ка изображения

Равномерный фон

Неравномерный фон

Площадные

Отдельно лежащие объекты одного типа

Пороговая сегментация

Методы математической морфологии

Объекты

Объекты, сопровождаемые объектами другого типа

Пороговая сегментация

Методы объединения областей

Протяженные объекты

Произвольное изображение

Пороговая сегментация

Морфологическая сегментация, основанная на утоньшении

Мелкие контрастные объекты

Произвольное изображение

Пороговая сегментация

Пороговая сегментация

Методы выделения контуров

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

Методы выделения контуров на изображении можно разделить на следующие основные классы:

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

Общим для всех этих методов является стремление рассматривать границы как область резкого перепада функции яркости изображения f(i,j); отличает же их вводимая математическая модель понятия границы и алгоритм поиска граничных точек. В соответствии с поставленными задачами к алгоритмам выделения контуров предъявляются следующие требования: выделенные контура должны быть утоньщенными, без разрывов и замкнутыми. Таким образом, процесс выделения контуров несколько усложняется в связи необходимостью применять алгоритмы утоньшения и устранения разрывов. Однако и это не всегда дает желаемого результата - в большинстве случаев контуры получаются незамкнутыми и, как следствие, непригодными для ряда процедур анализа.

Анимация. Приведение изображения к бинарному и определение его контура.

Вычисление морфометрических признаков

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

Пусть - количество пикселов, относящихся к объекту.

Все множество пикселов p(x,y), относящихся к объекту, обозначим . Тогда координаты центра масс объекта вычисляются как

Вычисляем несколько всомагательных величин:




Тогда длины максимальной и минимальной осей инерции вычисляются как:

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

Все эти параметры будут одаваться на входы нейронной сети.

Методы обучения нейронных сетей

Из существующих методов обучения нейронных сетей необходимо будет выбрать наиболее подходящие для решения поставленной задачи.

Таблица 2 - Алгоритмы обучения нейронных сетей

Парадигма

Обучающее правило

Архитектура

Алгоритм обучения

Задача

С учителем

Коррекция ошибки

Однослойный и многослойный перцептрон

Алгоритмы обучения перцептрона. Обратное распространение. Adaline и Madaline.

Классификация образов. Аппроксимация функций. Предскащание, управление.

Больцман

Рекуррентная

Алгоритм обучения Больцмана

Классификация образов

Хебб

Многослойная прямого распространения

Линейный дискриминантный анализ

Анализ данных. Классификация образов.

Соревнование

Соревнование

Векторное квантование

Категоризация внутри класса. Сжатие данных.

Сеть ART

ARTMap

Классификация образов

Без учителя

Коррекция ошибки

Многослойная прямого распространения

Проекция Саммона

Категоризация внутри класса. Анализ данных.

Хебб

Прямого распространения или соревнование

Анализ главных компонентов

Анализ данных. Сжатие данных.

Сеть Хопфилда

Обучение ассоциативной памяти

Ассоциативная память

Соревнование

Соревнование

Векторное квантование

Категоризация. Сжатие данных.

SOM Кохонена

SOM Кохонена

Категоризация. Анализ данных.

Сети ART

ART1, ART2

Категоризация

Смешанная

Коррекция ошибки и соревнование

Сеть RBF

Алгоритм обучения RBF

Классификация образов. Аппроксимация функций. Предсказание, управление.

Распознавание

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

Для реальных задач распознавания применяются, в основном, четыре подхода, использующие методы:

  • корреляционные, основанные на принятии решений по критерию близости с эталонами;
  • признаковые и синтаксические - наименее трудоемкие и нормализации, занимающие промежуточное положение по объему вычислений.

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

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

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

Главную трудность в признаковых методах составляет выбор признаков. При этом исходят из естественных правил:

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

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

Требования к системе

Система должна выполнять такие функции:

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

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

Литература

  1. Прэтт У.К. "Цифровая обработка изображений". - М.: Мир, 1982. T.I.-2. 792 с.
  2. Сойфер В.А. "Компьютерная обработка изображений". Часть 2 Методы и алгоритмы. - Соросовский образовательный журнал, №3, 1996.
  3. Ярославский Л.П. "Введение в цифровую обработку изображений". - М: Сов. радио, 1979 - 312 с., ил.
  4. Путятин Е.П., Аверин С.И. "Обработка изображений в робототехнике". - М.: Машиностроение, 1990. - 320 с.
  5. Рудаков П.И., Сафонов В.И. "Обработка сигналов и изображений". - М.: ДИАЛОГ-МИФИ, 2000. - 416 с.
  6. Андрианов Д.Е., и др. "Обработка УЗИ-снимков. Проблемы передачи и обработки информации в информационно-вычислительных сетях": Тез. докл. 7-го Международного научно-технического семинара / Под ред. И.А. Цветкова. М.: НИЦПрИС,– 1997. с. 50 – 51.

Главная :: Диссертация :: Библиотека :: Ссылки :: Поиск

Alex 2004