006X11图像的增强实验(程序)(精)

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); %将均衡化后图像的直方图显示%自主编程的直方图均衡化处理m,

2、n=size(I1);%计算图像数据矩阵大小p=zeros(1,256);%创建用于对图像灰度值计数的数组p for l=1:m %行循环lfor r=1:n %列循环rh=I1(l,r);%灰度值计数p(h+1)=p(h+1)+1;%灰度值计数endendR=p/(m*n);%将灰度级的概率分布存入数组Rsum=zeros(1,256);%创建累计直方图灰度级分布律存储数组sums=0;for c=1:256 %统计累计直方图灰度级分布律sum(c)=s+R(c);s=sum(c);end%更改每个需处理的像素的灰度值I3=I1;for c=1:256d=0;SS=round(sum(c)*

3、255);%映射后的灰度值数组RR=round(R(c)*255);if (abs(sum(c)-R(c)=0.05d=1;endif(d=0) %判断是否更改原灰度级for l=1:mfor r=1:nif(RR=I3(l,r) %查找出需更改灰度级的像素点I3(l,r)=SS-1; %更改灰度级endendendendend subplot(2,3,3);imshow(I3);title(自编直方图均衡化处理);q=zeros(1,256);%创建用于对自编直方图处理后图像灰度值计数的数组q forl=1:m %行循环lfor r=1:n %列循环rh=I3(l,r);%灰度值计数q(h+

4、1)=q(h+1)+1;%灰度值计数endendx=1:256;subplot(2,3,6);plot(x-1,q); %将自编程序均衡化处理后图像的直方图显示滤波实验clearall;I=imread(D:pcdzx1.jpg);%读入原图像文件subplot(2,4,1);imshow(I);%显示原图像title(原图像);I1=rgb2gray(I);%转变成灰度图像subplot(2,4,2);imshow(I1);%显示灰度图像title(灰度图像);I2=imnoise(I1,gaussian,0.1,0.002);subplot(2,4,3);imshow(I2);%显示加入高

5、斯噪声之后的图象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的矩阵I3=x.*x1; %模板上系数与模板下系数对应相乘I3=I3(1,:),I3(2,:),I3(3,:);a=uint8(sum(I3); %将所有乘积相加I4(l,r)=a; %取和值赋给所在元素endendsubplot(2,4,5);imshow(I4);%显示3X3均值滤波后的

6、图像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,:);%转换为1X9的矩阵I7=sort(I7); %按由小到大顺序排列I5(l,r)=I7(5); %取和值赋给所在元素endendsubplot(2,4,6);imshow(I5);%显示3X3中值滤波后的图像title(3X3中值滤波后的图像);%5X5均值滤波x=double(ones(5)/25);I8=I2;for l=3:(m-2)fo

7、r r=3:(n-2)x1=double(I8(l-2:l+2,r-2:r+2); %从噪声图像中取出5X5的矩阵I3=x.*x1; %模板上系数与模板下系数对应相乘I3=I3(1,:),I3(2,:),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(l7); %按由小到大顺序排列I9(l,r)=I7(13); %取和值赋给所在元素endendsubplot(2,4,8);imshow(l9);%显示5X5中值滤波后的图像title(5X5中值滤波后的图像);图像的增强实验滤波实验匣图蟻抿度阳牌加鸟,离涉嗓走之用的西第熬 3 均也滤漑帝的圏像和 3 中也溥澈帝的图像5均 0韓械拓的