Исследование зависимости результатов обработки радиолокационных данных ДЗЗ от параметров обработки. Основные результаты обработки методом PS данных РСА PALSAR/ALOS по территории Московского региона

Размещено на http://www.allbest.ru/

ИССЛЕДОВАНИЕ ЗАВИСИМОСТИ РЕЗУЛЬТАТОВ ОБРАБОТКИ РАДИОЛОКАЦИОННЫХ ДАННЫХ ДЗЗ ОТ ПАРАМЕТРОВ ОБРАБОТКИ. ОСНОВНЫЕ РЕЗУЛЬТАТЫ ОБРАБОТКИ МЕТОДОМ PS ДАННЫХ РСА PALSAR/ALOS ПО ТЕРРИТОРИИ МОСКОВСКОГО РЕГИОНА

А.А. Феоктистов,

А.И. Захаров,

М. А. Гусев,

П.В. Денисов

Введение

В данной работе продолжены исследования практических возможностей метода постоянных рассеивателей с помощью его программной реализации в виде модуля PS в программном пакете SARscape. В статьях [1, 2] изложены первые результаты обработки тестового набора из 25 фрагментов радиолокационных снимков территории вблизи города Лас Вегас, США, полученных радиолокатором с синтезированной апертурой (РСА) ASAR/ENVISAT (C-диапазон, длина волны излучения 5, 6 см). Особенность обрабатывавшегося тестового массива изображений была в том, что снятая территория является полупустынной территорией (влияние атмосферных фазовых искажений и осадков предположительно относительно невелико; вследствие положительных температур снежный покров и процессы замерзания-оттаивания, приводящие к нежелательной динамике отражающей поверхности, практически не наблюдаются).

На новом этапе нами проведены экспериментальные исследования возможностей и ограничений метода PS в случае обработки длинной временной серии из 29 радиолокационных изображений, полученных в период с 18.06.2006 по 14.02.2011 по территории г. Москвы и Московской области с помощью РСА PALSAR/ALOS (L-диапазон, длина волны излучения 23, 6 см). Некоторые результаты анализа группы близко расположенных рассеивателей приведены в [3], где на примере гостиницы Дельта (район метро Партизанская) проведена оценка влияния условий измерения смещений частей корпуса гостиницы с помощью метода PS на точность измерений. Вследствие малого расстояния между постоянными рассеивателями влияние атмосферных фазовых искажений на измерения взаимного положения PS считалось несущественным [3]. В ходе анализа истории смещений группы компактно расположенных PS сделано заключение о том, что основные ошибки измерений динамики отражающей поверхности порождены: (1) ошибками в знании высоты рельефа отражающей поверхности или положения постоянных рассеивателей по высоте, (2) шумами приемной системы, (3) изменением свойств отражающей поверхности под воздействием атмосферных осадков и температуры.

В данной работе сообщается о результатах дополнительных исследований для произвольного количества PS с произвольным взаимным положением на поверхности, выполненных с использованием модуля PS программного пакета SARscape и опорной цифровой модели рельефа (ЦМР) SRTM V4. Ограничения, связанные с компактностью расположения анализируемых PS, здесь отсутствуют, а потому проблеме атмосферных фазовых искажений уделялось основное внимание. Исследования проводились в пределах границ тестового фрагмента размером 14х15 км². Границы области "покрытия" земной поверхности каждым из 29 снимков РСА, а также границы тестового фрагмента, показаны на рис. 1 красным цветом. Московская кольцевая автомобильная дорога отмечена желтым цветом. Значения основных параметров алгоритмов обработки с помощью модуля PS полностью соответствуют значениям из [2].

Рисунок 1 - Границы области покрытия земной поверхности изображениями РСА PALSAR/ALOS и тестового фрагмента

Замечание. Метод PS разработан для анализа постоянных рассеивателей (точечных объектов типа квазиуголковых отражателей); в рамках метода PS “работает” разреженная сетка пикселей; возможность формирования пространственно "плотных" выходных продуктов реализована в рамках метода малых базовых линий [4, 5].

Основные результаты

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

Рисунок 2 - Усредненное изображение тестового фрагмента

Один из вариантов схемы распределения массива PS при , полученных с помощью пакета SARScape при указанных выше параметрах обработки, приведен на рис. 3. Средняя плотность PS по поверхности в этом районе города PS/км².

