Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Теория признаков распознавания образов на основе стохастической геометрии и функционального анализа

..pdf
Скачиваний:
5
Добавлен:
12.11.2023
Размер:
12.16 Mб
Скачать

2.3. Некоторые интегральные формулы

51

Эти интегралы берутся по всем прямым Gk и Ge, пересекающим F. На основании этого

q dG \ Л А ... Л dGjy

GiC\F=p0

N ( N - 1)

A ... A dGjy

qke dG\ Л

GtnF^0

= N (N - \ ) F SL u- 2.

 

Таким образом, для N случайных прямых, пересекающих вы­ пуклое множество F, среднее число точек пересечения, лежащих внутри F,

N (N -

l)nS

Mg:

(2.26)

L2

 

Аналогично можно найти Mg2. Исполвзуя обозначение интеграла от второй степени хорд (2 .20), запишем окончательный результат:

о

( N \ S , Л1 2 f N \ S 2 , Л1 f N \ h

(2-27)

Mg2

J2 + 24^ 2

~Гч-

 

3 L3

 

Выше мы решали пример и установили, что для окружности радиу­ са R с центром в начале координат / 2 = (16/3) 7гД3 с помощью (2.27) можно определить Мд2 для этого объекта:

Мд2

+ 167Г

2

N

3

 

 

 

Определим теперь среднее число областей г,

 

на которые N случай­

ных прямых делят F. Для нахождения Mr понадобится интерпретация

деления с помощью понятий теории

графов: хорды G/ П F

и грани­

ца 8F образуют плоский граф с г

областями

и g + 2^V

вершинами,

из которых g вершин лежат внутри

F и 2N

вершин —

на

границе.

Можно считать, что через каждую внутреннюю вершину проходят четыре ребра и через каждую вершину на границе — три ребра. Такой подсчет основан на том, что можно не учитывать случаи, когда более двух прямых пересекаются в одной точке или когда две из них пересекаются на границе. Общее число ребер плоского графа тп= 1/2 (4g + 6./V) = 2д + 3N, поскольку каждое ребро соединяет две вершины. Используем подстановку этого выражения в формулу Эйлера g —m + г = 1 :

г = q + N + 1.

Заменяя в этом выражении g средним значением Мд по форму­ ле (2.26), получаем выражение

M r = N { N - 1){TTS /L 2) + N + 1.

(2.28)

Определим среднее число сторон областей. Обозначим пц число сторон области Qi (i = 1 ,2 ,..., г). Как было установлено выше, общее число ребер m = 2q + 3N, причем 2N из них принадлежат границе

4*

52

Гл. 2. Траектории сканирования

и 2q + N являются внутренними ребрами. Заметим, что каждое внут­ реннее ребро является стороной двух областей и каждое ребро границы есть сторона лишь одной области, поэтому

ms = ^ то* = 2(2q + N) + 2N = 4q + 4N.

1

 

В результате

 

Mm.£ = М пц1 = 4N (N —1)7гД + 4N.

(2.29)

L2

 

Таким же образом определяется сумма периметров всех областей:

 

г

N

 

 

 

и = ^ 2 Ui = 2 X ] 9i + L.

 

 

i

i

 

 

Вместе с тем имеет место зависимость

= 7rS/L, справедливость

ее устанавливается из сопоставления (2.1) и (2.3). С учетом этого

 

Mw = 2NMg + L =

+ L .

(2.30)

Число

сторон областей,

на которые

случайные

