stringtranslate.com

Полиномы Цернике

Первые 21 полином Цернике, упорядоченные по вертикали в соответствии с радиальной степенью и по горизонтали в соответствии с азимутальной степенью

В математике полиномы Цернике представляют собой последовательность полиномов , ортогональных на единичном круге . Названные в честь физика-оптика Фрица Цернике , лауреата Нобелевской премии по физике 1953 года и изобретателя фазово-контрастной микроскопии , они играют важную роль в различных разделах оптики , таких как лучевая оптика и визуализация. [1] [2]

Определения

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

(четная функция по азимутальному углу ), а нечетные полиномы Цернике определяются как

(нечетная функция по азимутальному углу ), где m и n — неотрицательные целые числа с n ≥ m ≥ 0 ( m = 0 для сферических полиномов Цернике), — азимутальный угол , ρ — радиальное расстояние , а — радиальные полиномы, определенные ниже. Полиномы Цернике обладают свойством быть ограниченными диапазоном от −1 до +1, т. е . Радиальные полиномы определяются как

для четного числа nm , в то время как для нечетного числа nm он равен 0. Особое значение —

Другие представления

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

.

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

для nm четно.

Обратное отношение расширяется для фиксированного в

с рациональными коэффициентами [3]

для четных .

Фактор в радиальном многочлене может быть разложен в базисе Бернштейна для четных или умножен на функцию для нечетных в диапазоне . Радиальный многочлен может быть выражен конечным числом многочленов Бернштейна с рациональными коэффициентами:

Последовательные индексы Нолла

Приложения часто включают линейную алгебру, где интеграл по произведению полиномов Цернике и некоторого другого фактора строит матричные элементы. Чтобы перечислить строки и столбцы этих матриц по одному индексу, Ноллом было введено обычное отображение двух индексов n и l в один индекс j . [4] Таблица этой ассоциации начинается следующим образом (последовательность A176988 в OEIS ).

Правило следующее.

Стандартные индексы OSA/ANSI

Одноиндексные полиномы Цернике OSA [5] и ANSI с использованием:

Индексы Fringe/Университета Аризоны

Схема индексации Fringe используется в коммерческом программном обеспечении для оптического проектирования и оптического тестирования, например, в фотолитографии . [6] [7]

где - знаковая или signum функция . Ниже перечислены первые 20 чисел полосы.

индексы Вайанта

Джеймс С. Уайант использует схему индексации «Fringe», за исключением того, что она начинается с 0 вместо 1 (вычитание 1). [8] Этот метод широко используется, включая программное обеспечение для анализа интерферограмм в интерферометрах Zygo и программное обеспечение с открытым исходным кодом DFTFringe.

Формула Родригеса

Они удовлетворяют формуле Родригеса

и может быть связана с полиномами Якоби как

.

Характеристики

Ортогональность

Ортогональность в радиальной части имеет вид [9]

или

Ортогональность в угловой части представлена ​​элементарными

где (иногда называемый фактором Неймана , поскольку он часто появляется вместе с функциями Бесселя) определяется как 2, если и 1, если . Произведение угловой и радиальной частей устанавливает ортогональность функций Цернике относительно обоих индексов, если интегрировано по единичному кругу,

где — якобиан круговой системы координат, причем и — оба четные.

Преобразование Зернике

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

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

Коэффициенты Цернике тогда можно выразить следующим образом:

В качестве альтернативы можно использовать известные значения фазовой функции G на круговой сетке для формирования системы уравнений. Фазовая функция извлекается с помощью взвешенного произведения неизвестного коэффициента с (известными значениями) полинома Цернике по единичной сетке. Следовательно, коэффициенты также могут быть найдены путем решения линейной системы, например, путем обращения матрицы. Быстрые алгоритмы для вычисления прямого и обратного преобразования Цернике используют свойства симметрии тригонометрических функций , разделимость радиальных и азимутальных частей полиномов Цернике и их вращательные симметрии.

