Merge b9b90a4 into 65582b2

.gitignore

@@ -10,6 +10,12 @@ test/media/
 # media files generated by scripts
 scripts/media/
 
+# ipynb checkpoints in pythonscripts
+scripts/pythonscripts/.ipynb_checkpoints
+
+# idm acc trace files stored in lmfit
+scripts/lmfit/
+
 # Files generated by invoking Julia with --code-coverage
 *.jl.cov
 *.jl.*.cov

scripts/experiments.jl

@@ -51,6 +51,25 @@ for i in 1:length(s)
         JLD.save("media/upper_$i.jld","m",m,"p",p,"md",md,"veh_id_list",veh_id_list,"ts",ts,"te",te)
 end
 
+#*****************Generate idm params for lmfit************
+s = scenarios_upper
+f = FilteringEnvironment()
+name = "upper"
+
+for i in 1:length(s)
+        veh_id_list = s[i][1]
+        ts = s[i][2][1]
+        te = s[i][2][2]
+        filename = "media/upper_$i.jld";
+        id_list,ts,te = JLD.load(filename,"veh_id_list","ts","te");
+        scene_list_true = replay_scenelist(f,id_list=id_list,ts=ts,te=te);
+
+        dirname = "$(name)_$(i)"
+        mkdir("lmfit/$(dirname)")
+        feat_dict = scenelist2idmfeatures(f,scene_list_true,id_list=id_list,
+                        scenario_name=name,scenario_number=i);
+end
+
 
 #************Docstring example code from multiscenarios_pf in helpers.jl*****
 # USEFUL to keep around for making tikz plots later on by rerunning this script
@@ -87,6 +106,12 @@ coll_mat_list = [coll_mat_idm,coll_mat_cidm,coll_mat_lmidm,coll_mat_pf];
 coll_barchart(coll_mat_list,filename = "media/coll_barchart_upper.svg");
 
 
+#***********Histogram plot****************
+# Full example in docstring of Filtering/metrics.jl vel_distribution
+a = PGFPlots.Axis([veh_hist_true,v_hist_idm,v_hist_cidm,v_hist_lmidm,v_hist_pf,
+Plots.Command(raw"\legend{true,idm,cidm,lmidm,pf}")
+]);
+
 #********************Train upper test lower******************
 # We need to show a variability in the generated scenarios
 # So we need to combine the particle sets of different vehicles together
@@ -113,4 +138,4 @@ coll_barchart(coll_mat_list,filename = "media/coll_barchart_upper.svg");
 #         push!(test_tdata,Dict("t"=>t))
 # end
 
-# curve_fit(test_lmfit,test_tdata,ydata,[0.5,0.5])
+# curve_fit(test_lmfit,test_tdata,ydata,[0.5,0.5])