图像放大并进行BiCubic插值

BiCubic插值原理：

构造BiCubic函数：

其中，a取-0.5.

[source: http://en.wikipedia.org/wiki/Bicubic_interpolation]

BiCubic函数具有如下形状：

[source:  R. Keys, (1981). "Cubic convolution interpolation for digital image processing". IEEE Transactions on Signal Processing, Acoustics, Speech, and Signal Processing 29 (6): 1153–1160.]

对待插值的像素点(x,y)（x和y可以为浮点数），取其附近的4x4邻域点(xi,yj), i,j = 0,1,2,3。按如下公式进行插值计算：

实现代码：

#include "opencv2/imgproc/imgproc.hpp"
#include "opencv2/highgui/highgui.hpp"
#include <iostream>
#include <cmath>
#include <fstream>

using namespace cv;
using namespace std;

#define PI 3.14159265

float BiCubicPoly(float x);
void MyScaleBiCubicInter(Mat& src, Mat& dst, float TransMat[3][3]);
/**
 * @function main
 */
int main( int argc, char** argv )
{
  // load image
  char* imageName = "images/Lenna_256.png";

  Mat image;
  image = imread(imageName,1);
  if(!image.data)
  {
	  cout << "No