Фазовая декорреляция в дискретных когерентных ансамблях для уменьшения спекл-шума

Бюллетень науки и практики / Bulletin of Science and Practice

УДК 535.42: 535.317                              

Проблема спекл-шума является фундаментальным ограничением для всех систем активного оптического зондирования, использующих лазерные источники [1].

В классической постановке для полностью развитой спекл-картины, формируемой большим числом независимых рассеивателей, распределение интенсивности подчиняется негауссовой статистике с контрастом C = а/ /( Г) =1. Однако в современных приложениях, связанных с машинным зрением, голографической микроскопией и адаптивной оптикой, все чаще используются решетчатые ансамбли источников — например, матрицы вертикально-излучающих лазеров (VCSEL) с числом эмиттеров N от 16 до 256.

Переход от непрерывного распределения рассеивателей к дискретному ансамблю с конечным N коренным образом меняет статистику поля. Как показано в работах [2, 3], при N ^ 100 центральная предельная теорема выполняется плохо, и распределение интенсивности описывается обобщенным распределением Райса с параметром формы, зависящим от числа степеней свободы.

В частности, при синфазном сложении полей (что характерно для работы VCSEL-матриц в режиме максимальной мощности) контраст может достигать аномальных значений C >1, что делает невозможным применение стандартных алгоритмов фильтрации.

Существующие методы подавления спеклов делятся на три класса: (а) пространственное усреднение за счет движения диффузора [4-5]; (б) спектральное усреднение за счет перестройки длины волны [6-8]; (в) поляризационное усреднение.

Однако все они либо требуют подвижных механических частей, либо ограничены по скорости перестройки. Альтернативой является управление начальными фазами отдельных эмиттеров. Теоретические оценки показывают [9-10], что перебор M различных фазовых масок за время экспозиции эквивалентен увеличению числа эффективных независимых реализаций поля. Несмотря на интуитивную ясность данного подхода, до сих пор отсутствует строгая замкнутая теория, связывающая число фазовых состояний M , пространственную геометрию решетки (период Λ, апертуру a ) и результирующий контраст C .

Цель статьи — построить аналитическую модель, которая позволяет рассчитать C ( M ) и определить критерий оптимальности M _опт без привлечения эмпирических подгоночных параметров.

Теоретическая модель и вывод основного выражения для контраста. Рассмотрим плоскость источника z =0, в которой расположены N точечных монохроматических излучателей с комплексными амплитудами A _n и начальными фазами φ _n . Пространственная конфигурация задается векторами г = (х_п,у_п). Тогда комплексное поле в плоскости источника имеет вид:

N

U(r; 0) = У An 8(r — rn)ei(l>n.

n=1

Распространение света в свободном пространстве на расстояние d описывается интегралом Френеля (приближение параксиальной оптики):

U(.P)=-m\ U(r; 0) exp|^|p —r|2]dr, lAd              2d

где р = (% ,Г — координаты в плоскости наблюдения, k=2п / к . Подставляя (1) в (2), получаем:

eikd у—,                                   ^

U(P) =HdyAnei^n exp[l^n(p)],^n(p) =^|P-rn|2. n=1

Интенсивность  I(р )=| U(р )|2  раскладывается на сумму квадратов амплитуд и интерференционные члены:

N      ²

l(p) = yiAnl + У AnAmexp[l(фn — Фт) + Wn — Vm)] n=1      n^m

Рассмотрим процесс усреднения по всему множеству масок. Предположим, что мы

можем генерировать M независимых реализаций вектора начальных фаз φ ^{( m )} = ( φ ^{( m )} Усредненная интенсивность:

,...

_{, φ N} ( m ) _).

м

