Метод оптической обработки пространственной информации для целей наведения и посадки космических летательных аппаратов на малые тела солнечной системы

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

Метод оптической обработки пространственной информации для целей наведения и посадки космических летательных аппаратов на малые тела солнечной системы

Лобанов Александр Анатольевич

кандидат технических наук 

доцент, кафедра инструментального

и прикладного программного обеспечения,

МИРЭА - Российский технологический университет 

Аннотация

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

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

Ключевые слова: оптическая обработка, пространственная информация, посадка КЛА, наведение КЛА, малые тела, распознование образов, вычислительные комплексы, компьютерное зрение, оптический вычислитель, нейтрализация эталона

Abstract

In this article the authors focus on the issues that relate to using the promising technique of optical correlation for onboard guidance and landing a space probe on a surface of a solar system small body. The authors use the ability of optical systems to perform the Fourier transformation which allows to construct optical-electronic devices of high-speed and high-efficiency computing systems. Distribution of illumination (energy) inplane or through the lens of the entrance pupil is an input data flow for an optical computing system that has been constructed on the basis of the aforesaid principle. In a general way, this information will function as a reference standard tranmission according to two (x, y) coordinates. To analyze the aforesaid aspects, the authors of the article have used mathematical analysis methods, in particular, integral calculus for creating a mathematical model of an optical correlation bearing finder with a revolving the reference standard. The authors suggest to modulate an optical signal going through a transparent reference standard both upon the amplitude and phase. The advantage of the method is a superspeed of information processing that is restricted by the light velocity and happens almost at instance taking into account small sizes of devices. The authors give a theoretical justification of the efficiency of the optical polar correlation method. The authors demonstrate that in order to neutralize the reference standard when guiding and landing a space probe it is better to use a revolving of a reference image. The authors also analyze the method of revolving of a reference image to find the function extremum of the correlation signal and neutralizing the reference standard. 

Keywords: neutralizing the reference standard, optical computer, computer vision, computational complexes, pattern recognition, small bodies, landing a space probe, guidance a space probe, spatial information, optical processing 

Краткий анализ проблем наведения и посадки КЛА на малые тела солнечной системы.

Освоение космического пространства, изучение Солнечной системы, ее возникновения, формирования и строения, а также вопросы появления жизни на Земле в значительной степени связано с изучением малы тел солнечной системы: комет, астероидов и малых спутников планет. Большинство имеющихся научных работ основывается на материалах дистанционного зондирования этих тел. Дальнейшее изучение требует работы автоматических аппаратов на поверхности этих тел. В этом направлении уже сделаны первые шаги. Аппарат японского агентства аэрокосмических исследований (JAXA) Хаябуса взял пробы грунта астероида Итакава. Европейское космическое агентство с помощью спускаемого аппарата Филы надеялось получить сведения о комете P67 Чурюмова - Герасименко. Однако аппарат отклонился от заданной точки посадки, попал в тень и опрокинулся на бок. Исследования были выполнены в сокращенном объеме. Заметим, что точность наведения и посадки КЛА при исследовании малых тел имеет очень большое значение.

Другой научной проблемой, в решении которой важную роль играет наведение и посадка КЛА на поверхность малых тел солнечной системы, является отведение от орбиты Земли угрожающих космических объектов естественного происхождения. Анализ работ, посвященных изучению и отведению от Земли угрожающих малых тел солнечной системы, а также других способов защиты от астероидной опасности [1-5] показывает, что для решения задачи необходима посадка КЛА, либо наведение КЛА вплотную к поверхности МТСС и последующий взрыв [6]. В настоящее время координаты точки посадки КЛА рассчитываются предварительно [4], а посадка производится автономно по траекторным и баллистическим расчетам [7-10], однако, такой метод может не обеспечить требуемой точности, что в свою очередь приведёт к срыву отдельных задач или миссии в целом. Так спускаемый аппарат Филы при посадке на поверхность кометы P67/Чурюмова-Герасименко, отклонился на 100 метров от точки ожидаемой посадки. Отклонение имело крайне отрицательные последствия, аппарат попал в темную расщелину и кроме того, из-за неровностей поверхности опрокинулся на бок [10]. В случае защиты от астероидной опасности, такие промахи не допустимы. Отклонение от точки приложения силы на 100 метров при диаметре астероидов в 500 метров вполне может привести к срыву миссии [4].

