Skip to content

Add central transit accessibility calculation #375

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft
wants to merge 4 commits into
base: olusanya
Choose a base branch
from
Draft

Add central transit accessibility calculation #375

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
f8c73b7
Add central transit accessibility calculation
Sep 29, 2021
635dc2a
Add school tour cost parameter and save array
Oct 1, 2021
1e0dc01
Calculate weighted average over all home-based purposes
Oct 4, 2021
b89d763
Merge branch 'olusanya' into feat/malpakka
Oct 4, 2021
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
24 changes: 24 additions & 0 deletions Scripts/models/logit.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -457,6 +457,30 @@ def calc_accessibility(self, impedance):
Type (time/cost/dist) : numpy 2-d matrix
Impedances
"""
# Calculate cbd transit access
time_param = self.dest_choice_param["transit"]["impedance"]["time"]
try:
cost_param = self.dest_choice_param["transit"]["impedance"]["cost"]
except KeyError:
# School tours do not have a constant cost parameter
# Use value of time conversion from CBA guidelines instead
cost_param = -0.274445953
try:
time_access = time_param * impedance["transit"]["time"][:, 0]
except ValueError:
time_access = time_param[0] * impedance["transit"]["time"][:, 0]
try:
cost_access = cost_param * impedance["transit"]["cost"][:, 0]
money_utility = 1 / cost_param
except ValueError:
cost_access = cost_param[0] * impedance["transit"]["cost"][:, 0]
money_utility = 1 / cost_param[0]
central_access = pandas.Series(
time_access + cost_access, self.purpose.zone_numbers)
self.resultdata.print_data(
central_access, "central_transit_access.txt", self.purpose.name)
self.purpose.central_access = money_utility * central_access

mode_expsum = self._calc_utils(impedance)

# Calculate sustainable and car accessibility
Expand Down
4 changes: 4 additions & 0 deletions Scripts/modelsystem.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -360,6 +360,7 @@ def _calculate_accessibility_and_savu_zones(self):
logsum = 0
sust_logsum = 0
car_logsum = 0
central_access = 0
for purpose in self.dm.tour_purposes:
if (purpose.area == "metropolitan" and purpose.orig == "home"
and purpose.dest != "source"
Expand All @@ -370,6 +371,7 @@ def _calculate_accessibility_and_savu_zones(self):
logsum += weight * purpose.access
sust_logsum += weight * purpose.sustainable_access
car_logsum += weight * purpose.car_access
central_access += weight * purpose.central_access
pop = self.zdata_forecast["population"][bounds]
self.resultdata.print_line(
"\nTotal accessibility:\t{:1.2f}".format(
Expand All @@ -383,6 +385,8 @@ def _calculate_accessibility_and_savu_zones(self):
self.resultdata.print_data(
sust_logsum, "sustainable_accessibility.txt", "all")
self.resultdata.print_data(car_logsum, "car_accessibility.txt", "all")
self.resultdata.print_data(
central_access, "central_transit_access.txt", "all")
savu = numpy.searchsorted(zone_param.savu_intervals, sust_logsum) + 1
self.resultdata.print_data(
pandas.Series(savu, zone_numbers), "savu.txt", "savu_zone")
Expand Down