Построение цифровых моделей рельефа по результатам интерферометрической обработки радиолокационных изображений ALOS PALSAR

Ю.Б. Баранов, Ю.И.Кантемиров,  Е.В. Киселевский,  М.А. Болсуновский

Задача получения цифровых моделей рельефа (ЦМР) на территории объектов добычи, транспорта и хранения нефти и газа в настоящее время является одной из наиболее актуальных для ОАО «Газпром» в части картографического обеспечения деятельности. Это объясняется многими причинами. Готовятся к разработке новые месторождения, вводятся в разработку более глубокие продуктивные горизонты уже разрабатываемых месторождений, проектируются новые объекты наземного обустройства месторождений и нефтегазотранспортной инфраструктуры. Для ряда месторождений выполняется мониторинг смещений земной поверхности, вызванных подработкой залежей углеводородов, интерферометрическим методом, для чего также необходим высокоточный опорный рельеф. Наконец, цифровые модели рельефа используются в научных исследованиях, для моделирования геодинамических напряжений, анализа линеаментов, структурного и геоморфологического дешифрирования и т.д.

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

В данной статье, авторы предлагают подробнее рассмотреть технологию получения ЦМР на основе метода радиолокационной интерферометрии. Очевидными преимуществами получения модели современного рельефа на основе интерферометрической обработки радиолокационных данных являются:

— всепогодность радиолокационной съемки (нет зависимости от облачности и времени суток в отличие от оптической стерео-съемки);

— относительно низкая стоимость и широкая полоса захвата по сравнению с авиационной стереосъемкой и лазерным сканированием;

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

— в случае использования данных радиолокатора ALOS PALSAR с гарантированной субпиксельной точностью геопозиционирования возможно построение ЦМР без наземных контрольных точек, либо с минимальным их количеством;

— высокая точность и детальность получаемых результатов (в частности, для данных ALOS PALSAR по экспертной оценке авторов, результирующий рельеф близок рельефу топографической карты масштаба 1:25 000 при размере одной сцены порядка 70 х 70 км).

Кроме перечисленных преимуществ, отработка технологии получения ЦМР интерферометрическим методом представляет интерес в связи с планируемым запуском группировки спутников дистанционного зондирования Земли ОАО «Газком», среди которых будут и радарные.

В головном научном институте ОАО «Газпром» – ООО «ВНИИГАЗ» – в настоящее время выполняется полномасштабный анализ потенциала радиолокационных данных для решения различных задач нефтегазовой отрасли (мониторинг смещений земной поверхности, получение ЦМР, обнаружение нефтяных пятен в акваториях, дешифрирование геологических структур и т.д.).

Рассмотрим подробнее технологию получения ЦМР по радарным данным ALOS PALSAR на конкретном примере. В настоящее время авторами выполняется мониторинг смещений земной поверхности многопроходным радарным дифференциальным интерферометрическим методом на Заполярном газоконденсатонефтяном месторождении. Для такого мониторинга необходима многопроходная интерферометрическая цепочка радарных космоснимков и опорный входной рельеф. Выбор в пользу получения опорного рельефа также интерферометрическим методом (но по другим исходным снимкам) был сделан исходя из следующих соображений.

Чем выше точность опорного рельефа по высоте, тем увереннее определяются смещения и отфильтровываются фазовые помехи. Для максимально корректного разделения фазы опорного рельефа и фазы смещений земной поверхности необходима цифровая модель местности (ЦММ) по состоянию наиболее близкая времени интерферометрической съемки, поскольку в таком случае выделение высотной составляющей фазы будет более точным. В связи с этим ЦМР, построенная по результатам оцифровки топокарты, не была бы лучшим решением. С другой стороны, район Заполярного месторождения расположен севернее 600 северной широты, и, следовательно, ЦММ, полученная по результатам миссии SRTM, для этого района не доступна (и, в целом, ее пространственное разрешение – 90 м – недостаточно для использования ее в качестве опорного рельефа для отслеживания смещений). ЦМР GTOPO-30 (пространственное разрешение 1 км) также совершенно не подходит для решения поставленной задачи.

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

Выбор исходных данных

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

 Таблица 1

Диапазоны радиоволновой области электромагнитного спектра

Диапазон Частоты, ГГц Длины волн, см Спутниковые системы
X 5.20 – 10.90 2.75 – 5.77

(2.4 – 3.8)

USGS SLAR, TerraSAR-X,

