Функция конвертации геодезических координат в планетоцентрированные планетные координаты
p = llatoecef(lla); p = llatoecef(lla, f, Re);
lla - массив из трех вещественных чисел, координат некоторой точки в системе координат LLA (широта [град], долгота[град], высота[м]); либо матрица, сформированная несколькими точками построчно;
f - вещественное число, уплощение эллипса планеты;
Re - вещественное число, экваториальный радиус планеты;llatoecef(lla) – функция производит пересчет геодезических координат, заданных в виде массива lla широты [град], долготы [град] и высоты [м] в планетоцентрированные Земные координаты в виде (x,y,z)[м] с использованием модели Земли WGS84;
llatoecef(lla, f, Re) - функция производит пересчет геодезических координат, заданных в виде массива lla широты [град], долготы [град] и высоты [м] в планетоцентрированные планетные координаты в виде (x,y,z)[м] с использованием модели планеты, заданной экваториальным радиусом Re и уплощением f.
Если на вход подается одна точка в системе координат LLA (вектор из трех вещественных чисел), то функция возвращает одну точку (вектор из трех значений) в системе координат ECEF.
Если на вход подается несколько точек в системе координат LLA (матрица, где каждая строка это точка), то функция возвращает соответствующее число точек (матрица, где каждая строка это точка) в системе координа ECEF.
lla1 = [45 90 2000]; // точка 1 lla2 = [0 45 1000]; // точка 2 f = 1/196.877360; // уплощение некоторой планеты Re = 3397000; // экв. радиус некоторой планеты p1 = llatoecef(lla1); // земноцентрированные координаты точки1 (ECEF)) p2 = llatoecef(lla1, f, Re); // планетоцентрированные координаты точки1 p3 = llatoecef([lla1,lla2]); // земноцентрированные координаты точек 1 и 2 (ECEF)) p4 = llatoecef([lla1,lla2],f,Re); // планетоцентрированные координаты точек 1 и 2 (ECEF))
В примере производится пересчет координат из системы LLA в планетоцентрированные планетные координаты Земли и планеты с заданными характеристиками эллипса.