// 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;

}