Для совершения посадки КЛА на небесное тело с большей точностью необходимо специальное устройство навигации и посадки КЛА [11]. В миссиях «Фобос-Грунт» и спускаемом аппарате «Скиапарелли» европейского космического агентства для этой цели предусматривались системы технического зрения [12]. Анализ систем указанных выше аппаратов показывает, что эти системы были созданы методом расширения функционала других бортовых систем. Основными требованиями к такой системе являются обеспечение высокой точность наведения (с точностью в пределах разрешения исходного изображения, на котором производится выбора места посадки), помехозащищенность, а также стойкость к перспективным и ресурсным искажениям. Невысокие требования к бортовой вычислительной системе КЛА являются крайне важными для описываемого комплекса. Перечисленным требованиям в полной мере отвечают устройства оптической обработки информации, основанная на принципе оптической корреляции [13]. Рассмотрим далее некоторые теоретические аспекты возможности создания такого устройства.

Оптическая корреляция при наведении и посадки КЛА на малое тело солнечной системы.

Задача автономного наведения космического летательного аппарата (КЛА) по просматриваемой сцене (изображению поверхности МТСС) сводится, по существу, к совмещению двух функций - F1 (X, Y ) иF2 (X, Y ), формула (1). Первая из этих функций представляет собой распределение яркости на поверхности малого тела солнечной системы, вторая степень затемнения эталонного изображения поверхности объекта [14,15]. Если обе эти функции отличаются друг от друга ориентацией, т.е. одна из них переходит в другую при выполнении с ее аргументами преобразования, содержащего сдвиг начала координат на величину Y0, X0 , и угол поворота, ц :

 (1)

то критерием их совмещения можно считать достижение максимума величины

 (2)

Равенство в этом выражении достигается в единственном случае, когда

 (3)

Вычисление интеграла вида (2) аналоговым методом сравнительно просто удается выполнить при помощи оптического некогерентного коррелятора. Устройства подобного рода описаны в работе [16]. Некоторой сложностью в применении оптических корреляторов является различия в изображении одного и того же объекта. Различия эти связаны с целым рядом факторов:

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

· Фотосъемка и работа коррелятора будут производиться при разных дальностях до объекта посадки КЛА. При различных условиях освещения, возможно так же изменение условий видимости. Это приводит к тому, что контраст изображения отличается от контраста эталона, причем, контраст изображения может меняться при подлете КЛА к малому телу солнечной системы.

· Эталон и изображение имеют различную четкость. Связано это с тем, что объективы камеры и коррелятора имеют разную функцию передачи модуляции. Возможна потеря четкости изображения.

· Масштабирование и ракурсные искажения вследствие съемки места предполагаемой посадки и работы коррелятора во время самой посадки КЛА;

· Рабочее изображение и эталон различаются из-за разной дисторсии объектива камеры и объектива коррелятора;

· Искажения передачи яркости на эталоне ввиду логарифмической зависимости чувствительности матричного фотоприемника камеры КЛА;

В ряде работ [14-16] показано, что, несмотря на искажения перечисленные выше, по-прежнему, можно руководствоваться величиной интеграла (2) при довольно общих условиях, накладываемых на случайную функцию Ш(X, Y ), определяемую из соотношения:

 (4)

Неидентичность изображений, связанную с переменной контрастностью, указанную ранее, можно выразить следующим соотношением между функциями F 1 (X , Y ) и F 2 (X , Y ) [17]:

 (5)

Где g - фактор контрастности, e - постоянная составляющая.

C учетом зависимости (5), можно получить правило для совмещения двух изображений: совмещение изображений, описываемых функциями F 1 (X , Y ) и F 2 (X , Y ), будет оптимальным при таком выборе параметров Y 0 , ц и X 0 , когда величина

 (6)

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

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

Оптический коррелятор с вращающимся эталоном. Оптическое вычислительное устройство полярной корреляции с вращающимся эталоном.

Максимум величины R (6) совпадает с максимумом корреляционного интеграла (2) в одном частном случае, когда рабочее изображение отличается от эталона только поворотом на угол ц , а областью интегрирования является круг с центром в начале координат. В этом частном случае в выражении (6) корреляционный интеграл от зависит только угла поворота ц

 (7)

где A=const и B=const . При полном совпадении изображений (ц= 0), выражение имеет максимальное значение, которое достигается одновременно с максимальным значением корреляционного интеграла.

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

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

 (8)

Тогда зависимость фототока J от времени t , представляющую собой корреляционный интеграл можно записать в виде:

 (9)

где R - радиус кадрового окна, щ - частота вращения кадрового окна, K - коэффициент пропорциональности.

