#include <vector>
#include <iostream>
#include <algorithm>
#include <opencv/highgui.h>
#include <opencv/cv.h>
#include <boost/numeric/ublas/matrix.hpp>
#include <opencv2/nonfree/nonfree.hpp>
using namespace boost::numeric::ublas;

// Homography matrix
cv::Mat H;

// Point of images save in cv::Mat
cv::Mat srcPoints;
cv::Mat disPoints;

std::vector<IplImage*> images;

int globalMaxX, globalMaxY, globalMinX, globalMinY;
// Multiplication width convert to boost::matrix
vector<double> matrixMultiplication(const cv::Mat* H, int x, int y)
{
       matrix<double> m(3,3);  // <----------------------------------------------------------------------
        vector<double> pt(3);

        pt(0) = x; pt(1) = y; pt(2) = 1;

        // if there isn't matrix of homography return original point
        if (H == 0) return pt;

        // Copy to new matrix m
        for (unsigned i = 0; i < m.size1(); ++ i)
                for (unsigned j = 0; j < m.size2(); ++ j)
                        m(i, j) = (*H).at<double>(i, j);
       // Multiplication from boost            
       pt = prod(m, pt);
       // Normalization homogeneous coordinates
       pt = pt/pt(2);
       return pt;
}