Например, Бобцов

МЕТОДЫ ВЫДЕЛЕНИЯ И ВЫЧИСЛЕНИЯ ПРИЗНАКОВ СТРУКТУР НА СНИМКАХ ПРОМЫШЛЕННЫХ ИЗДЕЛИЙ

10
УДК 004.936
А. А. ОРЛОВ
МЕТОДЫ ВЫДЕЛЕНИЯ И ВЫЧИСЛЕНИЯ ПРИЗНАКОВ СТРУКТУР НА СНИМКАХ ПРОМЫШЛЕННЫХ ИЗДЕЛИЙ
Показана актуальность повышения оперативности и достоверности результатов обработки снимков промышленных изделий. Разрабатывается система методов анализа таких снимков, основанного на интегральных преобразованиях по сегментам линейчатых структур. Ключевые слова: автоматический анализ снимков изделий, определение дефектов, обработка изображений.
Введение. В настоящее время повышение качества продукции машиностроительных и металлообрабатывающих предприятий за счет внедрения новых технологий (в частности, технологий дефектоскопии) является очень актуальным. Необходимо повысить оперативность и достоверность процедур определения дефектов промышленных изделий по снимкам их микро- и наноструктур, но отсутствуют современные методы и алгоритмы их выявления.
На сегодняшний день в промышленности ручной контроль качества все еще преобладает над автоматизированным. Главная причина — отсутствие более совершенных методов автоматической обработки снимков (изображений), получаемых в большинстве способов контроля качества. В настоящее время для визуального контроля качества стали применяться автоматизированные системы анализа изображений (АСАИ).
Большой вклад в разработку систем и методов компьютерного анализа снимков промышленных изделий внесли Д. Мери, Р. Силва, Н. Насреддин, К. Демандт, Р. Гонсалес, И. Като, Г. Падуа, Х. Шафик, В. В. Клюев, М. В. Филинов, В. А. Лопухин и др. При обработке большинства изображений, как правило, требуется выделять и анализировать их линейчатые структуры, которыми насыщены микро- и наноснимки промышленных изделий, в частности, это дефекты-трещины, риски, царапины, швы, контуры, волокна, линейчатые вкрапления и др. Для оценки линейчатых образов на снимках измеряются геометрические характеристики.
ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 2

Методы выделения и вычисления признаков структур на снимках промышленных изделий 11

Линейчатые структуры могут быть локально прямолинейными и представленными в виде по-

лос с известным профилем и варьируемой шириной. При этом они часто довольно сильно

зашумлены, размыты, малоконтрастны, насыщены пересекающимися образами подобных или

других типов. Анализ таких структур на многих изображениях довольно сложен. Сложность

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

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

вием шумов и размытий, появившихся при формировании снимка, а также образов, не значи-

мых для исследования объектов и их теней.

АСАИ состоят из цифрового микроскопа и компьютерной системы, реализующей мето-

ды обработки и анализа изображений. Однако в них применяются довольно простые методы,

которые не позволяют выполнить более детальный разбор линейчатых по форме дефектов и с

заданной точностью определить геометрические характеристики дефектов. В связи с этим

уровень автоматической обработки и анализа дефектоскопических изображений очень низок.

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

контроля качества. Вклад в разработку методов выделения и анализа изображений линейча-

тых структур внесли Д. Канни, А. Рад, M. Хюккель, П. Хаф, Ф. Кларк, Д. Пао, K. Канатани,

С. Наяр, С. Кониши, С. В. Абламейко, В. А. Сойфер, В. В. Сергеев, Н. Г. Федотов, Я. А. Фур-

ман и многие другие. Однако существующие методы не позволяют выделять линейчатую

структуру, характеризующуюся вариациями кривизны, яркости, ширины, пересечениями и

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

металлографических снимков. Поэтому на практике для выделения линейчатой структуры

часто применяются следующие методы сегментации: пороговая, основанная на областях или

морфологическом анализе. Однако данные методы ориентированы на выделение площадных

образов и применимы лишь для обработки простых изображений линейчатых структур. Соз-

дание новых методов, основанных на анализе линейчатых структур, в рамках настоящего

