% Second project of Image Processing 
% Athor : Dariush Semnani

%Read image and show it
[filename,pathname]=uigetfile('e:\*.*','select a file:');
I=imread([pathname filename]);

Fmat=im2double(rgb2gray(I));
[M N]=size(Fmat);

h2=im2col(Fmat,[5 5],'distinct');
h3=mean(h2);
[k p]=size(h3);
for i=1:25
    for j=1:p
        h4(i,j)=h3(j);
    end
end
Fmat2=col2im(h4,[5 5],[M N],'distinct');

thick=5; %(mm)
r1=max(Fmat2);maximum=max(r1);
r2=min(Fmat2);minimum=min(r2);
Fmat3=Fmat2-minimum;
Fmat4=Fmat3/(maximum-minimum);
H=Fmat4*thick*0.25;
r1=max(H);maximum=max(r1);
r2=min(H);minimum=min(r2);
K=(H+thick-maximum);
K(M/2,N/2)=0;
mesh(K(round(M/3):round(2*M/3),round(N/3):round(2*N/3)));
Refplate(1:M,1:N)=mean2(H);