阿爾法均值濾波

1、圖像復(fù)原的目的是從觀測到的退化圖像重建原始圖像，它是圖像處理、 模式識別、機器視覺等的基礎(chǔ)，在天文學(xué)、遙感成像、醫(yī)療圖像等領(lǐng)域獲 得了重要應(yīng)用。代碼：I=imread(g.png);m,n=size(I);P=imnoise(I,salt & pepper,0.04);Q=imnoise(I,gaussian,0,0.02);R=imnoise(P,gaussian,0,0.02);%加椒鹽噪聲的三種情況A=I;%d=0d=0;for i=2:m-1for j=2:n-1img_turn = P(i-1:i+1,j-1:j+1);img_turn = sort(img_turn(:);img_

2、decrese = img_turn(floor(d/2)+1:9-floor(d/2);A(i,j) = sum(img_decrese(:)/(9-d);endendB=I;%d=2d=2;for i=2:m-1for j=2:n-1img_turn = P(i-1:i+1,j-1:j+1);img_turn = sort(img_turn(:);img_decrese = img_turn(floor(d/2)+1:9-floor(d/2);B(i,j) = sum(img_decrese(:)/(9-d);endendC=I;%d=4d=4;for i=2:m-1for j=2:n-1

3、img_turn = P(i-1:i+1,j-1:j+1);img_turn = sort(img_turn(:);img_decrese = img_turn(floor(d/2)+1:9-floor(d/2);C(i,j) = sum(img_decrese(:)/(9-d);end end%加高斯噪聲的三種情況D=I;%d=0d=0;for i=2:m-1for j=2:n-1img_turn = Q(i-1:i+1,j-1:j+1);img_turn = sort(img_turn(:);img_decrese = img_turn(floor(d/2)+1:9-floor(d/2);

4、D(i,j) = sum(img_decrese(:)/(9-d);endendE=I;%d=2d=2;for i=2:m-1for j=2:n-1img_turn = Q(i-1:i+1,j-1:j+1);img_turn = sort(img_turn(:);img_decrese = img_turn(floor(d/2)+1:9-floor(d/2);E(i,j) = sum(img_decrese(:)/(9-d);endendF=I;%d=4d=4;for i=2:m-1for j=2:n-1img_turn = Q(i-1:i+1,j-1:j+1);img_turn = sort

5、(img_turn(:);img_decrese = img_turn(floor(d/2)+1:9-floor(d/2);F(i,j) = sum(img_decrese(:)/(9-d);endend%加混合噪聲的三種情況G=I;%d=0d=0;for i=2:m-1for j=2:n-1img_turn = R(i-1:i+1,j-1:j+1);img_turn = sort(img_turn(:);img_decrese = img_turn(floor(d/2)+1:9-floor(d/2);G(i,j) = sum(img_decrese(:)/(9-d);endendH=I;%d

6、=2d=2;for i=2:m-1for j=2:n-1 img_turn = R(i-1:i+1,j-1:j+1);img_turn = sort(img_turn(:);img_decrese = img_turn(floor(d/2)+1:9-floor(d/2);E(i,j) = sum(img_decrese(:)/(9-d);endendJ=I;%d=4d=4;for i=2:m-1for j=2:n-1img_turn = R(i-1:i+1,j-1:j+1);img_turn = sort(img_turn(:);img_decrese = img_turn(floor(d/2

7、)+1:9-floor(d/2);J(i,j) = sum(img_decrese(:)/(9-d);endendfigure;imshow(I);title(原圖);figure;imshow(P);title(椒鹽噪聲圖);figure;imshow(Q);title(高斯噪聲圖);figure;imshow(R);title(混合噪聲圖);figure;imshow(uint8(A);title( d=0 時去椒鹽噪聲圖);figure;imshow(uint8(B);title( d=2 時去椒鹽噪聲圖);figure;imshow(uint8(C);title( d=4 時去椒鹽噪聲圖);figure;imshow(uint8(D);title( d=0 時去高斯噪聲圖);figure;imshow(uint8(E);title( d=2 時去高斯噪聲圖);figure;imshow(uint8(F);title( d=4 時去高斯噪聲圖);figure;imshow(uint8(G);title( d=0 時去混合噪聲圖);figure;imshow(uint8(H);title( d=2 時去混合噪聲圖);figure;imshow(uint8(J);title( d=4 時去混合噪聲圖);效果:原圖椒鹽噪聲圖混合噪聲圖d=0時去椒鹽噪