數(shù)字圖像處理-實(shí)驗(yàn)

數(shù)字圖像處理——實(shí)驗(yàn)實(shí)驗(yàn)一、計(jì)算圖象統(tǒng)計(jì)參數(shù)在matlab軟件編寫(xiě)程序：讀取圖像（文件名為‘cameraman.tif’）;最大值最小值均值直方圖（要求不能調(diào)用imhist函數(shù)，只能用該函數(shù)來(lái)驗(yàn)證自編直方圖程序塊的正確性）。實(shí)驗(yàn)二灰度和直方圖變換圖像數(shù)據(jù)讀出計(jì)算并分析圖像直方圖（imhist）利用直接灰度變換法對(duì)圖像進(jìn)行灰度變換（imadjust）直方圖均衡化處理及顯示（histeq）空間濾波給定圖像添加噪聲（高斯噪聲和椒鹽噪聲）（imnoise）使用平均濾波器、中值濾波器對(duì)不同強(qiáng)度的高斯噪聲和椒鹽噪聲，進(jìn)行濾波處理（fspecial,filter2,medfilt2）實(shí)驗(yàn)二——參考答案f=imread('medicine_pic.jpg');g=imhist(f,256);%顯示其直方圖g1=imadjust(f,[01],[10]);%灰度轉(zhuǎn)換，實(shí)現(xiàn)明暗轉(zhuǎn)換(負(fù)片圖像)figure,imshow(g1)%將0.5到0.75的灰度級(jí)擴(kuò)展到范圍[01]g2=imadjust(f,[0.50.75],[01]);figure,imshow(g2)g=imread('point.jpg');h=log(1+double(g));%對(duì)輸入圖像對(duì)數(shù)映射變換h=mat2gray(h);%將矩陣h轉(zhuǎn)換為灰度圖片h=im2uint8(h);%將灰度圖轉(zhuǎn)換為8位圖figure,imshow(h)J=histeq(f);%對(duì)原圖像進(jìn)行直方圖均衡化處理Figure,imshow(J);%顯示均衡化圖像%對(duì)直方圖均衡化后圖像進(jìn)行屏幕控制；作一幅子圖，并排兩幅圖的第1幅figure;subplot(1,2,1);imhist(I,64);%將原圖像直方圖顯示為64級(jí)灰度title(‘原圖像直方圖’);%給原圖像直方圖加標(biāo)題名subplot(1,2,2);%作第2幅子圖imhist(J,64);%將均衡化后圖像的直方圖顯示為64級(jí)灰度title(‘均衡變換后的直方圖’);%給均衡化后圖像直方圖加標(biāo)題名實(shí)驗(yàn)二——參考答案I=imread('electric.tif');J=imnoise(I,'gauss',0.02);%添加高斯噪聲J=imnoise(I,'salt&pepper',0.02);(注意空格)

%添加椒鹽噪聲ave1=fspecial('average',3);%產(chǎn)生3×3的均值模版ave2=fspecial('average',5);%產(chǎn)生5×5的均值模版K=filter2(ave1,J)/255;%均值濾波3×3L=filter2(ave2,J)/255;%均值濾波5×5M=medfilt2(J,[33]);%中值濾波3×3模板N=medfilt2(J,[44]);%中值濾波4×4模板imshow(I);figure,imshow(J);figure,imshow(K);figure,imshow(L);figure,imshow(M);figure,imshow(N);實(shí)驗(yàn)三1、利用MATLAB實(shí)現(xiàn)數(shù)字圖像的傅里葉變換讀入并顯示圖test1.tif，作該圖的二維FFT變換F，將其直流分量移到頻譜中心F1（fftshift()函數(shù)），計(jì)算其實(shí)部RR、虛部II，用兩種方法計(jì)算幅值A(chǔ)1=abs(F1)和A2=sqrt(RR.^2+II.^2)，分別顯示A1和A2，并加以比較。test1.tifi=imread(‘test1.tif');figure,imshow(i);title(‘test1.tif');f=fft2(i);f1=fftshift(f);rr=real(f1);ii=imag(f1);a1=abs(f1);a2=sqrt(rr.^2+ii^2);aa1=log(1+a1);aa2=log(1+a2);figure,imshow(aa1,[]);title('a1abs');figure,imshow(aa2,[]);title('a2sqrt');實(shí)驗(yàn)三——參考答案實(shí)驗(yàn)三2、近似沖擊函數(shù)二維傅里葉變換

A=zeros(99,99);

A(49:51,49:51)=1;作A的二維傅里葉變換B，將B直流分量移到頻譜中心B1，分別用函數(shù)imshow和mesh顯示A和B1模的對(duì)數(shù)（log(1+abs(B1))）A=zeros(99,99);A(49:51,49:51)=1;B=fft2(A);B1=fftshift(B);figure,imshow(A);title('沖擊函數(shù)imsow');figure,mesh(A);title('沖擊函數(shù)mesh');figure,imshow(log(1+abs(B1)));title('沖擊函數(shù)頻域模mesh');figure,mesh(log(1+abs(B1)))title('沖擊函數(shù)頻域模mesh');實(shí)驗(yàn)三——參考答案實(shí)驗(yàn)三3、空間濾波與頻域?yàn)V波將圖test2.tif分別進(jìn)行空間與頻域?yàn)V波?？臻g濾波：用fspecial產(chǎn)生9*9，標(biāo)準(zhǔn)差為2的高斯濾波器w，用函數(shù)imfilter對(duì)f進(jìn)行空間濾波，得到濾波后的圖像fi1。頻域?yàn)V波:將上述高斯濾波器w用函數(shù)freqz2獲得其256*256頻域形式W，用W在頻域?yàn)V波圖像f（注意W的直流分量在頻譜中心），得到Ff，求其傅里葉逆變換（ifft2），得到濾波后的圖像fi2。用函數(shù)imshow顯示w，用函數(shù)mesh顯示W(wǎng)，用mshow顯示f，fi1，F(xiàn)f的模的對(duì)數(shù)，fi2。比較fi1與fi2。test2.tifg=imread(test2.tif');G=fft2(g);w=fspecial('gaussian',9,2);fi1=imfilter(g,w);W=freqz2(w,256,256);W1=ifftshift(W);Ff=G.*W1;fi2=real(ifft2(Ff));

figure,imshow(w);title('