《数字图像处理》实验报告02.docVIP

广东海洋大学学生实验报告书（学生用表） 实验名称 图像灰度变换与直方图修正 课程名称 数字图像处理 课程号 学院(系) 理学院数信系 专业 信息与计算科学 班级 1112 学生姓名 麦杰洪 学号 201111921215 实验地点 钟05008 实验日期 13.10.11 一、实验目的 1.能熟练运用Matlab对图像进行灰度变换； 2.能熟练运用Matlab对图像进行直方图修正； 3.能熟练运用Maltab对图像进行空域平滑去噪。 二、实验内容 1.对图像’office_2.jpg’分别进行灰度线性变换、对数变换、幂次变换，使图像更清晰； 2.对图像’office_2.jpg’进行直方图均衡化； 以下3题任选1个 3.自编一个直方图均衡化的函数，输入输出格式：function J=myhisteq(I) 并用自编的函数对‘tire.tif’进行直方图均衡化； 4.对图像’coins.png’添加高斯噪声，自编代码用掩膜平滑法实现去噪； 5.对图像’coins.png’添加椒盐噪声，自编代码用中值滤波法实现去噪。 三、程序清单与运行结果 1. 灰度线性变换 clear all,clc; im=imread(office_2.jpg);% office_2.jpg为图片，在m文件所在目录下 s=size(im); R=im(:,:,1); G=im(:,:,2); f=im(:,:,3); R=reshape(R,[s(1),s(2)]); G=reshape(G,[s(1),s(2)]); B=reshape(B,[s(1),s(2)]); 编写程序： function J=my(I) [M,N]=size(I); figure,imhist(I); j=zeros(M,N); I=double(I); f=double(f); k1=min(min(I)); k2=max(max(I)); a=k1+50; b=k2-50; c=k1-30; d=k2+20; for i=1:M for j=1:N if (I(i,j)=a)(I(i,j)=b) f(i,j)=((d-c)/(b-a))*(I(i,j)-a)+c; end if (I(i,j)) f(i,j)=c; end if(I(i,j)=b) f(i,j)=d; end end end figure; subplot(121); imshow(I,[]); subplot(122); imshow(f,[]); 先对红色部分进行处理： J=my(R) r=f; 再对绿色部分进行处理： J=my(G) g=r; 最后对蓝色部分进行处理： J=my(B) b=f; 再把三个部分组合起来： I=cat(3,r,g,b); I=uint8(I); figure,subplot(121),imshow(im); subplot(122),imshow(I); 2. 灰度对数变换 clear all,clc; im=imread(office_2.jpg);% office_2.jpg为图片，在m文件所在目录下 s=size(im); R=im(:,:,1); G=im(:,:,2); f=im(:,:,3); R=reshape(R,[s(1),s(2)]); G=reshape(G,[s(1),s(2)]); B=reshape(B,[s(1),s(2)]); 先对红色部分进行处理： figure, imhist(R); %观察直方图，确定调整方案 %以下作分段线性变换 R=double(R); r=zeros(size(R)); idx1=find(R30); r(idx1)=R(idx1)*4; idx2=find(R=31); k=(255-121)/(255-31); r(idx2)=121+k*(R(idx2)-31); r=uint8(r); R=uint8(R); %还原成无符号8位整数 figure, subplot(2,2,1),imshow(R),subplot(2,2,2),imhist(R) subplot(2,2,3),imshow(r),subplot(2,2,4),imhist(r) 再分别对绿色部分和蓝色部分进行处理, 最后把三个部分组合起来： I=cat(3,r,g,b); I=uint8(I); figure,subplot(121),imshow(im); subplot(12