How to calculate euclidean distance between pair of rows of a numpy array

Rashmi Singh picture Rashmi Singh · Apr 12, 2017 · Viewed 10.9k times · Source

I have a numpy array like:

import numpy as np
a = np.array([[1,0,1,0],
             [1,1,0,0],
             [1,0,1,0],
             [0,0,1,1]])

I would like to calculate euclidian distance between each pair of rows.

from scipy.spatial import distance
for i in range(0,a.shape[0]):
    d = [np.sqrt(np.sum((a[i]-a[j])**2)) for j in range(i+1,a.shape[0])]
    print(d)

[1.4142135623730951, 0.0, 1.4142135623730951]

[1.4142135623730951, 2.0]

[1.4142135623730951]

[]

Is there any better pythonic way to do this since i have to run this code on a huge numpy array?

Answer

comendeiro picture comendeiro · Apr 12, 2017

In terms of something more "elegant" you could always use scikitlearn pairwise euclidean distance:

from sklearn.metrics.pairwise import euclidean_distances
euclidean_distances(a,a)

having the same output as a single array.

array([[ 0.        ,  1.41421356,  0.        ,  1.41421356],
       [ 1.41421356,  0.        ,  1.41421356,  2.        ],
       [ 0.        ,  1.41421356,  0.        ,  1.41421356],
       [ 1.41421356,  2.        ,  1.41421356,  0.        ]])