five的圖象處理大作業(yè)

憨憨選了個圖象處理課，結(jié)果老師上課瘋狂吹逼，具體方法只有很少一點，然后布置了一個大作業(yè)，大概就是讀取圖像然后加個高斯白噪聲，然后再想辦法把噪聲濾除，以及做一個銳化。這里直接上一個代碼

I=imread('lena512.bmp'); imfinfo('lena512.bmp'); figure,imshow(I); y = uint8(wgn(512,512,60)); figure,imshow(y); I1=I+y; for i=1:1:512for j=1:1:512if I1(i,j)>256I1(i,j)=256;endend end for i=1:1:512for j=1:1:512if I1(i,j)<0I1(i,j)=0;endend end figure,imshow(I1); k1=filter2(fspecial('average',3),I1)/255; %進行3*3模板平滑濾波 k2=filter2(fspecial('average',5),I1)/255; %進行5*5模板平滑濾波 k3=filter2(fspecial('average',7),I1)/255; %進行7*7模板平滑濾波 k4=filter2(fspecial('average',9),I1)/255; %進行9*9模板平滑濾波 subplot(141),imshow(k1);title('3*3 模板平滑濾波'); subplot(142),imshow(k2);title('5*5 模板平滑濾波'); subplot(143),imshow(k3);title('7*7 模板平滑濾波'); subplot(144),imshow(k4);title('9*9 模板平滑濾波');p1 =fspecial('sobel');%銳化 figure(); add1 =imfilter(I1,p1); out1=add1+I1; imshow(out1); title('sobel銳化結(jié)果'); % % figure,imhist(I); % % I2=imnoise(I,'gaussian',0,0.1); % % figure,imshow(I2);

第一步

讀取圖像，matlab有自帶的圖像讀取函數(shù)imread，讀取出來如果你是512512的黑白圖像就會用個512512的數(shù)組來儲存，如果是rgb就是5125123的數(shù)組，這下我門就完成了第一步，那個imfinfo是可以看圖像具體數(shù)據(jù)的函數(shù)，imread應(yīng)該可以讀很多格式的圖像

第二步

這里本來可以直接用imnoise函數(shù)加噪聲，但是由于老師要求單獨的高斯噪聲圖像這里使用的是wgn函數(shù)，wgn函數(shù)具體用法：
wgn（行，列，噪聲功率）
當(dāng)然也可以用rand（），之后還要乘一個常數(shù)再加到原圖像上去，兩者之間有個換算關(guān)系：wgn(m,n,p)=sqrt*(10^(p/10))*rand(m,n)
然后得到一個噪聲圖像和加了噪聲的圖像

線性平滑濾波

簡單地通過matlab的函數(shù)實現(xiàn)。。。
filter2(fspecial(‘a(chǎn)verage’,3),I1)/255，這里的3可以替換成其他奇數(shù)，一般越大越模糊這里給一組對比圖：

第四步銳化

用fspecial函數(shù)選擇sobel算子來做銳化，簡單的一批，關(guān)于銳化的原理，說白了提取邊緣分量，怎么提？用一個算子一個一個滑動過去，至于算子怎么設(shè)計？下次再開個專題來詳細(xì)說說有哪些算子，是怎么設(shè)計的，現(xiàn)在就給個成果圖：
總之先寫到這里，有弟弟和我說C++的opencv做這個會很快，我覺得matlab已經(jīng)夠快了，等我學(xué)了C++再出一期openCV做圖像處理的合集。

分割線-------------
一年過去了我還是沒有學(xué)會C++。?？佳悬h是人下人。

總結(jié)

以上是生活随笔為你收集整理的matlab图像处理Lena大作业的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。