У меня есть набор координат, которые я хочу повернуть в угол и выполнить анализ Anova. Для этого я использовал функцию maptools trackAzimuth, но я продолжаю терять одну или две точки. Я предполагаю, что он рассчитывает углы между двумя координатами. Но я хочу 1 точку (угол) на 1 широту и долготу, например на [-80,222, 30,555 до 45]. Поскольку моя долгота находится в Северной Америке, она имеет отрицательное значение, нужно ли преобразовывать ее в положительную, прежде чем преобразовать в угол? Пожалуйста, просветите меня об этом. Ниже приведены мои данные:
Longitude Latitude
-104.952 39.71478
-104.952 39.7149
-104.54 39.7148
-104.955 39.70441
-104.966 39.7175
Мои коды:
setwd("C:/Users/data")
install.packages("maptools")
library(maptools)
data <- read.table("data.csv", header = T, sep = ",")
dfr<-data.frame(data[3:4])
angle<-data.frame(trackAzimuth(as.matrix(dfr)))
Мои результаты: Угол -81,001 -95,57075 -175,254 -32,628 Здесь я потерял 1 широту и долготу. И мои углы отрицательные, я хочу положительные углы на этом. Как это сделать? Пожалуйста, помогите.
Спасибо