Кашкин В.Б., Сухинин А.И. Дистанционное зондирование Земли из космоса. Цифровая обработка изображений

Подождите немного. Документ загружается.

мил каждого канала квантуется на

256

уровней, скорость формирования

информационного потока

Мбит/с.

Если бы для каждого канала применялся один фотоприемник,

то при указанных скоростях сканирования не удалось бы обеспечить

указанного разрешения. Столь высокое разрешение сканеров достиг-

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

ной вдоль направления движения спутника, и последовательного счи-

тывания информации с элементов линейки.

Спутник SPOT-3 (Франция). Орбита почти круговая высотой

820

км

с наклонением 98,7°, период обращения 101 мин.

На

спутнике установ-

лены два сканера Н

VR с

линейной разверткой, фотоприемниками слу-

жат 1728-элементные ПЗС-линейки, ориентированные перпендику-

лярно движению спутника, цветоделение осуществляется с помощью

призм. Возможен панхроматический режим, перекрывающий зеле-

ный, красный

ближний ИК-диапазоны (0,51—0,73 мкм). При этом ис-

пользуются

линейки, число пикселов в

них

достигает

6000,

разреше-

ние в надире

AL = 10

м.

В многоспектральном режиме реализуются три канала: 0,5—0,59,

0,61—0,68,0,79-0,89

мкм,

по две

линейки на канал,

AL = 20 м. В

панхро-

матическом режиме сигнал квантуется на 64 уровня, в многоспект-

вальном

—

на 256 уровней для каждого канала. Полоса обзора равна

км

обоих случаях. Если включены оба сканера HVR, то полоса об-

I Юра

равна

117

км с

3-километровым

перекрытием. Периодичность на-

?блюдения заданного района составляет

раз в

дней. Камера может

отклоняться от направления в

надир до

+27°,

поэтому полоса обзора мо-

жет смещаться вправо и влево

до

475 км от трассы спутника. При этом

периодичность наблюдения заданного района может составлять

1—4

су-

ток. Особенностью спутника SPOT является возможность получать

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

того же

участка на двух последовательных витках.

Спутник ERS (Европейское космическое агентство). Высота орбиты

798x782 км с наклонением 98,54° и периодом обращения 100,67 мин.

В состав

бортовой аппаратуры включена радиолокационная станция ми-

кроволнового зондирования AMI (Active Microwave Instrument), кото-

рая обеспечивает три режима работы.

Режим построения радиолокационных изображений подстилающей

поверхности с использованием синтезированной апертуры антенны.

(AMI-SAR image mode) применяется при наблюдении береговой зоны,

полярных

льдов,

при определении состояния поверхности моря, выяв-

лении особенностей геологического строения земной поверхности,

Изучении растительного покрова. Сигналы, отраженные от поверхно-

сти Земли, могут приниматься двумя антеннами, расположенными од-

на над

другой.

По разности фаз их сигналов (интерферометрический ме-

тод измерения) можно определять высоту наземных объектов с точно-

стью 10 м.

Аппаратура AMI имеет следующие характеристики: мощность пе-

редатчика 1270 Вт, частота 5,3 ГГц, длительность импульса

37,1

мкс, ши-

рина спектра излучаемых сигналов 15,5 МГц, вертикальная линейная

поляризация излучаемых и принимаемых волн, пространственное раз-

решение 30 и 100 м, ширина полосы обзора 100 км при угле падения из-

лучаемых электромагнитных волн на поверхность Земли 23° в центре по-

лосы. На кодирование каждого отсчета отводится 16-битное слово.

Режим изучения морских волн с использованием синтезирован-

ной апертуры антенны (AM1-SAR wave mode) обеспечивает определе-

ние направления и длины морских волн. Данный режим программно

включается каждые 200—300 км, позволяя получать изображения раз-

мером 6x6 км, по которым можно оценивать характеристики морских

волн.

этом режиме излучаемая мощность 540 Вт, частота 5,3 ГГц, по-

ляризация излучаемых и принимаемых волн вертикальная, точность из-

мерения длины волны ±25% вдиапазоне 100-1000 м, точность опреде-

ления направления волн ±20° вдиапазоне 0—180°.

Режим трехлучевого скаттерометра (AMI Scatterometer mode) пред-

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

ветров.

этом режиме три передающие антенны формируют три луча,

