图像实验分析和总结.docxVIP

实验一 图像像素处理 提取图像 woman 的像素在 之间的像素点，并对提取处理的像素作如下处理： 对其提取的像素求平均值且将重新赋值于原像素； 将提取的像素全部置为 0； 将提取的像素全部置为 255； 要求显示处理前后的图像（4 幅图显示在同一窗口，且图像标题分别为“原图”、“置 0 图”、“置 255 图”、“置平均值图”）。 clc clear all load woman T=X; subplot(2,2,1) imshow(X,map) title(原图) b=find(X=50X=100); X(b)=mean(X(b)); subplot(2,2,2) imshow(X,map) title(置平均值图) X=T; b=find(X=50X=100); X(b)=0; subplot(2,2,3) imshow(X,map) title(置 0 图) X=T; b=find(X=50X=100); X(b)=255; subplot(2,2,4) imshow(X,map) title(置 255 图) 对图像 detfingr 进行阈值处理，阈值函数如下： ，要求显示处理前后的图像（2 幅图显示在同一窗口，且图像标题分别为“处理前图像”、“处理后图像”）。 clc clear all load detfingr subplot(1,2,1) imshow(X,map) title(处理前图像) b=find(X=127); X(b)=0; c=find(X127); X(c)=255; subplot(1,2,2) imshow(X,map) title(处理后图像) 实验二 图像灰度处理 1、分别对图像 woman 进行如下灰度变换 反转变换 灰度级数 8bit，函数取s ? L ?1? r ，要求显示处理前后的图像（显示在同一窗口）； 对数变换 函数取 s ? log(1? r) ，要求显示处理前后的图像（显示在同一窗口）； 幂次变换 函数取 s ? r? （ 分别取 0.9 和 1.1），要求显示处理前后的图像（显示在同一窗口）。 clc clear all load woman [m,n]=size(X) subplot(3,2,1) imshow(X,map) title( 原 图 ) L=256; T=X; for i=1:1:m for j=1:1:n X(i,j)=L-1-X(i,j); end end subplot(3,2,2) imshow(X,map) title(反转变换) X=T; for i=1:1:m for j=1:1:n X(i,j)=log(1+X(i,j)); end end subplot(3,2,3) imshow(X,map) title(对数变换) X=T; for i=1:1:m for j=1:1:n X(i,j)=X(i,j)^(0.9); end end subplot(3,2,4) imshow(X,map) title(幂次变换 1) for i=1:1:m for j=1:1:n X(i,j)=X(i,j)^(1.1); end end subplot(3,2,5) imshow(X,map) title(幂次变换 2) 实验三 图像直方图处理 对图像 woman 进行灰度级（即 8bit 转换为 3bit），并对转换后的图像进行直方图均衡化处理、且分别显示出每个灰度级的像素个数； clc load woman subplot(1,3,1) imshow(X,map) X=round(7*X/255); subplot(1,3,2) imshow(X,map) [m,n]=size(X); MN=m*n; for i=0:7 ss=find(X==i); l=i+1; nk(l)=length(ss); end prrk=nk/MN; sum=0; for i=1:8 sum=7*prrk(i)+sum; s(i)=round(sum); end for i=0:7 if s(i+1)==i X(find(X==i))=i; else end  end X(find(X==i))=s(i+1); subplot(1,3,3) imshow(X,map) 实验 4 1、 对图像 woman 进行均值滤波，自己编写均值滤波函数，用 3*3 模板。（注：不得用 matlab 中的滤波函数） clc clear all load woman T=X; subplot(121) imshow(X,map) title( 原 图 ) [m,n]=size(X); for i=2:m-1 for j=2:n-1 Y(i,j)=round((X(i-1,j-1