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

《《数字图像处理》实验报告02》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、GDOU-B-11-112广东海洋大学学生实验报告书（学生用表）实验名称图像灰度变换与直方图修正课程名称数字图像处理课程号19242504学院(系)理学院数信系专业信息与计算科学班级1112学生姓名麦杰洪学号201111921215实验地点钟05008实验日期13.10.11一、实验目的1.能熟练运用Matlab对图像进行灰度变换；2.能熟练运用Matlab对图像进行直方图修正；3.能熟练运用Maltab对图像进行空域平滑去噪。二、实验内容1.对图像’office_2.jpg’分别进行灰度线性变换、对数变换、幂次变换，使图像更清晰；2.对图像’office_2.jpg’进行直方图均衡化；

2、以下3题任选1个3.自编一个直方图均衡化的函数，输入输出格式：functionJ=myhisteq(I)并用自编的函数对‘tire.tif’进行直方图均衡化；4.对图像’coins.png’添加高斯噪声，自编代码用掩膜平滑法实现去噪；5.对图像’coins.png’添加椒盐噪声，自编代码用中值滤波法实现去噪。三、程序清单与运行结果1.灰度线性变换clearall,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

3、,[s(1),s(2)]);G=reshape(G,[s(1),s(2)]);B=reshape(B,[s(1),s(2)]);编写程序：functionJ=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;fori=1:Mforj=1:Nif(I(i,j)>=a)&(I(i,j)<=b)f(i,j)=((d-c)/(b-a))*(I(i,j)-a)+c;endif(I(

4、i,j))f(i,j)=c;endif(I(i,j)>=b)f(i,j)=d;endendendfigure;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.灰度对数变换clearall,clc;im=imread('

5、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(R<30);r(idx1)=R(idx1)*4;idx2=find(R>=31);k=(

6、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(122),imshow(I);3.灰度

7、幂次变换分别对红、绿、蓝三种颜色进行处理：红色部分：k=0.5;r=255*(double(R)/255).^k;r=uint8(r);figuresubplot(221),imshow(R);subplot(222),imhist(R);subplot(223),imshow(r);subplot(224),imhist(r);再分别对绿色和蓝色进行处理：最后把三个部分组合起来：I=cat(3,r,g,b);I=uint8(I);