基于matlab哈哈镜制作.pdfVIP

计算机视觉课程设计实验报告

一、题目：哈哈镜的制作

三、实验目的：

1、掌握几何变换的原理。

2、读入待制作的图片，实现图片的缩小放大、水平外凹、水平外凸、梯形变形、三角形S

形和波浪型。

四、实验原理：

变形处理整体都是一个思路，把原图每一行经过缩放拷到新图的相应行。

每次变形都是在原图的基础上，根据想要变换的图象形状，计算出每一行的起始位，再计算

出每一行新的总长度。

1）对于放大、缩小可以根据最近邻插值或者双线性插值来进行变化，放大我们采用了双线

性方法，缩小采用了最近邻插值方法。

2）对于水平内凹和水平外凸，以周期为2倍高度的正弦函数来进行变化，变化后远图象的

高度坐标不变，宽度坐标以一定比例变化。

3）对于三角形和梯形，原图象高度坐标不变，宽度坐标以线性来进行变化。

4）对于宽度不变的变形，如波浪形，直接一次变形拷贝一整行。对于高度不变的图形，如

S形，也可以用这个方法，或者说是把起始点确定后，将坐标点进行平移。

5）除放大效果外，缩小、水平外凹、水平外凸、梯形变形、三角形S形和波浪型均采用

最近邻插值方法进行变化。

五、实验步骤（附程序主代码）：

1、打开图象

clear;clc;%清空工作空间

[filename,pathname]=uigetfile({*.*,AllFiles(*.*)});

file=fullfile(pathname,filename);

t1=imread(file);

t2=rgb2gray(t1);%将打开的彩色图像转换为灰度图

imshow(t2)

savetempt2;

2、缩小

loadtempt2%加载打开的图片

[m,n]=size(t2);

x2=double(t2);

fori=1:m

forj=1:n

w=round(i*0.5);%宽度缩小为原来一半

h=round(j*0.5);%高度缩小为原来一半

x3(w,h)=x2(i,j);%最近邻插值，象素值赋给新的图象

end

end

figure,imshow(mat2gray(x3))

3、放大

x3=zeros(fix(m*2),fix(n*2));%开辟空间

fori=2:fix((m-1)*2)

forj=2:fix((n-1)*2)

p=fix(i/2);

q=fix(j/2);

x3(i,j)=[x2(p+1,q)-x2(p,q)]*(i/2-fix(i/2))+[x2(p,q+1)-x2(p,q)]*(j/2-

fix(j/2))+[x2(p+1,q+1)+x2(p,q)-x2(p,q+1)-x2(p+1,q)]*(i/2-fix(i/2))*

(j/2-fix(j/2))+x2(p,q);%双线性插值

end

end

4、水平内凹

t3(m,n)=0;%开辟新图象空间

fori=1:m

a=(2*pi/m)*i;%i转化为弧度形式

p=round(abs(30*sin(a/2))+0.5);

q=round(-abs(30*sin(a/2))-0.5);%正弦函数变换，周期为2M

forj=p:q+n

x=i;%变化后高度坐标不变

y=round((j-p)*n/((q+n-p+1)));%宽度坐标以(j-p)*n/((q+n-p+1)))比例变化

if(x0x=m)(y0y=n)

t3(i,j)=t2(x,y);

else

t3(i,j)=255;

end

end

end

5、水平外凸

a=round(sin(pi/2));

t3(m,n+2*a)=0;

grade=50;

fori=1:m

b=2*pi/m*i;%i转化为弧度形式

j1=round(abs(sin(2*pi-b/2)*grade+grade)+0.5);

j2=round(sin(b/2)*grade-grade-0.5);%正弦函数变换，周期为2M

f