图像拼接---图片柱面投影简单实现

算法思想参考：http://blog.csdn.net/weixinhum/article/details/50611750

柱面投影是图片拼接的前期的一部分工作，以下代码只是简单的实现了投影，还可以优化，

柱面半径设置位图片宽度的一半，即 R = width/2

代码运算流程是对于dst图片上的每一个像素点，通过公式计算出src上对应的位置(hnum,wnum)，把src上这个位置的像素值赋值给dst。

#include <iostream>
#include <opencv2/opencv.hpp>using namespace std;
using namespace cv;/*
实现一个简单的图像投影实例，
算法根据csdn博客http://blog.csdn.net/weixinhum/article/details/50611750
该代码对目录下的test图片进行柱面投影
*/
int main()
{Mat src = imread("test.jpg");Mat dst(src.rows,src.cols,src.type());int width = src.cols, height = src.rows;double x, y;double R = width / 2;int drcpoint;for (int hnum = 0; hnum < height;hnum++){for (int wnum = 0; wnum < width;wnum++){double k = R / sqrt(R*R + (wnum - width / 2)*(wnum - width / 2));x = (wnum - width / 2) / k + width / 2;y = (hnum - height / 2) / k + height / 2;if (0 < x && x < width && 0 < y &&y < height){dst.at<Vec3b>(hnum, wnum) = src.at<Vec3b>(int(y), int(x));}}}imshow("origin Image",src);imshow("柱面投影图",dst);waitKey(0);return 0;
}

运行结果：