《(完整word版)(绝对可运行)加密图像中的可逆数据隐藏算法matlab代码(基于张新鹏教授的算法).doc》由会员分享，可在线阅读，更多相关《(完整word版)(绝对可运行)加密图像中的可逆数据隐藏算法matlab代码(基于张新鹏教授的算法).doc（5页珍藏版）》请在金锄头文库上搜索。
1、(完整word版)(绝对可运行)加密图像中的可逆数据隐藏算法matlab代码(基于张新鹏教授的算法)%This is the program of Reversible Data Hiding in Encrypted Domain clc;clear;=Read image=original_p=rgb2gray(imread(LENA。tif)）;m,n=size(original_p）；%=Image encryption=random_bits=rand(m，n*8）=0.5; 伪随机序列for i=1：m for j=1:n s=0； for k=0:7 b（k+1）=mod（fix
2、（double（original_p（i，j）)/(2k），2）； eb(k+1）=xor(b(k+1)，random_bits（i,8j+k-7）; s=s+eb（k+1）*（2k)； end encrypted_image(i，j)=s； endendfigure（1）;subplot(1,2,1);imshow(original_p）；hold on；title(Original image,fontsize，16)；subplot(1,2,2）；encrypted_image=uint8(encrypted_image）；imshow（encrypted_image）；title（En
3、crypted image,fontsize,16)；%=Data embedding=size=8； %分块大小N=6;%置乱周期number1=3;number2=Nnumber1；lim_row=fix(m/size)； %图像的分块数lim_col=fix（n/size)； %图像的分块数bitts=3;watermessage=rand(lim_row，lim_col)0.5; %水印信息 %=pseudo_randomly=for i=1:lim_row for j=1:lim_col block_imagei,j=encrypted_image（size*isize+1):siz
4、ei，(size*j-size+1)：size*j）； 分块 for k=1:size for t=1:size index=arnold(k,t，size，number1)；arnold is permutateion function； arno_block_imagei,j（index(1)+1,index(2）+1）=block_imagei，j（k,t)； end end 数据嵌入 sum0=zeros（size/2，size）； if watermessage（i，j）=0 for k=0：7 bbbk+1=mod（fix（double（arno_block_imagei，j（1：
5、size/2,：）/（2k）),2); if kbitts bbbk+1=bbbk+1; end sum0=sum0+bbbk+1（2k); end data_imagei，j=sum0；double（arno_block_imagei，j(size/2+1：size，:)）； end sum1=zeros（size/2,size); if watermessage(i,j）=1 for k=0:7 bbbk+1=mod（fix（double(arno_block_imagei,j(size/2+1：size,：）)/（2k）),2）； if kbitts bbbk+1=bbbk+1; end
6、 sum1=sum1+bbbk+1（2k）； end data_imagei，j=double(arno_block_imagei,j（1:size/2,:)；sum1； end for k=1:size for t=1：size index=arnold（k，t，size,number2）;%arnold is permutateion function； re_data_imagei，j(index(1）+1,index(2)+1)=data_imagei,j（k,t）； end end 合成嵌入数据后的图像 data_en_image(（sizei-size+1):sizei,(size
7、j-size+1)：sizej)=re_data_imagei,j; endend%=Encrypt image which contains data=for i=1：lim_row*size for j=1：lim_colsize s=0； for k=0：7 b（k+1)=mod(fix(data_en_image(i，j)/(2k），2)； eb（k+1）=xor（b(k+1)，random_bits（i，8j+k-7）; s=s+eb（k+1）(2k）; end de_data_image(i，j)=s; endendfigure（2）；hold on；subplot（1，2,1)；
8、imshow（original_p);title(Original image，fontsize,16)；subplot（1，2，2）;imshow（uint8（de_data_image）);title（Decrypted image contains data,fontsize,16）；直接解密图像的峰值信噪比sumsss=0；for i=1：lim_rowsize for j=1：lim_col*size sumsss=（de_data_image(i，j)double(original_p（i，j）)）2+sumsss； endendmse=sumsss/(（lim_rowsize)(
9、lim_row*size)；psnr=10log10(2552/mse)；=Extract data and recover image=for i=1：lim_row for j=1：lim_col block_de_imagei，j=de_data_image（sizeisize+1）:size*i,(sizejsize+1):sizej)； 分块 for k=1:size for t=1：size index=arnold（k，t，size,number1）;%arnold is permutateion function; ar_block_de_imagei,j（index（1）+1
10、,index（2)+1)=block_de_imagei，j（k，t）; end end sum0=zeros（size/2,size); sum1=zeros(size/2，size）； for k=0：7 bbb0k+1=mod(fix（ar_block_de_imagei，j(1：size/2,：）/（2k），2); bbb1k+1=mod(fix(ar_block_de_imagei，j（size/2+1：size，：）/（2k），2）； if kbitts bbb0k+1=bbb0k+1; bbb1k+1=bbb1k+1; end sum0=sum0+bbb0k+1*（2k）; su
11、m1=sum1+bbb1k+1*（2k）; end ar_H0_imagei，j=sum0；ar_block_de_imagei，j（size/2+1:size,:）； ar_H1_imagei,j=ar_block_de_imagei,j(1：size/2，：);sum1； for k=1：size for t=1:size index=arnold(k,t,size,number2);%arnold is permutateion function； H0_imagei,j(index（1）+1,index（2)+1)=ar_H0_imagei,j(k,t）; H1_imagei,j(index(1）+1，index(2)+1)=ar_H1_imagei,j(k，t）; end end f0=0； f1=0； for u=2:size1 for v=2：size1
