《006X11图像的增强实验(程序)(精).doc》由会员分享，可在线阅读，更多相关《006X11图像的增强实验(程序)(精).doc（3页珍藏版）》请在三一文库上搜索。
1、图像的增强实验 clear all; I=imread(D:pcpic.jpg);% 读入原图像文件 subplot(2,3,1); imshow(I);% 显示原图像 title( 原图像 ); I1=rgb2gray(I);% 转变成灰度图像 subplot(2,3,4); imshow(I1);% 显示灰度图像 title( 灰度图像 ); I2=histeq(I1)% 对灰度图像进行直方图均衡化处理 subplot(2,3,2); imshow(I2);% 显示直方图均衡化处理后的图像 title( 直方图均衡化处理 ); subplot(2,3,5); Imhist(I2); % 将
2、均衡化后图像的直方图显示 %自主编程的直方图均衡化处理 m,n=size(I1);% 计算图像数据矩阵大小 p=zeros(1,256);% 创建用于对图像灰度值计数的数组 p for l=1:m % 行循环 l for r=1:n % 列循环 r h=I1(l,r);% 灰度值计数 p(h+1)=p(h+1)+1;% 灰度值计数 end end R=p/(m*n);% 将灰度级的概率分布存入数组 R sum=zeros(1,256);% 创建累计直方图灰度级分布律存储数组 sum s=0; for c=1:256 % 统计累计直方图灰度级分布律 sum(c)=s+R(c); s=sum(c)
3、; end %更改每个需处理的像素的灰度值 I3=I1; for c=1:256 d=0; SS=round(sum(c)*255);% 映射后的灰度值数组 RR=round(R(c)*255); if (abs(sum(c)-R(c)=0.05 d=1; end if(d=0) % 判断是否更改原灰度级 for l=1:m for r=1:n if(RR=I3(l,r) % 查找出需更改灰度级的像素点 I3(l,r)=SS-1; % 更改灰度级 end end end end end subplot(2,3,3); imshow(I3); title( 自编直方图均衡化处理 ); q=zer
4、os(1,256);% 创建用于对自编直方图处理后图像灰度值计数的数组 q for l=1:m % 行循环 l for r=1:n % 列循环 r h=I3(l,r);% 灰度值计数 q(h+1)=q(h+1)+1;% 灰度值计数 end end x=1:256; subplot(2,3,6); plot(x-1,q); % 将自编程序均衡化处理后图像的直方图显示 滤波实验 clear all; I=imread(D:pcdzx1.jpg);% 读入原图像文件 subplot(2,4,1); imshow(I);% 显示原图像 title( 原图像 ); I1=rgb2gray(I);% 转变
5、成灰度图像 subplot(2,4,2); imshow(I1);% 显示灰度图像 title( 灰度图像 ); I2=imnoise(I1,gaussian,0.1,0.002); subplot(2,4,3); imshow(I2);% 显示加入高斯噪声之后的图象 title( 加入高斯噪声之后的图象 ); m,n=size(I1);% 计算图像数据矩阵大小 %3X3 均值滤波 x=double(ones(3)/9); I4=I2; for l=2:(m-1) for r=2:(n-1) x1=double(I4(l-1:l+1,r-1:r+1); % 从噪声图像中取出 3X3 的矩阵 I
6、3=x.*x1; % 模板上系数与模板下系数对应相乘 I3=I3(1,:),I3(2,:),I3(3,:); a=uint8(sum(I3); % 将所有乘积相加 I4(l,r)=a; % 取和值赋给所在元素 end end subplot(2,4,5); imshow(I4);% 显示 3X3 均值滤波后的图像 title(3X3 均值滤波后的图像 ); %3X3 中值滤波 I5=I2; for l=2:(m-1) for r=2:(n-1) I6=I5(l-1:l+1,r-1:r+1); % 从噪声图像中取出 3X3 的矩阵 I7=I6(1,:),I6(2,:),I6(3,:);% 转换为
7、 1X9 的矩阵 I7=sort(I7); % 按由小到大顺序排列 I5(l,r)=I7(5); % 取和值赋给所在元素 end end subplot(2,4,6); imshow(I5);% 显示 3X3 中值滤波后的图像 title(3X3 中值滤波后的图像 ); %5X5 均值滤波 x=double(ones(5)/25); I8=I2; for l=3:(m-2) for r=3:(n-2) x1=double(I8(l-2:l+2,r-2:r+2); % 从噪声图像中取出 5X5 的矩阵 I3=x.*x1; % 模板上系数与模板下系数对应相乘 I3=I3(1,:),I3(2,:),
8、I3(3,:),I3(4,:),I3(5,:); a=uint8(sum(I3); % 将所有乘积相加 I8(l,r)=a; % 取和值赋给所在元素 end end subplot(2,4,7); imshow(I8);% 显示 5X5 均值滤波后的图像 title(5X5 均值滤波后的图像 ); %5X5 中值滤波 I9=I2; for l=3:(m-2) for r=3:(n-2) l6=l9(l-2:l+2,r-2:r+2); %从噪声图像中取出 3X3的矩阵 17=16(1,:),16(2,:),16(3,:),16(4,:),16(5,:);% 转换为 1X9 的矩阵 I7=sort
9、(l7); %按由小到大顺序排列 I9(l,r)=I7(13); %取和值赋给所在元素 end end subplot(2,4,8); imshow(l9);%显示5X5中值滤波后的图像 title(5X5中值滤波后的图像); 图像的增强实验 滤波实验 匣图蟻 抿度阳牌 加鸟,离涉嗓走之用的西第 熬 3 均也滤漑帝的圏像 和 3 中也溥澈帝的图像 5均 0韓械拓的囲像 5X5 中谊滤械启的囹愫 自購宜方蹲均働化加理 直方图均命上处煙 :* 1： 冶；1 fep i iksTF*- I jjii-Fi-tj-Fin-：匸可： J 7 I灯蛊=1 ! LJtaMd li 严六 lie b _ -EZ Id J 計 423*-