rationallogist.m

来自「有理bezier曲线应用于图像加密的源代码（Matlab）」· M 代码 · 共 54 行

clc;
clear all;

Ximage=double(imread('e:\tu\house.jpg'));
key0=0.1;
key1=4;
[M,N]=size(Ximage);
figure(1)
Ximage=uint8(Ximage);
imshow(Ximage,[]);
Ximage=double(Ximage);

for(i=1:M)
    for(j=1:N)
       key0=key1*key0*(1-key0);
       a(i,j)=key0;
    end
end

key3=0.2;
key2=4;
for(i=1:M)
    for(j=1:N)
        key3=key2*key3*(1-key3);
        b(i,j)=key3;
    end
end

key4=0.3;
key2=4;
for(i=1:M)
    for(j=1:N)
        key4=key2*key4*(1-key4);
        c(i,j)=key4;
    end
end

t=0.4;
w0=0.2;
w1=0.5;
w2=0.3;
w=(1-t)^2*w0+2*t*(1-t)*w1+t^2*w2;
for(i=1:M)
   for(i=j:N)  
       P(i,j)=(1-t)^2*a(i,j)*w0+2*t*(1-t)*b(i,j)*w1+t^2*c(i,j)*w2;
       d(i,j)=P(i,j)/w;
    end
end
% plot(d);
figure(2)
d=uint8(d);
imshow(d,[]);
d=double(d);