def main():
parser = ArgumentParser('mrt temperature visualizer')
parser.add_argument('--dir', type=str, dest='dir', default='.',
help='path to directory containing Icarus run data')
parser.add_argument('--log', type=str, dest='log', default=None,
help='path to file to save the process log; not saved by default')
parser.add_argument('--level', type=str, dest='level', default='info',
choices=('notset', 'debug', 'info', 'warning', 'error', 'critical'),
help='verbosity of the process log')
args = parser.parse_args()
handlers = []
handlers.append(log.StreamHandler())
if args.log is not None:
handlers.append(log.FileHandler(args.log, 'w'))
log.basicConfig(
format='%(asctime)s %(levelname)s %(filename)s:%(lineno)s %(message)s',
level=getattr(log, args.level.upper()),
handlers=handlers
)
database = SqliteUtil('database.db', readonly=True)
kind = 'mrt'
map_mrt_temperature(database, kind)
def main():
parser = ArgumentParser('mag abm parser')
parser.add_argument('--folder', type=str, dest='folder', default='.')
parser.add_argument('--log', type=str, dest='log', default=None)
parser.add_argument('--level',
type=str,
dest='level',
default='info',
choices=('notset', 'debug', 'info', 'warning', 'error',
'critical'))
args = parser.parse_args()
handlers = []
handlers.append(log.StreamHandler())
if args.log is not None:
handlers.append(log.FileHandler(args.log, 'w'))
log.basicConfig(
format='%(asctime)s %(levelname)s %(filename)s:%(lineno)s %(message)s',
level=getattr(log, args.level.upper()),
handlers=handlers)
path = lambda x: os.path.abspath(os.path.join(args.folder, x))
home = path('')
log.info('Running MAG ABM parsing tool.')
log.info(f'Loading run data from {home}.')
database = SqliteUtil(path('database.db'))
config = ConfigUtil.load_config(path('config.json'))
trips_file = config['population']['trips_file']
persons_file = config['population']['persons_file']
households_file = config['population']['households_file']
if not ready(trips_file, households_file, persons_file):
log.warning('Dependent data not parsed or generated.')
exit(1)
elif complete(database):
log.warning(
'Population data already parsed. Would you like to replace it? [Y/n]'
)
if input().lower() not in ('y', 'yes', 'yeet'):
log.info(
'User chose to keep existing population data; exiting parsing tool.'
)
exit()
try:
log.info('Starting population parsing.')
parse_abm(database, trips_file, households_file, persons_file)
except:
log.exception('Critical error while parsing population; '
'terminating process and exiting.')
exit(1)
def main():
parser = ArgumentParser('mrt temperature parser', add_help=False)
parser.add_argument('--help', action='help', default=SUPPRESS,
help='show this help menu and exit process')
parser.add_argument('--dir', type=str, dest='dir', default='.',
help='path to directory containing Icarus run data')
parser.add_argument('--log', type=str, dest='log', default=None,
help='path to file to save the process log; not saved by default')
parser.add_argument('--level', type=str, dest='level', default='info',
choices=('notset', 'debug', 'info', 'warning', 'error', 'critical'),
help='verbosity of the process log')
args = parser.parse_args()
handlers = []
handlers.append(log.StreamHandler())
if args.log is not None:
handlers.append(log.FileHandler(args.log, 'w'))
log.basicConfig(
format='%(asctime)s %(levelname)s %(filename)s:%(lineno)s %(message)s',
level=getattr(log, args.level.upper()),
handlers=handlers
)
path = lambda x: os.path.abspath(os.path.join(args.dir, x))
home = path('')
log.info('Running mrt temperature parsing tool.')
log.info(f'Loading run data from {home}.')
config = ConfigUtil.load_config(path('config.json'))
database = SqliteUtil(path('database.db'))
path = config['network']['exposure']['mrt_dir']
try:
log.info('Starting mrt temperature parsing.')
parse_mrt(
database,
path,
src_epsg=4326,
prj_epsg=2223,
bounds=50,
steps=96
)
except:
log.exception('Critical error while running mrt temperature '
'parsing; cleaning up and terminating.')
handlers.append(log.StreamHandler())
if args.log is not None:
handlers.append(log.FileHandler(args.log, 'w'))
log.basicConfig(
format='%(asctime)s %(levelname)s %(filename)s:%(lineno)s %(message)s',
level=getattr(log, args.level.upper()),
handlers=handlers)
path = lambda x: os.path.abspath(os.path.join(args.folder, x))
home = path('')
config = ConfigUtil.load_config(path('config.json'))
log.info('Running daymet exposure analysis tool.')
log.info(f'Loading run data from {home}.')
database = SqliteUtil(path('database.db'))
exposure = Exposure(database)
if not exposure.ready():
log.error('Dependent data not parsed or generated; see warnings.')
exit(1)
elif exposure.complete():
log.warn(
'Exposure analysis already run. Would you like to run it again? [Y/n]')
if input().lower() not in ('y', 'yes', 'yeet'):
log.info(
'User chose to keep existing exposure analysis; exiting analysis tool.'
)
exit()
try:
ax = plt.subplot(111)
ax.set_title('Maricopa County', fontsize=14)
graph = nx.Graph()
print('Adding nodes to the graph.')
for node_id, point in nodes:
x, y = xy(point)
graph.add_node(node_id, pos=(x, y))
print('Adding links to the graph.')
for _, source_node, terminal_node in links:
graph.add_edge(source_node, terminal_node)
print('Drawing the graph.')
pos = nx.get_node_attributes(graph, 'pos')
nx.draw_networkx(graph, pos=pos, ax=ax, with_labels=False, node_size=0)
print('Saving the graph.')
plt.tight_layout()
plt.savefig('result/network_usage.png', dpi=600)
if __name__ == '__main__':
database = SqliteUtil('database.db')
map_network_usage(database)
def main():
parser = ArgumentParser()
main = parser.add_argument_group('main')
main.add_argument('file',
type=str,
help='file path to save the exported routes to')
main.add_argument('--epsg',
dest='epsg',
type=int,
default=2223,
help='epsg system to convert routes to; default is 2223')
common = parser.add_argument_group('common')
common.add_argument(
'--folder',
type=str,
dest='folder',
default='.',
help='file path to the directory containing Icarus run data'
'; default is the working directory')
common.add_argument(
'--log',
type=str,
dest='log',
default=None,
help=
'file path to save the process log; by default the log is not saved')
common.add_argument(
'--level',
type=str,
dest='level',
default='info',
help='verbosity level of the process log; default is "info"',
choices=('notset', 'debug', 'info', 'warning', 'error', 'critical'))
common.add_argument(
'--replace',
dest='replace',
action='store_true',
default=False,
help='automatically replace existing data; do not prompt the user')
args = parser.parse_args()
handlers = []
handlers.append(log.StreamHandler())
if args.log is not None:
handlers.append(log.FileHandler(args.log, 'w'))
log.basicConfig(
format='%(asctime)s %(levelname)s %(filename)s:%(lineno)s %(message)s',
level=getattr(log, args.level.upper()),
handlers=handlers)
path = lambda x: os.path.abspath(os.path.join(args.folder, x))
home = path('')
log.info('Running link export tool.')
log.info(f'Loading run data from {home}.')
database = SqliteUtil(path('database.db'), readonly=True)
try:
export_links(database, args.file, 2223, args.epsg)
except:
log.exception('Critical error while exporting routes:')
exit(1)
database.close()
def main():
parser = ArgumentParser('daymet air temperature parser')
parser.add_argument('--dir',
type=str,
dest='dir',
default='.',
help='path to directory containing Icarus run data')
parser.add_argument(
'--log',
type=str,
dest='log',
default=None,
help='path to file to save the process log; not saved by default')
parser.add_argument('--level',
type=str,
dest='level',
default='info',
choices=('notset', 'debug', 'info', 'warning', 'error',
'critical'),
help='verbosity of the process log')
args = parser.parse_args()
handlers = []
handlers.append(log.StreamHandler())
if args.log is not None:
handlers.append(log.FileHandler(args.log, 'w'))
log.basicConfig(
format='%(asctime)s %(levelname)s %(filename)s:%(lineno)s %(message)s',
level=getattr(log, args.level.upper()),
handlers=handlers)
path = lambda x: os.path.abspath(os.path.join(args.dir, x))
home = path('')
config = ConfigUtil.load_config(path('config.json'))
database = SqliteUtil(path('database.db'))
tmin_files = config['network']['exposure']['tmin_files']
tmax_files = config['network']['exposure']['tmax_files']
day = config['network']['exposure']['day']
steps = config['network']['exposure']['steps']
log.info('Running roads parsing tool.')
log.info(f'Loading run data from {home}.')
if not ready(database, tmin_files, tmax_files):
log.error('Process dependencies not met; see warnings and '
'docuemntation for more details.')
exit(1)
if complete(database):
log.info('All or some of this process is already complete. '
' Would you like to proceed? [Y/n]')
valid = ('y', 'n', 'yes', 'no', 'yee', 'naw')
response = input().lower()
while response not in valid:
print('Try again; would you like to proceed? [Y/n]')
response = input().lower()
if response in ('n', 'no', 'naw'):
log.info('User chose to terminate process.')
exit()
try:
log.info('Starting road parsing.')
parse_temperatures(database, tmin_files, tmax_files, steps, day, 4326,
2223)
except:
log.exception('Critical error while parsing roads; '
'terminating process and exiting.')
exit(1)
def main():
parser = ArgumentParser('road network parser')
parser.add_argument('--dir',
type=str,
dest='dir',
default='.',
help='path to directory containing Icarus run data')
parser.add_argument(
'--log',
type=str,
dest='log',
default=None,
help='path to file to save the process log; not saved by default')
parser.add_argument('--level',
type=str,
dest='level',
default='info',
choices=('notset', 'debug', 'info', 'warning', 'error',
'critical'),
help='verbosity of the process log')
args = parser.parse_args()
handlers = []
handlers.append(log.StreamHandler())
if args.log is not None:
handlers.append(log.FileHandler(args.log, 'w'))
log.basicConfig(
format='%(asctime)s %(levelname)s %(filename)s:%(lineno)s %(message)s',
level=getattr(log, args.level.upper()),
handlers=handlers)
path = lambda x: os.path.abspath(os.path.join(args.dir, x))
home = path('')
log.info('Running roads parsing tool.')
log.info(f'Loading run data from {home}.')
database = SqliteUtil(path('database.db'))
networkpath = path('input/network.xml.gz')
if not ready(networkpath):
log.warning('Dependent data not parsed or generated.')
log.warning(
'Roads parsing dependencies include network generation as well '
'as exposure and regions parsing.')
exit(1)
elif complete(database):
log.warning(
'Roads already parsed. Would you like to replace it? [Y/n]')
if input().lower() not in ('y', 'yes', 'yeet'):
log.info(
'User chose to keep existing roads; exiting parsing tool.')
exit()
try:
log.info('Starting roads parsing.')
parse_roads(database, networkpath, 2223, 2223)
except:
log.exception('Critical error while parsing roads; '
'terminating process and exiting.')
exit(1)
persons = database.fetch_rows()
total = {}
data = []
for _, age, exposure in persons:
adj = (age // 5) * 5
if adj in total:
total[adj][0] += exposure
total[adj][1] += 1
else:
total[adj] = [exposure, 1]
for age, (total, count) in total.items():
data.append((age, total / count))
data = pd.DataFrame(data,
columns=('age (years)', 'average exposure (°C·sec)'))
axes = sns.barplot(x=data['age (years)'],
y=data['average exposure (°C·sec)'],
color='royalblue')
axes.set_title('Exposure By Age')
plot = axes.get_figure()
plot.savefig('result/age_exposure.png', bbox_inches='tight')
plot.clf()
if __name__ == '__main__':
database = SqliteUtil('database.db', readonly=True)
main(database)