📄 volkan2imge.m
字号:
clc; clear all; close all;
base_image= imread('projlab1.jpg');
base_image1=rgb2gray(base_image);
base_image2=double(base_image1);
% base_image2=base_image2(5:115,5:171);
base_image3 = imread('projlab2.jpg');
base_image4=rgb2gray(base_image3);
base_image5=double(base_image4);
% base_image5=base_image5(5:115,5:171);
base_image6 = imread('proj1lab3.jpg');
base_image7=rgb2gray(base_image6);
base_image8=double(base_image7);
% base_image8=base_image8(5:115,5:171);
base_image9 = imread('proj1lab4.jpg');
base_image10=rgb2gray(base_image9);
base_image11=double(base_image10);
% base_image11=base_image11(5:115,5:171);
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+10)
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:120
for j=1:176
if f5(i,j)==f
m=i; n=j;
m=x-m; n=y-n;
end
end
end
for i=1:120
for j=1:176
if f6(i,j)==g
a=i; b=j; b=y-b;
end
end
end
for i=1:120
for j=1:176
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/(102*141); mu2=toplam2/(102*141); mua=mu1-mu2;
mu3=toplam4/(101*104); mu4=toplam3/(101*104); mub=mu4-mu3;
mu5=toplam5/(100*86); mu6=toplam6/(100*86); muc=mu6-mu5;
imge(60:179,1:y)=base_image2(:,:);
imge(60-m-1:176-m-1,y-n+3:(2*y)-n+1)=base_image5(1:117,2:y)+mua;
imge(60-m+a-1:179-1-m+a,2*(y+1)-n:(3*y)+2-n-b)=base_image8(:,b:y)-mub;
imge(60-m+a-3-c:179-3-m+a-c,3*(y+1)-n-b:4*y+3-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))
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -