کد:
function [A,output]=getA(imput,m,n)
mnsize=3;
len1=zeros(1,(mnsize*2+1)^2);
len2=len1;
len3=len2;
output=zeros(m,n);
minmed=zeros(1,3);
imput1=zeros(m+2*mnsize,n+2*mnsize,3);
imput1(1:mnsize,mnsize+1:mnsize+n,=imput(1:mnsize,:,;
imput1(m+mnsize+1:end,mnsize+1:n+mnsize,=imput(end-mnsize+1:end,:,;
imput1(mnsize+1:m+mnsize,mnsize+1:n+mnsize,=imput(:,:,;
imput1(:,1:mnsize,=imput1(:,mnsize+1:mnsize*2,;
imput1(:,n+mnsize+1:end,=imput1(:,n+1:n+mnsize,;
for i=1+mnsize:m+mnsize
for j=1+mnsize:n+mnsize
len1=imput1(i-mnsize:i+mnsize,j-mnsize:j+mnsize,1);
len2=imput1(i-mnsize:i+mnsize,j-mnsize:j+mnsize,2);
len3=imput1(i-mnsize:i+mnsize,j-mnsize:j+mnsize,3);
a1=min(len1);
a2=min(len2);
a3=min(len3);
minmed=[a1,a2,a3];
output(i-mnsize,j-mnsize)=min(minmed);
end
end
[h,indx]=max(output);
[s,indy]=max(h);
A=imput(indx(indy),indy,;
figure,imshow(output);