СПЕКТРАЛЬНЫЙ АНАЛИЗ ФОРМУЛ РАСШИФРОВКИ МЕТОДОМ ПОШАГОВОГО ФАЗОВОГО СДВИГА

Спектральный анализ формул расшифровки методом пошагового фазового сдвига

В.И. Гужов, С.П. Ильиных, Г.А. Поздняков Новосибирский государственный технический университет, Новосибирск, Россия

Аннотация: Для оценки волнового фронта интерферограммы в интерферометрии с фазовым сдвигом (PSI) часто используют формулы расшифровки, построенные на основе алгоритмов фазового сдвига. Все Галгоритмы PSI используют в качестве входных данных смещенные по фазе интерферограммы и вычисляют фазу интерференционных полос. Наиболее известным алгоритмом PSI, специально разработанным для уменьшения источника систематической ошибки, был 5-точечный алгоритм Швидера-Харихарана. С тех пор были опубликованы десятки алгоритмов, предназначенные для уменьшения конкретных источников ошибок при определении фазы. В электротехнике функция Фурье преобразования частоты (ПФ) при конструировании линейных фильтров является стандартным инструментом проектирования. В последнее время ПФ также используется для разработки новых алгоритмов. В этой статье мы рассматриваем метод определения спектральных характеристик по числителю и знаменателю известных алгоритмов. Приведены АЧХ для формул расшифровки с различным количеством точек.

Ключевые слова: фазовый сдвиг, формула расшифровки, спектральный анализ, амплитудно-частотная характеристика.

ВВЕДЕНИЕ

При регистрации интерференционной картины ее интенсивность может быть описана как

I

В этом выражении три неизвестных: I (х, y) -

средняя интенсивность, V(x, y) - контраст, а

