Сегментация изображений гистологических объектов

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

Актуальность темы диссертации

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

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

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

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

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

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

Связь работы с крупными научными программами, темами

Диссертационная работа выполнена в лаборатории обработки и распознавания изображений Института технической кибернетики НАН Беларуси в сотрудничестве с лабораторией методов морфологических исследований МГМИ в рамках тем: Интеллект 17 “Создание теоретических основ разработки, применения и распознавания изображений с использованием искусственного интеллекта” – постановление президиума НАН Беларуси № 88 от 23.11.95; “Разработать рабочий проект АРМ клинического морфолога и сдать в опытную эксплуатацию в МГМИ” в соответствии с “Программой информатизации Республики Беларусь на 1991-1995 годы и на период до 2000 года”, утвержденной Постановлением Совета Министров Республики Беларусь N 444 от 27.11.92; “Разработать платформонезависимые инструментальные средства обеспечения удалённого сбора, хранения и визуализации биомедицинской информации в НИКС” в соответствии с Постановлением Совета Министров Республики Беларусь от 22 октября 1998 г. №1609 “О развитии в республике работ по созданию единой научно-информационной компьютерной сети”.

Цель и задачи исследования

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

Для достижения поставленной цели потребовалось:

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

- разработать алгоритмы сегментации волокон и сосудов;

- разработать алгоритмы сегментации клеток;

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

- программно реализовать алгоритмы сегментации гистологических объектов и проверить их на конкретных примерах.

Объект и предмет исследования

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

Гипотеза

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

Методология и методы проведённого исследования

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

Научная новизна полученных результатов

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

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

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

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

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

Экономическая и практическая значимость полученных результатов

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

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

Алгоритмы и программный комплекс анализа и обработки медицинских изображений внедрены и используются в учебных, научных и диагностических процессах двух медицинских вузов (Минском государственном медицинском институте, Гродненском государственном медицинском институте) и трёх НИИ (Бел НИИ эпидемиологии и микробиологии (г.Минск), Бел НИИ онкологии медицинской радиологии (г.Минск), Институте биохимии НАН РБ (г.Гродно)) Республики Беларусь, а также в двух медицинских вузах (Саратовской государственной медицинской академии, Ярославльской государственной медицинской академии), Институте физиологии детей и подростков РАМН (г.Москва) и Всероссийском центре пластической хирургии глаза МПЗП РФ (г.Уфа) Российской Федерации.

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

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

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

- алгоритм полутонового утоньшения объектов на слабоконтрастных изображениях гистологических объектов;

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

- алгоритм морфологической сегментации отдельных клеток;

- алгоритм сегментации клеток со сложным фоном, основанный на объединении областей;

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

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

Личный вклад соискателя

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

Апробация результатов диссертации

Основные результаты работы докладывались и обсуждались на конференциях и симпозиумах: научно-технической конференции по компьютерной графике и анимации (Минск, 1993), 8^th International Symposium on Diagnostic Quantitative Pathology. (Amsterdam, The Netherlands, 1994), Третей научной конференции по распознаванию и анализу изображений (Минск, 1995), I конгрессе Международной ассоциации патологоанатомов (г.Москва, 1995), республиканской научной конференции молодых ученых и студентов “Актуальные проблемы современной медицины” (Минск, 1997), System and signals in Intelligent Technologies (Минск, 1998), V международной конференции “Компьютерный анализ данных и моделирование” (8-12 июня 1998, Минск), Fifth International Conference Pattern Recognition and Information Processing PRIP99 (Минск, 18 - 20 мая 1999).

Опубликованность результатов

По материалам проведённых исследований опубликовано 12 научных работ, в том числе:

1 статья в научном журнале;

3 статьи в сборниках научных трудов;

8 докладов на международных конференциях.

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

Диссертация состоит из введения, четырёх глав, заключения, приложения и списка литературы из 152 наименований. Диссертация изложена на 121 странице, включая 72 иллюстрации.

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

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

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

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

2. Автоматический анализ изображений гистологических объектов слабо освещён в литературе.

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

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

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

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

Вторая глава посвящена сегментации гистологических объектов на полутоновых изображениях.

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

Таблица 1

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

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

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

1) p2ANDp6³x AND (p1ANDp4³x ORp3ANDp0³x ORp0³x ANDp4³x);

2) p6AND p2³x AND (p5ANDp0³x ORp7ANDp4³x ORp4³x ANDp0³x);

3) p4AND p0³x AND (p3AND p6³x ORp5AND p2³x ORp2³x AND p6³x);

4) p0AND p4³x AND (p7AND p2³x ORp1AND p6³x ORp6³x ANDp2³x),

где номера пикселей соответствуют рис. 1.

Рис. 1 Окрестность пикселя X

В работе предлагаются три варианта изменения пикселя, удовлетворяющего вышеуказанным условиям: