Методы радиолокационной интерферометрии в исследовании характеристик земных покровов тема автореферата и диссертации по физике, 01.04.03 ВАК РФ

На правах рукописи

ЗАХАРОВА Людмила Николаевна

МЕТОДЫ РАДИОЛОКАЦИОННОЙ ИНТЕРФЕРОМЕТРИИ К ИССЛЕДОВАНИИ ХАРАКТЕРИСТИК ЗЕМНЫХ ПОКРОВОВ

4052092

Специальность 01.04.03 — «Радиофизика»

АВТОРЕФЕРАТ диссертации на соискание учёной степени кандидата физико-математических наук

1 1 АВГ 2011

Фрязино - 2011

4852092

Работа выполнена в Учреждении Российской академии наук Институте радиотехники и электроники им. В. А. Котельникова РАН (Фрязинский филиал)

Научный руководитель: кандидат технических наук

Захаров Александр Иванович

Официальные оппоненты: доктор технических наук, профессор

Неронский Леон Богуславович

кандидат технических наук Павельсв Александр Геннадьевич

Ведущая организация: Учреждение РАН Институт

космических исследований РАН

Защита диссертации состоится «14» октября 2011 года, в 10-00 на заседании диссертационного совета Д 002.231.02 при Учреждении Российской академии наук Институте радиотехники и электроники им. В.А.Котельникова РАН по адресу: 125009, Москва, ул. Моховая, д.11,к.7.

С диссертацией можно ознакомиться в библиотеке ИРЭ им. В.А. Котельникова РАН.

Автореферат разослан «15» и>о/иу 2011 года.

Учёный секретарь диссертационного совета, доктор физико-математических наук

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

Актуальность исследований

Исследования Земли с помощью радиолокатора с синтезированной апертурой, расположенного на борту самолета или спутника — это современный высокопроизводительный метод дистанционного зондирования, позволяющий получать двумерные распределения интенсивности отраженного сигнала, а также изучать отражательные свойства и структуры рельефа поверхности. Преимуществами радиолокаторов перед оптическими сенсорами является возможность съёмки независимо от солнечного освещения и от наличия/отсутствия облачности. Прием отраженного сигнала на две антенны на борту носителя (или же, в схеме с одной антенной, съемка выбранного участка поверхности с близких траекторий носителя) и дальнейшая совместная обработка позволяют при помощи анализа фазовой составляющей эхо-сигнала добавить третье измерение к двумерным радиолокационным изображениям, выявляя, в зависимости от поставленной задачи и начальных условий, относительные высоты рельефа, степень декорреляции или изменение положения отражающей поверхности.

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

• извлечение информации о рельефе поверхности Земли;

• детектирование площадных подвижек поверхности;

• обнаружение изменений поверхности за время между съёмками;

• оценку скорости движущихся целей.

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

Цели и задачи работы

Целью диссертационной работы является исследование возможностей радиолокационной интерферометрии в изучении характеристик земных покровов при дистанционном зондировании Земли из космоса, демонстрация существующих методик обработки, а также их развитие: предложение модификации одного из методов двумерной развёртки фазы на интерферограмме, представление нового метода классификации естественных покровов на основе данных поляриметрической интерферометрии, введение и обоснование нового способа регистрации полной матрицы рассеяния при интерферомет-рической схеме съёмки. Для достижения поставленной цели выполнены следующие этапы:

1. Проведён анализ геометрических характеристик интерферометричс-ской съёмки, рассмотрены различные варианты организации интер-ферометрической съёмки, в том числе, конфигурации космических аппаратов.

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

3. Показаны измеряемые величины и приведены точности их оценки.

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

5. Изложены основы поляриметрической интерферометрии и проанализированы существующие на текущий момент её приложения.

6. Введён в рассмотрение новый метод классификации естественных покровов по поляриметрико-интерферометрическим данным на основе анализа поляриметрических сигнатур когерентности.

7. Предложен метод квази-поляриметрической съёмки с повторяющихся орбит, проведено моделирование искажений получаемой таким образом матрицы рассеяния, обработаны экспериментальные данные двух частотных диапазонов, показана работоспособность метода, проведён сравнительный анализ источников декорреляции.

Научная новизна

• Проведено оригинальное сравнительное исследование различных методов двумерной развёртки фазы на радиолокационных интерфе-рограммах.

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

• Предложен метод построения поляриметрических сигнатур когерентности и способ классификации на его основе. Проведён сравнительный анализ возможностей такой классификаций для данных Ь- и С-диапазона, а также для различных интервалов между съёмками.

• Предложен метод получения полной матрицы за два прохода носителя в схеме интерферометрической съёмки; промоделированы возможные искажения и предложен способ их коррекции.

Достоверность

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

Практическая значимость

Практическая значимость работы заключается в том, что:

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

• введённые поляриметрические сигнатуры интерферометрической когерентности являются новым признаком, позволяющим выполнять классификацию земных покровов на основе совместного анализа поляриметрической и интерферометрической информации;

• предложенный метод формирования полной матрицы рассеяния по частичным поляризационным измерениям расширяет возможности поляриметрических систем наблюдения поверхности Земли.

Положения, выносимые на защиту

1. Система критериев, обеспечивающая эффективное использование и сравнение методов устранения 2я-неоднозначности на интерферо-граммах.

2. Модифицированный метод вычетов, устраняющий 2л-неоднозначность на интерферограмме при условии пересечённого рельефа на снимаемой местности.

3. Поляриметрические сигнатуры интерферометрической когерентности как новый классификационный признак типов земных покровов.

4. Метод формирования полной матрицы рассеяния зондируемых сред по частичным поляризационным измерениям, выполняемым в двух-проходной интерферометрической схеме съёмки.

Апробация работы

Результаты работы, изложенные в диссертации, обсуждались и докладывались на научных семинарах Фрязинского филиала ИРЭ им. В.А.Котельникова РАН, конкурсах молодых учёных им. И.В. Анисимкина в ИРЭ им. В.А.Котельникова РАН, а также на 13 отечественных и 11 международных научных конференциях и симпозиумах в 2003-2010 гг.

Публикации

Основные результаты были опубликованы в 35 работах: из них 8 статей, в том числе 7 статей в журналах из рекомендованного перечня ВАК, 27 докладов на отечественных и международных конференциях.

Структура и объём диссертации

Диссертация состоит из введения, шести глав, заключения и списка литературы. Работа изложена на 213 страницах, содержит 51 рисунок, 6 таблиц и библиографический список из 182 наименований.

СОДЕРЖАНИЕ РАБОТЫ

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

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

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

В главе 3 подробно рассмотрен один из наиболее сложных при реализации этапов интерферометрической обработки — устранение фазовой 2л-неоднозначности на двумерном изображении («развёртка фазы»).

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

• Самая вычислительно ёмкая операция помимо интегрирования

• Независимость от фрагментации

• Тип маски

• Необходимость предварительная фильтрация

• Индивидуальные особенности

• Положительные стороны

• Отрицательные стороны

• Возможные модификации

Для сравнения были выбраны следующие локальные методы:

1. Метод вычетов

2. Потоковый метод разворота фазы

3. Частотный метод анализа фазовой неоднозначности интерферо-граммы

4. Метод восстановленных значений

5. Метод, основанный на локальной фильтрации интерферограммы

6. Метод определения граничных линий

Для одного из локальных методов, участвующих в сравнении (метода вычетов), была предложена модификация. Суть исходного метода состоит в следующем: на изображении разности фаз (интерферограмме), вообще говоря, присутствуют «особые» точки, названные, по аналогии с особыми точками в теории функций комплексного переменного, «вычетами». Замкнутый путь интегрирования, окружающий такую точку, приводит к ненулевому результату, что нарушает условие однозначности интегрирования по путям, проходящим через неё. В зависимости от знака результата интегрирования по минимальному замкнутому пути вокруг точки она называется «положительным» или «отрицательным» вычетом. В исходном методе ближайшие друг к другу отрицательный и положительный вычеты предлагается соединять «разрезом»: линией, которую не должны пересекать пути интегрирования. После построения всех разрезов интерферограмма становится однозначно интегрируемой. Модификация, предложенная автором диссертации, заключается в дополнительном использовании информации о степени декорреля-ции (интерферометрической когерентности), что позволяет более рационально использовать строить маску точек, запрещённых для путей интегрирования. Помимо вычисления вычетов, предлагается строить маску бракованных точек на основе анализа значений амплитуды комплексной интерферограммы, показывающей меру декорреляции отражённого сигнала (интерферомет-рическую когерентность). Затем подсчитывается алгебраическая сумма значений вычетов в связных компонентах этой маски («заряд» компоненты), и, при необходимости, компоненты с зарядами противоположных знаков соединяются друг с другом линиями разрезов, с тем чтобы итоговый баланс зарядов стал нулевым. Таким образом, маска становится распределённой по

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

Для сравнения глобальных методов были предложены следующие критерии, учитывающие специфику этого класса методов:

• Основа метода

• Необходимость взвешивания

• Наличие итераций

• Вычислительные особенности

• Ограничения на размеры обрабатываемого фрагмента

• Положительные стороны

• Отрицательные стороны

• Модифицируемые параметры метода

По этим параметрам было проведено сравнение четырёх глобальных методов:

1. Метод наименьших квадратов

2. Метод конечных элементов

3. Метод кусочно-линейного моделирования

4. Метод рекурсивного сглаживания

Помимо теоретической части, в рамках сравнения методов развёртки была проведена обработка экспериментальных данных. Для сравнения были выбраны три метода: модифицированный метод вычетов, метод частотного анализа и метод наименьших квадратов. На равнинной части все три метода работают одинаково хорошо: анализ разворота фазы по Баргузинской долине показал их согласованную работу: на квадратном участке изображения со стороной 512 элементов (около 10 км) среднее разности двух фазовых функций, развёрнутых двумя любыми методами, составляет 0.01 радиан при среднеквадратичном отклонении 0.05 радиан. В условиях пересечённого рельефа разница более заметна: в случае часто следующих интерференционных полос метод наименьших квадратов показал наихудшую работу, максимальная ошибка составила более 8 радиан. Вместе с тем, метод наименьших квадратов лишён главного недостатка локальных методов: в них маскированные области остаются «белыми пятнами», тогда как глобальные методы дают

решение задачи по всему полю изображения, без пропусков. Интересный результат показал спектральный анализ развёрнутой фазы (рис.1).

№---Р-®Я[Г-

*> г |

наим.кв. 1 наим.кв.« к

V. - '•)

вычет. вычет,

чаот.ан. чаег.ан. * ~

ю

Рис. 1. Спектральные характеристики интерферометрической фазы, развёрнутой тремя алгоритмами.

Следует отметить заметное подавление шумов на частотах более ^тах/Ю (где — частота следования отсчётов на интерферограмме) при предварительном сглаживании неразвернутой интерферограммы скользящим окном 5x5 и дальнейшем прямом интегрировании фазы с использованием маски недостоверных точек по методу частотного анализа. Ослабление шумов в данном случае достигает 6-И2 дБ по сравнению с другими методами (рис. 1, левый график). Уровень шумов для метода наименьших квадратов и метода вычетов на высоких частотах совпадает с точностью до долей децибела. Это значит, что метод наименьших квадратов не обладает свойствами низкочастотного фильтра. Для этого же метода на низких частотах видно превышение мощности спектральных компонент до 3 дБ (рис. 1, правый график). Это говорит о наличии крупномасштабных отклонений аппроксимации фазы, развернутой этим методом, от истинной.

Таким образом, показано, что устранение 2я-неоднозначности — один из наиболее сложных при реализации этапов интерферометрической обработки. Не существует универсального метода решения этой задачи: при выборе метода, оптимального для данных условий, следует учитывать особенности интерферограммы — характер распределения по изображению особых точек, их плотность, степень связности регулярных областей, — а также особенности постановки задачи: допустимость «белых пятен», допустимость локальных искажений с целью получения глобальной картины, и т.д.

Глава 4 посвящена обзору методов возникшей чуть более десяти лет назад отрасли интерферометрии — поляриметрической интерферометрии.

Глава носит справочный характер, представляя собой основу для следующего за ней материала глав 5 и 6: в ней введены понятия и величины, с которыми оперирует поляриметрическая интерферометрия, охарактеризованы используемые модели рассеяния, очерчен круг решаемых задач, проиллюстрированы основные направления исследований.

Глава 5 вводит в рассмотрение поляриметрические сигнатуры ин-терферометрической когерентности, сочетающие в себе информацию о рассеивающей поверхности/среде одновременно с интерферометрической и поляриметрической точки зрения. Сигнатуры когерентности имеют в качестве предшественника широко используемые в радиолокационной поляриметрии поляризационные сигнатуры. В современных публикациях поляризационная сигнатура — это графически представленная зависимость эффективной поверхности рассеяния (ЭПР) сигнала от параметров эллипса поляризации.

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

В интерферометрических работах принято вычислять когерентность по формуле

у= . {г(*г'-?1 :,где (1)

. г;.ч г;:;

21, 22 — значения коэффициента обратного рассеяния на первом и втором изображении;

Е() — операция взятия среднего по пространственной выборке;

* — операция комплексного сопряжения.

Графическое представление значений у для разных вариантов поляризации дают согласованную (для совпадающих поляризаций на излучении и приёме) или перекрёстную (для ортогональных поляризаций) сигнатуру интерферометрической когерентности.

Примеры сигнатур когерентности для согласованной и перекрестной комбинации поляризаций приведены на рис. 2. Вертикальная ось соответствует среднему значению (по некоторому пространственному ансамблю — в настоящей работе использовано скользящее окно усреднения размером 5x5 элементов разрешения) модуля интерферометрической когерентности, две горизонтальные оси — угол ориентации (ось направлена слева направо) и угол эллиптичности (ось трансверсальна плоскости страницы).

В данной работе основное внимание уделяется согласованным сигнатурам, которые соответствуют ситуации, в которой поляризации приёмника и источника совпадают. Везде далее слово «сигнатура» означает именно согласованную сигнатуру когерентности.

1т1 М

0.351 0.35

....... ■

¿у&Яфщ&АВМ+м ..■ву*

| ЙШШ ' ШЩШ

Щ щ ^ > ж

ЩШР Щу?

шШЩШЖШШ

0.25

О __ Л1Ж у о 95_ 0 _............_ 2п ¥

;п/4 - ""к/4

л/4 К п/4Л

а) б)

Рис. 2. Согласованная (а) и перекрёстная (б) сигнатуры когерентности для открытой поверхности (поле) в С-диапазоне.

Для иллюстрации свойств поляризационных сигнатур когерентности и классификации земных покровов, основанной на них, в настоящей работе использованы данные полнополяриметрического РСА 8111-С по территории юго-восточного Прибайкалья. Два повторяющихся прохода 7.10.1994 и 9.10.1994 позволяют вести интерферометрическую обработку, а два частотных диапазона (С-диапазон с длиной волны 5.6 см и Ь-диапазон с длиной волны 23 см) — видеть зависимость сигнатур от несущей частоты сигнала локатора. В пределах сцены расположены часть акватории Байкала, дельта реки Селенги и прилегающие районы. Характер подстилающих покровов достаточно разнообразен: имеется прибрежная полоса, культивируемые территории (поля), леса двух типов (смешанные и хвойные), необрабатываемые открытые пространства (луга, болота, пойменные земли).

Как правило, общий вид сигнатуры когерентности таков, что профили графиков при круговой поляризации (% = +Я" / 4) проходят ниже профиля линейной поляризации — 0 ), что хорошо видно на рис. 2. Глобальный максимум сигнатуры обыкновенно располагается также приблизительно на линии % = 0 , чаще всего на вертикальной ((// = тг/2) или горизонтальной (I// = 0 ну/ = 7Г) поляризации. Что касается глобального минимума, то он может быть как на профиле графика, соответствующем линейной поляриза-

ции {X = 0), так и на краю сигнатуры (на профилях, соответствующих круговым поляризациям % = ±7г/ 4). Кроме случая точечных максимумов и минимумов, возможен вариант максимальной (или минимальной) линии, как на рис. 2а. Выделяется также случаи хаотического распределения максимумов и минимумов; это бывает при сниженном среднем значении модуля когерентности по сигнатуре. Нулевая или околонулсвая интерферомстричсская когерентность свидетельствует о полной или почти полной декорреляции двух изображений, поэтому и сигнатура перестаёт иметь сколько-нибудь правильную форму.

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

Безлесные территории

Пример сигнатур на рис. 2 представляет именно этот тип поверхности. Среднее значение модуля когерентности по различным безлесным участкам для С-диапазоне составляет около 0.3, а для L-диапазона около 0.8.

Как было упомянуто выше, декорреляция изображении происходит по причинам различной природы. Мы не принимаем в рассмотрение так называемую «пространственную декорреляцию», которая возникает из-за различия в углах обзора для двух траекторий носителя, во-первых, потому что в данном случае эта составляющая несущественна, а во-вторых, потому что она медленно меняется от одного края сцены к другому, и для выбранных размеров тестовых участков не сказывается. Относительно декорреляции из-за соотношения сигнал-шум можно заметить следующее. График имеет максимальную прямую при нулевой эллиптичности. При этом традиционная поляризационная сигнатура, построенная по тому же участку изображения, имеет чётко выраженный минимум посередине сигнатуры ( % = 0,1// = тг/2 ), что соответствует вертикальной поляризации. В естественном предположении одинаковой шумовой составляющей для всех углов ориентации получаем, что соотношение сигнал-шум для нулевой эллиптичности зависит от угла ориентации, а модуль когерентности - - нет. Таким образом, декорреляция из-за соотношения сигнал-шум не вносит заметного вклада в значение модуля когерентности на линии ^ = 0 . Напротив, падение величины когерентности от линии нулевой эллиптичности к краям, то есть к линиям круговых по-

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

Водные поверхности

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

Леса

Леса являются одними из наиболее интересных объектов изучения поляриметрической интерферометрии благодаря наличию объёмного отражающего слоя. Тестовым в данной работе послужил лес у села Истомино. Это смешанный лес с неравномерным распределением пород деревьев. Из хвойных преобладают сосны, лиственницы и (более редкие) ели, из листвен-

ных — березы и осины. Богатый подлесок формирует нижний объемный отражающий слой.

В отличие от предыдущего типа поверхности, сигнатуры когерентности леса для Ь- и С-диапазонов не всегда имеют совпадающую форму для одного и того же участка леса. Близкие формы сигнатур показаны на рис. 3. Максимум соответствует линейной горизонтальной поляризации. Разница в среднем значении и вариации: и то, и другое в С-диапазоне снижено по сравнению с Ь-диапазоном. К этому типу можно отнести участки леса со множеством крепких горизонтальных ветвей, на которые не действуют ветровые колебания.

М м

0.23

0.75

0.17--п/4

2711,7

п/ЛХ

а) б)

Рис. 3. Сигнатуры когерентности леса в Ь-диапазоне (а) и в С-диапазоне (б)

Однако в С-диапазоне возникают и совсем другие сигнатуры: одни имеют максимум на вертикальной поляризации (рис. 4а), другие имеют максимумы на обеих согласованных поляризациях (рис. 46). Картину когерентности, показанную на рис. За), могут дать участки молодого леса с параллельными стволами и тонкими ветками. Тип сигнатуры, показанный на рис. 36), скорее всего, соответствует густой чаще, изобилующей стабильными отражателями и вертикальной, и горизонтальной ориентации. В соответствии с распределением максимумов сигнатуры можно выделить три типа леса. Каждый из них характеризуется определённой ориентацией наиболее стабильных отражателей в С-диапазонс. Среднее значение когерентности для всех трёх типов леса примерно одинаково и составляет около 0.2.

Кроме названных типов сигнатур когерентности, попадаются и другие - смешанные. Например, сигнатура, напоминающая поле, но с плавным максимумом (или двумя) на линии % = 0 . Чаще всего это означает, что в

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

0.23

-л/4

п/4^

о.:

0.19 -ti/4

ШШШ

к/К

а) б)

Рис. 4. Сигнатура когерентности леса второго (а) и третьего (б) типов в С-

диапазоне

Сравнение с картой лесных пород не дало основания считать, что три выделенные типа леса соответствуют определённым видам деревьев. Например, к типу 2 относится как участок сосен в центре леса, так и опушка, состоящая из берёз и осин. Общее свойство этих двух участков — частые стволы, тонкие горизонтальные ветви, неразвитый подлесок. К типу 1 относятся, как правило, участки смешанного леса, состоящие из деревьев разного возраста. Следует отметить также, что съёмка производилась в октябре, когда большинство лиственных деревьев теряет листву. Возможно, в летний период распределение сигнатур по площади леса имеет другой характер.

Шестая глава также относится одновременно к поляриметрии и интерферометрии, поскольку предлагает способ поэтапной регистрации полной матрицы рассеяния в схеме интерферометрической съёмки. Требования к орбитам носителя при этом остаются теми же, что и для обычной интерферометрической съёмки, оценка степени декорреляции проводится по формуле (1). Полученная квази-матрица рассеяния записывается в виде

Sq =

S S

1

нн 2

УН

S S

НУ 2

УУ

где нижние индексы указывают на поляризацию передачи/приёма, а верхние номер сеанса съемки.

Проведено моделирование возникающих искажений матрицы рассеяния в мультипликативном виде для четырёх базисных матриц. Отмечены искажения, возникающие вследствие амплитудной и фазовой составляющей искажающего множителя как по отдельности, так и совместно. На основании построенных искажённых базисных матриц предложен метод коррекции реальных квази-матриц, построенных по частичным поляризационным измерениям. Использован комплексный корректирующий множитель Б)^ / Буц для второй строки матрицы. Для оценки точности проведённой коррекции результаты сравнивались с исходной полной матрицей рассеяния, полученной в первом сеансе (рис. 5, 6).

• • • Коррекция Й-: ■ 5:г

Л" - А* Аг

▼А. о ч 1- 'чч""' ч ч у с

Рис. 5. Квази-матрица рассеяния, результат се коррекции и исходная полная матрица рассеяния, Ь-диапазон, тип поверхности: поле.

Зч -» Коррекция-^ -» *

Лг 1 л- \с Аг

шшШк . г

\ ИШ Шйшш- ЩрЦНр I

5- ■щш ШШ» _ о' 1°г

\ \Г

\ г 4 л

Рис. 6. Квази-матрица рассеяния, результат сё коррекции и исходная полная матрица рассеяния, С-диапазон, тип поверхности: поле.

В результате обработки экспериментальных данных поляриметрического радиолокатора БЖ-С продемонстрировано, что для однородных участков леса и поля (открытого пространства) — двух наиболее распространённых типов поверхности — квази-матрица рассеяния в Ь-диапазонс поддаётся

коррекции почти полностью (отличия в поляризационных сигнатурах составляют не более 5% полной вариации), в С-диапазоне в обоих случаях результат коррекции заметно хуже (отличия достигают 30%). Показано, что остаточный мешающий эффект связан с временной декорреляцией, которая в С-диапазоне проявляется сильнее при одном и том промежутке времени между съёмками.

В Заключении сформулированы результаты проделанной работы:

1. Систематизировано состояние радиолокационной интерферометрии к настоящему моменту:

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

б) проанализированы этапы интерферометрической обработки радиолокационных данных, а также приведены соотношения, указывающие пределы точности измеряемых в интерферометрии величин — радиального смещение и относительных высот рельефа (глава 2);

2. Детально рассмотрены методы развёртки фазы, для одного из которых предложена модификация, ускоряющая работу алгоритма в условиях областей декорреляции и сложного рельефа на снимаемой местности; проведено сравнение в классах локальных и глобальных методов по основным характеристикам методов, а также на примере обработки экспериментальных данных (глава 3);

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

а) приведены базовые соотношения величин, изложены основные существующие методы (глава 4);

б) введены в рассмотрение сигнатуры когерентности — предложенная автором визуализация поляриметрико-интерферометрических данных, могущая служить основой для классификации (глава 5);

в) предложен и обоснован способ регистрации полной матрицы рассеяния по частичным поляризационным измерениям в схеме интерферометрической съёмки с повторяющихся орбит носителя (глава 6).

Основные результаты работы заключаются в следующем:

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

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

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

проанализирована последовательность этапов иптерферо-метрической обработки и показана их взаимосвязь;

систематизированы основные направления поляриметрической интерферометрии:

приведены классы и конкретные примеры прикладных задач, решаемых с помощью поляриметрической интерферометрии.

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

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

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

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

ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ - статьи it журналах из рекомендованного перечня ВАК:

1. Захаров, А.И., Захарова, JI.II. Значимость информации о фазе отражённого сигнала при радиолокационном картировании земных покровов // Радиотехника. — 2003. — №12. — С.70-73.

2. Захарова, J1.I1., Захаров, А.И. Сравнение некоторых современных методов разворота разности фаз в радиолокационной интерферометрии // Радиотехника и электроника. — 2003. — Т. 48, №10. — С. 1208—1213.

3. Захарова, Л.П. Сигнатуры когерентности в поляриметрической интерферометрии для классификации земных покровов // Радиотехника. — 2005. — №8. — С. 45—50.

4. Захарова, Л.Н. Применение методов радиолокационной интерферометрии для исследования динамики берегов р. Волги в районе г. Ульяновска // Нелинейный мир. — 2007. —Т. 5, №5. — С. 273-274.

5. Захарова, Л.Н. Сравнение методов сигнатур когерентности и области когерентности для классификации земных покровов // Нелинейный мир. — 2008.

— Т. 6, №4. — С.245246.

6. Захарова, Л.Н.. Формирование полной матрицы рассеяния зондируемых сред по частичным поляризационным измерениям РСА // Нелинейный мир. — 2009. — Т.7, №3 — С.180—181.

7. Арманд, Н.А., Захаров, А.И., Захарова, Л. II. Космические радары с синтезированной апертурой в дистанционном зондировании Земли — современные системы и перспективные проекты // Исследование Земли из космоса. — 2010 —№2, —С. 3-13.

- труды и тезисы докладов, статьи и прочих журналах:

8. Захарова, Л.Н., Захаров, А.И. Проблема разворота разности фаз в радарной интерферометрии // Тр. XX и XXI Всероссийских Симп. «Радиолокационное исследование природных сред» (XXI Симпозиум, Санкт-Петербург, 15-17 апреля 2003 года). — СПб., 2003. — Вып. 3. — С. 35-47.

9. Zakharova, L. Comparison of Global and Local Approach to Phase Unwrapping For a Rugged Terrain // Proc. of Workshop on ERS SAR Interferometry (Fringe-2003), Frascati, Italy, December I 5. — 2003.

10. Захарова, Л.Н. Дистанционное зондирование земных покровов с помощью методов поляриметрической интерферометрии // Тр. XXII Всеросс. научн. конф. «Распространение радиоволн» (РРВ-22, RP XXII), Ростов-на-Дону, п. Лоо, 22-26 сентября 2008 г. — 2008. — С.32—35.

11. М. Raimadoya, et al. ALOS-Indonesia Pol-InSAR Experiment (AIPEX): Final Result / M. Raimadoya, B. Trisassongko, L. Zakharova, N. Sarbini // Abstract Book of the 4th Joint PI Symposium of ALOS Data Nodes for ALOS Science Program, Tokyo, November 15-17, 2010. — 2010. —P. 59.

12. Zakharova, L., Darizhapov, D., Kirbizhekova, I. Study of natural objects in Transbaikalia by means of polarimetry and polarimetric interferometry // Proc. of the 5lh European Conference on Synthetic Aperture Radar (EUSAR-2004), Ulm, Germany, May 25-27, 2004. — 2004. — Vol. 2. — P.733—736.

13. Zakharova, L. Polarimetric Coherence Signatures for Land Cover Study // Proc. of the 41'1 International Symposium on Retrieval of Bio- and Geophysical Parameters from SAR data for Land Applications, Innsbruk, Austria, November 16-19 2004.

— 2004. — Электронное издание на CD.

14. Захарова, JI.II. Классификация земных покровов по данным РСЛ с использованием методов поляриметрической интерферометрии // Тр. XXV""1 вес-росе. сими. «Радиолокационное исследование природных сред», Санкт-Петербург, 17-19 апреля 2007. — СПб, 2009. —Выпуск 7. —С. 83—90.

15. Захарова, JI.H. Новый метод анализа данных полноноляримстричсского РСД при использовании интерферометрии с повторяющихся орбит носителя // Сб. докл. XIII Международной науч.-техн. Коиф. «Радиолокация, навигация, связь» (RLNC-2007), Воронеж, 17-19 апреля 2007 года. — 2007. — С. 1761 — 1769.

16. Захарова, JI.II. Применение методов поляриметрической интерферометрии для исследования вертикальной структуры леса // Тр. XXII Сими, по радиолокационному зондированию природных сред, Санкт-Петербург, 18-20 апреля 2004. — 2004.

17. Zakliarova, L. Polarimelric Coherence for Land Covers Classification // Proc. of the 28lh General Assembly of URSI (Union Radio-Scientifique Internationale), New Delhi, India, October 22-29, 2005. — 2005.

18. Захарова JI.II. Применение методов поляриметрической интерферометрии для исследования растительных покровов // Сб. докл. Второй вссросс. научи, коиф. «Дистанционное зондирование земных покровов и атмосферы аэрокосмическими средствами», Санкт-Петербург, 16-18 июня 2004 года. — СПб., 2004. — Т. I. — С. 65- 69.

19. Захарова, Л.Н. Использование иптсрфсрометрической когерентности но данным полиополяриметричсского РСД для классификации земных покровов // Сборник трудов Пятой Юбилейной Открытой Всероссийской Конференции «Дистанционное зондирование Земли из Космоса», Москва, ИКИ РЛП, 1216 ноября 2007 года. — 2007.

20. Zakharova, L. Land cover classification based on polarimetric coherence signatures // Abstract book of the 3rd International Workshop on Science and Applications of SAR Polarimetry and Polarimetric Interfcrometry (POL1NSAR-2007), Frascati, Italy, January 22-26, 2007. - 2007. — P. 47.

21. Захарова, Л.Н. Возможности применения радиолокационной поляриметрической интерферометрии для исследования земных покровов// Тезисы докладов V конференции молодых учёных «Фундаментальные и прикладные космические исследования», Москва, ИКИ, 8-9 апреля 2008 года. — 2007. — С. 19.

22. Zakharova, L. Forest classification by means of two POLINSAR techniques // Proc. of the 5th Intl. Symp. on Retrieval of Bio- and Geophysical Parameters from SAR data for Land Application (BioGcoSAR-2007), Bari, Italy, September 2528, 2007. — 2007.

23. Zakharova, L. Polarimetric Coherence Signatures and Coherence Regions in Complex Unit Circle // Proc. of the 7lh European Conf. on Synthetic Aperture Radar (EUSAR-2008), Fricdrichshafcn, Germany, June 2-6, 2008. — 2008. — Vol. 3, —P. 309—312.

24. Zakharova, L., Zakharov, A. Quasi-scattering matrix registration in repeat pass mode // Proc. of the 4Ul Intl. Workshop on Science and Applications of SAR Polarimetry and Polarimetric Interfcrometry (POLINSAR-2009), Frascati, Italy, 26-30 January 2009. —2009.

25. Захарова, Л.Н., Захаров, A.II. Исследование подвижек почвы в районе г. Ульяновска методом радиолокационной интерферометрии // Тр. XXIV

Сими, по радиолокационному зондированию природных сред, С.-Петербург, 18-20 апреля 2006 г. — 2006. — С. 132-138.

26. Zakliarova, L. Interferometric study of Landslides Activity of Volga River Banks // Proc. of the 61'1 European Conference on Synthetic Aperture Radar (EUSAR-2006), Dresden, Germany, May 16-18, 2006. — 2006.

27. Захарова, JI.H. Наблюдение оползневых процессов методами радиолокационной интерферометрии па примере Ульяновска // Тр. XXII Всеросе. научн. конф. «Распространение радиоволн» (РРВ-22) RP XXII, Ростов-на-Дону, п. Лоо, 22-26 сентября 2008 г. — Ростов-на-Дону, 2008. — Т.З — С.36-39.

28. Захарова, JI.II., Захаров, А.И. Многолетние наблюдения динамики земных покровов в Ульяновске по данным космических радиолокаторов // Тр. IV Всеросс. научн. шк. и конф. "Радиофизические методы в дистанционном зондировании сред" (Муром'2009), Муром, 30 июня - 3 июля 2009 г. — 2009.

29. Лебедева, М.А., Захарова, Л.П., Захаров, А.И. Применение метода РСА-интерферометрии для изучения медленных смещений (на примере активных оползневых деформаций вблизи Северомуйского тоннеля БАМ)// Первая молодежи, тектонофизич. шк.-семпнар, Москва, ИФЗ им. О.Ю.Шмидта РАН, 21-25 сентября 2009 г. — 2009.

30. Захарова, Л.Н., и др. Применимость радиолокационной интерферометрии к мониторингу деформаций земной поверхности в районе Северомуйского участка БАМ // Тр. конф. «Всероссийские радиофизические научные чтения-конференции памяти H.A. Арманда. Всероссийская научно-практическая конференция «Космическая радиолокация»», Муром, 28 июня - 1 июля 2010 г. — 2010. — Электронное издание (CD).

31. Захаров, А.И., Захарова, Л.II., Лебедева, М.А. Применение РСА-интерферометрии для мониторинга транспортной инфраструктуры в зонах с опасной динамикой покровов// Журнал радиоэлектроники (электронный журнал) — 2010. — № 10 (приложение).— С. 32—38. — htlp://ire.cplire.ru/ire/library/Ulan-Udü-20l 0/pdffiles/cl 4.pdf

32. Лебедева, М.А., и др. Эндог енные и экзогенные деформации в зонах активных разломов Верхнеапгарско-Муиской междувпадинпой перемычки по данным дифференциальной РСА-иптерферометрии / Лебедева, М.А., Сань-ков, В.А., Захаров, А.И., Захарова, Л.II. // Журнал радиоэлектроники (электронный журнал). — 2010. -..... №10 (приложение). — С. 65—71. —

htlp://ire.cplire.ru/ire/library/Ulan-Ude-20l 0/pdfnies/cl 7.pdf

33. Захаров, А.И., и др. Методы космической радиолокационной интерферометрии в гляциологии / Захаров, А.И., Захарова, Л.II., Синило, В.П, Петрова, Л.Г. // Тр. XXV010 Всеросс. Сими. «Радиолокационное исследование природных сред», Санкт-Петербург, 17-19 апреля 2007 года. — СПб, 2009. — Вып. 7, —С. 91-94.

34. Zakharova, L., et al. On the Use of ERS INSAR Data in the Ecological Monitoring of the Baikal Region // Proc. of the Intl. Geoscience and Remote Sensing Symp. (IGARSS'03), Toulouse, France, 21-25 July 2003. — 2003. — Vol. 4. — P. 2477—2479.

35. Захаров, А.И., Захарова, Л.П. Применение интерферометрии для мониторинга районов добычи и транспортировки нефти и газа / Загл. и текст парал. рус., англ. // ROGTEC: Российские нефтегазовые технологии. — 2006. — Вып. 5,— С. 58-67.

Отпечатано в типографии ООО «Мещёра», г. Щёлково, ул. Свирская, д.8а. Тираж 100 экз. Заказ №435

СПИСОК ИСПОЛЬЗОВАННЫХ СОКРАЩЕНИЙ.

ВВЕДЕНИЕ.

1. ОРГАНИЗАЦИЯ ИНТЕРФЕРОМЕТРИЧЕСКОЙ СЪЁМКИ.

1.1. Геометрия интерферометрической съёмки.

1.1.1.' Интерферометрическая фаза.

1.1.2. Измерение высот рельефа.

1.1.3. Измерение смещений поверхности.

1.1.4. Критические величины орбитальных параметров.

1.2. Орбитальные конфигурации для реализации интерферометрической съёмки.

1.2.1. Полноактивные системы.

1.2.2. Полуактивные системы.

2. ЭТАПЫ ИНТЕРФЕРОМЕТРИЧЕСКОЙ ОБРАБОТКИ, ИЗМЕРЯЕМЫЕ ВЕЛИЧИНЫ И ИХ ТОЧНОСТИ.

2.1. Этапы интерферометрической обработки.

2.1.1. Совмещение.

2.1.2. Вычисление когерентности.

2.1.3. Вычисление разности фаз.

2.1.4. Фильтрация.

2.1.5. Интерпретация когерентности.

2.1.6. Устранение набега фазы по дальности.

2.1.7. Устранение 2л;-неоднозначности разности» фаз.

2.1.8. Устранение перепадов фазы, относящихся к рельефу.

2.1.9. Пересчет разности фаз в колебания высоты рельефа.

2.1.10. Пересчёт разности фаз в величину подвижек.

2.1.11. Визуализация.

2.1.12. Геокодирование (картографическая привязка).

2.2. Измеряемые величины и их точности.

2.2.1. Оценка степени декорреляции.

2.2.2. Оценка точности измерения высоты.

2.2.3. Точности интерферометрических измерений, достигаемые существующими радиолокационными системами.

2.2.4. Оценка точности локальных деформаций.

3. АНАЛИЗ МЕТОДОВ РАСКРЫТИЯ ФАЗОВОЙ НЕОДНОЗНАЧНОСТИ ПРИ ИНТЕРФЕРОМЕТРИЧЕСКИХ ИЗМЕРЕНИЯХ (МЕТОДЫ РАЗВОРОТА

ФАЗЫ).

3.1. Постановка задачи.

3.2. Локальные методы разворота фазы.

3.2.1. Метод вычетов и его модификация.

3.2.2. Потоковый метод разворота фазы.

3.2.3. Частотный метод анализа фазовой неоднозначности г интерферограммы.

3.2.4. Метод восстановленных значений.

3.2.5. Метод, основанный на локальной фильтрации интерферограммы.

3.2.6. Метод определения граничных линий.

3.3. Глобальные методы раскрытия фазовой неоднозначности.

3.3.1. Метод наименьших квадратов.

3.3.2. Метод рекурсивного сглаживания.

3.3.3. Метод кусочно-линейного моделирования.

3.3.4. Взвешенное развёртывание фазы методом конечных элементов.

3.4. Сравнение методов.

3.4.1. Сравнение локальных методов между собой.

3.4.2. Сравнение глобальных методов между собой.

3.4.3. Сравнение результатов работы локальных и глобальных методов на примере данных ERS.

4. ПОЛЯРИМЕТРИЧЕСКАЯ ИНТЕРФЕРОМЕТРИЯ.

4.1. Поляризационный базис, вектор рассеяния, матрица когерентности.

4.2. Векторная когерентность.

4.3. Преобразование когерентности при замене поляризационного базиса

4.4. Связь между поляриметрической и интерферометрической разностью фаз в поляриметрической интерферометрии.

4.4.1. Интерферограммы с разными поляризациями.

4.4.2. Интерферограмма для поляриметрических разностей фаз.

1 4.4.3. Дифференциальные интерферограммы вида XYUV-HHHH.

4.5. Модели рассеяния в скалярной и векторной (поляриметрической) интерферометрии.

4.5.1. Объёмное и поверхностное рассеяние в векторной интерферометрии

4.5.2. Идеи двух- и трёхуровневых моделей рассеяния для растительности

4.5.3. Линейная модель когерентности.

4.6. Поляриметрическая интерферометрия с двумя базами.

4.7. Современные области применения интерферометрии, решаемые задачи, измеряемые параметры.

4.7.1. Классификация.

4.7.2. Параметры высотной структуры растительности.

4.7.3. Топографические высоты и уклоны.

4.7.4. Мониторинг естественных процессов и чрезвычайных ситуаций

4.7.5. Обнаружение искусственных объектов.

5. ПОЛЯРИМЕТРИЧЕСКИЕ СИГНАТУРЫ ИНТЕРФЕРОМЕТРИЧЕСКОЙ КОГЕРЕНТНОСТИ.

5.1. Преобразование поляризационного базиса.

5.2. Сигнатуры когерентности: определение.

5.3. Общие свойства сигнатур.

5.4. Классификация земных покровов, основанная на сигнатурах когерентности.

5.5. Классификация леса, основанная на сигнатурах когерентности.

5.6. Области когерентности и сравнение классификации.

5.6.1. Замена поляризационного базиса.

5.6.2. Классификация на основе областей когерентности.

6. РЕГИСТРАЦИЯ МАТРИЦЫ РАССЕЯНИЯ ПО ЧАСТИЧНЫМ ПОЛЯРИЗАЦИОННЫМ ИЗМЕРЕНИЯМ В РЕЖИМЕ

ИНТЕРФЕРОМЕТРИЧЕСКОЙ СЪЁМКИ.

6.1. Моделирование квази-поляриметрической матрицы рассеяния.

6.2. Применение метода к данным радиолокатора SIR-C.

Актуальность исследований

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

Исследования Земли с помощью радиолокатора с синтезированной апертурой, расположенного на борту самолета или спутника — это современный высокопроизводительный метод дистанционного зондирования, позволяющий получать двумерные распределения интенсивности отраженного сигнала, а также изучать отражательные свойства и структуры рельефа поверхности. Преимуществами радиолокаторов перед оптическими- сенсорами является возможность съёмки независимо от солнечного освещения и от наличия/отсутствия облачности. Прием отраженного сигнала на две антенны на борту носителя (или же, в схеме с одной антенной, съемка выбранного участка поверхности с близких траекторий носителя) и дальнейшая совместная обработка позволяют при помощи анализа фазовош составляющей эхо-сигнала добавить третье измерение к двумерным радиолокационным изображениям, выявляя, в зависимости от поставленной задачи и начальных условий, относительные высоты рельефа, степень декорреляции или изменение положения отражающей поверхности.

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

• извлечение информации о рельефе поверхности Земли;

• детектирование площадных подвижек поверхности;

• обнаружение изменений поверхности за время между съёмками;

• оценку скорости движущихся целей.

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

Состояние исследований (обзор литературы)

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

Одной из первых публикаций, показывающей возможность количественной дистанционной оценки перепадов рельефа на поверхности Земли, является [1]. В этой работе 1974 года приведены первые результаты применения ин-терферометрического метода к радиолокационным данным, полученным с борта самолёта, для построения карты рельефа. Аппаратная реализация метода такова, что в качестве результата были получены линии постоянной высоты в координатах «азимут—дальность». Следующие публикации на эту тему относятся к середине 1980-х годов и также касаются измерений, полученных с самолёта [2]. В работе использованы данные радиолокатора Ь-диапазона (съёмка с самолёта, однопроходная, 11м жёсткой базы с антеннами на крыльях, точность измерения высоты-в пределах 10 м). Почти одновременно были продемонстрированы возможности отдельной методики интерферометрических измерений: в случае ориентации базы интерферометра вдоль трассы полёта становится возможным измерять скорости движущихся на поверхности объектов, например, океанских течений [3]. В настоящей работе эта тема остаётся в стороне, равно как и интерферометрическая съёмка с самолётных носителей, однако три названных публикации имеют безусловную ценность как первые исследования в области радиолокационных интерферометрических измерений.

Первая публикация, касающаяся интерферометрической обработки данных, полученных с космической орбиты, имела целью построения карты высот с повторяющихся орбит по данным радиолокатора SIR-B [4]. Также в этой работе уделяется внимание трудностям, возникающим при обработке данных с не параллельных орбит (пересекающихся под малым, но всё же заметным углом). Другой радиолокатор космического базирования, установленный на борту космического аппарата SeaSat, был использован для иллюстрации идеи дифференциальной интерферометрии в многопроходной схеме съёмки с компенсацией составляющей фазы, отражающей информацию о рельефе, за счёт дополнительной интерферограммы с другой пространственной базой [5, 6].

