- •Основные понятия курса. Оптическая и неоптическая голография
- •Что такое изображение
- •Методы восстановления изображений
- •Методы реконструкции изображений
- •Другие методы цифровой обработки изображений
- •Оптическая голография. Регистрация интерференционной картины.
- •Оптическая схема получения голограммы.
- •Неоптическая голография
- •Математический аппарат решения задач восстановления и реконструкции изображений
- •Дельта-функция
- •Свойства дельта-функции
- •Преобразование Фурье. Теорема о свёртке
- •Линейные системы. Импульсный отклик линейной системы
- •Прямые и обратные задачи. Уравнение Фредгольма
- •Решение уравнения типа свёртки. Частотная характеристика
- •Корректность решения обратной задачи. Существования решения
- •Единственность решения на примере уравнения типа свертки
- •Устойчивость решения
- •Регуляризация решени обратных задач
- •Регуляризация решения. Метод регуляризации Тихонова
- •Регуляризация решения уравнения типа свертки
- •Фильтр Тихонова. Невязка
- •Оптимальный фильтр Винера
- •Управляемая линейная фильтрация. Фильтр Бэйкуса-Гильберта
- •Гомоморфная фильтрация
- •Метод неопределенных коэффициентов
- •Пример решения обратной задачи
- •Коррекция искажений, вызванных равномерным прямолинейным движением объекта
- •Коррекция искажений, вызванных равномерным прямолинейным движением объекта. Учет граничных условий
- •Разрешающая способность систем формирования изображений
- •Понятие о разрешающей способности
- •Теоретическая оценка разрешающей способности на примере анализатора спектра
- •Представление Релея для монохроматических волн
- •Представление Релея для немонохроматических волн
- •Двойной физический смысл пространственной частоты
- •Частотная характеристика свободного пространства
- •Угловой спектр сферической волны
- •Импульсный отклик свободного пространства
- •Восстановление радиоголографических изображений
- •Алгоритм восстановления изображений в частотной области
- •Восстановление изображений в приближении Френеля
- •Азимутальное разрешение радиоголографической системы
- •Синтез апертуры сканированием одной антенной
- •Синтез апертуры сканирования двумя антеннами
- •Синтез радиоголограмм динамических объектов
- •Разрешающая способность в радиальном направлении
- •Многочастотная голография
- •Основы томографии
- •Прохождение плоскопараллельного пучка через среду с поглощением
- •Преобразование Радона
- •Преобразование Радона точечного объекта
- •Теорема о центральном сечении
- •Обратное преобразование Радона
- •Алгоритм обратного проецирования
- •Вычисление обратного преобразования Радона
- •Итерационные алгоритмы решения обратных задач
- •Понятие об итерационных алгоритмах решения обратных задач
- •Итерационные алгоритмы с ограничениями
- •Итерационное уравнение
- •Ряд Неймана
- •Итерационный оператор для уравнения типа свертки
Угловой спектр сферической волны
Выражение для комплексной амплитуды сферической волны
Рассмотрим, как разлагается по плоским волнам сферическая волна, комплексная амплитуда которой описывается выражением
,
где .
Такая волна является моделью, описывающей поле точечного источника излучения, и поэтому ее разложение важно для проведения некоторых дальнейших выводов. Уравнение колебаний плоской волны можно получить непосредственно из решения волнового уравнения для точечного источника колебаний.
В плоскости z = 0 распределение комплексных амплитуд будет описываться соотношением
,
где .
Вычисление углового спектра сферической волны
Подставляя выражение в соотношение , можно получить искомый угловой спектр. При этом необходимо вычислить интеграл
,
что можно сделать при помощи замены переменных , , и . После проведения вычислений можно получить, что
.
Полученное выражение представляет собой угловой спектр сферической волны. Анализируя , можно отметить, что при малых по сравнению с k значениях u1 и u2, что соответствует малым углам распространения плоских волн, спектр сферической волны практически постоянен. При увеличении углов распространения и приближении u1 и u2 к k спектр начинает возрастать и переходит в область неоднородных волн.
Подставляя выражение для углового спектра в представление Релея , получим искомое разложение:
Импульсный отклик свободного пространства
Выражение для импульсного отклика свободного пространства
Продолжая аналогию колебательных и волновых систем, выведем выражение для импульсного отклика свободного пространства.
В колебательных системах импульсным откликом называется реакция системы на единичное входное воздействие, или другими словами, импульсный отклик системы – это сигнал на выходе системы при подаче на ее вход бесконечно короткого импульса. В теории линейных систем показано, что импульсный отклик линейной системы и ее частотная характеристика связаны парой преобразований Фурье. Применим это правило к нашей волновой линейной системе и получим выражение для импульсного отклика как обратное преобразование Фурье от уже известной частотной характеристики:
.
Интеграл можно вычислить при помощи соотношения . Дифференцируя по z, получим
.
Сравнивая выражения и , можно получить, что
.
Таким образом, мы получили точное выражение для импульсного отклика свободного пространства.
Физический смысл импульсного отклика свободного пространства заключается в следующем. Представим некоторую монохроматическую волну, распределение комплексной амплитуды поля которой в плоскости z=0 представляет собой единичный импульс (дельта-функцию). Выражение представляет собой распределение комплексной амплитуды этого поля некоторой волны в плоскости z, которое образуется при прохождении этой волной свободного пространства длины z. Следует отметить, что поле такой волны физически нереализуемо.
Аналогично выражению , выражающему распределение комплексной амплитуды поля в некотором сечении через частотную характеристику, можно получить выражение для комплексной амплитуды поля через импульсный отклик свободного пространства:
.
Приближение Френеля. Принцип Гюйгенса-Френеля
Анализируя выражение , можно получить, что для больших (по сравнению с длиной волны) расстояний z импульсная характеристика записывается в виде
при .
Если предположить, что поперечные размеры плоскости, на которой анализируется комплексная амплитуда, также малы по сравнению с расстоянием до этой плоскости, выражение для импульсного отклика еще более упрощается:
при и .
Условия и являются условиями зоны Френеля, таким образом, приближенное выражение для импульсного отклика свободного пространства справедливо в этой зоне.
Вместе с тем можно придти к выводу, что в данном приближении импульсный отклик свободного пространства представляет собой не что иное, как сферическую волну (из сравнения с выражением ). Таким образом, мы приходим к принципу Гюйгенса-Френеля – фронт распространяющейся волны точно такой же, как если бы он был сформирован точечными источниками с амплитудами и фазами колебаний, соответствующими колебаниям поля на некотором другом волновом фронте.