Elevation Profiles in R

callan.csv

callan.R

gps <- read.csv("callan.csv", 
  header = TRUE)

# calculate moving average: (this section is optional)
library(TTR)
movingN <- 5 # define the n for the moving average calculations
gps$Altitude <- gps$Altitude * 3.281 # convert m to ft
gps$SMA <- SMA(gps$Altitude, 
  n = movingN)
gps <- gps[movingN:length(gps$SMA), ] # remove first n-1 points

# calculate distances, if not included in your dataset:
library(sp)
Dist <- 0
for(i in 2:length(gps$Longitude)) {
  Dist[i] = spDistsN1(as.matrix(gps[i,c("Longitude", "Latitude")]),
    c(gps$Longitude[i-1], gps$Latitude[i-1]),
    longlat = TRUE) / 1.609 # longlat so distances will be in km, then divide to convert to miles
}
gps$Dist <- Dist

DistTotal <- 0
for(i in 2:length(gps$Longitude)) {
  DistTotal[i] = Dist[i] + DistTotal[i-1]
}
gps$DistTotal <- DistTotal


# plot:
library(ggplot2)
ggplot(gps, aes(x = DistTotal)) +
geom_ribbon(aes(ymin = 600, # change this to match your min below
    ymax = SMA),
  fill = "#1B9E77") + # put your altitude variable here if not using moving averages
labs(x = "Miles",
  y = "Elevation") +
scale_y_continuous(limits = c(600,1200)) # change this to limits appropriate for your region

callanMetric.R

gps <- read.csv("callan.csv", 
  header = TRUE)

# calculate moving average: (this section is optional)
library(TTR)
movingN <- 5 # define the n for the moving average calculations
gps$SMA <- SMA(gps$Altitude, 
  n = movingN)
gps <- gps[movingN:length(gps$SMA), ] # remove first n-1 points

# calculate distances, if not included in your dataset:
library(sp)
Dist <- 0
for(i in 2:length(gps$Longitude)) {
  Dist[i] = spDistsN1(as.matrix(gps[i,c("Longitude", "Latitude")]),
    c(gps$Longitude[i-1], gps$Latitude[i-1]),
    longlat = TRUE) # longlat so distances will be in km
}
gps$Dist <- Dist

DistTotal <- 0
for(i in 2:length(gps$Longitude)) {
  DistTotal[i] = Dist[i] + DistTotal[i-1]
}
gps$DistTotal <- DistTotal


# plot:
library(ggplot2)
ggplot(gps, aes(x = DistTotal)) +
geom_ribbon(aes(ymin = 180, # change this to match your min below
    ymax = SMA),
  fill = "#1B9E77") + # put your altitude variable here if not using moving averages
labs(x = "Distance, km",
  y = "Elevation, m") +
scale_y_continuous(limits = c(180,360)) # change this to limits appropriate for your region