R+OMX County Flow

CountyFlow.r

modelPath = "C:\\Modelrun\\Model80\\Base\\"

source("http://raw.githubusercontent.com/osPlanning/omx/dev/api/r/omx.R")

ZtoC = A1201[,c("TAZ","COUNTY")] # This is the TAZ-County equivalence table


hbwpk = readMatrixOMX(paste(modelPath,"HBWPK.OMX",sep=""),"M1")

cnames = rbind('HAMILTO', 'CLERMON', 'BUTLER', 'WARREN', 'BOONE', 'CAMPBEL', 'KENTON', 'DEARBOR') #could probably also do rbind(unique(ZtoC$COUNTY))

hbwpk.cf = matrix(0,length(cnames),length(cnames))
for(ic in 1:length(cnames)){
  for(jc in 1:length(cnames)){
    hbwpk.cf[ic,jc] = sum(hbwpk[which(seq(1:nrow(hbwpk)) %in% ZtoC[which(ZtoC$COUNTY == cnames[ic]),"TAZ"]),
                                 which(seq(1:ncol(hbwpk)) %in% ZtoC[which(ZtoC$COUNTY == cnames[jc]),"TAZ"])])
  }
}

# Optional
rownames(hbwpk.cf) = cnames
colnames(hbwpk.cf) = cnames