Другие журналы

научное издание МГТУ им. Н.Э. Баумана

НАУКА и ОБРАЗОВАНИЕ

Издатель ФГБОУ ВПО "МГТУ им. Н.Э. Баумана". Эл № ФС 77 - 48211.  ISSN 1994-0408

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

# 06, июнь 2013
DOI: 10.7463/0613.0574140
Файл статьи: Artyuhova_P.pdf (1078.38Кб)
авторы: Артюхова О. А., Самородов А. В.

УДК 57.087

Россия, МГТУ им. Н.Э. Баумана

olga.artyukhova@mail.ru

avsbmstu@yandex.ru

 

Введение

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

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

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

 

Алгоритм сегментации флуоресцентно-микроскопических изображений

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

 

Untitled_01

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

 

Для сегментации цветных изображений обычно используются пороговые методы, методы кластеризации, наращивания областей, построения физических моделей формирования изображений и другие [1, 7]. Поскольку на рассматриваемых изображениях координаты цветности клеток и фона зачастую являются идентичными, для сегментации клеток был выбран метод пороговой обработки полутонового изображения – канала V (value) цветового пространства HSV[8] (рисунок 2).

 

Untitled_02

Рисунок 2 – Анализируемое изображение в цветовом канале V

 

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

 

 

(1)

 

где  и  – соответственно среднее значение и стандартное отклонение значений яркости изображения в некоторой окрестности рассматриваемой точки;

k – коэффициент, значение которого определяется экспериментально (в данном случае при сегментации клеток k = –1).

Результат автоматического выделения масок клеток с использованием бинаризации по методу Ниблэка приведен на рисунке 3. Размер стороны квадратного окна, по которому рассчитывается локальный порог бинаризации при сегментации клеток, соответствует утроенному среднему диаметру клетки и для изображений, зарегистрированных с шагом дискретизации 0,14 мкм/пиксель, равен 700 пикселей [10, 11].

 

Untitled_03

Рисунок 3 – Результат автоматического построения масок клеток по методу Ниблэка

 

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

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

 

 

,

(2)

 

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

 и  – концентрации красителей.

В соответствии с этим модель, описывающая формирование элемента изображения, имеет вид:

 

 

(3)

 

где , ,  – значения интенсивности элемента изображения в цветовых каналах R, G и Bсоответственно;

, ,  – коэффициенты пропорциональности, соответствующие вкладу в суммарную интенсивность флуоресценции для каналов R, G и Bсоответственно, i-го красителя.

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

 

 

,

(4)

 

где     ,     ,     .

Решение уравнения (4) методом наименьших квадратов обеспечивает нахождение профилей концентраций флуоресцирующих красителей:

 

 

.

(5)

 

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

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

 

Untitled_04

Рисунок 4 – Результат цветовой деконволюции анализируемого изображения, характеризующий профиль концентрации красителя ФИТЦ

 

Как правило, при приготовлении препаратов клеточных культур не происходит полного гашения фоновой флуоресценции, вызванной присутствием флуоресцирующего красителя не только в исследуемом биологическом микрообъекте, но и в целом на всей поверхности препарата (рисунок 4). В силу этого для построения масок хламидийных включений по результату цветовой деконволюции также была использована локальная адаптивная бинаризация по методу Ниблэка (рисунок 5). Эмпирически было установлено, что наилучшие результаты бинаризации достигаются при значении коэффициента kв соотношении (1), равного двум. Размер стороны квадратного окна, по которому рассчитывается локальный порог бинаризации для сегментации хламидийных включений, при шаге дискретизации изображения 0,14 мкм/пиксель, равен 500 пикселей [10, 11].

 

Untitled_06

Рисунок 5 – Результат автоматического построения масок хламидийных включений по методу Ниблэка

 

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

 

Untitled_09

Рисунок 6 – Итоговый результат сегментации анализируемого изображения

 

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

 

Разработка быстрого алгоритма бинаризации по методу Ниблэка

При использовании алгоритма бинаризации на основе метода Ниблэка изображение сканируется квадратным окном, по которому для каждого положения рассчитывается локальное среднее значение и стандартное отклонение интенсивности элементов изображения. Значительный размер окна сканирования приводит к существенным временным затратам при прямой реализации рассматриваемого алгоритма (несколько суток для ЭВМ с процессором Intel Core i5 с частотой 2,27 ГГц, объемом оперативной памяти 4,00 ГБ и 64-разрядной операционной системой) и невозможности его использования на практике. Для ускорения работы алгоритма бинаризации, т.е. для ускорения расчета локальных средних значений и стандартных отклонений интенсивности элементов изображения, был использован подход, основанный на построении интегрального изображения.

Интегральное изображение  строится по исходному изображению  следующим образом [15]:

 

 

(6)

 

С использованием интегрального изображения локальное среднее значение интенсивности изображения внутри прямоугольной области с центром в точке  и со сторонами  и , умноженное на число элементов этой области, может быть вычислено с использованием только трех операций сложения (рисунок 7):

 

 

(7)

 

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

 – число элементов рассматриваемой области.

 

а)

б)

