автореферат диссертации по авиационной и ракетно-космической технике, 05.07.12, диссертация на тему:Радиометрическая и радиационная коррекция информации, получаемой оптико-электронной камерой обзора Земли космического аппарата 17Ф118-О

НАУЧНО-ИССЛЕДОВАТЕЛЬСКИЙ ИНСТИТУТ МИКРОПРИБОРОВ НПАО"ЭЛАС"

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

УДК 551.51:629.78

РАДИОМЕТРИЧЕСКАЯ И РАДИАЦИОННАЯ КОРРЕКЦИЯ ИНФОРМАЦИИ, ПОЛУЧАЕМОЙ ОПТИКО- ЭЛЕКТРОННОЙ КАМЕРОЙ ОБЗОРА ЗЕМЛИ КОСМИЧЕСКОГО АППАРАТА 17Ф118-0.

05.07.12 - дистанционные аэрокосмические исследования

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

Москва 1996

Работа выполнена в Московском физико-техническом институте

Научный руководитель - доктор физико- математических наук,

профессор Кондранин Тимофей Владимирович

Официальные оппоненты - доктор физико-математических наук Козодеров Владимир Васильевич кандидат технических наук Широков Василий Васильевич

Ведущая организация - НПО им. С.А. Лавочкина г. Химки, Московской области

Защита состоится "Л-1" иЛОНСл 1996 г. в -10*" часов на заседании диссертационного совета Д142.06.02 при НИИ микроприборов по адресу: г. Москва 103460, НИИ МП.

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

Автореферат разослан " М "¿реЛзЬал* 1996 г.

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

/Широ Георгий Эдуардо]

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

Актуальность темы

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

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

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

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

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

Актуальность диссертационной работы связана с необходимостью разработки методического и программно алгоритмического обеспечения решения задачи по восстановлению спектральных отражающих свойств подстилающей поверхности Земли с использованием информации камеры ОЭК-ОЗ разработки НПЦ "ОПГЭКС", установленной на КА 17Ф118-0.

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

коэффициента отражения подстилающей поверхности (КО ПП) в рамках КН2 "ОБЗОР".

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

1) Научно методическое обоснование и разработка адекватного метод; восстановления КО ПП по сигналу, регистрируемому приемником на BepxHei границе атмосферы.

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

3) Расчетно-параметрическое исследование влияния различных факторо! системы "Земля - атмосфера - приемная аппаратура" (спектрального диапазона геометрии визирования, условий освещения, состояния атмосферы, отражающш характеристик подстилающей поверхности, калибровочных коэффициентов функций чувствительности спектральных каналов приемника) на точносп восстановлен™ КО по сигналу на входе приемной системы. При этом особеннс важно исследовать погрешности, возникающие при использовании CTatwapraoi атмосферы, и найти условия, когда погрешность минимальна, и способы сниженш погрешности, оптимизируя условия зондирования.

Научная новизна исследования заключается в том, что:

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

2) В результате численного моделирования установлены оптимальные условн* проведения съемки ПП оптико- электронными системами, которые обеспечиваю! минимальную погрешность восстановления КО ПП.

3) С использованием созданных программных модулей впервые проведена обработка, анализ и интерпретация данных международного подспутниковогс эксперимента FIFE (First International Field Experiment, Канзас, 1987г.).

4) С помощью разработанных методик и комплексов программ, включенных в состав штатного ПАО в/ч 32103, проведена полетная калибровка камеры ОЭК-ОЗ, что позволило получить достоверные данные по коэффициентам различных типов ПП.

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

1) всесторонним методическим обоснованием выбора исходной оптической модели атмосферы и отражающих свойств подстилающей поверхности;

2) физическим обоснованием и сведением к минимуму упрощающих предположений в модели, используемой при расчете переноса излучения в системе "подстилающая поверхность - атмосфера";

3) методическим исследованием внутренних свойств решения краевой задачи для уравнения переноса лучистой энергии (УПЛЭ), тестированием программ, а также сравнением с точными аналитическими и численными расчетами других авторов;

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

Практическая ценность работы состоит в том, что разработанные программно-реализованные алгоритмы радиационной и радиометрической коррекции включены, как составная часть, в комплекс первичной обработки спектрозональной видеоинформации СЗВИ фоно-целевой обстановки (ФЦО) на поверхности Земли, получаемой камерой ОЭК-ОЗ КА 17ФП8-0 (ПАО ЦТЕА.00011-01) в рамках КНЭ "ОБЗОР". Комплекс входит в состав ПАО центра обработки информации эксплуатирующей организации (в/ч 32103).

Разработанные методики, алгоритмы и пакеты прикладных программ могут быть использованы в организациях, занимающихся проблемами:

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

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

