实验说明

代码
%% prepare
% 获取图片
img = rgb2gray(imread('lena.jpg'));
width = size(img,2);
height = size(img,1);
% 创建新窗口,显示图片
figure;
imshow(img);
% 选取角点
moving_points = ginput(4);
hold on ;
plot(moving_points(:,1),moving_points(:,2),'ro');
% fixed_points矩阵
fixed_points = [0,0;
100,0;
0,200;
100,200];
%% use matlab function
% 计算变化矩阵tfom,获取相似变换对象
tfom = fitgeotrans(moving_points,fixed_points,'projective');
% 前置变换
X = moving_points(:,1);
Y = moving_points(:,2);
[x,y] = transformPointsForward(tfom,X(:),Y(:));
% 创建新窗口,显示网格
figure;
plot(x,y,'ro');
title('验证坐标点对齐')
grid on % 打开网格
% 计算所用的时间
tic; % 记录起始时刻
dst_img = imwarp(img,tfom); % 将几何变换应用于图像
t_sys = toc; % 计算时间差
% 创建新窗口,显示变换后的图片
figure;
imshow(dst_img);
title(['图像仿射变换后(系统函数),耗时(s):',num2str(t_sys)]) % 转换成字符串,输出时间
结果
