资源描述：
《数字图像处理实验报告(空间域图像增强)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、实验报告实验名称空间域图像增强课程名称数字图像处理姓名成绩班级学号日期地点备注：1.实验目的（1）了解空间域图像增强的各种方法（点处理、掩模处理）；（2）通过编写程序掌握采用直方图均衡化进行图像增强的方法；（3）使用邻域平均法编写程序实现图像增强，进一步掌握掩模法及其改进（加门限法）消除噪声的原理；（4）总结实验过程（实验报告，左侧装订）：方案、编程、调试、结果、分析、结论。2.实验环境（软件、硬件及条件）Windws7MATLAB6.xorabove3.实验方法对如图4.1所示的两幅128×128、256级灰度的数字图像fing_128.img和cell
2、_128.img进行如下处理：（1）对原图像进行直方图均衡化处理，同屏显示处理前后图像及其直方图，比较异同，并回答为什么数字图像均衡化后其直方图并非完全均匀分布。（2）对原图像加入点噪声，用4-邻域平均法平滑加噪声图像（图像四周边界不处理，下同），同屏显示原图像、加噪声图像和处理后的图像。①不加门限；②加门限T=(1/2)*avg(f(m,n)),其中avg（f(m,n)=(1/N^2)*f(i,j)）本次实验中的第一题，是对图像进行直方图统计和均衡化，在Matlab中有imhist()函数和histeq()函数直接调即可获得相应结果，代码如下：close
3、all;clearall;fid=fopen('cell_128.img','r');image1=fread(fid,[128,128],'uint8');image1=uint8(image1);fclose(fid);subplot(2,2,1);%显示原图像imshow(image1,[]);title('原图像');subplot(2,2,2);%统计图像直方图imhist(image1);title('原图像直方图');%直方图均衡化ima=histeq(image1);subplot(2,2,3);%显示均衡化后的图像imshow(ima);
4、title('直方图均衡化后图像');subplot(2,2,4);%显示直方图均衡化后的图像直方图imhist(ima);title('直方图均衡化后的直方图');本次实验的第二题，是对图像进行加噪后平滑，噪声为点噪声，分别采用不加门限的四邻域平均法和加门限的四邻域平均法进行平滑。主要思想为：对于不加门限的四邻域平均法，对加噪后图像的每一像素（边界不考虑）的灰度值用其周围（上下左右）四个像素的灰度值的平均值代替；对于加门限的四邻域平均法，是对先对图像每一像素的灰度值进行判别，若该像素灰度值与四周（上下左右）四个像素灰度值的平均值之差大于门限T，则用四周四
5、个像素灰度的平均值代替，否则保持该像素灰度值不变。对图像分别加如随机点噪声，并分别对其进行平滑，实验代码如下：closeall;clearall;fid=fopen('fing_128.img','r');image2=fread(fid,[128,128]);%image2=uint8(image2);fclose(fid);subplot(2,2,1);imshow(uint8(image2),[]);title('原图像');%为图像添加噪声snoise=0.1*randn(size(image2));J21=imadd(image2,snoise)
6、;%添加随机点噪声tmp21=J21;subplot(2,2,2);imshow(uint8(J21),[]);title('加随机点噪声');%对加随机点噪声的图像进行不加门限的四邻域法平滑tmp21=J21;fori=2:(size(J21,1)-1)forj=2:(size(J21,2)-1)tmp21(i,j)=(J21((i-1),j)+J21((i+1),j)+J21(i,(j-1))+J21(i,(j+1)))/4;endendsubplot(2,2,3);imshow(uint8(tmp21),[]);title('随机噪声不加门限');%
7、对加随机点噪声的图像进行加门限的四邻域法平滑N=128;f=J21;forx=1:size(J21,1)fory=1:size(J21,2)f=(1/N.*N)*(f+J21(x,y));endendT=(1/2)*f;tmp23=J21;fori=2:(size(J21,1)-1)forj=2:(size(J21,2)-1)tmp23(i,j)=(J21((i-1),j)+J21((i+1),j)+J21(i,(j-1))+J21(i,(j+1)))/4;if(J21(i,j)-tmp23(i,j))>Ttmp23(i,j)=tmp23(i,j);els
8、etmp23(i,j)=J21(i,j);endendendsub
