数字图像处理实验二告.杨海堂(U201012943)docx.docxVIP

数字图像处理实验报告 次数：第二次 姓名：杨海堂 班级：电信1006 学号：U201012943 指导老师：曹汉强 图像拼接 一．实验要求： 试从Fig1.bmp中提取出目标物体（飞机），并将其与背景图像background.bmp进行合成。 实验目的： 11熟悉并学会使用MATLAB中图像增强的相关函数 1掌握数字图像去除图像背景的方法及图像合成。 实验原理： 将飞机预想从原图中剥离出来，运用像素的差别进行判断，把属于飞机的像素部分保留下来，其余的像素置零（把背景变成黑色）。建立一个与背景图像大小相同的空图像zeros（）。先把截取的飞机图像复制进空图像中的相应位置，在把大小相同的图像进行覆盖，比较像素值把飞机周围的黑色背景用back.bmp中的像素填充。 流程图： 实验代码： f=imread(plane.png); figure,imshow(f,[]);%显示飞机原图 X=40;Y=50; fR=double(f(:,:,1)); fG=double(f(:,:,2)); fB=double(f(:,:,3)); for i=1:165 for j=1:273 if abs(fR(i,j)-65)40abs(fG(i,j)-60)40abs(fB(i,j)-200)=55 fR(i,j)=0; fG(i,j)=0; fB(i,j)=0; end end end for i=1:165 for j=1:273 if i160||j267 fR(i,j)=0; fG(i,j)=0; fB(i,j)=0; end end end f1=fR(1:2:end,1:2:end); f2=fG(1:2:end,1:2:end); f3=fB(1:2:end,1:2:end); ff(:,:,1)=f1; ff(:,:,2)=f2; ff(:,:,3)=f3; ff=mat2gray(ff); figure,imshow(ff);%得到去除背景的飞机图 g=imread(background.png); g=mat2gray(double(g)); figure,imshow(g);%显示彩色原图 ff1=zeros(210,274,3); ff1(X:X+82,Y:Y+136,1)=f1; ff1(X:X+82,Y:Y+136,2)=f2; ff1(X:X+82,Y:Y+136,3)=f3; ff1=mat2gray(ff1); figure,imshow(ff1);%对得到的背景图进行加工处理 for i=1:210 for j=1:274 for k=1:3 if ff1(i,j,k)==0 ff1(i,j,k)=g(i,j,k); end end end end figure,imshow(ff1,[]);%得到合成图片 显示结果图： 实验小结： 最后两张图片合成时候，应确保两张图片大小一致。在进行图像拼接的时候一定要注意两图像边缘的像素处理问题，需要进行平滑处理，否则拼接的图像边缘会出现明显的颜色失真。 二 Fig2.bmp是由灰度图像orig.bmp（256×256）经过空间线性变换得到的图像。试对Fig2.bmp进行几何校正，恢复出原始图像。 图 11经过线性变换扭曲的图像 关键代码设计思想 在数字图像处理中，几何变换由两个基本操作组成：（1）一个空间变换，它定义了图像平面上像素如何进行重新排列；（2）灰度级插补，它处理空间变换前后图像中像素灰度级的赋值。 空间变换设计思想 假设一幅图像f（x，y）经过几何失真后得到几何失真图像g（u，v），那么设 u=c1x+c2y+c3xy+c4 v=c5x+c6y+c7xy+c8 这时若已知原图像与目标图像中对应的四个点，代入上述式子当中，便可以求得c1至c8这些系数了。这样原图像到目标图像像素点便建立了映射的关系；同理，假设现在只知道目标图像，我们也可以反过来求得目标图像到原图的映射关系，再利用最近领域插值或双线性插值就可以求得原图像对应点的像素值。 灰度级插补设计思想 在图像缩放一节里我们已经接触到了最近领域插值与双线性插值两种插值方法，其实它们就是在做灰度级插补的工作，而缩放变换也是几何变换一个特例。这里我采取双线性灰度级插补。 流程图 实验代码： c