OpenCV - Как я могу рассчитать расстояние Римана с помощью C++?

Question

расчет расстояния Римана для двух матриц ковариации А и В в Matlab является:

C = sqrt(sum(log(eig(A,B)).^2)) 

Что является эквивалентом для OpenCV с C++?

Answer 1

Теперь я получил те же результаты, что и Matlab. Следит мое решение C++:

double calculateRiemannianDistance(cv::Mat A, cv::Mat B) { 
Eigen::MatrixXf e_A, e_B; 
cv::cv2eigen(A, e_A); 
cv::cv2eigen(B, e_B); 

/* compute generalized eigenvalues */ 
Eigen::GeneralizedEigenSolver<Eigen::MatrixXf> ges; 
ges.compute(e_A, e_B); 
cv::Mat gev; 
cv::eigen2cv(Eigen::MatrixXf(ges.eigenvalues().real()), gev); 

/* riemannian distance */ 
cv::Mat partial_res; 
cv::log(gev, partial_res); 
cv::pow(partial_res, 2, partial_res); 
double distance = sqrt(cv::sum(partial_res)[0]); 
return distance; 
}