Расчёт преломляющей поверхности для формирования изображения в виде линии Текст научной статьи по специальности «Математика»

РАСЧЕТ ПРЕЛОМЛЯЮЩЕЙ ПОВЕРХНОСТИ ДЛЯ ФОРМИРОВАНИЯ ИЗОБРАЖЕНИЯ В ВИДЕ ЛИНИИ

Дмитриев А.Ю.1, Досколович Л.Л.1'2 1 Институт систем обработки изображений РАН, 2 Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет)

Аннотация

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

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

Введение

Формирование изображений в виде линии является актуальным для задач лазерной оптики, радиофизики и светотехники. Задача расчёта оптического элемента для формирования изображения в виде линии состоит в определении формы преломляющей поверхности оптического элемента из условия, что преломлённые лучи приходят в точки на линии. При этом линия является каустикой семейства преломлённых лучей.

Методы решения задач данного класса в приближении геометрической оптики хорошо разработаны для дифракционных оптических элементов [14]. В этом случае задача ставится как задача расчёта эйконала светового поля на некоторой поверхности (обычно на плоскости), обеспечивающего заданные направления распространения лучей. Восстановление формы поверхности дифракционного элемента по функции эйконала предполагает использование приближения тонкого оптического элемента [5]. Использование данного приближения недопустимо при расчёте поверхностей, преломляющих лучи на большие углы (более 5...10°).

В работах [6-10] предложены итерационные методы расчёта преломляющих оптических элементов для формирования изображения в виде набора точек. Данные методы могут быть использованы для расчёта преломляющих поверхностей, формирующих сложные двумерные изображения, и в том числе изображения в виде линии [7]. Однако методы [6-10] являются вычислительно сложными из-за необходимости многократного решения прямой задачи на каждом шаге итерационного процесса.

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

1. Расчёт преломляющей поверхности для фокусировки в кривую Рассмотрим задачу расчёта поверхности Q для формирования изображения в виде линии при сферическом падающем пучке. Поверхность Q является границей раздела двух сред с показателями преломления п\ и п2. Сферический пучок формируется точечным источником, расположенным в начале координат (рис. 1). Линия задана векторной функцией одного аргумента

X(G) = (X(а), Y (а), г (а)) , а е [-П ], (1)

где о - натуральный параметр.

г ту

"2 /