прямые (?* (i =

= 1 ,

2 делят множество F, определяется как К = m s/r. Вели­

чина К — случайная, формула для нахождения ее среднего значения имеет сложный характер. Поэтому вместо среднего значения ограни­ чимся рассмотрением отношения средних значений:

М *К = 1 Г ^ = 4-Т7ЛГг

, \ 4/ , АГ

7ТТ9-

(2-31)

M r

N (N

—1)TTS + (N

+ 1)L2

 

По этой же причине вместо сложно вычисляемого среднего значе­ ния площади области MQ=S/r будем рассматривать отношение сред­

них значений величин 1:

 

M*Q = - f -

S L 2

(2.32)

М г

N ( N - \ ) TTS + {N + 1 )L2'

Вместо среднего периметра области M u = u /r воспользуемся отно­

шением средних значений:

 

Мм

2-KN S L + L3

М*м:

(2.33)

М ^ “

N ( N - \)TTS + (N + 1)L2'

Пуассоновский линейный процесс. Рассмотрим семейство выпук­ лых областей K(t), зависящее от параметра t, площадью S(t) и пери­ метром L(t), которое сканируется большим числом случайных прямых. В частном случае Ко — некоторый отрезок прямой длины /, содержа­ щийся в K(t). Поставим себе задачу определения вероятности того, что из N случайных прямых, пересекающих K(t), ровно п прямых пересечет и отрезок Ко. Для нахождения такой вероятности следует

1 Степень отличия этих величин можно оценить на следующем примере: при N = 2 для круга получается, что M *Q /M Q = 48/49.

оо
а среднее значение п Мгг = ^~^пР'* = IX, и оно не зависит от ориен-
к 0
тации и расположения отрезка K Q. В общем, линейный пуассоновский процесс означает, что прямые распределены случайным образом так, что в параметрическом пространстве им соответствует пуассоновское поле точек с плотностью XdpAdd. Прямые, соответствующие этому полю, делят плоскость на бесконечное число случайных многоуголь­ ников, параметры которых рассматриваются ниже. Из формул (2.26) и (2.27) следует

2.3. Некоторые интегральные формулы

53

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

Ри

1 )Ш

- !

(2.34)

 

 

Здесь применен результат (2.3), согласно которому вероятность того, что случайная прямая, пересекающая K(t), пересечет Ко, равна отношению длин выпуклых оболочек Ко и Kit). Отрезок прямой Ко можно рассматривать как вырожденный случай выпуклого множества, длина границы которого равна 21. С учетом этого вероятность рав­ на 21/L, она и фигурирует в формуле (2.34).

Рассмотрим теперь предельный случай, когда K it) распространяет­ ся на всю плоскость и число случайных линий сканирования возраста­

ет до бесконечности. При этомlim

т

А

t — Ю О

т

2

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

ти существует пуассоновская система прямых,

рассматриваемая

в стохастической геометрии. При таком процессе

 

Р* =

lim Рп

тп.-lx

(2.35)

t

— Ю О

п\

 

^

а

7гЛ2

7Г2 А4

lim М —=

4 •

16 ■

t — Ю О

S

