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

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

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

9.3. Распознавание дефектов сварных соединений системой

201

В этом случае для принятия решений о принадлежности распо­ знаваемого объекта к определённому классу целесообразно использо­ вать одну из разновидностей логического решающего правила CORAL [14, 19].

На основании данного правила построены логические высказыва­ ния S, отражающие факт попадания распознаваемого объекта а в об­ ласть I/, образованную границами подмножеств X iv, т. е. в гиперпарал­ лелепипед

Х\и х Х2и х . . . х Х п^ = J(a, X\jy^ X J(a, X%i^ X ... X */(ц, A )^),

где J(ai,Xj„) означает, что значение признака X j объекта а* принад­ лежит подмножеству А*„.

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

Принятие решения о принадлежности изображения одному из клас­ сов в распознающей системе происходит по информации о числе дефек­ тов w:

если w = 0, то изображение не содержит дефектов;

если w = 1, изображение содержит один дефект (сферическая пора, удлиненная пора, непровар, включение); для конкретизации класса распознаваемого изображения анализируются признаки, характеризу­ ющие форму дефекта, его расположение относительно центра изобра­ жения;

если w = 2, сварное соединение содержит два произвольно распо­ ложенных дефекта; для них определяются суммарная длина, площадь, расстояние между дефектами;

если w > 3, дефект принадлежит одному из следующих классов: цепочка дефектов, скопление дефектов, дефекты, расположенные вдоль прямой, произвольно расположенные дефекты; для принятия решения о принадлежности изображения одному из вышеперечисленных классов осуществляется проверка взаимного расположения дефектов и рассто­ яния между ними.

9.3. Распознавание дефектов сварных соединений системой

Результаты распознавания в «детальном режиме» тестовой сово­ купности рентгеновских изображений реальных соединений приведены в табл. 9.2.

202

Гл. 9. Распознавание дефектов сварных соединений

 

 

 

Общее

 

 

Таблица 9.2

 

 

Количество

Количество

Ошибки

Ошибки

Классы

число

верных

ошибочных

первого

второго

изобра­

 

 

ответов

ответов

рода

рода

 

 

жений

1

 

2

3

4

5

6

Класс 1

 

195

195

0

0

0

Нет дефектов

 

 

 

 

 

 

 

Класс 2

 

183

179

4

0,02

0

Сферическая пора

 

 

 

 

 

Класс 3

 

154

150

2(непров)

0,03

0,012

Включение

 

3(у/в)

Класс 4

 

79

77

2(вкл)

0,03

0,007

Включение или

 

 

 

 

 

удлиненная пора

 

 

 

 

 

Класс 5

 

73

71

2(вкл)

0,03

0,005

Непровар

 

 

 

 

 

 

 

Класс 6

 

87

87

0

0

0

Дефекты, располо­

 

 

 

 

 

женные на одной

 

 

 

 

 

прямой

 

 

 

 

 

 

Класс 7

 

65

65

0

0

0

Дефект в виде

 

 

 

 

 

цепочки пор или

 

 

 

 

 

включений

 

 

 

 

 

 

Класс 8

 

133

133

0

0

0

Несколько произ­

 

 

 

 

 

вольно располо­

 

 

 

 

 

женных дефектов

 

 

 

 

 

Класс 9

 

62

62

0

0

0

Дефект в виде

 

 

 

 

 

скопления пор

 

 

 

 

 

или включений

 

 

 

 

 

Итого:

 

1031

1019

13

Оценка информативности пространства признаков. Для оцен­ ки информативности пространства признаков используем подход, пред­ ложенный в [14], в основе которого лежит гипотеза компактности. Согласно гипотезе признаки одного класса сгущаются в некоторой области признакового пространства, образуя тем самым некий кластер (морфология кластеров исследована в работе [29]). Следовательно, для хорошего распознавания образов желательно, чтобы расстояния между точками — признаками одного класса были малыми, а расстоя­ ния между точками разных классов были бы большими. Показателем степени выполнимости этого свойства признакового пространства или

9.3. Распознавание дефектов сварных соединений системой

203

критерием информативности признаков является величина

Wi + Wj '

Ниже, в §11.2 поясняются вычисления составляющих этой фор­ мулы.

Результаты вычислений, характеризующие разнесённость распозна­ ваемых образов (классов) в пространстве признаков по критерию приведены в табл. 9.3.

 

 

 

Таблица 9.3

Классы

Пора

Включение

Включение или

Непровар

 

сферическая

 

пора удлинённая

 

Пора сферическая

1,84

83,33

21,67

229,92

Включение

83,33

12,06

84,78

131,42

Включение или

21,67

84,78

12,92

226,97

пора удлинённая

 

 

 

 

Непровар

229,92

131,42

226,97

25,26

Результаты вычислений, характеризующие информативность прост­ ранства признаков (при числе признаков, равном 26) по критерию J приведены в табл. 9.4.

 

 

 

Таблица 9.4

Классы

Пора

Включение

Включение или

Непровар

 

сферическая

 

пора удлинённая

 

Пора сферическая

5,99

1,47

8,49

Включение

5,99

3,39

3,52

Включение или

1,47

3,39

5,95

пора удлинённая

 

 

 

 

Непровар

8,49

3,52

5,95

Из табл. 9.4 видно, что рассчитанный критерий J превышает еди­ ницу, на основании этого применённый для распознавания дефектов сварных соединений набор признаков можно считать достаточно ин­ формативным.

Режим «годен-брак». Этот режим предназначен для быстрого разделения изображений на годные и содержащие дефекты. Он преду­ сматривает обнаружение дефектов и подсчёт их числа.

Алфавит классов (образов) в этом режиме состоит только из двух элементов: первый класс — изображения сварных соединений, не со­ держащие дефектов, второй класс — изображения сварных соединений, содержащие дефекты. Изображения, относящиеся ко второму классу, могут включать как один дефект, так и несколько дефектов одного или разных типов.

204 Гл. 9. Распознавание дефектов сварных соединений

Принятие решения о принадлежности изображения к одному из классов происходит в результате вычисления признака (8.4). В слу­ чае отсутствия дефектов этот признак равен нулю, при обнаружении дефектов он равен их количеству.

Проведённые эксперименты показали отсутствие ошибок при распо­ знавании системой годных соединений, что свидетельствует о высокой эффективности признака.

Программная реализация рассмотренных алгоритмов выполнена на языке C++. Разработанный программный комплекс для автоматиче­ ского распознавания дефектов сварных соединений по результатам радиографического контроля зарегистрирован в Роспатенте [56].

9.4. Триплетные признаки распознавания дефектов сварных соединений 1

П1: Т — подсчет числа пересечений изображения со сканирующей

прямой;

 

 

 

 

 

Р — подсчет

числа

групп

ненулевых

элементов: / = а^,

если

(dij Ф 0) & ( / = 0),

то Sj

= Sj +

1 , где ciij

трейс-матрицы i =

1, n,

j = \ ,k, n — число сканирующих линий, k — число поворотов решётки,

Sj и /

первоначально считаем равными 0;

0

— максимальное значение в полученном наборе: (*S'i, ... , S^}.

Соотношение геометрических размеров дефектов.

П2: Т — определение длины отрезка, высекаемого изображением на

сканирующей линии;

 

 

 

 

Р

нахождение

максимального

элемента в

столбце

Pj =

= max { c i i j , . . . , a nj } , j =

1 ,k, n — число строк трейс-матрицы,

k

число столбцов трейс-матрицы;

 

 

 

0

нахождение максимального m ax{Pi,. . . , Рф\ и минимального

min {P j,..., Pfc} значений полученного

набора чисел,

вычисление их

отношений max /m in.

 

 

 

 

Особенности контура дефекта.

ПЗ: Т — вычисление длины отрезка, высекаемого изображением на

сканирующей прямой;

 

 

 

Р

определить

максимальный элемент:

max;

если условие

(di-ij

> ciij + 1

+ 1 < dj+ij) выполняется

и

(max—ajj) > 3,

то Pj

=

П

 

 

 

J^ m a x —a,ij. Накопленную сумму делим на число ненулевых

 

i=1

 

 

 

элементов в столбце: Sj = S j + 1, если ац ф 0. Здесь Sj первоначально полагают равным 0;

0 — подсчет числа ненулевых элементов в данном наборе чисел.

1 Экстракция и тестирование признаков выполнены аспиранткой Т. В. Никифоровой.

9.4. Триплетные признаки распознавания дефектов

205

П4: Т - вычисление длины отрезка, высекаемого изображением на

сканирующей прямой;

 

 

 

 

 

 

 

Р

определить

максимальный

элемент:

шах;

если

условие

 

> a,ij + 1)&(a*j + 1 < сц+ij)

выполняется

и

(max—a^) > 3,

то P j

П

Накопленную сумму делим на число ненулевых

= J^ m a x —a,ij.

 

i= 1

 

 

 

Ф 0. Здесь

элементы

элементов в столбце: S j = S j + 1, если

трейс-матрицы i = 2,n — 1 ,

j = \ , k ,

п

— число строк, к

число

столбцов трейс-матрицы, S j

первоначально полагают равным 0;

0 — нахождение максимального значения среди полученных чисел { P i , ..., P k} , к — число столбцов в трейс-матрице.

П5: Т — вычисление длины отрезка, высекаемого изображением на сканирующей прямой;

Р — определить номер N среднего ненулевого элемента в столбце; если условия iy < а^, i = 2, N , j = \ ,к и a*_ij > aij, i = N + 2, n,

j = 1, к выполняются, то Pj = P j + 1, где первоначально P j = 0. Чтобы полученное значение не зависело от размера дефекта, результат вычис­ ления функционала следует разделить на число ненулевых элементов в столбце;

0 — нахождение среднего значения полученных чисел:

к

П6: Т — вычисление длины отрезка, высекаемого изображением на сканирующей прямой;

Р — определить номер Переднего ненулевого элемента в столбце;

если условия

iy < а^,

i = 2, N, j = 1, к и Oj_ij > ац, i = N + 2,п,

j = \,k выполняются, то

P j = P j + 1, где

элементы трейс-

матрицы, п — число строк, к — число столбцов трейс-матрицы, перво­ начально Pj = 0. Чтобы полученное значение не зависело от размера дефекта, результат вычисления функционала следует разделить на число ненулевых элементов в столбце;

0 — нахождение максимального значения среди полученных чисел:

{Ри . . . ,Рк}.

П7: Т — вычисление длины отрезка, высекаемого изображением на сканирующей прямой;

Р — нахождение максимального элемента: { a i j , ..., anj } , где j = = 1 ,к, п — число строк матрицы (число сканирующих линий), к — число столбцов матрицы (число поворотов решетки);

0 — нахождение дисперсии полученных чисел:

с

к(к — 1 )

206 Гл. 9. Распознавание дефектов сварных соединений

П8: Т — вычисление длины отрезка, высекаемого изображением на сканирующей прямой;

Р — определение номера Переднего ненулевого элемента в столбце;

если условия а*_\j <a„ij, i

= 2, N, j = 1, k и сцij > ад,-, i = N + 2, n,

j = \,k выполняются, то P j

= P j + 1, где первоначалвно P j = 0. Чтобы

полученное значение не зависело от размера дефекта, резулвтат вычисления функционала следует разделите на число ненулеввж элементов в столбце;

0 — подсчет числа ненулевых элементов в данном наборе чисел. П9: Т — вв1числение длины отрезка, высекаемого изображением на

сканирующей прямой; Р — подсчитываем число элементов в столбце трейс-матрицы, для

которых выполняется условие

iy > (4 j

+

1 , найденное число делим

на число ненулеввж элементов:

S j = S j

+

l, если а ^ ф 0. Здесь г =

= 2 , п 1 , j = 1, k. SQ первоначально полагают равным 0;

0 — подсчет числа ненулевых элементов в данном наборе чисел.

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

П10: Т — вычисление длины отрезка, высекаемого изображением на сканирующей прямой;

Р — сумма произведений полученных значений ад,- и расстояния между прямыми г представляет собой площадь многоугольника S , описанного около изображения:

П

г = 1

0 — среднее значение имеющихся чисел.

Длина минимального прямоугольника, описанного около де­ фекта.

П11: Т — подсчет числа пересечений изображения со сканирующей линией;

Р — подсчет числа ненулевых элементов P j = P j + 1, если ад,- ф О, где первоначально P j = 0;

0 — произведение максимального значения из полученного набора

чисел на расстояние между сканирующими прямыми.

Ширина минимального прямоугольника, описанного около де­ фекта.

П12: Т — подсчет числа пересечений изображения со сканирующей линией;

Р — подсчет числа ненулевых элементов P j = P j + 1, если ад,- ф О, где первоначально P j = 0;

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

Координаты центра дефекта.

П13: Т — подсчет числа пересечений сканирующей прямой с изоб­ ражением;

9.4. Триплетные признаки распознавания дефектов

207

 

 

No N1

номер

первого

Р — вычисляем N\ -------------- В, где N\

отличного от

нуля

элемента, а

N2 — номер последнего ненулевого

 

,

 

Л72 —Л7!

 

 

элемента в столбце, слагаемое — ----- округляем до целого, В — число

сканирующих прямых, деленное на два;

 

 

 

0 — функционал:

 

 

 

 

 

1. в = arccos

J = = = ^ , а

элемент из

полученного

набора

чисел, номер

которого N\

определяется порядком

следования; b —

элемент, номер которого N 2

находят по формуле N2

= P /4 + JVi, Р —

число прямых.

 

 

 

 

 

 

 

2. Если ab > 0

и а < 0,

то

ф = <р + в + 180;

ip — угол,

который

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

Если ab >

а > 0, тоф = р>+ в.

 

 

 

Если а& <

0и а < О, тоф = <р $+180.

 

 

 

Если аЪ <

а > 0, тоф =

— в.

 

 

 

3. х = cos

ф

■у/а2 + б2, у = sin ф ■у/а2 + Ь2. Окончательный резуль­

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

Расположение дефекта относительно центра изображения.

П14: Т — определение длины отрезка, высекаемого изображением на сканирующей линии;

Р — нахождение максимального значения; 0 — нахождение порядкового номера максимального элемента сре­

ди данного набора чисел.

П15: Т — подсчет числа пересечений изображения со сканирующей линией;

Р — подсчет числа ненулевых элементов: Pj = Pj + 1, если ац ф О, где первоначально Pj = 0;

0 — определить номер N минимального значения среди получен­ ного набора чисел; произведение угла поворота сканирующей решётки Ав на N — 1 представляет собой градусную меру угла, образованного осью Ох с направлением наибольшей протяжённости дефекта.

П16: Т — определение расстояния от края изображения до дефекта; Р — нахождение дисперсии полученных чисел:5

р _

5 Е al

~ ( Е aij)

i = 1

щ = 1

J

j ~

S ( S - 1 )

 

где S(S > 1) — число ненулевых элементов в столбце;

0 — определить номер N минимального значения среди полученно­ го набора чисел; произведение угла поворота сканирующей решетки Ав на N — 1 представляет собой градусную меру угла, образованного осью Ох с направлением наибольшей протяженности дефекта.

2 0 8

Гл. 9. Распознавание дефектов сварных соединений

П17: Т — определение расстояния от края изображения до дефекта; Р — нахождение дисперсии полученных чисел:

S ( S - 1)

где Ojj — элементы трейс-матрйцы, S (S > 1) — число ненулевых элементов в столбце;

0 — определить минимальное значение среди полученного набора чисел: min {P i,..., Рф\, к — число столбцов.

Взаимное расположение дефектов.

П18: Т — сумма длин отрезков, высекаемых изображением на сканирующей прямой;

Р — вычисление отношения: m axi/m ax2, где maxi — наибольшее число нулевых элементов, расположенных между ненулевыми элемен­ тами этого столбца трейс-матрицы, т а х 2 — наибольшее число нену­ левых элементов, расположенных между нулевыми элементами этого столбца;

0 — значение, соответствующее направлению сканирования, при котором прямые пересекают наибольшее число дефектов.

П19: Т — сумма длин отрезков, высекаемых изображением на

сканирующей прямой;

элементов: / = а^-, если

Р — подсчет

групп ненулевых

(о^ Ф 0) & ( / = 0),

то Sj = Sj + 1 , где Sj

и / первоначально считаем

равными 0;

 

 

0 — вычисление |Лф —Л^|, где N i — номер максимального элемен­ та, Лф — номер минимального элемента из полученного набора чисел.

П20: Т — сумма длин отрезков, высекаемых изображением на

сканирующей прямой;

элементов: / =

если

Р — подсчет

групп ненулевых

(a,ij ф 0) & ( / = 0),

то Sj = Sj + 1 , где Sj

и / первоначально

считаем

равными 0;

 

 

 

0 — нахождение минимального значения в полученном наборе чисел.

П21: Т — подсчет числа отрезков, высекаемых изображением на сканирующей прямой;

Р — нахождение максимального элемента { а ц , . . . , a n j } , где j = = фк;

0 — вычисление |N\ —Лф|, где N i — номер максимального элемен­ та, Лф — номер минимального элемента из полученного набора чисел.

Расстояние между двумя дефектами.

П22: Т — сумма длин отрезков, высекаемых изображением на сканирующей прямой;

Р — подсчёт числа нулевых элементов, расположенных между ненулевыми элементами столбца трейс-матрицы;

0 — определить максимальный элемент max в полученном наборе чисел, найти произведение m ax—1 на расстояние между прямыми.

9.4. Триплетные признаки распознавания дефектов

2 0 9

Наибольшее расстояние между соседними дефектами.

П23: Т — сумма длин отрезков, высекаемых изображением на сканирующей прямой;

Р — определение максималвного числа нулевв1х элементов, распо­ ложенных между ненулевыми элементами столбца трейс-матрицы;

0 — определите максимальный элемент max в полученном наборе чисел, найти произведение max—1 на расстояние между прямыми.

Суммарная длина.

П24: Т — сумма длин отрезков, высекаемых изображением на сканирующей прямой;

Р — нахождение максимального элемента: {ац, ■■■, anj};

0 — значение, соответствующее направлению сканирования, при котором прямые пересекают наименьшее число дефектов.

Размер наибольшего дефекта.

П25: Т — сумма длин отрезков, высекаемых изображением на

сканирующей прямой; Р — определение максимального числа ненулевых элементов, рас­

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

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

Расстояние между наиболее удалёнными дефектами.

П26: Т — сумма длин отрезков, высекаемых изображением на сканирующей прямой;

Р — определение разности между номерами первого и последнего ненулевых элементов столбца трейс-матрицы;

0 — определить максимальный элемент max в полученном наборе чисел, найти произведение max —1 на расстояние между прямыми.

14 Федотов Н. Г.

ГЛА ВА 10

РАСПОЗНАВАНИЕ

СЛОЖНОСТРУКТУРИРОВАННЫХ

ИЗОБРАЖЕНИЙ

10.1.Сложноструктурированные изображения

вмедицинской диагностике

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

Несмотря на значительное развитие в последние годы новых ме­ тодов медицинской диагностики, цитологический и гистологический анализ по-прежнему остаются «золотым стандартом» в области диагно­ стики онкологических заболеваний. Только этот вид исследования поз­ воляет установить точный диагноз с высокой степенью достоверности. Суть цитологического и гистологического анализа заключается в по­ лучении и подготовке препарата и рассмотрении его под микроскопом при различных увеличениях на предмет выявления морфологических признаков, характерных для онкологических заболеваний.

Цитологическое исследование относится к инвазивным методам диагностики, оно достаточно дорого и болезненно для пациента. Для ранней, неинвазивной диагностики рака используется УЗИ, но оно с трудом позволяет выявить очаги поражения малых размеров, а также определить характер опухоли. Для повышения эффективности ультра­ звуковой диагностики рака щитовидной железы, в рамках программы INTAS 1, научной группе под руководством автора была поставлена

1 П р о е к т

№ 0 4 - 0 7 - 7 0 3 6 IN T A S

2 0 0 5 - 2 0 0 7

гг.: « D isc lo sin g in trin ic

b etw een

u ltra so n ic

for

h isto lo g ic al im a g e s for

im p ro v in g

th y ro id c a n c e r d ia g n o s is

a fte r th e

C h e rn o b y l

R e a c to r A ccid en t» .