Многочастотная голография

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

Эффект улучшения радиального разрешения системы обуславливается специфическим видом радиальной функции рассеяния. На рисунке 19 показаны несколько различных радиальных функций рассеяния системы, рассчитанных для разных значений длин волн. Функции рассеяния имеют максимум в точке z = 0 и при удалении от этой точки затухают с колебаниями, период которых зависит от длины волны.

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

Рисунок 0.19 – Функции рассеяния радиоголографической системы при разных частотах облучения

Алгоритм восстановления изображений с использованием многочастотной информации достаточно прост. Его принцип заключается в формировании набора одночастотных распределений поля в плоскости объекта любым из методов и последующем их сложении.

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

Функция рассеяния для вписанной круговой апертуры определяется выражением

.

Это выражение получено, если восстановление изображения производится при помощи преобразования вида :

.

Мы будем использовать несколько другую форму записи преобразования Френеля:

.

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

.

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

В пределе, при увеличении числа частот, сумму в выражении можно заменить интегралом:

Вычислить интеграл непосредственно не представляется возможным из-за наличия в нем функции sinc(). Однако видно, что с учетом использованных приближений ( и ) аргумент функции sinc() будет величиной, значительно меньшей единицы. В этом случае значение функции sinc() будет примерно равно единице.

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

Запишем выражение для функции рассеяния с учетом этих соображений:

Запишем граничные значения волновых чисел в следующем виде:

.

Выражение в этом случае запишется в следующем виде:

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

.

где  – полоса длин волн излучения, f – полоса частот облучения, c – скорость распространения волн в среде.

В многочастотном случае радиальное разрешение может быть сопоставимо с азимутальным. Допустим, пусть в предыдущем примере (размер апертуры 0,5 м, расстояние до нее 5 м) используются волны c длинами от 1 до 1,1 см, что соответствует частотам от 27 ГГц до 30 ГГц) радиальное разрешение составит 0,1 м.

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

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

Блок-схема алгоритма восстановления изображений с использованием многочастотной информации представлена на рисунке 20.

Рисунок 0.20 – Блок-схема алгоритма восстановления изображения с использованием многочастотной информации