欧美性猛交XXXX免费看蜜桃,成人网18免费韩国,亚洲国产成人精品区综合,欧美日韩一区二区三区高清不卡,亚洲综合一区二区精品久久

雖然matlab里面有這些函數，但是攀立民老師要求自己編寫(xiě)，計算機視覺(jué)上有這個(gè)實(shí)驗，到網(wǎng)上找了半天才零散的找到一些碎片，還是整理以后發(fā)上來(lái)吧！

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列的所有元素
for i=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)開(kāi)始的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;
for i=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)開(kāi)始的n行n列元素,即模板(n×n的)
        e=c(1,:);      %是c矩陣的第一行
        for u=2:n
            e=[e,c(u,:)];     %將c矩陣變?yōu)橐粋€(gè)行矩陣   
        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);%計算圖象中心
for i=1:n
    for j=1:n
      b(i,j) =exp(-((i-n1)^2+(j-n1)^2)/(4*k))/(4*pi*k);
    end
end
%生成高斯序列b。
Img1=conv2(Img,b,'same'); %用生成的高斯序列卷積運算，進(jìn)行高斯濾波
d=uint8(Img1);

 

 

%此為程序主文件，包含主要功能單元，以及對子函數進(jìn)行調用
try
%實(shí)驗步驟一：彩色、灰度變換
h=imread('photo.jpg'); %讀入彩色圖片
c=rgb2gray(h); %把彩色圖片轉化成灰度圖片，256級
figure,imshow(c),title('原始圖象'); %顯示原始圖象
g=imnoise(c,'gaussian',0.1,0.002);  %加入高斯噪聲
figure,imshow(g),title('加入高斯噪聲之后的圖象');  %顯示加入高斯噪聲之后的圖象

%實(shí)驗步驟二：用系統預定義濾波器進(jìn)行均值濾波
n=input('請輸入均值濾波器模板大小\n');
A=fspecial('average',n);  %生成系統預定義的3X3濾波器
Y=filter2(A,g)/255;           %用生成的濾波器進(jìn)行濾波,并歸一化
figure,imshow(Y),title('用系統函數進(jìn)行均值濾波后的結果'); %顯示濾波后的圖象

%實(shí)驗步驟三:用自己的編寫(xiě)的函數進(jìn)行均值濾波
Y2=avefilt(g,n);     %調用自編函數進(jìn)行均值濾波，n為模板大小
figure,imshow(Y2),title('用自己的編寫(xiě)的函數進(jìn)行均值濾波之后的結果'); %顯示濾波后的圖象

%實(shí)驗步驟四:用Matlab系統函數進(jìn)行中值濾波
n2=input('請輸入中值濾波的模板的大小\n');
Y3=medfilt2(g,[n2 n2]);   %調用系統函數進(jìn)行中值濾波，n2為模板大小
figure,imshow(Y3),title('用Matlab系統函數進(jìn)行中值濾波之后的結果');  %顯示濾波后的圖象

%實(shí)驗步驟五:用自己的編寫(xiě)的函數進(jìn)行中值濾波
Y4=midfilt(g,n2);      %調用自己編寫(xiě)的函數進(jìn)行中值濾波，
figure,imshow(Y4),title('用自己編寫(xiě)的函數進(jìn)行中值濾波之后的結果');

%實(shí)驗步驟六:用matlab系統函數進(jìn)行高斯濾波
n3=input('請輸入高斯濾波器的均值\n');
k=input('請輸入高斯濾波器的方差\n');
A2=fspecial('gaussian',k,n3);      %生成高斯序列
Y5=filter2(A2,g)/255;              %用生成的高斯序列進(jìn)行濾波
figure,imshow(Y5),title('用Matlab函數進(jìn)行高斯濾波之后的結果');    %顯示濾波后的圖象

%實(shí)驗步驟七:用自己編寫(xiě)的函數進(jìn)行高斯濾波
Y6=gaussfilt(n3,k,g);  %調用自己編寫(xiě)的函數進(jìn)行高斯濾波，n3為均值，k為方差
figure,imshow(Y6),title('用自編函數進(jìn)行高斯濾波之后的結果');      %顯示濾波后的圖象

catch           %捕獲異常
    disp(lasterr);     %如果程序有異常，輸出
end