示例#1
0
from numpy import shape
import matplotlib.pyplot as plt
import numpy as np
import os.path
import pandas as pd
import win32com.client as com
import timeit
import copy

# Load visum
path = "E:\\Thenuwan\\Singapore_Calibration"
version_file = "7_HeadwayBased_remove_unwanted_lines_speed_up.ver"
visum_path = os.path.join(path, version_file)
Visum = com.Dispatch("Visum.Visum.170")

ocv.loadVisum(VisumComDispatch=Visum, verPath=visum_path)

observed_stop_point_df = pd.read_csv(
    "E:\\Thenuwan\\Singapore_Calibration\\data\\observed\\observed_stop_data_all_movements_2902200.csv"
)
observed_line_route_df = pd.read_csv(
    "E:\\Thenuwan\\Singapore_Calibration\\data\\observed\\observed_line_route_data_14032020.csv"
)
observed_line_route_df["LineName"] = observed_line_route_df["LineName"].astype(
    str)
observed_line_route_df["Name"] = observed_line_route_df["Name"].astype(str)

# Create dataframe and assign values
df_rmsn_columns = [
    'in_veh', 'transfer_walk', 'origin_wait', 'transfer_wait',
    'transfer_penalty', 'pax_trans_total_rmsn', 'pax_trans_walkb_rmsn',
import win32com.client as com
import timeit

# Load Visum Version and create a Network Object
path = "E:\\Thenuwan\\Singapore_Calibration"
#ver_file = "11_Remove_premium_0.0475_mean_headway_28042020.ver"
ver_file = "11.2_Remove_premium_0.01_mean_head_way_27042020.ver"
version_path = os.path.join(path, ver_file)
Visum = com.Dispatch("Visum.Visum.170")

# save results
save_result_path = "E:\\Thenuwan\\Singapore_Calibration\\data\\results\\run_15\\hp_set_13_spsa_10072020_run15_calibrated_model.csv"
save_results_all_path = "E:\\Thenuwan\\Singapore_Calibration\\data\\results\\run_15\\hp_set_13_spsa_all_19072020_run15_calibrated_model.csv"

# load visum file
ocv.loadVisum(VisumComDispatch=Visum, verPath=version_path)

# read observed data
observed_stop_df = pd.read_csv(
    "E:\\Thenuwan\\Singapore_Calibration\\data\\observed\\observed_stop_data_all_movements_total_transfers_g_5000_07102020.csv"
)
observed_line_route_df = pd.read_csv(
    "E:\\Thenuwan\\Singapore_Calibration\\data\\observed\\observed_line_route_data_24042020.csv"
)
observed_line_route_df["LineName"] = observed_line_route_df["LineName"].astype(
    str)
observed_line_route_df["Name"] = observed_line_route_df["Name"].astype(str)

# saving results from the calibration in a dictionary

results_columns = [