В 1991 году был запущен искусственный спутник Земли ERS-1 с радиолокатором на борту, ознаменовавший начало новой эры в радиолокационной интерферометрии, поскольку 35-дневный период повторения съёмок и охват съёмками практически всей поверхности планеты благодаря околополярной орбите позволили получать пары изображений, пригодные для интерферометрической обработки, в небывалом ранее количестве. Таким образом, начиная с 1992 года, количество публикаций, посвящённых интерферометрической тематике в радиолокации, стремительно увеличивается. На симпозиуме, посвящён-ным первым результатам работы ERS-1, а также на крупнейшем международном ежегодном симпозиуме по дистанционному зондированию IGARSS (International Geoscience And Remote Sensing Symposium), проведённым в 1992 году, появились доклады, посвящённые радиолокационной интерферометрии для картирования рельефа Земли. Также вышли работы, посвящённые теоретичеt ским основам и анализу ошибок интерферометрических измерений [7,8]. i f а i i

I1

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

1994 год ознаменовался съёмками радиолокаторов SIR-C/X-SAR (совместный американо-германо-итальянский проект, аппаратура производила съёмку в трёх частотных диапазонах X, С и L, причём в двух последних в поляриметрическом режиме). Несмотря на то, что съёмка производилась в течение всего двух коротких периодов 9.04.1994—20.04.1994 и 30.09.1994-11.10.1994, полученные данные представляют большую научную ценность, в том числе для развития интерферометрической и поляриметрической радиолокации. Помимо ежегодных симпозиумов, посвящённых дистанционному зондированию вообще и при использовании данных космического аппарата ERS-1 в частности — IGARSS, 2nd ERS-1 Symposium, Progress in Electromagnetics Research Symposium (PIERS), 1st Workshop on ERS-1 Pilot Projects — в 1994 году проводится первый международный научный симпозиум, посвящённый задачам собственно радиолокационной интерферометрии — 1st Workshop on SAR Interferometry. К этому времени тематику как обсуждений на конференциях, так и журнальных публикаций можно условно разбить на три группы: интерферометрия как средство измерения рельефа [13-15], интерферометрия как средство измерения динамики поверхности [16-18], проблемы интерферометрической обработки радиолокационных данных и направления их решения [19-23]. Среди последних следует особо выделить растущее количество работ, посвящённых двумерной развёртке фазы (устранения 27г-неоднозначности), которая является нетривиальной задачей и представляет собой одну из самых вычислительно ёмких процедур интерферометрической обработки.

В 1995 году были запущены сразу два ИСЗ с РСА на борту: европейский ERS-2 (с радиолокатором, имеющим те же характеристики, что и у ERS-1 и орбитой, организованной так, что съёмка каждого участка поверхности Земли проводится спустя ровно сутки после съёмки ERS-1) [24] и канадский RADARS AT с радиолокатором С-диапазона [25]. Расширяется сфера приложений радиолокационной интерферометрии: появляются работы, посвящённые исследованиям состояния лесов [26] и других типов растительности [27], вулканов [28] и береговых полос [29]. Усиливается внимание к точностям интерфе-рометрических измерений [30] и причинам искажений на интерферограммах (в том числе, влиянию атмосферы [31, 32]), рассматриваются различные аспекты интерферометрической обработки, в частности, по-прежнему большое внимание уделяется поиску оптимального метода развёртки фазы [33].

В 1996-1997 годах возрастает количество публикаций, посвящённых аспектам интерферометрической обработки для конкретных инструментов (SIR-C/X-SAR, ERS-1/2 Tandem, JERS-1, SIR-C и планируемая в тот момент миссия SRTM), продолжаются исследования возможностей интерферометрии для регистрации рельефа и подвижек земной поверхности, вместе с тем рассматриваются различные особенности обработки и коррекции интерферометрических изображений. Кроме того, появляются обзоры, иллюстрирующие потенциальные возможности интерферометрического направления с точки зрения научных и хозяйственных приложений. Материалы нескольких секций трёх крупных международных конференций (IGARSS, EuSAR, Fringe) посвящены радиолокационной интерферометрии. Всё большее внимание уделяется приложениям к изучению растительных покровов ([34]) и анализу точности измеряемых величин, особенно атмосферным искажениям [35].

В 1998 году возникает новая перспективная ветвь интерферометрических исследований: поляриметрическая интерферометрия, использующая поляриметрические данные PC А [36]. Особенности методов поляриметрической интерферометрии позволяют оценивать характеристики объёмных рассеивающих сред (главным образом, растительных), на основе интерферометрической обработки сигналов с различными комбинациями поляризаций на излучении/приёме.

В 1999 году впервые рассмотрена возможность однопроходной интерфе-рометрической съёмки, организованной при помощи полуактивной конфигурации спутников: активный излучающий инструмент установлен на одном из них, остальные (малогабаритные) снабжены только приёмными устройствами, регистрирующими сигнал, причём орбиты всех космических аппаратов расположены достаточно близко, чтобы можно было формировать интерферометрическое изображение [37].

В 2000-м году детальное внимание было уделено различным этапам ин-терферометрической обработки: геокодированию [38], высокоточному совмещению пары изображений [39], учитыванию тропосферных [40] и ионосферных [41] эффектов. Кроме того, впервые рассмотрена возможность интерферомет-рической обработки данных двух разных инструментов, несколько отличающихся по частотному диапазону [42]. В 2001 году продолжают появляться публикации по всем названным выше направлениям, ширится поток исследований, связанных с техникой постоянных отражателей (permanent scatterers), которая позволяет оценивать смещения поверхности по данным даже с высокой степенью декорреляции [43]'. В 2002 году впервые упомянуты уголковые отражатели в связи с интерферометрическими измерениями [44]'.

В 2003 году в научной печати появляются первые результаты интерферо-метрической обработки данных ENVISAT [45], который был выведен на орбиту в марте предыдущего, 2002 года. Всё шире рассматриваются возможности совместной обработки данных, полученных из различных источников, для достижения наилучших результатов [46, 47] Также рассматриваются результаты интерферометрии для L-диапазона (радиолокатор JERS-1), на основе которых происходит планирование интерферометрической составляющей миссии Тегга-SAR-L [48]. В 2004 году появляется первое упоминание о планируемой миссии COSMO-SkyMed [49] — первой реализованной системы из четырёх спутников, следующих по орбите друг за другом и несущих идентичные инструменты. Множество работ посвящены обнаружению деформаций земной поверхности при помощи интерферометрии, а также методологии и частным аспектам обработки. Всё чаще привлекаются сторонние источники информации для совместного анализа с интерферометрическими измерениями. Впервые упоминаются планируемые миссии 11ас1агза1;-2 и ЯаёагБа^З [50].

С середины первого десятилетия 21 века были выведены на околоземную орбиту несколько новых космических аппаратов с РСА на борту. Это АЬОБ (Япония, 2006) с радиолокатором РАТРАК. Ь-диапазона, семейство СОЯМО-БкуМес! (Италия, 2007—2010) из четырёх спутников, оснащённых РСА Х-диапазона, КАЕ)АЛ8АТ-2 (Канада, 2007) с РСА С-диапазона, а также группировка Тегга8АИ-Х и ТапЕ)ЕМ-Х (Германия, 2007 и 2010) снова Х-диапазона. Преобладающие темы публикаций в эти годы: деформации земной поверхности, получение точных цифровых карт земной поверхности, особенности различных этапов обработки, особенности миссий, способы выявления и устранения атмосферного влияния, поляриметрическая интерферометрия. Продолжают совершенствоваться системы обработки, отслеживаются различные динамические явления на земной поверхности, предлагаются методы коррекции атмосферной составляющей, предлагаются модели для оценки возможностей распознавания деформаций методом интерферометрии. Также большое внимание уделяется точности полученных интерферометрическими методами измерений, которая подтверждается альтернативными способами [51, 52], а таюке продолжаются работы по детектированию деформаций, в том числе зданий [53—55] Сравниваются возможности интерферометрической обработки данных различных диапазонов [56]. Выходит несколько книг, посвящённых особенностям обработки радиолокационных данных, в каждую из которых входят разделы, по-свящённые интерферометрии [57-60]. Интенсивно развивается поляриметрическая интерферометрия [61-63].

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

Цели и задачи работы

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

1. Проведён анализ геометрических характеристик интерферометрической съёмки, рассмотрены различные варианты организации интерферометрической съёмки, в том числе, конфигурации космических аппаратов.

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

3. Показаны ¡измеряемые величины и приведены точности их оценки.

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

5. Изложены основы поляриметрической интерферометрии и проанализированы существующие на текущий момент её приложения.

6. Введён в рассмотрение новый метод классификации естественных покровов по поляриметрико-интерферометрическим данным на основе анализа поляриметрических сигнатур когерентности.

7. Предложен метод квази-поляриметрической съёмки с повторяющихся орбит, проведено моделирование искажений получаемой таким образом матрицы рассеяния; обработаны экспериментальные данные двух

14 частотных диапазонов, показана работоспособность метода, проведён сравнительный анализ источников декорреляции.

Научная новизна.

- Проведено оригинальное сравнительное исследование различных методов двумерной развёртки фазы на радиолокационных интерферограммах.

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

- Предложен метод построения поляриметрических сигнатур когерентноI сти и способ классификации на его основе. Проведён сравнительный анализ возможностей такой классификаций для данных Ь- и С-диапазона, а также для различных интервалов между съёмками.

- Предложен метод получения полной матрицы за два прохода носителя в схеме интерферометрической- съёмки; промоделированы возможные искажения и предложен способ их коррекции.

Достоверность

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

Практическая значимость

Практическая значимость работы заключается в том, что

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

- введённые поляриметрические сигнатуры интерферометрической когерентности являются новым признаком, позволяющим выполнять классификацию земных покровов на основе совместного анализа поляриметриче1

I ской и интерферометрической информации;

- предложенный метод формирования полной матрицы рассеяния по частичным поляризационным измерениям расширяет возможности поляриметрических систем наблюдения поверхности Земли.

Положения, выносимые на защиту

- Система критериев, обеспечивающая эффективное использование и сравнение методов устранения 27г-неоднозначности на интерферограммах.

- Модифицированный метод вычетов, устраняющий 2тг-неоднозначность на интерферограмме при условии пересечённого рельефа на снимаемой местности.

- Поляриметрические сигнатуры интерферометрической когерентности как новый классификационный признак типов земных покровов.

- Метод формирования полной матрицы рассеяния зондируемых сред по частичным поляризационным измерениям, выполняемым в двухпроход-ной интерферометрической схеме съёмки.

Личный вклад автора

Основные результаты диссертации получены либо лично автором, либо при его непосредственном участии. Автором лично разработана модификация одного из методов развёртки фазы, предложен метод классификации на основе поляриметрических сигнатур интерферометрической I когерентности, а также проведено моделирование возможных искажений матрицы рассеяний, полученной по частичным поляризационным измерениям, предложен метод их коррекции, проведена обработка экспериментальных данных для демонстрации работы предложенных методов. Автором разработан фрагмент пакета программ интерферометрической обработки радиолокационных данных в части этапа устранения^ 2тс-неоднозначности (двумерной развёртки фазы), а также полностью модуль поляриметрической интерферометрии. Интерпретация научных результатов осуществлялась вместе с соавторами публикаций, которым автор глубоко благодарен за плодотворную совместную работу. с

Публикации f

По результатам диссертационной работы опубликовано 35 работ, из них 7 в журналах, рекомендованных ВАК (ссылки [64], [65], [102], [155], [166], [169], [171] в общем списке литературы и [1—7] в Приложении).