Замечание. При обработке набора данных РСА С-диапазона ASAR/ENVISAT (часть 2, раздел 4.2 [2]) при было получено значительно более высокое значение PS/км². Более высокий уровень пространственного разрешения данных РСА PALSAR/ALOS (14 м вместо 30 м у данных РСА ASAR/ENVISAT) должен был привести к дополнительному значительному увеличению . Обнаруженное нами падение до 85 PS/км², по-видимому, возникло вследствие превалирующего влияния другого фактора. Длина волны излучения РСА PALSAR/ALOS более чем в 4 раза превышает длину волны излучения ASAR/ENVISAT, что приводит к существенному снижению уровня эффективной площади рассеяния квазиуголковых отражателей типичной городской застройки, исчезновению некоторых из них и результирующему падению плотности PS.

Рисунок 3 - Схема распределения PS в восточной части города

В качестве варианта 3D-визуализации расположения постоянных рассеивателей на рис. 4 представлен фрагмент схемы распределения PS в районе спорткомплекса Олимпийский (его местоположение обозначено на рис. 3 звездочкой). Данная опция должна обеспечивать однозначную идентификацию элементов инфраструктуры, соответствующих отдельным PS(путем визуального анализа с разных ракурсов; опция рекомендована для использования разработчиками программного пакета SARscape; важна при решении задач оперативного мониторинга критически важных объектов); однако вследствие ошибок в геометрических преобразованиях текущая версия программного обеспечения работает недостаточно удовлетворительно.

Рисунок 4 - Пример отображения в 3D-формате с использованием интерфейса Google Earth

Анализ зависимости историй смещения PS от количества изображений K

Ввиду отсутствия каких-либо результатов подспутниковых измерений значительное внимание уделено анализу относительных точностных характеристик выходных продуктов модуля PS; в полной аналогии с проведенным ранее анализом данных РСА ASAR/ENVISAT (часть 2, раздел 9.1 [2]) исследована зависимость историй смещения трех PS, взятых на разных участках тестового фрагмента, от значений одного из ключевых параметров обработки - количества изображений временной серии (история смещения - совокупность значений величин смещений, рассчитанных на момент съемки парных изображений каждой интерферометрической пары; характеристики набора данных представлены в [3]). Дислокация PSобозначена цифрами 1, 2, 3 на рис. 2. Критерий выбора PS - однозначная идентификация на схемах распределения PS, рассчитанных с использованием временных серий из изображений; при проведении обработки последовательно исключались изображения с порядковыми номерами 14 (второй этап обработки) и 2 (третий этап) [3].

Анализ графических зависимостей. На рис. 5 представлены графические зависимости для историй смещения трех PS для каждой из трех временных серий; зависимости построены с использованием функции линейной аппроксимации программы Excel 2010 (величины смещений из выходного файла модуля PS). Цвета - красный, зеленый, синий при , соответственно. Цветные пунктирные линии - линии регрессии (линейные аппроксимирующие зависимости); скорости смещения PS определяются тангенсами угла наклона прямых.

Необходимо отметить наличие проблем, связанных с возникновением на графиках как одиночных "выбросов", так и серий "выбросов" с экстремальным изменением значений смещения. Разности значений смещений PS в соседних отсчетах (амплитуды) экстремальных "выбросов", выделены красными звездочками, примерно равны половине длины волны излучения РСА:

(1)

Это является серьезным аргументом в пользу того, что причиной данного артефакта могут быть сбои в работе процедуры развертки фазы при формировании историй смещений PS (некорректно восстановленное количество компонент фазы, кратных 2р) [6].

Рисунок 5 - Истории смещения трех PS для серий снимков разной длины

Замечание. Можно отметить также, что аналогичные эффекты возникали и ранее - при обработке данных РСА ASAR/ENVISAT (часть 2, раздел 9.1 [2]). Амплитуда экстремальных "выбросов" в этом случае падает примерно в четыре раза - что полностью соответствует отношению длин волн излучения РСА PALSAR/ALOS и ASAR/ENVISAT - 5, 6 и 23, 6 см, соответственно.

Сопоставление рассчитанных значений скорости смещений PS. В таблице 1 (столбцы 2-4) представлены рассчитанные модулем PS значения скорости смещения трех PS для каждой из трех временных серий изображений.

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