проекта позволит повысить оперативность и достоверность определения дефектов промыш-

ленных изделий по их снимкам. Целью настоящей статьи является разработка системы таких

методов.

Выделение линейчатых структур с варьирующейся кривизной. Введем следующее

интегральное преобразование:

∫ ∫h(x0, y0, θ, σ) =

1 sL

2

∞∞ −∞ −∞

f

( x,

y)sL (x



x0 ,

y



y0, θ, σ)dxdy

,

где f(x,y) — исходное изображение; h(x0,y0,θ,σ) — спектральная функция; sL(x,y,θ,σ) — детектор сегмента линейчатой структуры; σ — масштабный коэффициент; θ — угол поворота.
Для анализа образов с варьирующейся кривизной разработан метод выделения и вычис-

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

( )sL (x, y) =

1− x2

e−

x2

+ 2

y

2

.

Детектор линейчатой структуры в нашем случае будет задаваться функцией трех пере-

менных:

sL(x,

y,θ)

=

sL

⎢⎣⎡Tθ

⎛ ⎝⎜

x σ

,

y σ

⎞⎤ ⎠⎟⎥⎦

= ⎜⎝⎛⎜1−

(x cos θ − y σ2

sin

θ)2

⎞ ⎠⎟⎟


e

x2+ y 2σ2

2

,

где

T
θ

=

⎛ cos θ

⎜ ⎝

sin

θ

− sin θ⎞

cos θ

⎟ ⎠



оператор

(матрица) поворота

на

угол

θ.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 2

12 А. А. Орлов

Функция sL(x,y,θ) может быть разложена на три составляющие:

sL (x, y, θ) = s1(x, y) cos 2θ + s2 (x, y)sin 2θ + s3 (x, y) ,

где s1(x,y), s2(x,y) — нормированные гауссоидом гиперболические параболоиды:

s1(x, y)

=



x2 − y2 2σ2


e

x

2 + y2 2σ2

,

s2 (x, y)

=

xy σ2

e



x2+ y2 2σ2

,

а s3(x,y) — лапласиан гауссоида:

s3 (x,

y)

=

⎝⎜⎜⎛1 −

x2 + y2 2σ2

⎞ ⎟⎟⎠


e

x2+ y 2σ2

2

.

Спектральная функция h(x0,y0,θ) принимает максимальное значение в точке

θ0

=

1 2

<

(S1, S2 )

,

так что

h(θ0 ) = S12 + S22 + S3 ,
где S1, S2, S3 — свертки исходного изображения f(x,y) с функциями s1(x,y), s2(x,y) и s3(x,y) соответственно:

∞∞
∫ ∫Si = f (x, y)si (x − x0 , y − y0 )dxdy, i = 1, 2, 3. −∞ −∞
Преобразование исходного изображения f(x, y) в изображение g(x, y), так что

g(x, y) = max h(x,y,θ) = ( f ⊗ s1)2 + ( f ⊗ s2 )2 + f ⊗ s3 ,
θ

назовем выделением линейчатых структур без вращения их сегментов.

Метод выделения и вычисления признаков структур без вращения их сегментов пред-

полагает вычисление g(x, y) и θ(x, y) по данным формулам.

Выделение линейчатых структур с изменяющейся шириной. Для анализа образов с

изменяющейся шириной разработан метод выделения линейчатых структур различного мас-

штаба.

Назовем интегральным преобразованием по сегменту линейчатой структуры различного

масштаба преобразование Hθ, выполняющее фильтрацию функции f(x, y) по масштабу в пределах [σ1, σ2] вдоль каждого направления θ. При этом формируется спектральная функция

g(x, y, θ), такая что

( )∫G(ω,ν,θ)

=

F (ω, ν) Cs

σ2 σ1

SL

T⋅ (σω, ν)
θ

2 dσ , σ

Cs — нормирующий параметр, SL(ω, ν), G(ω, ν, θ), F(ω, ν) — фурье-спектры функций sL(x, y), g(x, y, θ) и f(x, y) соответственно.

Преобразование исходного изображения f(x, y) в изображение g(x, y), так что

