О географической привязке космических снимков
Досліджено задачу географічної прив’язки космічних знімків, що включає процедуру польотного калібрування знімальної апаратури космічного апарата дистанційного зондування Землі. Обґрунтовано можливість реалізації польотного калібрування за інформацією від єдиного знімка наземного полігону з відомими...
Gespeichert in:
| Datum: | 2014 |
|---|---|
| 1. Verfasser: | |
| Format: | Artikel |
| Sprache: | Russian |
| Veröffentlicht: |
Інститут кібернетики ім. В.М. Глушкова НАН України
2014
|
| Schriftenreihe: | Проблемы управления и информатики |
| Schlagworte: | |
| Online Zugang: | https://nasplib.isofts.kiev.ua/handle/123456789/207835 |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| Назва журналу: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| Zitieren: | О географической привязке космических снимков / Д.В. Лебедев // Проблемы управления и информатики. — 2014. — № 5. — С. 71-79. — Бібліогр.: 14 назв. — рос. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraine| id |
irk-123456789-207835 |
|---|---|
| record_format |
dspace |
| spelling |
irk-123456789-2078352025-10-15T00:21:07Z О географической привязке космических снимков Про географічну прив'язку космічних знімків On geographical coordinate determination of space images Лебедев, Д.В. Космические информационные технологии и системы Досліджено задачу географічної прив’язки космічних знімків, що включає процедуру польотного калібрування знімальної апаратури космічного апарата дистанційного зондування Землі. Обґрунтовано можливість реалізації польотного калібрування за інформацією від єдиного знімка наземного полігону з відомими точковими орієнтирами. Наведено результати комп’ютерного моделювання, що підтверджують ефективність запропонованих алгоритмічних розв’язків. The problem of geographical coordinate determination of space images, including the in-flight calibration procedure of opto-electronic equipment of remote sensing SC has been investigated. The possibility of implementing calibration information from a single image with known ground polygon point landmarks has been substantiated. The results of computer simulation confirm the effectiveness of the proposed algorithmic solutions. 2014 Article О географической привязке космических снимков / Д.В. Лебедев // Проблемы управления и информатики. — 2014. — № 5. — С. 71-79. — Бібліогр.: 14 назв. — рос. 0572-2691 https://nasplib.isofts.kiev.ua/handle/123456789/207835 330.4:658.012.32 10.1615/JAutomatInfScien.v46.i9.40 ru Проблемы управления и информатики application/pdf Інститут кібернетики ім. В.М. Глушкова НАН України |
| institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
| collection |
DSpace DC |
| language |
Russian |
| topic |
Космические информационные технологии и системы Космические информационные технологии и системы |
| spellingShingle |
Космические информационные технологии и системы Космические информационные технологии и системы Лебедев, Д.В. О географической привязке космических снимков Проблемы управления и информатики |
| description |
Досліджено задачу географічної прив’язки космічних знімків, що включає процедуру польотного калібрування знімальної апаратури космічного апарата дистанційного зондування Землі. Обґрунтовано можливість реалізації польотного калібрування за інформацією від єдиного знімка наземного полігону з відомими точковими орієнтирами. Наведено результати комп’ютерного моделювання, що підтверджують ефективність запропонованих алгоритмічних розв’язків. |
| format |
Article |
| author |
Лебедев, Д.В. |
| author_facet |
Лебедев, Д.В. |
| author_sort |
Лебедев, Д.В. |
| title |
О географической привязке космических снимков |
| title_short |
О географической привязке космических снимков |
| title_full |
О географической привязке космических снимков |
| title_fullStr |
О географической привязке космических снимков |
| title_full_unstemmed |
О географической привязке космических снимков |
| title_sort |
о географической привязке космических снимков |
| publisher |
Інститут кібернетики ім. В.М. Глушкова НАН України |
| publishDate |
2014 |
| topic_facet |
Космические информационные технологии и системы |
| url |
https://nasplib.isofts.kiev.ua/handle/123456789/207835 |
| citation_txt |
О географической привязке космических снимков / Д.В. Лебедев // Проблемы управления и информатики. — 2014. — № 5. — С. 71-79. — Бібліогр.: 14 назв. — рос. |
| series |
Проблемы управления и информатики |
| work_keys_str_mv |
AT lebedevdv ogeografičeskojprivâzkekosmičeskihsnimkov AT lebedevdv progeografíčnuprivâzkukosmíčnihznímkív AT lebedevdv ongeographicalcoordinatedeterminationofspaceimages |
| first_indexed |
2025-10-15T01:14:11Z |
| last_indexed |
2025-10-16T01:06:52Z |
| _version_ |
1846098591789088768 |
| fulltext |
© Д.В. ЛЕБЕДЕВ, 2014
Международный научно-технический журнал
«Проблемы управления и информатики», 2014, № 5 71
КОСМИЧЕСКИЕ ИНФОРМАЦИОННЫЕ
ТЕХНОЛОГИИ И СИСТЕМЫ
УДК 528.7
Д.В. Лебедев
О ГЕОГРАФИЧЕСКОЙ ПРИВЯЗКЕ
КОСМИЧЕСКИХ СНИМКОВ
Введение
Под географической привязкой снимков Земли из космоса будем понимать
процедуры обработки снимка земной поверхности, в результате которых каждому
элементу изображения ставятся в соответствие его географические координаты
(широта и долгота).
Стремительное развитие и совершенствование элементной базы оптико-элек-
тронного оборудования космических аппаратов (КА) дистанционного зондирова-
ния Земли (ДЗЗ), использование современных ПЗС-матриц, работающих по тех-
нологии временной задержки и накопления сигналов (time delay and integration),
стимулируют разработку таких информационных технологий привязки изображе-
ний земной поверхности, которые обеспечивали бы точность, адекватную точно-
сти пространственного разрешения съемочной аппаратуры. Чтобы иметь пред-
ставление, о каких уровнях точности пространственного разрешения идет речь,
отметим, что для съемочной аппаратуры КА очень высокого и сверхвысокого раз-
решения размеры наименьших объектов, различимых на изображении, составля-
ют 1–10 м и меньше 0,3–0,9 м соответственно [1].
Из всех возмущающих факторов наибольшее влияние на точность привязки
космических снимков оказывает остаточная неопределенность ориентации систе-
мы координат целевой аппаратуры относительно бортовых средств носителей
опорной системы координат, например, звездного датчика (ЗД). Как отмечено
в [2], отклонения взаимного положения упомянутых систем координат от начально-
го положения, определенного в результате предполетной геометрической калибров-
ки, обусловлены воздействием ряда возмущающих факторов, таких как термо-
и вибронагрузки, имеющие место при транспортировке, запуске КА и в процессе
его функционирования на орбите, а также невесомость, старение материалов конст-
рукции КА и др.
Высокая точность координатной привязки изображений возможна за счет
применения космических платформ с высокой стабильностью и улучшенной точ-
ностью определения ориентации спутника ДЗЗ. Под стабильностью платформы
здесь понимается возможность длительного поддержания допустимых значений
параметров, характеризующих взаимную ориентацию съемочной аппаратуры
и измерительных средств системы управления ориентацией и стабилизацией КА.
При отсутствии такой возможности возникает необходимость периодического
уточнения параметров взаимной ориентации указанного оборудования КА в про-
цессе его функционирования (полетная калибровка).
Известны, по крайней мере, два способа решения задачи полетной калибровки,
отличающиеся используемыми источниками внешней информации: калибровка по
звездам и калибровка в режиме наблюдения наземных полигонов с известными опор-
ными объектами (например, точечными ориентирами) [1]. Калибровка по наземным
72 ISSN 0572-2691
точечным ориентирам может проводиться как с использованием топографически при-
вязанных (топопривязанных) ориентиров, так и по априори неизвестным ориентирам
(см., например, [3–10]). В обоих случаях искомые параметры калибровки определяют-
ся из решения нормальных уравнений метода наименьших квадратов (МНК). Следует
отметить, что точность решения задачи полетной калибровки существенным образом
зависит от количества измерений координат наблюдаемых ориентиров (числа контак-
тов с ориентирами), влияющего на обусловленность системы уравнений МНК. Это об-
стоятельство может сказаться на оперативности и достоверности получения информа-
ции о качестве производимой съемки. В связи с этим перспективным направлением
исследований представляется разработка методов полетной калибровки оптико-элект-
ронного комплекса КА ДЗЗ при наличии минимального числа снимков земной по-
верхности, содержащих топопривязанные ориентиры.
Учитывая, что высокоточное и оперативное решение задачи полетной калиб-
ровки играет приоритетную роль в реализации качественной географической при-
вязки космических снимков, в этой работе сначала ставится и решается именно
эта задача. При этом используется информация из единственного снимка назем-
ного полигона с известными опорными точками. Затем реализуется процедура
географической привязки элементов кадра и оценивается ее точность.
1. Системы координат. Постановка задачи полетной калибровки
Пусть КА ДЗЗ обращается вокруг Земли по орбите, близкой к круговой. На
нем установлены съемочная аппаратура (камера) высокого разрешения, звездный
датчик и GPS-приемник. Не нарушая строгости, будем считать, что упомянутая
аппаратура установлена в центре масс О космического аппарата.
Введем необходимые далее правые ортогональные системы координат (СК)
и базисы, составленные из ортов осей этих СК:
XYZ (I) — инерциальная СК с вершиной в центре Земли, осью X, направ-
ленной в точку весеннего равноденствия, и осью Z, ориентированной по оси мира
в сторону Полярной звезды;
GGG ZYX (G) — географическая СК, в качестве которой принимается WGS-84
(World Geodetic System 1984). В WGS-84 нулевым меридианом считается IERS Refer-
ence Meridian, расположенный в 13,5 к востоку от Гринвичского меридиана. За осно-
ву взят эллипсоид с экваториальным радиусом 1R 6378137 м и полярным радиусом
2R 6356752,3142 м. Угловая скорость вращения Земли — 510292115,7 рад/c [11];
О123 (S) — СК, связанная со звездным датчиком так, что ось 3 совпадает
с оптической осью датчика;
Oxyz (K) — СК, связанная со съемочной камерой таким образом, что ее
вершина находится в центре проекции камеры, ось z перпендикулярна чувстви-
тельной площадке (плоскости изображения) и направлена по оптической оси ка-
меры в сторону, противоположную объекту съемки;
TTTT ZYXO (Т) — топоцентрическая СК, вершина которой находится в
точке TO земной поверхности с координатами и (широта и долгота соот-
ветственно). Ось TY совпадает с касательной к меридиану и направлена на север,
а ось TZ — вертикально вверх.
Связь между базисами I, G, S, K и T иллюстрирует схема на рис. 1, где ijC —
матрица вращения, осуществляющая преобразование векторов из базиса i в базис j.
Сформулируем задачу полетной геометрической калибровки. Пусть известны
ориентация топоцентрической СК относительно инерциального пространства — мат-
рица направляющих косинусов ,ITC местоположение и ориентация космического ап-
парата в инерциальной СК — радиус-вектор r и матрица ISС (рис. 2), а также радиу-
сы-векторы ir ),1( ni топографически привязанных точечных ориентиров iM
),1( ni , находящихся в поле зрения камеры. По имеющейся информации и единст-
Международный научно-технический журнал
«Проблемы управления и информатики», 2014, № 5 73
венному снимку наземного полигона, полученному с орбиты полета КА, не обходи-
мо уточнить значение матрицы ,SKC определяющей взаимную ориентацию базисов
S и K звездного датчика и камеры.
r
jVr
iMr ijr
jVe
iMe
O
TO jV
iM TY
TZ
TX
Рис. 2
2. Алгоритм полетной калибровки
Если матрицу направляющих косинусов KSС параметризовать координатами
вектора ориентации eθ },,{ 321 или вектора Эйлера (e — орт эйлеровой
оси вращения, θ — угол поворота относительно этой оси), то ее связь с ко-
ординатами вектора θ задается соотношениями [12]
),(
cos1
)(
sin 2
23 θθ
EСKS
,
0
0
0
)(
12
13
23
θ
где 3E — единичная матрица третьего порядка. При малых значениях угла , ко-
гда 2 значительно меньше 1, матрица KSС может быть аппроксимирована вы-
ражением ).(3 θ EСKS Тогда матрица SKC имеет вид
).(3 θ EСSK (1)
Из равенства (1) следует, что задача калибровки сводится к определению ко-
ординат вектора .θ
Из n известных точек в топоцентрической системе координат TZTTT YXO
выберем точку iM и произвольную точку ),,1( mjV j координаты которой
можно задавать произвольным образом, — виртуальную точку (рис. 2).
При известных векторах
iMrr, и
jVr вычислим орты
jViM ee , направлений
iOM и ,jOV а также вектор .
ij MVij rrr В основу алгоритма уточнения зна-
чения матрицы SKC положим условие компланарности векторов
iMij rr , и :
jVr
0)(T
jViMij eer (2)
(верхний индекс Т означает транспонирование).
Запишем условие (2) в системе координат .TTTT ZYXO Из показаний камеры
},,{ Fyx ii , где F — фокусное расстояние съемочной аппаратуры, сформируем
единичный вектор KiM ,e направления на i-й ориентир. Как следует из схемы на
рис. 1, в топоцентрической СК он определяется соотношением
.,, KiMKTTiM C ee (3)
STC KTC
GTC
ITC
SKC
S
IGC
G T
ISC
K I
Рис. 1
74 ISSN 0572-2691
Из этой же схемы находим выражения, необходимые для вычисления матрицы :KTC
,T
SKSTKT CCC .T
ISITST CCC
Тогда равенство (3) с учетом формулы (1) можно записать в виде
.,))((, 3 KiMSTTiM EC eθe (4)
При известных векторах r и
jVr орт
jVe направления из точки О на выбран-
ную виртуальную точку jV вычисляется по формуле
./)( rrrre
jjj VVV (5)
Векторное произведение в условии (2) с учетом равенств (4) и (5) определя-
ется соотношением
,),()(,)( θeeeeee STKiMSTVKiMSTVViM CCC
jjj
которое позволяет записать условие компланарности (2) в форме
.,)(
,),()(B
,B
T
T
KiMSTVij
STKiMSTVij
Cq
CC
q
j
j
eer
eer
θ
(6)
При обработке информации обо всех априори известных и виртуальных ориен-
тирах число уравнений типа (6) равно mn. В связи с этим для формирования уравне-
ний МНК используется схема, приведенная в [8]. Решая полученную с ее помощью
систему алгебраических уравнений относительно вектора θ , уточняем взаимную
ориентацию базисов S и K звездного датчика и съемочной камеры, принимая это
в качестве результата работы собственно алгоритма калибровки.
Будем считать известными статистические характеристики точности информа-
ции, используемой для реализации изложенной процедуры оценки матрицы .SKC Ме-
тодом статистических испытаний выполним расчет SN вариантов оценки координат
вектора θ , используя в каждом варианте расчета свою систему виртуальных точек.
Усредненное по множеству решений θ значение вектора θ (математическое ожида-
ния )θ̂ принимается в качестве решения задачи полетной калибровки оптико-
электронной аппаратуры КА ДЗЗ:
.
1
}ˆ,ˆ,ˆ{ˆ
1
321
θθ
SN
SN
(7)
Чтобы получить сведения о точности и надежности оценки (7), воспользуем-
ся системой доверительных интервалов )3,2,1(, iI i [13] для координат векто-
ра .θ̂ Для этого параллельно с расчетом векторов θ будем формировать их по-
грешность в виде разности
θθδ }{ ,i
с априори известным вектором рассогласования θ базисов S и K.
Характеристики случайных координат )3,2,1(, ii векторов δ — мате-
матическое ожидание im , и дисперсия iD , — неизвестны. Однако считаем, что
в процессе моделирования для них получены оценки
).3,2,1()~(
1
1~
,~ 2
,,
1
,,
1
1
,
im
N
DNm ii
N
S
ii
N
Si
SS
Международный научно-технический журнал
«Проблемы управления и информатики», 2014, № 5 75
Построим доверительный интервал iI , для математического ожидания ,,im
соответствующий доверительной вероятности . Согласно [13] он определяется
формулой
).~;~( ,,,,, iiiii tmtmI (8)
В соотношении (8) ./
~
,, Sii ND Величина t определяет для нормального
закона распределения число средних квадратических отклонений, которое следу-
ет отложить слева и справа от центра рассеивания, чтобы вероятность попадания
в полученный интервал была равна . Численное значение t зависит от и, на-
пример, для 999,0 равно 3,29.
Выполнив расчеты доверительных интервалов для каждой из координат вектора
δ при заданных значениях SN и доверительной вероятности , можно оценить диа-
пазоны возможных значений погрешности работы алгоритма полетной калибровки.
3. Географическая привязка элементов кадра. Точность привязки
Пусть на снимке земной поверхности в системе координат ,Oxyz связанной
с камерой, выбрана точка Р с координатами x, y и .F Необходимо вычислить
координаты этой же точки в географической СК, т.е. осуществить ее привязку.
Из координат точки Р сформируем единичный вектор .Ke После того как вы-
полнена процедура полетной калибровки, орт Pe направления на эту точку из точки
съемки О, представленный в системе ,GGG ZYX определяется равенством (рис. 1)
.))ˆ((},,{ 3 KSGKKGP ECCnml eθee (9)
Как и в работе [7], в качестве координат точки Р в системе WGS-84 примем ко-
ординаты точки пересечения прямой, проходящей через точку съемки O с координа-
тами
GGG ZYX ,, параллельно направляющему вектору ,Pe с эллипсоидом вра-
щения, принятым в WGS-84. Таким образом, искомые координаты — решение
системы уравнений
.1
,
2
2
2
2
1
2
2
1
2
***
R
Z
R
Y
R
X
n
ZZ
m
YY
l
XX
GGG
GGGGGG
(10)
Система алгебраических уравнений (10) имеет два решения:
.,,
,,,
2
22222
2
11111
qppZZYZX
qppZZYZX
GGGGG
GGGGG
(11)
В выражениях (11) введены следующие обозначения:
,,/,,/ ****
GGGG ZYnmZXnl
)).(/(,)(,)( 222
2
2
1
2
2
2
1
22 RRRDDRqDp
Из двух решений системы (10) выбирается то, которому соответствует точка,
ближайшая к точке O ).,,(
GGG ZYX Ее радиус-вектор обозначим
}.,,{ GGG ZYXE
Для вычисления географических координат PP , точки Р достаточно оп-
ределить значения координат орта EEe /},,{ ZYXG eee вектора Е и вос-
пользоваться равенствами
76 ISSN 0572-2691
.sin,coscos,coscos PZPPYPPX eee
Рассмотрим два варианта оценки точности координатной привязки космиче-
ского снимка.
Вариант 1. На снимке наземного полигона выберем один из известных ориенти-
ров и по методике, изложенной выше, вычислим его координаты в топоцентрической
СК. Из сравнения найденных координат ориентира с априори известными координа-
тами оценим точность его привязки. Выполнив такие расчеты со всеми известными и
виртуальными ориентирами при различных реализациях случайных погрешностей
бортовой аппаратуры, информация которой необходима для топопривязки, определим
статистические характеристики точности координатной привязки снимка.
Вариант 2. Видоизменим рис. 2,
приведя его к виду, представленному на
рис. 3. Введем векторы
ijV
i
eρ
)(
1
i)(
1 и
.
)(
2
i)(
2 iM
i
eρ Совместно с вектором jb
они характеризуют стереопару, содержа-
щую известный ориентир iM на земной
поверхности. Чтобы найти оценку
iMr̂ век-
тора ,
iMr воспользуемся приемом, изло-
женным в [9]. Согласно этому приему век-
тор
iMr̂ определяется по формулам
,2/)(ˆ )(
2
)(
1
ii
VM ji
ρρrrr
.,
)(
,
)( T
)(
2
T
)(
1
iji
iji
i
iji
iji
iji
ii
VM
VM
d
VM
Vjdi
VM
Mjdi
ee
ee
e
ee
ebe
ρ
ee
ebe
ρ
Аналогично варианту 1 характеристики точности привязки оцениваются ме-
тодом статистических испытаний.
4. Результаты моделирования
На первом этапе моделирования процесса географической привязки снимка
в соответствии с методикой, описанной в разд. 2, оцениваются параметры взаим-
ной ориентации базисов S и K, на втором этапе исследуется точность координатной
привязки снимка. Моделирование выполняется по технологии, изложенной в [9].
Система параметрических возмущений, оказывающих влияние на точность
полетной калибровки, формируется в виде совокупности нормально распределен-
ных центрированных случайных величин и включает:
— погрешности звездного датчика 3),2,1,( ii которые имитируются по-
воротами базиса S вокруг осей 1, 2, 3 в его расчетном положении и имеют средние
квадратические отклонения ;i
— аддитивные погрешности ,,, YX Z возникающие при формировании
координат центра масс КА по информации от GPS-приемника. Средние квадрати-
ческие отклонения погрешностей считаются одинаковыми и равными ;GPS
— аддитивные погрешности фокусного расстояния F съемочной аппарату-
ры (среднее квадратическое отклонение — );F
— ошибки считывания координат kk yx , k-го ориентира на плоскости изо-
бражения камеры из-за их пиксельного представления; они подчиняются равно-
мерному закону распределения.
Координаты вектора ,θ характеризующие начальную неопределенность
взаимной ориентации базисов S и K, принимались центрированными нормально
распределенными случайными величинами с одинаковыми среднеквадратически-
ми отклонениями .
O
ijVe
jb
r
jVr
iMr
iMe
TO
jV
iM TY
TZ
TX
Рис. 3
Международный научно-технический журнал
«Проблемы управления и информатики», 2014, № 5 77
Полетная калибровка. Предполагалось, что движение КА происходит по
круговой орбите со следующими параметрами: высота — 680 км, наклонение ор-
биты — 98 °, долгота восходящего узла — 142 °.
Моделирование выполнялось для случая, когда опорные точки размещались
на наземном полигоне размером м2200900 на высоте из диапазона м.10
Виртуальные точки выбирались на плоскости ;TTT YXO размер площадки —
м.60006000 Координаты n известных и m виртуальных ориентиров формиро-
вались с использованием системы псевдослучайных равномерно распределенных
чисел. Географические координаты вершины TO топоцентрической системы ко-
ординат — .33,45 Параметры съемочной аппаратуры: фокусное рас-
стояние — ,м25,2F среднее квадратическое отклонение — м,0,0075F
размер пиксела — мкм.99 Точность стабилизации углового движения КА во
время съемки характеризовалась среднеквадратическим отклонением в пределах
0,5 ° по каждой из степеней свободы.
При моделировании анализировалось влияние двух наборов погрешностей звезд-
ного датчика и GPS-приемника на точность калибровки: I — 20}2,{2,σ угл. с,
15GPS м; II — 0}1,1{1,σ угл. с, 5GPS м. Значение вектора σ в наборе I
отвечает характеристикам звездного датчика БОКЗ-М [14]. Использовалось оди-
наковое число известных (опорных) и виртуальных точек ).7( nm
Вначале генерировались координаты опорных точек на полигоне (или задава-
лись — при заранее известных координатах) и значения координат вектора ,θ кото-
рые оставались постоянными в процессе калибровки. Затем по алгоритму разд. 2 рас-
считывались варианты уточнения параметров взаимной ориентации базисов S и K.
Этот этап моделирования заканчивался расчетом оценки θ̂ вектора θ и фор-
мированием диапазонов возможных погрешностей процедуры калибровки.
Координаты вектора ошибки оценки ,θ̂ полученные в результате моделиро-
вания трех вариантов задания начального рассогласования вектора θ (его коор-
динаты генерировались при )01
T)3(T)2(T)1(
10,17]10,552,70[,10,55]3,33[25,24,1,64]18,848,93[ θθθ
и наборов I и II погрешностей звездного датчика и GPS-приемника, приведены
в табл. 1. Данные из таблицы получены при 10000SN и доверительной вероят-
ности .999,0 Они свидетельствуют, во-первых, о возможности реализации
процесса полетной калибровки по одному снимку подспутникового полигона с
топографически привязанными ориентирами и, во-вторых, о высокой точности
решения задачи калибровки. По результатам выполненных расчетов при началь-
ном рассогласовании базисов S и K в несколько десятков угловых минут его зна-
чение уточняется по координатам 1 и 2 с ошибкой порядка 01 (для набора I
характеристик камеры и ЗД) и около 7 (для набора II).
Таблица 1
№
варианта
Координаты вектора ошибки оценки ,θ̂ угл. с
I II
1 16,003,2 17,063,5 31,269,9 06,001,0 06,010,2 62,134,5
2 16,089,1 17,015,10 31,217,12 06,012,0 07,058,6 28,283,7
3 16,023,1 17,092,4 36,269,6 06,077,0 06,038,1 30,233,2
Координатная привязка снимка. Полученное в результате калибровки
уточненное значение θ̂ вектора θ (матрицы )SKC ) позволяет либо воспользо-
ваться формулой (9) для вычисления направляющего вектора Pe в варианте 1 то-
попривязки снимка, либо откорректировать орт
iMe — в варианте 2.
78 ISSN 0572-2691
Для большей наглядности точность топопривязки точки iM будем характери-
зовать не разностью между вычисленными значениями широты и долготы и их точ-
ными значениями, а вектором
iii MMM rrr ˆ в топоцентрической системе ко-
ординат
iMr̂( и
iMr — оценка радиуса-вектора точки iM и его точное значение).
Моделирование с целью определения точности топопривязки при фиксиро-
ванном значении вектора θ̂ было организовано следующим образом: каждый ва-
риант моделирования включал генерирование расположения систем опорных то-
чек и виртуальных ориентиров, формирование описанной выше совокупности
случайных параметрических возмущений и оценку координат опорных точек.
Чтобы судить об эффективности и точности процедуры топографической
привязки, приведем статистические характеристики (математические ожидания
Rm и средние квадратические отклонения )Rσ ошибок топопривязки.
Для случая, когда остаточная неопределенность во взаимной ориентации ба-
зисов S и K не уточнялась в процессе полета КА, в табл. 2 представлены резуль-
таты расчетов, выполненных для тех же значений вектора начального рассогласо-
вания ,θ что и в табл. 1. Радиус-вектор точки iM оценивался по методике вари-
анта 2 топопривязки.
Таблица 2
№
варианта
Координаты вектора mR математическо-
го ожидания ошибки топопривязки, м
Координаты вектора R средней квадра-
тической ошибки топопривязки, м
1 – 691,3 – 3005,0 6,2 293,1 260,9 10,6
2 3562,5 1288,9 9,6 322,7 354,6 10,4
3 – 621,8 938,8 0,8 104,6 99,8 3,4
Данные из табл. 2, полученные без использования процедуры полетной калиб-
ровки съемочной аппаратуры КА, свидетельствуют об очень грубой топографиче-
ской привязке снимка (по отдельным координатам погрешность превышает 4 км).
При моделировании процесса топопривязки с учетом данных полетной калиб-
ровки получены следующие результаты. Математические ожидания Rm и сред-
ние квадратические отклонения Rσ ошибок топопривязки для трех вариантов за-
дания значений вектора начального рассогласования θ одинаковые и имеют такие
численные значения (в метрах): ,]07,144,3[ TRm T]1,04,179,16[Rσ — для
набора характеристик I; ,]09,43,0[ TRm T]04,62,6[Rσ — для набора ха-
рактеристик II.
Эффективность топопривязки с использованием данных, полученных в ре-
зультате решения задачи полетной калибровки по информации из единственного
снимка подспутникового полигона, очевидна. Следует отметить почти троекрат-
ное повышение точности топопривязки при использовании звездного датчика
и GPS-приемника с улучшенными характеристиками (набор II).
Заключение
В работе исследуется задача топографической привязки элементов объекта
съемки, в которую включена процедура полетной калибровки съемочной аппара-
туры КА ДЗЗ. Обоснована возможность реализации процесса калибровки по ин-
формации из единственного снимка наземного полигона с топографически привя-
занными ориентирами, полученного с орбиты. Ключевой момент процедуры ка-
либровки — использование системы так называемых виртуальных ориентиров.
Математическое моделирование алгоритмов полетной калибровки и топо-
привязки, выполненное на основе метода статистических испытаний, свидетель-
ствует об эффективности предлагаемых алгоритмических решений.
Международный научно-технический журнал
«Проблемы управления и информатики», 2014, № 5 79
Д.В. Лебедєв
ПРО ГЕОГРАФІЧНУ ПРИВ’ЯЗКУ
КОСМІЧНИХ ЗНІМКІВ
Досліджено задачу географічної прив’язки космічних знімків, що включає про-
цедуру польотного калібрування знімальної апаратури космічного апарата дис-
танційного зондування Землі. Обґрунтовано можливість реалізації польотного
калібрування за інформацією від єдиного знімка наземного полігону з відоми-
ми точковими орієнтирами. Наведено результати комп’ютерного моделювання,
що підтверджують ефективність запропонованих алгоритмічних розв’язків.
D.V. Lebedev
ON GEOGRAPHICAL COORDINATE
DETERMINATION OF SPACE IMAGES
The problem of geographical coordinate determination of space images, including the
in-flight calibration procedure of opto-electronic equipment of remote sensing SC has
been investigated. The possibility of implementing calibration information from
a single image with known ground polygon point landmarks has been substantiated.
The results of computer simulation confirm the effectiveness of the proposed algo-
rithmic solutions.
1. Классификация снимков по пространственному разрешению. — http://mapexpert.com.ua.
2. Самойлов С.Ю. Способы уменьшения погрешности географической привязки снимков
с космических аппаратов дистанционного зондирования Земли // Вестник ФГУП НПО
им. С.А. Лавочкина. — 2012. — № 2. — С. 19–22.
3. Лебедев Д.В., Ткаченко А.И. Калибровка информационно-измерительного комплекса кос-
мического аппарата, предназначенного для съемки земной поверхности // Проблемы
управления и автоматики. — 2004. — № 1. — С. 101–120.
4. In-flight geometric calibration — an experience with Cartsat-1 and Cartsat-2 / T.P. Srinivasan,
B. Islam, S.K. Singh, Crishna B. Copala, P.K. Srivastava //Archives of the Photogrammetry and
Remote Sensing and Spatial Information Sciences. XXXVII. Part B1. Beijing, 2008. — P. 83–88.
5. Сомов Е.И., Бутырин С.А., Скирмунт В.К. Полетная геометрическая калибровка космиче-
ского телескопа и системы звездных датчиков // Сб. материалов юбилейной XV Санкт-
Петербургской международной конференции по интегрированным навигационным систе-
мам. СПб: ЦНИИ «Электроприбор», 2008.
6. Сомов Е.И., Бутырин С.А. Технология обработки сопровождающей измерительной инфор-
мации для высокоточной координатной привязки космических снимков // Известия Самар-
ского научного центра РАН. — 2009. — 11, № 5. — С. 156–163.
7. Пятак И.А. Задачи координатной привязки снимков, выполненных КА // Вісник Дніпропет-
ровського університету. Сер. «Ракетно-космічна техніка». — 2011. — Вип. 14. — С. 116–122.
8. Лебедев Д.В., Ткаченко А.И. Параметрическая юстировка комплекса «камера и звездный
датчик», установленного на низкоорбитальном космическом аппарате // Известия РАН.
Теория и системы управления. — 2012. — № 2. — С. 153–165.
9. Лебедев Д.В. Полетная геометрическая калибровка оптико-электронной аппаратуры кос-
мического аппарата наблюдения Земли по неизвестным ориентирам // Международный на-
учно-технический журнал «Проблемы управления и информатики». — 2013. — № 5. —
С. 114–125.
10. Ткаченко А.И. О полетной юстировке оптико-электронного комплекса космического аппа-
рата // Известия РАН. Теория и системы управления. — 2013. — № 6. — С. 122–130.
11. WGS84. — http://ru.wikipedia.org/wiki/WGS-84.
12. Бранец В.Н., Шмыглевский И.П. Введение в теорию бесплатформенных инерциальных на-
вигационных систем. — М. : Наука, 1992. — 280 с.
13. Вентцель Е.С. Теория вероятностей. — М. : Наука, 1964. — 576 с.
14. Основные характеристики звездных координаторов БОКЗ. — http://www.iki.rssi.ru/ofo/
bokz_spec.html.
Получено 30.05.2014
http://mapexpert.com.ua/
http://ru.wikipedia.org/wiki/WGS-84
|