Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета, 2007, №29

Научный журнал КубГАУ, №29(5), май 2007 года 
 

http://ej.kubagro.ru/2007/05/pdf/03.pdf
 

1

УДК 681.513 
 
UDC 681.513 
 
ПРИМЕНЕНИЕ ПРИНЦИПОВ ПОСТРОЕНИЯ 
СИСТЕМ МАШИННОГО ЗРЕНИЯ В ЗАДАЧЕ 
АНАЛИЗА ИЗОБРАЖЕНИЙ КЛЕТОЧНЫХ 
СТРУКТУР 
 

APPLICATION OF PRINCIPLES OF 
CONSTRUCTION OF MACHINE VISION 
SYSTEM IN A TASK OF  CELL STRUCTURE 
IMAGE ANALYSIS 

Ковригин Алексей Владимирович 
 аспирант 
 

Kovrigin Alexey Vladimirovich 
post- graduate student 

Политехнический институт Сибирского 
федерального университета 
 

Polytechnic Institute of Siberian Federal University, 
Russia 

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

The approach on the base of principles of construction 
of machine vision system to decide the tasks of cell 
structure image analysis are considered. The advantage 
of the offered approach consists of possibilities  
to overcome the problem of increasing of error 
classification under presence of intersections  
of objects with semitransparent  structure. 

Ключевые слова: МАШИННОЕ ЗРЕНИЕ, 
КОМПЬЮТЕРНОЕ ЗРЕНИЕ, 
АВТОМАТИЧЕСКАЯ КЛАССИФИКАЦИЯ, 
СЕГМЕНТАЦИЯ, КЛАСТЕРИЗАЦИЯ 
 

Key words: MACHINE VISION, COMPUTER 
VISION, AUTOMATIC CLASSIFICATION, 
SEGMENTATION, CLUSTERING. 

 

1. Введение 

Большинство медицинских методик ручного микроскопического 

анализа эксплуатирует свойства глаз человека распознавать соотношения 

формы, оттенков, текстуры сложных изменчивых объектов без их 

измерения, на основе качественных характеристик. Если исследуемых 

клеток много, а требуемый объем выборки невелик (100-200), то 

визуальный анализ (например, подсчет лейкоформулы) при наличии 

достаточного числа опытных врачей-лаборантов может выполняться в 

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

дифференциальным подсчетом более серьезных выборок клеток (около 

500), сопровождающимся большим количеством измерений (например, 

определение 
индекса 
овалоцитоза 
эритроцитов). 
Весьма 
трудны 

психологически и физически для человека задачи поиска редких клеток, 

подсчета сложных дифференциальных формул, аккуратного просмотра 

больших пространств препарата (например, подсчет формул при 

цитопении, поиск юных клеток в мазке крови, подсчет миелограммы в 

Научный журнал КубГАУ, №29(5), май 2007 года 
 

http://ej.kubagro.ru/2007/05/pdf/03.pdf
 

2

мазке костного мозга, подсчет патологических типов эритроцитов). Очень 

трудны 
и 
даже 
вредны 
для 
глаз 
врача 
анализы 
с 
подсчетом 

флуоресцирующих меток (например, подсчет субпопуляции меченых FITC 

лимфоцитов). Далеко не вся присутствующая в цитологическом препарате 

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

В связи с этим доминировавшая до последнего времени ручная 

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

субъективна и ограничена по своим возможностям.  

Развитие алгоритмов машинного зрения позволяет постепенно 

перейти от ручного способа анализа клеточных структур. Принципы 

построения системы машинного зрения подразумевают наличие трех 

уровней (фильтрация, сегментация, классификация), каждому из которых 

соответствует свой набор алгоритмов. Однако если существующий набор 

алгоритмов на нижнем уровне (фильтрация) позволяет качественно 

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

поставленной задачи, то решения, существующие на среднем уровне 

(сегментация), недостаточны. Так традиционные алгоритмы сегментации 

изображений, основанные на принципе близости точек объекта по цвету 

(яркости) 
и 
не 
учитывающие 
пересечения 
объектов, 
исключают 

