def get_poly_avg(input_poly_fc):
# as of 11/26/2019, each of these outputs are dictionaries
accdata = acc.get_acc_data(input_poly_fc, p.accdata_fc, p.ptype_area_agg, get_ej=False)
collision_data = coll.get_collision_data(input_poly_fc, p.ptype_area_agg, p.collisions_fc, 0)
mix_data = mixidx.get_mix_idx(p.parcel_pt_fc, input_poly_fc, p.ptype_area_agg)
intsecn_dens = intsxn.intersection_density(input_poly_fc, p.intersections_base_fc, p.ptype_area_agg)
bikeway_covg = bufnet.get_bikeway_mileage_share(input_poly_fc, p.ptype_area_agg)
tran_stop_density = trn_svc.transit_svc_density(input_poly_fc, p.trn_svc_fc, p.ptype_area_agg)
emp_ind_wtot = lubuff.point_sum(p.parcel_pt_fc, input_poly_fc, p.ptype_area_agg, [p.col_empind, p.col_emptot], 0)
emp_ind_pct = {'emp_ind_pct': emp_ind_wtot[p.col_empind] / emp_ind_wtot[p.col_emptot]}
pop_x_ej = lubuff.point_sum(p.parcel_pt_fc, input_poly_fc, p.ptype_area_agg, [p.col_pop_ilut], 0, p.col_ej_ind)
pop_tot = sum(pop_x_ej.values())
pct_pop_ej = {'pct_ej_pop': pop_x_ej[1] / pop_tot}
job_pop_dens = lubuff.point_sum_density(p.parcel_pt_fc, input_poly_fc, p.ptype_area_agg, \
[p.col_du, p.col_emptot], 0)
total_dens = {"job_du_dens_ac": sum(job_pop_dens.values())}
out_dict = {}
for d in [accdata, collision_data, mix_data, intsecn_dens, bikeway_covg, tran_stop_density, pct_pop_ej,\
emp_ind_pct, total_dens]:
out_dict.update(d)
return out_dict
def complete_streets_idx(fc_pclpt, fc_project, project_type, posted_speedlim,
transit_event_fc):
arcpy.AddMessage("getting complete street score...")
# get transit service density around project
tran_stops_dict = ts.transit_svc_density(fc_project, transit_event_fc,
project_type)
transit_svc_density = list(tran_stops_dict.values())[0]
lu_fac_cols = [p.col_area_ac, p.col_k12_enr, p.col_emptot, p.col_du]
lu_vals_cols = [p.col_k12_enr, p.col_emptot, p.col_du]
#get sums of the lu_fac_cols within project buffer area
lu_vals_dict = lubuff.point_sum(fc_pclpt, fc_project, project_type,
lu_fac_cols, p.cs_buffdist)
#dens_score = (student_dens + trn_svc_dens + job_dens + du_dens)
dens_score = sum(
[lu_vals_dict[i] / lu_vals_dict[p.col_area_ac]
for i in lu_vals_cols]) + transit_svc_density
csi = dens_score * (
1 - (posted_speedlim - p.cs_threshold_speed) * p.cs_spd_pen_fac)
out_dict = {'complete_strt_idx': csi}
return out_dict
def get_poly_avg(input_poly_fc):
# as of 11/26/2019, each of these outputs are dictionaries
pcl_pt_data = params.parcel_pt_fc_yr()
accdata = acc.get_acc_data(input_poly_fc, params.accdata_fc, params.ptype_area_agg, get_ej=False)
collision_data = coll.get_collision_data(input_poly_fc, params.ptype_area_agg, params.collisions_fc, 0)
mix_data = mixidx.get_mix_idx(pcl_pt_data, input_poly_fc, params.ptype_area_agg)
intsecn_dens = intsxn.intersection_density(input_poly_fc, params.intersections_base_fc, params.ptype_area_agg)
bikeway_covg = bufnet.get_bikeway_mileage_share(input_poly_fc, params.ptype_area_agg)
tran_stop_density = trn_svc.transit_svc_density(input_poly_fc, params.trn_svc_fc, params.ptype_area_agg)
emp_ind_wtot = lubuff.point_sum(pcl_pt_data, input_poly_fc, params.ptype_area_agg, [params.col_empind, params.col_emptot], 0)
emp_ind_pct = {'EMPIND_jobshare': emp_ind_wtot[params.col_empind] / emp_ind_wtot[params.col_emptot] \
if emp_ind_wtot[params.col_emptot] > 0 else 0}
pop_x_ej = lubuff.point_sum(pcl_pt_data, input_poly_fc, params.ptype_area_agg, [params.col_pop_ilut], 0, params.col_ej_ind)
pop_tot = sum(pop_x_ej.values())
key_yes_ej = max(list(pop_x_ej.keys()))
pct_pop_ej = {'Pct_PopEJArea': pop_x_ej[key_yes_ej] / pop_tot if pop_tot > 0 else 0}
job_pop_dens = lubuff.point_sum_density(pcl_pt_data, input_poly_fc, params.ptype_area_agg, \
[params.col_du, params.col_emptot], 0)
# total_dens = {"job_du_perNetAcre": sum(job_pop_dens.values())}
out_dict = {}
for d in [accdata, collision_data, mix_data, intsecn_dens, bikeway_covg, tran_stop_density, pct_pop_ej,\
emp_ind_pct, job_pop_dens]:
out_dict.update(d)
return out_dict
def get_singleyr_data(fc_project, projtyp, adt, posted_speedlim, out_dict={}):
pcl_pt_fc = p.parcel_pt_fc_yr(2016)
pcl_poly_fc = p.parcel_poly_fc_yr(2016)
accdata = acc.get_acc_data(fc_project, p.accdata_fc, projtyp, get_ej=False)
collision_data = coll.get_collision_data(fc_project, projtyp, p.collisions_fc, adt)
complete_street_score = {'complete_street_score': -1} if projtyp == p.ptype_fwy else \
cs.complete_streets_idx(pcl_pt_fc, fc_project, projtyp, posted_speedlim, p.trn_svc_fc)
truck_route_pct = {'pct_proj_STAATruckRoutes': 1} if projtyp == p.ptype_fwy else \
linex.get_line_overlap(fc_project, p.freight_route_fc, p.freight_route_fc) # all freeways are STAA truck routes
ag_acres = luac.get_lutype_acreage(fc_project, projtyp, pcl_poly_fc, p.lutype_ag)
pct_adt_truck = {"pct_truck_aadt": -1} if projtyp != p.ptype_fwy else truck_fwy.get_tmc_truck_data(fc_project, projtyp)
intersxn_data = intsxn.intersection_density(fc_project, p.intersections_base_fc, projtyp)
npmrds_data = npmrds.get_npmrds_data(fc_project, projtyp)
transit_data = trnsvc.transit_svc_density(fc_project, p.trn_svc_fc, projtyp)
bikeway_data = bnmi.get_bikeway_mileage_share(fc_project, p.ptype_sgr)
infill_status = urbn.projarea_infill_status(fc_project, p.comm_types_fc)
# total job + du density (base year only, for state-of-good-repair proj eval only)
job_du_dens = lu_pt_buff.point_sum_density(pcl_pt_fc, fc_project, projtyp,
[p.col_emptot, p.col_du], p.ilut_sum_buffdist)
comb_du_dens = sum(list(job_du_dens.values()))
job_du_dens['job_du_perNetAcre'] = comb_du_dens
# get EJ data
ej_data = lu_pt_buff.point_sum(pcl_pt_fc, fc_project, projtyp, [p.col_pop_ilut],
p.ilut_sum_buffdist, p.col_ej_ind, case_excs_list=[])
ej_flag_dict = {0: "Pop_NonEJArea", 1: "Pop_EJArea"} # rename keys from 0/1 to more human-readable names
ej_data = utils.rename_dict_keys(ej_data, ej_flag_dict)
ej_data["Pct_PopEJArea"] = ej_data["Pop_EJArea"] / sum(list(ej_data.values()))
accdata_ej = acc.get_acc_data(fc_project, p.accdata_fc, projtyp, get_ej=True) # EJ accessibility data
ej_data.update(accdata_ej)
# for base dict, add items that only have a base year value (no future year values)
for d in [accdata, collision_data, complete_street_score, truck_route_pct, pct_adt_truck, ag_acres, intersxn_data,
npmrds_data, transit_data, bikeway_data, infill_status, job_du_dens, ej_data]:
out_dict_base.update(d)
outdf = pd.DataFrame.from_dict(out_dict_base, orient='index')
return outdf
def complete_streets_idx(fc_pclpt, fc_project, project_type, posted_speedlim,
transit_event_fc):
'''Calculate complete street index (CSI) for project
CSI = (students/acre + daily transit vehicle stops/acre + BY jobs/acre + BY du/acre) * (1-(posted speed limit - threshold speed limit)*speed penalty factor)
'''
# col_area_ac
# col_k12_enr
# col_emptot
# col_du
# cs_buffdist
# cs_threshold_speed
# cs_spd_pen_fac
# don't give complete street score for freeway projects or if sponsor didn't enter speed limit
if project_type == params.ptype_fwy or posted_speedlim <= 1:
csi = -1
else:
# arcpy.AddMessage("Calculating complete street score...")
# get transit service density around project
tran_stops_dict = ts.transit_svc_density(fc_project, transit_event_fc,
project_type)
transit_svc_density = list(tran_stops_dict.values())[0]
lu_fac_cols = [
params.col_area_ac, params.col_k12_enr, params.col_emptot,
params.col_du
]
lu_vals_cols = [params.col_k12_enr, params.col_emptot, params.col_du]
# get sums of the lu_fac_cols within project buffer area
lu_vals_dict = lubuff.point_sum(fc_pclpt, fc_project, project_type,
lu_fac_cols, params.cs_buffdist)
print(lu_vals_dict)
#dens_score = (student_dens + trn_svc_dens + job_dens + du_dens)
if lu_vals_dict[params.col_area_ac] == 0:
dens_score = 0
else:
dens_score = sum([
lu_vals_dict[i] / lu_vals_dict[params.col_area_ac]
for i in lu_vals_cols
]) + transit_svc_density
csi = dens_score * (1 - (posted_speedlim - params.cs_threshold_speed) *
params.cs_spd_pen_fac)
out_dict = {'complete_street_score': csi}
return out_dict
ag_acres = luac.get_lutype_acreage(project_fc, p.parcel_poly_fc,
p.lutype_ag)
pct_adt_truck = {
"pct_truck_aadt": -1
} if project_type != p.ptype_fwy else truck_fwy.get_tmc_truck_data(
project_fc, project_type)
intersxn_data = intsxn.intersection_density(project_fc,
p.intersections_base_fc,
project_type)
npmrds_data = npmrds.get_npmrds_data(project_fc, project_type)
transit_data = trnsvc.transit_svc_density(project_fc, p.trn_svc_fc,
project_type)
bikeway_data = bnmi.get_bikeway_mileage_share(project_fc, p.ptype_sgr)
infill_status = urbn.projarea_infill_status(project_fc, p.comm_types_fc)
# total job + du density (base year only, for state-of-good-repair proj eval only)
job_du_dens = lu_pt_buff.point_sum_density(p.parcel_pt_fc, project_fc,
project_type,
[p.col_emptot, p.col_du],
p.ilut_sum_buffdist)
comb_du_dens = sum(list(job_du_dens.values()))
job_du_dens['job_du_perNetAcre'] = comb_du_dens
# get EJ data
ej_data = lu_pt_buff.point_sum(p.parcel_pt_fc,