<1(Р»м=1У1(т,(р')

т=1

Второе слагаемое в (4) содержит фактор:

м

Xnm^^ exP[ 1(фППП)) - ФтЛ m m-1

При статистической независимости и равномерном распределении фаз на [0,2 π ) для n≠m имеем χ → 0 при M →∞. Для конечного M остаточная корреляция определяется функцией:

г. x(M} = [

sln^/M} . n/M

]'■

что следует из дискретного суммирования по равномерной сетке. Контраст спеклов определяется как:

OM)^—^

\4M

После усреднения произведения двух интенсивностей и учёта статистической независимости фаз для различных масок, получаем основное выражение:

С 2_{(М) =} 1.^П=Лап^ _{(   1}}

^{G (M)} M  (%^N = ₁\ A „ I²')²    ( M) ^M "ростр'

где пространственный фактор когерентности:

Для регулярной прямоугольной решетки с периодом Л и апертурой а = VN Л , при условии a»Xd/A, фактор '_простр аппроксимируется экспоненциально:

^М простр

ехр (

2я ² а ² Л ²

A ² d ²

Подставляя (11) в (9) и полагая равные амплитуды An=A, получаем компактную формулу:

1/1              аЛ

c(M'N'a) = \мы + (1-м)е^2 а 'а = й-

Требуя C≤ε, из (12) находим:

м опт

= тах

{Ng2'

- _е -2л²а²

_— е —2^²а²

}

Анализ (13) выявляет два предельных случая:

При а«1/(2п) (сильная когерентность) М_опт ~ 1/(Ne²)

При а»1/(2п) (слабая когерентность) М_опт ~ 1/е².

Переходный режим определяется условием α=1/(2π), которое совпадает с критерием ван Циттерта–Цернике для дискретных источников. Спектральные свойства, границы применимости и сравнение с механическими методами. Пространственно-частотный спектр поля (3) имеет вид:

^^

u(f) = F{U(p)} =

^Апе1фп п

2я

ехр[ — i~^rn

• f] • const.

Тогда ФПМ системы (Функция передачи модуляции - критерий качества, показывающий, насколько хорошо оптическая система (или в нашем случае — процесс распространения и усреднения поля) передает контраст объекта при различных пространственных частотах), усреднённая по M фазовым маскам, выражается как:

(MW))M =-1У|Лп|2+-1У ЛМ)ехр[-12НГп-г,) •f]. п            п^1

При M→∞(χ→0) ФПМ стремится к постоянной величине, что соответствует полностью некогерентному источнику. При конечном M сохраняется узкий пик на нулевой частоте, амплитуда которого пропорциональна x(M)~1/M2. Относительная доля когерентной энергии:

PKorxx(MyN -N/M2.

Это означает, что подавление когерентной составляющей происходит быстрее, чем при механическом усреднении (где Р ког ~ 1/VM), и не влияет на высокочастотные компоненты, что важно для сохранения резкости краев. Предложенная теория основана на нескольких допущениях, которые необходимо учитывать при практической реализации.

Параксиальное приближение справедливо при d»a2/k Для типичных параметров (a=1 мм, Х=0.85 мкм) это условие выполняется при d»1.2 м. В ближней зоне требуется использование интеграла Рэлея–Зоммерфельда, однако структура формул сохраняется с заменой Ψ n на гиперболическую фазу.

• 1.    Независимость фаз предполагает отсутствие перекрестных тепловых связей между соседними эмиттерами VCSEL. В реальности такие связи уменьшают эффективное число степеней свободы Neff

• 2.    Временная когерентность не учитывалась. Ширина линии VCSEL (Av~50 МГц) даёт длину когерентности Lc~6 м. Для d«Lc временными эффектами можно пренебречь; в противном случае в (9) добавляется фактор ехр( — lt_n — t_m|/r_c).

• 3.    Негауссовость статистики для малых N<10 может приводить к отклонениям от формулы (9). Однако большинство практических систем имеют N≥16, где корреляционный подход второго порядка даёт хорошую точность.

Для вращающегося диффузора эффективное число независимых реализаций M мех ≈Ω T /(2 π ), где Ω — угловая скорость, T — время экспозиции. При T =1 мс и типичных скоростях Q~10 3 об/мин, M мех ~10. В то же время электронное переключение фаз (с временем ~100 нс на состояние) позволяет за 1 мс достичь M =10⁴. Теоретический выигрыш в контрасте (при одинаковой геометрии) составляет:

9мех - [Еээл - 71000 - З1.браз.

С эл   \ Ммех                   ^F

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

Заключение

В настоящей работе представлена аналитическая теория уменьшения спекл-шума путём управления начальными фазами дискретного ансамбля когерентных источников. Основные результаты:

• 1.    Выведено замкнутое выражение для контраста спеклов C(M, N, α) — формула (12), объединяющая статистические (M, N) и волно-геометрические (α=aΛ/λd) параметры системы.

• 2.    Доказано, что при M→∞ контраст спадает как 1/M в слабокогерентном режиме и как 1/MN в режиме сильной когерентности. Критерий перехода между режимами имеет вид α=1/(2π).

• 3.    Получен аналитический критерий оптимальности Mопт (13), позволяющий расчётно определить минимальное число фазовых состояний для достижения заданного порога контраста ε. Для типичных параметров (N=64, α=0.5) достаточно M=64, чтобы обеспечить C<0.08.

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

• 5.    Теоретически обосновано, что электронное управление фазами обеспечивает более чем на порядок более эффективное подавление шума по сравнению с механическими деспеклерами при равном времени экспозиции.

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