exponenta event banner

lla2flat

Перевести из геодезической широты, долготы и высоты в плоское положение Земли

Описание

пример

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.

Примеры

свернуть все

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

p = lla2flat( [ 0.1 44.95 1000 ], [0 45], 5, -100 )
p =
   1.0e+04 *

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

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

    0.5588   -0.3465   -0.0900
   -0.1373    1.7975   -0.1900

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

свернуть все

Геодезические координаты (широта, долгота и высота), заданные как массив m-by-3 в градусах метров. Значения широты и долготы могут быть любыми. Однако значения широты + 90 и -90 могут возвращать непредвиденные значения из-за сингулярности на полюсах .

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

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

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

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

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

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

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

Эллипсоидная модель планеты, указанная как 'WGS84'.

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

Сведение на каждом полюсе, указанное как скаляр.

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

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

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

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

свернуть все

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

Совет

  • Эта функция предполагает, что траектория полета и угол крена равны нулю.

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

Алгоритмы

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

dmi = λ мк0d,

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

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

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

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

dN = дмкатан (1RM),

и

dE = dtiatan (1RNcosмк0).

При преобразовании координат Север и Восток в плоские координаты Земля x и y преобразование имеет вид

[pxpy] = [cosü sinstart− singroupcos,] [NE],

где

(ψ)

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

Плоская ось Z Земли - это отрицательная высота минус опорная высота (href):

pz = h href.

Ссылки

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

См. также

Темы

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