а) расположение сканирующего окна на изображении; б) точки интегрального изображения, использующиеся для вычисления локальных характеристик изображения

Рисунок 7 – Использование интегрального изображения для вычисления локальных характеристик изображения

 

Основные потери времени при бинаризации изображений по методу Ниблэка связаны с вычислением локального стандартного отклонения значений интенсивности элементов изображения. В данной работе для ускорения указанных вычислений было предложено использовать не только интегральное изображение (6), но и интегральный квадрат изображения , т.е. интегральное изображение, построенное для поэлементно возведенного в квадрат исходного изображения:

 

 

(8)

 

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

 

 

(9)

 

Для расчета локального порога бинаризации при прямом вычислении среднего значения и стандартного отклонения в некоторой точке изображения, лежащей в центре квадратной области со стороной d, в соответствии с (1) используется  операций сложения, три операции умножения,  операций возведения в квадрат и одна операция извлечения квадратного корня. При использовании интегральных представлений (6) и (8) для вычисления локального порога число требуемых операций не зависит от размера сканирующего окна и в пересчете на один элемент изображения снижается в среднем до двенадцати сложений, четырех умножений, двух возведений в квадрат и одного извлечения корня.

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

 

Заключение

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

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

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

 

Список литературы

1. Wu Q., Merchant F A., Castleman K.R. Microscope Image Processing. Oxford: Academic Press, 2008. 548 p.

2. Imaging cellular and molecular biological functions / Shorte S.L., Frischknecht F. (Eds.). SpringerBerlinHeidelberg, 2007. 472 p. (Ser. Principles and Practice). DOI: 10.1007/978-3-540-71331-9

3. Световаямикроскопиявбиологии. Методы: пер. с англ. / Под ред. А. Лейси. М.: Мир, 1992. 464 с.

4. Лукичева Т.И. и др. Клинико-лабораторные аналитические технологии и оборудование: учеб. пособ. для студ. средн. проф. учеб. заведений / под ред. проф. В.В. Меньшикова. М.: Издательский центр «Академия», 2007. 240 с.

5. Гранитов В. М. Хламидиозы. М.: Медицинская книга; Н. Новгород: Изд-во НГМА, 2002. 192 с.

6. Якубович А.И., Корепанов А.Р. Урогенитальный хламидиоз. Иркутск: Полиграфический центр «РИЭЛ», 2007. 108 с.

7. Skarbek W., Koschan A. Colour image segmentation: a survey. Technical report № 94-32. Technical University of Berlin, 1994. 80 p.

8. Гонсалес Р., Вудс Р., Эддинс С. Цифровая обработка изображений в среде MATLAB : пер. с англ. М.: Техносфера, 2006. 616 с.

9. Trier O.D., Taxt T. Evaluation of Binarization Methods for Document Images // IEEE Transactions on Pattern Analysis and Machine Intelligence. 1995. Vol. 17, no. 3. P. 312-315. DOI: 10.1109/34.368197

10. Артюхова О.А. Автоматизированный анализ флюоресцентных изображений микробиологических препаратов //X МНТК «Физика и радиоэлектроника в медицине и экологии ФРЭМЭ’2012»: материалы. Владимир, 2012. Кн. 2. С. 111-113.

11. Артюхова О.А., Кост Е.А., Самородов А.В., Диденко Л.В. Разработка алгоритма автоматизированного анализа изображений хламидий в культуре клеток // V ВНПК «Цитоморфометрия в медицине и биологии: фундаментальные и прикладные аспекты»: материалы. М., 2012. С. 10-11.

12. Ruifrok A.C., Johnston D.A. Quantification of histochemical staining by color deconvolution //Anal Quant Cytol Histol. 2001. Vol. 23. P. 291-299.

13. Оптическая биомедицинская диагностика: учеб. пособие для вузов / пер. с англ. под ред. В.В. Тучина. В 2 т. Т. 1. М.: Физматлит, 2007. 560 с.

14. Digital Microscopy / Sluder G., Wolf D. (eds). 2nd  ed. of "Video Microscopy". Academic Press, 2003. 544 p. (Ser. Methods in Cell Biology, vol. 72). ISBN 13: 978-0-12-564169-2

15. Viola P., Jones M.J. Robust real-time face detection // International Journal of Computer Vision. 2004. Vol. 57, no. 2. P. 137-154. DOI: 10.1023/B:VISI.0000013087.49260.fb

Поделиться:
 
ПОИСК
 
elibrary crossref ulrichsweb neicon rusycon
 
ЮБИЛЕИ
ФОТОРЕПОРТАЖИ
 
СОБЫТИЯ
 
НОВОСТНАЯ ЛЕНТА



Авторы
Пресс-релизы
Библиотека
Конференции
Выставки
О проекте
Rambler's Top100
Телефон: +7 (915) 336-07-65 (строго: среда; пятница c 11-00 до 17-00)
  RSS
© 2003-2020 «Наука и образование»
Перепечатка материалов журнала без согласования с редакцией запрещена
 Тел.: +7 (915) 336-07-65 (строго: среда; пятница c 11-00 до 17-00)