flat2lla

Преобразуйте от плоского Наземного положения до массива геодезической широты, долготы и высотных координат

Синтаксис

lla = flat2lla(flatearth_pos, llo, psio, href) lla = flat2lla(flatearth_pos, llo, psio, href, ellipsoidModel) lla = flat2lla(flatearth_pos, llo, psio, href, flattening, equatorialRadius)

Описание

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

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

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

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

`flatearth_pos`	Плоские Наземные координаты положения, в метрах.
`llo`	Ссылочное местоположение, в градусах, широты и долготы, для источника оценки и источника плоской Наземной системы координат.
`psio`	Угловое направление плоской Земли x - ось (степени по часовой стрелке от севера), который является углом, в градусах используемым для преобразования плоской Земли x и координаты y к Северным и Восточным координатам.
`href`	Ссылочная высота от поверхности Земли к плоской Земле структурирует относительно плоского Наземного кадра в метрах.
`ellipsoidModel`	Определенная модель планеты эллипсоида. Эта функция поддерживает только `'WGS84'`. Значение по умолчанию: WGS84
`flattening`	Пользовательская планета эллипсоида задана путем выравнивания.
`equatorialRadius`	Планетарный экваториальный радиус, в метрах.

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

lla

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

Примеры

Оцените широту, долготу и высоту в заданной координате:

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+003 *

    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+003 *

    0.0001    0.0451   -0.0200
   -0.0000    0.0451   -4.3980

Алгоритмы

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

$[\begin{matrix} N \\ E \end{matrix}] = [\begin{matrix} потому что ψ & - \sin ψ \\ \sin ψ & потому что ψ \end{matrix}] [\begin{matrix} p_{x} \\ p_{y} \end{matrix}]$

где $(\bar{ψ})$ угол в градусах по часовой стрелке между x - ось и севером.

Чтобы преобразовать Северные и Восточные координаты в геодезическую широту и долготу, оценка использует радиус искривления в главной вертикали (_RN) и радиус искривления в меридиане (_RM). (_RN) и (_RM) заданы следующими отношениями:

$\begin{array}{l} R_{N} = \frac{R}{\sqrt{1 - (2 f - f^{2}) \sin^{2} μ_{0}}} \\ R_{M} = R_{N} \frac{1 - (2 f - f^{2})}{1 - (2 f - f^{2}) \sin^{2} μ_{0}} \end{array}$

где (R) является экваториальным радиусом планеты и $(\bar{f})$ выравнивание планеты.

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

$\begin{array}{l} d μ = atan (\frac{1}{R_{M}}) d N \\ d ι = atan (\frac{1}{R_{N} потому что μ}) d E \end{array}$

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

$\begin{array}{l} μ = μ_{0} + d μ \\ ι = ι_{0} + d ι \end{array}$

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

$h = - p_{z} - h_{r e f}$

Ссылки

Etkin, B., динамика атмосферного рейса. Нью-Йорк: John Wiley & Sons, 1972.

Стивенс, B. L., и Ф. Л. Льюис, Управление Самолетом и Симуляция, 2-й редактор Нью-Йорк: John Wiley & Sons, 2003.

Смотрите также

lla2flat

Документация