ГОСТ Р 59759-2021 Данные дистанционного зондирования Земли из космоса. Радиометрическая коррекция данных дистанционного зондирования Земли из космоса, получаемых с космических аппаратов оптико-электронного наблюдения в видимом и ближнем инфракрасном диапазоне. Требования к алгоритмам

ГОСТ Р 59759-2021

НАЦИОНАЛЬНЫЙ СТАНДАРТ РОССИЙСКОЙ ФЕДЕРАЦИИ

Данные дистанционного зондирования Земли из космоса

РАДИОМЕТРИЧЕСКАЯ КОРРЕКЦИЯ ДАННЫХ ДИСТАНЦИОННОГО ЗОНДИРОВАНИЯ ЗЕМЛИ ИЗ КОСМОСА, ПОЛУЧАЕМЫХ С КОСМИЧЕСКИХ АППАРАТОВ ОПТИКО-ЭЛЕКТРОННОГО НАБЛЮДЕНИЯ В ВИДИМОМ И БЛИЖНЕМ ИНФРАКРАСНОМ ДИАПАЗОНЕ

Требования к алгоритмам

Remote sensing data of the Earth from space. Radiometric correction for remote sensing data of the Earth from space obtained from optical-electronic observation satellites in the visible and near infrared range. Requirements for algorithms

ОКС 35.240.70

49.140

Дата введения 2022-05-01

Предисловие

1 РАЗРАБОТАН по заказу Государственной корпорации по космической деятельности "Роскосмос" Автономной некоммерческой организацией высшего образования "Университет Иннополис"

2 ВНЕСЕН Техническим комитетом по стандартизации ТК 321 "Ракетно-космическая техника"

3 УТВЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Приказом Федерального агентства по техническому регулированию и метрологии от 19 ноября 2021 г. N 1514-ст

4 ВВЕДЕН ВПЕРВЫЕ

Правила применения настоящего стандарта установлены в статье 26 Федерального закона от 29 июня 2015 г. N 162-ФЗ "О стандартизации в Российской Федерации". Информация об изменениях к настоящему стандарту публикуется в ежегодном (по состоянию на 1 января текущего года) информационном указателе "Национальные стандарты", а официальный текст изменений и поправок - в ежемесячном информационном указателе "Национальные стандарты". В случае пересмотра (замены) или отмены настоящего стандарта соответствующее уведомление будет опубликовано в ближайшем выпуске ежемесячного информационного указателя "Национальные стандарты". Соответствующая информация, уведомление и тексты размещаются также в информационной системе общего пользования - на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет (www.rst.gov.ru)

Введение

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

- относительная радиометрическая коррекция;

- абсолютная радиометрическая коррекция;

- атмосферная коррекция;

- коррекция влияния анизотропии земной поверхности.

Все перечисленные виды радиометрической коррекции изменяют значения пикселей данных дистанционного зондирования Земли из космоса, не затрагивая их геометрических свойств.

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

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

1 Область применения

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

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

2 Нормативные ссылки

В настоящем стандарте использована нормативная ссылка на следующий стандарт:

ГОСТ Р 53460 Глобальная справочная атмосфера для высот от 0 до 120 км для аэрокосмической практики. Параметры

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

3 Термины и определения

В настоящем стандарте применены следующие термины с соответствующими определениями:

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

3.2 энергетическая яркость: Отношение потока излучения к произведению телесного угла, в котором он распространяется, и проекции площади излучающего элемента на плоскость, перпендикулярную к направлению излучения.

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

3.4 спектральная плотность энергетической яркости: Отношение энергетической яркости, приходящейся на малый спектральный интервал, к ширине этого интервала.

3.5 спектральная плотность облученности: Отношение облученности, приходящейся на малый спектральный интервал, к ширине этого интервала.

3.6 эффективное значение спектральной плотности энергетической яркости: Средневзвешенное по функции спектральной чувствительности спектрального канала значение спектральной плотности энергетической яркости.