сканирующие в полосе шириной до 500 км, позволяя определять на-

правление и скорость ветра. Элементы разрешения размером 50x50 км

формируются с интервалом 25 км. Характеристики AMI при этом сле-

дующие: излучаемая мощность 540 Вт, частота 5,3 ГГц, точность опре-

деления скорости ветра ±2 м/с или 10% в диапазоне 2—24 м/с, точ-

ность определения направления ветра ±20° в диапазоне 0—360°, угол

падения лучей изменяется в пределах 27—58°, обеспечивая полосу об-

зора шириной 500 км. Скорость передачи информации 500 кбит/с.

Во всех трех режимах полоса обзора аппаратуры AMI смещена относи-

тельно трассы ИСЗ вправо по ходу движения.

В состав измерительной аппаратуры входит также радиолокаци-

онный высотомер RA(Radar Altimeter) на частоту 13,8 ГГц для опреде-

ления скорости ветра, измерения характерной высоты волн, топогра-

фии морской поверхности, ледяного покрова и поверхности суши,

построения контуров ледяных массивов, а также выявления границ

морских

льдов.

Высотомер может работать в режиме исследования оке-

ана (Ocean Mode), обеспечивая точность измерения скорости волн

2 м/с и точность измерения высоты волн 0,5 м в диапазоне 1-20 м

в пределах пятна размером

1,6—2,0

км, точность определения высоты

подъема поверхности моря 10 см.

режиме определения характеристик ледяного покрова (Ice Mode)

высотомер работает с более низким пространственным разрешением

(около

км). Этот режим используется при топографии ледяных покро-

вов,

определении типа льда и выявлении границ ледяного покрова.

Для периодических уточнений параметров орбиты, а также калиб-

ровки радиолокационного высотомера на ERS устанавливаются лазер-

ные уголковые отражатели.

Комплекс приборов ATSR (Along-Track Scanning Radiometer and

Microwave Sounder) включает радиометр оптического диапазона и двух-

канальное микроволновое устройство вертикального зондирования.

Радиометр предназначен для наблюдения поверхности моря и суши, из-

мерения их температуры, температуры верхней облачности и обеспечи-

вает прием излучения в спектральных каналах 0,65; 0,85; 1,27; 1,6; 3,7;

11 и 12 мкм с пространственным разрешением

км в надире.

Микроволновое устройство зондирования на частоты 23,8 и 36,5 ГГц

обеспечивает поступление информации дистанционного зондирования

о содержании аэрозолей в земной атмосфере, концентрации водяного

пара и облачных капель, а также о состоянии растительного покрова

Земли. Точность определения вертикального профиля концентрации во-

дяного пара в пределах пятна размером 20—25 км составляет 10%.

Спектрометр GOME (Global Ozone Monotoring Experiment) исполь-

зуется для построения вертикальных профилей концентрации озона

) и малых газовых компонентов (NO, N0

, BrO, Н

0) в тропосфере

и стратосфере, измерения потоков солнечного излучения, отражаемо-

го поверхностью Земли и рассеиваемого атмосферой. Прибор работает

в ультрафиолетовом диапазоне в спектральных каналах 0,24—0,295;

0,29—0,405; 0,4—0,605 и 0,59—0,79 мкм. Каждый канал содержит решет-

ку

детекторов из 1024 фотодиодов, температура которых поддерживает-

ся в пределах 39—41 °С термоэлектрическими охладителями. Верти-

кальное разрешение при определении концентрации озона 0

составляет

5 км. Полоса обзора изменяется от

120 до

960 км, а пространственное раз-

решение

—

от 40x40 до 40x320 км.

Аппаратура PRARE (Precise Range and Rate Equipment) обеспечива-

ет определение параметров орбиты спутника путем одновременной пе-

редачи двух радиосигналов с разной частотой на сеть специальных назем-

ных станций, измеренная разница времени прихода сигналов позволяет

выполнить коррекцию относительной дисперсии, которая обусловлена

влиянием ионосферы. Информация

дальности до спутника и его ради-

альной скорости передается обратно на борт спутника и накапливается

специальном бортовом запоминающем устройстве, а затем передается

в пункт приема информации при пролете над ним. Точность определе-

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

4—8

см.

Информация с аппаратуры AMI, работающей врежиме построения

радиолокационных изображений (AMI-SAR image mode), передается со

скоростью 105 Мбит/с на частоте 8,14 ГГц только в реальном масшта-

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

ского аппарата ERS, за исключением радиолокационной, объединяет-

ся в цифровой поток со скоростью около 1,1 Мбит/с и передается на

Землю на частоте 8,04

ГГц

либо записывается на бортовой магнитофон

для последующего воспроизведения в зоне радиовидимости станции

приема информации со скоростью 15 Мбит/с и передачи на той же ча-

стоте.

Спутник «Океан-О» (Россия—Украина). Высота орбиты 667 км, пери-

од обращения 97,98 мин, угол наклонения орбиты 98,03°. На спутнике

установлены две некогерентные РЛС БО, антенны которых направле-

ны соответственно влево и вправо от направления движения спутника.

Полоса обзора каждой РЛС 455 км, угол падения излучаемых электро-

магнитных волн на подстилающую поверхность 20,48°. Пространст-

венное разрешение в среднем 1,3 км поперек направления движения

спутника и 2,5 км вдоль него.

Многоканальный сканирующий микроволновый радиометр «Дель-

та-2Д» имеет 4 канала со средней длиной волны 0,8 мкм (1-й канал),

1,35 см (2-й канал), 2,25 см (3-й канал) и 4,3 см (4-й канал). Полоса об-

зора составляет 1126 км, размеры элементов разрешения на подстила-

ющей поверхности не более 17x22 км в 1-м канале, 28x37 км во 2-м,

49x65 км в 3-м и

х 120 км в 4-м.

Два микроволновых несканирующих радиометра на длины волнах

2,25 и 6 см направлены в надир, имеют полосу обзора около 130 км.

На спутнике «Океан-О» установлен ряд приборов, работающих

в оптическом диапазоне. Сканер МСУ-СК с конической разверткой

аналогичен сканеру МСУ-СК спутника «Ресурс-О». Многоканальное

сканирующее устройство высокого разрешения МСУ-В с полосой об-

зора 195 км имеет 8 спектральных каналов:

1-й-0,48-0,52

мкм;

2-й —0,54-0,61 мкм;

3-й

—

0,63-0,73 мкм;

4-й —0,78-0,92 мкм;

5-й-0,915-0,997 мкм;

6-й

—

1,47-1,62

мкм;

7-й-2,060-2,385 мкм;

8-й— 10,6-12,0 мкм.

Пространственное разрешение в 1—5-м каналах 50 м, в 6-м

—

100 м,

в 7-м

—

300 м, в 8-м

—

250 м.

Многоканальное сканирующее устройство малого разрешения

МСУ-М имеет полосу обзора 2000 км, пространственное разрешение

1,5x1,7

км в центре строки и 1,8х 1,7 км на краю строки.

Поляризационный спектрорадиометр «Трассер-О» в диапазоне

длин волн 411-809 нм имеет пространственное разрешение 45 км.

Для сбора информации с морских буев предназначена аппаратура

сбора и передачи информации «Кондор-2».

Кроме России, США, Франции, Украины, Европейского космиче-

ского агентства, спутники для дистанционных исследований Земли

запушены Японией, Индией, Китаем, Бразилией и некоторыми други-

ми странами.

1.4.4.

Искажения спутниковых сканерных изображений

Рассмотрим схему

формирования спутнико-

вых изображений

учетом

кривизны Земли. Пусть

спутник огибает поверх-

ность Земли по круговой

орбите на высоте Н. Это

равносильно движению

спутника по прямой,

но над поверхностью ци-

линдра. Кажущееся изоб-

ражение, которое строит

сканер спутника

цилин-

дрической или линейной

разверткой,

—

это проек-

ция на поверхность цилиндра

/-"(рис.

1.27).

В действительности же сканер отображает земную поверхность,

т.е.

внешнюю поверхность цилиндра G. Будем считать, что Земля

—

шар радиусом

6370 км; собственное вращение Земли за время по-

строения сканерного изображения учитывать не будем.

Пусть

L = SA —

расстояние от спутника Sдо некоторой точки зем-

ной поверхности при сканировании, AL

—

разрешение в надире,

Дер —

мгновенное поле зрения, а

—

угол сканирования (угол визирования),

под углом удуга АКвидна из центра

Земли.

Используя теорему синусов

и учитывая, что для всех значений угла а выполняется соотношение

(#+

)

ос <

Я, получаем

Рис.

1.27. Диаграмма формирования спутниковых

сканерных изображений

50/sin

Р =

/Ш/sin a; sin

р = (1

H/Ro)sin

а;

у =

к -

- а;

у= я

—

arcsin|(

+ H/R

)sm a] -a; L=

Rsin

y/sin a.

Длина

дуги

К,

отвечающей развертке внешней поверхности цилин-

дра С, равна у=АК=

R^y.

Длина дуги

А/Л",

соответствующей кажущему-

ся изображению, х

МК= На, откуда а

х/Н. Таким образом,

у=

R{n

- arcsin |(l+tf/R

) sin (x/H)] -x/H]. (1.15)

Используя (1.15), можно пересчитать координату х

вдоль

строки ка-

жущегося изображения в координату вдоль строки реального изобра-

жения.

Знание расстояния L позволяет уточнить пространственное разре-

шение. Минимальная дальность

min

при этом разрешение

AL =

АфЯ;

при максимальном угле сканирования а,

шх

дальность равна

тях

разре-

шение

AL,

= A(pL

max

вдоль направления движения спутника. Но поверх-

ность видна со спутника под углом

тах

поэтому поперек направления

движения

Д/,

AZ,|/cosa

max

Таким образом, на максимальной дально-

сти поперек направления движения пикселы оказываются сжатыми

в AL

/AL раза, вдоль

—

в AL^/AL раза. Кажущаяся полоса обзора

maK

Яа

тах

, реальная полоса С

тах

=2Л|З

тах

Для спутника NOAA Н = 870 км,

тах

55°, Лф = 1,26

•

10~

рад,

AL=\,\

км. Максимальная дальность L

max

= 1827 км, полоса обзора

G~ 3000 км, разрешение AL

2,31 км и AL

= 4,02 км. На максималь-

ной дальности поперек направления движения пикселы оказываются

сжатыми в 3,65 раза, вдоль него

—

2,1

раза. Две реки, текущие парал-

лельно направлению движения спутника

находящиеся на расстоянии

1,1

км друг

от друга,

надире были бы различимы, но на краю скана они

сливаются. Если не проводить коррекцию, то полоса обзора кажется

равной

F= 2

На

тш

1670 км.

На рис 1.28 показаны кажущееся и действительное изображения

двух сканов. За время одного скана спутник NOAA перемещается на

' " Кажущееся изображение

Реальное изображение

Рис.

1.28. Кажущееся и реальное изображения двух строк сканирования (спутник NOAA)

— 45 км -

- 60,5 км

1,1 км при разрешении в надире 1,1 км. На некотором расстоянии от на-

дира строки изображения начинают заметно перекрываться.

Определим область спутникового изображения с малыми искаже-

ниями как участок, где строки перекрываются не более чем на 10%, т.е.

примерно на 0,1 км. При этом Д1| = 1,3 км. Этот участок соответству-

ет кажущейся полосе обзора в ±475 км и действительной полосе обзо-

ра в +540 км от надира. Разрешение вдоль скана на границе участка со-

ставляет AL

= 1,54 км, вдоль скана пикселы вытянуты в 1,4 раза.

Для спутников с малой

полосой обзора (спутник

«Ресурс-01-3», сканер

МСУ-Э; спутник SPOT,

сканер HVR) при наблюде-

нии в надир геометричес-

кие искажения изображе-

ний за счет кривизны Земли

незначительны. Искажения

возникают при отклонении

оси сканирования от нади-

ра, которое может достигать

около 30° (перспективные

искажения). Это иллюстри-

рует рис. 1.29, построенный на основании (1.15), где в центре показа-

на полоса шириной 45 км

—

область кажущегося изображения при от-

клонении оси сканера МСУ-Э на

30°.

действительности изображение

несколько вытянуто вдоль направления сканирования (60,5 км).

При этом разрешение вдоль направления движения спутника состав-

ляет AL| =

м на максимальной дальности и AI, = 39 м на минималь-

ной, в то время как в надире AL\

35 м. В надире для МСУ-Э разреше-

ние вдоль линии сканирования Д£

= 45 м, при отклонении оси

сканирования на 30° AL

60 м.

Сканер МСУ-СК спутника

«Ресурс-01 -3»

имеет спирально-кониче-

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

ности Земли

дугу.

Так как расстояние от спутника до каждой точки дуги

одно и то же, то искажения сканерного изображения за счет кривизны

Земли не возникают. Строки на изображении также следует строить в ви-

де

дуг.

Однако такое изображение передается со спутника в том же фор-

мате, что и изображения от сканеров с цилиндрической

линейной раз-

верткой, у которых строки — прямые линии. Переход от линейного

формата к формату, учитывающему особенности сканера со спирально-

конической разверткой, при известных параметрах сканера несложен.

Направление

движения

спутника

Рис.

1.29. Искажение изображения со спутника

«Ресурс-01-3» при отклонении оси сканирования

на

30°

от надира

Рис.

1.30. Преобразование строк изображения, полу-

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

кой

Пусть

— строка

линейного формата, х

и у

—

координаты пик-

села М на этой строке

(рис.

1.30).

Пикселы ре-

ального изображения

должны лежать на дуге

AS в точке К

коорди-

натами и и v. Отрезок

х должен отобра-

жаться

на дугу

KB,

длина

которой равна Ла,

У?=

FB.

Таким образом,

х, откуда а

x/R.

Из

рис.

1.30 следует, что

и =

/?sina

/f sin x/R.

Поскольку уравнение окружности

+ (v —

,7o v=y+ (Л

—

- ы

)

1/2

. Чтобы выполнить преобразования, необходимо знать размер R,

который зависит от конструкции сканера и высоты полета спутника,

а также максимальное значение угла

тах

связанное с полосой обзора

соотношением

2/?sina

max

Радиолокационные станции бокового обзора наблюдают поверх-

ность Земли под углом к надиру, существенно превышающим углы от-

клонения оси сканирования оптических сканеров МСУ-Э

др.

Таким

образом, радиолокационным изображениям в еще большей степени

присущи те же перспективные искажения, что и упомянутым скане-

рам. На радиолокационных изображениях предметов на местности

также возникают некоторые специфические искажения. При фор-

мировании сканерных изображений в оптическом диапазоне источ-

ник энергии (Солнце) находится в одной точке небесной сферы,

а приемник на борту спутника

— в

другой. При использовании РЛС

БО и источник энергии (передатчик), и приемник находятся в одной

точке (на спутнике), наблюдение производится под углом визирова-

ния в 20—30° к горизонту. Может случиться так, что радиоволна упа-

дет на склон горы под прямым углом. В этом случае изображение

ближней части склона сожмется почти до нуля, а дальней окажется

в зоне радиотени и растянется (это можно заметить на рис.

1.15).

Полностью устранить такие искажения путем обработки на ЭВМ не

представляется возможным. В то же время на изображении, полу-

ченном оптическим сканером, оба склона имеют размеры, более близ-

кие к реальным.

Рассмотренные причины геометрических искажений спутниковых

изображений не единственно возможные. Солнечно-синхронные ор-

биты природоведческих спутников не пересекают ось вращения Зем-

ли,

а наклонены относительно нее. Поэтому если спутник движется с се-

вера на юг (нисходящий виток орбиты), то вверху изображения будет

не север, как на карте, а, например, север-северо-восток. К тому же во

время сеанса приема спутниковой информации Земля поворачивается

на некоторый угол (на

0,25°

за

мин ). Кроме «свежих» сканерных изо-

бражений, приходится обрабатывать архивные изображения, для кото-

рых неизвестны и орбитальные данные, и угол отклонения оси скани-

рования от надира. В этих случаях необходим другой метод

геометрической коррекции (см. п. 3.1.2).

КОНТРОЛЬНЫЕ ЗАДАНИЯ

На

какую длину

волны приходится максимум теплового излучения

от

лес-

ного пожара, если горящий лес имеет температуру 600 "С?

Используя рис. 1.5, оцените вегетационный индекс для кукурузы и сои.

Найдите на рис. 1.6 линии поглощения озона и углекислого газа в ИК-

диапазоне.

Пусть оптическая толщина атмосферы

т =

0,2.

Во

сколько раз отличают-

ся интенсивности принимаемого аппаратурой спутника излучения при на-

блюдении в надир и под углом

5 =

30"?

Каково мгновенное поле зрения сканеров спутников

NOAA

(высота ор-

биты Н= 870 км) и «Ресурс-01-3» (//= 650 км)?

6. Карте какого масштаба соответствует панхроматическое изображение со

спутника SPOT с разрешением

м?

Каковы период обращения

орбитальная скорость японского спутника

исследования природных

ресурсов

Jers-1,

если высота

орбиты

составляет

568

км?

8. С какой скоростью движется по поверхности Земли проекция спутника,

если высота орбиты 900 км?

9. Солнечная постоянная на орбите 1,4 кВт/м

, КПД солнечных батарей

0,15.

Какова мощность новых солнечных батарей спутника NOAA?

10.

Подсчитайте скорость спутника NOAA на орбите.

11.

Сколько надо делать сканов в секунду, чтобы обеспечить разрешение

в 50 м с полосой обзора 2400 км? Какова должна быть скорость передачи ин-

формации в этом случае?

12.

Какова скорость спутника «Ресурс» на орбите и скорость движения

проекции спутника на земную поверхность? С какой частотой следует качать

зеркало сканера, чтобы получить разрешение

в 45 м

при использовании скане-

ра с качанием зеркала вместо ПЗС-линейки?

13.

Какого размера зеркальный объектив следует установить на геостаци-

онарном спутнике, чтобы получить разрешение в

м в зеленом участке спе-

ктра?

14.

Подсчитайте скорость передачи информации со сканера МСУ-Э спут-

ника «Ресурс-01» при разрешении 35x45 м.

15.

Подсчитайте число элементов в фотоприемной линейке каждого кана-

ла сканеров MSS и ТМ спутника LANDSAT.

16.

Орбитальные элементы в формате NORAD спутника OCEANSAT (Ин-

дия) имеют вид

1 25758U 9Ч029С 99351.20895206.00001054 00000-0 27840-3 0 1586

2 25758 98.2916 85.3537 0001002 80.8953 279.2357 14.50869716 29712

К какому моменту времени относятся орбитальные элементы, каковы угол

наклонения орбиты, долгота восходящего узла, эксцентриситет орбиты, аргу-

мент перигея, средняя аномалия, среднее движение в оборотах за сутки? Вы-

числите период обращения спутника в минутах и высоту орбиты в километрах.

17.

Какой размер должна иметь антенна для приема потока информации

в 100 Мбит/с со спутника, находящегося на расстоянии 2000 км, если шумовая

температура МШУ составляет 100 К, а мощность передатчика 10 Вт?

18.

Каков максимальный угол сканирования прибора МСУ-СК спутника

«Ресурс-01-3», если полоса обзора составляет 600 км? Чему равно разрешение

в надире и на границе полосы обзора?

19.

Чему равно разрешение, даваемое прибором MSS спутника LAN DSAT

на границе полосы обзора?

20.

К каким дополнительным искажениям изображения приводит враще-

ние Земли в течение 15-минутного сеанса приема информации со спутника

NOAA на средних широтах?

ЛИТЕРАТУРА

Кондратьев

К.Я.,

Тимофеев

Ю.М. Метеорологическое зондирование атмо-

сферы из космоса. Л.: Гидрометеоиздат, 1978. 279 с.

Зуев В.Е., Креков

Г.М.

Оптические модели атмосферы.

J1.:

Гидрометеоиздат,

1986.256 с.

Хргиан А.Х. Физика атмосферы. М.: Изд-во МГУ, 1988. 327 с.

Гарбук СВ.,

Гершензон

В.Е. Космические системы дистанционного зонди-

рования Земли. М.: Сканэкс, 1997. 296 с.

Киенко Ю.П. Введение в космическое природоведение и картографирова-

ние.

М.: Картгецентр— Геодезиздат, 1994. 214 с.

Дистанционное зондирование: количественный подход: Пер. с англ. / Под

ред.

А.С. Алексеева. М.: Недра, 1983. 415 с.

Брюханов

А.

В.,

Господинов

Г.

В.,

Книжников Ю.Ф. Аэрокосмические методы

в географических исследованиях. М.: Изд-во МГУ, 1982. 231 с.

КацЯ.Г.,

ТевелевА.В., Полетаев А.И. Основы космической геологии. М.: Не-

дра, 1988.236 с.

Гонин

Г.

Б.

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

Л.:

Недра, 1989.320 с.