Преобразуйте азимут на эллипсоиде в направление на карте
th = vfwdtran(lat,lon,az)
th = vfwdtran(mstruct,lat,lon,az)
[th,len] = vfwdtran(...)
th = vfwdtran(lat,lon,az)
преобразует угол азимута в заданных точках широты и долготы на сфере в проекционное пространство. Отображаемая в данный момент проекция карты используется для определения проекционного пространства. Углы входа должны быть в тех же модулях, что и заданные текущей проекцией карты. Входами могут быть скаляры или матрицы равного размера. Угол в проекционном пространстве определяется как положительный против часовой стрелки от оси X.
th = vfwdtran(mstruct,lat,lon,az)
использует проекцию карты, заданный вход mstruct
для вычисления проекции карты.
[th,len] = vfwdtran(...)
также возвращает длину вектора в проективной системе координат. Значение 1 указывает на отсутствие искажений шкалы.
Направление севера легко определить на трехмерной сфере, но сложнее на двумерной карте. Для цилиндрических проекций в нормальном аспекте север всегда в положительном направлении y. Для конических проекций север может быть слева или справа от оси Y. Эта функция преобразует любой угол азимута в сфере в соответствующий угол в проективных координатах бумаги.
Выборочные расчеты:
axesm('eqdconic','maplatlim',[-10 45],'maplonlim',[-55 55]) gridm; framem; mlabel; plabel quiverm([0 0 0],[-45 0 45],[0 0 0],[10 10 10],0) quiverm([0 0 0],[-45 0 45],[10 10 10],[0 0 0],0)
vfwdtran([0 0 0],[-45 0 45],[0 0 0]) ans = 59.614 90 120.39 vfwdtran([0 0 0],[-45 0 45],[90 90 90]) ans = -30.385 0.0001931 30.386
Это преобразование ограничено областью, заданной пределами системы координат в текущем определении карты.
Измеряют угол географического азимута по часовой стрелке с севера. Измеряют угол проекционного пространства против часовой стрелки от оси X.
Эта функция использует метод конечного различия. Географические координаты слегка возмущаются в разных направлениях и проецируются. Небольшая величина ошибки введена численным расчетом производных и вариацией параметров искажения карты.