Файл: Иваницкий Г.Р. Исследование микроструктуры объектов методами когерентной оптики.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 22.06.2024
Просмотров: 67
Скачиваний: 0
После промежуточных преобразовании, положив
<г(х, y)z(x', |
y ')> = R,(Ах, |
Ay), |
|
где Ах = х'—х, Ау = у'— у, |
можно |
прийти |
к выражению |
2Lх 2L |
|
|
|
° J = z i r J |
j1 |
A y ) - |
оо
- < г ( х , |
(148) |
Несложно показать [Л. 39], что предполагая отсутст вие корреляции в распределениях плотностей эталона и шума, (148) можно представить в виде
2I.x 21у
4 = 1 ^ 7 ( f Яш (А-v. ДУ) ( 1 - ^ - ) ( 1 - 2 1 7 ) dAx
оо
(149)
где Е — энергия эталонной функции; Rm(Ax, Ay) — функ ция корреляции шума.
Таким образом, дисперсия оценок корреляционной функции зависит от протяженности анализируемых вы
борок и функции |
корре |
|
|
|
|
|
|||||
ляции шума. |
|
величина |
|
|
|
|
|
||||
Предельная |
|
|
|
|
|
||||||
Lx, определяющая |
грани |
|
|
|
|
|
|||||
цы |
существования |
вход |
|
|
|
|
|
||||
ного |
распределения |
(в |
|
|
|
|
|
||||
направлении |
х, |
рис. |
14) |
|
|
|
|
|
|||
в общем случае ограниче |
|
|
|
|
|
||||||
на |
параметрами оптиче |
|
|
|
|
|
|||||
ской системы и частотным |
|
|
|
|
|
||||||
составом |
изображения |
|
|
|
|
|
|||||
(выражение для Ьу иден |
F hc. 14. |
Соотношение |
между |
про |
|||||||
тично) |
|
|
|
|
тяженностью |
входного распреде |
|||||
|
|
|
|
ления, |
максимальным |
углом |
ди |
||||
Lx^ A —2FtgQ, |
(ISO) |
||||||||||
фракции |
и апертурой |
системы |
|||||||||
где |
А — диаметр |
|
линзы; |
LX^ A — 2F |
— - 1 |
|
|
||||
Л — ее фокусное |
расстоя- |
|
|
> 4тс2 — А2о)- |
|
||||||
ине; |
0 — максимальный |
|
|
|
|
|
угол дифракции, соответствующий наивысшей простран
ственной^ частоте, пропускаемой |
системой при |
данном |
Lx. В свою очередь sin 0=A co/2jt, |
где К — длина |
волны |
g* |
|
83 |
используемого освещения; со — пространственная |
ча- |
||||
стота. |
|
|
|
|
|
Таким образом: |
jA |
|
|
|
|
Lx< |
V-№ — l-w- |
2F. |
(151) |
||
2 F |
|||||
Отношение Aj2F имеет смысл относительной апер- |
|||||
туры системы. Очевидно, что при |
|
|
|||
|
_Л |
Хсо |
|
|
|
|
21- > Y 4я2 -- Л2С02 |
|
|
оптическая система не ограничивает полосы частот вход ного распределения. Предельное значение угла дифрак ции, еще не ограничиваемого оптической системой, удов
летворяет соотношению |
|
|
|
А _ |
Хсо |
(152) |
|
IgO = - F |
Y - I tz2 — \ 2ы 2 |
||
|
|||
При этом входное распределение стягивается в точку. |
|||
Дальнейшее уменьшение |
отношения A/2F приводит |
||
к сокращению полосы частот входа. |
|
Анализ статистических характеристик шумовых со ставляющих изображений объекта показал, что коэффи циент корреляции гш(Дх) = /?ш(Дх)/а2ш достаточно хоро
шо аппроксимируется функцией вида |
|
||
|
Гш(А.т) =ехр[— а|Ах|]. |
(153) |
|
Подставив выражение для коэффициента корреляции |
|||
в (149) и рассматривая |
случай одной переменной, |
окон |
|
чательно получим: |
|
|
|
|
2 |
|
|
= |
~TL | ( 1 ~ |
ехР (— а I Аа' I ) dhx = |
|
= |
-Щ - |
i l1exp ( - 2*LX)]}. |
(154) |
Коэффициент затухания корреляционной функции а является величиной, обратно пропорциональной интер
валу корреляции Дхкор, |
который в свою очередь связан |
||
с полосой частот шумового процесса |
соотношением |
||
&хКОр = 2п/А(ох. Поэтому |
выражение (154) |
можно |
пере |
писать так: |
|
|
|
2 |
А<Дд-L н |
(155) |
|
gr |
71 |
. |
|
|
|
||
84 |
|
|
|
Характер |
изменения |
ди |
|
|||||
сперсии в функции произве |
(8т/см2-ср)2 |
|||||||
дения |
Au>xLx |
показан |
на |
|
||||
рис. 15. Анализ зависимости |
|
|||||||
ozR=f(AioxLx) |
показывает, |
|
||||||
что при заданной протяжен |
|
|||||||
ности |
объекта |
во |
входной |
|
||||
плоскости наибольшее укло |
|
|||||||
нение оценок функции кор |
|
|||||||
реляции |
наблюдаются |
для |
|
|||||
низкочастотных |
составляю |
|
||||||
щих |
шумовых |
компонентов |
|
|||||
изображений объекта, т. е. |
|
|||||||
для сильнокоррелированных |
Рис. 15. График изменения ди |
|||||||
изменений формы |
объекта. |
сперсии оценок корреляцион |
||||||
ной функции в зависимости от |
||||||||
Расширение |
полосы частот, |
|||||||
произведения 'Дм.*/.*. |
||||||||
занимаемой |
шумовыми |
со |
|
ставляющими, вызывает уменьшение дисперсии по зако ну рис. 15. Полагая, что величина ошибки не может быть больше некоторой заданной величины, можно опре делить значение допустимой дисперсии в пределах класса.
Полученные выводы могут быть также использованы для оптимального с точки зрения минимизации диспер сии оценок корреляционной функции синтеза эталонов.
Например, в структуре анафазных шапок (анафазная стадия деления клетки обычно используется для ис следования хромосомных повреждений) медленные, сильнокоррелированные изменения формы, как правило, связаны с участками, обращенными к центру митотиче ского аппарата и образованными центромерами хромо сом, движущимися вдоль нитей веретена. Такого же рода изменения наблюдаются на боковых сторонах шап ки. Быстрые, слабокоррелированные изменения формы чаще всего связаны с областями, формируемыми хроматндными нитями. Поэтому для обеспечения минималь ной дисперсии оценок корреляционной функции в осно ву классификации целесообразно положить признаки, определяющие устойчивость сильно коррелированных изменений в пределах класса. Другими словами, форми рование классов нужно производить таким образом, чтобы в пределах класса сохранялось максимальное по добие элементов, связанных с сильнокоррелированными изменениями формы. Эталоны, синтезированные для фор-
85
В качестве таких признаков при 'классификации формы анафазных шапок были использованы следующие приз наки (рис. 19): 1) площадь; 2) протяженность стороны шапки, обращенной в сторону полюса 'митотического ап парата и образованной центромерами хромосом, движу щимися вдоль нитей веретена; 3) углы, образованные стороной, обращенной к полюсу, н боковыми сторонами; 4) отношение высоты шапки к ее ширине (на уровне половины высоты); 5) характер «бахромы», образован ной хроматидами. Классификация по указанным приз накам позволила выявить устойчивые формы, доста точно распространенные в исследуемой популяции клеток.
Однако с точки зрения проведенного анализа диспер сия оценок корреляционной функции может быть умень шена, а границы класса несколько расширены, если в пределах класса сохранять в первую очередь устойчи вость сильнокоррелированных изменений формы.
Представляет интерес исследование вопроса о харак тере изменения дисперсии оценок корреляционной функ ции за счет флюктуаций формы при регулировке дискри минирующих свойств фильтра.
Целесообразность такого рода преобразований ста новится очевидной, если иметь в виду, -что этот путь по зволяет наиболее рационально использовать информаци онную емкость двумерного согласованного фильтра. Как правило, подобная операция осуществляется изменением коэффициента передачи оптической системы для двумер ной согласованной фильтрации. Так, уменьшение селек тивных свойств фильтра может быть осуществлено вве дением пространственного интегрирования (коэффици ент передачи системы уменьшается при этом обратно пропорционально пространственным частотам со* и соу) или просто блокированием верхних пространственных частот. Последнее эквивалентно введению в тракт фильт ра -низких пространственных частот. Обратная операция, имеющая целью увеличить селективные свойства фильт ра, сводится к ослаблению удельного веса низких про странственных частот. Она может быть выполнена вве дением в оптическую систему дополнительного фильтра верхних пространственных частот или осуществлению пространственного дифференцирования.
Сигнал на выходе такой системы можно представить в виде свертки аппаратной функции системы и входной реализации
88
со
я * (АХ. A|/) = J J h ( x - x ' , y - y ' ) f BK( x ’ , y ' ) d x ’ d y'. (156)
—00
При этом аппаратная функция представляет произведение аппа ратной функции согласованного фильтра и аппаратной функции дополнительного фильтра
Л(х, у) = lic .ty[х, |
</)Лд.ф'(*- У)- |
(‘157) |
. Можно показать, что дисперсия |
оценок корреляционной |
функции |
с учетом дополнительного фильтра определяется выражением [Л . 40]
00
|
|
|
° к = |
|
л®ня(д*. *у)~ |
|
|
|||
|
|
|
—00 |
|
|
|
|
|
|
|
|
|
— |
< f n x ( x ', |
y ' ) f э т (х ', |
у ') |
> |
2] d h x dAy, |
(158) |
||
где Rn(Ax, A y ) — |
автокорреляционная |
функция аппаратной |
функции |
|||||||
системы; |
R(Ax, |
Ay) — |
корреляционная |
|
функция |
произведения |
||||
[fn x]*, у) fat (.v, (/)], т. е. |
если |
дисперсия |
оценок |
корреляционной |
||||||
функции при неограниченной полосе частот оптической |
системы |
|||||||||
определяется |
корреляционной |
функцией |
произведения |
у) X |
||||||
Х 1ат{х, |
(/)] |
[выражение |
(148)], то при |
введении дополнительного |
фильтра она описывается интегральным выражением типа свертки (158). В этом выражении в качестве одной из свертываемых функ ций используется автокорреляция аппаратной функции системы, а вторая свертываемая функция пропорциональна функции корре ляции произведения входного и эталонного распределений.
Рассмотрим более детально каждый из сомножителей в урав нении (158). В предположении независимости распределений плот ности эталона и шума корреляционная функция произведения вход ного и эталонного распределений R(<Ах, А у) определяется корреля ционной функцией шума.
Величина < f Bx'(x, y ) f B?(x, у ) > 2 эквивалентна энергии эталонной функции.
В ряде случаев полоса пространственных частот системы огра ничена апертурой используемых линз или специально вводимыми
вспектральную плоскость диафрагмами.
Вкачестве примера ниже рассмотрена автокорреляция аппа ратной функции оптической системы при изотропном ограничении
полосы пространственных частот.
Коэффициент передачи такого фильтра нижних частот описы
вается выражением |
|
|
|
|
|
Н |
(со*, |
сон) = g | при |
со2, + со2 | ^ |
\ |
(158а) |
где а — радиус |
зоны |
прозрачности |
в частотной |
плоскости. |
|
В этом случае зона прозрачности представляет собой функцию— круг. Поэтому аппаратная функция подобной системы определяется преобразованием Фурье от функции, описываемой выражением (158а). Как известно [Л . 50], преобразование Фурье такой функции определяется как
£ ( Р ) = ^ Д |
(159) |
89