- разработки перспективных оптико-электронных систем наблюдения (обеспечение необходимой точности калибровки, выбор спектральных каналов и т.д.).

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

1) Методика восстановления спектрального альбедо ПП по значениям яркости восходящего излучения на верхней границе атмосферы. Методика позволяет достигать приемлемой точности восстановления спектрального альбедо при минимальных затратах машинного времени.

2) Результаты обработки, анализа и интерпретации данных международного подспутникового эксперимента FIFE (Канзас, 1987г.),.которые позволили провести верификацию всех разработанных алгоритмов и программ.

3) Разработанное ПАО радиометрической, атмосферной коррекций и полетной калибровки, предназначенное для обработки С31?И, получаемой аппаратурой ОЭК-ОЗ КА 17Ф118-0. Разработанный комплекс программ входит в состав ПАО ! (ОИ АНИУК-Э эксплуатирующей организации.

4) Результаты полетной калибровки камеры ОЭК-ОЗ КА 17Ф118-0 и значения альбедо, восстановленные по данным СЗВИ, получаемой этой камерой.

Работа выполнялась на кафедре "Системы, устройства и методы геокосмической физики" (СУМГФ) Московского физико-технического института (МФТИ) в период 1991-1995 гг. Представленные в работе исследования проводились в соответствии с планами работ кафедры СУМГФ МФТИ и ОКР "СТРОЙ-О" (КНЭ "ОБЗОР") выполняемой в НПЦ "ОПТЭКС" НПО "ЭЛАС".

Апробация работы. Результаты работы докладывались на 15 Международной конференции по лазерному зондированию в г. Томске (1990г.) на Всероссийской конференции в ВИА им. Можайского (1995 г.), ежегодных научных, конференциях МФТИ.

Объем работы

Диссертационная работа изложена на 150 страницах машинописного текста, состоит из: Введения, 3-х глав, Заключения и Приложения; включает 32 рисунка и список основной испрльзованной литературы из 148 наименований отечественных и зарубежных авторов на 15 страницах.

СОДЕРЖАНИЕ ДИССЕРТАЦИИ

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

В главе 1 излагается метод расчета альбедо ГШ по данным спектрозональной информации в видимом и ближнем ИК диапазоне, учитывающий особенности оптических свойств системы "атмосфера- подстилающая поверхность". Дается обоснование основных упрощающих предположений и оценивается точность вычисления ОПФА для нескольких методов.

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

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

<Рд падает параллельный пучок солнечных лучей, создающих на верхней границе атмосферы спектральную освещенность площадки, перпендикулярной солнечном) излучению, равной Индекс длины волны X в дальнейшем будем опускать. Е

среде, ограниченной снизу горизонтально неоднородным Ламбертовым дном I альбедо я(х,у), происходят процессы многократного рассеяния и истинногс поглощения с произвольной нормированной индикатрисой рассеяния Р("у,г) показателем ослабления а(г) и рассеяния а(г), вероятностью выживания квантг (альбедо однократного светорассеяния) зависящих от высоты ъ

Направление вышедшего из атмосферы света будем определять углом агссо5г|=6 между рассеянным излучением и нормалью к слою, а также азимутом ф отсчитываемым от плоскости солнечного вертикала.

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

1) эффектами рефракции и поляризации при рассеянии солнечного излучения е атмосфере можно пренебречь;

2) о тражение от поверхности происходит по закону Ламберта;

3) эффекты, связанные со сферичностью Земли, не проявляются;

4) ослаблением контрастов яркостей горизонтально- неоднородной подстилающей поверхности в зависимости от пространственной частоты можно пренебречь.

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

¿ = + — (яе'^+дАо) (1)

1 п

где О - яркость атмосферной дымки; пЕ - освещенность земной поверхности усредненная по горизонтальным координатам при q=0; Е - освещенность земной

поверхности, усредненная по горизонтальным координатам; А^ - диффyзнa^ составляющая коэффициента пропускания; с„ - сферическое альбедо атмосферы.

На практике, вместо яркости излучения удобно использовать эквивалентное шанетарное альбедо q (альбедо системы "атмосфера-подстилающая поверхность"). 1т о связь сданными измерений яркости в i-м канале L. задается уравнением

а* =—--. =а0,- +al.--q + ali-q, (2)

cos(©0) Е0(Л*)

где cl - расстояние между Землей и Солнцем, выраженное в астрономических единицах; i- номер канала; Л* - эффективная длина волны в i-м канале;

(з)

cos(©o) £-0(А,) d1 ЕЛХ)У

a\j = —--------, (4)

cos(©0) Eq (Л,- )

mi1 E (A*)А0(Л*)

al: =------2---. (5)

cos(0o) £0(Д.)

Коэффициенты aO, a! и a2 зависят от альбедо ПП (через £ ), однако для

