In ahorawzy/HVS: This package is designed for Highway Virtual Traffic Survey Station project.

rm(list=ls())

library(HVS)
library(magrittr)
library(knitr)
library(plyr)

load("D:/R/packages/HVS/data-raw/hlj_roadnet.RData")
load("D:/R/packages/HVS/data-raw/indextable.RData")
load("D:/R/packages/HVS/data-raw/hljsf21.RData")

标准过程

result <- dijkstra_matrix(weight_mat = hlj_roadnet)
mileage <- result[[1]]
path <- result[[2]]

hljsf21$ENSEQ <- indextoseq(hljsf21$ENSTATIONINDEX)
hljsf21$EXSEQ <- indextoseq(hljsf21$EXSTATIONINDEX)
hljsf21 <- hljsf21[hljsf21$ENSEQ != 0 & hljsf21$EXSEQ != 0,]

hljsf21$path <- mapply(function(O,D) path[[O]][[D]],hljsf21$ENSEQ,hljsf21$EXSEQ)

hljsf21$link <- lapply(hljsf21$path,extract_link)

alllink <- all_link(hlj_roadnet)

library(parallel)
cores <- detectCores()
cluster <- makePSOCKcluster(cores)
l <- parLapply(cluster,alllink,caculate_natureflow,hljsf21$link)
names(l) <- alllink

d <- data.frame(traffic_volumn = simplify2array(l))
d$seqlink <- rownames(d)
d$indexlink <- seqlinktoindexlink(rownames(d))
rownames(d) <- 1:nrow(d)
d <- d[,c(3,2,1)]
kable(head(d,20))

system.time(caculate_natureflow_all(hljsf21))

x <- caculate_natureflow_all(hljsf21)
head(x,20)

分车型（VT）

system.time(caculate_natureflow_by(hljsf21,"VT"))

flowbyvt <- caculate_natureflow_by(hljsf21,"VT")

head(flowbyvt,20)

分车等级(VC)

system.time(caculate_natureflow_by(hljsf21,"VC"))

flowbyvc <- caculate_natureflow_by(hljsf21,"VC")

head(flowbyvc,20)

ahorawzy/HVS documentation built on May 29, 2019, 1:52 a.m.