示例#1
0
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
0
#########################################################################
# 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!')