Симметрии

Отражения тригонометрических функций приводят к тому, что четность относительно отражения вдоль оси x равна

для l ≥ 0,
для l < 0.

Сдвиги тригонометрических функций на π приводят к тому, что четность относительно отражения точки в центре координат равна

где можно было бы также записать, потому что четные числа являются единственными случаями получения неисчезающих полиномов Цернике. (Если n четное, то l также четное. Если n нечетное, то l также нечетное.) Это свойство иногда используется для классификации полиномов Цернике на четные и нечетные полиномы с точки зрения их угловой зависимости. (также можно добавить еще одну категорию с l = 0, поскольку она имеет особое свойство отсутствия угловой зависимости.)

Радиальные многочлены также являются либо четными, либо нечетными в зависимости от порядка n или m :

Эти равенства легко увидеть, поскольку при нечетном (четном) m содержит только нечетные (четные) степени ρ (см. примеры ниже).

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

Рекуррентные соотношения

Полиномы Цернике удовлетворяют следующему рекуррентному соотношению, которое не зависит ни от степени, ни от азимутального порядка радиальных полиномов: [10]

Из определения видно, что и . Следующее трехчленное рекуррентное соотношение [11] позволяет тогда вычислить все остальные :

Приведенное выше соотношение особенно полезно, поскольку производную можно вычислить из двух радиальных полиномов Цернике смежной степени: [11]

Дифференциальное уравнение гипергеометрической функции Гаусса эквивалентно

Примеры

Радиальные многочлены

Первые несколько радиальных многочленов:

Полиномы Цернике

Ниже показаны первые несколько мод Цернике с различными индексами. Они нормализованы так, что: , что эквивалентно .

Приложения

Результат первых 21 полиномов Цернике (как указано выше), введенных как аберрации на плосковершинном пучке. Пучок отображается линзой, осуществляющей преобразование Фурье, интенсивность которого представлена ​​на этом рисунке

Функции являются базисом, определенным над круговой опорной областью, обычно плоскостями зрачка в классическом оптическом изображении на видимых и инфракрасных длинах волн через системы линз и зеркал конечного диаметра. Их преимуществами являются простые аналитические свойства, унаследованные от простоты радиальных функций и факторизации в радиальных и азимутальных функциях; это приводит, например, к замкнутым выражениям двумерного преобразования Фурье в терминах функций Бесселя. [12] [13] Их недостатком, в частности, если задействованы высокие n , является неравномерное распределение узловых линий по единичному диску, что вносит эффекты звона вблизи периметра , что часто приводит к попыткам определить другие ортогональные функции по круговому диску. [14] [15] [16]

В прецизионном оптическом производстве полиномы Цернике используются для характеристики ошибок более высокого порядка, наблюдаемых в интерферометрических анализах. В датчиках наклона волнового фронта, таких как Шак-Гартман , коэффициенты Цернике волнового фронта могут быть получены путем подгонки измеренных наклонов с производными полинома Цернике, усредненными по субапертурам выборки. [17] В оптометрии и офтальмологии полиномы Цернике используются для описания аберраций волнового фронта роговицы или хрусталика от идеальной сферической формы, которые приводят к ошибкам рефракции . Они также широко используются в адаптивной оптике , где их можно использовать для характеристики атмосферных искажений . Очевидными приложениями для этого являются ИК- или визуальная астрономия и спутниковые изображения .

Другое применение полиномов Цернике можно найти в расширенной теории дифракции и аберраций Нийбура–Цернике.

