volkankesilmis.m

If you have not registered, Please [regist first].You should upload at least five sourcecodes/docume

clc; clear all; close all;

base_image= imread('projlab1.jpg');
base_image1=rgb2gray(base_image);
base_image1=base_image1(5:115,5:171);
base_image2=double(base_image1);

base_image3 = imread('projlab2.jpg');
base_image4=rgb2gray(base_image3);
base_image4=base_image4(5:115,5:171);
base_image5=double(base_image4);

base_image6 = imread('proj1lab3.jpg');
base_image7=rgb2gray(base_image6);
base_image7=base_image7(5:115,5:171);
base_image8=double(base_image7);

base_image9 = imread('proj1lab4.jpg');
base_image10=rgb2gray(base_image9);
base_image10=base_image10(5:115,5:171);
base_image11=double(base_image10);

toplam1=0; toplam2=0; toplam3=0; toplam4=0; toplam5=0; toplam6=0;
[x y]=size(base_image5); 
for i=1:(x*2)
    for j=1:(y*2+20)
        imge(i,j)=255;
    end
end 
f1=fft2(base_image1); f2=fft2(base_image4);
f3=fft2(base_image7); f4=fft2(base_image10);

f5=abs(ifft2((f1.*conj(f2))./abs(f1.*conj(f2))));
f6=abs(ifft2((f2.*conj(f3))./abs(f2.*conj(f3))));
f7=abs(ifft2((f3.*conj(f4))./abs(f3.*conj(f4))));
f=max(max(f5)); g=max(max(f6)); h=max(max(f7));
for i=1:111
    for j=1:167
        if f5(i,j)==f
            m=i; n=j;
           m=x-m; n=y-n; 
        end
    end
end
for i=1:111
    for j=1:167
        if f6(i,j)==g
            a=i; b=j; b=y-b;
        end
    end
end
for i=1:111
    for j=1:167
        if f7(i,j)==h
            c=i; d=j; c=x-c; d=y-d;
        end
    end
end
for i=1:x-m
    for j=n:y
        toplam1=base_image2(i,j)+toplam1;
    end
end
for i=m+1:x
    for j=1:y-n+1
        toplam2=base_image5(i,j)+toplam2;
    end
end
for i=m-a+1:x
    for j=b:y
        toplam3=base_image5(i,j)+toplam3;
    end
end
for i=1:x-(m-a+1)
    for j=1:y-b+1
        toplam4=base_image8(i,j)+toplam4;
    end
end
for i=1:x-(m-a+3+c)
    for j=d:y
        toplam5=base_image8(i,j)+toplam5;
    end
end
for i=m-a+3+c:x
    for j=1:y-d+1
        toplam6=base_image11(i,j)+toplam6;
    end
end
mu1=toplam1/((x-m-1)*(y-n)); mu2=toplam2/((x-m-1)*(y-n)); mua=mu1-mu2;
mu3=toplam4/((x-(m-a+2))*(y-b)); mu4=toplam3/((x-(m-a+2))*(y-b)); mub=mu4-mu3+mua;
mu5=toplam5/((x-(m-a+4+c))*(y-d));  mu6=toplam6/((x-(m-a+4+c))*(y-d)); muc=mu6-mu5;

imge(60+2:170+2,1:y)=base_image2(:,:)-mua;
imge(60-m:170-m,y:(2*y)-n)=base_image5(:,n:y)-mua;
imge(60-m+a-1:170-m-1+a,2*y+1-n:(3*y)+1-n-b)=base_image8(:,b:y)+mub;
imge(60-m-3+a-c:170-m-3+a-c,3*y+2-n-b:4*y+2-n-b-d)=base_image11(:,d:y)-muc;
% imge(180-m:180,2*(y+1)-n:(3*y)+2-n-b)=base_image8(,y-b:y)
imshow(uint8(imge))