encryption.m

数字图象加密程序（编写的数字图象加密程序）

clear all;

Image=double(imread('e:\work\lena.jpg'));

figure(1)
subplot(2,2,1)
imshow(Image,[]);%显示原始图像
title('原始图像');
Image=uint8(Image);
figure(1)
subplot(2,2,2)
imhist(Image);
title('原始图像的直方图');

Ih=size(Image,1);%图像的长度
Iw=size(Image,2);%图像的宽度



key=double(input('请输入加密密钥:','s'));%输入任意位密钥，为安全起见，密钥长度不得少于四位

len=length(key);

while(mod(len,4)~=0)%密钥位数不是4的倍数则程序自动取固定值补足末几位
    len=len+1;
    key(len)=15;%程序自动补足默认值为15
end

N=len/4;%四个密钥为一组，共N组
for(i=1:Ih)
    for(j=1:Iw)
        Imagebak(i,j)=Image(i,j);
    end
end

for(n=1:N)
    select=mod(key((n-1)*4+1),4);%四种参数取法，每种取法有两个独立参数
    switch(select)
    case 0
        a=1;
		b=mod(key((n-1)*4+2),Ih);   
		c=mod(key((n-1)*4+3),Ih);   
		d=mod((b*c+1),Ih);
    case 1
		a=mod(key((n-1)*4+2),Ih);
		b=1;
		c=mod(a*mod(key((n-1)*4+3),Ih)-1,Ih);
		d=mod(key((n-1)*4+3),Ih);
   	case 2
        a=mod(key((n-1)*4+2),Ih);
		b=mod(a*mod(key((n-1)*4+3),Ih)-1,Ih);
		c=1;
		d=mod(key((n-1)*4+3),Ih);
    case 3
		a=mod(key((n-1)*4+2),Ih)*mod(key((n-1)*4+3),Ih);
		b=mod(key((n-1)*4+2),Ih);
		c=mod(key((n-1)*4+3),Ih);
		d=1;
    end%这四种取法使猫映射为保面积映射，即映射矩阵行列式为1
    
    iterative=mod(key((n-1)*4+4),17);%每组迭代次数
    
    if(iterative<10)
        iterative=iterative+10;%为防止迭代次数过多或过少，取模17后的值，但若小于10则要加上10
    end	
    
    for(i=1:Ih)
        for(j=1:Iw)
            absc=i;
            ordi=j;
            for(k=1:iterative)
                absc1=mod((a*absc+b*ordi),Ih)+1;
				ordi1=mod((c*absc+d*ordi),Iw)+1;
				absc=absc1;
				ordi=ordi1;
% 				imagebak[i][j]=(image[i][j]+G*abs(absc-ordi))%256;
            end
            Enimage(absc1,ordi1)=Imagebak(i,j);
        end
    end
    for(i=1:Ih)
        for(j=1:Iw)
            Imagebak(i,j)=Enimage(i,j);
        end
    end
end


figure(1)
subplot(2,2,3)
imshow(Enimage,[]);
title('加密图像');
figure(1)
subplot(2,2,4)
imhist(Enimage);
title('加密图像的直方图');

Enimage=uint8(Enimage);

imwrite(Enimage,'Entest.bmp','bmp');