start PT car legs fig

U-Shift · Nov 28, 2023 · 19aebc2 · 19aebc2
1 parent 0b8daf5
commit 19aebc2
Show file tree

Hide file tree

Showing 4 changed files with 32 additions and 0 deletions.
diff --git a/paper/PaperCEUS/PTcar_legs.R b/paper/PaperCEUS/PTcar_legs.R
@@ -0,0 +1,31 @@
+# see heat-prep_car.R to see where these data comes from
+
+carlegs_2_or = OD_2legs_lts3 |> st_as_sf(sf_column_name = "origins")
+carlegs_2_de = OD_2legs_lts3 |> st_as_sf(sf_column_name = "destinations")
+
+# mapview::mapview(carlegs_2)
+# id = 5123
+
+carlegs_2_od = st_coordinates(OD_2legs_lts3 |> st_as_sf(sf_column_name = "origins")) |> as.data.frame() |> 
+  cbind(st_coordinates(OD_2legs_lts3 |> st_as_sf(sf_column_name = "destinations")) |> as.data.frame()) |> 
+  cbind(OD_2legs_lts3$id)
+
+carlegs_2_od = stplanr::od_coords2line(carlegs_2_od)
+
+mapview::mapview(carlegs_2_od)
+id_route = 45563
+
+
+###### LEG 2 in PT
+mapview::mapview(carlegs_2_od[669,])
+
+
+
+carlegs_13_od = OD_13legs_lts4
+
+###### LEG 2 in Bike
+mapview::mapview(carlegs_13_od[carlegs_13_od$id_route == 45563,])
+
+
+# Try to run the routing model WITH grometry to get car routing
+
diff --git a/paper/PaperCEUS/PaperCEUS.Rmd b/paper/PaperCEUS/PaperCEUS.Rmd
@@ -297,6 +297,7 @@ In particular, for the urban train and tram -- with 100% electric traction -- on
 The conversion of avoided emissions into avoided welfare loss and respective monetary valuation was based on the EU Guide to Cost-benefit Analysis [@EuropeanCommission2014] and the best up-to-date reference values for the various gases [@EuropeanCommission2014; @bickel2006; @UNITE]: 8.44 €/ton for CO, 2,867.85 €/ton for NO~X~, 340,969.27 €/ton for PM10, 7,169.62 €/ton for VOC and 35.85 €/ton for CO~2~eq.
 The social impacts are in avoided premature mortality. This result is finally monetized using the *Statistical Value of Life* for Portugal: €3,055,358/fatality [@ANSR2021].
 We updated all the monetary reference values of the literature based on the annual inflation rate in Portugal for 2022[^7], and our 10-years estimations assumed a discount rate of 5% and inflation of 3%. 
+See [Research Data](#research-data) for all the input values we used.
 
 [^7]: See [Statistics Portugal tool for inflation rate estimates between years](https://www.ine.pt/xportal/xmain?xpid=INE&xpgid=ipc).
 

diff --git a/paper/PaperCEUS/PaperCEUS_files/figure-latex/mododist-1.pdf b/paper/PaperCEUS/PaperCEUS_files/figure-latex/mododist-1.pdf
diff --git a/paper/PaperCEUS/PaperCEUS_files/figure-latex/mododist-2.pdf b/paper/PaperCEUS/PaperCEUS_files/figure-latex/mododist-2.pdf