python输入两个点的坐标,在python中获取两个（X，Y）坐标之间所有点的最快方法...

So I have a shapely LineString:

print np.round(shapely_intersecting_lines.coords).astype(np.int)

>>> array([[ 1520, -1140],

[ 1412, -973]])

This can be interpreted as a numpy array as well as seen above.

I want to get all the points in between, that is I want to get the points of the line in between as integer values. The output should be something like this:

array([[ 1520, -1140],

[ 1519, -1139],

[ 1519, -1138],

...,

[ 1413, -975],

[ 1412, -974],

[ 1412, -973]], dtype=int32)

I posted this earlier in gis.stackexchange hoping there was a solution in shapely that was efficient. The solution was good at first, however, the solution is now too slow as I run this over 50000 times in my code. On my computer each loop takes about 0.03s resulting in over a day of running. It is too slow