实验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。