均值濾波與matlab編程

MatLab自編的均值濾波、中值濾波、高斯濾波圖像處理函數。

%自編的均值濾波函數。x是需要濾波的圖像,n是模板大小(即n×n)

function d=avefilt(x,n)

a(1:n,1:n)=1; %a即n×n模板,元素全是1

p=size(x); %輸入圖像是p×q的,且p>n,q>n

x1=double(x);

x2=x1;

%A(a:b,c:d)表示A矩陣的第a到b行,第c到d列的所有元素

fori=1:p(1)-n+1

for j=1:p(2)-n+1

c=x1(i:i+(n-1),j:j+(n-1)).*a; %取出x1中從(i,j)開始的n行n列元素與模板相乘

s=sum(sum(c)); %求c矩陣(即模板)中各元素之和

x2(i+(n-1)/2,j+(n-1)/2)=s/(n*n); %將模板各元素的均值賦給模板中心位置的元素end

end

%未被賦值的元素取原值

d=uint8(x2);

%自編的中值濾波函數。x是需要濾波的圖像,n是模板大小(即n×n)

function d=midfilt(x,n)

p=size(x); %輸入圖像是p×q的,且p>n,q>n

x1=double(x);

x2=x1;

fori=1:p(1)-n+1

for j=1:p(2)-n+1

c=x1(i:i+(n-1),j:j+(n-1)); %取出x1中從(i,j)開始的n行n列元素,即模板(n×n的)

e=c(1,:); %是c矩陣的第一行

for u=2:n

e=[e,c(u,:)]; %將c矩陣變為一個行矩陣

end

mm=median(e); %mm是中值

x2(i+(n-1)/2,j+(n-1)/2)=mm; %將模板各元素的中值賦給模板中心位置的元素end

end

%未被賦值的元素取原值

d=uint8(x2);

%自編的高斯濾波函數，S是需要濾波的圖象，n是均值，k是方差

function d=gaussfilt(k,n,s)

Img = double(s);

n1=floor((n+1)/2);%計算圖象中心

fori=1:n

for j=1:n

b(i,j) =exp(-((i-n1)^2+(j-n1)^2)/(4*k))/(4*pi*k);

總結

以上是生活随笔為你收集整理的matlab均值编程,均值滤波与matlab编程的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。