С учетом этого дисперсия D § = М ^2 - (М

2 —>0.

Из этого факта в соответствии с теорией вероятностей вытекает,

что с вероятностью, равной единице,

 

lim |

(2.36)

t—юо Ь

4

а это есть среднее число вершин на единицу площади.

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

54 Гл. 2. Траектории сканирования

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

Если сфера имеет единичный радиус, а число секущих плоско­ стей N, то при N = 1 число различных областей равно 2, при N = =2 оно равно 4, при N = 3 равно 8. Применив метод математической

индукции, можно доказать,

что число

областей равно 2 + N (N — 1),

ибо на каждом шаге увеличения N

на единицу прибавляется 2N

областей. По аналогии для

N = 2,3,4

число сторон областей равно

числу высекаемых дуг окружностей (соответственно 4, 12, 24). По индукции получаем, что для любого N оно равно 2N (N — 1). Отсюда при увеличении N среднее число сторон будет равно пределу отноше-

4iV(iV —1)

, _

ния 2 [_ jv(iV— 1) ’ т0 есть

“ числителе этого выражения записано

4N (N —1), а не 2N (N —1), поскольку каждая дуга является стороной для двух многоугольников и должна учитываться дважды. Средняя площадь в пределе стремится к 4TT/ N 2, общий периметр 2 + N (N —1) областей равен 4TTN, поэтому средний периметр области в пределе стремится к 4лN. Далее, при больших N рассматривается распреде­ ление областей в малом круге радиуса R на поверхности сферы. На эти области круг рассекается большими окружностями, число которых равно NR. Для упрощения анализа предполагается, что эти области также являются кругами. Затем находят среднюю площадь MQ обла­ стей на плоскости, образованных пуассоновским полем с плотностью dp Add. Она равна 4/лА2, а средний периметр М и равен 4/А.

В стохастической геометрии рассматриваются полученные вероят­ ностными методами средние значения и некоторых других параметров случайных многоугольников, образующихся при пуассоновском линей­ ном процессе [17, 31, 105].

2.4.Развёртка в виде случайных отрезков прямой

ипризнаки распознавания

Сканирование случайно ориентированными и случайно располо­ женными в поле изображения отрезками прямой фиксированной дли­ ны (называемыми для простоты случайными отрезками) оказывается информативнее, чем сканирование случайными линиями. При скани­ ровании случайными линиями единственный признак — интеграл от

2.4. Развёртка в виде случайных отрезков прямой

55

отрицательной степени хорд Д_i — обладает чувствительностью к уг­ лам (/_ 1 является сходящимся при отсутствии углов и расходящимся, если они имеют место на изображении объектов) и может служить как индикатор наличия углов на изображении объекта. Сканирова­ ние отрезками линий со случайными параметрами дает возможность измерять углы у объектов, причем информацию об углах содержит легко определяемый признак — число пересечений изображения с разверткой. Поэтому информационная ценность такого сканирования выше. Изучение свойств развертки в виде случайно ориентированных и случайно расположенных отрезков прямой длиной I начнем с рас­ смотрения сканирования модельных объектов, представляющих собой геометрические элементы: множества, ломаные линии, углы, фигуры, системы параллельных линий.

Пересечение отрезка с геометрическими элементами. Рассмот­ рим пересечение ориентированного отрезка К прямой с выпуклым множеством. Определим меру множества отрезков длиной I, пере­ секающих выпуклое множество Ко, площадь которого равна So, а периметр До. Для этого применим формулу (1.32) для кинематиче­ ской плотности dK. Зафиксируем прямую G как опору, на которой находится отрезок К, обозначая через g длину хорды, высекаемой множеством Ко на прямой:

 

 

 

1 — а

ц{К\ К С К о ф 0 )

dG* Л dt

dG*

dt =

КС\Коф 0

ОПКоф 0

 

— ( 1+а)

= 2

+ /) dG = 2TTS0 + 21L0. (2.37)

 

апкоф0

 

 

Переход от последнего интеграла в этой цепи равенств к оконча­ тельному результату осуществляется в соответствии с формулами (2 .1) и (2.3).

Рассмотрим вырожденный случай, когда выпуклое множество Ко представляет собой отрезок прямой длины /о; периметр такого мно­ жества будет равен 21Q. Применяя вышеприведенную формулу, получа­ ем меру

/х(Д; К П К 0 = 0 ) = 4//0.

(2.38)

Этот результат можно распространить на ломаную линию. Если Ко — ломаная линия, длина которой До. то, применяя последнюю формулу для каждого ее звена и складывая результаты, получаем

Р-{К; К П Ко ф 0 )

n d K =41LQ,

(2.39)

где п — число звеньев ломаной К 0, которые пересекаются с К при данном положении отрезка. Определим меру множества отрезков дли-

2.4. Развёртка в виде случайных отрезков прямой

57

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

М= sm а

sin в sin + в) d6 - [1 + (тг - a ) c t g a ] .

(2.41)

 

о

 

