adapt1.m

来自「Adaptive Filter for image processing, a 」· M 代码 · 共 40 行

clc;clear all;close all;
I=imread('C:\Davood\Arshad_2\Image Process\Folder1\images_chapter\images_chapter_05\Fig5.07(a).jpg');
J = imnoise(I,'gaussian');
imshow(I);title('original');
figure;imshow(J);title('niose image')
%**************** var noise**************
J1=imcrop(J,[100 248 100 40]);figure;imshow(J1);
[n m]=size(J1);
V=[];
mnoise=sum(J1(:))/(m*n);
for i=1:n
    for j=1:m
     V(i,j)=(J1(i,j)-mnoise)^2;
    end
end
vnoise=sum(V(:))/(m*n);
%*****************************************
d1=7;
d=(d1-1)/2;
[n m]=size(J);
for i=(d+1):(n-d);
    for j=(d+1):(m-d);
   locmat=J((i-d):(i+d),(j-d):(j+d));
   [n1 m1]=size(locmat);
   ml=sum(locmat(:))/(m1*n1);
   V1=[];
   for k=1:n1
    for l=1:m1
     V1(k,l)=(locmat(k,l)-ml)^2;
    end
   end
vl=sum(V1(:))/(m1*n1);
est(i,j)=J(i,j)-(vnoise/vl)*(J(i,j)-ml);
    end
end
figure;imshow(est);title('adaptive')