vinvtran
Азимут на эллипсоиде от направляющего угла в плоскости карты
Синтаксис
az = vinvtran(x,y,th)
az = vinvtran(mstruct,x,y,th)
[az,len] = vinvtran(...)
Описание
az = vinvtran(x,y,th) преобразовывает угол в пробел проекции в точке, заданной x и y в угол азимута в географических координатах. Проекция карты, в настоящее время отображенная, используется, чтобы задать пробел проекции. Входные углы должны быть в тех же модулях, как задано текущей проекцией карты. Входные параметры могут быть скалярами или матрицами равного размера. Угол в углу пробела проекции th задан как положительный против часовой стрелки от оси X.
az = vinvtran(mstruct,x,y,th) использует проекцию карты, заданную входом struct, чтобы вычислить проекцию карты.
[az,len] = vinvtran(...) также возвращает длину вектора в географической системе координат. Значение 1 не указывает ни на какое искажение шкалы для того угла.
Фон
В то время как векторы вдоль оси Y всегда указывают на север в цилиндрической проекции в нормальном аспекте, они могут указать к востоку или к западу от севера на теории конических сечений, azimuthals, и других проекций. Эта функция вычисляет географический азимут для углов на спроектированном пробеле.
Примеры
Демонстрационные вычисления:
axesm('eqdconic','maplatlim',[-10 45],'maplonlim',[-55 55])
gridm; framem; mlabel; plabel
[x,y] = mfwdtran([0 0 0],[-45 0 45]);
quiver(x,y,[ .2 .2 .2],[0 0 0],0)
quiver(x,y,[0 0 0],[ .2 .2 .2],0)
vinvtran(x,y,[ 0 0 0])
ans =
57.345 90.338 124.98
vinvtran(x,y,[ 90 90 90])
ans =
331.99 0 28.008Ограничения
Это преобразование ограничивается областью, заданной пределами кадра в текущем определении карты.
Советы
Географический угол азимута измеряется по часовой стрелке от севера. Угол пробела проекции измеряется против часовой стрелки от оси X.
Эта функция использует метод конечной разности. Географические координаты встревожены немного в различных направлениях и спроектированы. Небольшое количество ошибки введено численным расчетом производных и изменением параметров искажения карты.