Этой формуле можем придатв вероятностный смысл, если рассмотретв случай, когда отрезок длиной I, пересекающий угол а, пересекает и стороны болвшего угла [3 с общей вершиной в точке О, расположен­ ною вне угла а. Вероятности этого геометрического события

1 + (-7Г- ( 3 ) ctg 1 3

(2.42)

1 + a) ctg а

Рассмотрим сканирование случайными отрезками выпуклого мно­ гоугольника. Предположим, что Ко — выпуклый многоугольник, име­ ющий углы оу при вершинах, площадь SQ, периметр LQ. Обозначим через К ориентированный отрезок длиной /, который не может пересе­ кать две несмежные стороны многоугольника Ко.

Далее, пусть /х* (7 = 0,1,2) есть мера множества положений отрез­ ка К, в которых он имеет г общих точек с контуром многоугольника Ко (Мо мера множества отрезков, размещенных внутри Ко). Тогда пред­ шествующие формулы (2.37), (2.39) и (2.41) можно записать в виде

Mo + Mi + М2 —27rSo + 2/Lo,

Mi Т 2м2 = 4ILQ,

М2 = (/2/ 2) ^ [ 1 + (Tr-oy)ctg оу].

°ч

Из этих уравнений получаем

Mo = 27rS'o - 2ILo + (l2/ 2 ) ^ 2 П + (^ - ay) ctg ay],

(2.43)

М2 = (/2/ 2 ) ^ [ l + (д —ay) ctg ay].

Для приложений к распознаванию образов нас будет интересо­ вать вероятностный смысл этих соотношений. Чтобы выявить его, рассмотрим случай, когда многоугольник Ко располагается внутри выпуклой фигуры Ф, имеющей площадь S<p и периметр Ьф (роль такой фигуры играет сетчатка распознающей системы). В такой ситуации

вероятность пересечения случайным отрезком К длины I выпуклого многоугольника

(2.44)

58

Гд. 2. Траектории сканирования

 

Вероятность пересечения его контура

 

р

_

Mi + 2/Ц2 _

4/Lp

(2.45)

 

 

Мо(Ф)

27г5ф + 21Ьф

 

 

 

Вероятность двукратных пересечений

 

р _

М2

{I2/ 2 ) Е П + (7 r - a J')ctgai ]

 

_ _____ «з__________________

(2.46)

/хо(Ф)

 

27г5ф + 21Ьф

 

 

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

В заключение рассмотрим пересечение ориентированного отрезка длиной I с параллельными линиями на плоскости. В главе 1, пред­ лагая идею распознающей системы, мы оттолкнулись от этой задачи, называемой также задачей Бюффона об иголках. Итак, напомним ее постановку: иголка (т. е. ориентированный отрезок длиной /) бросается случайно на плоскость, на которой начерчены параллельные прямые на единичном расстоянии друг от друга. Требуется определить вероят­ ность того, что иголка пересечет эти прямые. Рассмотрим случай, когда I < 1. Это означает, что возможно только одно пересечение. Обозначим через х расстояние от центра иголки до ближайшей параллели и через 9 угол, составленный иголкой с этой параллелью. Величины х и 9 полностью определяют положение иголки. Тогда, поскольку случайное бросание означает, что х распределена равномерно на (—1/2, 1/2) и 9 распределена равномерно на (0,27г), искомая вероятность

2тг

 

1

/| sin

2/

 

d9

 

 

 

о

 

 

Если I больше единицы, то вероятность одного пересечения

 

arcsin I

2

 

Р

1

——arcsin l + 1

(1 —Z|sin#|) d9

 

 

Если иголка бросается N раз и при п бросках получается хотя бы одно пересечение, тогда Р = n N ~ l является несмещенной оценкой для Р с дисперсией N ~ lP( 1 —Р).

Можно рассматривать вероятность того, что брошенная случайно иголка пересечет прямоугольную решетку, образованную двумя мно­ жествами параллельных линий, пересекающих друг друга под прямым углом. Если линии разделены промежутками а и Ь и I < а, I < Ъ, то

2.5. Сканирование по случайным криволинейным траекториям

59

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

При Ъ—>оо получаем известный результат для вероятности пересе­ чения параллельных линий одного множества 21/жа.

Задача Бюффона обобщена на случай бросания иголки в трех­ мерном пространстве, в котором размещено бесконечное множество параллельных плоскостей на единичном расстоянии друг от друга. Пусть I < 1. Считая, что ось координат направлена перпендикулярно плоскостям, применим полярные координаты для определения направ­ ления иголки. В этих координатах вероятностным элементом является

47г- 1 sin в dO dtp. После интегрирования по его можно записать как (l/2)d(cos0). По этой причине длина проекции lz иголки на ось, пер­ пендикулярную плоскостям, равномерно распределена на (0,1). Так как положение проекции случайно, то вероятность пересечений для всякой данной ориентации равна lz, а это после усреднения по распределению ориентации дает для полной вероятности пересечений значение 0,5/.

Рассмотренные теоретические положения о свойствах пересечений отрезков и выпуклых фигур на плоскости можно распространить на случай пересечения тел отрезками в трехмерном пространстве. Это да­ ет возможность делать выводы о телах на основе свойств пересечений. Приведем в завершение данной темы одно из таких положений. Мера множества отрезков К длиной I, пересекающих выпуклое тело Ко,

y ( K n K o = 0) = 2n2(4Vo + lSo),

где Vo — объем и 5о — площадь поверхности тела.

Приведенные теоретические положения применены при построении распознающих систем, которые рассматриваются ниже в главе 12 и приложении А.

2.5. Сканирование по случайным криволинейным траекториям

Сканирование объектов по криволинейным траекториям со слу­ чайными параметрами (называемых для простоты случайными кри­ волинейными траекториями) ценно тем, что для некоторых видов траекторий дает возможность упростить конструкцию распознающего устройства. Пример такого устройства приводится в этом параграфе ниже. Кроме того, можно показать, что при определенной форме ска­ нирующей кривой (например, двухзвенная ломаная с прямым углом между звеньями) возрастает точность измерения параметров распозна­ ваемого объекта, т. е. повышается надежность распознавания. Соответ­ ственно при заданной точности уменьшается число шагов развертки, необходимое для ее достижения, из-за лучшей сходимости процесса стохастического измерения. Анализ такого эффекта приведен в главе 4, здесь же рассмотрены геометрические свойства пересечений кривых и изображений объектов.

60 Гл. 2. Траектории сканирования

Пересечения кривых. Рассмотрим ситуацию, когда объекты пред­ ставляют собой ломаную линию Гд длиной I/O, а подвижная кривая, случайным образом брошенная на плоскость и пересекающая объект, также есть некоторая m-звенная ломаная Г[ длиной L\. На основании соотношения (2.39), примененного к каждому звену ломаной Гь

щ d K = 4ZjLo.

 

 

Здесь щ — число точек пересечения звена k с ломаной Го.

 

Такие формулы можно получить для

всех звеньев

подвижной ло-

маной Г;. Сложим их левые и правые

части. Считая,

что

т

J2ni = n,

приходим к выводу

 

 

1

 

 

 

n d K = 4LLQ.

 

(2.47)

Здесь п — общее число точек пересечения ломаных Го и Г), инте­ грирование ведется по всем возможным положениям Г). Возможность сложения предшествующих формул для каждого звена следует из того, что на основании свойства инвариантности кинематической меры относительно выбора подвижной системы координат (см. § 1.2) можем выбрать одну и ту же систему координат для всех звеньев k (i = = 1 , 2 , ... , то).

Рассмотрим пересечение объекта со случайными кривыми. Пока­ жем, что формула (2.47) носит универсальный характер и применима к кривым более общего класса, чем ломаные линии. Предположим,

имеется некоторая фиксированная спрямляемая 1 кривая Го длиной Lo, образованная конечным числом дуг, имеющих касательную в каждой точке. Будем полагать, что в некоторой системе координат (0 ,х,у) она определяется уравнениями

х = х(С0); у = у(С0),

где параметр Со является длиной дуги кривой Го.

Обозначим через Г[ подвижную спрямляемую кривую, задаваемую

в подвижной системе координат (А, X, Y)

уравнениями

Х = Х ( С \ ) ; Y =

y(C i),

где С\ — длина дуги Г[.

Для анализа пересечений наиболее удобной формой кинематиче­

ской плотности является

следующая:

 

dK

= |sin.'у| dCo A dC\ A dj,

(2.48)

где 7 — угол между касательными к кривым Го и Г[ в точке пересече­ ния, определяемой значениями параметров Со и С\.

Эта формула очень полезна для приложений. Из нее, в частности, следует, что вероятность того, что угол между кривой Го и случайной

кривой

Г[ в любой точке их пересечения лежит между 7 и 7 + dj

(О < 7 ^

7г), равна 0,5 sin7 d7 , среднее значение угла М 7 = тг/2.

Спрямляемая кривая —это кривая, имеющая конечную длину.