Структура и объём диссертации

Диссертация состоит из введения, шести глав, заключения и списка литературы. Работа изложена на 213 страницах, содержит 51 рисунок, 6 таблиц и библиографический список из 182 наименований.

Основные результаты работы заключаются в следующем:

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

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

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

- проанализирована последовательность этапов интерферометрической обработки и показана их взаимосвязь;

- систематизированы основные направления поляриметрической интерферометрии:

- приведены классы и конкретные примеры прикладных задач, решаемых с помощью поляриметрической интерферометрии.

2. Предложена оригинальная система критериев сравнения методов устранения 27Г-неоднозначности на интерферограммах, проведён сравнительный анализ на её основе.

3. Предложена и обоснована модификация одного из локальных методов устранения фазовой 2л;-неоднозначности, учитывающая условия наличия областей низкой когерентности на интерферограмме.

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

ЗАКЛЮЧЕНИЕ

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

1. Систематизировано состояние радиолокационной интерферометрии к настоящему моменту: а) подробно освещены такие аспекты, как организация интерферомет-рической съёмки, приведены геометрические соотношения, иллюстрирующие связь разности фаз сигналов с топографическими характеристиками поверхности и величиной подвижек (глава 1), б) проанализированы этапы интерферометрической обработки радиолокационных данных, а также приведены соотношения, указывающие пределы точности измеряемых в интерферометрии величин — радиального смещение и относительных высот рельефа (глава 2);

2. Детально рассмотрены методы развёртки фазы, для одного из которых предложена модификация, ускоряющая работу алгоритма в условиях областей декорреляции и сложного рельефа на снимаемой местности; проведено сравнение в классах локальных и глобальных методов по основным характеристикам методов, а также на примере обработки экспериментальных данных (глава 3);

3. Представлены методы поляриметрической интерферометрии — возникшей чуть более десяти лет назад новой отрасли дистанционного зондирования, в рамках которой: а) приведены базовые соотношения величин, изложены основные существующие методы (глава 4); б) введены в рассмотрение сигнатуры когерентности — предложенная автором визуализация поляриметрико-интерферометрических данных, могущая служить основой для классификации (глава 5); в) предложен и обоснован способ регистрации полной матрицы рассеяния по частичным поляризационным измерениям в схеме интерфе-рометрической съёмки с повторяющихся орбит носителя (глава 6).

1. Graham, L.C. Synthetic interferometer radar for topographic mapping // Proc. of the IEEE. — 1974. — Vol.62, №6. — P. 763-768.

2. Zebker, H.A., Goldstein, R.M.Topographic Mapping From Interferometric Synthetic Aperture Radar Observations // Journ. of Geophysical Research. — 1986.

3. Vol. 91, № B5. — P. 4993-4999.

4. Zebker, H.A., Goldstein, R.M. Interferometric radar measurement of ocean surface current // Nature. — 1987. — vol. 328, №20. — P. 707-709.

5. Gabriel, A.K., Goldstein, R.M. Crossed orbit interferometry: theory and experimental results from SIR-B // International Journ. of Remote Sensing. — 1988.1. Vol. 9, №5. — P. 857-872.

6. Gabriel, A.K., Goldstein, R.M., Zebker, H.A. Mapping Small Elevation Changes Over Large Areas: Differential Radar Interferometry // Journ. of Geophysical Research. — 1989. — Vol. 94, № B7. — P. 9183—9191.

7. Li, F.K., Goldstein, R.M. Studies of Multibaseline Spaceborne Interferometric Synthetic Aperture Radars // IEEE Trans, on Geoscience and Remote Sensing.1990. — Vol. 28, № 1. — P. 88—97.

8. Rodriges, E., Martin, J. Theory and design of interferometric synthetic aperture radars // IEE Proc. F — 1992. — Vol. 139, № 2. — P. 147—159.

9. Zebker, H., Villasenor, J. Decorrelation in interferometric radar echoes // IEEE Trans, on Geoscience and Remote Sensing. — 1992. — Vol.30, №5. — P. 950—959.

10. Massonnet, D., et al. The displacement field of the Landers earthquake mapped by radar interferometry // Nature. — 1993. — Vol. 364, №8. — P. 138—142.

11. Goldstein, R.M., et al. Satellite radar interferometry for monitoring ice sheet motion: Application to an antarctic ice stream // Science. — 1993. — Vol. 262.1. P. 1525—1530.

12. Hagberg, J.O., Ulander, L.M.H. On the optimization of interferometric SAR for topographic mapping // IEEE Trans, on Geoscience and Remote Sensing. — 1993. — Vol. 31, №1. — P.303—306.

13. Massonnet, D., Rabaute, T. Radar interferometry: Limits and potential // IEEE Trans, on Geoscience and Remote Sensing. — 1993. — Vol. 31, №2. — P. 455—464.

14. Prati, C., Rocca, F., Monti Guarnieri, A. Topographie Capabilities of SAR exemplified with ERS-1 // Geo-Information-Systems. — 1994. — Vol. 7, №1. — P.17—22.

15. Zebker, H.A., et al. Mapping the world's topography using radar interferometry: The TOPSAT mission // Proc. of the IEEE. — 1994. — Vol. 82, №12. — P. 1774—1786.

16. Zebker, H.A., et al. Accuracy of topographic maps derived from ERS-1 interfe-rometric data // IEEE Trans, on Geoscience and Remote Sensing. — 1994. — Vol.32, №4. — P.823—836.

17. Massonnet, D., et al. Radar interferometric mapping of deformation in the year after the Landers earthquake // Nature. — 1994. Vol. 369. — P. 227—230. — May 19 1994.

18. Peltzer, G., Hudnut, K.W., Feigl, K.L. Analysis of coseismic surface displacement gradients using radar interferometry: New insights into the Landers earthquake // Journ. of Geophysical Research. — 1994. — Vol. 99, № Bll. — P.21971—21981.

19. Zebker, H.A., et al. On the derivation of coseismic displacement fields using differential radar interferometry: The Landers earthquake // Journal of Geophysical Research. — 1994. — Vol. 99, №B10. — P. 19617—19634.

20. Gateiii, F., et al. The wavenumber shift in SAR interferometry // IEEE Trans, on Geoscience and Remote Sensing. — 1994. — Vol. 32, №4. — P. 855—865.

21. Just, D., Bamler, R. Phase statistics of interferograms with applications to synthetic aperture radar // Applied Optics. 1994. - Vol. 33, №20. - P.4361-4368.

22. Lee, J.S., et al. Intensity and phase statistics of multilook Polarimetrie and interferometric sar imagery // IEEE Trans, on Geoscience and Remote Sensing. — 1994. —Vol. 32, №5. P. 1017—1028.

23. Ghiglia, D.C., Romero, L.A. Robust two-dimensional weighted and unweighted phase unwrapping that uses fast transforms and iterative methods // Jour, of the Opt. Soc. of America A. — 1994. — Vol. 11, №1, :107-117, January 1994.

24. Pritt, M.D., Shipman, J.S. Least-squares two-dimensional phase unwrapping using FFT's // IEEE Trans, on Geoscience and Remote Sensing. — 1994. — Vol. 32, №3.—P. 706—708.

25. Duchossois, G., Martin, P. ERS-1 and ERS-2 Tandem Operations // ESA Bulletin. — 1995. —№83, P. 54—60. —August 1995.

26. Vachon, P.W., et al. ERS-1 synthetic aperture radar repeat-pass interferometry studies: Implications for Radarsat // Canadian Journ. of Remote Sensing. — 1995. — Vol. 21, №4. — P. 441—454.

27. Hagberg, J.O., Ulander, L.M.H., Askne, J. Repeat-pass SAR interferometry over forested terrain // IEEE Trans, on Geoscience and Remote Sensing. — 1995. — Vol. 33, №2. — P.331—340.

28. Wegmiiller, U., et al. Land-surface analysis using ERS-1 SAR interferometry // ESA Bulletin. — 1995. —№81. — P.30-37.

29. Mouginis-Mark, PJ. Preliminary observations of volcanoes with the SIR-C radar // IEEE Trans, on Geoscience and* Remote Sensing. — 1995. — Vol. 33, №4. —P.934—941.

30. Wright, P A., Cordey, R.A., Rogers, G. Coastal application of ERS-1 SAR interferometry// Proc. of Intl. Geoscience and Remote Sensing Symp. (IGARSS), Florence, Italy, 10-14 July 1995. — 1995. — P. 547—549.

31. Carrasco, D., Alonso, J., Broquetas, A. Accuracy assessment of SAR interferometry using the ERS-1// Proc. of Intl. Geoscience and Remote Sensing Symp., Florence, Italy, 10-14 July 1995. — 1995. — P.781—783.

32. Goldstein, R. Atmospheric limitations to repeat-track radar interferometry // Geophysical Research Letters. — 1995. — Vol. 22, №18. — P.2517—2520.

33. Zivanovic, S.S., Foster, J.R., Welch, W.J. A new method for improving the in-terferometric resolution by compensating for the atmospherically induced phase shift // Radio Science. — 1995. — Vol. 30, № 4. — P. 877—884.

34. Spagnolini, U. 2-D phase unwrapping and instantaneous frequency estimation // IEEE Trans, on Geoscience and Remote Sensing. — 1995. — Vol 33, № 3. — P. 579—589.

35. Wegmüller, U., Werner, C. Retrieval of vegetation parameters with sar interfe-rometry // IEEE Trans, on Geoscience and Remote Sensing. — 1997. — Vol. 35, № 1. —P. 18—24.

36. Zebker, H.A., Rosen, P.A., Hensley, S. Atmospheric effects in interferometric synthetic aperture radar surface deformation and topographic maps // Journ. of Geophysical Research. — 1997. — Vol. 102, № B4. — P. 7547—7563.

37. Cloude, S.R., Papathanassiou, K.P. Polarimetrie SAR Interferometry // IEEE Trans, on Geoscience and Remote Sensing. — 1998. — Vol. 36, №5. — P. 1551-1565.

38. Massonnet, D. Capabilities and limitations of the interferometric cartwheel // Proc. of CEOS SAR Workshop, ESA-CNES, Toulouse, France, 26-29 October 1999. — 1999. —P.13.

39. Hellwich, O., Ebner, H. Geocoding SAR interferograms by least squares adjustment // ISPRS Journ. of Photogrammetry and Remote Sensing. — 2000. — Vol: 55-. — P. 277—288.

40. Rignot, E.J.M. Effect of Faraday rotation on L-band interferometric and Polarimetrie synthetic-aperture radar data // IEEE Trans, on Geoscience and Remote Sensing. — 2000. Vol. 38, №1. — P. 383—390.

41. Bamler, R. Interferometric stereo radargrammetry: Absolute height determination from ERS-ENVISAT interferograms // Proc. of Intl. Geoscience and Remote Sensing Symp., Honolulu, Hawaii, 24-28 July, 2000. — 2000. — Vol. 2. —-P. 742—745.

42. Ferretti, A., Prati, C, Rocca, F. Permanent scatterers in SAR interferometry // IEEE Trans, on Geoscience and Remote Sensing. — 2001. — Vol. 39, №1. — P. 8—20.