3.7 коэффициент спектральной плотности энергетической яркости: Отношение эффективного значения спектральной плотности энергетической яркости излучения, отраженного от изучаемой поверхности, к эффективному значению спектральной плотности энергетической яркости излучения, отраженного от совершенного рассеивателя, при тех же условиях облучения и наблюдения.

3.8 совершенный рассеиватель: Поверхность, которая полностью отражает падающий поток равномерно во всех направлениях.

3.9 ламбертовский рассеиватель: Поверхность, которая отражает падающий поток равномерно во всех направлениях.

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

3.11 анизотропия земной поверхности: Зависимость отражательных свойств земной поверхности от направления облучения и направления наблюдения этой поверхности.

3.12 зенитный угол объекта: Отсчитываемый от зенита угол между направлением на объект из точки наблюдения и направлением отвесной линии в этой точке.

3.13

3.14

4 Сокращения

В настоящем стандарте применены следующие сокращения:

5 Требования к алгоритму относительной радиометрической коррекции

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

5.2 Входными данными для алгоритма относительной радиометрической коррекции являются:

- исходные данные ДЗЗ из космоса для каждого спектрального канала;

- результаты калибровки ФПУ КА, представленные в виде файла(ов) справочных таблиц (LUT);

- метаданные, содержащие служебную информацию о текущем состоянии ФПУ КА.

5.3 Значения пикселей исходных данных ДЗЗ из космоса должны быть представлены цифровыми отсчетами, полученными АЦП на борту КА путем преобразования аналогового сигнала ФПУ.

5.4 Файл(ы) с результатами калибровки ФПУ может (могут) содержать следующие сведения:

- коды, информирующие о работоспособности каждого цифрового детектора;

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

- коэффициенты, описывающие погрешность линейности каждого цифрового детектора в каждом спектральном канале;

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

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

- значения цифровых отсчетов, соответствующие темновому сигналу;

- коэффициенты, определяющие зависимость спектральных коэффициентов преобразования каждого цифрового детектора от температуры фокальной плоскости ФПУ КА;

- опорное значение температуры фокальной плоскости ФПУ КА;

- иные параметры калибровки ФПУ КА, зависящие от его конструктивных особенностей.

5.5 Метаданные, содержащие служебную информацию, должны включать в себя следующие сведения о текущем состоянии ФПУ КА:

- время проведения съемки в системе UTC;

- температуру фокальной плоскости ФПУ КА во время проведения съемки;

- иные параметры, зависящие от конструктивных особенностей ФПУ КА.

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

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

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

T - текущая температура фокальной плоскости;

5.9 Скорректированные значения цифровых отсчетов считаются ненадежными:

- для сбойных пикселей, возникших при передаче данных ДЗЗ из космоса с борта КА на наземную инфраструктуру приема и обработки данных ДЗЗ из космоса;

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

- пикселей, значения цифровых отсчетов которых выходят за пределы динамического диапазона АЦП на борту КА;

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

5.10 В результате проведения относительной радиометрической коррекции должен быть получен набор данных ДЗЗ из космоса, содержащий:

- данные ДЗЗ из космоса для каждого спектрального канала, значения пикселей которых представлены исправленными цифровыми отсчетами;

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

6 Требования к алгоритму абсолютной радиометрической коррекции

6.1 Абсолютная радиометрическая коррекция проводится с целью преобразования значений цифровых отсчетов пикселей данных ДЗЗ из космоса, прошедших относительную радиометрическую коррекцию, к эффективным значениям СПЭЯ либо значениям КСПЭЯ, отнесенных к уровню верхней границы атмосферы.

Примечание - Справочный солнечный спектр Мирового центра радиационных данных ВМО в Давосе рекомендован к использованию при обработке данных ДЗЗ из космоса подгруппой по оптическим датчикам видимого и инфракрасного диапазонов группы по калибровке и валидации CEOS.

6.4 Каждый раз перед преобразованием эффективных значений СПЭЯ пикселей, отнесенных к верхней границе атмосферы, к значениям КСПЭЯ пикселей, отнесенных к тому же уровню, должны быть рассчитаны:

- зенитные углы Солнца для каждого пикселя данных;

- расстояние между Землей и Солнцем.

6.5 Алгоритм абсолютной радиометрической коррекции при расчете зенитных углов Солнца должен аппроксимировать поверхность Земли эллиптической поверхностью, поднятой на среднюю фиксированную высоту местности над общеземным эллипсоидом. В качестве общеземного эллипсоида рекомендуется использовать эллипсоид GRS80, одобренный Международным геодезическим и геофизическим союзом [1].

где d - расстояние между Землей и Солнцем, выраженное в астрономических единицах;

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

- данные ДЗЗ из космоса для каждого спектрального канала, значения пикселей которых представлены эффективными значениями СПЭЯ либо значениями КСПЭЯ;

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

7 Требования к алгоритму атмосферной коррекции

7.1 Атмосферную коррекцию проводят с целью преобразования эффективных значений СПЭЯ либо значений КСПЭЯ, отнесенных к уровню верхней границы атмосферы, к значениям тех же величин, отнесенных к уровню поверхности Земли.

7.2 Перед проведением атмосферной коррекции должна быть проведена классификация данных ДЗЗ из космоса в пределах обрабатываемой сцены с созданием маски облачности и маски облачных теней.

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

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

- решение уравнения атмосферной коррекции для каждого пикселя данных.

7.4 Требования к алгоритму составления уравнения атмосферной коррекции

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

(9)

Примечание - Физический смысл слагаемых, входящих в уравнение (7), заключается в том, что первое слагаемое отвечает за солнечное излучение, отраженное внутри атмосферы в направлении КА. Второе слагаемое главным образом связано с солнечным излучением, напрямую отраженным от изучаемого элемента поверхности в направлении КА. Третье слагаемое вызвано солнечным излучением, отраженным в направлении КА элементами поверхности, соседними с изучаемым элементом.

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

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

7.4.3 Алгоритм составления уравнения атмосферной коррекции должен извлекать объемные спектральные коэффициенты поглощения из международных баз данных о спектральных свойствах газов.

Примечание - Наиболее известной и используемой базой данных о спектральных свойствах газов является HITRAN.

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

Примечание - Наиболее известными операторами выходных данных реанализов являются Европейский центр среднесрочного прогнозирования (ECMWF) и Национальное управление океанических и атмосферных исследований (NOAA).

7.4.5 Алгоритм составления уравнения атмосферной коррекции должен вычислять коэффициент внутреннего отражения атмосферой солнечного излучения в направлении КА за счет молекулярного (рэлеевского) рассеяния через молекулярную оптическую толщину, которая зависит от высотного профиля концентрации рассеивателей (молекул газов), а также абсолютной высоты изучаемого элемента поверхности. Угол рассеяния следует рассчитывать по зенитным углам и разности азимутов Солнца и КА. При этом следует придерживаться требований, указанных в 6.5.

7.4.6 Алгоритм составления уравнения атмосферной коррекции должен вычислять высотный профиль концентрации молекул атмосферных газов с помощью климатических моделей высотных профилей давления и температуры воздуха. Рекомендуется использовать климатическую модель глобальной справочной атмосферы по ГОСТ Р 53460.

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

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

Пример - Модель смеси аэрозолей городского ландшафта предполагает 17% пылевидных аэрозолей, 61% водорастворимых аэрозолей и 22% сажи.

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

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

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

Примечание - Продукты обработки данных ДЗЗ из космоса, содержащие оптическую толщину аэрозолей, создаются на основе применения методики Кауфмана или методики "Глубокого синего". Методика Кауфмана оценивает оптическую толщину аэрозолей над темными объектами (влажная растительность, пашня и др.) исходя из корреляции значений СПЭЯ этих объектов в канале прозрачности среднего инфракрасного диапазона (как правило, 2,1 мкм) со значениями СПЭЯ этих объектов в красном и синем каналах видимого диапазона. Методика "Глубокого синего" позволяет оценить оптическую толщину аэрозолей над яркими поверхностями в результате анализа наблюдений в каналах сине-фиолетовой области видимого диапазона (как правило, 0,412 или 0,470 мкм), для которых имеет место сильное аэрозольное рассеяние, но слабое отражение от земной поверхности.

7.4.11 Алгоритм составления уравнения атмосферной коррекции должен вычислять спектральный коэффициент направленного пропускания и спектральный коэффициент диффузного пропускания через молекулярную и аэрозольную оптические толщины, а также зенитные углы КА или Солнца.

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

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

- коэффициент внутреннего отражения атмосферы солнечного излучения в направлении КА за счет молекулярного (рэлеевского) и аэрозольного рассеяния;

- спектральный коэффициент направленного пропускания;

- спектральный коэффициент диффузного пропускания;

- сферическое альбедо атмосферы.

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

Таблица 1 - Рекомендуемые максимальный шаг сетки и диапазон значений параметров условий съемки в файлах справочных таблиц

7.4.15 Алгоритм составления уравнения атмосферной коррекции должен обеспечивать интерполяцию элементов уравнения атмосферной коррекции от значений параметров условий съемки, заданных в узлах сетки справочных таблиц, к текущим значениям параметров условий съемки, полученных из внешних источников согласно рекомендациям, указанным в 7.4.2-7.4.10.

7.5 Требования к алгоритму решения уравнения атмосферной коррекции

7.5.1 Алгоритм решения уравнения атмосферной коррекции, составленного для каждого пикселя данных ДЗЗ из космоса, должен состоять из следующих шагов:

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

7.5.3 Значения КСПЭЯ, отнесенные к уровню поверхности Земли, для пикселей, попавших в маски облачности и маски облачных теней, а также для пикселей, для которых оптическая толщина аэрозолей более 1,5 или зенитный угол Солнца более 70°, считаются ненадежными.

8 Требования к алгоритму коррекции влияния анизотропии земной поверхности

8.1 Коррекция влияния анизотропии земной поверхности проводится с целью исключить из временных рядов значений КСПЭЯ, отнесенных к уровню поверхности Земли, вариации, вызванные изменениями геометрических условий облучения этой поверхности Солнцем и наблюдения этой поверхности ФПУ КА.

8.2 Коррекцию влияния анизотропии земной поверхности следует проводить для данных ДЗЗ из космоса низкого и сверхнизкого пространственного разрешения.

8.3 Для проведения коррекции влияния анизотропии земной поверхности значения КСПЭЯ каждого пикселя, отнесенные к уровню поверхности Земли, следует представлять в виде:

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

8.4 Алгоритм проведения коррекции влияния анизотропии земной поверхности должен вычислять численные коэффициенты по формуле (12) для каждого спектрального канала и каждого пикселя данных ДЗЗ из космоса на основании набора значений КСПЭЯ, полученных для одной и той же поверхности, но при разных углах облучения и наблюдения. Рекомендуемый период сбора данных составляет 30 дней.

8.5 По найденным значениям численных коэффициентов следует вычислить по формуле (12) для каждого пикселя и каждого спектрального канала среднее для периода сбора данных значение КСПЭЯ, отнесенное к уровню поверхности Земли и к заданным геометрическим условиям облучения и наблюдения.

8.6 Геометрические условия облучения и наблюдения, к которым приводятся значения КСПЭЯ в 8.5, должны оставаться фиксированными для всей коллекции данных ДЗЗ из космоса, прошедших коррекцию влияния анизотропии земной поверхности.

Приложение А

(обязательное)

Справочный солнечный спектр Мирового центра радиационных данных Всемирной метеорологической организации в Давосе

Таблица А.1

Продолжение таблицы А.1

Продолжение таблицы А.1

Продолжение таблицы А.1

Продолжение таблицы А.1

Окончание таблицы А.1

Библиография

[1] Resolution N 7 of the XVII General Assembly of the IUGG in Canberra, December 1979