def hmm_graph():
global hmm_graph_
if hmm_graph_ is None:
if not os.path.exists(hmm_graph_fname):
if graph_type == 'simple':
hmm_graph_ = model.createHMMGraphFromNetwork(net, mode_counts=mode_counts())
else:
# Complex model not implemented
assert False
else:
tic("Reading completed hmm graph from %s"%hmm_graph_fname)
hmm_graph_ = pickle.load(open(hmm_graph_fname,'r'))
return hmm_graph_
def hmm_graph():
global hmm_graph_
if hmm_graph_ is None:
if not os.path.exists(hmm_graph_fname):
if graph_type == 'simple':
tic("creating empty hmm graph", experiment_name)
hmm_graph_ = model.createHMMGraphFromNetwork(net, mode_counts=mode_counts())
tic("done creating empty hmm graph", experiment_name)
tic("saving hmm graph pickle", experiment_name)
pickle.dump(hmm_graph(),open(hmm_graph_fname,'w'))
tic("done saving hmm graph pickle", experiment_name)
else:
# Complex model not implemented
assert False
else:
tic("Reading completed hmm graph from %s"%hmm_graph_fname)
hmm_graph_ = pickle.load(open(hmm_graph_fname,'r'))
tic("done reading completed hmm graph from %s"%hmm_graph_fname)
return hmm_graph_
traj_conv_param = experiment_design['trajectory_conversion']['params']
traj_conv = createTrajectoryConversion(graph_type=graph_type,
process=experiment_design['trajectory_conversion']['process'],
params=traj_conv_param,
network=net,
max_nb_mixture=traj_conv_param['max_n_modes'])
traj_conv_one_mode = createTrajectoryConversion(graph_type=graph_type,
process='mixture_auto',
params=traj_conv_param,
network=net,
max_nb_mixture=1)
# mode_counts = dict([(link_id,1) for link_id in net.keys()])
if graph_type == 'simple':
hmm_graph = model.createHMMGraphFromNetwork(net, mode_counts=traj_conv.modeCounts())
hmm_graph_one_mode = model.createHMMGraphFromNetwork(net, mode_counts=traj_conv_one_mode.modeCounts())
# hmm_graph = model.createHMMGraphFromNetwork(net, mode_counts=mode_counts)
else:
# Complex model not implemented
assert False
tt_graph = createTravelTimeGraph(hmm_graph, radius=2e-4)
tt_graph.checkInvariants()
tt_graph_one_mode = createTravelTimeGraph(hmm_graph_one_mode, radius=2e-4)
tt_graph_one_mode.checkInvariants()
gmrf = emptyValues(tt_graph)
gmrf_one_mode = emptyValues(tt_graph)