Функции щ изображений Ц1( с , ц ) и Ц2( с , ц ) периодичны по углу и допускают разложение в ряд Фурье по переменной ц . При n = m для выражения тока (10) справедливы следующие соотношения:

 (10)

Последнее выражение представляет собой разложение сигнала в ряд Фурье по временной координате.

Вводя обозначения:

соотношение (10) можно представить в виде:

 (11)

где

 (12)

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

 (13)

а также справедливы следующие соотношения:

 (14)

Тогда величины Bn в формулах (12) равны нулю, а вместе с ними равны нулю и фазы Шn всех гармонических составляющих сигнала, что приводит к возникновению симметричных пиков, максимумы которых соответствуют по времени совпадению изображений. При поступательном смещении изображений (нарушении соотношения (14)), каждая гармоническая составляющая приобретает отличную от нуля фазу, что приводит к исчезновению сигнала - пика. Приведенные выкладки показывают, что на сигнал коррелятора с вращающимся эталоном не влияют различия изображений, описываемые формулой (13) из-за разности контраста изображений. Таким образом, предлагаемая модель пеленгатора обладает свойством, которое называют нейтрализацией эталона. «Нейтрализация эталона» является важным свойством предлагаемого устройства наведения и навигации КЛА, поскольку посадка может выполнятся в различных условиях освещения и контрастности объекта.

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

Обсуждение результатов, полученных в статье.

Дальнейшее изучение малых тел солнечной системы актуальная научная задача, которая требует сбора и анализа состава вещества из которого состоят эти тела, уточнения параметров их движения и других исследований, которые могут быть выполнены в полном объеме при посадке на них КЛА. Посадка космического летательного аппарата на МТСС так же необходима в большинстве предложенных способов защиты Земли от астероидной опасности, отведения угрожающих небесных тел. Процедура посадки, как правило, должна выполняться на достаточно большом расстоянии, которое не позволяет управлять посадкой с Земли вручную в реальном масштабе времени. Очевидно, что посадка должна выполняться автономно в полностью автоматическом режиме. Для повышения точности этой процедуры, на взгляд авторов, целесообразно использовать отдельный специальный бортовой комплекс навигации и наведения КЛА. В настоящей статье предложен возможный вариант практической реализации дифференциального корреляционно-экстремального устройства наведения и навигации космического летательного аппарата, входящего в состав такого комплекса. В целом, дифференциальные оптические корреляционно-экстремальные устройства навигации и наведения обладают необходимыми для этого свойствами, в том числе помехоустойчивостью и автономностью. Наведение в таких системах осуществляется в местной системе координат путем сравнения с эталонным изображением. Однако, на данный момент наибольшее распространение получили системы с электронным вычислителем. Такая реализация предполагает наличие мощной вычислительной аппаратуры, которой, как правило, не располагают космические летательные аппараты. В отличие от ресурсоемких систем с электронным вычислителем, предлагаемое в данной статье устройство с корреляцией в оптическом тракте предъявляет минимальные требования к вычислительным ресурсам, что особенно ценно при реализации бортового комплекса навигации и наведения КЛА. Отсутствие высоких требований объясняется тем, что корреляционный алгоритм (преобразование Фурье) выполняется оптической подсистемой (аналоговый оптический вычислитель). Кроме того, оптические корреляционно-экстремальные устройства в составе бортового комплекса навигации и наведения КЛА обладают рядом других важных потенциальных достоинств. Оптические корреляционные устройства способны выполнять вычисления со скоростью до 1012 бит/с, что в настоящее время недостижимо другими экономически оправданными средствами [15]. Оптические корреляционно-экстремальные устройства непосредственно обрабатывают двухмерную пространственную информацию без дополнительных преобразований, в отличие от систем с электронным вычислителем, которые работают только по одному параметру - времени, а, следовательно, требуют дополнительных преобразований. Также весьма важным для бортового навигационно-посадочного комплекса КЛА являются такие свойства оптических корреляционно-экстремальных устройств, как высокая надежность, простота и компактность оптических деталей и узлов, а также малый вес и габариты аппаратуры.