реальных оптических толщин атмосферы т=0.1-0.6 сферическое альбедо принимает значения с^ =0.05-0.17, поэтому можно рассчитывать Е при некотором среднем

альбедо (например, q=0.3). Простые оценки показывают, что такая аппроксимация приводит к погрешности, не превышающей 2-3% при q=0.0-0.6. Среднее альбедо вычисляется по формуле

J Î+I /+(

=7?7+пт s (б)

(2/ + 1) k^i-t l=l-i

Как показало сравнение с точными расчетами, такое предположение приводит к хорошим результатам, если правильно выбран размер "окна" I, который зависит от размера пикселя, параметров атмосферы в момент съемки, спектрального интервала, пространственных частот самого изображения, и, по результатам различных авторов, находится в пределах от 100 до 500 метров.

При известных L, aO, al, а2 уравнения (1) и (6) позволяют однозначно вычислять альбедо подстилающей поверхности (ПП). Решение проводилось итерационным методом:

(7)

. i '+' 7+1 . q* -а„ +а7 -<7?"'

II]arn>1 ^ =7^7,7 I ' °а J - W

(2/ + 1) k*i-n=j-t ai

Проведенный анализ показал возможным использовать п=2 для аппаратуры высокого и среднего пространственного разрешения (г<500м) и ограничиться учетом вклада атмосферной дымки (n=l, t=0) для аппаратуры низкого пространственного разрешения (г>500м).

Так как в (7), (8) входят только простейшие операции, время обработки изображения определяется скоростью получения коэффициентов aO, al, а2.

В части 2 первой главы описываются различные методы расчета оптически? передаточных функций атмосферы (и, следовательно, зависящих от них аО, а1 и а2; и их сравнение с данными точных расчетов.

В результате анализа литературных источников были выбраны и программне реализованы следующие алгоритмы расчета оптических передаточных функций атмосферы (и следовательно зависящих от них аО, а1 и а2):

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

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

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

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

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

В части 3 главы I приведено численное сравнение реализованных алгоритмов с данными тестовых расчетов, выполненных для реальной атмосферы, взятых из опубликованных данных, которое показало, что большинство из используемых ни практике методов расчетов дают погрешность расчетов яркости более 50-100% и освещенностей более 10-20%. В результате был предложен метод расчета ОПФА. сочетающий высокую скорость (сравнимую со скоростью ввода исходных данных с программу) с приемлемой точностью (1-3% для освещенности и 5-15% при вычислении яркости атмосферной дымки). Моделирование процесса расчета КС НП по данным измеренной яркости показало, что погрешность восстановления КС 1111, вызванная методическими погрешностями расчета оптических передаточных функций атмосферы, будет находиться в большинстве случаев в пределах 10-30%..

Часть 4 главы 1 посвящена анализу чувствительности расчетной схемы алгоритма восстановления альбедо ПП к вариациям входных параметров.

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

Одним из основных источников погрешности восстановления КО ПП является неточное знание типа аэрозоля в момент съемки. Если тип аэрозоля известен, то погрешность практически во всех важных случаях не превышает 10%. возрастая при больших зенитных углах Солнца 0^>5ОО, малых КО ПП я< 10-30% н

азимутальных углах ф=0-50° . Ошибка возрастает, если тип аэрозоля не определен, так погрешность, вызванная использованием при расчетах оптической модели континентального аэрозоля вместо городского, не опускается ниже 10%, достигая 20 - 40% и более. Ошибка несколько ниже (в 2- 3 раза), если континентальный аэрозоль используется вместо морского. Таким образом, если известен тип аэрозоля (например, при проведении атмосферной коррекции далеко от крупных промышленных городов, в глубине материка можно с достаточной точностью использовать континентальный аэрозоль), при проведении атмосферной коррекции

юзможно использование стандартного набора вектора оптических параметров, :оответствук>щего метеорологической дальности видимости 20 км и высоте орозольного слоя 2 км. Получаемая при этом погрешность восстановления КО ПП 1ЛЯ достаточно ярких поверхностей (q>0.1) будет в пределах 10-20%, что достаточно i большинстве прикладных задач.

Наблюдается сложная зависимость погрешности восстановления КО ПП от зенитного угла солнца. До 0^=60" погрешность может: плавно возрастать;

монотонно убывать, достигая минимума в районе 20-40° с последующим монотонным возрастанием Или иметь еще более сложную зависимость. После 60° погрешность во всех случаях монотонно возрастает.

Максимальная погрешность восстановления КО ПП монотонно убывает с увеличением КО ПП в тех случаях, когда имеется неопределенность в типе аэрозоля. В случае определенного состава аэрозоля погрешность достигает минимума, при этом положение минимума может перемещаться от q=0.15 до q=0.55 в зависимости от геометрии визирования.

Отмечена сильная зависимость погрешности восстановления КО ПП от азимута визирования, отсчитываемого относительно плоскости солнечного вертикала (ф). Зависимость уменьшается, во-первых при увеличении отражательной способности зондируемой поверхности q, во- вторых, при стремлении зенитных углов к нулю. В некоторых случаях изменение азимута от 0° до 90° или 180° , при одинаковых остальных условиях приводит к возрастанию погрешности восстановления КО ПП от 2 - 3 до 8 раз. Анализ результатов, полученных в диссертации, показал, что для зондирования ПП оптимальными являются азимуты ©0=40-150°. При таких азимутах погрешность восстановления КО ПП может в

отдельных случаях быть меньше погрешности для надирных условий наблюдения.

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

дымки в сигнал, регистрируемый приемником на верхней границе атмосферы. Показано, что целесообразно иметь аппаратурные погрешности на уровне 1-5%. Если аппаратурные погрешности будут больше, то они могут стать основным источником погрешности даже при восстановлении КО ПП для стандартной атмосферы.

В главе 2 дано описание международного подспутникового эксперимента FIFE (Канзас 1987). Приведены результаты сравнения коэффициентов отражения, рассчитанных с помощью разработанной в главе 1 мегодики по данным, полученным AVHRR NOAA9.10, сданными наземных измерений (MMR).

В частях 1-2 главы 2 описываются цели и задачи международного эксперимента; описывается аппаратура, использованная в ходе эксперимента. В части 3 обсуждаются полученные результаты, которые позволяют сделать вывод о том, что разработанная модель позволяет достаточно точно восстанавливать КО ПП даже при низких углах солнца над горизонтом и ненадирных условиях съемки.

Кроме того сопоставление позволило выявить следующие закономерности. 1) Хотя зависимость коэффициента корреляции (г) от зенитного угла Солнца (0Ä} (после

которого производилась отбраковка спутниковых измерений) несколько различается для данных с NOAA-9 и NOAA-IO, можно отметить