g(x, y) = max Hθ[f(x, y)],
θ

назовем выделением линейчатых структур различного масштаба.

Метод выделения линейчатых структур различного масштаба предполагает вычисление

G(ω,ν,θ) и g(x, y) по данным формулам, при этом

σ1 ≤ σ ≤ σ2. Преобразование заключается в поиске угла наклона детектора, который соответствует

максимуму в спектре g(x, y, θ).

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 2

Методы выделения и вычисления признаков структур на снимках промышленных изделий 13

( )Если sL (x, y) =

1− x2




e

x2+ y2 2

,

то

G(ω, ν, θ)

=

1 Cs

σ2

F

(ω,

ν)

ψ

L

[T(σω,
θ

lν)]

σ1

,

где ψL (ω, ν) = (1 − ω2 )e−ω2−ν2 .

Воспроизведение линейчатых структур. Для выделения образов на снимках без иска-

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

многомасштабной фильтрации профилей линейчатых структур. Для осуществления много-

масштабной фильтрации профилей линейчатых структур использованы два преобразования:

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

разование для воспроизведения линейчатой структуры, подробно описанные в работе [1].

Метод воспроизведения линейчатых структур по максимуму в откликах воспроизводя-

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

повторить часть образа на снимке.

Вычисление кривизны линейчатых структур. Для анализа линейчатых структур раз-

личной формы разработан метод вычисления их кривизны образующих.

Значения координат образующей (x, y), угол наклона нормали к ней θ и масштаб σ воз-

можно найти, определив максимумы спектральной функции, т.е. в результате решения задачи

оптимизации:

h(x, y, θ, σ) → max,

0 ≤ θ < 2π,

σ1 ≤ σ ≤ σ2. В статье [2] показано, что кривизна линии яркостного перепада в точке (x, y) на плоско-

сти xOy определяется как абсолютное значение дивергенции векторного поля градиента яр-

кости в этой точке:

k(x, y) = div (cosθ,sinθ)

для дискретного пространства

( )k(i,

j)



1 2

Xi+1, j − X ij + Yi, j+1 − Yij + Xi−1, j − X ij + Yi, j−1 − Yij

,

где X, Y — матрицы значений составляющих нормали к линейчатой структуре по оси x и y

соответственно.

Заключение. Рассмотренные методы реализованы и апробированы на реальных изо-

бражениях промышленных объектов. Результаты экспериментальных исследований под-

тверждают правильность теоретических выводов [1—4]. Таким образом, разработана система

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

вариацией кривизны, яркости, ширины, пересечениями, разветвлениями, наблюдающиеся на

снимках промышленных изделий.

СПИСОК ЛИТЕРАТУРЫ
1. Садыков С. С., Орлов А. А., Ермаков А. А. Теория, алгоритмы и методика обработки линейчатых образов на дефектоскопических снимках // Изв. вузов. Приборостроение. 2009. Т. 52, № 2. С. 11—16.
2. Орлов А. А. Реализация системы обработки изображений линейчатых объектов // Программные продукты и системы. 2007. № 4. С. 61—62.
3. Орлов А. А. Компьютерный рентгенографический анализ качества сварных соединений и выделение линейчатых объектов в них // Автоматизация и современные технологии. 2009. № 6. C. 5—7.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 2

14 С. С. Садыков, С. В. Савичева

4. Орлов А. А., Ермаков А. А. Фильтрация полосовых образов прямоугольного профиля // Науч.-техн. ведомости СПбГПУ. Информатика. Телекоммуникации. Управление. 2008. № 3(60). С. 52—56.

Алексей Александрович Орлов

Сведения об авторе — д-р техн. наук, доцент; Муромский институт Владимирского государ-
ственного университета им. А. Г. и Н. Г. Столетовых, кафедра физики и прикладной математики; E-mail: AlexeyAlexOrlov@rambler.ru

Рекомендована Юго-Западным государственным университетом

Поступила в редакцию 24.10.11 г.

ИЗВ. ВУЗОВ. ПРИБОРОСТРОЕНИЕ. 2012. Т. 55, № 2