Полиномы Цернике широко используются в качестве базисных функций моментов изображения . Поскольку полиномы Цернике ортогональны друг другу, моменты Цернике могут представлять свойства изображения без избыточности или перекрытия информации между моментами. Хотя моменты Цернике существенно зависят от масштабирования и перемещения объекта в области интереса (ROI), их величины не зависят от угла поворота объекта. [18] Таким образом, их можно использовать для извлечения признаков из изображений, которые описывают характеристики формы объекта. Например, моменты Цернике используются в качестве дескрипторов формы для классификации доброкачественных и злокачественных масс молочной железы [19] или поверхности вибрирующих дисков. [20] Моменты Цернике также использовались для количественной оценки формы линий клеток рака остеосаркомы на уровне отдельных клеток. [21] Более того, моменты Зернике использовались для раннего выявления болезни Альцгеймера путем извлечения дискриминационной информации из изображений МРТ группы людей с болезнью Альцгеймера, с легкими когнитивными нарушениями и группы здоровых людей. [22]

Более высокие измерения

Концепция переносится в более высокие измерения D , если многочлены в декартовых координатах преобразуются в гиперсферические координаты , , умноженные на произведение многочленов Якоби угловых переменных. В измерениях угловые переменные являются сферическими гармониками , например. Линейные комбинации степеней определяют ортогональный базис, удовлетворяющий

.

(Обратите внимание, что здесь фактор включен в определение R , тогда как при нормализации он выбирается немного иначе. Это в значительной степени дело вкуса, зависящее от того, желаете ли вы сохранить целочисленный набор коэффициентов или предпочитаете более строгие формулы, если задействована ортогонализация.) Явное представление имеет вид [3]

для четных , в противном случае тождественно нулю.

Смотрите также

