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

數(shù)字圖像處理—實(shí)驗(yàn)三一．實(shí)驗(yàn)內(nèi)容：⑴圖像的銳化：使用Sobel，Laplacian，unsharpmasking，high-boostfiltering算子分別對(duì)圖像進(jìn)行運(yùn)算，觀察并體會(huì)運(yùn)算結(jié)果。⑵綜合練習(xí)：對(duì)需要進(jìn)行處理的圖像分析，正確運(yùn)用所學(xué)的知識(shí)，采用正確的步驟，對(duì)圖像進(jìn)行各類處理，以得到令人滿意的圖像效果。二．實(shí)驗(yàn)?zāi)康模簩W(xué)會(huì)用Matlab中的以下函數(shù)對(duì)輸入圖像按實(shí)驗(yàn)內(nèi)容進(jìn)行運(yùn)算；感受各種不同的圖像處理方法對(duì)最終圖像效果的影響。imfilter；fspecial；imadjust；imadd；immultiply.三．實(shí)驗(yàn)步驟：仔細(xì)閱讀Matlab幫助文件中有關(guān)以上函數(shù)的使用說明，能充分理解其使用方法并能運(yùn)用它們完成實(shí)驗(yàn)內(nèi)容。將Fig3.41(c).jpg圖像文件讀入Matlab，使用imfilter函數(shù)分別采用不同的算子對(duì)其作銳化運(yùn)算，顯示運(yùn)算前后的圖像。算子的輸入可采用直接輸入法。其中Sobel，Laplacian，unsharpmasking也可用fspecial函數(shù)產(chǎn)生。各類算子如下： High-boostfilteringUnsharpmasking對(duì)角線LaplacianSobel垂直梯度High-boostfilteringUnsharpmasking對(duì)角線LaplacianSobel垂直梯度將Fig3.46(a).jpg圖像文件讀入Matlab，按照以下步驟對(duì)其進(jìn)行處理：用帶對(duì)角線的Laplacian對(duì)其處理，以增強(qiáng)邊緣。用imadd函數(shù)疊加原始圖像?？梢钥闯鲈肼曉鰪?qiáng)了，應(yīng)想法降低。獲取Sobel圖像并用imfilter對(duì)其進(jìn)行5×5鄰域平均，以減少噪聲。用immultiply函數(shù)獲取Sobel，Laplacian增強(qiáng)圖像，噪聲得以減少。最后用imadjust函數(shù)做power-law灰度變換，以增強(qiáng)dynamicrange。四．實(shí)驗(yàn)報(bào)告要求： 用imshow函數(shù)顯示經(jīng)每一步處理后的圖像，打印后提交。本次實(shí)驗(yàn)后要求認(rèn)真寫出三次實(shí)驗(yàn)的心得體會(huì)及提出改良實(shí)驗(yàn)的一些建議。實(shí)驗(yàn)報(bào)告上的其他內(nèi)容，按常規(guī)實(shí)驗(yàn)報(bào)告要求辦。%2i2=imread('Fig3.41(c).jpg');i2=im2double(i2);subplot(3,2,1);imshow(i2);title('orgimg');h1=[-1-2-1000121];h2=[1111-81111];h3=[1111-91111];A=-1;h4=[1111-A-81111];subplot(3,2,2);g1=imfilter(i2,h1);g1=imadd(i2,g1);imshow(g1);title('g1');subplot(3,2,3);g2=imfilter(i2,h2);g2=imadd(i2,g2);imshow(g2);title('g2');subplot(3,2,4);g3=imfilter(i2,h3);g3=imadd(i2,g3);imshow(g3);title('g3');subplot(3,2,5);g4=imfilter(i2,h4);g4=imadd(i2,g4);imshow(g4);title('g4');%5i5=imread('Fig3.46(a).jpg');i5=im2double(i5);subplot(3,2,6);imshow(i5);title('5orgimg');%51figure;subplot(3,2,1);g51=imfilter(i5,h2);imshow(g51);title('51');%52subplot(3,2,2);i52=imadd(i5,g51);imshow(i52);title('52');subplot(3,2,3);h52=fspecial('average',[33]);%生成一3×3鄰域平均窗函數(shù)g52=filter2(h52,i51);%求鄰域平均imshow(g52);title('g52');%53subplot(3,2,4);H=fspecial('sobel');SOBEL=imfilter(i5,H);i53=imadd(i5,SOBEL);imshow(i53);title('53');subplot(3,2,5);H=fspecial('average',[55]);%生成一5×5鄰域平均窗函數(shù)g53=filter2(H,i53);%求鄰域平均imshow(g53);title('g53');%54subplot(3,2,6);i54l=immultiply(i52,i52);imshow(i54);title('54l');figure;subplot(3,2,1);i54s=immultiply(i53,i53);imshow(i54);title('54s');%55subplot(3,2,2);i55l=imadjust(i52,[],[],0.2);imshow(i55l);title('55l');subplot(3,2,3);i55=imadjust(i54l,[],[],0.1);imshow(i55);title('g55l');subplot(3,2,4);i55s=imadjust(i53,[],[],0.1);imshow(