如何在python中创建LogEuclidean算法？

Question

我有一个矩阵（它是DTMRI的一个）作为输入，我想找出像素点的测地距离；因为每个像素有6个点。。。然后将其输出到一个数组，其中对角线将读取所有零（因为从一个点到它自身的距离是零）。我该怎么做？在

我要知道矩阵的计算结果。在

Dij=dTLE（T1，T2）=

sqrt（trace（（logm（T1））−（logm（T2））^2），其中T1和T2是两个DT，logm是矩阵对数，通过分解T=U∧U^T定义为logm（T）=Udiag（log（diag（∧）））U^T

任何关于如何开始编码的指示将不胜感激-