数字图像实验二：几何变换与变形c++实现

实验2.1：图像缩放

实现一个图像缩放函数，可以对输入图像进行任意倍数的缩放；
采用双线性插值进行重采样；
X,Y方向的缩放倍数参函数参数的形式传入；
可以只考虑输入图像为3通道，8位深度的情况；
不能调用图像处理库的缩放函数来完成；
实验2.2：图像变形
记[x’, y’]=f([x, y])为像素坐标的一个映射，实现f所表示的图像形变。f的逆映射为：

2.1代码：

#include<opencv2/core/core.hpp>
#include<opencv2/highgui/highgui.hpp>
#include<iostream>
using namespace cv;
using namespace std;

Mat scale(Mat SImage, double sx, double sy)
{
    int x = (int)(SImage.rows*sx);
    int y = (int)(SImage.cols*sy);
    Mat DImage(x, y, CV_8UC3);
        for (int i = 0; i < DImage.rows; i++)
        {
            for (int j = 0; j < DImage.cols; j++)
            {
                int tempx = (int)(i / sx);
                int tempy = (int)(j / sy);
                if (tempx >= SImage.rows - 1 || tempy >= SImage.cols - 1)
                {
                    DImage.ptr(i, j)[0] = SImage.ptr(SImage.rows - 1, SImage.cols - 1)[0];
                    DImage.ptr(i, j)[1] = SImage.ptr(SImage.rows - 1, SImage.cols - 1)[1];
                    DImage.ptr(i, j)[2] = SImage.ptr(SImage.rows - 1, SImage.cols - 1)[2];
                }
                else { 
                    double b = i / sx - tempx;
                    double a = j / sy - tempy;
                    DImage.ptr(i, j)[0] = SImage.ptr(tempx, tempy)[0] * (1 - a)*(1 - b) + SImage.ptr(tempx, tempy + 1)[0] * a*(1 - b)
                        + SImage.ptr(tempx + 1, tempy)[0] * b*(1 - a) + SImage.ptr(tempx + 1, tempy + 1)[0] * a*b;
                    DImage.ptr(i, j)[1] = SImage.ptr(tempx, tempy)[1] * (1 - a)*(1 - b) + SImage.ptr(tempx, tempy + 1)[1] * a*(1 - b)
                        + SImage.ptr(tempx + 1, tempy)[1] * b*(1 - a) + SImage.ptr(tempx + 1, tempy + 1)[1] * a*b;
                    DImage.ptr(i, j)[2] = SImage.ptr(tempx, tempy)[2] * (1 - a)*(1 - b) + SImage.ptr(tempx, tempy + 1)[2] * a*(1 - b)
                        + SImage.ptr(tempx + 1, tempy)[2] * b*(1 - a) + SImage.ptr(tempx + 1, tempy + 1)[2] * a*b;
                }               
            }
        }
    return DImage;
}


int main()
{
    Mat SImg = imread("god.jpg");
    Mat DImg = scale(SImg, 1.5, 2.3);
    namedWindow("原图");
    imshow("原图", SImg);
    namedWindow("缩放后");
    imshow("缩放后", DImg);
    waitKey();

}

2.2代码：

#include<opencv2/core/core.hpp>
#include<opencv2/highgui/highgui.hpp>
#include<iostream>
#include<math.h>

using namespace cv;
using namespace std;

int main()
{
    Mat img = imread("god.jpg");
    Mat dest(img.rows, img.cols, CV_8UC3);
    int x = dest.rows;
    int y = dest.cols;
    for (int i = 0; i < x; i++)
    {
        for (int j = 0; j < y; j++)
        {
            double tempx = i / (0.5*x) - 1;
            double tempy = j / (0.5*y) - 1;
            double srcx;
            double srcy;
            double r = sqrt(pow(tempx, 2.0) + pow(tempy, 2.0));
            double angle = pow(1 - r, 2.0);
            if (r>=1)
            {
                srcx = tempx;
                srcy = tempy;
            }
            else {
                    srcx = cos(angle)*tempx - sin(angle)*tempy;
                    srcy = sin(angle)*tempx + cos(angle)*tempy;
            }
            int sx = (srcx + 1)*0.5*x;
            int sy = (srcy + 1)*0.5*y;
            dest.ptr(i, j)[0] = img.ptr(sx, sy)[0];
            dest.ptr(i, j)[1] = img.ptr(sx, sy)[1];
            dest.ptr(i, j)[2] = img.ptr(sx, sy)[2];
        }
    }
    imshow("原图", img);
    imshow("变形后", dest);
    waitKey();
}

注意：运行代码前你需要配置好opencv。