ф(x, y) = (рр (x, y) - ( (x, y) - фазовая разность

между интерферирующими фронтами. В [1] показано, что для однозначного определения разности фаз, необходимо не менее трех взаимно независимых уравнений. Взаимная независимость достигается, например, если в каждое из уравнений (1) добавить постоянный фазовый сдвиг g. В этом случае можно

определить значения разности фаз ((x, y) из

системы уравнений, если значения средней яркости и контраста не меняются в точке (x, y) .

Iiix,y) = I0

В этом уравнении измеряются значения интенсивностей I ( x, y) при известных фазовых

сдвигах g. По этим значениям с помощью

формул расшифровки определяются значения фазового сдвига. Найти вид этих формул можно, непосредственно решая

тригонометрические уравнения вида (2), или с помощью использования матричных выражений описанных в [2].

Общий вид формул расшифровки будет иметь вид

ф( x, y ) = arctg

( N-l \

IК I * ^ )

i=0

где S. и C

!

V i=0 y

соответствующие значения

синусов и косинусов при заданом угле сдвига 8 , а коэффициенты, которые определятся при

получении соотвествующих формул. Если 8 меняются

как

g = 0o, g=—, g = —, K ,g = 2л , то

1 , 2 n ' 3 N ' ' "

I Ii sing

ф( x, y) = arctan

I Iicos gi

Выражение (7) можно переписать в виде

ZI sin I 2л—

' V N

ф( x, y) = arctan N--^-

I Ii cos I 2л

N

Данное выражение совпадает с формулой вычисления фазы первой гармоники преобразования Фурье.

Угол фазового сдвига можно задавать с помощью модуляторов света, простейщим из которых является зеркало, закрепленное на пьезокерамике. Метод определения разности фаз при последовательном задании известного сдвига (PSI - Phase Shift Interferometry) впервые предложен в работе [3] и привлек внимание специалистов возможностью получения 2019, № 4 (30) 51

i=0

i=0

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

К сожалению, значения фазовых сдвигов задаются с некоторой ошибкой, измеряются значения интенсивности в каждой точке тоже с некоторой погрешностью. Это сказывается на погрешности при определении фазового сдвига.

1. ИСТОЧНИКИ ОСНОВНЫХ ПОГРЕШНОСТЕЙ ВОЗНИКАЮЩИХ ПРИ ОПРЕДЕЛЕНИИ РАЗНОСТИ ФАЗ МЕТОДОМ ПОШАГОВОГО ФАЗОВОГО СДВИГА

Как видно из выражения (3), погрешность измерения фазы зависит от погрешности задания фазового сдвига и погрешности измерения интенсивности.

Дж. Швидер [4] получил выражение для погрешности определения фазы от абсолютной ошибки при установке сдвига при условии малости ошибок задания сдвига и измерения интенсивности

Аф « Фг- и AI « Il :

^ А - C cos 2ф - S sin 2ф

Аф = arctg

l=1

m - C sin 2ф - S cos 2ф

где C = £Аф cos 2^. , S = sin 2^. .

i=1

i=1

В [5] Б.В. Карпюком получено выражение для погрешности определения фазы от суммарных ошибок при измерении интенсивности и установки сдвига при тех же условиях

2. АНАЛИЗ ФОРМУЛ РАСШИФРОВКИ НА ОСНОВЕ СПЕКТРАЛЬНОГО АНАЛИЗА

Спектральный анализ в радио- и электротехнике давно стал одним из основных методов анализа, в частности, использование Фурье преобразования является стандартным способом анализа и синтеза цифровых фильтров [6]. Спектральный анализ в теории алгоритмов пошагового фазового сдвига (PSI) впервые использован в 1990-х годах, и основан на работе Freishlad и Koliopoulos [7]. Применение спектрального анализа к синтезу алгоритма PSI можно также встретить в [8]. В 1996 году Surrel [9] разработан алгебраический подход к анализу алгоритмов, которые основан на так называемом характеристическом полиноме, связанным с фильтром квадратурных составляющих.

На Рис. 1 показан способ формирования исследуемого сигнала. Сигнал формируется из последовательности отсчетов значений интенсивности I

интерференционной картины при изменении фазового сдвига §..

Рис. 1. Способ формирования сигнала

Аф*£

i=1

дф

ds;

(±АЗ M

i=i

дф

д!

дф _ cos ôfj + sin ôlSd

^ C )2+( Sd )2 дф = sinôiCd - cosôiSd dIi " (Cd )2+ (Sd )2

где m

m

число фазовых сдвигов, а

т

С =^Мг 008^. , Б, = £Д/, . 1 =1 1=1 В реальных системах погрешность задания сдвига и измерения интенсивности может составлять 20%-30%. В этом случае условие (4) не выполняется.

Такой подход основан на представлении формулы расшифровки в виде цифрового фильтра с конечной импульсной характеристикой

ф = arctg

I

где - (*) означает операцию свертки, ^ (п) и

к0 (п) - импульсная характеристика синфазной

и квадратурной составляющих, т.е. числителя и знаменателя в формуле (9).

Представим интерференционный сигнал в следующем виде

1.(х, у, /) = 10(х, у,/)[1 + У(х, у)со8(8. + ю0 (х, у)/)]'

здесь t

- некоторое произвольное

вещественное число, такое, что г е Я, щ -

несущая пространственная частота в радианах. Тогда уравнение расшифровки (9) примет вид

ф = arctg

К С0Я

_ 10, йШ

Импульсная характеристика определяется коэффициентами числителя и знаменателя формулы (3) и с учетом соотношения

*) Z ckS(t - nT) = h(t )e

k

h(t -10) = eiao(t-t*) может быть записана в следующем виде

hN (t) = [8(t - t0)cos(v0t0) + S(t - t1)cos(ffloti) + K

K + S(t - tn-1)cOs(v0tn-1) + S(t - tn )cOs(v0tn )]

hD(t) = i[S(t - t0)sin(v0t0) + 8(t -i1)sin(«0i1) + K

K + S(t - tn-1)sin(«0tn-1) + S(t - tn )sin(v0tn )]

или

hN(t) = ZS(t -1*)ck

k

hD(t) = iZS(t -1*)sk

где и Sk - соответствующие коэффициенты

формулы (5), I - мнимая единица. Отсюда

h

Импульсная и амплитудно-частотная передаточная характеристики связаны между собой преобразованием Фурье H (а>) = 3[h(t)].

Поэтому можно записать

HN

HD() = i[s0 ei(t-to) + s ei(t-t,)+K + sn-1 ei(t-t-) + sn ei(t-t-)]

или

Hn V) = Z ck ei(v-vk) к

HD (v) = iZ sk e

i(m-mt) .

Амплитудно-частотную и фазо-частотную характеристики можно получить из выражений

H

ф(ю) = tg

Re

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

Рассмотрим получение частотных характеристик для известных формул расшифровки [10-12].

Пусть имеется формула расшифровки с шестью фазовыми сдвигами

8=| 0, -, —, ж, — —Т . Тогда

3 3 3 3

получим следующую формулу расшифровки

tsA= 1 210 + I1 - I2 - 2I3 - I4 + I5 .

g S Ii + I2 -14 - I5 Для неё находится импульсная характеристика

hN

-2S

Hn

-28

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

Я / \ л Ню-ох) , i

N

-2ei

HD

На рис. 2 показан модуль амплитудно-частотной характеристики шеститочечной формулы расшифровки (26).

\Н{ф)\

7 .6 -5 -4 -S -2 -1 о 1 2 3 4 5 6 7 8 а/0