Ссылки

  1. ^ Зернике, Ф. (1934). «Beugungstheorie des Schneidenverfahrens und Seiner Verbesserten Form, der Phasenkontrastmethode». Физика . 1 (8): 689–704. Бибкод : 1934Phy.....1..689Z. дои : 10.1016/S0031-8914(34)80259-5.
  2. ^ Борн, Макс и Вольф, Эмиль (1999). Принципы оптики: Электромагнитная теория распространения, интерференции и дифракции света (7-е изд.). Кембридж, Великобритания: Cambridge University Press. стр. 986. ISBN 9780521642224.(см. также в Google Books)
  3. ^ ab Mathar, RJ (2009). «Базис Зернике для декартовых преобразований». Сербский астрономический журнал . 179 (179): 107–120. arXiv : 0809.2368 . Bibcode :2009SerAJ.179..107M. doi :10.2298/SAJ0979107M. S2CID  115159231.
  4. ^ Нолл, Р. Дж. (1976). «Полиномы Цернике и атмосферная турбулентность» (PDF) . J. Opt. Soc. Am . 66 (3): 207. Bibcode : 1976JOSA...66..207N. doi : 10.1364/JOSA.66.000207.
  5. ^ Thibos, LN; Applegate, RA; Schwiegerling, JT; Webb, R. (2002). «Стандарты отчетности об оптических аберрациях глаз» (PDF) . Журнал рефракционной хирургии . 18 (5): S652-60. doi :10.3928/1081-597X-20020901-30. PMID  12361175.
  6. ^ Лумис, Дж., «Компьютерная программа для анализа интерферометрических данных», Оптические интерферограммы, обработка и интерпретация, ASTM STP 666, под ред. А. Х. Гюнтера и Д. Х. Либенберга, Американское общество по испытаниям и материалам, 1978, стр. 71–86.
  7. ^ Генберг, В. Л.; Михельс, Г. Дж.; Дойл, КБ (2002). «Ортогональность полиномов Цернике». Optomechanical design and Engineering 2002. Proc SPIE. Vol. 4771. pp. 276–286. doi :10.1117/12.482169.
  8. ^ Эрик П. Гудвин; Джеймс С. Вайант (2006). Полевое руководство по интерферометрическому оптическому тестированию . стр. 25. ISBN 0-8194-6510-0.
  9. ^ Лакшминараянан, В.; Флек, Андре (2011). «Полиномы Зернике: руководство». J. Mod. Opt . 58 (7): 545–561. Bibcode :2011JMOp...58..545L. doi :10.1080/09500340.2011.554896. S2CID  120905947.
  10. ^ Хонарвар Шакибаи, Бармак (2013). «Рекурсивная формула для вычисления радиальных полиномов Цернике». Opt. Lett . 38 (14): 2487–2489. Bibcode :2013OptL...38.2487H. doi :10.1364/OL.38.002487. PMID  23939089.
  11. ^ ab Kintner, EC (1976). «О математических свойствах полиномов Цернике». Opt. Acta . 23 (8): 679–680. Bibcode :1976AcOpt..23..679K. doi :10.1080/713819334.
  12. ^ Татулли, Э. (2013). «Преобразование коэффициентов Цернике: метод на основе Фурье для масштабированных, смещенных и повернутых апертур волнового фронта». J. Opt. Soc. Am. A. 30 ( 4): 726–32. arXiv : 1302.7106 . Bibcode : 2013JOSAA..30..726T. doi : 10.1364/JOSAA.30.000726. PMID  23595334. S2CID  23491106.
  13. ^ Janssen, AJEM (2011). "Новые аналитические результаты для полиномов окружности Цернике из базового результата в теории дифракции Нийбура-Цернике". Журнал Европейского оптического общества: Rapid Publications . 6 : 11028. Bibcode :2011JEOS....6E1028J. doi : 10.2971/jeos.2011.11028 .
  14. ^ Баракат, Ричард (1980). «Оптимально сбалансированные аберрации волнового фронта для радиально-симметричных распределений амплитуд: обобщения полиномов Цернике». J. Opt. Soc. Am . 70 (6): 739–742. Bibcode : 1980JOSA...70..739B. doi : 10.1364/JOSA.70.000739.
  15. ^ Янссен, AJEM (2011). «Обобщение полиномов круга Цернике для прямых и обратных задач в теории дифракции». arXiv : 1110.2369 [math-ph].
  16. ^ Матар, Р. Дж. (2018). «Ортогональная базисная функция над единичной окружностью со свойством минимакса». arXiv : 1802.09518 [math.NA].
  17. ^ Аконди, Вьяс; Дубра, Альфредо (22 июня 2020 г.). «Средний градиент полиномов Цернике по полигонам». Optics Express . 28 (13): 18876–18886. Bibcode : 2020OExpr..2818876A. doi : 10.1364/OE.393223 . ISSN  1094-4087. PMC 7340383. PMID 32672177  . 
  18. ^ Тахмасби, А. (2010). Эффективная система диагностики опухолей молочной железы с использованием моментов Зернике . 17-я Иранская конференция по биомедицинской инженерии (ICBME'2010). Исфахан , Иран : IEEE . стр. 1–4. doi :10.1109/ICBME.2010.5704941.
  19. ^ Tahmasbi, A.; Saki, F.; Shokouhi, SB (2011). «Классификация доброкачественных и злокачественных масс на основе моментов Цернике». Компьютеры в биологии и медицине . 41 (8): 726–735. doi :10.1016/j.compbiomed.2011.06.009. PMID  21722886.
  20. ^ Rdzanek, WP (2018). «Звуковое излучение вибрирующей упруго поддерживаемой круглой пластины, встроенной в плоский экран, пересмотренное с использованием полиномов круга Цернике». J. Sound Vib . 434 : 91–125. Bibcode : 2018JSV...434...92R. doi : 10.1016/j.jsv.2018.07.035. S2CID  125512636.
  21. ^ Ализаде, Элахех; Лайонс, Саманта М; Касл, Джордан М; Прасад, Ашок (2016). «Измерение систематических изменений в форме инвазивных раковых клеток с использованием моментов Цернике». Интегративная биология . 8 (11): 1183–1193. doi :10.1039/C6IB00100A. PMID  27735002.
  22. ^ Горджи, Х. Т. и Дж. Хаддадния. «Новый метод ранней диагностики болезни Альцгеймера на основе псевдо-момента Цернике из структурной МРТ». Neuroscience 305 (2015): 361–371.

Внешние ссылки