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

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- 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 к Северным и Восточным координатам. Преобразование имеет форму

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

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

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

RN=R1(2ff2)sin2μ0RM=RN1(2ff2)1(2ff2)sin2μ0

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

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

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

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

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

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

h=pzhref

Ссылки

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

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

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

Введенный в R2011a
Для просмотра документации необходимо авторизоваться на сайте