Информационные технологии
УДК 621.517, 681.142.36
ВЛИЯНИЕ ОПТИЧЕСКОГО ФАКТОРА НА ОЦЕНКУ ШЕРОХОВАТОСТИ ПОВЕРХНОСТИ ОПТИКО-ЭЛЕКТРОННЫМ КОМПЛЕКСОМ
А.Д. Абрамов
Самарский государственный технический университет 443100, г. Самара, ул. Молодогвардейская, 244
E-mail: esib@samgtu.ru
Рассматривается метод компенсации дополнительной погрешности измерения оптико-электронными средствами параметров микрорельефа поверхности, которая возникает вследствие отклонения уровня светового потока от его номинального значения. Метод основывается на определении автокорреляционных функций по изображениям анализируемых поверхностей.
Ключевые слова: метод, оптико-электронный комплекс, измерение, преобразователь, поверхность, световой поток, изображение, погрешность, компенсация, автокорреляция.
Введение. Использование оптико-электронных измерительных преобразователей совместно со средствами компьютерной обработки сигналов для оценки параметров микрорельефа механически обработанных поверхностей в настоящее время находит все более широкое применение в машиностроительном производстве [1]. При этом наибольшее распространение получили такие средства, как атомносиловая микроскопия, оптическая микроинтерферометрия, специализированные оптические микроскопы и рефлектометры. Атомно-силовая микроскопия и современная компьютеризированная микроинтерферометрия характеризуются очень высокой разрешающей способностью, которая позволяет измерять высотные параметры микрорельефа от долей нанометров до единиц микрометров. Отмеченное определяет ее основное использование при производстве высокоточных оптических приборов и изделий микроэлектроники. При этом аппаратура, реализующая эти методы, используется только в лабораторных условиях, так как она очень чувствительна к вибрациям, отличается высокой стоимостью и имеет значительные габариты и вес.
При производстве многих механических изделий в различных отраслях машиностроения такая высокая точность при измерении параметров микрорельефа не требуется. В этом случае применяются наряду с контактными профилометрами оптические приборы, использующие методы светового сечения, теневой проекции, растровый и рефлектометрический методы [2], позволяющие измерять высотные параметры микрорельефа (шероховатости) исследуемой поверхности в диапазоне 0,01 - 0,5 мкм по параметру Ra - среднему арифметическому отклонению профиля поверхно-
Алексей Дмитриевич Абрамов (к.т.н., доц.), доцент каф. электронных систем и информационной безопасности.
сти от средней линии. При производстве таких изделий на первое место выходят такие характеристики измерительной аппаратуры, как способность оперативно оценивать параметры микрорельефа непосредственно в ходе технологического процесса и в труднодоступных местах изготавливаемых изделий (внутренние полости штампов и прессформ, различные детали подшипников, отверстия небольшого диаметра и т. п.). Кроме того, аппаратура, обладающая такими характеристиками должна быть малогабаритной, обеспечивать простоту эксплуатации, иметь возможность цифровой обработки измерительной информации; обеспечивать энергосберегающий режим использования измерительных средств и иметь низкую стоимость, соизмеримую с щуповыми приборами - профилографами-профилометрами. Оптические приборы, использующие методы светового сечения, теневой проекции и растровый метод, предполагают использование специализированных микроскопов и применяются также только в лабораторных условиях. При этом требуются значительные временные затраты при производстве измерений, что полностью исключает использование данных приборов непосредственно в производственных условиях для оперативного контроля параметров шероховатости.
Отмеченные недостатки во многом устранены в рефлектометрических приборах, которые основаны на непосредственном определении мощности отраженного от исследуемой поверхности светового потока [2, 3]. Однако и этот метод имеет существенный недостаток, а именно зависимость отраженного светового потока не только от микрорельефа исследуемой поверхности, но и от условий эксплуатации прибора.
Выполненные исследования. Для проведения исследований влияния микрорельефа поверхности и внешних условий эксплуатации измерительной аппаратуры рефлектометрического типа на отраженный световой поток авторами был создан и применен измерительный комплекс, который включает в себя оптическую систему, видеокамеру, персональный компьютер и специально разработанное программное обеспечение. Данный комплекс образует информационно-измерительную систему (ИИС), которая предполагает определение таких параметров интегрального характера, как среднее арифметическое отклонение профиля исследуемой поверхности от средней линии Яа, высота неровностей профиля по конечному числу точек Яг и т. п. [4].
С помощью рассматриваемой ИИС участок анализируемой поверхности размером L1 х L2 считывают по строкам и столбцам строчно-кадровым растром телевизионной передающей камеры с ПЗС-матрицей и записывают в память ЭВМ. Таким образом, в памяти ЭВМ формируется кадр исходного полутонового изображения анализируемой поверхности формата К1 х К2 точек дискретизации -пикселей.
Для исследований были изготовлены методом бесцентрового шлифования три образца из стали ШХ15 с различной шероховатостью поверхности. Шероховатость определялась на профилографе модели SJ-201P: образец № 1 имел Яа = 0,56мкм, образец № 2 - Яа = 0,13мкм и образец № 3 - Яа = 0,084мкм.
Оптико-электронная система комплекса была настроена таким образом, что анализируемая поверхность эталонных образцов имела размер 3*2,5 мм. Световой поток падал на исследуемую поверхность под углом 45°. Формат видеокадра, записываемого в память компьютера, составлял К1 х К2 = 320*240 пикселей [4].
Отметим, что при использовании черно-белой видеокамеры на каждый пиксель изображения поверхности в памяти компьютера отводится один байт. Следовательно,
сигнал яркости может принимать значения, лежащие в диапазоне 0-255 относительных единиц. Анализ полученных видеосигналов показал, что с уменьшением шероховатости поверхности наблюдается существенное уменьшение амплитуды переменной составляющей сигнала и и, кроме того, ее сильная
зависимость от освещенности исследуемой поверхности.
Для выбранных образцов с заданным микрорельефом были выполнены исследования влияния изменений интенсивности опорного светового потока Фот на традиционные характеристики видеосигнала, такие как его среднее значение Вср по всему кадру изображения, средняя амплитуда переменной составляющей видеосигнала иср и средний период его колебания Тср на уровне Вср. Анализ полученных результатов показал, что наиболее информативной характеристикой из них, по которой можно идентифицировать (распознавать) микрорельефы с различными параметрами шероховатости, является средняя амплитуда переменной составляющей видеосигнала и . Г рафик зависимости и = / (Фо) представлен на рис. 1.
200 400 600 800 ЮООФо-Ю'лм
Рис. 1. Влияние изменения интенсивности опорного светового потока на среднюю амплитуду переменной составляющей видеосигнала и ср :
ср
1-Ка = 0,56мкм, 2 - Ка = 0,13мкм, 3 - Ка = 0,084мкм
Примем за номинальное значение опорного светового потока величину Фон =600-10~3лм, а для величин Фо . , Фо- значения 200-10~3лм и
ОН о 111111 о 11ЮХ
1000 -10~3лм соответственно. Выбранный диапазон значений вполне соответствует рабочим значениям светового потока, встречающимся в практике применения рассматриваемых ИИС в производственных условиях. Отмеченные изменения светового потока могут происходить по разным причинам, например вследствие колебаний напряжения питания источников света, изменения прозрачности атмосферы на рабочем месте при выполнении технологического процесса изготовления изделия, неконтролируемых бликов на исследуемой поверхности от посторонних источников света и т. п.
Отметим, что величина отраженного светового потока Фот от микрорельефа исследуемой поверхности является величиной промежуточной, на основании которой определяется параметр микрорельефа Ка. На основе полученных эксперимен-
тальных данных был построен график зависимости иср = /(Яа) для разных значений Фо =сош^ приведенный на рис. 2. Его анализ показал, что если в результате измерения величина иср приняла значение 24 отн. ед., то при опорном световом потоке ФощП =200-10~3лм этому значению иср будет соответствовать значение Яа = 0,48мкм. При номинальном световом потоке Фон =600 10 ~3лм - значение Яа = 0,35мкм и при значении светового потока Ф____________ах =1000 10 3 лм будет соответ-
О ^ОаХ »
ствовать значение Яа = 0,29мкм. Таким образом, изменение светового потока от нижней границы рабочего диапазона до его верхней границы под действием влияющей величины создает абсолютную погрешность в определении Яа, равную А = 0,19мкм . При этом относительная погрешность 8 , определяемая как отношение (А /ЯаНОМИН )• 100 %, составит величину 8 =54,3%, где ЯаНОМИН = 0,35мкм.
Рис. 2. Зависимость иСР = /(Яа, мкм) для различных значений опорного светового потока: 1 - Фо = 200-10 3лм , 2 - Фон = 600-10 3лм, 3 - ФО = 1000-10 3лм
Таким образом, определение высотных параметров микрорельефа на основе непосредственного измерения средней амплитуды видеосигнала иср может привести
к значительным погрешностям. Эти непредусмотренные воздействия на освещенность исследуемой поверхности можно охарактеризовать с помощью функции влияния /вл (АФ) [5]. Для определения доминирующего характера этой функции были
вычислены отношения величин иср для разных параметров Яа при значениях светового потока Ф_________. , Ф ,,, Ф_________ . Обозначив эти отношения как
о ^ош он ^ о ^ОаХ
К1,2 = иСР(Яа=0,56мкм) / иСР(Яа=0,13мкм) , К1,3 = иСР(Яа=0,56мкм) / ^СР(Яа=0,084мкм) и
К 2,3= и СР( Яа=0,13 мкм)/ 11 СР (Яа=0,084 мкм), полуЧим их соответствуюЩие значения: при
световом потоке Фот.п К12= 1,31, К13=2,09, К23 =1,59; при световом потоке Фон К12= 1,32, К13= 2,10, К23=1,б0 и при световом потоке Фотах К12= 1,33, К13= 2,19 и К2 3 =1,60. Полученные значения указанных отношений позволяют 48
обоснованно утверждать, что функция влияния /вл (АФ) нестабильности опорного
светового потока и соответственно дополнительная погрешность, вызываемая действием данного оптического фактора, имеют мультипликативный характер.
Наличие мультипликативного характера функции /вл (АФ), воздействующей на
выходную величину оптико-электронного преобразователя (видеокамеры), диктует выбор логометрического способа устранения рассматриваемой дополнительной погрешности [6]. Однако применение логометрического способа предусматривает введение в состав рассматриваемой аппаратуры дополнительного осветительного канала и второй ПЗС-матрицы, содержащих компенсационные величины. Требуемое введение дополнительных устройств в измерительную аппаратуру существенно увеличивает ее габариты и усложняет ее использование в производственных условиях. Поэтому применительно к используемой разновидности измерительного преобразования был выдвинут иной принцип коррекции дополнительной погрешности, основанный на применении специальной математической обработки исследуемого изображения микрорельефа поверхности.
Реализация предлагаемого принципа коррекции осуществлена с привлечением теории оптимальной линейной фильтрации сигналов известной формы [7, 8]. Известно, что если на входе приемника-фильтра действует двумерная величина х(п1, п2), по отношению к которой этот фильтр согласован, то выходной сигнал у(п1, п2) будет пропорционален автокорреляционной функции двумерного входного сигнала и, кроме того, будет достигнуто максимальное отношение сигнала к помехе на выходе фильтра.
Нахождение двумерной автокорреляционной функции в рамках рассматриваемой задачи осуществлялось согласно известному выражению [8]
N -ш2 -1
ЕЁ (м(п1’ П2 ) - ти ) • (х(п1 - k1, П2 - k2 ) - тх )
Гш (*1, *2>= -----------------------------------------------------> (!)
СТ1 • °2
где через и(п1, п2) обозначен фрагмент изображения (эталон), который расположен внутри зоны поиска х(п1, п2) , а а1 и сг2- средние квадратические отклонения величин и(п1, п2) и х(п1, п2) от их математических ожиданий ти и тх. Анализируя выражение (2), можно отметить, что его числитель и знаменатель представляют собой суммы произведений сигналов текущего фрагмента изображения х(п1, п2) и эталона и(п1, п2) . Следовательно, эти сигналы будут подвержены, как было показано выше, мультипликативному воздействию функции влияния /вл (АФ). Тогда выражение
для нормированного коэффициента корреляции (2) с учетом воздействия функции влияния примет вид
^ -Ш2 -1
ЕЕ (и(п1’ п2) - ти ) • /вл (АФ) • (х(п1 - ^ п2 - *2) - тх ) • /вл (АФ)
щ =0 «2 =0
СТ1 • /вл (АФ) ^2 • /вл (АФ)
Как видно из приведенного выражения для гхи(к1,к2), функции влияния /вл (АФ), стоящие в числителе и знаменателе, сокращаются. Таким образом, ис-
Гхи (% к2/
пользуемая автокорреляционная функция обладает искомым компенсационным свойством, не требуя при этом введения дополнительной аппаратуры.
Для вычисления двумерной автокорреляционной функции в кадре формата К1 х К2 точек дискретизации по специально разработанной программе начиная с первой строки выделялась полоса шириной Ы2 пикселя. По центру этой полосы задавался эталон размером Ы1 х Ы2 пикселей. Затем эталон начиная с крайней левой позиции исходного изображения перемещался по выделенной полосе с шагом в 1 пиксель. Схема формирования и перемещения эталона по исходному изображению исследуемой поверхности приведена на рис. 3. При каждом совмещении эталона и(п1,п2) и текущего фрагмента полутонового изображения и(п1,п2) подсчитывался коэффициент корреляции по формуле (1). В связи с тем, что сравниваемые фрагменты и(п1,п2) и и(п1,п2) берутся из одного и того же изображения поверхности, гии (Ар А2) является коэффициентом автокорреляции.
сдвиг
ЭТАЛОНА
ТЕКУЩИЙ
ФРАГМЕНТ
N1
N1
г I
г
о
о
1
и
К1
Рис. 3. Схема формирования и перемещения эталона
По окончании вычисления коэффициентов автокорреляции в первой полосе задавалась следующая полоса того же формата в исходном изображении, но смещенная вниз на один пиксель. В этой полосе по центру задавался новый эталон, но с теми же размерами, что и предыдущий, и выполнялись те же самые действия для вычисления коэффициента автокорреляции и т. д. После обработки всего кадра исходного изображения по отмеченной методике в запоминающем устройстве ИИС формировалась матрица размерностью М1 х М2 коэффициентов автокорреляции, то есть двумерная автокорреляционная функция. Для этой функции определялось среднее значение переменной составляющей иср .
Заключение. Экспериментальная проверка предложенного метода коррекции дополнительной погрешности показала, что средний уровень переменной составляющей двумерной автокорреляционной функции остается постоянным для поверх-50
ности с Ra = const при изменении светового потока от Фо . до Фо . График
а о о max ^ ^
зависимости Ucp = f (АФ) представлен на рис. 4.
Ucp 20
15
Ю
5
О 200 400 600 800 ЮОО
Фо-Ю ЛЛЛ
Рис. 4. Влияние изменения опорного светового потока на величину Ucp :
1 - Яа = 0,56 мкм, 2 - Яа = 0,13 мкм, 3 - Яа = 0,084 мкм
Сам принцип определения неизвестной шероховатости по рассматриваемому методу заключается в следующем. Для эталонных образцов с известной шероховатостью, определенной по гостовским методикам, например с помощью профилографа, строится зависимость Яа = f (иср) и для заданной вероятности распознавания Яа
определяется доверительный интервал Iр = f (иср) . Для неизвестной, исследуемой шероховатости по рассмотренному алгоритму обработки изображения вычисляют иср. Используя зависимость для доверительного интервала, находят исрш1п и исршах. Затем для найденных значений иср, исрш1п и исршах по уравнению Яа = f (иср) определяют Яа , ЯаШп и Яашах для исследуемой поверхности. Практическая реализация корреляционного метода при оценке шероховатости цилиндрической поверхности бомбинированного ролика подшипника ступицы колеса автомобиля ВАЗ приведена в [4]. На данный метод определения шероховатости поверхности получен патент Российской Федерации № 2413179 «Способ контроля шероховатости поверхности».
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Уайтхауз Д. Метрология поверхностей. Принципы, промышленные методы и приборы: Научное издание / Д. Уайтхауз. - Долгопрудный: Интеллект, 2009. - 472 с.
2. Дунин-Барковский И.В., Карташова А.Н. Измерение и анализ шероховатости, волнистости и не-круглости поверхности. - М: Машиностроение, 1987. - 232 с.
3. Ким К.Ю. Исследование и разработка оптического метода бесконтактного контроля шероховатости поверхностей. Дисс. ... канд. техн. наук. - М.: Гос. техн. ун-т электронной техники, 2009. 150 с.
4. Абрамов А.Д. Определение микрогеометрии поверхности деталей машин и механизмов на основе компьютерной обработки их видеоизображений. - М.: Вестник компьютерных и информационных технологий, 2008. - № 11. - С. 18-25.
5. ГОСТ 8.009-72. Нормируемые метрологические характеристики. Средства измерения.
6. Измерение электрических и неэлектрических величин: Учеб. пособие для вузов / Н.Н. Евтихеев, Я.А. Купершмидт, В.Ф. Папуловский, В.Н. Скугоров. - М.: Энергоатомиздат, 1990. - 352 с.
7. Прэтт У. Цифровая обработка изображений. Том 2. - М.: Мир, 1982. - 790 с.
8. Даджион Д., Мерсеро Р. Цифровая обработка многомерных сигналов. - М.: Мир, 1988. - 486 с.
Статья поступила в редакцию 13 марта 2012 г.
INFLUENCES OF OPTICAL FAKTOR ON ESTIMATE OF SURFACE ROUGHESS BY OPTICO-ELECTRONIC COMPLEX
A.D. Abramov
Samara State Technical University
244, Molodogvardeyskaya st., Samara, 443100
This paper presents a method of compensation of the supplementary measurement error by the optic-electronic means of the surface microgeometry parameters, which appears as a result of the value of reference incident light deviation. The Method is based on determining autocorrelative functions by the vision of surfaces being analyzed.
Keywords: method, optic-electronic complex, measurement, convertes, surface, reference incident light, image, error, compensation, autocorrelation.
Aleksey D. Abramov (Ph.D. (Techn.)), Associate Professor.