Моделирование алгоритмов обработки сигналов двухпроходного авиационного интерферометрического РСА
П.Е. Шимкин, А.И. Баскаков
Национальный исследовательский университет «МЭИ»
Аннотация
В работе проведено имитационное моделирование алгоритмов интерферометрической обработки сигналов в пакете прикладных программ MATLAB и обоснован выбор основных параметров одного из возможных вариантов построения интерферометрического РСА.
Ключевые слова: интерферометрический РСА с «мягкой базой», разность фаз, трехмерное изображение, рельеф поверхности.
Abstract
The simulation algorithms interferometric signal processing in the application package MATLAB was made in the paper and choice of the basic parameters of one of the possible variants of interferometric SAR was justified.
Key words: interferometric SAR with a "soft base", the phase difference, three-dimensional images, surface relief.
Введение
В последнее время, развитие теоретических принципов и технических возможностей построения радиолокаторов с синтезируемой апертурой антенны (РСА) позволило получать радиолокационные изображения (РЛИ) с разрешающей способностью порядка единиц метров и выше. Но дальнейшее повышение геоинформативности авиационных и космических средств дистанционного зондирования Земли (ДЗЗ) связано с развитием теории и техники получения детального рельефа местности и формированием трехмерного изображения [1]. Поэтому, на сегодняшний день, основной задачей является создание перспективных комплексов, способных с высокой точностью и высокой разрешающей способностью формировать кроме привычных двумерных РЛИ также и трехмерные изображения земной поверхности. Одним из решений построения такого комплекса является двухпроходный интерферометрический РСА (ИРСА). Он может быть выполнен как в авиационном, так и в космическом варианте построения.
Такой комплекс является сложной системой, поэтому реальная проверка алгоритмов его работы затруднена из-за финансовых или технических возможностей. Также, реальные эксперименты могут выполняться только в небольшом диапазоне вариаций условий наблюдения, что недостаточно для полного анализа системы. В связи с этим, создание имитационной модели и проведение компьютерного моделирования являются одним из основных и эффективных методов изучения такой сложной системы.
Геометрия визирования
Моделируется двухпроходный авиационный ИРСА, геометрия визирования которого изображена на рис. 1, где синим цветом условно указан первый проход носителя, красным - второй. Наблюдение за исследуемой поверхностью происходит в режиме бокового обзора с помощью одного носителя РСА, но за два прохода, так называемый режим «мягкой базы». В двух проходах носитель РСА движется вдоль координаты с постоянной скоростью , формируя синтезированную апертуру , но с различными положениями по высоте и направлению поперечному движению. На рис. 1 это проявляется в наличии сдвигов во втором проходе и . В результате образуется интерферометрическая пара антенн, разделенных базой интерферометра . Следует отметить, что в общем случае наличие сдвигов и приводит к наклону базы к горизонтали, обозначим этот угол . В большинстве случаев принимает значения (строго горизонтальное смещение) и (строго вертикальное смещение). В статье рассматривается случай строго вертикального смещения. На рис. 1 - размеры элемента разрешения по дальности и по азимуту, и - координаты центра кадра; и - размеры кадра по соответствующим координатам.
Рис. 1. Геометрия визирования двухпроходного ИРСА.
Структура модели
Процесс имитационного моделирования может быть представлен в виде блок-схемы (рис. 2), где ЦМР - цифровая модель рельефа, РЛИ - радиолокационное изображение, ИРФ - интерферометрическая разность фаз.
Рис. 2. Структурная схема имитационного моделирования ИРСА
Цифровая модель рельефа
Происходит формирование радиолокационного рельефа исследуемой поверхности согласно феноменологической модели [1]. Каждый элемент разрешения на поверхности Земли представляется набором отдельных парциальных отражателей, распределенных по нормальному закону, на которые накладываются известные из экспериментальных результатов условия рассеяния.
В качестве исходных данных для имитационной модели двухпроходного интерферометра возьмем модель рельефа, изображенную на рис. 3. Согласно рис. 3 центр кадра имеет координаты м, м. Размеры кадра определяются пятном облучения диаграммы направленности антенны на поверхности.
Рис. 3. Цифровая модель рельефа поверхности
В качестве модели антенны используем фазированную антенную решетку (ФАР), вид и ДНА которой приведены на рис. 4.
Рис. 4. Внешний вид и ДНА
Ширина луча в обоих плоскостях составляет примерно . При дальностях порядка 7,5 км и высоте полета 5 км получим кадр размером м. Поэтому, элементарные отражатели, формирующие распределенную цель, будут расположены на поверхности в интервале горизонтальных дальностей м и азимута м. Тип поверхности представляет собой сухую почву, по виду - овражно-холмистая с среднеквадратической ординатой мелкой шероховатости см и перепадами неровностей рельефа от до 80 м.
Основной параметр, характеризующий распределенную цель - удельная эффективная поверхность рассеяния (УЭПР). В настоящее время для приближенности к реальному миру разработчики используют экспериментально полученные натурные модели УЭПР, хорошо зарекомендовавшие себя на протяжении многих лет. Одной из таких моделей является, разработанная научно-исследовательским институтом Georgia Tech (Georgia Tech Research Institute, Atlanta, Georgia, USA) модель УЭПР для различного типа поверхностей [2]. Данная модель учитывает среднеквадратическую ординату мелкой шероховатости поверхности , угол падения на поверхность , отсчитываемый от вертикали и длину волны , и имеет следующий вид [3] (1):
,(1)
где - коэффициенты эмпирической модели. В [3, 4] приведены значения для этих констант в частотном диапазоне частот от 3 до 95 ГГц для некоторых типов поверхностей.
На рис. 5 представлена зависимость УЭПР от угла падения для сухой почвы на частоте 10 ГГц.
Рис. 5. Зависимость УЭПР от угла падения для сухой почвы
Параметры двухпроходного ИРСА
Рабочая частота передатчика РСА ГГц (соответствующая длина волны см). Например, для получения разрешения по координате наклонной дальности м и, следственно, горизонтальной дальности м используется зондирующий ЛЧМ сигнал с шириной спектра МГц.
Разрешение вдоль координаты , линии пути (азимуту) достигается за счет применения техники синтезированной апертуры. Например, для достижения разрешающей способности по азимуту примерно равной разрешающей способности по горизонтальной дальности необходимо синтезировать апертуру длиной м.
Длительность импульса можно выбрать из скважности . Пусть , тогда с учетом того что мкс, мкс.
Период повторения зондирующих импульсов самолетного РСА выбирается из условия однозначности формирования РЛИ как по азимуту, так и по дальности (2):
,(2)
где , скорость полета носителя порядка 250 м/с, типичный размер антенны авиационной РЛС м, мкс - время работы антенного переключателя. Тогда: мкс. Пусть мкс.
Выбор оптимальной базы интерферометра
Из рис. 1 может быть получено соотношение между параметрами отклонения базы , угла визирования и разности фаз :
(3)
Из (3) видно, что что ордината элемента разрешения z является функцией параметров: высоты носителя , наклонной дальности азимутального угла , размера базы и разности фаз . Поэтому общую погрешность измерения высоты элемента разрешения, с предположением о некоррелированности отдельных составляющих, можно записать в виде [5]:
,(4)
где - дисперсия оценки ординаты элемента разрешения из-за погрешности оценки разности фаз; - дисперсия оценки ординаты элемента разрешения из-за погрешности задания угла наклона базы; - дисперсия оценки ординаты элемента разрешения из-за погрешности измерения высоты орбиты; - дисперсия оценки ординаты элемента разрешения из-за погрешности измерения расстояния до элемента разрешения; - дисперсия оценки ординаты элемента разрешения из-за погрешности неточного знания размера базы.
Выбор наиболее выгодного режима интерферометрического зондирования сводится к определению параметров съемки, минимизирующих вышеперечисленные погрешности. Например, погрешность, связанная с погрешностью оценки разности фаз, определяется следующим выражением (5):
(5)
Так как потенциальные точностные характеристики определяются только через , то наибольший интерес представляет анализ погрешности, связанной именно с оценкой ИРФ .
Согласно [6], оценку ИРФ можно провести использую метод «оценки максимального правдоподобия (ОМП)». Тогда, среднеквадратическое отклонение ИРФ можно оценить следующим выражением (6):
(6)
где - кратность некогерентного накопления, - коэффициент корреляции.
Интерферометрическая когерентность или коэффициент корреляции между двумя снимками, составляющими интерферометрическую пару, в общем случае состоит из четырех составляющих [7]: декорреляция за счет пространственного разнесения двух антенн - ; декорреляция, вследствие наличия мелкой шероховатости поверхности - ; временная декорреляция, связанная с физическими изменениями на поверхности между двумя наблюдениями - ; декорреляция за счет теплового шума в системе РСА - . Для простоты .
Результирующий коэффициент корреляции представляется в виде произведения вышеперечисленных составляющих:
.(7)
Получены аналитические выражения каждого из факторов. При расчетах было принято предположение о представлении реальной поверхности в виде распределённой радиолокационной цели, для которой используется феноменологическая модель из набора независимых парциальных отражателей внутри элемента разрешения. Также было использовано описание плотности распределения ординат парциальных отражателей внутри элемента разрешения в виде нормального закона. Результаты расчетов каждой из составляющих (7):
;(8)
;(9)
;(10)
;(11)
где - разрешающая способность по наклонной дальности; - среднеквадратическая ордината мелкой шероховатости поверхности; - разрешающая способность по азимуту; - отношение сигнал/шум.
Получив результирующий коэффициент корреляции в виде произведения указанных составляющих и, подставив их в (6) можем получить выражение для СКО оценки ИРФ. Затем подставляя (6) в (5) получим зависимость СКО оценки ординат рельефа через СКО оценки ИРФ. В связи с громоздкостью этих выражений, их полное отображение в документе не представляется возможным.
Для выбранных параметров комплекса ИРСА и моделируемой поверхности, имеем следующую зависимость потенциальной точности оценивания рельефа поверхности от базы интерферометра (рис. 6).
Рис. 6. СКО оценивания рельефа в зависимости от размера базы при параметрах: МГц; ; ; км; дБ; ; см.
Из рис. 6 видно, что оптимальным является смещение носителя во втором проходе на 7,5 м, при котором потенциальное СКО оценки рельефа примерно 2 м.
Моделирование траекторного сигнала
двухпроходный авиационный интерферометрический радиолокатор
Для начала необходимо определить некоторые параметры, характеризующие сигналы на входе процессора регистрации и обработки траекторного сигнала РСА [8].
Имеем:
см и ГГц - длина волны РСА и несущая частота зондирующего сигнала;
мкс - длительность зондирующего импульса;
МГц - ширина спектра зондирующего ЛЧМ сигнала;
мкс и кГц - частота и период повторения импульсов;
МГц; нс - частота и период дискретизации принятых сигналов;
- число отсчетов на длительности импульса;
м/с - скорость полета носителя;
с - время синтеза;
- число синтезируемых импульсов;
м - шаг дискретизации по наклонной дальности (масштаб голограммы и РЛИ);
м и м - наклонные дальности первого и последнего отсчета сигналов в строке радиолокационных данных;
- число отсчетов в стробе приема по дальности.
Излучаемые, а также принимаемые двумерные сигналы рассматривают раздельно в координатах «быстрого времени» (по наклонной дальности) и «медленного времени» по азимуту , которое можно выразить через номер зондирующего импульса . Излучаемый сигнал представляет собой последовательность зондирующих импульсов, каждый из которых характеризуется в «быстром времени» выражением (12):
,(12)
где при .
Принимаемый комплексный видеосигнал на выходе фазового детектора приемника является суперпозицией сигналов, отраженных от элементарных отражателей, расположенных в пределах кадра формируемой поверхности (13):
(13)
где - амплитудный множитель, нормированный на излученную передатчиком мощность и коэффициент усиления антенны; -УЭПР отражателя; ; - текущая наклонная дальность до цели; - номер отсчета дальности для переднего фронта принятого сигнала от цели. В общем случае зависит от положения РСА по азимуту (номера строки дальности ), что вызывает эффект «миграции» дальности.
Обработка траекторного сигнала и формирование РЛИ
Обработка траекторного сигнала РСА происходит по алгоритму прямой свертки принимаемого сигнала с комплексно сопряженной двумерной дискретной опорной функцией.
Для определения вида опорной функции, аппроксимируют временной закон текущей наклонной дальности квадратичным членом [9]. Азимутальные отсчеты дальности и фазы определяются формулами:
(14)
Тогда, двумерная опорная функция будет иметь вид (15) [8]:
(15)
Прямая свертка двумерных дискретных сигналов выполняется в два этапа - сначала производится свертка по дальности, а затем по азимуту. Наиболее распространенным и одним из самых производительных алгоритмов вычисления свертки является алгоритм быстрой свертки. Преимущества быстрой свертки реализуются при использовании для дискретного преобразования Фурье алгоритма быстрого преобразования Фурье (БПФ, Fast Fourier Transform - FFT). Таким образом, свертки по дальности и азимуту (16):