matlab实现Jsteg隐写分析

clc; clear all; gray=imread('lena.jpg'); temp=double(gray); [m,n]=size(temp); T=dctmtx(8);%产生二维DCT变换矩阵 Q_table=[16 11 10 16 24 40 51 61 12 12 14 19 26 58 60 55 14 13 16 24 40 57 69 56 14 17 22 29 51 87 80 62 18 22 37 56 68 109 103 77 24 35 55 64 81 104 113 92 49 64 78 87 103 121 120 101 72 92 95 98 112 100 103 99]; B=blkproc(temp,[8,8],'P1*x*P2',T,T');%blkproc分块处理函数 B=blkproc(B,[8,8],'round(x./P1)',Q_table);%round四舍五入函数，得到DCT系数 C=zeros(1,m*n-m*n/64);%创建一个1行m*n-m*n/64列的矩阵,存储量化后的DCT系数 K=1; for i=1:m for j=1:n if (mod(i,8)~=1)||(mod(j,8)~=1)%去掉DC(直流)值，每一块左上角的第一个数是DC系数。 C(1,K)=B(i,j);%提取AC值到C矩阵中 K=K+1; end end end %计算AC系数不等于0且不等于1的个数 v=0; for i=1:m*n-m*n/64%m*n-m*n/64是指用所有像素个数减去DC值，就剩下AC值 if(C(1,i)~=0)&&(C(1,i)~=1)%系数不等于0且不等于1 v=v+1; end end rate=0.8; A=floor(v*rate);%向下取整，取一部分 M=rand(1,A)>0.5;%随机产生矩阵 %从C中取一不大于A的矩阵 j=0; for i=1:m*n-m*n/64 if(C(1,i)~=0)&&(C(1,i)~=1) j=j+1; end if j==A break; end end%在长度为A的隐写率下，C应该取到第i位 S=zeros(1,i);%构造C的一个子集，i=j+AC位0或1的个数 %隐写过程 l=1; for r=1:i S(1,r)=C(1,r); end for r=1:i if (S(1,r)~=0)&&(S(1,r)~=1) if (S(1,r)>0)%正数隐写 S(1,r)=S(1,r)-mod(S(1,r),2)+uint8(M(1,l));%隐写过程，M是随机产生的秘密信息 elseif (S(1,r)<0)%负数隐写 if (mod(S(1,r),2)==1)%负奇数项 if (M(1,l)==0) S(1,r)=S(1,r)-1; else S(1,r)=S(1,r); end else%负偶数项 if (M(1,l)==1) S(1,r)=S(1,r)+1;%-4变-3，-2变-1 else S(1,r)=S(1,r); end end end l=l+1; if l>A break;%隐写结束 end end end for r=1:i C(1,r)=S(1,r);%隐写之后的C矩阵 end p = JSteg(C,m,n,rate);