lla2flat

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

Синтаксис

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

Описание

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

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

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

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

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

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

flatearth_pos

Плоские Наземные координаты положения, в метрах.

Примеры

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

p = lla2flat( [ 0.1 44.95 1000 ], [0 45], 5, -100 )

p =

  1.0e+004 *

    1.0530   -0.6509   -0.0900

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

p = lla2flat( [ 0.1 44.95 1000; -0.05 45.3 2000 ], [0 45], 5, -100, 'WGS84' )

p =

  1.0e+004 *

    1.0530   -0.6509   -0.0900
   -0.2597    3.3751   -0.1900

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

f = 1/196.877360;
Re = 3397000;
p = lla2flat( [ 0.1 44.95 1000; -0.05 45.3 2000 ], [0 45], 5, -100,  f, Re )

p =

  1.0e+004 *

    0.5588   -0.3465   -0.0900
   -0.1373    1.7975   -0.1900

Советы

Эта функция принимает, что курс полета и угол банка являются нулем.
Эта функция принимает, что плоская Земля z - ось нормальна к Земле только в начальной геодезической широте и долготе. Эта функция имеет более высокую точность по маленьким расстояниям от начальной геодезической широты и долготы. Это также имеет более высокую точность на расстояниях ближе до экватора. Функция вычисляет долготу с более высокой точностью, когда изменения в широте меньше. Кроме того, долгота сингулярна в полюсах.

Алгоритмы

Оценка начинается путем нахождения небольших изменений в широте и долготе от выходной широты и долготе минус начальная широта и долготе.

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

Чтобы преобразовать геодезическую широту и долготу к Северным и Восточным координатам, оценка использует радиус искривления в главной вертикали (_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) является экваториальным радиусом планеты и $f$ выравнивание планеты.

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

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

С преобразованием Северных и Восточных координат к плоской Земле x и координаты y, преобразование имеет форму

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

где

$(ψ)$

угол в градусах по часовой стрелке между x - ось и севером.

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

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

Ссылки

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

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

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

flat2lla

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