def test_point_projects_to_edge(self):
# p = (114.83299055, 26.8892277)
p = (121.428387, 31.027371)
a = time.time()
edges, segments = self.sg.point_projects_to_edges(p, 0.01)
print(time.time() - a)
if self.show_plots:
plt.figure()
s2g.plot_lines(MultiLineString(segments), color='orange') # original roads
for i in range(0, len(edges)):
s, e = edges[i]
sxy = self.sg.node_xy[s]
exy = self.sg.node_xy[e]
plt.plot([sxy[0], exy[0]], [sxy[1], exy[1]], color='green') # graph edges
plt.plot(p[0], p[1], color='red', markersize=12, marker='o') # bridges
plt.show()
评论列表
文章目录