双三次插值C代码（利用opencv）

双三次插值C代码（利用opencv）

phasecubic2.cpp

D:\文件及下载相关\文档\Visual Studio 2010\Projects\phasecubic2\phasecubic2

// phasecubic2.cpp : 定义控制台应用程序的入口点。
//
#include "stdafx.h"
#include <opencv2/opencv.hpp>
#include <opencv2/highgui/highgui.hpp>
#include <opencv2/imgproc/imgproc.hpp>
#include <string>
using namespace cv;
using namespace std;

int _tmain(int argc, _TCHAR* argv[])
{
    cv::Mat matSrc, matDst1, matDst2;
    matSrc = cv::imread("003.png",  | );
    matDst1 = cv::Mat(cv::Size(, ), matSrc.type(), cv::Scalar::all());
    matDst2 = cv::Mat(matDst1.size(), matSrc.type(), cv::Scalar::all());

    double scale_x = (double)matSrc.cols / matDst1.cols;
    double scale_y = (double)matSrc.rows / matDst1.rows;
    int iscale_x = cv::saturate_cast<int>(scale_x);
    int iscale_y = cv::saturate_cast<int>(scale_y);

    for (int j = ; j < matDst1.rows; ++j)
    {
        float fy = (float)((j + 0.5) * scale_y - 0.5);
        int sy = cvFloor(fy);
        fy -= sy;
        sy = std::min(sy, matSrc.rows - );
        sy = std::max(, sy);

        const float A = -0.75f;

        float coeffsY[];
        coeffsY[] = ((A*(fy + ) - *A)*(fy + ) + *A)*(fy + ) - *A;
        coeffsY[] = ((A + )*fy - (A + ))*fy*fy + ;
        coeffsY[] = ((A + )*( - fy) - (A + ))*( - fy)*( - fy) + ;
        coeffsY[] = .f - coeffsY[] - coeffsY[] - coeffsY[];

        short cbufY[];
        cbufY[] = cv::saturate_cast<short>(coeffsY[] * );
        cbufY[] = cv::saturate_cast<short>(coeffsY[] * );
        cbufY[] = cv::saturate_cast<short>(coeffsY[] * );
        cbufY[] = cv::saturate_cast<short>(coeffsY[] * );

        for (int i = ; i < matDst1.cols; ++i)
        {
            float fx = (float)((i + 0.5) * scale_x - 0.5);
            int sx = cvFloor(fx);
            fx -= sx;

            if (sx < ) {
                fx = , sx = ;
            }
            if (sx >= matSrc.cols - ) {
                fx = , sx = matSrc.cols - ;
            }

            float coeffsX[];
            coeffsX[] = ((A*(fx + ) - *A)*(fx + ) + *A)*(fx + ) - *A;
            coeffsX[] = ((A + )*fx - (A + ))*fx*fx + ;
            coeffsX[] = ((A + )*( - fx) - (A + ))*( - fx)*( - fx) + ;
            coeffsX[] = .f - coeffsX[] - coeffsX[] - coeffsX[];

            short cbufX[];
            cbufX[] = cv::saturate_cast<short>(coeffsX[] * );
            cbufX[] = cv::saturate_cast<short>(coeffsX[] * );
            cbufX[] = cv::saturate_cast<short>(coeffsX[] * );
            cbufX[] = cv::saturate_cast<short>(coeffsX[] * );

            for (int k = ; k < matSrc.channels(); ++k)
            {
                matDst1.at<cv::Vec3b>(j, i)[k] =
                    abs(
                    (matSrc.at<cv::Vec3b>(sy-, sx-)[k] * cbufX[] * cbufY[]
                    + matSrc.at<cv::Vec3b>(sy, sx-)[k] * cbufX[] * cbufY[]
                    + matSrc.at<cv::Vec3b>(sy+, sx-)[k] * cbufX[] * cbufY[]
                    + matSrc.at<cv::Vec3b>(sy+, sx-)[k] * cbufX[] * cbufY[]
                    + matSrc.at<cv::Vec3b>(sy-, sx)[k] * cbufX[] * cbufY[]
                    + matSrc.at<cv::Vec3b>(sy, sx)[k] * cbufX[] * cbufY[]
                    + matSrc.at<cv::Vec3b>(sy+, sx)[k] * cbufX[] * cbufY[]
                    + matSrc.at<cv::Vec3b>(sy+, sx)[k] * cbufX[] * cbufY[]
                    + matSrc.at<cv::Vec3b>(sy-, sx+)[k] * cbufX[] * cbufY[]
                    + matSrc.at<cv::Vec3b>(sy, sx+)[k] * cbufX[] * cbufY[]
                    + matSrc.at<cv::Vec3b>(sy+, sx+)[k] * cbufX[] * cbufY[]
                    + matSrc.at<cv::Vec3b>(sy+, sx+)[k] * cbufX[] * cbufY[]
                    + matSrc.at<cv::Vec3b>(sy-, sx+)[k] * cbufX[] * cbufY[]
                    + matSrc.at<cv::Vec3b>(sy, sx+)[k] * cbufX[] * cbufY[]
                    + matSrc.at<cv::Vec3b>(sy+, sx+)[k] * cbufX[] * cbufY[]
                    + matSrc.at<cv::Vec3b>(sy+, sx+)[k] * cbufX[] * cbufY[]
                    ) >> );
            }
        }
    }
    //cv::imwrite("cubic_1.jpg", matDst1);
    //cv::imwrite("cubic_2.jpg", matDst2);
    const char SrcWinNa[] = "SrcWindow";
    const char ResWinNa[] = "ResWindow";
    cvNamedWindow(SrcWinNa, CV_WINDOW_AUTOSIZE);
    imshow(SrcWinNa, matDst1);
    cv::resize(matSrc, matDst2, matDst1.size(), , , );
    cvNamedWindow(ResWinNa, CV_WINDOW_AUTOSIZE);
    imshow(ResWinNa, matDst2);
    cvWaitKey();
    return ;
    //getchar();
    //return 0;
}