Значения среднеквадратичных отклонений скорости смещения трех PS, , равны 0, 58; 0, 38; 1, 19 мм/год, соответственно; среднее значение . Поскольку значения скорости смещения пикселя с минимальным значением индекса рассеяния амплитуды для этих трех наборов оказались различными (столбец 5), выполнена "ручная" коррекция рассчитанных при значений скорости смещения всех трех PS (часть 1, раздел 3 [1]). Скорректированные значения скорости смещения трех PS представлены в столбцах 6-8; значения скорости смещения равны 0, 41; 0, 53; 0, 84 мм/год, соответственно; среднее значение . Таким образом, проведение коррекции привело к снижению с 0, 71 мм/год до 0, 6 мм/год, т.е. примерно на 16%.

Таблица 1 - Значения скорости смещений PS по данным модуля PS, мм/год

Далее также в полной аналогии с проведенным ранее анализом данных РСА ASAR/ENVISAT (часть 2, раздел 9.1 [2]) выполнены контрольные расчеты значений скорости смещений PS с использованием функции линейной аппроксимации программы Excel 2010. Предварительно в выходном файле модуля PS (со значениями смещений PS) исключены экстремальные"выбросы".

Полученные результаты представлены в таблице 2. Значения скорости смещения PS1, PS2 и PS3 равны, соответственно, 0, 34; 1, 32; 0, 40 мм/год. Среднее значение . Важно отметить - значения по данным программы Excel 2010 и модуля PS оказались близкими (0, 69 и 0, 6 мм/год, соответственно).

Таблица 2 - Значения скорости смещений PS по данным функции линейной аппроксимации, реализованной в программе Excel 2010, мм/год

Разница между значениями средней скорости смещения PS1, PS2 и PS3 по данным Excel и SARscape (последние после "ручной" коррекции значений скорости смещения) оказалась равной, соответственно, 3, 12; 2, 65 и 2, 85 мм/год. Это расхождение может быть объяснено двумя возможными причинами: (1) необходимостью введения порогового ограничения на амплитуду отклонения смещений PS от линии регрессии при проведении расчетов с использованием Excel и (2) различием алгоритмов расчета скорости смещений в SARscape и Excel (линейный/нелинейный; часть 2, раздел 9.1 [2]).

Априорные оценки точности расчетов () значений скорости смещения PS (часть 1, раздел 2.1.5 [1] и часть 2, раздел 6 [2]) при приведены в таблице 3. Можно констатировать, что: (1) среднее значение по данным модуля PS и Excel значительно меньше априорных оценок точности расчетов значений скорости смещения (различие примерно в три раза), (2) зафиксированные различия между значениями средней скорости смещения PS1, PS2, PS3 по данным Excel и SARscape лишь незначительно превышают уровень априорных оценок точности .

Таблица 3 - Априорные оценки точности расчетов () значений скорости смещения PS, мм/год

Анализ зависимости значений остаточных высот PS от количества изображений K

В таблице 4 представлены рассчитанные с помощью модуля PS значения остаточных высот для трех PS при ; значения () равны 0, 96; 0, 91; 1, 04 м, соответственно; среднее значение . Таким образом, вариабельность значений остаточных высот PS в данном случае выросла на два порядка в сравнении с результатами проведенного ранее анализа данных РСА ASAR/ENVISAT (часть 2, раздел 9.2 [2]).

Априорные оценки точности расчетов () значений остаточных высот PS (часть 1, раздел 2.1.5 [1] и часть 2, раздел 6 [2]) при приведены в таблице 5. Можно констатировать, что среднее значение достаточно близко к уровню априорных оценок точности , превышая их примерно на 50%.

Таблица 4 - Значения остаточных высот PS, м

Таблица 5 - Априорные оценки точности расчетов () значений остаточных высот PS, м

Остаточные атмосферные искажения