Предлагается метод вращения эталонного изображения для отыскания экстремума функции корреляционного сигнала. Применение этого метода исключает необходимость «нейтрализации эталона», что в свою очередь, существенно упрощает конструкцию системы, значительно повышают помехозащищенность и другие характеристики оптического вычислительного корреляционного устройства обработки пространственной информации для целей навигации и наведения КЛА на малые тела солнечной системы. Показано, что предлагаемый пеленгатор с вращающимся эталоном может работать при различной четкости рабочего изображения и эталона. Это свойство пеленгатора существено важно при посадке на небесные тела, имеющие шлейф из газа и пыли, например, кометы, при этом разрешающая способность объектива пеленгатора может быть ограничена без ущерба для процесса отыскания максимума функции взаимной корреляции.

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

1. William Napier. Hazards from comets and asteroids, Oxford University Press, 2008. 289 p.

2. Медведев Ю.Д., Астероидно-комметная опасность / Ю.Д. Медведев, М.Л. Свешников, А.Г. Сокольский и др.-СПб.: Изд-во ИТА-МИПАО, 1996 - 244с.

3. Бурков В.Д., Есаков В.А., Куфаль Г.Э. и др. Проблема противодействия астероидной опасности космическими средствами. //Лесной вестник. 2011 № 5 С.157 - 169

4. Эйсмонт Н. А., Боярский М. Н., Ледков А. А., Назиров Р. Р., Данхэм Д., Шустов Б. М.. О возможности наведения малых астероидов на опасные небесные объекты с использование гравитационного маневра // Астрономический вестник, 2013, том 47, №4, С. 1-9.

5. Gehrels T (Ed), Hazards due to Comets and Asteroids, University of Arizona Press, 1994, ISBN 0-8165-1505-0 This covers all aspects of NEOs, comprehensively.

6. Coppinger R. NASA plans 'Armageddon' spacecraft to blast asteroid [Электронный ресурс] / Rob Coppinger // FlightGlobal. - 2007. - Режим доступа: https://www.flightglobal.com/news/articles/nasa-plans-armageddon-spacecraft-to-blast-asteroid-215924/ - (дата обращения 20.03.2018).

7. Карпенко С.А., Малые тела Солнечной системы посадка зонда NEAR на поверхность 433 Эроса [Электронный ресурс] // Новости космонавтики 2004. URL: http://galspace.spb.ru/index346.html (Дата обращения: 20.03.2018).

8. Takahiro Hiroi et al. Developing space weathering on the asteroid 25143 Itokawa // Nature. 2006. V. 443. P. 56-58.

9. ESA confirms the primary landing site for Rosetta. [Электронный ресурс] // ESA. 2015. URL:http://www.esa.int/Our_Activities/Space_Science/Rosetta/ESA_confirms_the_primary_landing_site_for_Rosetta. (дата обращения: 20.03.2018).

10. NASA Rosetta Instrument Reignites Debate on Earth's Oceans . [Электронный ресурс] // ESA. 2015. URL:http://www.jpl.nasa.gov/news/news.php?release=2014-423. (дата обращения: 20.03.2018).

11. Дубов С.С., Зельков К.В., Конопихин А.А., Краснопевцева Б.В., Лобанов А.А., Рожнев И.Ю., Шингарева К.Б. Выбор поверхности относимости и системы координат для картографо-геодезического обеспечения посадки на Фобос КА «ФОБОС-ГРУНТ» // Приложение к журналу Известия вузов. Геодезия и аэрофотосъемка. Сборник статей по итогам научно-технической конференции. 2009. № 2-2. С. 98-101.

12. Белинская, Е.В. История и перспективы использования систем технического зрения для управления процессом посадки на планеты и малые тела солнечной системы / Е.В. Белинская // Механика, управление и информатика (см. в книгах). - 2009. - №1. - С. 268-278.

13. Боголюбов, И.А. Применение оптического пеленгатора для целей посадки космического летательного аппарата на малые тела солнечной системы / И.А. Боголюбов, А.А. Лобанов, А.С. Филонов // Научный альманах. - 2017. №5-3(31). - С. 45-48.

14. Щербинин В.В. Построение инвариантных корреляционно-экстремальных систем навигации и наведения летательных аппаратов / В. В. Щербинин. - М. : Изд-воМГТУим. Н. Э. Баумана, 2011. - 230 с.

15. Якушенков Ю.Г. Проектирование оптико-электронных приборов. 2-е изд. - М.: Логос, 2000. - 489 с.

16. Сырямкин В.И., Шидловский В.С. Корреляционно-экстремальные радионавигационные системы. Томск: Изд-во Том. ун-та, 2010. - 316 с.

17. Благородов А.М., Ванециан Р.А. и др. Оптический коррелятора. - Заявка на изобретение №1567398/18-10.

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