Computing Mahalanobis Distance Between Set of Points and Set of Reference Points

Question

Royi Avital 2015-8-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/232943-computing-mahalanobis-distance-between-set-of-points-and-set-of-reference-points

编辑： Royi Avital 2015-8-15

采纳的回答： Royi Avital

在 MATLAB Online 中打开

Hello,

I have an n x p matrix - mX which is composed of n points in R^p.

I have another m x p matrix - mY which is composed of m reference points in R^p.

I would like to create an n x m matrix - mD which is the Mahalanobis Distance matrix.

D(i, j) means the Mahalanobis Distance between point j in mX, mX(j, :) and point i in mY, mY(i, :).

Namely, is computes the following:

mD(i, j) = (mX(j, :) - mY(i, :)) * inv(mC) * (mX(j, :) - mY(i, :)).';

Where mC is the given Mahalanobis Distance PSD Matrix.

It is easy to be done in a loop, is there a way to vectorize it?

Namely, is the a function which its inputs are mX, mY and mC and its output is mD and fully vectorized without using any MATLAB toolbox?

Thank You.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Royi Avital 2015-8-10

Anyone? Any Vectorization wizard? Thank You.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Royi Avital 2015-8-15

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/232943-computing-mahalanobis-distance-between-set-of-points-and-set-of-reference-points#answer_189402

编辑：Royi Avital 2015-8-15

在 MATLAB Online 中打开

Here is the solution fully vectorized (Though uses much more multiplications than needed):

mA = reshape(bsxfun(@minus, permute(mY, [1, 3, 2]), permute(mX, [3, 1, 2])), [(m * n), p]);
mD = reshape(diag(A* inv(mC) * A.'), [m, n]);

If anyone has faster way (Not necessarily fully vectorized) I'd be happy to see.

Thank You.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Image Analyst 2015-8-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/232943-computing-mahalanobis-distance-between-set-of-points-and-set-of-reference-points#answer_188739

Use *mahal(* ) in the Statistics and Machine Learning Toolbox. I haven't used it yet so I don't have any demo for you. Why do you need it?

2 个评论
显示无隐藏无

Royi Avital 2015-8-9

Hi, First I would like a solution without any toolbox as I noted above. Second, the function you suggested doesn't allow predefined weighting PSD Matrix (It infers it from the data). There's also `pdist`, yet I want to be toolboxes independent. Thank You.

Image Analyst 2015-8-9

What's m, n, and p? Your code is already at least partially vectorized. If m, n, and p are less than a few million, then it probably won't take much time at all even if it's not vectorized.

请先登录，再进行评论。