Opencv下双线性插值法进行图像放缩

关于图像放缩的算法有很多，本文主要介绍双线性插值法进行图像放缩，本文参考了： http://www.cnblogs.com/funny-world/p/3162003.html

我们设源图像src的大小为m*N，目标图像的大小为a*b。无论是放大还是缩小，我们所遵循的原则就是对于目标图像的像素点 f(i,j)，我们在源图像当中找到其对应点 f(x,y)。

其算法为：i/x=m/a  j/y=n/b 因此我们可以得到：x=i*(m/a)     y=j*(n/b)

但是为了使放缩的效果更好，可以使两个图片的几何中心重合，所以我们可以：x=(i+0.5)*(m/a)-0.5   y=(j+0.5)*(n/b)-0.5

通过这样算法，我们便可以在源图像当中找到目标图像的对应点了，但是这样计算 x,y 很可能是float型，因此我们需要继续对 x,y进行处理

我们可以找到 于P（所求的点）最为临近的四个点 Q₁₁ ,Q₁₂ ,Q₂₁,Q₂₂ ，通过下面的算法进行处理：

f(R₁)=(X₂-X)/(X₂-X1)*f(Q₁₁)+(X-X₁)/(X₂-X₁)*f(Q₂₁)

f(R₂)=(X₂-X)/(X₂-X₁)*f(Q₁₂)+(X-X₁)/(X₂-X₁)*f(Q₂₂)

f(i,j)=(Y₂-Y)/(Y₂-Y1)f(R₁)+(Y-Y₁)/(Y₂-Y₁)*f(R₂)

通过这样的处理，我们便可以计算出目标图像法 f(i,j)的像素值了，从而可以求出目标图像.

下面是Opencv的代码：

#include<opencv2\opencv.hpp>
#include<opencv2\core\core.hpp>
#include<opencv2\highgui\highgui.hpp>

using namespace std;
using namespace cv;

int main(){
    Mat src;
    src=imread("G:/save1.jpg",);    //输入图片;
    if(src.empty()){
        return -;
    }
    Mat dist=Mat((*src.rows),(*src.cols),CV_32FC3);        //目标图片，将原图放大两倍;

    //遍历目标图片的点，将目标图片的归溯到原图的点当中，若结果为浮点型，则按照算法计算；
    //归溯算法，原图大小为a*b，目标图大小m*n，i/x=a/m，j/y=b/n;

    float x,y;
    int Q11[],Q21[],Q12[],Q22[];          //Q11 Q12 Q21 Q22分别对应着与P点最为临近的四个点，P点为目标图点还原到原图的点的对应位置;
    float R1[],R2[];
    for(int i=;i<dist.rows-;i++){           //这里减2是因为对于(src.rows-1,src.cols-1)的位置对应的是（dist.rows-2,dits.cols-2),再大一点就会超过内存。
        for(int j=;j<dist.cols-;j++){
            x=(i+0.5)*0.5-0.5;               //通过加减0.5，将像素点移到中心点的位置，使放缩效果更好。
            y=(j+0.5)*0.5-0.5;

            Q11[]=(int)x;
            Q12[]=(int)x;
            Q21[]=(int)x+;
            Q22[]=(int)x+;

            Q11[]=(int) y;
            Q21[]=(int) y;
            Q12[]=(int) y+;
            Q22[]=(int) y+;

            // 插值法算法如下：
            // X1=(int)x; X2=(int)x+1; Y1=(int)y; Y2=(int)y+1;
            // f(R1)=(X2-X)/(X2-X1)*f(Q11)+(X-X1)/(X2-X1)*f(Q21);
            // f(R2)=(X2-X)/(X2-X1)*f(Q12)+(X-X1)/(X2-X1)*f(Q22);
            // f(i,j)=(Y2-Y)/(Y2-Y1)f(R1)+(Y-Y1)/(Y2-Y1)*f(R2);
            for(int m=;m<;m++){
                R1[m]=((Q21[]-x)*src.at<Vec3b>(Q11[],Q11[])[m])+((x-Q11[])*src.at<Vec3b>(Q21[],Q21[])[m]);
                R2[m]=((Q21[]-x)*src.at<Vec3b>(Q12[],Q12[])[m])+((x-Q11[])*src.at<Vec3b>(Q22[],Q22[])[m]);
                dist.at<Vec3f>(i,j)[m]=(Q12[]-y)*R1[m]+(y-Q11[])*R2[m];
            }
        }
    }
    imwrite("G:/save.jpg",dist);          //输出图片;
}