blob: ecf3273ba301053a76e0666155164e731bcec5d0 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
|
//The input to the function are the corrupted image ‘a’ and the dimension
function [Out_Imag] = Func_medianall(a,N)
a=double(a);
[m n]=size(a);
Out_Imag=a;
if(modulo(N,2)==1)
Start=(N+1)/2;
End=Start;
else
Start=N/2;
End=Start+1;
end
if(modulo(N,2)==1)
limit1=(N-1)/2;
limit2=limit1;
else
limit1=(N/2)-1;
limit2=limit1+1;
end
for i=Start:(m-End+1),
for j=Start:(n-End+1),
I=1;
for k=-limit1:limit2,
for l=-limit1:limit2,
mat(I)=a(i+k,j+l);
I=I+1;
end
end
mat=gsort(mat); //Sort the elements to find the median
if(modulo(N,2)==1)
Out_Imag(i,j)=(mat(((N^2)+1)/2));
else
Out_Imag(i,j)=(mat((N^2)/2)+mat(((N^2)/2)+1))/2;
end
end
end
|