Pca compact trick

PCA PCA_trick(const Mat& pcaset, int maxComponents)
{
     int n = pcaset.rows, p = pcaset.cols;
     cout << "\tcalculating 'true' means, varance and standard deviation..." << endl;
     Mat means(1, p, CV_32FC1);
     Mat variance(1, p, CV_32FC1);
     for (size_t i = 0; i < p; i++)
     {
       float avg = mean(pcaset.col(i)).val[0];
       means.at<float>(0, i) = avg;
       Mat p2 = Mat(1, n, CV_32F);
       for (size_t j = 0; j < n; j++)
         p2.at<float>(0, j) = pow(pcaset.at<float>(j, i), 2);
       variance.at<float>(0, i) = (1 / (float)n) * sum(p2).val[0] - pow(avg, 2);
     }

     //covariance matrix, AA', not the A'A like usual
     Mat M;
     Mat centred(n, p, CV_32FC1);
     for (size_t i = 0; i < n; i++)
       centred.row(i) = (pcaset.row(i) - means) / variance;
     mulTransposed(centred, M, 0);

     //compute eigenvalues and eigenvectors
     PCA pca;
     pca = PCA(M, cv::Mat(), CV_PCA_DATA_AS_ROW, maxComponents);

     //this is the compact trick
     pca.mean = means;
     pca.eigenvectors = pca.eigenvectors * centred;
     pca.eigenvectors = pca.eigenvectors.rowRange(Range(0, maxComponents));

     return pca;
}

Source: stackoverflow.com

Similar codes

0 vote198 views1 answer
What does map.count() return in c++
C++
0 vote207 views1 answer
Waiting in a serial as the spool reflect the queue operation. Demonstrate Printer Behavior in context of Queue.Subject to the Scenario implement the Pop and Push Using C++.
C++
4 votes189 views1 answer
Git branch in my bash prompt
C++
0 vote165 views4 answers
C++ to c converter
C++
0 vote210 views1 answer
How to seek to the start of afile in c++
C++
0 vote215 views1 answer
How to use qdebug c++
C++
0 vote204 views1 answer
Print number with leading zeros
C++
0 vote200 views1 answer
How to use mersenne_twister_engine in c++ to generate random numbers
C++
1 vote211 views1 answer
Ue4 ftext to int
C++
4 votes191 views2 answers
C++ function as pointer
C++