Влияние фотограмметрической и внутренней рефракций на определение линейных элементов внешнего ориентирования
В.И. Куштин, Н.Ф. Добрынин, Т.М. Пимшина
Ростовский государственный университет путей сообщения
Аннотация: Объектом исследований в представленной работе является влияние фотограмметрической и внутренней рефракций на определение линейных элементов внешнего ориентирования. Также предложен более простой по сравнению со строгим метод решения обратной пространственной фотограмметрической засечки, проверка правильности которого, выполненная методом моделирования, показала его достаточно высокую точность.
Ключевые слова: фотограмметрия, засечка, обратная пространственная засечка, рефракция, элементы внутреннего и внешнего ориентирования, опорные точки, метод моделирования.
Определение элементов внутреннего и внешнего ориентирования и поправок в координаты точек снимка по опорным точкам основано на решении обратной пространственной фотограмметрической засечки. Строгое решение обратной пространственной засечки по трем опорным точкам и их изображениям на снимке (рис.1) приводит к решению уравнения четвертой степени [1].
Исследование этого метода решения задачи с использованием моделей снимка и местности показало, что решение является довольно сложным и недостаточно эффективным для практического использования. Решение задачи осложняется тем, что проектирующие лучи вследствие оптической рефракции не являются прямыми, в общем случае они являются пространственными кривыми [2].
При выводе формул фотограмметрии предполагается, что световой луч от объекта до центра проекции является прямой линией [3-9]. В действительности луч света, проходя через слои атмосферы, имеющие различные плотности, искривляется. Это искривление называется рефракцией.
Б
3
Рис.1 - Обратная пространственная засечка по трем опорным точкам
Летательный аппарат, в котором находится приемник излучения, перемещаясь в атмосфере, создает вблизи себя поток обтекания, атмосферная среда в котором по своим оптическим характеристикам отличается от атмосферной среды в свободной атмосфере. Оптические лучи, проходя через поток обтекания, получают дополнительное искривление, называемое рефракцией в потоке обтекания. Затем оптические лучи проходят через иллюминатор и воздушную среду внутри летательного аппарата. Так как оптические характеристики иллюминатора и атмосферы внутри летательного аппарата отличаются от атмосферных условий в потоке обтекания, то возникает искривление и преломление оптических лучей в системе: поток обтекания - иллюминатор - атмосфера внутри летательного аппарата, называемое внутренней рефракцией [1].
Вопросы учета рефракции при решении обратной пространственной фотограмметрической засечки рассматривались в [1], но в этой работе было
1
сделано предположение о равенстве углов г/и гё (рис.1), что не соответствует действительности. Поэтому перед решением обратной пространственной фотограмметрической засечки целесообразно в координаты точек снимка или местности ввести поправки за фотограмметрическую рефракцию. Используем для этой цели достаточно строгие формулы, полученные в [10]:
дх = - г/
<( /2
х( /'2 + х2 + у2)
/ £2 х
р-/2 р-/2 = го у( /2 + х2 + у2) = г^
У ~ / ^ т2 ~ ^ г2
0 Р'0 р■/
где х, у - координаты точек на горизонтальном снимке,
/ - фокусное расстояние аэрофотоаппарата (АФА),
(1)
ГГ = С0
Я
(Р, - Ра )
- Н )
Ра
т
С0 =Р-10
- 6 ОД 0
Р0
(2)
(3)
где ЯС,,0 - удельная газовая постоянная и ускорение силы тяжести,
Р,, Ра - давление в снимаемой точке на земной поверхности и у АФА,
Ту $ = Т — 0,138—^ — виртуальная температура у АФА, еа - влажность воздуха у АФА,
На, Н,, - высоты центра проектирования Б и точки на земной поверхности, Ы0 - индекс преломления газовой среды (воздуха) при виртуальной температуре Ту0 и давлениир0, р= 206265''.
В тех случаях, когда оптический луч проходит последовательно через атмосферу, поток обтекания, иллюминатор и атмосферу кабины летательного аппарата, полагая, что границы перечисленных сред являются плоскими и параллельными друг другу (горизонтальными) [10]
г1 = г^ф
]к1
где
Г/к со'
Я ( Р, - Ра ) (- н&)
&
Рк_
т
(4)
где
с,,
N0 Т
Ро
р• 10
- 6
И0к - индекс преломления газовой среды в кабине при давлении р0 и температуре Т0. Если газовый состав в кабине летательного аппарата не отличается от газового состава внешней атмосферы, то Ск = С0.
В тех случаях, когда снимок является наклонным, то целесообразно поправки за рефракцию вводить в высоты точек на местности. На рис.2, учитывая = имеем
42,= = ^
= -Г г
3 3
81Пф
008 ф
2, + 42, - 23
= Г г
(X, - Х8)2 +( - )2 +(2, - 28)2 _ (2, +42, - 2з )2
(5)
Рис.2 - Поправки за рефракцию в высоты точек для наклонного снимка
Так как в формулу (5) входят определяемые величины ХБ, Уб, и Л2\, задача решается методом приближений.
Выясним, при каких углах наклона можно пользоваться формулами, полученными для горизонтальных снимков. Дифференцируя формулы (1) по переменным х, у, получим
dSx
P'
1 +
3x2 + y
2 Л
dx +
2 xy
Т2
dy
(6)
dSy
Р
1 +
x2 + 3 y
2 Л
2 xy dy + —г dx
Т
Максимальные значения йдх и йду будут при |х| = |у| в углах снимка. При г/0 = 60", |х| = |у| = 80 мм, / = 100 мм при дх =ду, получим |^дх| = \dSyl =0,0014|оХ|. При допустимом значении dдx = dдy = 0,010 мм имеем
сЬг = с1ху*2 = 10,10 мм, которые будут
;;.y-= / .1т мм, а
0JDS14
соответствовать максимальному искажению радиус-вектора на снимке
Sr
r а0
р , откуда r 2 113
Следовательно, при углах наклона а < 4,50 ошибки в Sx, Sy, вычисленные по формулам (1) при отмеченных выше значениях, не будут превышать 10 мкм. Поэтому формулами (1) можно пользоваться не только для горизонтальных, но и для плановых снимков.
Как отмечено выше, строгое решение обратной пространственной засечки является неэффективным, поэтому целесообразно разработать более простой метод решения этой задачи.
На рис.2 по теореме косинусов в треугольнике S12 имеем
S2 + S22 - 2SS cose - А2 = 0
или
а =
SrmJp 10.10 • 100 • 57.3
= 4.50
(7)
r
0
r
1
(8)
((1 - X,)2 + (( - 7,)2 + (71 - )2 + (X2 - X, )2 + (72 - 7,)2 +
+ (72 - 7,)2 - 2{[(( - X,)2 +(71 - 7,)2 + (71 - 7,)2] х
1
х [(( - X,)2 + (72 - 7,)2 + (72 - 7,)2 ]}* созД - Д2.
Примем Х3 = Ух= У/ + 23 = Щ + Подставив эти
значения в (8), после преобразований, отбрасывая члены второго и более высокого порядка малости, для трех точек имеем систему уравнений
+ МБ + С + 1г = 0 ,
+ + ¡2 = 0, ■+ Ьлд¥в ■+ ■+ ;а - О,
где
, = (( - X,0) + ((+1 - X,0)', Ьг =(7, - 7,)а + (7+1 - 7,)а', = (7, - 7,°) + (7,+1 - 7,° )£', /, — 1 (Д2 - Д2) = 1 (( - Д)((, + Д0),
(9)
а
0 — 1 - д, созД , 0 — 1 -
Ч,
,0
— +1 , 4 =
, = VX2 + У? + /2 , созв
Х X,+1 + У У+1 + /2
^ 2 + У2 + /2)(21+У+1 + /2)
Значения Б, определяют по формуле (7), в которую вместо , подставляют Задача решается методом приближений.
Проверка правильности полученных формул, выполненная методом моделирования по данным, приведенным в табл. 1, показала, что уже второе приближение позволяет определить значения X,, 7,, 7, с ошибками 0,135 м; 0,032 м и 0,502 м соответственно, т.е. с достаточно высокой точностью.
Для количественной оценки влияния рефракции на результаты решения обратной пространственной засечки предварительно в координаты точек снимка введем поправки 8х, 8у по формуле (1).
Таблица 1
Координаты точек на снимке и соответствующие им координаты точек на местности. Х$ = Уб = 0, Дб = 353815.0 м, / = 140 мм, X/ = 14700.00 м, У/ = -9819.35 м, 2° = 348319.00 м
№ точки х, мм у, мм X, м У, мм 2„ мм
1 73,73582 82,90761 196229,74 199939,31 -368,83
2 -89,69884 97,87368 -203754,14 203708,18 -618,38
3 -94,10511 -69,20215 -208153,80 -203195,47 -452,45
Положим, что фотографирование выполнено из кабины летательного аппарата, в которой давление рк = 1020 гПа, Тк = 290,0 К. Давление pg = 980 гПа, рa = 0. Для зеленой части спектра (X = 0,55 мкм) имеем С0 = 16.297 угл.с
я
/ гПа, С0— = 0.47702 угл.с / гПа. Газовый состав в кабине принят
,0
соответствующим газовому составу атмосферы. Результаты вычислений и сравнение полученных с учетом рефракций координат центра проектирования с их значениями, вычисленными без учета рефракции, показывают, что разности ЛХБ = - 4,36 мм, ЛУб = - 0,78 м, Л2б = 158,88 м.
Следовательно, фотограмметрическая и внутренняя рефракция в основном приводит к искажению координаты или высоты
фотографирования, тогда как искажение плановых координат является сравнительно небольшим.
Литература
1. Куштин И.Ф., Бруевич П.Н., Лысков Г.А. Справочник техника -фотограмметриста. - М.: Недра, 1988., 320 с.
2. Куштин В.И. Разработка и исследование методов аналитического трансформирования снимков и их использование при решении научно-технических задач. Диссертация на соискание степени к.т.н. - Ростов-на-Дону.: РГСУ, 1999, 226 с.
3. Dengsheng Lu, Emilio Moran, Scott Hetrick. Detection of impervious surface change with multitemporal Landsat images in an urban- rural frontier // ISRRS J. Photogramm. and Retome Sens, 2011. 66, № 3 pp. 298-306.
4. Balci Murat, Foroosh Hassan. Subpixel estimation of shifts directly in the Fourier domain //IEEE Trans. Image Process, 2006. 15, - № 7. - pp. 1965-1972.
5. Brandstatter G. Fundamentals of Algebro-Projective Photogrammetry URL: planet-austria.at/0xc1aa500d_0x00022343.pdf.
6. Куштин В.И., Ревякин А.А., Добрынин Н.Ф., Щербак П.Н., Петренко А.М. Разработка аналитических методов трансформирования и их применение для обработки космических снимков // Инженерный вестник Дона, 2017, №4 URL: ivdon.ru/ru/magazine/archive/n4y2017/4539.
7. Добрынин Н.Ф., Пимшина Т.М. Взаимное ориентирование снимков с новым сочетанием угловых элементов в стереопаре // Инженерный вестник Дона, 2014, №2 URL: ivdon.ru/ru/magazine/archive/n2y2014/2347.
8. Добрынин Н.Ф., Пимшина Т.М. Использование космических средств позиционирования при обработке аэро- и космической информации Инженерный вестник Дона, 2013, №3 URL: ivdon.ru/magazine/archive/n3y2013/1835.
9. Бирюков B.C. Обработка цифровых снимков в фотограмметрии: Монография / Москва: ВИУ, 2001. - 194 с.
10. Куштин И.Ф. Справочник техника - геодезиста. - М.: Недра, 1993., 320 с.
References
1. Kushtin I.F., Bruevich P.N., Lyskov G.A. Spravochnik tehnika -fotogrammetrista [Reference book of the photogrammetric technician]. M. Nedra,
1988., 320 p.
2. Kushtin V.I. Razrabotka i issledovanie metodov analiticheskogo transformirovanija snimkov i ih ispol'zovanie pri reshenii nauchno-tehnicheskih zadach. Dissertacija na soiskanie stepeni k.t.n. [Development and research of methods of analytical transformation of images and their use in solving scientific and technical problems. Thesis for the degree of Candidate of Technical Sciences.] Rostov-na-Donu. RGSU, 1999, 226 p.
3. Lu Dengsheng, Emilio Moran, Scott Hetrick. ISRRS J. Photogramm. and Retome Sens, 2011. 66, № 3. pp. 298-306.
4. Balci Murat, Foroosh Hassan. IEEE Trans. Image Process, 2006. 15, № 7. pp. 1965-1972.
5. Brandstatter G. Fundamentals of Algebro-Projective Photogrammetry. URL: planet-austria.at70xc1aa500d_0x00022343.pdf.
6. Kushtin V. I., Revyakin A. A., Dobrynin N. F., Shcherbak P. N., Petrenko A. M. Inzenernyj vestnik Dona (Rus), 2017, №4 URL: ivdon.ru/ru/magazine/archive/n4y2017/4539.
7. Dobrynin N.F., Pimshina T.M. Inzenernyj vestnik Dona (Rus), 2014, №2 URL: ivdon.ru/ru/magazine/archive/n2y2014/2347.
8. Dobrynin N.F., Pimshina T.M. Inzenernyj vestnik Dona (Rus), 2013, №3 URL: ivdon.ru/magazine/archive/n3y2013/1835.
9. Birjukov B.C. Obrabotka cifrovyh snimkov v fotogrammetrii. [Processing digital images in photogrammetry]. Monografija. Moskva: VIU, 2001. 194 p.
10. Kushtin I.F. Spravochnik tehnika - geodezista [Handbook of a surveyor]. M.
Nedra, 1993, 320 p.