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-на-3 в [степенях метров]. Значения широты и долготы могут быть любым значением. Однако значения широты + 90 и -90 могут вернуть неожиданные значения из-за особенности на полюсах.

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

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

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

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

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

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

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

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

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

Уплощение на каждом полюсе, заданное как скаляр.

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

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

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

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

свернуть все

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

Совет

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

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

Алгоритмы

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

dμ=μμ0dι=ιι0.

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

RN=R1(2ff2)sin2μ0,

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

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

dN=dμatan(1RM),

и

dE=dιatan(1RNcosμ0).

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

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

где

(ψ)

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

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

pz=hhref.

Ссылки

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

См. также

Темы

Введенный в R2011a