实验二-空间域图像增强

实验二空间域图像增强一、实验目的1、了解空间域图像增强的各种方法〔点处理、掩模处理〕；2、通过编写程序掌握采用直方图均衡化进行图像增强的方法；3、使用邻域平均法编写程序实现图像增强，进一步掌握掩模法及其改良〔加门限法〕消除噪声的原理；二、实验环境及开发工具Windws2000/XP、MATLAB6.x、VisualC++、VisualBasic或其它三、实验方法对如下图的两幅128×128、256级灰度的数字图像fing_128.img和cell_128.img进行如下处理：指纹图fing_128.指纹图fing_128.img显微医学图像cell_128.img〔1〕对原图像进行直方图均衡化处理，同屏显示处理前后图像及其直方图，比拟异同，并答复为什么数字图像均衡化后其直方图并非完全均匀分布。〔2〕对原图像参加点噪声，用4-邻域平均法平滑加噪声图像〔图像四周边界不处理，下同〕，同屏显示原图像、加噪声图像和处理后的图像。①不加门限；②加门限，〔其中〕四、实验结果及分析1、直方图均衡化处理程序如下：clc;fid=fopen('F:\数字图像\图像\img\fing_128.img','r');f=fread(fid,[128,128],'uchar');subplot(2,1,1);imshow(f,[0,255]);q=zeros(1,256);forx=1:128fory=1:128q(f(x,y)+1)=q(f(x,y)+1)+1;endends=q./(128*128);X=0:255;subplot(2,1,2);bar(X,s');figure;t=zeros(1,256);t(1)=s(1);fori=2:256t(i)=t(i-1)+s(i);endsubplot(2,1,1);bar(X,t');t0=floor(255*t+0.5);subplot(2,1,2);bar(X,t0');figure;t1=zeros(1,256);fori=1:256t1(t0(i)+1)=s(i)+t1(t0(i)+1);endsubplot(2,1,1);bar(X,t1');f1=zeros(128,128)forx=1:128fory=1:128f1(x,y)=t0(f(x,y)+1);endendsubplot(2,1,2);imshow(f1,[0,255]);运行结果：〔1〕、指纹均衡化处理〔2〕、细胞的均衡化处理2、对原图像参加点噪声程序代码不加门限：clc;fid=fopen('F:\数字图像\图像\img\fing_128.img','r');f=fread(fid,[128,128],'uchar');subplot(2,2,1);imshow(f,[0,255]);forx=1:128fory=1:128ifx==yf(x,y)=255;elseifx+y==100f(x,y)=0;endendendsubplot(2,2,2);imshow(f,[0,255]);f0=f;forx=2:127fory=2:127f0(x,y)=(f((x-1),y)+f((x+1),y)+f(x,(y-1))+f(x,(y+1)))./4;endendsubplot(2,2,3);imshow(f0,[0,255]);t=fft2(f);T=t(1,1)/128;f1=f;forx=2:127fory=2:127h=(f((x-1),y)+f((x+1),y)+f(x,(y-1))+f(x,(y+1)))./4;ifabs(f(x,y)-h)>Tf1(x,y)=h;elsef1(x,y)=f(x,y);endendendsubplot(2,2,4);imshow(f1,[0,255]);运行结果：加门限：clc;fid=fopen('F:\数字图像\图像\img\fing_128.img','r');fg=fread(fid,[128,128],'uchar');subplot(2,2,1);imshow(fg,[0,255]);a=randn(128,128);f=a.*20+fg;subplot(2,2,2);imshow(f,[0,255]);f0=f;forx=2:127fory=2:127f0(x,y)=(f((x-1),y)+f((x+1),y)+f(x,(y-1))+f(x,(y+1)))./4;endendsubplot(2,2,3);imshow(f0,[0,255]);t=fft2(f);T=t(1,1)/128;f1=f;forx=2:127fory=2:127h=(f((x-1),y)+f((x+1),y)+f(x,(y-1))+f(x,(y+1)))./4;ifabs(f(x,y)-h)>Tf1(x,y)=h;elsef1(x,y)=f(x,y);endendendsubplot(2,2,4);imshow(f1,[0,255]);运行结果：细胞参加噪声及去噪3、结果分析〔1〕直方图均衡化处理后