實驗二空間域圖像增強一、實驗目的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〕直方圖均衡化處理后