outpath = "{}/_output".format(mainpath)
junkpath = "{}/junk".format(outpath)
nw70 = "{}/lines_for_network_1870.shp".format(outpath)
nw90 = "{}/lines_for_network_1890.shp".format(outpath)
centroids = "{}/county_centroids.shp".format(outpath)
out70 = "{}/_odmat_1870.csv".format(outpath)
gp = GeoProcess()
###########################################################################
###########################################################################
###########################################################################
print('computing origin-destination matrix, 1870')
odmat70 = gp.odmat_nn(nw70,
centroids,
'NHGISNAM',
tolerance=10,
criterion='speed',
speed_field='speed',
default_speed=6.966857)
print('output to csv, 1870')
gp.output_csv(odmat70, out70)
print('DONE!')
#########################################################################
# 2) centroids and distance to coast
#########################################################################
#########################################################################
#########################################################
# Fix geometries
#########################################################
print('fixing geometries, countries')
fixgeo_countries = gp.fix_geometry(adminin)
#########################################################
# Centroids
#########################################################
print('finding country centroids')
country_centroids = gp.centroids(fixgeo_countries)
#########################################################
# Closest coordinates from points to lines
#########################################################
result = gp.find_closest_coordinates_from_points_to_lines(country_centroids, coastin, intermediate_output_folder=junkpath,
input_id='ne_10m_adm', input_name='cent', line_name='coast',
verbose=True, output_object='memory:')
#########################################################
# Centroids
#########################################################
print('writing to csv')
gp.output_csv(result, csvout)
print('DONE!')