Рис. 2.АЧХ шеститочечной формулы расшифровки

Из Рис. 2 следует, что данная формула расшифровки подавляет 2, 3, 4, 5 и 6 и кратные им гармоники.

Для пятиточечной формулы расшифровки Schwider-Hariharan с фазовыми сдвигами

ж ж Т

0, —, ж 2'

8<=\-я, -

2

I0 - 2I2 +14

k

k

Для этой формулы

hN

hD

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

HN (©) = 2

si

Яд

На рис. 3 показан модуль частотной характеристики пятиточечного алгоритма.

Рис. 3. АЧХ пятиточечной формулы расшифровки

Из Рис. 3 следует, что данная формула расшифровки подавляет 2, 3, 4 и кратные им гармоники, в то же врем хорошо выделяет первую гармонику.

Представляет интерес сравнение АЧХ различных формул расшифровки. На Рис. 4 показаны АЧХ пятиточечной [1] (кривая 1) и шеститочечной (кривая 2) формул

расшифровки.

Рис. 4 - Сравнение АЧХ пятиточечной (кривая 1) и шеститочечной (кривая 2) формул расшифровки

Из Рис. 4 следует, что АЧХ шеститочечной формулы расшифровки имеет меньшую амплитуду боковых лепестков. Поэтому этот метод лучше подавляет шумы. Кроме того, чем уже лепесток для первй гармоники, тем ближе этот алгоритм к выделению правильного значения фазовой разности.

ВЫВОДЫ

В работе рассмотрен метод спектрального анализа формул расшифровки используемых в методе пошагового фазового сдфига. Известные методы анализа погрешности основаны на малости ошибок измерения интенсивности и задания фазовых сдвигов и не могут быть применены для реальных измерительных

систем. В то же время методы спектрального

анализа показывают качество используемых

алгоритмов.

ЛИТЕРАТУРА

[1] Гужов В.И., Ильиных С.П. Оптические измерения. Компьютерная интерферометрия. -Москва: Изд-во Юрайт, 2018. (ISBN: 978-5-53406855-9) - 258с.

[2] Generic algorithm of phase reconstruction in phase-shifting interferometry /Guzhov V., Ilinykh S., Kuznetsov R., Haydukov D.// Optical Engineering, -2013.-Vol.52(3) - pp. 030501-1 - 030501-2.

[3] P. Carré, "Installation etutilisation du comparateurphotoélectrique et interférentiel du Bureau International des Poids et Mesures," Metrologia 2, 13-23 (1966).

[4] Schwider J., Burow R.,Elsner K.-E.,Grzanna J., Spolaczyc R.,Mezkel K. Digital wave-front measuring interferometry: some systematic error sources// Applied Optics.- 1983.- V.22.- No.21.-P.3421-3432.

[5] Карпюк Б.В., Солодкин Ю.Н. Анализ погрешностей измерения фазы интерферометров с управляемым фазовым сдвигом // Автометрия.-1992.-№6.-С.16-21

[6] M. Servin, J. C. Estrada, and J. A. Quiroga, "Spectral analysis of phase shifting algorithms," Opt. Express 17(19), 16423-16428 (2009).Wyant J.C. Interferometric optical metrology: basic system and principles // Laser Focus.- 1982.- pp.65-67.

[7] K. Freischlad, and C. L. Koliopoulos, "Fourier description of digital phase-measuring interferometry," J. Opt. Soc. Am. A 7(4), 542-551 (1990).

[8] D. W. Phillion, "General methods for generating phase-shifting interferometry algorithms," Appl. Opt. 36(31), 8098-8115 (1997).

[9] Y. Surrel, "Design of algorithms for phase measurements by the use of phase stepping," Appl. Opt. 35(1), 51-60 (1996).

[10] Wyant J.C., Creath K. Recent advances in interferometric optical testing // Laser Focus.- 1985. -pp.118-132.

[11] Wyant J.C. Interferometric optical metrology: basic system and principles // Laser Focus.- 1982.- pp.6567.

[12] D. Malacara, ed., Optical Shop Testing, 3rd ed. (John Wiley and Sons, New York, 2007).

Владимир Иванович Гужов

- профессор кафедры ССОД Новосибирского Государственного Технического университета, профессор, доктор

технических наук. Автор более 250 научных работ. Область научных интересов: высокоточные интерференционные измерения.

E-mail: vigguzhov@gmail.com

Григорий Александрович Поздняков - аспирант кафедры систем сбора и обработки данных НГТУ. Является автором более 10 научных работ. Область научных интересов -компьютерные методы гологрфии.

20 июля 2019 г.

