МЕТОДИКА ОБРАБОТКИ ЛИНЕЙЧАТЫХ ОБРАЗОВ НА ДЕФЕКТОСКОПИЧЕСКИХ СНИМКАХ
Методика обработки линейчатых образов на дефектоскопических снимках
11
УДК 004.936
С. С. САДЫКОВ, А. А. ОРЛОВ, А. А. ЕРМАКОВ
МЕТОДИКА ОБРАБОТКИ ЛИНЕЙЧАТЫХ ОБРАЗОВ НА ДЕФЕКТОСКОПИЧЕСКИХ СНИМКАХ
Дается определение полосового (линейчатого) образа и строится его модель. Разрабатываются метод и алгоритм многомасштабной фильтрации полосовых образов. Представлена методика обработки дефектоскопических снимков, содержащих линейчатые образы дефектов. Приводится пример обработки реального изображения.
Ключевые слова: пространственно-распределенный объект, топологическое отношение, муниципальные геоинформационные системы, матрица шестнадцати пересечений.
В настоящее время прослеживается тенденция развития систем технического зрения и методов обработки изображений в цифровой дефектоскопии. Цифровой анализ позволяет быстро обнаруживать, измерять и классифицировать характер дефектов для обеспечения неразрушающего контроля различных материалов и изделий. Существует множество публикаций, связанных с анализом дефектоскопических снимков различной природы путем компьютерной обработки изображений. Основные подходы к решению таких задач представлены в статьях [1, 2].
Многие дефекты имеют специфические сложные образы. Например, очень распространены дефекты (трещины, риски, царапины), имеющие протяженную линейчатую форму (полосы). В существующих системах технического зрения, применяемых в дефектоскопии, в основном используются классические методы цифровой обработки изображений [3], не позволяющие выполнить более детальный разбор сложных по форме дефектов. Поэтому актуальным является создание, исследование и применение новых, более эффективных, специальных методов обработки изображений полосовых дефектов.
Из теории сигналов известен так называемый фильтр воспроизведения полезного сигнала [3]. Этот фильтр позволяет выделить из шума и помех сигнал заданной формы. Однако такой фильтр используется в основном для обработки одномерных сигналов. В нашем случае
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 2
12 С. С. Садыков, А. А. Орлов, А. А. Ермаков
к исходным данным для анализа относится изображение, которое рассматривается как двумерный сигнал. Предполагается использовать специфику формы образов — их полосовую структуру (производить фильтрацию вдоль полос).
Цель настоящей работы заключается в разработке специальной технологии обработки дефектоскопических снимков, содержащих линейчатые образы дефектов.
Математическая модель полосового образа. Пусть x = x(t), y = y(t) — параметрические уравнения некоторой кривой (рис. 1, пунктирная кривая) на плоскости R2. Множество точек, расположенных вдоль кривой x = x(t), y = y(t) на расстоянии, не превышающем значения σ,
S = {(x0 , y0 ) (x0 − x(t))2 + ( y0 − y(t))2 ≤ σ2} ,
назовем областью полосы, заданной этой кривой. Кривую x = x(t), y = y(t) будем называть образующей полосы S.
y
2l
θ 2σ
x
Рис. 1
Если значение σ (полуширина сегмента) является функцией σ = σ(t), то полосу S назовем полосой с переменной шириной, так что 2σ(t) — ширина этой полосы, изменяющаяся в соответствии с образующей кривой x = x(t), y = y(t).
Сцена (изображение) полосы задается характеристической функцией
χ(x0 ,
y0 )
=
⎧1, ⎩⎨0.
(x0 ,
y0 ) ∈
S,
Профилем полосы S в точке (x, y) будем называть сечение поверхности z = χ(x, y) нормальной плоскостью в пространстве R3 к образующей кривой x = x(t), y = y(t) в точке (x, y). Профиль изображения полосы S будет иметь прямоугольную форму.
Обобщим понятие образа полосы, видоизменив ее профиль. Для этого будем полагать, что полосовое изображение состоит (точнее, интегрируется) из бесконечного множества профилей, расположенных вдоль кривой x = x(t), y = y(t).
Пусть A(x) — профиль полосы, так что
lim A(x) = 0 ,
x→±∞
тогда характеристическая функция изображения полосы профиля A(x) определится как криволинейный интеграл:
χ(x0 , y0 ) = ∫ A(| (x, y) − (x0 , y0 ) |) ds . x= x (t ), y= y(t)
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 2
Методика обработки линейчатых образов на дефектоскопических снимках
13
Примеры изображения кривой (а), полосы с прямоугольным профилем (б) и полосы, где A(x) — гауссиан (в), приведены на рис. 2.
а) б)
в)
Рис. 2
Прямоугольник, аппроксимирующий часть полосы, будем называть сегментом полосы.
Обозначим полудлину сегмента как l, а угол наклона сегмента (равный углу наклона нормали
к образующей полосы) — θ (см. рис. 1).
Воспроизведение одномерных сигналов. Рассмотрим одномерный случай. Фурье-
спектр фильтра воспроизведения полезного сигнала определяется как
Ψ
(ω)
=
Ws
Ws (ω) (ω) + Wn
(ω)
,
где Ws(ω) и Wn(ω) — энергетический спектр (плотность мощности) полезного сигнала и помех. Под помехами (шумами) будем понимать все сигналы вида, отличного от s(x).
Энергетические спектры определяются как
Ws(ω) = |S(ω)|2 и Wn(ω) = |N(ω)|2,
где S(ω), N(ω) — фурье-спектры полезного сигнала s(x) и помех n(x) соответственно. Полезный сигнал (в нашем случае это функция s(x)) содержится в исходном сигнале f(x)
в сумме с шумами и помехами n(x): f(x) = s(x) + n(x).
Пусть ψ(x) — фильтр воспроизведения полезного сигнала, g(x) — функция восстанов-
ленного сигнала (g(x) = f(x)⊗ψ(x)). С помощью фильтра ψ(x) производится минимизация средней квадратической ошибки, которая характеризует степень отклонения функции восстановленного сигнала g(x) от полезного s(x):
ε = M[(g(x) − s(x))2 ] = M[( f (x) ⊗ ψ(x) − s(x))2 ] → min ,
где M — оператор среднего.
Пусть теперь s(x) — многомасштабная функция, т.е. s = s(x/σ), так что σ ∈[σ1, σ2] (σ1, σ2 — минимальное и максимальное значения масштаба соответственно). Для фильтрации сигналов различного масштаба в равной степени необходимо, чтобы средняя квадратическая
ошибка была одинаковой для каждого значения σ из [σ1, σ2]. В качестве критерия оценки точности функционирования фильтра примем среднее значение ошибки для различных масштабов:
ε = Mσ{M[( f (x) ⊗ ψ(x) − s(x / σ))2 ]} → min ,
где Mσ — оператор среднего по масштабу σ. Можно доказать, что такой критерий обеспечивается фильтром при
Ψ(ω)
=
Mσ[Ws (σω)] Mσ[Ws (σω)] + Wn
(ω)
,
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 2
14 С. С. Садыков, А. А. Орлов, А. А. Ермаков
так что
∫M σ [Ws
(σω)]
=
σ2
1 −
σ1
σ2
Ws
σ1
(σω)d σ
.
Двумерная фильтрация вдоль полосы. Функцию s(x, y), которая максимально соот-
ветствует изображению полосы χ(x, y) в каждой точке области ее сегмента, назовем детектором сегмента полосы (или детектором полосы).
В данном случае детектор сегмента полосы будет определять полезный двумерный сиг-
нал. Можно показать, что если полоса имеет прямоугольный профиль, то
s(x, y)
=
⎪⎧e− y2 ⎨
2,
−1 ≤ x ≤ 1,
⎪⎩0.
Пусть S(ω, ν), N(ω, ν) — фурье-спектры детектора сегмента полосы и функции шума
соответственно.
Назовем воспроизводящим интегральным преобразованием по сегменту полосы
(ВИПСП) преобразование Hθ[⋅], выполняющее фильтрацию в пределах [σ1, σ2] функции f(x, y) вдоль каждого направления θ. Спектральная функция g(x, y, θ) формируется по правилу
G(ω, ν, θ)
=
F (ω, ν)
Ms
+
Ms N (ω, ν)
2
,
где F(ω, ν), F(ω, ν, θ) — фурье-спектры функций f(x, y) и g(x, y, θ) соответственно,
∫M s
=
σ2
1 − σ1
σ2 σ1
2
S(Rot(σω, lν))
θ
dσ .
Преобразование исходного изображения f(x, y) на основе ВИПСП в изображение g(x, y),
при котором выполняется условие
g(x, y) = max Hθ[f(x, y)],
θ
назовем воспроизведением полосы (ВП).
С помощью ВП можно найти максимумы в откликах ВИПСП, тем самым среди сегмен-
тов разной ориентации определив сегмент, который наилучшим образом соответствует части
полосы.
Приведем алгоритм воспроизведения полосы:
1) формирование фурье-спектра F(ω, ν) исходного изображения f(x, y);
2) вычисление спектра шума N(ω, ν) по области изображения, которая не содержит по-
лезного сигнала;
3) выполнение воспроизводящего интегрального преобразования
G(ω, ν, θ) = Hθ[ F(ω, ν)]; 4) обратное преобразование фурье-спектра G(ω, ν, θ) в g(x, y, θ);
5) поиск максимумов в функции g(x, y, θ)
g(x, y) = max g(x, y, θ).
θ
Методика обработки дефектоскопических изображений. Методика обработки и ана-
лиза дефектоскопических изображений, содержащих полосовые образы, включает следую-
щие этапы:
1) воспроизведение полосовых образов с помощью разработанного алгоритма;
2) бинаризация фильтрованного изображения;
3) вычисление площади выделенных областей;
4) устранение малых (незначащих) областей;
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 2
Методика обработки линейчатых образов на дефектоскопических снимках
15
5) определение признаков (характеристик) выделенных дефектных областей; 6) классификация дефектов. Фильтрация позволяет выделить образы протяженной формы, которые, вероятнее всего, являются дефектными, и подавить фоновые шумы и помехи, вызванные неоднородной структурой оцениваемой поверхности. Пороговая бинаризация изображения необходима для отделения дефектных областей от фона. В результате бинаризации возможно выделение малых зашумленных областей, которые не являются дефектными. Для устранения из рассмотрения таких областей выполняется их отсечение по площади. Порог бинаризации необходимо выбирать в зависимости от уровня значимости обнаруживаемых дефектов. Этап вычисления геометрических признаков объектов основан на использовании известных алгоритмов обработки монохромных (бинарных) изображений, вычисляются их площадь, длина, ширина и др. На этапе классификации образов выполняются распознавание и классификация областей. В зависимости от их ширины определяется значимость дефекта. В результате получается характеристическое (маркированное) изображение с отмеченными областями значимых дефектов.
Экспериментальные результаты. На рис. 3 приведен пример воспроизведения образа трещины как полосы прямоугольного профиля на дефектоскопическом изображении (здесь а — исходное изображение, б — результат фильтрации; в, г — их бинарные образы соответственно). Обработка выполнена при параметрах, которые соответствуют выделяемому объек-
ту: угол наклона сегмента полосы θ ∈ [0, π), полуширина полосы σ ∈ [0, 5, 3], полудлина сег-
мента полосы l = 3. Спектр помех N(ω, ν) определен по области изображения, которая не содержит трещины.
а) б)
в) г)
Рис. 3
Проведена бинаризация исходного и обработанного изображений. Порог бинаризации выбран таким образом, чтобы образ трещины не содержал разрывов. Из рис. 3, в, г видно, что с помощью фильтрации возможно довольно четко отделить фон от изображения, следовательно, возможно выделять изображения дефектов такого рода для дальнейшей автоматической оценки качества поверхности анализируемого изделия.
Заключение. Таким образом, в настоящей работе построена математическая модель обработки дефектоскопических изображений, учитывающая множество структурных особенностей полосовых объектов. Это обеспечивает качественный результат обработки изображений,
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 2
16 В. В. Гордиенко, В. М. Довгаль, Р. А. Пузына
который подтверждается экспериментальными исследованиями. Обработка образов полос на основе предлагаемой методики обеспечивает их непрерывность на малоконтрастных изображениях. Появляется возможность дальнейшего автоматического анализа.
Следует отметить возможность и целесообразность применения метода в промышленности для автоматического анализа различного рода сцен, которые содержат полосовые образы.
СПИСОК ЛИТЕРАТУРЫ
1. Mery D., Berti M. A. Automatic Detection of Welding Defects Using Texture Features // Insight. 2003. Vol. 45, N 10. P. 676—681.
2. Silva R. R., Caloba L. P., Siqueira M. H. S. Evaluation of the Relevant Characteristic Parameters of Welding Defects and Probability of Correct Classification Using Linear Classifiers // Insight. 2002. Vol. 44, N 10. P. 197—218.
3. Методы компьютерной обработки изображений / Под ред. В. А. Сойфера. М.: Физматлит, 2003. 784 с.
Султан Садыкович Садыков Алексей Александрович Орлов
Александр Андреевич Ермаков
Сведения об авторах — д-р техн. наук, профессор; Муромский институт (филиал) Владимирско-
го государственного университета, кафедра информационных систем — канд. техн. наук, доцент; Муромский институт (филиал) Владимир-
ского государственного университета, кафедра информационных систем; E-mail: AlexeyAlexOrlov@rambler.ru — аспирант; Муромский институт (филиал) Владимирского государственного университета, кафедра информационных систем
Рекомендована кафедрой информационных систем
Поступила в редакцию 12.09.08 г.
УДК 681.3
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 2
11
УДК 004.936
С. С. САДЫКОВ, А. А. ОРЛОВ, А. А. ЕРМАКОВ
МЕТОДИКА ОБРАБОТКИ ЛИНЕЙЧАТЫХ ОБРАЗОВ НА ДЕФЕКТОСКОПИЧЕСКИХ СНИМКАХ
Дается определение полосового (линейчатого) образа и строится его модель. Разрабатываются метод и алгоритм многомасштабной фильтрации полосовых образов. Представлена методика обработки дефектоскопических снимков, содержащих линейчатые образы дефектов. Приводится пример обработки реального изображения.
Ключевые слова: пространственно-распределенный объект, топологическое отношение, муниципальные геоинформационные системы, матрица шестнадцати пересечений.
В настоящее время прослеживается тенденция развития систем технического зрения и методов обработки изображений в цифровой дефектоскопии. Цифровой анализ позволяет быстро обнаруживать, измерять и классифицировать характер дефектов для обеспечения неразрушающего контроля различных материалов и изделий. Существует множество публикаций, связанных с анализом дефектоскопических снимков различной природы путем компьютерной обработки изображений. Основные подходы к решению таких задач представлены в статьях [1, 2].
Многие дефекты имеют специфические сложные образы. Например, очень распространены дефекты (трещины, риски, царапины), имеющие протяженную линейчатую форму (полосы). В существующих системах технического зрения, применяемых в дефектоскопии, в основном используются классические методы цифровой обработки изображений [3], не позволяющие выполнить более детальный разбор сложных по форме дефектов. Поэтому актуальным является создание, исследование и применение новых, более эффективных, специальных методов обработки изображений полосовых дефектов.
Из теории сигналов известен так называемый фильтр воспроизведения полезного сигнала [3]. Этот фильтр позволяет выделить из шума и помех сигнал заданной формы. Однако такой фильтр используется в основном для обработки одномерных сигналов. В нашем случае
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 2
12 С. С. Садыков, А. А. Орлов, А. А. Ермаков
к исходным данным для анализа относится изображение, которое рассматривается как двумерный сигнал. Предполагается использовать специфику формы образов — их полосовую структуру (производить фильтрацию вдоль полос).
Цель настоящей работы заключается в разработке специальной технологии обработки дефектоскопических снимков, содержащих линейчатые образы дефектов.
Математическая модель полосового образа. Пусть x = x(t), y = y(t) — параметрические уравнения некоторой кривой (рис. 1, пунктирная кривая) на плоскости R2. Множество точек, расположенных вдоль кривой x = x(t), y = y(t) на расстоянии, не превышающем значения σ,
S = {(x0 , y0 ) (x0 − x(t))2 + ( y0 − y(t))2 ≤ σ2} ,
назовем областью полосы, заданной этой кривой. Кривую x = x(t), y = y(t) будем называть образующей полосы S.
y
2l
θ 2σ
x
Рис. 1
Если значение σ (полуширина сегмента) является функцией σ = σ(t), то полосу S назовем полосой с переменной шириной, так что 2σ(t) — ширина этой полосы, изменяющаяся в соответствии с образующей кривой x = x(t), y = y(t).
Сцена (изображение) полосы задается характеристической функцией
χ(x0 ,
y0 )
=
⎧1, ⎩⎨0.
(x0 ,
y0 ) ∈
S,
Профилем полосы S в точке (x, y) будем называть сечение поверхности z = χ(x, y) нормальной плоскостью в пространстве R3 к образующей кривой x = x(t), y = y(t) в точке (x, y). Профиль изображения полосы S будет иметь прямоугольную форму.
Обобщим понятие образа полосы, видоизменив ее профиль. Для этого будем полагать, что полосовое изображение состоит (точнее, интегрируется) из бесконечного множества профилей, расположенных вдоль кривой x = x(t), y = y(t).
Пусть A(x) — профиль полосы, так что
lim A(x) = 0 ,
x→±∞
тогда характеристическая функция изображения полосы профиля A(x) определится как криволинейный интеграл:
χ(x0 , y0 ) = ∫ A(| (x, y) − (x0 , y0 ) |) ds . x= x (t ), y= y(t)
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 2
Методика обработки линейчатых образов на дефектоскопических снимках
13
Примеры изображения кривой (а), полосы с прямоугольным профилем (б) и полосы, где A(x) — гауссиан (в), приведены на рис. 2.
а) б)
в)
Рис. 2
Прямоугольник, аппроксимирующий часть полосы, будем называть сегментом полосы.
Обозначим полудлину сегмента как l, а угол наклона сегмента (равный углу наклона нормали
к образующей полосы) — θ (см. рис. 1).
Воспроизведение одномерных сигналов. Рассмотрим одномерный случай. Фурье-
спектр фильтра воспроизведения полезного сигнала определяется как
Ψ
(ω)
=
Ws
Ws (ω) (ω) + Wn
(ω)
,
где Ws(ω) и Wn(ω) — энергетический спектр (плотность мощности) полезного сигнала и помех. Под помехами (шумами) будем понимать все сигналы вида, отличного от s(x).
Энергетические спектры определяются как
Ws(ω) = |S(ω)|2 и Wn(ω) = |N(ω)|2,
где S(ω), N(ω) — фурье-спектры полезного сигнала s(x) и помех n(x) соответственно. Полезный сигнал (в нашем случае это функция s(x)) содержится в исходном сигнале f(x)
в сумме с шумами и помехами n(x): f(x) = s(x) + n(x).
Пусть ψ(x) — фильтр воспроизведения полезного сигнала, g(x) — функция восстанов-
ленного сигнала (g(x) = f(x)⊗ψ(x)). С помощью фильтра ψ(x) производится минимизация средней квадратической ошибки, которая характеризует степень отклонения функции восстановленного сигнала g(x) от полезного s(x):
ε = M[(g(x) − s(x))2 ] = M[( f (x) ⊗ ψ(x) − s(x))2 ] → min ,
где M — оператор среднего.
Пусть теперь s(x) — многомасштабная функция, т.е. s = s(x/σ), так что σ ∈[σ1, σ2] (σ1, σ2 — минимальное и максимальное значения масштаба соответственно). Для фильтрации сигналов различного масштаба в равной степени необходимо, чтобы средняя квадратическая
ошибка была одинаковой для каждого значения σ из [σ1, σ2]. В качестве критерия оценки точности функционирования фильтра примем среднее значение ошибки для различных масштабов:
ε = Mσ{M[( f (x) ⊗ ψ(x) − s(x / σ))2 ]} → min ,
где Mσ — оператор среднего по масштабу σ. Можно доказать, что такой критерий обеспечивается фильтром при
Ψ(ω)
=
Mσ[Ws (σω)] Mσ[Ws (σω)] + Wn
(ω)
,
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 2
14 С. С. Садыков, А. А. Орлов, А. А. Ермаков
так что
∫M σ [Ws
(σω)]
=
σ2
1 −
σ1
σ2
Ws
σ1
(σω)d σ
.
Двумерная фильтрация вдоль полосы. Функцию s(x, y), которая максимально соот-
ветствует изображению полосы χ(x, y) в каждой точке области ее сегмента, назовем детектором сегмента полосы (или детектором полосы).
В данном случае детектор сегмента полосы будет определять полезный двумерный сиг-
нал. Можно показать, что если полоса имеет прямоугольный профиль, то
s(x, y)
=
⎪⎧e− y2 ⎨
2,
−1 ≤ x ≤ 1,
⎪⎩0.
Пусть S(ω, ν), N(ω, ν) — фурье-спектры детектора сегмента полосы и функции шума
соответственно.
Назовем воспроизводящим интегральным преобразованием по сегменту полосы
(ВИПСП) преобразование Hθ[⋅], выполняющее фильтрацию в пределах [σ1, σ2] функции f(x, y) вдоль каждого направления θ. Спектральная функция g(x, y, θ) формируется по правилу
G(ω, ν, θ)
=
F (ω, ν)
Ms
+
Ms N (ω, ν)
2
,
где F(ω, ν), F(ω, ν, θ) — фурье-спектры функций f(x, y) и g(x, y, θ) соответственно,
∫M s
=
σ2
1 − σ1
σ2 σ1
2
S(Rot(σω, lν))
θ
dσ .
Преобразование исходного изображения f(x, y) на основе ВИПСП в изображение g(x, y),
при котором выполняется условие
g(x, y) = max Hθ[f(x, y)],
θ
назовем воспроизведением полосы (ВП).
С помощью ВП можно найти максимумы в откликах ВИПСП, тем самым среди сегмен-
тов разной ориентации определив сегмент, который наилучшим образом соответствует части
полосы.
Приведем алгоритм воспроизведения полосы:
1) формирование фурье-спектра F(ω, ν) исходного изображения f(x, y);
2) вычисление спектра шума N(ω, ν) по области изображения, которая не содержит по-
лезного сигнала;
3) выполнение воспроизводящего интегрального преобразования
G(ω, ν, θ) = Hθ[ F(ω, ν)]; 4) обратное преобразование фурье-спектра G(ω, ν, θ) в g(x, y, θ);
5) поиск максимумов в функции g(x, y, θ)
g(x, y) = max g(x, y, θ).
θ
Методика обработки дефектоскопических изображений. Методика обработки и ана-
лиза дефектоскопических изображений, содержащих полосовые образы, включает следую-
щие этапы:
1) воспроизведение полосовых образов с помощью разработанного алгоритма;
2) бинаризация фильтрованного изображения;
3) вычисление площади выделенных областей;
4) устранение малых (незначащих) областей;
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 2
Методика обработки линейчатых образов на дефектоскопических снимках
15
5) определение признаков (характеристик) выделенных дефектных областей; 6) классификация дефектов. Фильтрация позволяет выделить образы протяженной формы, которые, вероятнее всего, являются дефектными, и подавить фоновые шумы и помехи, вызванные неоднородной структурой оцениваемой поверхности. Пороговая бинаризация изображения необходима для отделения дефектных областей от фона. В результате бинаризации возможно выделение малых зашумленных областей, которые не являются дефектными. Для устранения из рассмотрения таких областей выполняется их отсечение по площади. Порог бинаризации необходимо выбирать в зависимости от уровня значимости обнаруживаемых дефектов. Этап вычисления геометрических признаков объектов основан на использовании известных алгоритмов обработки монохромных (бинарных) изображений, вычисляются их площадь, длина, ширина и др. На этапе классификации образов выполняются распознавание и классификация областей. В зависимости от их ширины определяется значимость дефекта. В результате получается характеристическое (маркированное) изображение с отмеченными областями значимых дефектов.
Экспериментальные результаты. На рис. 3 приведен пример воспроизведения образа трещины как полосы прямоугольного профиля на дефектоскопическом изображении (здесь а — исходное изображение, б — результат фильтрации; в, г — их бинарные образы соответственно). Обработка выполнена при параметрах, которые соответствуют выделяемому объек-
ту: угол наклона сегмента полосы θ ∈ [0, π), полуширина полосы σ ∈ [0, 5, 3], полудлина сег-
мента полосы l = 3. Спектр помех N(ω, ν) определен по области изображения, которая не содержит трещины.
а) б)
в) г)
Рис. 3
Проведена бинаризация исходного и обработанного изображений. Порог бинаризации выбран таким образом, чтобы образ трещины не содержал разрывов. Из рис. 3, в, г видно, что с помощью фильтрации возможно довольно четко отделить фон от изображения, следовательно, возможно выделять изображения дефектов такого рода для дальнейшей автоматической оценки качества поверхности анализируемого изделия.
Заключение. Таким образом, в настоящей работе построена математическая модель обработки дефектоскопических изображений, учитывающая множество структурных особенностей полосовых объектов. Это обеспечивает качественный результат обработки изображений,
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 2
16 В. В. Гордиенко, В. М. Довгаль, Р. А. Пузына
который подтверждается экспериментальными исследованиями. Обработка образов полос на основе предлагаемой методики обеспечивает их непрерывность на малоконтрастных изображениях. Появляется возможность дальнейшего автоматического анализа.
Следует отметить возможность и целесообразность применения метода в промышленности для автоматического анализа различного рода сцен, которые содержат полосовые образы.
СПИСОК ЛИТЕРАТУРЫ
1. Mery D., Berti M. A. Automatic Detection of Welding Defects Using Texture Features // Insight. 2003. Vol. 45, N 10. P. 676—681.
2. Silva R. R., Caloba L. P., Siqueira M. H. S. Evaluation of the Relevant Characteristic Parameters of Welding Defects and Probability of Correct Classification Using Linear Classifiers // Insight. 2002. Vol. 44, N 10. P. 197—218.
3. Методы компьютерной обработки изображений / Под ред. В. А. Сойфера. М.: Физматлит, 2003. 784 с.
Султан Садыкович Садыков Алексей Александрович Орлов
Александр Андреевич Ермаков
Сведения об авторах — д-р техн. наук, профессор; Муромский институт (филиал) Владимирско-
го государственного университета, кафедра информационных систем — канд. техн. наук, доцент; Муромский институт (филиал) Владимир-
ского государственного университета, кафедра информационных систем; E-mail: AlexeyAlexOrlov@rambler.ru — аспирант; Муромский институт (филиал) Владимирского государственного университета, кафедра информационных систем
Рекомендована кафедрой информационных систем
Поступила в редакцию 12.09.08 г.
УДК 681.3
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2009. Т. 52, № 2