практически монотонное убывание г с увеличением &*,акс. Кроме того отмечен-неравномерное падение коэффициента корреляции с увеличением плавно

уменьшение до некоторого критического значения (43° - для АУНЯЯ ЫОАА-9 и 60 - для АУНЯЯ ЫОАА-Ю); последующее резкое падение в пределах 5° изменени фмакс и медленное дальнейшее уменьшение, по крайней мере, до углов 65-70

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

Солнца ©^й^с=70° коэффициент корреляции оставался достаточно высоки!

(г=0.81).

2) Отмечена более высокая корреляция данных АУНЯЯ ЫОАА-Ю с данным] наземных измерений (г^оллю-0-^ тКО/Ы9-°-81)- Д™ данных с ЫОАА-Ю такж отмечено меньшее расхождение между альбедо по абсолютным значениям. Так коэффициенты линейной регрессии {Ццш< =а+Ь* Чыолл > гДе альбедо определяется процентах) находятся в следующих пределах --1,88ч-0.81; ат/1/ПВ--4.2^-2.2

Ъыолм -23: Ьноллю =0-91 +0.97.

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

В главе 3 излагается научно- методическое и программно- алгоритмическо обеспечение решения задач по обработке спектральной видеоинформации в рамка: КНЭ "ОБЗОР", получаемой камерой ОЭК-ОЗ КА, функционирующей с августа 199 года.

Для повышения эффективности использования СЗИ, получаемо! аппаратурой ОЭК-ОЗ, в том числе и для обеспечения восстановления абсолютны: значений коэффициента отражения (КО) природных образований, наряду предполетной метрологической аттестацией осуществляется текущий полетиьп контроль энергетических и частотно-контрастных характеристик ОЭК-ОЗ. С это1 целью разработаны, научно и методически обоснованы с учетом специфически: особенностей аппаратуры алгоритмы радиометрической (РК) и радиационно! (атмосферной) коррекции (АК), а также полетной калибровки (ПК) ОЭК-ОЗ Алгоритмы реализованы в виде соответствующих программных модулей, входящи: в состав штатного математического обеспечения первичной обработки СЗИ н; АППИ.

В части 1 главы 3 описаны задачи, решаемые в рамках КНЭ "ОБЗОР", I показана работа автора, результатом которой было обоснование, разработка I реализация программ радиометрической, атмосферной коррекций и методик1 полетной калибровки.

В части 2 главы 3 описана методика радиометрической коррекцш изображения.

Как в большинстве оптико- электронных систем, в ОЭК-ОЗ наблюдаете! линейная зависимость сигнала на выходе АЦП (О) от яркости излучения на вход' приемного устройства, выраженной в стандартных энергетических единицах I [Вт/(м2-стсрмкм)]:

