6

我有两个清单

available_points = [[2,3], [4,5], [1,2], [6,8], [5,9], [51,35]]

和

solution = [[3,5], [2,1]]

我试图弹出一个点available_points并将其附加到solution从该点到所有点的欧几里得距离之和solution最大的地方。

所以，我会得到这个

solution = [[3,5], [2,1], [51,35]]

我能够像这样选择最初的 2 个最远点，但不知道如何继续。

import numpy as np
from scipy.spatial.distance import pdist, squareform

available_points = np.array([[2,3], [4,5], [1,2], [6,8], [5,9], [51,35]])

D = squareform(pdist(available_points)
I_row, I_col = np.unravel_index(np.argmax(D), D.shape)
solution = available_points[[I_row, I_col]]

这给了我

solution = array([[1, 2], [51, 35]])

4

4 回答 4

2

你可以使用cdist-

In [1]: from scipy.spatial.distance import cdist

In [2]: max_pt=available_points[cdist(available_points, solution).sum(1).argmax()]

In [3]: np.vstack((solution, max_pt))
Out[3]: 
array([[ 3,  5],
       [ 2,  1],
       [51, 35]])
于 2018-01-12T06:29:59.367 回答
2

既然你标记numpy 

import numpy as np 

solution=np.array(solution)
available_points=np.array(available_points)
l=[]
for x in solution:
    l.append(np.linalg.norm(available_points-x, keepdims=True,axis=1))


np.append(solution,[available_points[np.argmax(np.array(l).sum(0))]],axis=0)
Out[237]: 
array([[ 3,  5],
       [ 2,  1],
       [51, 35]])
于 2018-01-12T05:19:22.197 回答
2

您可以使用 max 函数在“available_points”列表中找到最大值，然后将“available_points”列表的最大值附加到“解决方案”列表中！我还附上了输出的截图！

available_points = [[2,3], [4,5], [1,2], [6,8], [5,9], [51,35]];
solution = [[3,5], [2,1]]
solution.append(max(available_points));
print(solution);

输出

于 2018-01-12T05:08:40.410 回答
0

我想通了使用cdist

from scipy.spatial.distance import cdist

d = cdist(solution, available_points)

distances = []

for q in range(len(available_points)):
    y = d[:,q]
    distances.append(sum(y))

# Largest Distance
max_point = available_points[distances.index(max(distances))]

# Update datasets
solution = np.append(solution, [max_point], axis=0)
universe = np.delete(available_points, max_index, 0)
于 2018-01-12T06:33:37.290 回答