2064
.pdf
Замена исходного изображения его «проекциями» позволяет решать задачи по измерению параметров используемых признаков с удовлетворительной точностью.
Использование результатов теории на практике пока в основном ограничено следующей схемой. Конструируются два параллельно работающих обнаружителя. Один настроен на обнаружение образа с круговой симметрией, второй – протяженного. В качестве первого может быть использован рассмотренный выше «полосовой» фильтр.
Таблица 2.1
Оценка эффективности метода «проекций»
% от размеров поля |
|
|
|
Форма объекта |
|
|
|
|
|
|||
изображения |
протяжен |
протяжен |
круглый |
|
круглый |
круглы |
||||||
|
ный, |
ный, |
, 50% |
|
, 70% |
|
й, |
|||||
|
5% |
20% |
|
|
|
|
|
90% |
||||
Отношение с/п |
1 |
2 |
1 |
2 |
1 |
2 |
|
1 |
2 |
1 |
|
2 |
Безошибочно |
|
|
|
|
|
|
|
|
|
|
|
|
измеренные параметры |
100 |
100 |
100 |
100 |
0 |
22 |
|
7 |
76 |
35 |
|
80 |
признаков |
|
|
|
|
|
|
|
|
|
|
|
|
Объекты, измеренные с |
0 |
0 |
0 |
0 |
75 |
0 |
|
60 |
0 |
15 |
|
0 |
погрешностью более 20 |
|
|
||||||||||
Отклонение |
|
|
|
|
|
|
|
|
|
|
|
|
измеренной площади |
0 |
0 |
0 |
0 |
─ |
15 |
|
─ |
4 |
─ |
|
2,2 |
объектов от |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
действительной |
|
|
|
|
|
|
|
|
|
|
|
|
Более распространен вариант пространственного фильтра, согласованного с изображением распознаваемого образа минимальных размеров, а использование потенциальных возможностей увеличения отношения сигнал/шум реализуется с помощью параллельных линейно-разностных фильтров во временной области (рис. 2.6).
Рис. 2.6. Структурная схема обнаружения объектов: 1 – информация о характеристиках объекта,
ПФ – пространственный фильтр, ВЛРФ – «временный» линейно-разностный фильтр, РУ – решающее устройство
41
Для локальных образов в аналитических исследованиях используется модель гауссоиды
|
|
x |
2 |
y |
2 |
||
|
|
|
|
|
|||
|
|
|
2 |
|
|||
B(x, y) B0 exp |
|
2c |
|
|
|||
|
|
|
|
|
|
||
и весовая функция пространственного фильтра
|
|
x |
2 |
y |
2 |
|
||
h(x, y) h |
|
|
|
|
, |
|||
|
|
|
|
|
||||
exp |
|
|
|
2 |
|
|
||
0 |
|
|
2a |
|
|
|||
|
|
|
|
|
|
|
||
где 
2a и 
2c – эффективные радиусы гауссоид, в которых сосредоточено 63% энергии от объекта и отклика фильтра на сигнал
f (t) B0h0 |
2a2c2 |
|
2t2 |
|
|
|
|||||
|
|
|
|
exp |
|
|
|
|
|
, |
|
a |
2 |
c |
2 |
2(a |
2 |
c |
2 |
|
|||
|
|
|
|
|
|
) |
|||||
где – скорость движения апертуры; t – время.
В момент совпадения центров апертуры и объекта
f t |
|
B h a2 |
exp |
2t2 |
|
|
с |
|
2 |
, |
|||
|
0 0 |
|
4a |
|||
|
|
|
|
|
|
|
а дисперсия шума и отношение сигнал/шум:
2 |
h2a2 2 |
|
4 a2c4 |
|
B2 |
B2 |
||
; ОСШ |
|
|
0 |
M |
0 |
, |
||
|
|
|
||||||
П |
0 |
|
(a2 c2)2 |
|
2 |
2 |
||
|
|
|
|
|||||
где σ – среднеквадратичное значение помех исходного изображения.
Зависимость ОСШ f (a/c) приведена на рис. 2.7.
42
|
Рис. 2.7. Зависимость ОСШ от величины апертуры |
|
при фиксированном значении размера объекта |
Эффект |
от пространственной фильтрации наблюдается при |
a c. Для |
более крупных изображений образов возможности |
пространственной фильтрации не реализуются. Однако выходной сигнал пространственного фильтра содержит информацию, позволяющую повышать ОСШ с помощью «временного» линейноразностного фильтра (ВЛРФ).
ВЛРФ применяется для исключения тренда входной реализации. Его характеристика при настройке на искомый сигнал
1,если |
|
|
|
t |
|
|
T / 4; |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
T / 4 |
|
t |
|
T / 2; |
|
||||||
|
|
|
||||||||||
t 1,если |
|
|
(2.6) |
|||||||||
|
|
|
|
|
|
|
|
|
||||
0,если |
|
t |
|
T / 2, |
|
|
|
|||||
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
||||
где Т – постоянная времени фильтра.
При поступлении на вход фильтра (2.6) суммы прямоугольного импульса (полезный сигнал f t ) длительностью Т /2, аддитивного шума на фоне линейного тренда l t реакция фильтра в момент совпадения центров сигнала и фильтра характеризуется
следующими оценками: АТ/2, 2 kT, ОСШвых= ОСШвхT /4τk , |
где |
τk – интервал корреляции входных шумов. |
по |
При увеличении (уменьшении) длительности сигнала |
сравнению с постоянной фильтра Т его реакция резко снижается. Поэтому предложено использовать параллельно включенные ЛРФ, число которых определяется из следующих соображений. Первый фильтр настраивается на изображение образа минимального
размера min, т.е. его постоянная времени выбирается из условия
43
T1 |
|
|
Для остальных она растет |
по закону |
T |
|
gT |
, |
где |
|||
|
|
2 min . |
1 i |
i |
|
|||||||
g |
4G 1 |
|
|
, G – допускаемый |
|
|
|
|||||
|
8G 1 |
проигрыш |
в |
повышении |
||||||||
|
|
|
||||||||||
|
|
|
2 |
|
|
|
|
Тогда общее |
||||
ОСШвых по сравнению с согласованным фильтром. |
||||||||||||
число запараллеливаемых ЛРФ определяется диапазоном размеров
объектов |
(r |
r ) |
и величиной G: M loga rmax /rmin . |
Для |
|
r |
/r |
min |
max |
N = 6. |
|
= 2 |
G 1,1 |
|
|||
max |
min |
|
|
|
|
Алгоритм распознавания протяженных образов также основан на совмещении операций «обнаружение – опознавание». Для реализации этого алгоритма используется пространственный фильтр в виде вращающего штриха. При совмещении штриха с изображениями образа выходной сигнал принимает наибольшее значение. Он отражает информацию: где находится образ (обнаружение) и «этот образ имеет протяженную форму (опознавание)».
При теоретических исследованиях распределения интенсивности сигнала протяженного образа и весовую функцию
штрихового |
|
анализатора |
|
|
описывают |
|
несимметричными |
|||||||||||||||||||
гауссоидами: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
B x, y B |
|
|
|
x |
2 |
|
y |
2 |
|
h x, y h |
|
|
|
|
x |
2 |
|
y |
2 |
|
||||||
exp |
|
|
|
|
|
|
; |
|
exp |
|
|
|
|
|
|
, |
||||||||||
|
|
2 |
|
|
2 |
|
|
|
2 |
|
|
2 |
||||||||||||||
c |
0 |
|
|
2c |
|
2d |
|
(2.7) |
0 |
|
|
|
2a |
|
2b |
|
||||||||||
|
|
|
|
|
|
|
|
|
c и |
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
где |
|
d – |
длины полуосей |
||||||||||
эллипса |
равной |
интенсивности, |
на |
границе |
которого |
|
уровень |
|||||||||||||||||||
сигнала |
снижается |
в |
e 2,718раз; |
aиb– эффективные размеры |
||||||||||||||||||||||
апертуры, |
определяемые аналогично полуосям c |
|
и |
|
d; В0 – |
|||||||||||||||||||||
максимальная |
интенсивность |
сигнала; |
h0 |
|
– |
|
максимальная |
|||||||||||||||||||
чувствительность анализатора. Уравнения (2.7) соответствуют случаю совпадения координат сигнала и пространственного фильтра. При развертке (движении центра объекта вдоль оси x) с постоянной скоростью
B0 x, y,t |
|
|
|
|
|
|
|
|
2 |
|
x t sin ycos |
2 |
|
B |
хexp |
x t cos ysin |
|
|
, |
|
|
|
|
||||
0 |
|
2c2 |
|
2d2 |
|
|
|
|
|
|
|
||
(2.8)
где γ – угол между осью ординат и наибольшей осью изображения дефекта; t – время.
44
Совместим центр штрихового анализатора с началом координат и зададим ему вращение с постоянной скоростью. Тогда его описание записывается в виде
h x, y,t h |
exp |
xcos ysin 2 |
xsin ycos 2 |
|
, |
|
0 |
|
|
2a2 |
2b2 |
|
|
|
|
|
|
|
||
(2.9)
где 2 t /TB – угол поворота штриха; TB – период вращения. Реакция фильтра (2.9) на сигнал (2.8)
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
t |
2 |
sin |
2 |
|
|
1 |
|
1 2 |
|
||||||||||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
2 |
|
|||||||||||||||||
f t B h |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c |
|
d |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
exp 2t2A |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
0 0 |
|
4 A B C D E2 |
|
|
|
|
|
|
|
|
|
|
|
4 C D |
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
2 |
|
2 |
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
t |
|
|
[4AC D ] sin cos E |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
d2 |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
4C D [4 A B C D E2 |
] |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
B h |
|
2 |
abcd |
exp |
2t2R |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2.10 |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
0 |
|
0 |
|
|
y |
|
|
4C D y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
где y a2 |
c2 b2 |
d2 b2 |
a2 d2 c2 sin2 ; |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
2 |
2 |
2 |
|
2 |
|
1 |
|
|
1 |
|
|
|
|||
R a |
b |
c |
d |
|
|
|
|
|
|
4AD |
|
|
|
4BC |
|
|
2 |
d |
2 |
2 |
b |
2 |
|||||||||
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
c |
|
a |
|
|
||||||
|
2 |
|
2 |
|
2 |
|
2 |
|
1 |
1 |
2 |
|
1 |
|
1 |
2 |
||
sin |
|
cos |
|
sin |
|
cos |
|
|
|
|
|
|
|
|
|
|
|
; |
|
|
|
|
|
b2 |
|
d2 |
|||||||||||
|
|
|
|
|
|
|
|
a2 |
|
|
c2 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
45
A |
cos2 |
|
sin2 |
|
|
B |
cos2 |
|
sin2 |
|
||||||||||
|
|
|
|
|
|
|
, |
|
|
|
|
|
|
|
|
; |
|
|||
2c2 |
|
|
|
|
|
|
|
|
|
2b2 |
|
|
||||||||
|
|
|
|
2d2 |
|
|
2a2 |
|
|
|
|
|||||||||
C sin2 /2c2 cos2 /2d2 ; |
|
D sin2 /2a2 cos2 /2b2 ; |
||||||||||||||||||
|
|
|
|
1 |
|
|
1 |
|
|
|
|
1 |
|
1 |
|
|||||
F sin cos |
|
|
|
|
|
|
|
sin cos |
|
|
|
|
|
. |
||||||
|
|
|
d2 |
|
|
b2 |
||||||||||||||
|
|
|
c2 |
|
|
|
|
|
|
a2 |
|
|
||||||||
Сигнал на выходе вращающегося штрихового анализатора (2.10) значительно отличается по форме и длительности от сигналов, получаемых при сканировании локальных образов анализатором симметричной формы.
Наличие периодических пульсаций, обусловленных вращением анализатора, позволяет выйти на алгоритм распознавания:
разложение отрезков f t в ряд Фурье и вычисление гармоники с максимальной амплитудой. Если эта амплитуда превышает установленный порог принятия решений, объект считается опознанным. Во временной плоскости максимальная амплитуда сигнала находится из выражения
2 abcd
0 max B0h0 
a2 c2 b2 d2
при условии , а минимальная (при взаимноперпендикулярном положении штриха относительно дефекта: / 2
2 abcd
0 min B0h0 
a2 d2 b2 c2 .
46
Разность |
F 0 max 0 min |
характеризует амплитуду |
полезного сигнала. |
|
|
Доказано, что при сканировании изображения с равномерной |
||
спектральной плотностью и дисперсией шумов 2 двумя взаимноперпендикулярными анализаторами дисперсия шума на выходе вычитающего устройства определяется формулой
b2 |
2h2 |
|
ab |
|
|
2 ab 1 |
|
. |
|||
a2 b2 |
|||||
|
|
|
|
Тогда отношение сигнал/шум на входе решающего звена
ОСШВ |
в02 |
|
2 abc2d2 a2 b2 |
|
||||||||||
|
a2 |
b2 ab |
|
|||||||||||
в2 |
||||||||||||||
|
|
1 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|||
a2 |
c2 b2 d2 |
|
|
|
|
|
|
|||||||
|
|
|
|
a2 d2 b2 c2 |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
При совпадении размеров апертуры и объекта в результате согласованной фильтрации увеличение отношения сигнал/шум составляет
ab a b 4
2 a2 b2 ab а2 b2 .
Проиллюстрируем эффективность согласованной пространственной фильтрации на примере.
Объект смоделирован в виде прямоугольной щели длиной 5 мм с одинаковыми размерами по глубине раскрытия и ширине –
наименьший, |
выявляемый образ: (0,1 0,1 5)мм3, |
а также |
(0,2 0,2 5)мм3 |
и (0,4 0,4 5)мм3. Результаты |
обработки |
изображений этих образов с помощью пространственных фильтров с круговой симметрией приведены в табл. 2.2, 2.3.
Таблица
2.2
Амплитуда сигнала от распознаваемого образа на выходе
47
пространственного фильтра с круговой симметрией, ед. оптич. плотности
Размеры апертуры, |
мм2 |
|
Размеры образа, |
мм2 |
|
|
|
|
|
||
|
|
0,1 5 |
|
0,2 5 |
0,4 5 |
0,1 0,1 |
|
0,01 |
|
0,02 |
0,04 |
0,2 0,2 |
|
0,005 |
|
0,02 |
0,04 |
0,4 0,4 |
|
0,0025 |
|
0,01 |
0,04 |
Таблица
2.3
ОСШ на выходе пространственного фильтра с круговой симметрией
Размеры апертуры, |
мм |
|
Размеры образа, |
мм |
||
|
|
|
|
|
||
|
|
0,1 5 |
|
0,2 5 |
|
0,4 5 |
0,1 0,1 |
|
0,21 |
|
0,83 |
|
3,30 |
0,2 0,2 |
|
0,13 |
|
2,04 |
|
8,16 |
0,4 0,4 |
|
0,052 |
|
0,83 |
|
13,22 |
Как следует из приведенных данных, уверенно можно зарегистрировать только объекты размерами (0,4 5)мм2. Протяженные образы с меньшими размерами по ширине обнаружить проблематично.
Иной вывод вытекает из анализа результатов обработки тех же данных штриховым анализатором (табл. 2.4 и 2.5).
Таблица
2.4
Амплитуда сигнала от образа на выходе штрихового анализатора единиц оптич. плотности
Размеры апертуры, мм2 |
|
Размеры образов, мм2 |
|
|||
|
|
|
|
|
|
|
|
|
0,1 5 |
|
0,2 5 |
|
0,4 5 |
0,1 5 |
0,5 |
|
1 |
|
2 |
|
0,2 5 |
0,125 |
|
0,5 |
|
1 |
|
0,4 5 |
0,031 |
|
0,125 |
|
0,5 |
|
|
|
|
|
|
|
Таблица |
ОСШ на выходе штрихового анализатора |
2.5 |
|||||
|
||||||
Размеры апертуры, мм2 |
|
|
Размеры образов, мм2 |
|
||
|
0,1 5 |
|
0,2 5 |
|
0,4 5 |
|
|
|
|
|
|||
0,1 5 |
|
12,38 |
|
49,5 |
|
198 |
0,2 5 |
|
6,01 |
|
96,1 |
|
384 |
0,4 5 |
|
2,37 |
|
37,9 |
|
606 |
|
48 |
|
|
|
|
|
Поскольку дисперсия шумов на выходе фильтра увеличивается в меньшей степени, чем амплитуда сигнала, итоговое ОСШ возрастает настолько (табл. 2.5), что проблема надежного обнаружения протяженных образов снимается с повестки дня. Однако следует иметь в виду потери, которые сопряжены с применением таких фильтров. В каждой точке растра изображения необходимо найти свертки анализируемого снимка и штрихового анализатора при его разных положениях относительно декартовой системы координат. Имеющиеся в литературе данные дают разумную оценку увеличения временных затрат на обработку – не менее чем в 20 раз. Возможности современных ПЭВМ не ограничивают реализацию таких фильтров, работающих в реальном масштабе времени.
В качестве второй частной задачи рассмотрим комплексную биометрическую систему по защите информационных ресурсов от внутренних угроз.
На сегодняшний день неавторизованный доступ можно считать одной из самых серьезных угроз информационной безопасности.
По данным Information Security Breaches Survey 2010, угроза неавторизованного доступа стоит на третьем месте в списке самых опасных и губительных для бизнеса угроз информационной безопасности.
Проблемы ограничения доступа к информационным ресурсам и терминалам информационно-вычислительных сетей или изолированным ПЭВМ могут быть решены введением биометрических систем идентификации пользователя. Существуют системы идентификации человека по особенностям геометрии лица, руки, рисунку кровеносных сосудов на лице, тыльной стороне кисти, по радужной оболочке глаз. Однако системы аутентификации, построенные на базе «открытых биометрических образов», характеризуются незащищенностью от подделок предъявляемых образов. При удаленной аутентификации пользователь недоступен для прямого наблюдения за ним службой безопасности, нельзя выявить и пресечь попытки предъявления муляжа «легитимного пользователя». Другая проблема биометрических систем связана с неприемлемыми ошибками первого и второго рода, которыми характеризуются предлагаемые для использования разработки. Перспективным направлением ее
49
решения считается разработка комплексных систем идентификации. Поэтому рассмотрим несколько независимых признаков, описывающих особенности двигательных навыков человека, которые проявляются при выполнении операций на подсознательном уровне, а именно особенности набора паролей на клавиатуре ПЭВМ, речевой пароль, параметры траектории и давления пера на планшет при написании автографа (пароля). Структурная схема комплексной системы приведена на рис. 2.8.
Рис.2.8. Схема работы системы: AB – априорные вероятности
В качестве признака для идентификации по речевому паролю используется преобразованный амплитудный спектр речевого сигнала, полученный путем сглаживания текущего спектра F(vj):
j n
F * ( j ) F ( i ) /(2n 1)
i j n
где n – эмпирически подбираемая величина, с последующим интегрированием в окрестности экстремумов.
|
|
i x |
|
|
|
F * j |
|
F ** ( j |
) |
j i x |
|
2x 1 |
|
||
|
|
|
50