43. Xia, Y., Kaufmann, H., Guo, X. Differential SAR interferometry using corner reflectors // Proc. of the Intl. Geoscience and Remote Sensing Symp., Toronto, Canada, 24-28 June 2002. — 2002. — Vol.2. — P. 1243—1246.

44. Monti Guarnieri, A., et al. Multi-mode ENVISAT AS AR interferometry: techniques and preliminary results // IEE Proc. on Radar Sonar Navigation. — 2003.

45. Vol. 150, №3. — P. 193—200.

46. Nishimura, T., et al. A preliminary fault model of the 2003 July 26, M 6.4 northern Miyagi earthquake, northeastern Japan, estimated from joint inversion of GPS, leveling, and InSAR data // Earth Planets Space. — 2003. — Vol. 55.1. P. 751—757.

47. Daito, K., et al. L-band PS analysis: JERS-1 results and TerraSAR-L predictions// Proc. of Third International Workshop on ERS SAR Interferometry, (Fringe'03), Frascati, Italy, 1-5 Dec 2003. — 2004.

48. Bamler, R., Holzner, J. ScanSAR interferometry for RADARSAT-2 and RADARSAT-3 // Canadian Journ. of Remote Sensing. — 2004. — Vol. 30, №3.1. P.437—447.

49. Ferretti, A., et al. Sub-millimeter accuracy of InSAR time series: Experimental validation //IEEE Trans, on Geoscience and Remote Sensing. — 2007. — Vol. 45, №5. — P. 1142—1153.

50. Marinkovic, P., et al. InSAR quality control: Analysis of five years of corner reflector time series // Proc. of the Fifth International Workshop on ERS/Envisat SAR Interferometry (Fringe'07), Frascati, Italy, 26-30 Nov 2007. — 2007.

51. Sandwell, D.T., et al. Accuracy and Resolution of ALOS Interferometry: Vector Deformation Maps of the Father's Day Intrusion at Kilauea // IEEE Trans, on Geoscience and Remote Sensing. — 2008. — Vol. 46, №11. — P. 3524—3534.

52. Bombrun, L., et al. DEM Error Retrieval by Analyzing Time Series of Differential Interferograms // IEEE Geoscience and Remote Sensing Letters. -— 2009. — Vol. 6, №4. P. 830—834.

53. Stramondo, S., et al. X-, C-, and L-Band DInSAR Investigation of the April 6, 2009, Abruzzi Earthquake // IEEE Geoscience and Remote Sensing Letters. — 2011. —Vol. 8, №1. — P. 49—53.

54. Franceschetti, G., Lanari, R. Synthetic Aperture Radar Processing. — CRC Press. Boca Raton, L, NY, 1999. 328 pages.

55. Elachi, C., van Zyl, J. Introduction to the Physics and Techniques of Remote Sensing. —- Wiley Series in Remote Sensing. John Wiley & Sons, Inc. — 2nd edition. — 2006. — 552 pages.

56. Woodhouse, I.H. Introduction to Microwave Remote Sensing. — Taylor & Francis, Boca Raton, FL, 2006. — 370 pages.

57. Cumming, I., Wong, F. Digital Processing Of Synthetic Aperture Radar Data: Algorithms And Implementation. — Artech House Publishers, New York, 2005. — 632 pages.

58. Garestier, F., Le Toan, T. Estimation of the Backscatter Vertical Profile of a Pine Forest Using Single Baseline P-Band (Pol-)InSAR Data // IEEE Trans, on Geoscience and Remote Sensing. — 2010. — Vol. 48, № 9. — P. 3340—3348.

59. Garestier, F., Le Toan, T. Forest Modeling For Height Inversion Using SingleBaseline InSAR/Pol-InSAR Data // IEEE Trans, on Geoscience and Remote Sensing.—2010.—Vol. 48, №3, Part 2'. —P. 1528—1539.

60. Захаров, А.И:, Захарова, JI.H. Значимость информации о фазе отражённого сигнала при радиолокационном картировании земных покровов // Радиотехника. — 2003 . — №12, — С.70-73.

61. Арманд, Н:А., Захаров, А:И., Захарова, JI. Н. Космические радары с синтезированной апертурой в дистанционном зондировании Земли — современные системы и перспективные проекты // Исследование Земли из космоса.i2010 — №2. — С. 3-13.

62. Colesanti, С., et al! Generation of DEM with1 sub-metric, vertical accuracy from 30 ERS-ENVISAT pairs // Proc. of Third International Workshop on ERS SAR Interferometry, (Fringe'03), Frascati, Italy, 1-5 Dec 2003. — 2003.

63. Arnaud, A., et al. ASAR ERS Interferometric Phase Continuity // Proc. of the Intl. Geoscience and Remote Sensing Symp. (IGARSS'03), Toulouse, France, 21-25 July 2003. — 2003. — Vol. 2. — P. 1133—1135.

64. Krieger, G., et al. TanDEM-X: Mission Concept, Product Definition and Pertliformance Prediction // Proc. of the 5 European Conference on Synthetic Aperture Radars (EUSAR 2006). — 2006.

65. Hajnsek, I., Moreira, A. TanDEM-X: Mission and Science Exploration duringththe Phase A Study // Proc. of the 5 European Conference on Synthetic Aperture Radars (EUSAR 2006). — 2006.i a

66. Candela, L., Caltagirone, R. COSMO-SkyMed: mission definition, main applications, and products // Proc.of the Workshop on applications of SAR Polarime-try and Polarimetrie Interferometry (POLINSAR-2003), Frascati, Italy, January 14-16, 2003. —2003.

67. Geudtner, D., Seguin, G. Capabilities of Canada's planned RADARSAT conthstellation. Proc. of the 5 European Conference on Synthetic Aperture Radars (EUSAR 2006). —2006.

68. Zink, M., Krieger, G., Amiot, T. Interferometric Performance of a Cartwheel Constellation for TerraS AR-L // Proc. of Third International Workshop on ERS SAR Interferometry, (Fringe'03), Frascati, Italy, 1-5 Dec 2003. — 2003.

69. Massonnet, D. Capabilities and Limitations, of the Interferometric Cartwheel // IEEE Trans, on Geoscience and Remote Sensing. — 2001. — Vol. 39, №3. — P. 506—520.

70. Massonnet, D. The interferometric cartwheel: a constellation of passive satellites to produce radar images to be coherently combined // Int. J. of Remote Sensing. — 2001. —Vol. 22, № 12. —P. 2413—2430.

71. Fiedler, H., et al'. Analysis of bistatic configurations for spaceborne SAR InterfetVi —.rometry // Proc. of the 5 European Conference on Synthetic Aperture Radar (EUSAR 2002), 4-6 June 2002, Cologne, Germany. — 2002.

72. Ramongassié, S., et al. High resolution SAR based on Micro-satellite. Technological step versus performances for the next 10 years // Proc.of 4th European Conference on Synthetic Aperture Radars (EUSAR 2004). — 2004: — Vol. 1. —P. 55—58.

73. Li, Z., Bethel, J. Image Coregistration in SAR Interferometry // The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences. — Beijing, 2008. — Vol. XXXVII, Part Bl. — P.433—438.

74. Eineder, M., et al. Filtering of interferometric SRTM X-SAR data // Proc. of the Intl. Geoscience and Remote Sensing Symposium, Toronto, Canada, 24-28 June 2002. —2002. —Vol. 1. —P. 164—166.

75. Huangyin, Y, et al. A SAR interferogram filter based on the empirical mode decomposition method // Proc. of the Intl. Geoscience and Remote Sensing Symp., Sydney, Australia, 9-13 July 2001. — 2001. — Vol. 5 — P.2061-2063.

76. Ciuc, M., et al. Amplitude-driven coherence filtering in complex interferograms // Proc: of the Intl. Geoscience and; Remote Sensing Symposium, Toronto, Canada, 24-28 June 2002. — 2002. — Vol.6i — P: 3453—3455.

77. Trouve, E., Nicolas, J.-M., Maitre, H; Improving phase unwrapping techniques by the use of local frequency estimates // IEEE Trans, on Geoscience and Remote Sensing. — 1998. — Vol. 36, №6. —P. 1963—1972.

78. Goldstein, R.M., Werner, C.L. Radar interferogram filtering for geophysical applications // Geophys. Res. Letters. 1998.-Vol. 25, №21. - P. 4035-403883. Baran, I., et al. A Modification to the Goldstein Radar Interferogram Filter //

79. EE Trans, on Geoscience and Remote Sensing. — 2003. — Vol. 41, №9. — P. 2114—2118.

80. Lee, J.-S., Grünes, M.R., de Grandi, G. Polarimetrie SAR speckle filtering and its implication for classification // IEEE Transactions on Geoscience and Remote Sensing, vol; 37, no. 5, September 1999, P. 2363-2373.

81. Hajnsek, I., et al: The potential, of surface parameter; estimation from interfereth • metric SAR. // Proc. of the 4 European Conference on Synthetic Aperture Radar (EUSAR 2002), 4 6 June 2002, Cologne, Germany. — 2002.

82. Madsen, S.N.,, Martin, J.M., Zebker, H.A. Analysis and Evaluation of the NASA/JPL TOPS AR Across-Track Interferometric SAR System // IEEE Trans, on Geoscience and Remote Sensing. —1995. — Vol: 33^ №2. — P. 383—-391.

83. Bamler, R., Just, D. Phase statistics and decorrelation in SAR interferograms // Proc. of the International Geoscience and Remote Sensing Symposium (IGARSS'93), Tokyo, Japan, 18-21 August 1993. — 1993. — P. 980—984.

84. Shiping, S. DEM Generation Using ERS-1/2 Interferometric SAR Data // Proc. of the Intl. Geoscience and Remote Sensing Symposium, 2000 (IGARSS 2000), Hawaii, USA, 24-28 July, 2000. — 2000. — Vol.2. — P. 788—790:

85. Dall, J., et al. Using airborne SAR interferometry to measure the elevation of a ' Greenland ice cap. Proc. of the Geoscience and Remote Sensing Symposium,2000 (IGARSS 2000), Hawaii, USA, 24-28 July, 2000. — 2000. — Vol. 3. — P. 1125—1127.

86. Mouginis-Mark, P.J., et al. Topographic change on volcanoes from SRTM and other interferometric radars // Proc. of the Geoscience and Remote Sensing Symposium, 2001 (IGARSS'01), 9-13 July 2001, Sydney, Australia. — 2001.1. Vol. 2. — P. 757—758.

87. Koch, A., Heipke, C., Lohmann, P. Analysis of SRTM DTM: Methodology and Practical Results // Proc. of Geospatial Theory, Processing and Applications ISPRS Commission IV Symposium 2002, Ottawa, Canada, July 9-12, 2002.

88. Colesanti, C., et al. Monitoring known seismic faults using the permanent scat-terers (PS) technique. Proc. of the International Geoscience and Remote Sensingi

89. Usai, S., et al. Modelling terrain deformations at the Phlegrean Fields with INSAR. Proc. of the International Geoscience and Remote Sensing Symposium (IGARSS 2000), Hawaii, USA, 24-28 July, 2000. — 2000. — Vol. 5. — P. 2245—2247.