/ Щ ^ст \г(а1//\У У/ X

О

А

Аг Х(ст) А2

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

Поверхность Q предлагается определить в виде огибающей семейства поверхностей по параметру о, каждая из которых фокусирует падающий пучок от источника в точку X(а) линии. Уравнение поверхности несложно получить из принципа Ферма (принцип равенства оптических путей от источника к точке на кривой) в виде [5]:

-|^ ±|X(a)-^ = У(а): п

(2)

где M = (х, у, z) , п = п2 / п1 - относительный показатель преломления, у(о) - константа при фиксированном о. Уравнение (2) определяет поверхность в

неявном виде в декартовых координатах. Знаки «±» в левой части соответствуют случаям действительного и мнимого фокуса. Поверхность (2) является картезианским овалом [5]. Уравнение огибающей определяется уравнением (2) и его производной по параметру с [10]:

( Х(ст) - М dX(ст)А |Х(ст) - М ' dст

¿у(ст)

dст

(3)

Таким образом, поверхность для формирования изображения в виде линии (1) определяется системой уравнений (2), (3).

Так как с является натуральным параметром, то из уравнения наклонов [1] можно показать, что

ау(ст) < 1. Таким образом, уравнение (3) является

¿ст

уравнением кругового конуса с единичным векто ром оси

¿Х(ст)

а(а) =-

¿ст

(4)

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

'¿у(ст)

а(ст) = arccos

¿ст

(5)

Введём единичный вектор q(ст, ф), являющийся образующей конуса лучей (3), где ф - некоторый параметр, описывающий положение вектора на конусе. Запишем уравнение преломляющей поверхности в параметрическом виде:

М(ст, ф) = (х(ст, ф), у(ст, ф), 2(<5, ф)) = = Х(ст) + q(ст, ф) • /(ст, ф), где / (ст, ф) - расстояние по направлению q(ст, ф) от точки линии Х(ст) до поверхности. Подставив (6) в (2), получим для / (ст, ф) следующее уравнение:

(6)

— |Х(ст) + /(ст,ф^(ст, ф)| ±/(ст,ф) = у(ст).

(7)

Легко видеть, что расчёт /(ст, ф) из уравнения (7) сводится к решению квадратного уравнения.

Получим ещё одно представление огибающей семейства картезианских овалов. Представим картезианский овал (2) в виде:

М(и;ст) = е(и) • /(и;ст), (8)

где и = (и, у) - некоторые криволинейные координаты, е(и) - единичный вектор направления (вектор луча от точечного источника), /(и; ст) - расстояние по направлению луча от источника до поверхности. Подставив (8) в (2), получим для /(и; ст) следующее уравнение:

1

-/(и;ст)±|Х(ст)-/(и;ст)е(и)| = у(ст).

(9)

Расчёт /(и; ст) из уравнения (9) сводится к решению квадратного уравнения. Уравнение огибающей семейства картезианских овалов, записанных в параметрическом виде (8), определяется уравнениями (8), (9) и условием [11]: д/ (и;ст) =

дст

= 0.

(10)

Последнее условие также определяет уравнение конуса (3), где точка поверхности М определена в виде (8).

Рассмотрим случай, когда размеры поверхности являются малыми по сравнению с расстоянием до линии, то есть при / / |Х| << 1. Разлагая корень в (9), получим:

У1(ст)

/ (и;ст) = -

(11)

1+п ((е (и), р (ст))) '

где у Дст) = п (у (ст) + |Х(ст)|) , р (ст) = Х(ст)/ |Х(ст)| .

Без потери общности рассмотрим случай знака «-» в знаменателе в (11). При условии (11) уравнение (10) принимает вид:

(

е(и, у),

V (ст) IV (ст)|

Л

= С(ст),

где

V(ст) = I р(ст)

¿у^ст) ёр(ст)

¿ст

¿ст

У1 (ст)

С(ст) = ^^

п аст /

к(ст)

(12)

(13)

(14)

Уравнение (12) соответствует уравнению кругового конуса с вершиной в начале координат (точка расположения источника излучения) и осью v(o). Точно такое же уравнение конуса возникает в задаче формирования диаграммы направленности, представленной функцией р(о)=Х(о)/|Х(о)| [12, 13]. Таким образом, в рассмотренном приближении (|Х| >> /) лучи, приходящие от источника на преломляющую поверхность, которые после преломления попадают в точку Х(с), лежат на поверхности кругового конуса (12).

Представление поверхности в виде (8), (9), (12) является приближённым, но в ряде случаев позволяет существенно упростить расчёт. Как показано далее, при формировании изображения в виде отрезка указанное представление позволяет свести расчёт преломляющей поверхности к интегрированию двух дифференциальных уравнений, разрешённых относительно производной.

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

Рассмотрим задачу формирования диаграммы направленности (ДН) в виде прямоугольника. Приближённым решением данной задачи является преломляющая поверхность, формирующая изображение в виде отрезка.

п

Будем рассматривать формирование мнимого изображения отрезка, так как в этом случае потери энергии, связанные с полным внутренним отражением лучей от поверхности, будут значительно меньше, чем при формировании изображения действительного отрезка. Геометрия задачи представлена на рис. 1. Поверхность для формирования мнимого изображения в виде отрезка А\А2 является огибающей семейства картезианских овалов (2), каждый из которых формирует сферический пучок, исходящий из точки Х(с), лежащей на отрезке А1А2. Угол у(с) определяет размер ДН по оси Оу и зависит от расстояния h от начала координат до отрезка А^2. Уравнение огибающей поверхности определяется выражениями (3), (6), (7).

Аналогично подходу, представленному в работах [12, 13], предлагается расчёт трёхмерной поверхности огибающей проводить через расчёт цилиндрического профиля г (ст), являющегося сечением поверхности плоскостью у = 0. Функция г (с) определяется из дифференциального уравнения вида [12, 13]:

Аг (ст) ¿ст

= -г (ст)

sm(ст-P(ст))

(п1 / п2) - cos (ст - Р(ст))

(15)

Функция Р(с) определяет распределение энергии вдоль отрезка. Для расчёта функции Р(ст) будем использовать дифференциальное уравнение первого порядка, полученное в [13] для случая формирования ДН в виде отрезка. Данное дифференциальное уравнение получено при условии, что световой поток, заключенный между конусами (12), переходит в соответствующий элемент ДН. Указанное условие соответствует приближению |Х| >> /, рассмотренному в пункте 1. Авторы не приводят в данной статье вид дифференциального уравнения для функции Р(ст) из работы [13] в силу его громоздкости. Однако следует отметить, что дифференциальное уравнение в [13] является разрешённым относительно производной. Интегрирование этого уравнения является тривиальной задачей и может быть выполнено с использованием стандартного метода Рунге-Кутты. После расчёта функции Р(ст) дальнейший расчёт огибающей поверхности может быть провёден по формуле (6), где

q(ст, ф) =

- -Р(ст)

sin ф sin I — -Р(ст)

(16)

sinI — -Р(ст) |cosф

( , ч • , ч к-г(ст)^(ст)А г (ст^ш(ст) -

Х(ст) =

с1ё (Р(ст))

-к

(17)

У(ст) = пг(ст) +|М0(ст) - Х(ст)| , М0(ст) = (г(ст)sin(ст),0,г(ст)соз(ст)) .

(18) (19)

При этом функция / (с, ф) может быть получена из квадратного уравнения (7) в виде:

/(ст, ф) = (-( п2у-( Х, q)) +

1 '' 2 -(X,q))2 -(1-1

1 - п2

п2у

-п2у2)).

(20)

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

Использование данного приближения при расчёте функции Р(ст) обуславливает приближённый характер решения задачи формирования ДН в виде прямоугольника.

На рис. 2 представлена поверхность для формирования ДН в виде прямоугольника с угловыми размерами 120° по оси Ох и 40° по оси Оу, рассчитанная по полученным формулам. Данная ДН соответствует формированию мнимого изображения отрезка в плоскости г = 10 мм. Максимальные размеры поверхности вдоль осей координат Ох, Оу, Ог составляют 16,57 мм, 10,52 мм и 10 мм, соответственно.

Рис. 2. Поверхность для формирования ДН в виде прямоугольника с угловыми размерами 120°*40°

На рис. 3 представлен результат моделирования работы рассчитанной поверхности в коммерческой программе по светотехнике ТгасеРго [14]. Результат моделирования показывает хорошее качество формируемой ДН. Среднеквадратичное отклонение полученного распределения интенсивности от равномерного составило около 15%. Световая эффективность оптического элемента (доля излучённого светового потока, попадающая в области заданной ДН) составила около 67%. Такие потери объясняются полным внутренним отражением лучей. Действительно, максимально возможный угол поворота луча при преломлении на поверхности в случае п = 1,5, п2 = 1 составляет менее 49°. При формировании ДН с рис. 3 максимальный угол поворота луча составляет

2

2

0

90°, что и обуславливает значительные энергетические потери. ffrm 80U

te

40-

Ш:

SOL 40L 0 40R 80R -80 -40 0 40 80 Degrees Degrees

Рис. 3. ДН в виде прямоугольника, формируемая поверхностью с рис. 2

Заключение Получено общее представление преломляющей поверхности для формирования изображения в виде линии, соответствующее огибающей семейства картезианских овалов. Предложено приближённое решение задачи формирования прямоугольной ДН соответствующей преломляющей поверхности, формирующей мнимое изображение отрезка. Проведён расчёт преломляющей поверхности для формирования прямоугольной ДН с угловыми размерами 120°*40°. Результаты моделирования работы поверхности показали высокое качество формируемой ДН.

Благодарности Работа выполнена при финансовой поддержке Минобрнауки РФ (государственный контракт № 07.514.11.4060).

Литература

1. Дифракционная компьютерная оптика / под ред. В.А. Сойфера - М.: Физматлит, 2007. - 736 с.

2. Гончарский, А.В. Введение в компьютерную оптику /

A.В. Гончарский, В.В. Попов, В.В. Степанов - М.: Изд-во МГУ, 1991. - 309 с.

3. Данилов, В.А. Теория когерентных фокусаторов /

B.А. Данилов, Б.Е. Кинбер, А.Е. Шилов // Компьютерная оптика. - 1987. - № 1. - С. 40-52.

4. Дмитриев, АЮ. Геометрооптический расчёт дифракционных оптических элементов для фокусировки в плоскую кривую в непараксиальном случае / А.Ю. Дмитриев, Л.Л. Досколович, С.И. Харитонов // Компьютерная оптика - 2009. - Т. 33, № 4. - C. 420-426.

5. Борн, М. Основы оптики / М. Борн, Э. Вольф. - М.: Наука, 1973.

6. Oliker, V.I. Mathematical Aspects of Design of Beam Shaping Surfaces in Geometrical Optics / V.I. Oliker // Trends in Nonlinear Analysis. - 2002. - P. 191-222.

7. Michaelis, D. Cartesian oval representation of freeform optics in illumination systems/ D. Michaelis, P. Schreiber, and A. Brauer // Optics Letters. - 2011. - V. 36, N 6. -P. 918-920.

8. Wang, L. Design of compact freeform lens for application specific light-emitting diode packaging / Kai Wang, Fei Chen, Zongyuan Liu // Optics Express. - 2010. -V. 18, N 2. - P. 413-425.

9. Fournier, F.R. Fast freeform reflector generation using source-target maps / F.R. Fournier, W.J. Cassarly, J.P. Rolland // Optics Express. - 2010. - V. 18, N 5. - P. 52955304.

10. Moiseev, M.A. Design of high-efficient freeform LED lens for illumination of elongated rectangular regions/ M.A. Moiseev L.L. Doskolovich, N.L. Kazanskiy, // Optics Express. - 2011. - Vol. 19, N S3. - P. A225-A233.

11. Математическая энциклопедия. Т.3. - М.: «Советская энциклопедия», 1982.

12. Дмитриев, АЮ. Расчёт преломляющих поверхностей для формирования диаграммы направленности в виде отрезка линии / А.Ю. Дмитриев, Л.Л. Досколович // Компьютерная оптика - 2010. - Т. 34, № 3. - C. 297-301.

13. Дмитриев, А.Ю. Расчёт преломляющих поверхностей для формирования диаграммы направленности в виде отрезка / А.Ю. Дмитриев, Л.Л. Досколович // Компьютерная оптика. - 2010. - Т. 34, № 4. - C. 476-480.

14. http://www. lambdares.com/software_products/tracepro

References

1. Diffractive Computer Optics. Chapter 3 / ed. by V.A. Soi-fer. - Moscow: Fizmatlit Publisher, 2007. - (In Russian).

2. Goncharsky, A.V. Introduction to computer optics /

A.V. Goncharsky, V.V. Popov, V.V. Stepanov - Moscow: Moscow State University Publisher, 1991. - 309 p. - (In Russian).

3. Danilov, V.A. Theory of coherent focusers / V.A. Danilov,

B.E. Kinber, A.E. Shilov // Computer optics. - 1987. -Vol. 1. - P. 40-52. - (In Russian).

4. Dmitriev, A.Yu. Geometric-optics design of diffractive optical elements to focus into a plane line / A.Yu. Dmitriev, L.L. Doskolovich // Computer Optics. - 2009. -V. 33, N 4. - P. 297-301. - (In Russian).

5. Born, M. Principles of optics / M. Born, E. Wolf -Moscow: "Nauka" Publisher, 1973. - (In Russian).

6. Oliker, V.I. Mathematical Aspects of Design of Beam Shaping Surfaces in Geometrical Optics / V.I. Oliker // Trends in Nonlinear Analysis. - 2002. - P. 191-222.

7. Michaelis, D. Cartesian oval representation of freeform optics in illumination systems / D. Michaelis, P. Schreiber and A. Brauer // Optics Letters. - 2011. - V. 36, N 6. -P. 918-920.

8. Wang, L. Design of compact freeform lens for application specific light-emitting diode packaging / Kai Wang, Fei Chen, Zongyuan Liu // Optics Express. - 2010. -V. 18, N 2. - P. 413-425.

9. Fournier, F.R. Fast freeform reflector generation using source-target maps / F.R. Fournier, W.J. Cassarly, J.P. Rolland // Optics Express. - 2010. - V. 18, N 5. - P. 52955304.

10. Moiseev, M.A. Design of high-efficient freeform LED lens for illumination of elongated rectangular regions/ M.A. Moiseev L.L. Doskolovich, N.L. Kazanskiy, // Optics Express. - 2011. - Vol. 19, N S3. - P. A225-A233.

11. Mathematical encyclopaedia. V. 3. - Moscow: "Sovet-skaya encyclopaedia" Publusher, 1982. - (In Russian).

12. Dmitriev, A.Yu. Designing of reafracting surface to gena-rate a line-shaped directivity diagram / A.Yu. Dmitriev, L.L. Doskolovich // Computer Optics. - 2010. - V. 34, N 3. - P. 297-301. - (In Russian).

13. Dmitriev, A.Yu. Design of refracting surface to generate a line-segment directivity diagram / A.Yu. Dmitriev, L.L. Doskolovich // Computer Optics. - 2010. - V. 34, N 4. - P. 476-480. - (In Russian).

14. http://www. lambdares.com/software_products/tracepro/

DESIGN OF REFRACTING SURFACE TO GENERATE A LINE-SHAPED IMAGE

A. Yu. Dmitriev1, L.L. Doskolovich1'2 1 Image Processing Systems Institute of the RAS, 2 S.P. Korolyov Samara State Aerospace University (National Research University)

Abstract

A method of designing a refracting surface to generate a line-shaped image is presented. A relationship for the surface is derived as an envelope of family of cartesian ovals. The refracting surface that generates the imaginary picture of the line-segment is presented as an approximate solution of the problem of generating the rectangular directivity diagram.

Key words: geometric optics, refractive freeform surface, directivity diagram, rays.

Сведения об авторах

Дмитриев Антон Юрьевич, кандидат физико-математических наук, научный сотрудник Федерального государственного бюджетного учреждения Российской академии наук Института систем обработки изображений РАН. Область научных интересов: геометрическая оптика, дифракционная оптика.

E-mail: tonydm@yandex.ru .

Anton Yurievich Dmitriev, a researcher of Image Processing Systems Institute of the RAS. Academic degree: Candidate of Physico-mathematical Sciences. Research interests: geometric optics, diffraction optics.

Сведения об авторе Досколович Леонид Леонидович - см. стр.340 этого номера.

Поступила в редакцию 12 мая 2012 г.