for ii in range(3):
seeds, seedweight = splitclusters(datapointwts, seeds)
print('cluster split: ', time.time() - start)
gedges = coocurematrix(datapointwts, seeds)
print('coocurence matrix computed: ', time.time() - start)
gedges = prunegraph(gedges, seeds)
print('graph pruning. number of edges = ', len(gedges),
time.time() - start)
getgeojson(gedges, seeds)
printedges(gedges, seeds, datapointwts)
printclusters(seeds)
# 2. Run the map fusion code.
start = time.time()
# Work on the static map: OSM, etc.
G = OSMReader.build_road_network_from_osm_shapefile(basemap_path)
print('Generated OSM graph from shapefiles in:', time.time() - start)
seedsosm, neighborsosm = getosmmap(base_map=G)
seedsosm, neighborsosm = densifymap(seedsosm, neighborsosm, 20)
print ('osm map read/densified: ', time.time()-start, len(seedsosm), \
sum([len(neighborsosm[xx]) for xx in neighborsosm]))
clusters, clusterneighbour = getclusters(seedsosm, neighborsosm)
deadclusters = detectdeadclusters(densified_data, clusters,
clusterneighbour)
print('dead clusters found: ', time.time() - start, len(deadclusters))
clusters1, clusterneighbour1 = deletenodesfromgraph(
clusters, clusterneighbour, deadclusters)
# Work on Kharita
seedskha, neighborskha = getourmap()
seedskha, neighborskha = stitchtwomaps(seedskha, neighborskha, clusters1,
clusterneighbour1)
