R/caculate_freight.R

In ahorawzy/Mreport: Monthly Report of traffic condition

#' @export
caculate_travel_volume <- function(section_volume,mileage,atts){
  l <- split(section_volume,section_volume[[atts]])
  m <- split(mileage,mileage[[atts]])
  l2 <- lapply(l, `[`, 3:6)
  m2 <- lapply(m, `[`, 2:5)
  k2 <- list(length = length(l2))
  for (i in 1:length(l2)) {
    k2[[i]] <- l2[[i]] * matrix(as.numeric(m2[[names(l2[i])]]),nrow = nrow(l2[[i]]),ncol = 4,byrow=T)
  }
  names(k2) <- names(l2)
  k3 <- lapply(k2, rowSums, na.rm = T)
  for(i in 1:length(k3)){
    names(k3[[i]]) <- l[[1]][[1]]
  }
  return(data.frame(k3))
}

#' @export
caculate_freight_trip <- function(section_volume,mileage,atts){
  l <- split(section_volume,section_volume[[atts]])
  m <- split(mileage,mileage[[atts]])
  l2 <- lapply(l, `[`, 3:6)
  m2 <- lapply(m, `[`, 2:5)
  k2 <- list(length = length(l2))
  for (i in 1:length(l2)) {
    k2[[i]] <- l2[[i]] * matrix(as.numeric(m2[[names(l2[i])]]),nrow = nrow(l2[[i]]),ncol = 4,byrow=T)
  }
  names(k2) <- names(l2)
  k3 <- k2
  for (i in 1:length(k3)) {
    k3[[i]][1] <- k2[[i]][1]/182
    k3[[i]][2] <- k2[[i]][2]/160
    k3[[i]][3] <- k3[[i]][3]/182
    k3[[i]][4] <- k3[[i]][4]/160
  }
  k4 <- lapply(k3, rowSums, na.rm = T)
  for(i in 1:length(k4)){
    names(k4[[i]]) <- l[[1]][[1]]
  }
  return(data.frame(k4))
}