exponenta event banner

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-by-2 в градусах, для начала оценки и начала плоской системы координат Земли.

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

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

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

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

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

Эллипсоидная модель планеты. 'WGS84' является единственным вариантом.

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

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

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

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

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

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

свернуть все

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

Алгоритмы

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

[NE] = [cosstart− sinü sinü cosstart] [pxpy],

где (δw) - угол в градусах по часовой стрелке между осью X и севером.

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

RN = R1 (2f f2) син2мк0,

и

RM = RN1 (2f f2) 1 − (2f f2) син2мк0,

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

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

= atan (1RM) dNd, = atan (1RNcosμ) dE

и

= атан (1RNcosμ) дЭ.

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

мкм = мк0 + мк

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

h = pz href

Ссылки

[1] Эткин, Б., Динамика атмосферного полета. Нью-Йорк: John Wiley & Sons, 1972.

[2] Стивенс, Б. Л., и Ф. Л. Льюис, Управление и моделирование летательных аппаратов, 2-й ред. Нью-Йорк: Джон Уайли и сыновья, 2003.

См. также

Представлен в R2011a