Дуда Р., Харт П. - Распознование образов и анализ сцен (1033979), страница 64
Текст из файла (страница 64)
Разработайте другой оператор для обнаружения текстур типа «елочхн», как на твидовой ткани. 6. Сформулируйте критерий для процедуры объединения областей, который способствовал бы созданию и сохранению областей с прямолинейными границами. 7. Постройте алгоритм прослеживания контуров, пригодный для выделения контуров 8.связных дискретных объектов. Глава В АНАЛИЗ ПРОСТРАНСТВЕННЫХ ЧАСТОТ эгй ВВЕДЕНИЕ В этой главе мы хотим определить операцию двумерного преобразования Фурье и обсудить некоторые ее применения в анализе сцен.
Наш интерес к этому преобразованию вызван двумя причинами: оно дает одновременно и хороший аппарат для исследования некоторых теоретических аспектов обработки изображений, и удобные средства для реализации таких операций, как сглаживание, повышение резкости и сравнение с эталоном. Прежде чем переходить к деталям, мы должны сделать два замечания. Во-первых, мы предполагаем, что читатель хотя бы слегка знаком с одномерным преобразованием Фурье, которое широко используется в теории связи. Во-вторых, читателя, склонного к цифровым методам, мы спешим уверить, что, хотя наше обсуждение будет проводиться полностью на языке аналоговых функций интенсивности, имеются весьма эффективные цифровые методы вычисления результатов преобразования Фурье для дискретного изображения, Вообще говоря, для выполнения этого преобразования имеются как хорошие аналоговые, так и хорошие цифровые методы.
Одновременное развитие оптических средств для изображений в аналоговой форме и так называемого быстрого преобразования Фурье для изображений в цифровой форме удваивает наш интерес к обсуждению возможных применений частотных методов в анализе сцен. Наше обсуждение мы начнем с нескольких фундаментальных определений. Пусть имеется аналоговая функция интенсивности у(х, у), заданная на бесконечной плоскости. Определим преобразование Фурье У, создающее спектр Фурье 6(1„, 1„), следующим образом: 6(~„, ~„) =еГ(у(х, у))= = ) ) у (х, у) ехр( — 2я(((„х+ ~„у)]Ых Ыу. (1) Определим также обратное лреобразоеание Фурье а ', переводящее спектр 6(г„, Г„) в исходную функцию д(х, у): у (х у) = У' '(6(1.. 1,)) = = ~ ~ 6 (1„, ~„) ехр[2п((Г„х+ ~„у)]б~„г()„, (2) Гл.
а. Анализ ираспранственньи «а«тат Оставим пока в стороне вопросы математического характера о существовании определенных интегралов и посмотрим, в чем заключается смысл этих определений. Можно считать, что формула (2) описывает разложение функции интенсивности я(х, у) в ряд, который представляет собой «обобщенную сумму» экспоненциальных функций комплексного аргумента. Для каждой пары значений пространственных частот г„и га в обобщенной сумме имеется одна экспонента, которая умножается на весовой коэффициент 6(~„, 1„).
Как же получаются эти весовые коэффициентыу Согласно фундаментальному результату анализа Фурье, называемому обратной теоремой, весовые коэффициенты определяются по формуле (1). Таким образом, спектр Фурье функции я(х, у) можно считать просто набором весовых коэффициентов ее разложения в ряд, представляющий собой (обобщенную) сумму экспоненциальных функций.
Вернемся на короткое время к вопросу о существовании пары преобразований Фурье (1) и (2). Ответ на этот вопрос может показаться не слишком элегантным по форме. Он заключается в следующей совокупности обычно принимаемых достаточных условий: 1) функция д должна быть абсолютно интегрируемой на всей плоскости изображения; 2) функция д должна иметь конечное число разрывов и конечное число максимумов и минимумов в любой прямоугольной области конечных размеров; 3) функция у не должна иметь бесконечных разрывов.
Хотя вопрос о существовании преобразования Фурье интересен с математической точки зрения, мы в дальнейшем им заниматься не будем, поскольку проблемы, связанные с существованием определенных интегралов, на практике, как правило, никаких трудностей не вызывают. Нам хотелось бы рассмотреть несколько более детально, какие следствия вытекают из определения пары преобразований Фурье и в'чем заключается их смысл. В частности, будет показано, что высокие пространственные частоты соответствуют резким изменениям интенсивности на изображении. Мы видели, что спектр функции интенсивности определяет весовые коэффициенты ее разложения в ряд, составленный из комплексных экспонент; поэтому для понимания смысла преобразования нужно хорошо представлять себе, как выглядит экспоненциальная функция ехр 12п((1„х+ ~„у)1 для заданных значений 1„и га, К сожалению, эта функция двух переменных х и у принимает комплексные значения, и ее трудно изобразить графически.
Некоторое представление о ее виде можно получить, если пометить на плоскости (Х, У) области, в которых эта функция вещественна и положительна. Эти так называемые области нулевой фазы можно найти, приравняв показатель экспоненты 8Л. Вн«дениз величине 2п(п для всех целых и. Это приводит к уравнению 2п( (Г„к+Гау) = 2п(п, или л у= — — х+ —. « Еа Таким образом, как показано на рис. 8.1 '), области нулевой фазы представляют собой прямые параллельные линии.
Наклон каждой Рис. 8,1. Геометрическое место точек нулевой фазы. (Из книги Дж. Гудмэна «Введение в фурье-оптикун Авторские права 1968 г. издательства МсС«гатт-НИ1, 1пс. Воспроизводится с разрешения Мсбгатч-Н1П Воо)«Сошрапу.) прямой равен — (Г„чг), а их общая нормаль ориентирована под углом О=иге(и(ГД ). Расстояние Е между этими прямыми называется пространственным периодом и определяется формулой 1 'г' 1'+~~ Следовательно, чем выше пространственные частоты, тем теснее расположены линии нулевой фазы.
С каждой точкой (~„, Г„) на плоскости пространственных частот связаны два геометрических параметра: угол ориентации О н интервал Е. т) Из книги Гудмэна (1968). 322 Гл. 8. Анализ прастаанстззннслх частот Попробуем теперь понять, что означает появление пространственной частоты с большой амплитудой в спектре Фурье некоторого изображения. Пусть, например, у нас есть изображение у(х, у), спектр которого СЦ„, 1а) достигает большой величины, скажем, в точке (и, о). Поскольку значение 6(и, о) велико, слагаемое 6(и, о) ехр (2п ((их+оу)! вносит значительный вклад в обобщенную сумму экспонент в формуле (2). Легко показать, что, поскольку функция у(х, у) вешественная, справедливо равенство С(~н, Да)= =6" ( — ~„, — га), где звездочка обозначает комплексно сопряженйую величину.
Так как модули комплексно сопряженных величин равны, то величина 6 ( — и, — о) равна по модулю 6(и, о), и поэтому слагаемое 6( — и, — о) ехр [2пс( — их, — оу)] также вносит значительный вклад в сумму экспонент формулы (2). Предположим теперь, что функция 6 Д„, )з) мала по модулю всюду, кроме точек (и, о) и ( — и, — о). Тогда изображение достаточно хорошо аппроксимируется формулой 6 (и, о) ехр [2пс (их+ оу)) + 6 ( — и, — о) ехр (2пс ( — их — оу)), поскольку зто слагаемое преобладает в сумме экспонент выражения (2). Легко показать, что это слагаемое представляет собой вещественную величину; если изобразить ее графически, она будет выглядеть, как волнистая поверхность синусоидальной формы, гребни которой образуют параллельные линии, подобные изображенным на рис.
8.1. Следовательно, каждая симметричная пара пространственных частот (и, о) и ( — и, — о) вызывает появление в обобщенной сумме (2) слагаемого, представляющего изображение параллельных полос с изменением интенсивности по сииусоидальному закону. Чем больше модуль спектра в точке (и, о), тем значительней вклад, вносимый этим слагаемым, Пусть, например, у нас есть функция интенсивности д(х, у), состоящая из вертикальных темных и светлых полос. Если интенсивность при переходе от темных полос к светлым меняется по синусоидальному закону, спектр Фурье СД„, )а) не равен нулю только в двух точках плоскости пространственных частот, причем обе точки лежат на осн г"„. Если перепады интенсивности более резкие, спектр Фурье будет отличен от нуля более чем в двух точках; однако вне оси ~„он по-прежиему равен нулю, так как полосы на изображении вертикальны.
Как и в случае одномерного преобразования Фурье, чем резче перепады интенсивности, тем больше будет величина 6()„, О) для больших значений ~„. В общем случае края на изображении создают пространственные частоты вдоль линии на плоскости комплексных частот, перпендикулярной краю, Чем резче край, тем дальше от начала координат мы должны уйти по этой линии, чтобы величина весовых коэффициен- 323 8.2. Теорема ошомаов тов стала незначительной. Таким образом, наш-качественный вывод заключается в том, что высокие пространственные частоты соответствуют резким краям, низкие — отсутствию краев (т.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
