projestandartsapma.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
sigma1=std2(base_image2); sigma2=std2(base_image5);
sigma4=std2(base_image11); sigma3=std2(base_image8);
alpha3=sigma4/sigma3;
alpha1=sigma2/sigma1; 
for i=1:111
    for j=1:167
        beta1=mean2(base_image5(i,j)-alpha1*mean2(base_image2));
        beta3=mean2(base_image11(i,j)-alpha3*mean2(base_image8));
    end
end
base_image2=(base_image2-beta1)/alpha1;
base_image5=(base_image5-beta1)/alpha1;  
base_image11=(base_image11-beta3)/alpha3;

sigma3=std2(base_image8); sigma2=std2(base_image5);
alpha2=sigma3/sigma2;
for i=1:111
    for j=1:167
        beta2=mean2(base_image8(i,j)-alpha2*mean2(base_image5));
    end
end
base_image8=(base_image8-beta2)/alpha2;






imge(60+2:170+2,1:y)=base_image2(:,:);
imge(60-m:170-m,y:(2*y)-n)=base_image5(:,n:y);
imge(60-m+a-1:170-m-1+a,2*y+1-n:(3*y)+1-n-b)=base_image8(:,b:y);
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);
% imge(180-m:180,2*(y+1)-n:(3*y)+2-n-b)=base_image8(,y-b:y)
imshow(uint8(imge))