Cosmo-SkyMed

 

C 3.9 – 6.2 3.8 – 7.6 ERS-1,2, ENVISAT-ASAR, RADARSAT-1,

RADARSAT-2

L 0.39 – 1.55 19.3 – 76.9

(15 – 30)

SIR-A,B, ALOS
P 0.225 – 0.391 40.0 – 76.9

(30 – 100)

AIRSAR

 

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

Первый вариант – использование результатов тандемной съемки с небольшой временной базой (например, архив снимков, сделанных спутниками ERS-1 и ERS-2). Плюсы: идеальная интерферограмма в большинстве случаев; высокая когерентность; пространственное разрешение 20-25 м – приблизительно в 4 раза выше, чем у ЦМР SRTM. Минусы: неполное покрытие территории РФ; отсутствие свежих архивных данных (это минус, в случае если необходима современная ЦМР); резкое падение когерентности в районе лесных массивов, характерное в целом для С-диапазона радиоволновой области электромагнитного спектра.

Второй вариант – использование результатов радиолокационной съемки со спутника ALOS (бортовым радиолокатором с синтезированной апертурой PALSAR). Основным плюсом использования данных ALOS PALSAR для получения информации о рельефе земной поверхности представляется длина зондирующей волны этого радиолокатора – 23 см, соответствующая L-диапазону. Соответственно в большинстве случаев снимается проблема временной декорреляции интерферометрических пар радарных снимков, вызванной влиянием растительности. В связи с этим интерферометрические пары снимков ALOS PALSAR характеризуются высокой когерентностью, даже, несмотря на то, что минимально возможный временной интервал между ними составляет 46 дней. При этом режим съемки PLR (полная поляризация излучения) позволяет дополнительно повысить когерентность и оптимизировать интерферограмму за счет съемки с полной поляризационной матрицей. В свою очередь, данные ALOS PALSAR FBS (однополяризационный режим съемки) характеризуются высоким пространственным разрешением – 7 м, что в 3-4 раза превышает разрешение данных тандемной съемки ERS-1 и 2.

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

Поиск  архивных данных ALOS PALSAR

 Подбор архивных данных с учетом характеристик, необходимых для обеспечения результатов работ был выполнен официальным российским дистрибьютором данных со спутника ALOS – компанией «Совзонд». На участок интереса – район месторождения Заполярное – было найдено несколько комплектов данных, подходящих для интерферометрической обработки. Соответствующая карта покрытия приведена на рисунке 1.

Отметим, что данные ALOS PALSAR FBD (двойная поляризация) покрывают район интереса на 80 % двумя интерферометрическими парами снимков, в то время как одна интерферометрическая пара снимков ALOS PALSAR FBS (одинарная поляризация) полностью покрывает район интереса и, кроме того, характеризуется более высоким пространственным разрешением (7 м). Интерферометрические данные ALOS PALSAR PLR (полная поляризационная матрица) также присутствуют в архиве на рассматриваемый участок, но они отсняты в весенний период – по опыту авторов не самый благоприятный для этой территории.

В итоге для построения рельефа была выбрана интерферометрическая пара снимков от 09.06.2006 и 25.07.2006 гг., сделанных в режиме FBS (поляризация H/H).

Рис. 1. Схема покрытия района Заполярного месторождения (красный многоугольник) интерферометрическими данными ALOS PALSAR. Два перекрывающихся прохода от 18.07.2007 и от 02.09.2007 сделаны в режиме FBD (двойная поляризация), а два перекрывающихся прохода от 09.06.2006 и от 25.07.2006 сделаны в режиме FBS (одинарная поляризация).
Рис. 1. Схема покрытия района Заполярного месторождения (красный многоугольник) интерферометрическими данными ALOS PALSAR. Два перекрывающихся прохода от 18.07.2007 и от 02.09.2007 сделаны в режиме FBD (двойная поляризация), а два перекрывающихся прохода от 09.06.2006 и от 25.07.2006 сделаны в режиме FBS (одинарная поляризация)

Краткие теоретические основы получения рельефа интерферометрическим методом

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

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

