Фотограмметрія Дорожинський

Розміщення матеріальної точки і визначається вектором η, фіксується такими величинами^ η - геоцентричним радіусом, at - прямим сходженням, δί -

Екваторіальна система координат безперервно змінює своє розміщення через вплив прецесії (рух осі власного обертання по коловій конічній поверхні; тоді точка весняного рівнодення зміщується; загальна річна прецесія для Землі становить 50,3") та нутації (коливання земної осі чи осі Світу, зумовлені збуренням Сонця і Місяця). Для того, щоб така система стала інерційною, її необхідно зафіксувати на певну епоху Т0, якою прийнята епоха Т0 = 2000,0. Інерційна система координат

встановлюється каталогами координат зірок. У космічній фотограмметрії використовують католог Смітсоніанської астрофізичної обсерваторії (SAOC), створений на підставі спостережень більше ніж 250 тис. зірок.

Для фіксації положення точок земної поверхні інерційна система незручна, бо треба вказати ще й час, якого стосуються ці координати. Адже через вплив гравітаційних сил полюси Землі змінюють своє положення на поверхні планети. Тому для Землі прийнята Грінвіцька геоцентрична система координат (рис. 8.9). Тут початок 0 збігається з центром мас Землі, вісь oz скерована в середній Північний полюс Землі, вісь ох лежить у площині геоцентричного Грінвіцького меридіана та площині екватора. Вісь оу доповнює цю систему до просторової прямокутної, як показано на рис. 8.9.

У літературі наводяться формули, за якими можна від істинних інерційних координат α, δ зафіксованих на епоху t, перейти до грінвіцьких координат.

300

Рис. 8.9, Грінвіцька система координат (Північний полюс - середнє значення у 1900—1905 pp., площина екватора - середнє положення у 1900-1905 pp.)

У фотограмметрії (як і в геодезії) широко використовуються геодезичні системи координат. Початок геодезичної системи, як правило, не збігається з центром мас Землі, проте осі координат майже паралельні до осей грінвіцької системи. Розбіжність становить 1-2" і встановлюється на підставі високоточних астрономічних спостережень. Для конкретної математичної моделі еліпсоїда для Землі - т. зв. референц-еліпсоїда - геодезичні координати точки задаються

широтою В , довготою L та висотою Η (див. п. 2.3).

Для вивчення планет (побудова опорних мереж та картографування) використовують планетоцентричні та планетографічні системи координат. Вони зв'язані з прийнятою для планети референц-поверхнею.

Упланетоцентричній системі початок розміщений у центрі мас планети, вісь

ζзбігається з положенням середньої осі обертання планети, площина ху

збігається з площиною середнього екватора. Тут бачимо повну аналогію з рис. 8.9, проте замість Грінвіцького меридіана треба взяти щось своє. Таким нульовим меридіаном є меридіан, що проходить через наперед вибрану точку планети. Наприклад, для Меркурія це кратер Хун Кал; для Марса - це кратер - Ері- О; для Юпітера та Сатурна нульові меридіани виводять з вимірів магнітних полів.

301

Планетоцентричні широти φ і довготи λ як сферичні координати зв'язані з просторовими прямокутними координатами рівняннями, аналогічними до (8.7).

ζ= r sin φ.

Улітературі подаються формули для переходу від прямокутних координат для стандартної епохи до прямокутних координат в описаній планетоцентричній системі. Для цього використовують дані МАС (Міжнародного астрономічного союзу) для кожної планети: координати північного полюсу обертання та напрямку нульового меридіана (на епоху В 1950 або на епоху 2000).

Для картографування поверхні планет використовують планетографічну систему. В ній положення точки задається широтою Φ\ довготою Л' та висотою Η .

Планетографічна широта φ' - кут між нормаллю до референц-поверхні та площиною екватора планети.

Планетографічна довгота Л' - кут, вимірюваний по екватору від нульового меридіана до меридіана точки.

Планетографічна висота Η - віддаль по нормалі від точки до поверхні референц-еліпсоїда.

Планетоцентричні та планетографічні координати обчислюють так:

λ= 360° - Л' (для планет з прямим обертанням);

λ= Λ' (для планет з оберненим обертанням);

е - а ~ с - полярне стиснення;

а

а- екваторіальний радіус;

с- полярний радіус.

Згідно з рекомендаціями МАС для планет із прямим обертанням довготи ростуть від 0° до 360° від нульового меридіана праворуч, а для зворотного обертання - ліворуч.

302

8.3.4. Визначення планетоцентричних координат точок планети із застосуванням кадрових систем

Кадрова знімальна система дає змогу отримати знімок (образ) як центральну проекцію, в якій точка зображення визначається або прямокутними координатами (х,у- / ) , або сферичними координатами (a9j39r) (рис. 8.9).

Отримуючи формули, що зв'язують точку т знімка з точкою Μ поверхні планети, виходять з наявності вихідних даних, які можуть бути у різних комбінаціях. Наприклад, можуть бути відомі кутові елементи орієнтації знімальної камери в системі космічного апарата, або ж кутове орієнтування камери та елементи орбіти в планетоцентричній системі координат на конкретну епоху прецесійні елементи як відхилення епохи t від стандартної епохи 2000, чим визначається положення нульового меридіану та Північного полюса планети.

Як видно з рис. 8.9, для точок S , т , Μ справедлива умова колінеарності; тому за повною аналогією з аерозніманням запишемо класичні рівняння:

системи координат камери до планетоцентричної системи:

північного полюса обертання планети у поточний момент часу t; обчислюють за рекомендованими даними МАС, що уже відзначалося раніше;

космічного апарата на момент знімання; ці елементи визначають з формул руху небесних тіл, вони виникають під впливом реального гравітаційного поля, впливу

303

Сонця та планет атмосферного гальмування, магнітного та теплового полів, сонячної радіації;

Al0,All9Al2 - матриці, обчислені за кутовими елементами орієнтування координатної системи космічного апарату у геоекваторіальній системі на момент знімання;

Л13,Л14,Л15 - матриці, обчислені за кутовими елементами орієнтування камери у координатній системі космічного апарата;

А16 - матриця, що описує кутовий перехід θ від планетоцентричної орбітальної системи до планетоцентричної екваторіальної системи (формула (8.3),

де rs - віддаль від центра планети до космічного апарата (КА) (рис. 8.9);

визначається за параметрами орбіти на момент знімання.

Віддаль г як віддаль від космічного апарата до точки на поверхні планети обчислюється, якщо відомі дані про рельєф планети. Якщо ж таких даних немає, то беруть віддаль від КА до точки на референц-поверхні. У деяких випадках на КА може працювати лазерний біддалемір, що безпосередньо вимірює таку віддаль.

Використання рівнянь колінеарності (8.8) у космічній фотограмметрії таке саме, як і в аерофотограмметрії: вони становлять базу для космічної фотогріангуляції (див. п. 8.4).

8.3.5. Визначення планетоцентричних координат точок планет із використанням сканерних систем

Для отримання математичних залежностей необхідно брати до уваги геометрію формування зображення. На рис. 8.10 маємо:

кут сканування; г - віддаль від сканера до точки планети. Координати точки визначаються згідно з рис. 8.10 так:

X "0

У = rsin/?

ζ- r c o s / ?

304

Тому для загального випадку знімання сканерною системою з лінійною розгорткою з (8.14) маємо:

Рис. 8.10. Дія сканера з лінійною розгорткою

Віддаль rs до планети визначається так само, як і для рівнянь (8.14).

Сучасні сканерні знімальні системи по-різному орієнтовані під час знімання. Найпоширенішими є два способи, які, до речі, дають змогу отримати стереопари:

-системи з постійною орієнтацією центрального променя на центр мас планети;

-системи з постійним кутовим положенням сканера у просторі.

Упершому випадку головний промінь сканера розміщений у площині орбіти КА і скерований на центр мас планет, а площина сканування перпендикулярна до площини орбіти КА.

Координати точки планети обчислюють так:

305

(Нагадуємо, що θ - кут, названий аномалією орбіти; це один із параметрів, що характеризує орбіту КА; rs - довжина планетоцентричного радіус-вектора КА).

8.3.6. Визначення планетоцентричних координат точок планети з даних радіолокаційного знімання

В останні роки набуло широкого поширення раділокаційне знімання поверхні планети, зокрема Землі. При такому виді знімання є беззаперечні переваги, головна з яких - це незалежність від атмосферного стану (задимленість, хмарність). Радіопромені вільно проникають через атмосферу, а відбитий сигнал від поверхні дає інформацію про радіояскравість елементарної поверхні. Положення точки земної поверхні у просторі визначається віддаллю від радіолокатора до планети, складовою швидкості та елементами осколювальної орбіти на момент знімання.

При зніманні планети в такому випадку вісь ζ р скерована в центр мас

планети, а антена радіолокатора відхилена від площини орбіти КА на заданий кут ліворуч або праворуч.

З певних геометричних міркувань (рис. 8.11) та (8.14) можна отримати формули для обчислення просторових координат точки поверхні планети:

ζrcosap sin βρ

де кути а р , β р показано на рис. 8.11.

306

Рис. 8.11. Сферичні координати (орієнтація координатної системи)

при радіолокаційному зніманні

Взагалі кажучи, для створення коректної математичної моделі опрацювання радіолокаційних зображень необхідно знати геометричні засади формування цих образів. Наприклад, у фундаментальній монографії проф. Ю. Тюфліна [19] згадується, що під час радіолокального знімання поверхні Венери були задіяні два різні способи.

Перший, названий "смуговим зніманням", грунтувався на тому, що відбиті від поверхні планет сигнали у межах вузької поверхневої смуги створювали один рядок панорамного образу. Для кожного піксела смуги були відомі віддаль г , кут β р та момент часу знімання t. За цими даними можна обчислити координати

точки планети в планетоцентричній орбітальній системі.

У другому способі для кожного елемента поверхні планет фіксується час запізнення імпульсу τ (проходження радіоімпульсу до об'єкта і назад) та доплерівське зміщення частоти ν, які зв'язані з віддаллю г та величиною складової швидкості Vr по напрямку на точку поверхні:

307

де c - швидкість світла; λ - довжина хвилі.

Ці дані дають змогу визначити координати точки планети у планетоцентричній орбітальній системі координат.

Подамо приклад-схему знімання з використанням радара бокового бачення із синтезованою апертурою SAR (англ. Synthetc Aperture Radar) на рис. 8.12 а, б.

На рис. 8.12, а показана геометрична схема отримання зображення із супутника. Антена радара має розміри 10м χ їм, довшою стороною зорієнтована в напрямку польоту. Висилає вузькоспрямований пучок радіохвиль під кутом 23°. Роздільна здатність на Землі становить 30x30 м. Віддаль між витками (в проекції на поверхню Землі) становить 300 км, а ширина смуги становить 100 км. Це система безперервної дії.

На рис. 8.12, б показана схема т. зв. хвильового способу (SAR WM - від англ. SAR Wave Mode), який використовується для реєстрації невеликих територій

оцінки океанічних хвиль за допомогою вимірювання напрямку поширення та їхньої довжини.

308

8.4.Космічна фототріангуляція

8.4.1.Математичні моделі космічної фототріангуляції

Основним завданням космічної фототріангуляції є створення мережі опорних точок планети, а в низці випадків уточнення орбітальних та навігаційних параметрів руху космічного апарату.

Для побудови мереж найбільш придатні кадрові знімки. Проте це не виключає застосування сканерних та радіолокаційних знімальних систем.

Залежно від типу знімальної системи вхідними фотограмметричними величинами є координати точки на знімку та елементи внутрішнього орієнтування (для кадрових систем), координати точок на панорамному зображенні та масштабні коефіцієнти (для сканерних систем), орієнтація радіолокаційного вектора - променя - кути відхилення променя від планетоцентричного радіуса - вектора та віддаль до точки планети (для радіолокаційних систем).

Види опорних даних у космічній фототріангуляції доволі різноманітні. Частково про це згадувалось раніше, а зараз ще раз виокремимо найважливіші:

- координати центрів проекцій знімальної системи на момент знімання, визначені в планетоцентричній геоекваторіальній системі певної епохи, отримані з диференційних рівнянь руху космічного апарата;

-кутові елементи зовнішнього орієнтування зображень, визначені за допомогою знімків зоряного неба та при відомих параметрах взаємного розташування на космічному апараті знімальної системи та зоряної камери (що фотографує зоряне небо);

-дані радіопрофілювання поверхні планети з космічного апарату;

-дані віддалемірних вимірювань від КА до поверхні планети (лазерні або радіовіддалеміри);

-координати точок поверхні планети в планетоцентричній системі координат, зафіксованої відносно тіла планети.

Можливі різні комбінації з цих опорних даних, і все залежить від конкретної

ситуації: від постановки задачі, мети, якої треба досягти, наявного апаратного комплексу тощо.

Загальна математична модель достатньо добре описується рівняннями, поданими в розділі 2, проте в кожній конкретній ситуації є свої особливості. Покажемо це на прикладі фототріангуляції з одночасним визначенням початкових умов руху супутника.

В основу фототріануляції покладено умову колінеарності для точки поверхні планети Μ, точки зображення т та центру знімання S (рис. 8.9). Ця умова описується рівняннями (8.11). Якщо задача сформульована так, як вказано вище, то невідомими параметрами є:

309