opencv透視變換

關於透視投影的幾何知識，以及求解方法，可以參考

http://media.cs.tsinghua.edu.cn/~ahz/digitalimageprocess/chapter06/chapt06_ahz.htm

http://blog.csdn.net/xiaowei_cqu/article/details/26471527

此外，opencv1那本書183頁也有講

這里，實現一個貼圖，主要參考

http://www.cnblogs.com/tiandsp/p/4033071.html

那篇博客用的matlab，這里使用opencv2實現

原圖

代碼

#include <cv.h>
#include <highgui.h>
#include <iostream>
#include <vector>
#include <fstream>

using namespace std;
using namespace cv;

Mat rawImg, dstImg,boardImg,dstboard,maskboard;
vector<Point2f>srcQuad(4);
vector<Point2f>dstQuad(4);

void on_mouse(int event, int x, int y, int flags, void *ustc)//event鼠標事件代號，x,y鼠標坐標，flags拖拽和鍵盤操作的代號
{
    Point pt;//坐標點;
    char coordinateName[16];
    static int n=0;

    if (event == CV_EVENT_LBUTTONDOWN)//左鍵按下，讀取坐標，並在圖像上該點處划圓
    {
        pt = Point(x, y);
        cout << x << " " << y << endl;
        dstQuad[n] = pt;
        n++;
        circle(boardImg, pt, 2, Scalar(255, 0, 0, 0), CV_FILLED, CV_AA, 0);//划圓
        sprintf(coordinateName, "(%d,%d)", x, y);
        putText(boardImg, coordinateName, pt, FONT_HERSHEY_SIMPLEX, 0.5, Scalar(0, 0, 0, 255), 1, 8);//在窗口上顯示坐標
        if (n >= 4)
        {
            //imshow("board", boardImg);
            cvDestroyAllWindows();
        }
    }
}
int main()
{
    rawImg = imread("lena.jpg",0);//圖片路徑
    boardImg = imread("board.jpg",0);
    dstboard = boardImg.clone();
    cvNamedWindow("board",0);
    setMouseCallback("board",on_mouse,0);
    imshow("board",boardImg);
    waitKey(0);
    imshow("raw", rawImg);
    imshow("board", boardImg);
    int imgHeigth = rawImg.rows;
    int imgWidth = rawImg.cols;
    int dstHeigth = boardImg.rows;
    int dstWidth = boardImg.cols;
    srcQuad[0] = Point2f(0,0);//左上 右上 左下 右下
    srcQuad[1] = Point2f(imgWidth-1, 0);
    srcQuad[2] = Point2f(0, imgHeigth - 1);
    srcQuad[3] = Point2f(imgWidth - 1, imgHeigth - 1);
    /*dstQuad[0] = Point2f(imgWidth*0.05, imgHeigth*0.33);
    dstQuad[1] = Point2f(imgWidth*0.9, imgHeigth*0.25);
    dstQuad[2] = Point2f(imgWidth*0.2, imgHeigth*0.7);
    dstQuad[3] = Point2f(imgWidth*0.8, imgHeigth*0.9);*/
    for (int i = 0; i < 4; i++)
        cout << dstQuad[i] << endl;

    Mat warpMatrix = getPerspectiveTransform(srcQuad, dstQuad);
    cout << warpMatrix << endl;
    warpPerspective(rawImg, dstImg, warpMatrix, Size(dstWidth, dstHeigth));

    threshold(dstImg, maskboard, 0, 1, THRESH_BINARY);
    cout << "begin a test of cout to file." << endl;
    // 保存cout流緩沖區指針
    streambuf* coutBuf = cout.rdbuf();
    ofstream of("out.txt");
    // 獲取文件out.txt流緩沖區指針
    streambuf* fileBuf = of.rdbuf();
    // 設置cout流緩沖區指針為out.txt的流緩沖區指針
    cout.rdbuf(fileBuf);
    cout << maskboard << endl;

    of.flush();
    of.close();
    // 恢復cout原來的流緩沖區指針
    cout.rdbuf(coutBuf);
    //imshow("rawboard", dstboard);
    cout << "Write Personal Information over..." << endl;
    for (int i = 0; i < dstHeigth - 1; i++)
        for (int j = 0; j < dstWidth - 1; j++)
            dstImg.at<uchar>(i, j) = dstImg.at<uchar>(i, j)*maskboard.at<uchar>(i, j)
            + dstboard.at<uchar>(i, j)*(1-maskboard.at<uchar>(i, j));
    imshow("result", dstImg);
    waitKey();
    return 0;
}

結果