90. Ferretti, A., et al. SAR analysis of building collapse by means of the permanent scatterers technique. Proc. of the International Geoscience and Remote Sensing Symposium (IGARSS 2000), Hawaii, USA, 24-28 July, 2000. — 2000. — Vol. 7.—P. 3219—3221.

91. Захарова, Л.Н., Захаров, А.И. Сравнение некоторых современных методов разворота разности фаз в радиолокационной интерферометрии // Радиотехника и электроника. — 2003. — Т. 48, №10. — С. 1208—1213.

92. Costantini, М. A Novel Phase Unwrapping Method Based on Network Programming // IEEE Trans, on Geoscience and Remote Sensing. — 1998. — Vol. 36,№3.—P.813—821.

93. Ahuja, R.K., Magnanti, T.L., Orlin, J.B. Network Flows: Theory, Algorithms, and Applications // Prentice-Hall, Englewood Cliffs. —N.-J., 1993.

94. Eineder, M., Holzner, J. Phase unwrapping of low coherence differential interfe-rograms // Proc. of the International Geoscience and Remote Sensing Symposium (IGARSS'99), June 28 July 2 1999, Hamburg, Germany. — 1999. — Vol.3. —P. 1727—1730.

95. Zakharov, A.I., Tugarinov, P.V. On the Use of FFT for Phase Unwrapping of Interferogram With Poor Quality // Proc. of the 3rd European Conference on Synthetic Aperture Radars (EUSAR-2000), 23-25 May 2000, Munich. — 2000. — P. 509—512.

96. Suksmono, A. B., Hirose, A. A study of Interferometric SAR Image Restoration Using Complex-Valued Neural Networks and its Application to Phase Unwrapping Problem // Proc. of CEOS SAR Workshop. — 2001. — P. 29—33.

97. Zebker, H.A., Lu, Y. Phase Unwrapping Algorithms for Radar Interferometry: Residue-Cut, Least-Squares, and Synthesis Algorithms // Journ. of Opt. Soc. Of America A. — 1997. — Vol.15, №3. — P. 586—598.

98. Kadono, H., Takei, H., Toyooka, S. A Noise-immune Method of Phase Unwrapping in Speckle Interferometry // Optics and Lasers in Engineering. — 1997. — №26. — P. 151—164.

99. Li, X., et al. Phase Unwrapping of SAR Interferogram Based on Dyadic Wavelets // Proc. of the Intl. Geoscience and Remote Sensing Symp., Toronto, Canada, 24-28 June 2002. — 2002. — Vol. 4. — P. 2537—2539.

100. Hellwich, O. SAR Phase Unwrapping: Implications of Terrain Shape and Smoothing // Proc. of the 2nd European Conf. on synthetic aperture radars (EUSAR'98), 25-27 May, Friedrichshafen, Germany. — 1998 — P.51—56.

101. Tarayre-Oriot, H., Massonet, D. New methods of phase unwrapping in SAR interferometry // Proc. of the Workshop on ERS SAR Interferometry (Fringe'96), Zurich, 30 Sept. 2 Oct. 1996.

102. Smits, P.C., et al. Iterative Model Reconstruction for phase Unwrapping // Proc. 3rd ERS Symposium "Space at the service of our Environment", Florence, Italy, 17-21 March 1997. — 1997. — P. 1707—1710.

103. Fornaro, G., et al. Interferometric SAR Phase Unwrapping Using the Finite Elements Method// IEE Proc. Radar, Sonar and Navigation. — 1997. —Vol. 144, №5. — P.266—274.

104. Zakharova, L. Comparison Of Global And Local Approach To Phase Unwrapping For a Rugged Terrain // Proc. of Workshop on ERS SAR Interferometry (Fringe-2003), Frascati, Italy, December 1-5. — 2003.

105. Boerner, W.-M., et al. Chapter 5 "Polarimetry at Remote Sensing" in the Manual of Remote Sensin. 3rd Edition, ASPRS publishing, Bethesda, MD; 1997.

106. Strozzi T. et al; Measurement of slow uniform surface displacement with mm/year accuracy // Proc. of Intl. Geoscience and Remote Sensing Symp., Honolulu, Hawaii, 24-28 July, 2000. — 2000. — Vol. 5. — P. 2239—2241.

107. Papathanassiou, K.P., Reigber, A., Cloude, S.R. Vegetation and ground parameters estimation' using Polarimetrie Interferometry. Parts I, II. // Proc. of ESA CEOS Workshop, 1999. Toulouse, France. — 1999.

108. Dali, J., Papathanassiou, K., Skriver, H. Polarimetrie SAR Interferometry applied to land ice. First result // Proc. of Intl. Geoscience and Remote Sensing Symp. (IGARSS'03), Toulouse, France, 21-25 July 2003. — 2003. — Vol. 3. — P. 1432—1434.

109. Treuhaft, R, Siqueira, P. Vertical structure of vegetated land surfaces from in-terferometric and Polarimetrie radar // Radio Science. — 2000. — Vol 35, №1.1. P.141—177.

110. Cloude, S. R., Papathanassiou, K.P. A 3-Stage Inversion Process for Polarime-trie SAR Interferometry // Proc. of the 4 European Conference on Synthetic Aperture Radar (EUSAR 2002), 4-6 June 2002, Cologne, Germany. — 2002.

111. Cloude, S.R., Papathanassiou, K.P. Three-stage inversion process for Polarimetrie SAR interferometry // IEE Proc. Radar, Sonar and Navigation 2003. — Vol. 150, №3. —P. 125—134.

112. Cloude, S.R, et al. Multi-Frequency Polarimetrie SAR Interferometry for Vegetation Structure Extraction// Proc. of the Intl. Geoscience and Remote Sensing Symposium, 2000 (IGARSS 2000), Hawaii, USA, 24-28 July, 2000. — 2000.1. Vol.1. —P.129—131.

113. Gabriel, J, et al. Tree Height Extraction Using Polarimetrie SAR Interferometry // Proc. of the Intl. Geoscience and Remote Sensing Symp. (IGARSS'99), 28 June-2 July 1999, Hamburg, Germany. — 1999. — Vol. 4. — P. 2131—2133.

114. Wallington, E.D., et al. Assessing L-band SAR modes for commercial forest management // Proc. of the Intl. Geoscience and Remote Sensing Symp. (IGARSS'03), Toulouse, France, 21-25 July 2003. — 2003. — Vol. 4. — P. 2541—2543.

115. Colin, E., Titin-Schnaider, C., Tabbara, W. A new parameter for IFPOL coherence optimization methods // Proc. of the Intl. Geoscience and Remote Sensing Symp. (IGARSS'03), Toulouse, France, 21-25 July 2003. — 2003. — Vol. 6. — P. 3979—3981.

116. Shimada, M., et al. Tree height estimation using an airborne L-band Polarimetrie interferometric SAR // Proc. of the Intl. Geoscience and Remote Sensing Symp., Sydney, Australia, 9-13 July 2001. — 2001. — Vol. 3. — P.1430—1432.

117. Mette, T., et al. Forest biomass estimation using Polarimetrie SAR interferome-try // Proc. of the Intl. Geoscience and Remote Sensing Symp., Toronto, Canada, 24-28 June 2002. — 2002. — Vol.2. — P. 817—819.

118. Cloude, S.R., Papatanassiou, K.P., Reigber, A. Polarimetrie SAR Interferometry at P-band for Vegetation Structure Extraction // Proc. of the 3rd European Conference on Synthetic Aperture Radars (EUSAR-2000), 23-25 May 2000, Munich.2000.—P. 249—252.

119. Breit, H., et al. SRTM X-SAR DEM of Europe-Results and algorithmic improvements // Proc. of the Intl. Geoscience and Remote Sensing Symp., Toronto,

120. Canada, 24-28 June 2002: — 2002. — Vol.1. — P. 155—157.

121. Quartulli, M, Datcu, M. 3D scene reconstruction and structure recognition from metric resolution InSAR // Proc. of the 4th European Conference on Synthetic Aperture Radar (EUSAR 2002), 4-6 June 2002, Cologne, Germany. — 2002.

122. Wang Chao et al. Retrieval of DEM from SIR-C Data in Kunlun Mountain, West China // Proc. of the Intl. Geoscience and Remote Sensing Symp. (IGARSS'99), 28 June 2 July 1999, Hamburg, Germany. — 1999. — Vol. 2.1. P. 1357—1359.

123. Schuler, D.L., Lee, J.S., de Grandi, G. Measurement of Topographic Surface Slopes Using Polarimetrie SAR Images // Proc. of the Intl. Geoscience and Remote Sensing Symp. (IGARSS'94), 8-12 August 1994, Pasadena, CA, USA. — 1994. —Vol.3. —P. 1467—1469.

124. Li, X., et al. Generation and Error Analysis of DEM Using Spaceborne Polarimetrie SAR Interferometry Data // Proc. of the Intl. Geoscience and Remote Sensing Symp., Toronto, Canada, 24-28 June 2002. — 2002. — Vol.5. — P. 2705—2707.

125. Woodhouse, I.H., et al. Polarimetrie Interferometry in the Glen Affric Project: Results & Conclusions // Proc. of the Intl. Geoscience and Remote Sensing Symp., Toronto, Canada, 24-28 June 2002. — 2002. — Vol.2. — P. 820—822.

126. Balzter, H., et al. Potential of Polarimetrie SAR Interferometry for Forest Carbon Accounting // Proc. of the Intl. Geoscience and Remote Sensing Symp. (IGARSS'03), Toulouse, France, 21-25 July 2003. — 2003. — Vol. 3. — P. 1945—1947.

127. Lukowski, T.I., Charbonneau, F.J. Synthetic Aperture Radar and Search and Rescue // Proc. of the Intl. Geoscience and Remote Sensing Symposium, 2000 (IGARSS 2000), Hawaii, USA, 24-28 July, 2000. — 2000. — Vol.5. — P.2374—2376.

128. Guillaso, S., et al. Analysis of built-up areas from Polarimetrie interferometric SAR images // Proc. of the Intl. Geoscience and Remote Sensing Symp. (IGARSS'03), Toulouse, France, 21-25 July 2003. — 2003. — Vol. 3. — P. 1727—1729.

129. Sagues, L., et al. Wide-band Polarimetrie SAR Interferometry for Buried Mine Detection // Proc. of the 3rd European Conference on Synthetic Aperture Radars (EUSAR-2000), 23-25 May 2000, Munich. — 2000. — P. 527—530.

130. Захарова, JI. Сигнатуры когерентности в поляриметрической интерферометрии для классификации земных покровов // Радиотехника. — 2005. — №8. — С. 45—50.

131. Захарова, Л'. Применение методов поляриметрической интерферометрии для исследования вертикальной структуры леса // Тр. XXII Симп. по радиолокационному зондированию природных сред, Санкт-Петербург, 18—20 апреля,2004. — 2004.

132. Zakharova, L. Polarimetrie Coherence for Land Covers Classification // Proc. ofiL ^^^the 28 ' General Assembly of URS! (Union Radio-Scientifique Internationale), New Delhi, India, October 22-29, 2005. — 2005.

133. Спб., 2004. — Т. 1. — С. 65—69.