Преобразование из плоского положения Земли в массив геодезических координат
оценивает массив геодезических координат, lla = flat2lla(flatearth_pos,llo,psio,href)lla, из массива плоских координат Земли, flatearth_pos. Эта функция оценивает lla значение относительно ссылку местоположения, которое вы задаете с помощью llo, psio, и href.
оценивает координаты для определенной эллипсоидной планеты.lla = flat2lla(___, ellipsoidModel)
оценивает координаты для пользовательской эллипсоидной планеты, заданные как lla = flat2lla(___,flattening,equatorialRadius)flattening и equatorialRadius.
Оценка начинается с преобразования плоских x Земли и y координат в северные и восточные координаты. Преобразование имеет форму
где - угол в степенях по часовой стрелке между осью x и севером.
Для преобразования координат севера и востока в геодезическую широту и долготу в оценке используется радиус кривизны в основной вертикали (RN) и радиус кривизны в меридиане (RM). (RN) и (RM) определяются следующими отношениями:
и
где (R) - экваториальный радиус планеты и - сплющивание планеты.
Небольшие изменения широты и долготы аппроксимируются из небольших изменений в положениях Север и Восток
и
Выход широта и долгота - это начальные широта и долгота плюс небольшие изменения широты и долготы.
Высота является отрицательным значением плоской Земли z -оси минус эталонная высота (href).
[1] Эткин, Б., Динамика атмосферного рейса. Нью-Йорк: John Wiley & Sons, 1972.
[2] Stevens, B. L., and F. L. Lewis, Aircraft Control and Simulation, 2nd ed. New York: John Wiley & Sons, 2003.