ztong上机3

  

二、实验名称：数字图像处理matlab上机

三、实验学时：2学时

四、实验目的：（详细填写）

• 掌握几何变换
• 掌握插值
• 理解配准的概念

五．实验内容

• （1）首先自己写一个对图像进行旋转和缩放的复合变换程序；
• （2）然后通过基于点的配准找到一个变换矩阵，以将前面变换后的图像再变换回去。

六、实验器材（设备、元器件）：

硬件：计算机 软件：windows XP ,matlab7.0

七、实验数据及结果分析：

内容一程序：

function out = LinearT(
In,theta,sx, sy)
im
= imread(In);
degree = theta;
imm = im2double(im);
cos_d = cos(degree);
sin_d = sin(degree);
[width, height] = size(im);
R
= [cos_d -sin_d 0
sin_d cos_d 0
0
0 1];
T
= [1 0 -width/2.0
0
1 -height/2.0
0
0 1];
S=[sx,0,0;
0,sy,0;
0,0,1];
R=inv(T)*S*R*T;
[W, H, ws, we, hs, he] =
ComputeOutImSize(im, R);

inR = inv(R);
ii=0;

for i = ws:we%row
ii=ii+1; jj=0;
for j = hs:he%col
jj=jj+1;
locate_r = [i j 1]';
locate_s = inR*locate_r;
if(((locate_s(1)<=1)||(locate_s(1)>=width))||((locate_s(2)<=1)||(locate_s(2)>=height)))
V(ii,jj) = 0;
else
V(ii,jj) =
BilinearInterplate(imm,locate_s);

end
end
end

out= im2uint8(V);
end

调用程序
In='pre.bmp';

theta=pi/8;
sx=2;

sy=4;

out = LinearT( In,theta,sx,
sy);
imshow(out)

原始图像

  

处理后的图像

  

  

  

内容二程序

s=1;

f1=imread('pre.bmp');

f2=imread('flop.bmp');

figure(1)
subplot(121),imshow(f1)
subplot(122),imshow(f2)
N=4;
imtool(f2)
[w1,h1]=size(f1);
[w2,h2]=size(f2);
a=[0 0]';b=[w1-1 0]';c=[0 h1-1]';d=[w1-1
h1-1]';
a1=[4 93]';b1=[224 1]';c1=[113
359]';d1=[333 267]';
x=[a b c d];
y=[a1 b1 c1 d1];
H=zeros(2);
xmean=[sum(x(1,:))/N;sum(x(2,:))/N];

ymean=[sum(y(1,:))/N;sum(y(2,:))/N];

xx=[x(1,:)-xmean(1);x(2,:)-xmean(2)];

yy=[y(1,:)-ymean(1);y(2,:)-ymean(2)];

for i=1:N
H=H+xx(:,i)*(yy(:,i))';
end
[U,S,V]=svd(H);
r=det(U*V);
W=[1 0;0 r];
R=V*W*U';
s=0;

for i=1:N
s1=s+R*xx(:,i).*yy(:,i);
end

for i=1:N
s2=s+R*xx(:,i).*xx(:,i);
end

s=s1/s2;
t=xmean-s*R*ymean;
imm=im2double(f2);
ii=0;
for i
=1:w1
ii=ii+1; jj=0;
for j =
1:h1
jj=jj+1;
SS=inv(R)*([i j]'-t);
if(((SS(1)<=1)||(SS(1)>=w2))||((SS(2)<=1)||(SS(2)>=h2)))

S(ii,jj) = 0;
else
S(ii,jj) = BilinearInterplate(imm,SS);

end
end

end

S0=
im2uint8(S);
figure(2)
imshow(S0)

  

  

  

八、总结及心得体会：

• 掌握了几何变换
• 掌握了插值
• 理解了配准的概念

  

  

报告评分：

指导教师签字：