Bilinear Approximation to Increase the Size of a Numerical Sequence in Projection

Methods

V.I. Guzhov, S.P. Ilinykh, G.A. Pozdnyakov

Novosibirsk State Technical University, Novosibirsk, Russia

Abstract. To evaluate the wavefront of interferograms in phase shift interferometry (PSI), decryption formulas based on phase shift algorithms (PSA) are often used. All PSAs use phase-shifted interferograms as input and calculate the phase of the interference fringes. The first and best-known PSA specifically designed to reduce the source of bias was the 5-point Shveder-Hariharan PSA (SH-PSA). Since then, dozens of PSAs have been published, designed to reduce specific sources of error in phase determination. In electrical engineering, the Fourier function of frequency conversion (FT) in the design of linear filters is their standard design tool. Recently, FT is also used to develop PSA. In this article, we consider a method for determining spectral characteristics from the numerator and denominator of the PSA. The frequency response for decoding formulas with a different number of points is given.

Key words: phase shift, decoding formula, spectral analysis, amplitude-frequency characteristic..

REFERENCES

[1] Guzhov V.I., Il'inyh S.P. Opticheskie izmereniya. Komp'yuternaya interferometriya. - Moskva: Izd-vo YUrajt, 2018. (ISBN: 978-5-534-06855-9) - 258s.

[2] Generic algorithm of phase reconstruction in phase-shifting interferometry /Guzhov V., Ilinykh S., Kuznetsov R., Haydukov D.// Optical Engineering, -2013.-Vol.52(3) - pp. 030501-1 - 030501-2.

[3] P. Carré, "Installation etutilisation du comparateurphotoélectrique et interférentiel du Bureau International des Poids et Mesures," Metrologia 2, 13-23 (1966).

[4] Schwider J., Burow R.,Elsner K.-E.,Grzanna J., Spolaczyc R.,Mezkel K. Digital wave-front measuring interferometry: some systematic error sources// Applied Optics.- 1983.- V.22.- No.21.-P.3421-3432.

[5] Guzhov V.I. Metody izmerenija 3D profilja obektov. Kontaktnye, trianguljacionnye sistemy i metody strukturirovannogo osveshhenija.-Novosibirsk: Izd-vo NGTU, 2015. — 82s..

[6] M. Servin, J. C. Estrada, and J. A. Quiroga, "Spectral analysis of phase shifting algorithms," Opt. Express 17(19), 16423-16428 (2009).Wyant J.C. Interferometric optical metrology: basic system and principles // Laser Focus.- 1982.- pp.65-67.

[7] K. Freischlad, and C. L. Koliopoulos, "Fourier description of digital phase-measuring interferometry," J. Opt. Soc. Am. A 7(4), 542-551 (1990).

Сергей Петрович

Ильиных - доцент кафедры Вычислительная техника Новосибирского Государственного Технического

университета, кандидат технических наук, доцент. Является автором более 130 научных трудов. Область научных

интересов: разработка алгоритмов анализа

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

Статья поступила

[8] D. W. Phillion, "General methods for generating phase-shifting interferometry algorithms," Appl. Opt. 36(31), 8098-8115 (1997).

[9] Y. Surrel, "Design of algorithms for phase measurements by the use of phase stepping," Appl. Opt. 35(1), 51-60 (1996).

[10] Wyant J.C., Creath K. Recent advances in interferometric optical testing // Laser Focus.- 1985. -pp.118-132.

[11] Wyant J.C. Interferometric optical metrology: basic system and principles // Laser Focus.- 1982.-pp.65-67.

[12] D. Malacara, ed., Optical Shop Testing, 3rd ed. (John Wiley and Sons, New York, 2007).

Vladimir Guzhov - Professor, Department of Special Purpose Traffic Control, Novosibirsk State Technical University, Professor, Doctor of Technical Sciences. The author of more than 250 scientific papers. Research interests: high-precision interference

measurements.

E-mail: vigguzhov@gmail.com Novosibirsk, Prosp. of K. Marksa, h.20, NSTU.

Sergey Petrovich Ilyinykh -

Associate Professor,

Department of Computer Engineering, Novosibirsk State Technical University, Candidate of Technical Sciences, Associate

Professor. He is the author of more than 130 scientific papers. Research interests: development of image analysis algorithms in optical measuring systems. Novosibirsk, Prosp. of K. Marksa, h.20, NSTU.

Grigory Aleksandrovich

Pozdnyakov - graduate student of the Department of Data Acquisition and Processing Systems, NSTU. He is the author of more than 10 scientific papers. The area of scientific interests is computer methods of holography. Novosibirsk, Prosp. of K. Marksa, h.20, NSTU. The paper has been received on 20/07/2019.