flat2lla

Преобразование из плоского положения Земли в массив геодезических координат

Описание

пример

lla = flat2lla(flatearth_pos,llo,psio,href) оценивает массив геодезических координат, lla, из массива плоских координат Земли, flatearth_pos. Эта функция оценивает lla значение относительно ссылку местоположения, которое вы задаете с помощью llo, psio, и href.

lla = flat2lla(___, ellipsoidModel) оценивает координаты для определенной эллипсоидной планеты.

lla = flat2lla(___,flattening,equatorialRadius) оценивает координаты для пользовательской эллипсоидной планеты, заданные как flattening и equatorialRadius.

Примеры

свернуть все

Вычислите широту, долготу и высоту в одной координате:

lla = flat2lla( [ 4731 4511 120 ], [0 45], 5, -100)
lla =
    0.0391   45.0441  -20.0000

Вычислите широты, долготы и высоты в нескольких координатах с помощью WGS84 эллипсоидной модели:

lla = flat2lla( [ 4731 4511 120; 0 5074 4498 ], [0 45], 5, -100, 'WGS84' )
lla =
   1.0e+03 *

    0.0000    0.0450   -0.0200
   -0.0000    0.0450   -4.3980

Оцените широты, долготы и высоты в нескольких координатах с помощью пользовательской эллипсоидной модели:

f = 1/196.877360;
Re = 3397000;
lla = flat2lla( [ 4731 4511 120; 0 5074 4498 ], [0 45], 5, -100,  f, Re )
lla =
   1.0e+03 *

    0.0001    0.0451   -0.0200
   -0.0000    0.0451   -4.3980

Входные параметры

свернуть все

Плоские координаты положения Земли, заданные как 3-элементный вектор, в метрах.

Типы данных: double

Широта и долгота опорного местоположения, заданные как массив m-на-2 в степенях, для начала оценки и источника плоской системы координат Земли.

Типы данных: double

Угловое направление плоской оси x Земли, заданное как скаляр. Угловое направление - степени по часовой стрелке с севера, это угол в степенях, используемый для преобразования плоских x Земли и y координат в северные и восточные координаты.

Типы данных: double

Эталонная высота от поверхности Земли до плоской системы координат Земли относительно плоской системы координат Земли, заданная в виде скаляра, в метрах.

Типы данных: double

Модель эллипсоидной планеты. 'WGS84' является единственной опцией.

Типы данных: char | string

Уплощение планеты, заданное как двойной скаляр.

Типы данных: double

Планетарный экваториальный радиус, заданный как скаляр, в метрах.

Типы данных: double

Выходные аргументы

свернуть все

Геодезические координаты (широта, долгота и высота), возвращенные как m-на-3 массив, в [степенях, степенях, метрах].

Алгоритмы

Оценка начинается с преобразования плоских x Земли и y координат в северные и восточные координаты. Преобразование имеет форму

[NE]=[cosψsinψsinψcosψ][pxpy],

где(ψ¯) - угол в степенях по часовой стрелке между осью x и севером.

Для преобразования координат севера и востока в геодезическую широту и долготу в оценке используется радиус кривизны в основной вертикали (RN) и радиус кривизны в меридиане (RM). (RN) и (RM) определяются следующими отношениями:

RN=R1(2ff2)sin2μ0,

и

RM=RN1(2ff2)1(2ff2)sin2μ0,

где (R) - экваториальный радиус планеты и(f¯) - сплющивание планеты.

Небольшие изменения широты и долготы аппроксимируются из небольших изменений в положениях Север и Восток

dμ=atan(1RM)dNdι=atan(1RNcosμ)dE

и

dι=atan(1RNcosμ)dE.

Выход широта и долгота - это начальные широта и долгота плюс небольшие изменения широты и долготы.

μ=μ0+dμι=ι0+dι

Высота является отрицательным значением плоской Земли z -оси минус эталонная высота (href).

h=pzhref

Ссылки

[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.

См. также

Введенный в R2011a