clear
A=imread('meili.BMP')
subplot(331)%1行1列第一个位置
imshow(A)
A=double(A)/255;
[M,N]=size(A)
C=imread('suolong.BMP')
subplot(332)
imshow(C)
C=double(C)/255;
[M,N]=size(A)
B(M,N)=0;
for i=1:M
for j=1:N
B(i,j)=C(i,j)+A(i,j);
end
end
subplot(333)
imshow(B);
M=600,N=600,K=3;
X(M,N,K)=0;
X(:,:,1)=1;
X(:,:,2)=1;
X(:,:,3)=1;
X1(M,N,K)=0;
X1(:,:,1)=0.2;
X1(:,:,2)=0.5;
X1(:,:,3)=0.6;
subplot(121)
imshow(X);
subplot(122)
imshow(X1);
Frames=5;
Y(M,N,K,Frames)=0;
x0=300;y0=300;
L=40;Sigama=200;V=10;
Sigama2=Sigama*Sigama*2;
for t=1:Frames
for i=1:M
for j=1:N
rr=(i-x0)*(i-x0)+(j-y0)*(j-y0);
r=sqrt(rr);
alpha=exp(-rr/Sigama2)*((cos(2*pi*(r-V*(t-1))/L)+1)/2);
for k=1:K
Y(i,j,k,t)=alpha*X(i,j,k)+(1-alpha)*X1(i,j,k);
end
end
end
end
clc
for t=1:Frames
pause(0.1)
imshow(Y(:,:,:,t))
end
A=imread('shi.BMP')
subplot(334)
imshow(A)
A=double(A)/255;
[M,N,K]=size(A)
C=imread('shamo.BMP')
subplot(335)
imshow(C)
C=double(C)/255;
B(M,N,K)=0;
for i=1:M
for j=1:N
for k=1:3
B(i,j,k)=C(i,j,k)+A(i,j,k);
end
end
end
subplot(336)
imshow(B);
for i=1:M
for j=1:N
for k=1:3
B(i,j,k)=1-A(i,j,k);
end
end
end
figure
imshow(B)
subplot(337)
晕,程序都有了,合成也要问?删减,改参数就完了!