Ц (х,}>) = С11(х)-П(х,у) + С01(х), (9)

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

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

В соответствии с методикой разработана программа, на которую оформлена документация в соответствии с ЕСПД, и которой присвоен регистрационный номер: ХА.00043-01 34 05 ЛУ. Программа устраняет дефекты ("полосатость" изображения), обусловленные разбросом чувствительности элементов линейки ФПУ, а также рассчитывает коэффициенты перевода цифровых отсчетов скорректированного изображения в величины альбедо подстилающей поверхности.

Часть 3 главы 3 посвящена методике атмосферной коррекции изображения, используемой при обработке информации, получаемой камерой ОЭК-ОЗ. Методика описана в главе 1. На основании методики создана программа, на которую оформлена документация в соответствии с ЕСПД, и которой присвоен регистрационный номер: ХА.00043-01 34 08 ЛУ. Программа устраняет влияние атмосферы путем учета коэффициента пропускания атмосферы, вычитания засветки атмосферной дымки, устранения влияния боковых подсветов от соседних элементов.

Часть 4 главы 3 посвящена научно- методическому и программно-алгоритмическому обеспечению решения задачи полетной калибровки аппаратуры ОЭК-ОЗ КА 17Ф118-0.

Процедура полетной калибровки аппаратуры заключается в получении калибровочных коэффициентов, необходимых для проведения радиометрической коррекции получаемого изображения. Так как на камере ОЭК-ОЗ КА 17Ф118-0 отсутствует аппаратура для проведения калибровки, единственная возможность получения калибровочной информации заключается в использовании расчетов яркости на входе приемного устройства над поверхностями с известными отражающими свойствами и построении на основании этих расчетов калибровочной кривой.

Решение этой задачи проводится в три этапа:

1) Находятся коэффициенты для относительной калибровки: коэффициенты СО и С1, устраняющие разброс фоточувствительностей между различными детекторами.

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

3) Вычисляются с помощью полученных величин калибровочные коэффициенты СО и С1 - коэффициенты перевода сигнала в разрядах АЦП. в яркость излучения на уровне приемника излучения.

Для нахождения темновых токов было предлржено использовать результаты ночных съемок. Для анализа возможности такого пути нахождения темновых токов был исследован вопрос о влиянии отражающих свойств ППЗ,.наличия облачности и влияния луны на величину яркости выходящего из атмосферы, излучения в сумеречное и ночное время. Исследования проводились путем математического моделирования с помощью пакета программ LOWTRAN-7. Результаты моделирования позволили сделать вывод о том, что при зенитных углах Солнца больше 93-95° ночные снимки .могут использоваться в качестве черного тела.

Для.: подтверждения • возможности использования ночной съемки для определения темновых токов была проанализирована зависимость сигнала на выходе аппаратуры AVHRR NOAA-9 и NOAA-IO от зенитного угла Солнца для вечерних и ночных снимков, полученных в рамках FIFE (см. рис. 1). Такая проверка важна, потому что в расчетах яркости при низких положениях Солнца над горизонтом сложно учесть все факторы, начинающие играть важную роль в формировании излучения, такие как сферичность земли, эффекты рефракции, учет горизонтальной неоднородности атмосферы и др..

Анализ полученных зависимостей показал, что при зенитных углах солнца меньше 90° сигнал на выходе приемника AVHRR NOAA-IO не превышал 42 ед. АЦП, что хорошо согласуется с данными, полученными по предполетной калибровке и сигналом при направлении камеры в открытый космос: 38-40 ед. АЦП. Однако при углах более 100° AVHRR . NOAA-IO начинает показывать заниженные значения сигнала (D<20-30). Кроме того отмечено противофазное изменение яркостей в подспутниковой точке и при визировании в открытый космос. Полученная закономерность по видимому объясняется нелинейными аппаратурными эффектами, возникающими при прохождении терминатора: Таким образом необходимо с осторожностью использовать данные ночных съемок пp^ калибровке аппаратуры, и проводить исследования диапазона зенитных углои солнца, при которых не возникает указанных эффектов, т.е. оценить максимальный зенитный угол Солнца, при котором возможно производить калибровку аппаратуры. В частности для AVHRR NOAA-IO углы Солнца, при которьо возможно использовать ночные съемки для калибровки, находятся в диапазоне 90 100°: '

Для получения темновых токов камеры была произведена серия ночны; съемок пустынь Руб-эль-Хали и Сахары. Анализ гистограмм полученньп изображений показал, что средний сигнал во время ночных снимков был на уровш 0-3 ед. АЦП. На основании данных предполетной калибровки ожидалось, что это: сигнал будет на уровне десятков едениц АЦП. Полученные значения показали, чт< величина темновых токов существенно изменилась в процессе эксплуатации ОЭ1С ОЗ. Так как полученный по данным ночных снимков сигнал находится на уровт шумов, то использование ночных съемок не позволяет сделать вывод о величин! темновых токов (возможно С0<0).

Таким образом в случае ОЭК-ОЗ КА 17Ф118-0 использование ночных съемо! не позволяет извлечь информацию о коэффициенте СО.

Анализ снимков, полученных при использовании камеры ОЭК-ОЗ, позвол1и сделать следующие выводы:

80 90 100 110 120

зенитный угол Солнца (градусы)

не. 1. Сигнал (в ед. АЦП) полученный аппаратурой AVHRR NOAA-IO полученный з время эксперимента FIFE (Канзас, 1987 г.) вечернее и ночное время в зависимости I' зенитного угла Солнца (1,2) и сигнал полученный при направлении приемник в гкрытый космос (3,4): 1,3- канал 1; 2,3 - канал 2.

ис.2. Фрагмент изображения от 20.02.1995, включающий полигон White Sands, Нью [ексико (виток 2769), полученный с помощью камеры ОЭК-ОЗ МКА-ЭКО зениты й угол Солнца 0=56.4 , 1-й канал.

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

2) Изучение сигнала над однородными (квази-однородными) поверхностно показал плавные изменения сигнала вдоль линии сканирования.

Для удаления указанных эффектов был использован метод относительно! калибровки приемников по однородным поверхностям, в качестве которых былс принято решение использовать снежные поверхности Антарктиды и пустынь. Былг произведена серця съемок камерой ОЭК-ОЗ КА 17Ф118-0 для обеспечения решена задачи относительной калибровки. После визуального анализа изображений был! отобраны изображения, имеющие фрагменты с однородными поверхностями пустыня Рубь-Эль-Хали, Антарктида.

Для исследования степени однородности приведенных изображений бьип вычислень! средние яркости и стандартные отклонения яркости по столбцам i строкам этих изображений. Полученные результаты показывают сильное различи! стандартных отклонений яркости, вычисленных для строк и столбцов исходны; изображений (от 3-4 до 60 раз в зависимости от изображения и номер: спектрального канала). Анализ стандартных отклонений, усредненных по строка сигналов, рассчитанных по исходному изображению, подтверждают предположена об однородности исходных фрагментов изображений. Так, средняя яркость п< строкам изображений тестовых полигонов меняется в пределах 0.J-3.1% О зависимости от изображения и номера канала). Таким образом величин! стандартного отклонения средней яркости в строках изображения от средней п< снимку величины может являться хорошей количественной характеристике! степени однородности снимаемого объекта. А так как отсутствуют выделенные направления на исходных изображениях, этой величиной можно оценить i горизонтальную неоднородность изображения.

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

После анализа литературных источников было принято решени использовать следующие типы поверхностей для абсолютной калибровк: аппаратуры ОЭК-ОЗ:

- гипсовый песок - полигон White Sands;-Нью-Мексико (см. рис.2);

- темные объекты на всех имеющихся в наличии изображениях.

Для расчетов яркости на верхней границе атмосферы был использован паке LOWTRAN-7, потому что в отличии от .метода, приведенного в главе I, oi позволяет учитывать высоту полигона'над уровнем моря, оптические свойств Пустынного аэрозоля и моделировать влияние разных типов облачности н регистрируемый сигнал. Результаты' расчетов яркости для условий в момен зондирования и соответствующие им^игналы на выходе АЦП приведены на рис. 2 Чтобы оценить погрешность, которые возникают из-за неточного знания состояни атмосферы в момент съемки, на рисунке приведен разброс яркости, рассчитаны; для различной прозрачности атмосферы для конкретного полигона.

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

а) Канал 1 ед. АЦП 6) Канал 2 ед. АЦП

200 180 160 1 140 120 100 80 60 40

Тт1 1 ✓

1

I ■ 1 ■ • 1

100

200

в) Канал 3

300 ед. АЦП

400

500

600

300

г) Канал

Рнс.З. Калибровочные характеристики ОЭК-ОЗ КА 17Ф118-0: а) - 1-й, б) - 2-й, в) - 3-й, г) - 4-й канал.

зависимости между сигналами в ед. АЦП (О) и яркостью на верхней границ! атмосферы, рассчитанной при метеорологической дальности видимости 23 км Соответствующие кривые приведены на рис. 3.

Анализ полученных результатов позволяет сделать следующие выводы:

1) Для ярких объектов (типа гипсовых песков с коэффициентами отражения 0.6-0.7 погрешность определения яркости будет на уровне 3-4%. Учет всех источнике] погрешности (таких как неоднородность полигона, методические погрешности неучет неортотропности) дает погрешность не более 8-15%.

2) Точность определения яркости темных объектов будет несколько ниже: 10- 50%,; в некоторых случаях до 100%.

Высокая погрешность определения яркости темных объектов связана неопределенностью оптического состояния атмосферы, поэтому приходилоа учитывать возможность метеорологической дальности видимости У=4 км Численные расчеты показали, что если бы метеорологическая дальность видимосп лежала в пределах 10-50 км, то точность расчета яркости была бы на уровне 10%.

С другой стороны анализ яркостей и сигналов, вычисленных для водно1 поверхности, показывает их тесную связь между собой: увеличение яркости пр: увеличении сигнала при уменьшении зенитного угла солнца. Последнее возможно случае, когда состояние атмосферы в момент получения указанных четыре; изображений было похожим. Так как съемка производилась в разных местах ; разное время при полном отсутствии или минимальной облачности и изображени: отличаются высокой контрастностью, то вероятность того, что во всех четыре; случаях низкая прозрачность атмосферы, невысока. Поэтому точность определени яркости темных объектов по спутниковым измерениям будет на уровне 10-15%.

Показаны пути увеличения точности калибровки. Разработанный комплек программ в виде составной части комплекса первичной обработки изображена поверхности земли, получаемых камерой ОЭК-ОЗ КА 17Ф118-0, передан эксплуатирующую организацию.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

Сформулируем основные результаты, полученные в работе:

I. Разработана методика восстановления спектрального альбедо ПП по значения» яркости восходящего излучения на верхней границе атмосферы. Методик позволяет достигать приемлемой точности восстановления спектрального альбед! при минимальных затратах машинного времени. При этом для расчетов ОПФА бы. предложен метод расчета ОПФА,-сочетающий высокую скорость (сравнимую а скоростью ввода исходных данных в программу) с приемлемой точностью (1-3% дд освещенности и 5-15% при вычислении яркости атмосферной дымки).

II. Проведены расчеты предложенным методом поля рассеянного излучения дл системы: атмосфера- подстилающая поверхность, на основе которых установлен! закономерности распределений погрешности восстановления альбедо ПП п< излучению, попадающему в апертуру приемника, и дано их объяснение. Сред] полученных закономерностей можно отметить следующие.

1) Одним из основных источников погрешности восстановления КО ПГ является неточное знание типа аэрозоля в момент съемки. Если тип аэрозол известен, то погрешность практически во всех важных случаях не превышает 107!

возрастая при больших зенитных углах Солнца ©0>50", малых КО ПП ц< 10-30% |

азимутальных углах ф=0-50°. Ошибка возрастает, если тип аэрозоля не определен

Гак, погрешность, вызванная использованием при расчетах оптической модели сонтинентального аэрозоля вместо городского, не опускается ниже 10%, достигая 20 - 40% и более. Ошибка несколько ниже (в 2- 3 раза), если континентальный ирозоль используется вместо морского. Таким образом, если известен тип аэрозоля например, при проведении атмосферной коррекции далеко от крупных лромышленных городов, в глубине материка можно с достаточной точностью «пользовать континентальный аэрозоль), при проведении атмосферной коррекции юзможно использование стандартного набора вектора оптических параметров, ;оответствуюiцего метеорологической дальности видимости 20 км и высоте 1эрозольного слоя 2 км. Получаемая при этом погрешность восстановления КО ПП 1ля достаточно ярких поверхностей (q>0.1) будет в пределах 10-20%, что достаточно s большинстве прикладных задач.

2) Наблюдается сложная зависимость погрешности восстановления КО ПП в

¡ависимости от зенитного угла солнца. До 0Q=6O° погрешность может: плавно

юзрастать; монотонно убывать, достигая минимума в районе 20-40° с последующим монотонным возрастанием или иметь еще более сложную зависимость. После 60° югрешность во всех случаях монотонно возрастает.

3) Максимальная погрешность восстановления КО ПП монотонно убывает с Увеличением КО ПП в тех случаях, когда имеется неопределенность в типе 1эрозоля. В случае определенного состава аэрозоля погрешность достигает минимума, при этом положение минимума может перемещаться от q=0.15 до q=0.55 ! зависимости от геометрии визирования.

4) Отмечена сильная зависимость погрешности восстановления КО ПП от 1зимута визирования ф. Зависимость уменьшается, во-первых при увеличении >тражательной способности зондируемой поверхности q, во- вторых, при премлении зенитных углов, к нулю. В некоторых случаях изменение азимута от 0° 10 90° или 180°, при одинаковых остальных условиях приводит к возрастанию югрешности восстановления КО ПП от 2 - 3 до 8 раз. Анализ результатов,

юлученных в диссертации, показал, что для зондирования ПП при больших 0()

штимальными являются азимуты ф=40-150°. При таких азимутах погрешность ¡осстановления КО ПП может в отдельных случаях быть меньше погрешности для шдирных условий наблюдения.

II. Проанализировано влияние аппаратурных погрешностей на точность ^становления КО ПП. На основании анализа фундаментальных выражений делан вывод об усилении аппаратурных погрешностей атмосферой, причем оэффициент усиления к^ возрастает при увеличении вклада яркости атмосферной

(.ымки в сигнал, регистрируемый приемником на верхней границе атмосферы. 1оказано, что целесообразно иметь аппаратурные погрешности на уровне 1-5%. •ели аппаратурные погрешности будут больше, то они могут стать основным [сточником погрешности даже при восстановлении КО Г1П для стандартной тмосферы.

V. Проведено сопоставление альбедо поверхности, рассчитанных по измерениям NOAA, и альбедо, измеренных на Земле в ходе международного эксперимента rIFE (1987). Полученные результаты позволили верифицировать все разработанные лгоритмы и программы по радиометрической и радиационной коррекции и делать вывод о том, что разработанная модель позволяет с приемлемой точностью

восстанавливать альбедо ГШ даже при низких углах солнца над горизонтом i нснадирных условиях зондирования.

V. Разработана, обоснована, реализована в виде программы методик; радиометрической коррекции изображения. На программу оформлен; документация в соответствии с ЕСПД, и ей присвоен регистрационный номер ХА.00043-01 34 05 ЛУ. Программа устраняет дефекты ("полосатость" изображения) обусловленные разбросом чувствительности элементов линейки ФПУ, а такж' переводит цифровые отсчеты исходного сигнала в яркости излучения на вход приемного устройства.

VI. Разработана, обоснована, реализована в виде программы методик; атмосферной коррекции изображения. На программу оформлена документация i соответствии с ЕСПД, и ей присвоен регистрационный номер: ХА.00043-01 34 0! ЛУ. Программа устраняет влияние атмосферы путем учета коэффициент; пропускания атмосферы, вычитания засветки атмосферной дымки, устранени: влияния боковых подсветов от соседних элементов.

VII. Разработана методика полетной калибровки аппаратуры видимого и БИ! областей спектра. Проведена относителытя и абсолютная калибровка камер; ОЭК-ОЗ на основе разработанной методики. Проведенный анализ показал, чт! точность относительной калибровки находится на уровне 1-3%, а абсолютной 8 15%. Показаны'пути увеличения точности калибровки.

Разработанный комплекс программ в виде составной части комплекс первичной обработки изображений поверхности земли, получаемых камерой ОЭК ОЗ КА 17Ф118-0, передан в эксплуатирующую организацию.

Основные результаты диссертации опубликованы в следующих работах:

1. Borin O.V.,Bushnev S.V.,Grebennik A.R.,Kondranin T.V., Myasniko E.V.,Yurin D.V.. Applied Methods for Modeling the Backscattering Characteristics in th Airborne Laser Remote Sensing of the Ocean. II 15th International Laser Rada Conference. July 23-27 1990, Tomsk. -Vol.2. -P.42-43.

2. Baranov Yu.F., S.S.Bogdanov, O.V.Borin, T.V.Kondranin, E.V. Myasnikov V.A. Filippov, I.P. Chernyavskyi - The experimental Investigation of the Lida Backscatter Signal (Black Sea, August, 1989) - 15 th Internat. Laser Radar Conf. July 23 27 1990, Tomsk. -Vol.2, p. 32-36.

3. О.В.Борин, С.В.Бушнев. Применение методов атмосферной коррекции npi обработке космической спектрозональной информации по подегилающа поверхности. Отчет НПЦ "ОПТЭКС" инв. 0Ю/"ОПТЭКС"/30.03.94г'.

4. О.В.Борин, ПО-СЗВИ-ФЦО ПРОГРАММЫ ОБРАБОТКА СПЕКТРОЗОНАЛЬНОЙ ВИДЕОИНФОРМАЦИИ ПО ФЦО, Радиометрйческа коррекция изображений по данным предполетной калибровки ФПУ. Руководств оператора. ХА.00043-01 34 05-ЛУ 1994г.

5. О.В.Борин, ПО-СЗВИ-ФЦО ПРОГРАММЫ ОБРАБОТКА СПЕКТРОЗОНАЛЬНОЙ ВИДЕОИНФОРМАЦИИ ПО ФЦО, Программ атмосферной коррекции СЗИ. Руководство оператора. ХА.00043-01 34 08 -ЛУ 199 г.

6. О.В.Борин, С.В. Бушнев, Т.В. Кондранин, Е.В.Мясников. Погрешност определения био-оптических параметров верхнего слоя океана по измерения! спектра, коэффициента яркости.// Прикладные задачи аэромехайники геокосмической.физики, Междуведомственный сборник научных трудов, Москв 1992 г.