使用pdist的Python中的字符串距离矩阵

时间:2017-09-27 16:18:22

标签: python string jaro-winkler pdist

如何计算Python中字符串的Jaro Winkler距离矩阵?

我有大量手工输入的字符串(姓名和记录号),我试图在列表中找到重复项,包括可能在拼写上有轻微变化的重复项。 response to a similar question建议使用Scipy的pdist函数和自定义距离函数。我试图用Levenshtein包中的jaro_winkler函数实现这个解决方案。这个问题是jaro_winkler函数需要字符串输入,而pdict函数似乎需要2D数组输入。

示例:

import numpy as np
from scipy.spatial.distance import pdist
from Levenshtein import jaro_winkler

fname = np.array(['Bob','Carl','Kristen','Calr', 'Doug']).reshape(-1,1)
dm = pdist(fname, jaro_winkler)
dm = squareform(dm)

预期输出 - 像这样:

          Bob  Carl   Kristen  Calr  Doug
Bob       1.0   -        -       -     -
Carl      0.0   1.0      -       -     -
Kristen   0.0   0.46    1.0      -     -
Calr      0.0   0.93    0.46    1.0    -
Doug      0.53  0.0     0.0     0.0   1.0

实际错误:

jaro_winkler expected two Strings or two Unicodes

我假设这是因为jaro_winkler函数看到的是ndarray而不是字符串,而且我不确定如何在pdist函数的上下文中将函数输入转换为字符串。 / p>

是否有人建议允许此功能?提前谢谢!

3 个答案:

答案 0 :(得分:10)

你需要包裹距离函数,就像我在下面的例子中展示的Levensthein距离

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

# my list of strings
strings = ["hello","hallo","choco"]

# prepare 2 dimensional array M x N (M entries (3) with N dimensions (1)) 
transformed_strings = np.array(strings).reshape(-1,1)

# calculate condensed distance matrix by wrapping the Levenshtein distance function
distance_matrix = pdist(transformed_strings,lambda x,y: distance(x[0],y[0]))

# get square matrix
print(squareform(distance_matrix))

Output:
array([[ 0.,  1.,  4.],
       [ 1.,  0.,  4.],
       [ 4.,  4.,  0.]])

答案 1 :(得分:0)

对于任何有类似问题的人 - 我刚发现的一个解决方案是从pdist函数中提取相关代码并在jaro_winkler函数输入中添加[0]以从numpy数组中调用字符串。

示例:

X = np.asarray(fname, order='c')
s = X.shape
m, n = s
dm = np.zeros((m * (m - 1)) // 2, dtype=np.double)

k = 0
for i in xrange(0, m - 1):
    for j in xrange(i + 1, m):
        dm[k] = jaro_winkler(X[i][0], X[j][0])
        k = k + 1

dms = squareform(dm)

即使这个算法有效,我仍然想知道是否有一个“正确”的计算机技巧方法来使用pdist函数。谢谢,希望这有助于某人!

答案 2 :(得分:0)

这是一个简洁的解决方案,不需要numpy或scipy:

from Levenshtein import jaro_winkler
data = ['Bob','Carl','Kristen','Calr', 'Doug']
dm = [[ jaro_winkler(a, b) for b in data] for a in data]
print('\n'.join([''.join([f'{item:6.2f}' for item in row]) for row in dm]))

  1.00  0.00  0.00  0.00  0.53
  0.00  1.00  0.46  0.93  0.00
  0.00  0.46  1.00  0.46  0.00
  0.00  0.93  0.46  1.00  0.00
  0.53  0.00  0.00  0.00  1.00