Рис. 2. Схема космосъёмки интерферометрической пары изображений. H1 и H2 – высота радиолокатора; R1 и R2 – дальность (путь зондирующей волны); θ1 и θ2 – углы между дальностью и высотой; B┴ и B║ - перпендикулярная и параллельная составляющие базовой линии
Рис. 2. Схема космосъёмки интерферометрической пары изображений.
H1 и H2 – высота радиолокатора; R1 и R2 – дальность (путь зондирующей волны); θ1 и θ2 – углы между дальностью и высотой; B┴ и B║ — перпендикулярная и параллельная составляющие базовой линии

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

Каждый радиолокационный снимок интерферометрической пары (или цепочки) содержит в себе амплитудный и фазовый слой. Амплитудный слой более пригоден для визуального анализа. Результирующая фаза Ф, полученная в ходе интерферометрической обработки фазовых слоев снимков интерферометрической пары, в общем случае состоит из следующих компонентов:

Ф = Фtopo + Фdef + Фatm + Фn; где:

Фtopo – фазовый набег за счет обзора топографии под двумя разными углами;

Фdef – фазовый набег за счет смещения поверхности в период между съемками;

Фatm – фазовый набег за счет различия длин оптических путей из-за преломления в среде распространения сигнала;

Фn – вариации фазы в результате электромагнитного шума.

Непосредственно интерферометрическая обработка пары радиолокационных изображений с целью получения ЦМР в общем случае состоит из нескольких базовых шагов:

  • Совмещение основного и вспомогательного радиолокационных изображений интерферометрической пары (в автоматическом режиме либо с ручным вводом контрольных точек);
  • Генерация интерферограммы, являющейся результатом комплексного поэлементного перемножения основного изображения и изображения, комплексно-сопряжённого к вспомогательному;
  • Разделение компонентов фазы Фtopo и Фdef за счет синтеза фазы рельефа;
  • Фильтрация интерферограммы, позволяющая в определенной степени уменьшить фазовый шум (помехи) за счет загрубления выходной ЦМР;
  • Получение файла когерентности для области перекрытия двух снимков, составляющих интерферометрическую пару, в значениях от 0 до 1 для каждой пары соответствующих друг другу пикселей;
  • Развертка фазы (процедура перехода от относительных значений фазы к абсолютным);
  • Коррекция базовой линии;
  • Преобразование абсолютных значений фазы в относительные либо абсолютные высотные отметки в метрах с получением на выходе ЦМР;

Практическая реализация (программные решения)

 Авторы, на основе анализа мирового рынка программных продуктов, позволяющих выполнять интерферометрическую обработку радиолокационных снимков, сделали выбор в пользу специализированного программного комплекса SARscape (разработчик – SARMAP, Швейцария), распространяемого в качестве дополнительного модуля к программному комплексу ENVI (разработчик – ITT Visual Solutions, США). Эксклюзивным дистрибьютором обоих программных комплексов на территории Российской Федерации и стран СНГ является компания «Совзонд». ООО «ВНИИГАЗ» является пользователем обоих пакетов программ (ENVI + SARscape) с 2007 г. Программный комплекс SARscapeпозволяет проводить полную комплексную обработку материалов радиолокационной съемки. В нем реализованы наиболее современные алгоритмы обработки и поддержка всех видов радиолокационных данных, доступных пользователям на сегодняшний день.

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

  • Процесс совмещения основного и вспомогательного радиолокационных изображений интерферометрической пары радарных снимков с помощью двухэтапного автоматизированного алгоритма (первый этап – совмещение на пиксельном уровне, второй этап – на субпиксельном) с последующим расчетом интерферограммы и заданием параметров некогерентного накопления;
  • Процесс разделения компонентов фазы Фtopo и Фdef выполняется на основе синтеза фазы рельефа (с помощью имеющейся более грубой ЦМР либо набора наземных контрольных точек с известными абсолютными высотами, либо – в случае отсутствия таких данных – с помощью модели эллипсоида) и последующего ее вычитания из общей фазы с получением на выходе так называемой дифференциальной интерферограммы;
  • Процесс фильтрации интерферограммы (усредняющим либо адаптивными фильтрами), позволяющий в определенной степени уменьшить фазовый шум за счет загрубления выходной ЦМР, выполняется одновременно с получением файла когерентности для области перекрытия двух снимков, составляющих интерферометрическую пару;
  • Развертка фазы (процедура перехода от относительных значений фазы к абсолютным) по алгоритму растущей области с отсечением по порогу когерентности;
  • Коррекция базовой линии выполняется с помощью наземных контрольных точек либо специальным адаптивным алгоритмом на основе сопоставления развернутой фазы, фильтрованной интерферограммы, карты когерентности и уже имеющейся более грубой ЦМР;
  • Преобразование абсолютных значений фазы в относительные либо абсолютные высотные отметки в метрах с получением на выходе ЦМР;

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

Получение ЦМР интерферометрическим методом на примере Заполярного месторождения

Рассмотрим конкретный пример получения ЦМР интерферометрическим методом по паре снимков Заполярного месторождения, сделанных аппаратурой ALOS PALSAR 09.06.2006 и 25.07.2006 гг.

В качестве исходных данных были использованы два вышеуказанных снимка ALOS PALSAR уровня обработки 1.0 (сырые данные). Их фокусировка (приведение сырых данных к виду Single Look Complex) была выполнена в дополнительном модуле программного комплекса SARscape – SARscape Focusing module.

Сфокусированные снимки были успешно совмещены в полуавтоматическом режиме (необходимые параметры процесса автоматического совмещения задаются пользователем). Далее была рассчитана интерферограмма (рис. 3). Отметим высокое качество интерферограммы – наблюдаются непрерывные интерферометрические полосы. Зашумленная область в правом верхнем углу интерферограммы соответствует водной поверхности (р. Таз).

Для выделения топографической компоненты фазы из интерферограммы были использованы две имеющиеся на район исследований более грубых ЦМР. Они были построены по горизонталям и высотным отметкам, оцифрованным с топокарт масштабов   1 : 500 000 и 1 : 100 000 соответственно. На рисунках 4а и 4б представлены две соответствующие дифференциальные интерферограммы (прошедшие процедуру незначительной адаптивной фильтрации) на фрагмент участка исследований. Как и следовало ожидать, рельефная компонента более полно исключена из дифференциальной интерферограммы, полученной с использованием стотысячной топокарты (рис. 4б). На дифференциальной интерферограмме, полученной с помощью пятисоттысячной топокарты, рельефная компонента все еще довольно ярко выражена (рис. 4а).

Отметим, что перпендикулярная составляющая базовой линии для рассматриваемой пары снимков составляет около 2006 м. При таком значении базовой линии рельефная компонента выражена довольно сильно, и, чтобы извлечь ее полностью, даже стотысячной топокарты оказалось недостаточно. А в целом, для получения рельефа интерферометрическим методом по радарным снимкам в L-диапазоне, необходимы базовые линии примерно на порядок большие, чем для С-диапазона. Так, если для радарных спутников С-диапазона (ERS-1 и 2, ENVISAT, RADARSAT-1) приемлемыми для получения рельефа являются перпендикулярные базы в первые сотни метров, то для L-диапазона (ALOS PALSAR) таковыми являются базы в первые тысячи метров.

Рис. 3. Район Заполярного месторождения. Интерферограмма, рассчитанная по паре снимков ALOS PALSAR от 09.06.2006 и 25.07.2006
Рис. 3. Район Заполярного месторождения. Интерферограмма, рассчитанная по паре снимков ALOS PALSAR от 09.06.2006 и 25.07.2006
Рис. 4а. Район Заполярного месторождения. Дифференциальная интерферограмма, полученная путем вычитания из интерферограммы фазы опорного рельефа, синтезированной с помощью топокарты масштаба 1 : 500 000.  Рис. 4б. Район Заполярного месторождения. Дифференциальная интерферограмма, полученная путем вычитания из интерферограммы фазы опорного рельефа, синтезированной с помощью топокарты масштаба 1 : 100 000
Рис. 4а. Район Заполярного месторождения. Дифференциальная интерферограмма, полученная путем вычитания из интерферограммы фазы опорного рельефа, синтезированной с помощью топокарты масштаба 1 : 500 000.
Рис. 4б. Район Заполярного месторождения. Дифференциальная интерферограмма, полученная путем вычитания из интерферограммы фазы опорного рельефа, синтезированной с помощью топокарты масштаба 1 : 100 000

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

Построенная карта когерентности, говорит о достаточно высоких ее значениях, несмотря на разницу в 46 дней между съемками в летний период. Это на практике подтверждает теоретическую предпосылку о просвечивании растительности радио-лучами L-диапазона. На рисунке 5 представлена карта когерентности для рассматриваемой пары. Красным выделены пиксели со значением когерентности свыше 0,15. Отметим, что для интерферометрической пары снимков ENVISAT на эту же территорию и за приблизительно такой же период (35 дней) средняя когерентность ниже примерно в три раза.