Считается [1, 2], что основной объем атмосферных фазовых искажений связан с низкочастотными пространственными компонентами, которые эффективно устраняются при проведении атмосферной коррекции с использованием атмосферных фазовых экранов (APS), рассчитываемых в рамках процедуры высокочастотной временной и низкочастотной пространственной фильтрации (при проведении расчетов использовано задаваемое по умолчанию значение линейных размеров низкочастотного фильтра (часть 2, раздел 2 [2]).

В данном разделе описаны результаты проведенного визуального анализа остаточных атмосферных искажений (после проведения атмосферной коррекции). Анализировались только "корректные" истории смещения; истории смещения, содержащие вызванные сбоями в работе процедуры развертки фазы экстремальные "выбросы", из рассмотрения исключены. В качестве "опорной" графической зависимости выбран график для PS1 при К=29 (синий график на рис 5.1). На рис. 6 представлены графические зависимости для историй смещения трех групп из шести PS. Все расчеты выполнены для временной серии из изображений.

Первая группа включает рассеиватель PS1 и случайную выборку из пяти PS, находящихся в локальной окрестности PS1. Вторая группа включает PS1 и случайную выборку из других пяти PS, находящихся на расстоянии порядка нескольких сотен метров друг от друга. Наконец, третья группа кроме PS1 включает случайную выборку из пяти PS, расположенных на расстоянии порядка полутора км друг от друга (что превышает линейные размеры низкочастотного пространственного фильтра ). Графические зависимости для первой группы PS -на рис 6.1, второй - на рис. 6.2, третьей - на рис. 6.3.

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

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

И, наконец, для PS, находящихся на расстоянии порядка полутора км друг от друга, рис. 6.3, явная коррелированность сохранилась только в двух позициях (синие звездочки), причем в этих случаях заметна кросс-корреляция с соответствующими фрагментами графических зависимостей на рис. 6.1. и 6.2. Необходимо подчеркнуть, что ситуации, когда сохраняется коррелированность для PS, находящихся друг от друга на расстояниях, превышающих линейные размеры низкочастотного фильтра , возникают вследствие недостаточной эффективности используемой процедуры низкочастотной пространственной фильтрации.

Количественная оценка эффективности алгоритма атмосферной коррекции

При проведении количественной оценки эффективности атмосферной коррекции использован эвристический параметр - эффективность алгоритма атмосферной коррекции , равный отношению:

(2)

- оценка среднеквадратичного отклонения () для атмосферных фазовых искажений, устраненных на этапе атмосферной коррекции:

(3)

Рисунок 6 - Графические зависимости для историй смещения PS1 и трех случайных выборок PS

где - значений фазы (в радианах) для пяти точек , рассчитанных по набору атмосферных фазовых экранов ; ; . Положение пяти точек указано цифрами на рис. 7, на котором в качестве примера приведено изображение одного из атмосферных фазовых экранов (после выполнения процедур низкочастотной пространственной и высокочастотной временной фильтрации, часть 2, раздел 7 [2]); минимальное и максимальное значения фазы равны 0, 22 и 1, 83 рад, что согласно (4) соответствует смещениям 4, 1 и 34, 4 мм.

Переход от значений фазы в радианах к смещениям в мм выполнен в (1) с использованием коэффициента:

(4)

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

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

Рисунок 7 - Пример атмосферного фазового экрана РСА

Результаты проведенной количественной оценки эффективности алгоритма атмосферной коррекции представлены в первой строке таблицы 6. Для сравнения во второй строке приведены аналогичные результаты расчетов по данным ASAR/ENVISAT (часть 2, раздел 7 [2]).

Таблица 6 - Результаты оценки эффективности атмосферной коррекции

Можно отметить два момента. Во-первых, полученные по результатам обработки данных РСА PALSAR/ALOS и ASAR/ENVISAT значения отношения оказались очень близкими (0, 045 и 0, 061) - иными словами, уровень остаточных атмосферных искажений оказался близким к значению - для обоих упомянутых РСА, длины волн которых резко различаются (5, 6 и 23, 6 см, соответственно). Возможно, это значение характеризует предельные возможности варианта алгоритма атмосферной коррекции, реализованного в модуле PS программного пакета SARscape. Разумеется, чтобы однозначно сделать столь серьезное заключение, необходимо провести дополнительные экспериментальные исследования.

Во-вторых, следует обратить внимание на достаточно сильное различие значений параметра эффективности алгоритма атмосферной коррекции , рассчитанных по данным РСА PALSAR/ALOS и ASAR/ENVISAT (3, 14 и 1, 31).

Возможная интерпретация: массив данных РСА PALSAR/ALOS получен по территории г. Москвы и Московской области, для которой характерны достаточно сильные атмосферные искажения, тогда как в тестовом массиве РСА ASAR/ENVISAT (по полупустынной территории в США) уровень атмосферных фазовых искажений должен быть существенно более низким (лишь незначительно превышая тот уровень, ниже которого алгоритм коррекции становится недостаточно эффективным).

Выводы

Представлены основные результаты дополнительных экспериментальных исследований практических возможностей метода PS на примере модуля PS программного пакета SARscape и данных РСА PALSAR/ALOS по территории Московского региона.

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

Установлено, что значения отношения остаточных атмосферных искажений к длинам волн излучения РСА PALSAR/ALOS и ASAR/ENVISAT достаточно близки - 0, 045 и 0, 061, соответственно (т.е. остаточных атмосферных искажений примерно равно для обоих РСА, длины волн которых существенно различаются - 5, 6 и 23, 6 см). Возможно, это значение характеризует предельные возможности реализованного в модуле PS варианта алгоритма атмосферной коррекции - необходимо проведение дополнительных экспериментальных исследований.

Установлено также существенное различие между значениями параметра эффективности алгоритма атмосферной коррекции , рассчитанными по данным РСА PALSAR/ALOS и ASAR/ENVISAT (3, 14 и 1, 31). Возможная интерпретация - массив данных РСА PALSAR/ALOS получен по территории г. Москвы и Московской области, где характерно достаточно сильное влияние атмосферных факторов, тогда как для полупустынной территории в США (РСА ASAR/ENVISAT) характерен существенно более низкий уровень атмосферных фазовых искажений, лишь незначительно превышающий тот уровень, ниже которого алгоритм коррекции становится недостаточно эффективным.

Установлено, что наиболее вероятной причиной возникновения экстремальных "выбросов" на историях смещений PS являются сбои в работе процедуры развертки фазы на этапе формирования историй смещения PS (амплитуда выбросов равна примерно ). Отмечается, что аналогичные эффекты возникали и ранее - при обработке данных РСАASAR/ENVISAT (амплитуда экстремальных "выбросов" в этом случае падает примерно в четыре раза - что полностью соответствует отношению длин волн излучения РСА PALSAR/ALOS и ASAR/ENVISAT).

Построены схемы распределения PS по территории тестового фрагмента; зафиксировано резкое падение средней плотности PS по сравнению с аналогичными результатами проведенной ранее обработки данных РСА ASAR/ENVISAT. Падение плотности PS объясняется четырехкратным различием длин волн РСА PALSAR/ALOS и ASAR/ENVISAT, приводящим к соответствующему снижению уровня эффективной площади рассеяния квазиуголковых отражателей, формируемых элементами типичной городской инфраструктуры.

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

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

рассеиватель программный пакет московский

Литература

1. Феоктистов А.А., Захаров А.И., Гусев М.А., Денисов П.В. Исследование зависимости результатов обработки радиолокационных данных ДЗЗ от параметров обработки. Часть 1. Ключевые моменты метода постоянных рассеивателей // Журнал радиоэлектроники [электронный журнал]. 2014. №12. URL: http://jre.cplire.ru/jre/dec14/5/text.html

2. Феоктистов А.А., Захаров А.И., Гусев М.А., Денисов П.В. Исследование зависимости результатов обработки радиолокационных данных ДЗЗ от параметров обработки. Часть 2. Экспериментальные результаты // Журнал радиоэлектроники [электронный журнал]. 2014. №12. URL: http://jre.cplire.ru/jre/dec14/6/text.html

3. Захаров А.И., Феоктистов А.А., Гусев М.А., Денисов П.В. Влияние условий проведения радарной съемки на точность измерений динамики урбанизированной территории методом постоянных отражателей. Труды VII Всероссийской научной конференции «Радиофизические методы в дистанционном зондировании сред» в рамках VI Всероссийских Армандовских чтений «Муром`2016». Муром, 31.05-02.06. 2016 г., Муром. Электронное издание: CD-ROM, с. 389-396.

4. Феоктистов А.А., Захаров А.И., Гусев М.А., Денисов П.В. Исследование возможностей метода малых базовых линий на примере модуля SBAS программного пакета SARscape и данных РСА ASAR/ENVISAT и PALSAR/ALOS. Часть 1. Ключевые моменты метода // Журнал радиоэлектроники [электронный журнал]. 2015. №9. URL: http://jre.cplire.ru/jre/sep15/1/text.html.

5. Феоктистов А.А., Захаров А. И., Гусев М.А., Денисов П.В. Исследование возможностей метода малых базовых линий на примере модуля SBAS программного пакета SARscape и данных РСА ASAR/ENVISAT и PALSAR/ALOS. Часть 2. Экспериментальные результаты // Журнал радиоэлектроники [электронный журнал]. 2015. №9. URL: http://jre.cplire.ru/jre/sep15/2/text.html.

6. Захаров А.И., Яковлев О.И, Смирнов В.М. Спутниковый мониторинг Земли. Радиолокационное зондирование поверхности // М.: КРАСАНД, 2012. - 248 С.

Размещено на Allbest.ru