возможность 
последующей 
классификации 
с 
приемлемой 
для 

рассматриваемой задачи ошибкой.  

В работах Ковригина А.В., Хлопова А.Н. [1, 2] рассматривался 

синтез 
алгоритмов 
автоматической 
классификации 
объектов 
на 

изображении. На основе синтезированных алгоритмов и принципов 

построения системы машинного зрения, в настоящей статье предложен 

подход  к анализу изображений, содержащих полупрозрачные структуры 

клеток, позволяющий снизить ошибку классификации, в том числе и при 

наличии геометрического пересечения. 

2. Постановка задачи исследования 

Научный журнал КубГАУ, №29(5), май 2007 года 
 

http://ej.kubagro.ru/2007/05/pdf/03.pdf
 

3

Задана исходная выборка 
)
,1
,
,
,
(
n
i
y
x
V
i
i
i
=
=
α
, представляющая собой 

цифровой снимок клеточных структур, где 
i
i
i y
x
α
,
,
 – соответственно 

абсцисса, ордината и яркость i -ой точки изображения. 

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

3. 
Анализ 
изображений, 
содержащих 
полупрозрачные 

структуры клеток 

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

классификации изображения на принципе построения системы машинного 

зрения.  

На первом уровне выполняется предобработка исходной выборки, с 

целью повышения качества распознавания. Предобработка включает в себя 

пересчет значений  яркостей 
i
α  на основе выравнивания гистограммы 

яркостей исходной выборки (1). 

На втором уровне используется принцип алгоритма автоматической 

классификации цветных изображений [1], основанный на построении R  

окрестности  вокруг точек 
)
,
,
(
i
i
i y
x
α
 выбранного класса Ω и проверки 

гипотезы (2) о принадлежности точек 
)
,
,
(
i
i
i y
x
α
, попавших в 
R  

окрестность,  классу 
m
Ω .  

Третий 
уровень 
включает в себя логический контроль на 

однородность класса 
m
Ω  (4)  и, в случае неоднородности, построение 

карты минимальных расстояний по точкам класса 
m
Ω  (5) с последующим 

их перераспределением между дополнительно выявленными классами [2]. 

 

Предлагаемая методика: 

Уровень 1: 

1. Вычислить число точек выборки соответствующих каждому из 

уровней яркости 

∑
=
=
n

i

i
l
l
h
h
1
)
,
(α
,  

Научный журнал КубГАУ, №29(5), май 2007 года 
 

http://ej.kubagro.ru/2007/05/pdf/03.pdf
 

4

 
где  
255
,0
=
l
- уровень яркости, 

 
≠
=
=
l
l
l
h
i

i
i
α
α
α
,0
,1
)
,
(
- единичная функция, определяющая равенство 

значения яркости i -ой точки выборки V  и уровня яркости l . 

2. Вычислить «выровненные» значения яркости 

∑
=
=

i

j
j
i
h

α
α

1
256
*
256
255
 
 
 
 
 
 
 
 
(1) 

3. Сформировать 
из 
выборки 
V  
и 
значений 
i
α  
выборку 

)
,1
,
,
,
(

__
n
i
y
x
V
i
i
i
=
=
α
. 

Уровень 2: 

1. Выбрать i -ую точку из выборки 

__
V  и отнести её к классу 
m
Ω . 

Включить точку 
)
,
,
(
i
i
i y
x
α
 в подвыборку 
m
V
m -го класса,  с 

последующим её удалением из выборки 

__
V . 
2. Определить R  окрестность класса 
m
Ω  из условия попадания в нее 

соседних точек 
)
,
,
(
j
j
j y
x
α
. 

3. Проверить 
все 
точки 
)
,
,
(
j
j
j y
x
α
 
из 
R  
окрестности 
на 
принадлежность их классу 
m
Ω , путем применения решающего 
правила 

>
−
∀
Ω
∉

≤
−
∀
Ω
∈
=
ζ
α
α
α

ζ
α
α
α

j
i
m
i
i
i

j
i
m
i
i
i

y
x

y
x
x
m
)
,
,
(

)
,
,
(
)
(
1
,  
 
 
 
(2) 

где  

i
α - значение яркости точки 
)
,
,
(
i
i
i y
x
α
, принадлежащей классу 

m
Ω , 

j
α - значение яркости точки 
)
,
,
(
j
j
j y
x
α
, принадлежащей R  
окрестности, 

ζ
α
α
α
α
≤
−
=
−
j
i
j
i
3
- мера близости между точками 
)
,
,
(
i
i
i y
x
α
 

и 
)
,
,
(
j
j
j y
x
α
, 

5
=
ζ
 - коэффициент размытости изображения. 
 
При условии принадлежности точки 
)
,
,
(
j
j
j y
x
α
 из R  окрестности 

классу 
m
Ω  включить 
)
,
,
(
j
j
j y
x
α
 в выборку 
m
V
m -го класса,  с 

последующим ее удалением  из 

__
V . 

Научный журнал КубГАУ, №29(5), май 2007 года 
 

http://ej.kubagro.ru/2007/05/pdf/03.pdf
 

5

 
4. При отсутствии точек из R  окрестности, которые можно отнести к 
классу 
m
Ω , следует увеличить счетчик классов на единицу 
1
+
= m
m
 

При отсутствии точек из выборки 

__
V , которые можно отнести к 

классу 
m
Ω  при заданном коэффициент размытости изображении ζ , 

следует перейти к этапу 5, в противном же случае перейти к этапу 1. 

5. Установить индекс текущего класса 
1
=
k
 
6. Построить R  окрестность вокруг точек класса 
k
Ω  и проверить 
гипотезу о принадлежности каждой точки 
)
,
,
(
j
j
j y
x
α
 из этой 
окрестности классу 
k
Ω , путем применения решающего правила 

−
≥
−
∀
Ω
∈

−
<
−
∀
Ω
∈
=

Ω
Ω

Ω
Ω

)
(
)
(
)
,
,
(

)
(
)
(
)
,
,
(
)
(
2
α
α
α
α
α

α
α
α
α
α

r
k

r
k

M
M
y
x

M
M
y
x
x
m
j
j
r
j
j
j

j
j
k
j
j
j
  
(3) 

где  
∑
Ω
Ω
Ω
=

k

k
k

k
M
α
α
1
)
(
- среднее значение яркости по классу 
k
Ω

∑
Ω
Ω
Ω
=

r

r
r

r
M
α
α
1
)
(
- среднее значение яркости по классу 
r
Ω

r - индекс класса, к которому принадлежит точка 
)
,
,
(
j
j
j y
x
α
 
При условии перехода точки из одного класса в другой по правилу 
(3), необходимо пересчитать средние показатели этих классов. Для 
этого следует применять рекуррентные формулы. 

n
n
n

i

i
n
M
n
n
n
n
M
k
k
α
α
α
1
1
)1
(
1
1

1
+
=
+
−
−
=
Ω

−

=
Ω
∑
 , 

n
n
n

i

i
n
M
n
n
n
n
M
r
r
α
α
α
1
1
)1
(
1
1

1
−
=
−
−
−
=
Ω

−

=
Ω
∑
, 

где  
k
Ω - класс, к которому перешла точка 
)
,
,
(
j
j
j y
x
α
, 
r
Ω - класс, 
который она покинула. 
7. При выполнении условия 
m
k
≤
+1
, увеличить индекс текущего класса 
на 1 и перейти к этапу 6, в противном же случае перейти к этапу 8. 
8. Увеличить коэффициент размытости изображения ζ  на величину 

ζ
d и перейти к этапу 5.  
Величину 
ζ
d  выбирать равной 3, пока ζ не достигнет 20, далее 
рекомендуется ζ
d  выбирать равным 1. 
9. Если значение ζ достаточно велико (
[
]
35
,
30
≈
ζ
) следует перейти к 
выполнению уровня 3. 
Уровень 3:

Научный журнал КубГАУ, №29(5), май 2007 года 
 

http://ej.kubagro.ru/2007/05/pdf/03.pdf
 

6

1. Для каждого класса 
k
Ω , 
m
k
,1
=
, вынести решение об однородности 

путем применения решающего правила (4) 

Ω

Ω
>
=
н
неоднороде
класс
иначе

однороден
класс
x
m

k

k
i
,

,
*
)
(

2
2

3
σ
λ
σ
, 
 
 
 
 
(4) 

где  
2
i
σ  - значение дисперсии точек, принадлежащих классу 
m
Ω , 

2
σ - среднее по всем классам 
k
Ω , где 
m
k
,1
=
,  значение 

дисперсии, 

8.0
≈
λ
 - заранее выбранное пороговое значение. 

При условии, что все классы однородны, следует прекратить 

выполнение алгоритма. 

2. Для каждой i -ой точки выборки 
k
V , где k - индекс неоднородного 

класса 
k
Ω , вычислить евклидово расстояние до всех точек, не 

принадлежащих выборке 
k
V ,  и выбрать минимальное 

 
2
2
)
(
)
(
min
j
i
j
i
i
y
y
x
x
d
−
+
−
=
, 
 
 
 
 
 
(5) 

 
где  
i
i y
x ,
- координаты точки 
k
i
i
i
V
y
x
∈
)
,
,
(
α
, 

 
 
j
j y
x ,
- координаты точки
k
j
j
j
V
y
x
∉
)
,
,
(
α
. 

3. Сформировать 
из 
выборки 
k
V  
и 
значений 
id  
выборку 

)
,1
,
,
,
,
(
1
'
n
i
d
y
x
V
i
i
i
i
k
=
=
α
. 

4. Установить начальное пороговое значение 
i
d
T
max
=
 и счетчик 

дополнительно выявленных классов 
1
=
p
. 

5. Выбрать i -ую точку из выборки 
'
k
V , имеющую значение
T
i
d
≥
, 

отнести её к классу 
p
m+
Ω
. 

6. Определить R  окрестность как набор точек 
)
,
,
,
(
i
j
j
j
d
y
x
α
 из выборки 

'
k
V , имеющих значение 
T
d j ≥
 и непрерывно связанных с классом 

p
m+
Ω
. 

Научный журнал КубГАУ, №29(5), май 2007 года 
 

http://ej.kubagro.ru/2007/05/pdf/03.pdf
 

7

7. Отнести все точки 
)
,
,
,
(
j
j
j
j
d
y
x
α
из R  окрестности к классу 
p
m+
Ω
 и 

увеличить счетчик дополнительно выявленных классов на единицу 

1
+
= p
p
. 

8. При наличии точек из выборки 
'
k
V , которые можно отнести к классу 

p
m+
Ω
 следует вернуться к пункту 12. 

9. Уменьшить пороговое значение на единицу 
1
−
=T
T
. 

10.  Для каждого из классов 
l
Ω , где 
p
m
m
l
+
=
,
, определить 
R  

окрестность как набор точек 
)
,
,
,
(
i
j
j
j
d
y
x
α
 из выборки 
'
k
V , имеющих 

значение 
T
d j ≥
 и граничащих с классом 
l
Ω . В случае если точка 

)
,
,
,
(
i
j
j
j
d
y
x
α
 из выборки 
'
k
V  граничит одновременно с несколькими 

классами, ее следует удалить из выборки 
'
k
V  

11. Отнести 
все 
точки 
)
,
,
,
(
j
j
j
j
d
y
x
α
 
из 
R  
окрестностей 
к 

соответствующим классам с последующим удалением этих точек из 

выборки 
'
k
V . 

12. При наличии точек из выборки 
'
k
V , которые можно отнести к одному 

из классов 
l
Ω , где 
p
m
m
l
+
=
,
, следует вернуться к пункту 17. 

13. При наличии точек в выборке 
'
k
V  с значением 
T
d j ≥
, следует 

вернуться к пункту 12. 

14. При условии 
0
>
T
 и наличии  не расклассифицированных точек из 

выборки 
'
k
V  следует вернуться к пункту 16. 

15. Удалить класс 
k
Ω  из списка найденных классов и уменьшить 

счетчик классов на единицу 
1
−
= m
m
. 

16.  При отсутствии неоднородных классов 
k
Ω , где 
m
k
,1
=
, следует 

прекратить выполнение алгоритма, в противном случае следует 

вернуться к пункту 2. 

4. Экспериментальные данные 

Научный журнал КубГАУ, №29(5), май 2007 года 
 

http://ej.kubagro.ru/2007/05/pdf/03.pdf
 

8

С целью проверки эффективности предлагаемого подхода, в качестве 

тестовых примеров, были использованы цифровые снимки мазка 

периферической крови человека.  

На сухое обезжиренное предметное стекло наносили 4 мкл 

периферической крови. Мазок производили шлифованным стеклом под 

углом 30-45°, высушивали на воздухе. Фиксировали мазки в метаноле в 

течении 5 мин., высушивали на воздухе. 

Окраску производили по методу Романовского-Гимза. Высушенные 

мазки помещали в кювету с рабочим раствором краски (1:40) на 45 мин. 

Промывали проточной водой и высушивали на воздухе. 

Далее разные участки полученного мазка фотографировали камерой 

Sony SSC-M370CE через тринокулярный микроскоп МБС-9 с увеличением 

x64 (40x1,6), полученные изображения размером 255x255, 8 bit(grey) 

обрабатывались созданными программными средствами. 

Анализ на основе тестовой выборки, включающей в себя 30 

изображений, показал, что оценка вероятности ошибки составила 0.07. 

а 
б 

 

в 
г 

Научный журнал КубГАУ, №29(5), май 2007 года 
 

http://ej.kubagro.ru/2007/05/pdf/03.pdf
 

9

 

д 
е 

 
 

Рисунок 1 – Примеры работы предложенного подхода (а, в, д – 

исходные фотографии мазков крови; б, г, е – результат работы 

предложенного подхода).

5. Заключение 

Анализ предложенного подхода на основе принципов построения 

системы машинного зрения показал преимущества перед использованием 

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

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

отказаться от метода ручного микроскопического анализа. Важное 

достоинство данного подхода заключено в снижении степени зависимости 

Научный журнал КубГАУ, №29(5), май 2007 года 
 

http://ej.kubagro.ru/2007/05/pdf/03.pdf
 

10

качества решения задачи от цвета, прозрачности, формы и расположения в 

пространстве объектов. 

Литература 

1. Ковригин А.В., Хлопов А.Н. Синтез и анализ алгоритмов автоматической 
классификации цветных изображений // Электронный журнал "Исследовано в 
России". 2005. c. 2417-2422. http://zhurnal.ape.relarn.ru/articles/2005/236.pdf 
2. Ковригин А.В. Многоуровневый алгоритм автоматической классификации 
изображений в условиях отсутствия четких границ между объектами // Вестник 
КрасГАУ. Красноярск, 2006, выпуск № 12, с. 407-411 
3. Форсайт Д., Понс Ж. Компьютерное зрение. Современный подход. : Пер. с англ.М.: Издательский дом «Вильямс», 2004. -928 с.: ил. – Парал.тит.англ. 
4. Непараметрические системы классификации / А.В. Лапко, В.А. Лапко, М.И. 
Соколов, С.В. Ченцов. – Новосибирск: Наука, 2000. – 240 с. 
5. Лапко А.В. Непараметрические методы классификации и их применение. - 
Новосибирск: ВО «Наука», 1993.-152 с 
6. Автоматизированная подсистема распознавания и оконтуривания клеток / А.А. 
Коков, Е.В. Меркулова, В.Г. Адамов, В.Я. Спорыхин // Наукові праці ДонДТУ, 
Донецьк, 2001, випуск 25, с. 110-116 
7. Спорыхин В.Я. Автоматизация процесса исследования функционального 
состояния клеток / В.Я. Спорыхин, Е.В. Адамова // Наукові праці ДонДТУ, 
Донецьк, 2001, випуск 25, с. 110-116 
8. Спорыхин В.Я. Автоматизированная система определения жизнеспособности 
клеток / В.Я. Спорыхин, Е.В. Адамова. // Наукові праці ДонНТУ, Донецьк, 2002, 
випуск 38, с. 110-114