Далее была развернута фаза и выполнена коррекция базовой линии. Для коррекции использовался файл наземных контрольных точек, состоящий из 10 точек. Из них 8 точек соответствовали высотным отметкам стотысячной топокарты, а 2 точки соответствовали положениям GPS-станций. После коррекции базовой линии построена результирующая ЦМР. Фрагмент этой ЦМР представлен на рисунке 6б в сравнении с ЦМР, полученной по результатам оцифровки горизонталей и высотных отметок стотысячной топокарты (рис. 6а). Визуально детальность ЦМР, полученной по данным ALOS PALSAR FBS, по информативности значительно превосходит стотысячный масштаб и по нашей оценке соответствует масштабу 1 : 25 000.

Наконец, была выполнена проверка точности по высоте полученной интерферометрическим методом ЦМР. Для этого в программном комплексе ENVI было выполнено вычитание из интерферометрической ЦМР стотысячного опорного рельефа. Классифицированная карта погрешности по высоте интерферометрической ЦМР (по сравнению с ЦМР, полученной из стотысячной топокарты) приведена на рисунке 7. Как видно из рисунка 7, масштабу 1 : 25 000 по высоте соответствует порядка 70 % результирующей ЦМР. Систематическое ухудшение точности по высоте – до точности соответствующей масштабу 1 : 100 000, наблюдается в северо-западной части полученной ЦМР, что может быть связано с отсутствием опорных наземных точек в этом районе и соответственно не совсем правильными результатами коррекции базовой линии на этом участке.

Рис. 5. Район Заполярного месторождения. Карта когерентности, рассчитанная по паре снимков ALOS PALSAR от 09.06.2006 и 25.07.2006. Красным цветом выделены пиксели с значение когерентности выше 0,15
Рис. 5. Район Заполярного месторождения. Карта когерентности, рассчитанная по паре снимков ALOS PALSAR от 09.06.2006 и 25.07.2006. Красным цветом выделены пиксели с значение когерентности выше 0,15
Рис. 6а. Фрагмент ЦМР, построенной по результатам оцифровки горизонталей и высотных отметок топокарты масштаба 1 : 100 000 Рис. 6б. Фрагмент ЦМР, построенной по результатам интерферометрической обработки пары снимков ALOS PALSAR от 09.06.2006 и 25.07.2006 с использованием в качестве опорного рельефа ЦМР, показанной на рис. 6а
Рис. 6а. Фрагмент ЦМР, построенной по результатам оцифровки горизонталей и высотных отметок топокарты масштаба 1 : 100 000
Рис. 6б. Фрагмент ЦМР, построенной по результатам интерферометрической обработки пары снимков ALOS PALSAR от 09.06.2006 и 25.07.2006 с использованием в качестве опорного рельефа ЦМР, показанной на рис. 6а
Рис. 7. Карта погрешности по высоте ЦМР, полученной интерферометрическим методом по данным ALOS PALSAR FBS. Красный цвет – погрешность по высоте в пределах 5 м, зеленый цвет – в пределах 10 м, синий цвет – в пределах 20 м, белый и черный цвета – разрыв фазы
Рис. 7. Карта погрешности по высоте ЦМР, полученной интерферометрическим методом по данным ALOS PALSAR FBS. Красный цвет – погрешность по высоте в пределах 5 м, зеленый цвет – в пределах 10 м, синий цвет – в пределах 20 м, белый и черный цвета – разрыв фазы
Рис. 8. Район Заполярного месторождения. 3D-визуализация в увеличенном вертикальном масштабе цифровой модели рельефа, полученной по результатам обработки интерферометрической пары снимков ALOS PALSAR от 09.06.2006 и 25.07.2006 (оценочный масштаб 1 : 25 000), совмещенной с детальным оптическим снимком со спутника Quickbird (пространственное разрешение 60 см). По центру рисунка  – крупный выпученный ледяной массив
Рис. 8. Район Заполярного месторождения. 3D-визуализация в увеличенном вертикальном масштабе цифровой модели рельефа, полученной по результатам обработки интерферометрической пары снимков ALOS PALSAR от 09.06.2006 и 25.07.2006 (оценочный масштаб 1 : 25 000), совмещенной с детальным оптическим снимком со спутника Quickbird (пространственное разрешение 60 см). По центру рисунка – крупный выпученный ледяной массив

 

Комментарии