imagehide.m

直接4bit替换法的图像伪装技术以及对第四bit考察的matlab代码

%参数说明：
%cover是载体图像的地址
%message是秘密图像的地址
%stegocover是隐藏后图像的地址
%level是作为载体的机体层，R为1，G为2，B为3
%data是隐藏后图像的矩阵
function data=imagehide(cover,message,stegocover,permission,level)
%提取图像信息并分层
cover=imread(cover,permission);
data=cover;
msg=imread(message,permission);
[row,col]=size(cover);
cover1=cover(:,:,level);
%置载体图像R层的低四位为0
for i=1:row
    for j=1:col/3
        cover1(i,j)=bitand(cover1(i,j),240);
    end
end

%置秘密图像的低四位为0
takemsg4=bitand(msg,240);
%将秘密图像的高四位右移四位
shiftmsg4=bitshift(takemsg4,-4);
%图像隐藏
for i=1:row
    for j=1:col/3
        cover1(i,j)=bitor(cover1(i,j),shiftmsg4(i,j));
    end
end

%写回并保存
data(:,:,level)=cover1;
imwrite(data,stegocover,permission);

%提取秘密图像信息，检测隐藏效果
data=imread(stegocover,permission);
[row,col]=size(data);
A=data(:,:,level);
for i=1:row
    for j=1:col/3
        A(i,j)=bitand(A(i,j),15);
    end
end
A=bitshift(A,4);
%显示结果
subplot(221);imshow(cover);title('两幅图像中所选的载体图像');
subplot(222);imshow(message);title('两幅图像中所选的秘密图像');
subplot(223);imshow(data);title('将秘密图像隐藏在载体图像后所得到的的图像');
subplot(224);imshow(A);title('从处理后的载体图像中提取出的秘密图像');
% figure;
% t=double(cover(:,:,level))-double(data(:,:,level));
% imshow(t,[]);





%参数说明：
%cover是载体图像的地址
%message是秘密图像的地址
%stegocover是隐藏后图像的地址
%level是作为载体的机体层，R为1，G为2，B为3
%count是替换表
%thershold是阈值
function count=fourthbitcmp(cover,message,permission,level,threshold)
%提取图像信息并分层
cover=imread(cover,permission);
data=cover;
msg=imread(message,permission);
cover1=cover(:,:,level);
%对cover和msg的第四位进行处理
tempc=cover1;
tempm=msg;
tempc=bitand(tempc,8);
tempm=bitand(tempm,8);
temp=bitxor(tempm,tempc);
[row,col]=size(temp);

%记录图像每个分块的n值
k1=0;
k2=0;
a=row*col/64;
count=zeros([1 a]);
for i=1:a
    for m=1:8
        for n=1:8
            if temp(8*k1+m,8*k2+n)==0
                count(1,i)=count(1,i)+1;
            end
        end
    end
    k2=k2+1;
    if k2*8==col
        k2=0;
        k1=k1+1;
    end
end

%计算每块的u值与阈值进行比较
count=count/64;
for i=1:a
    if count(i)>=threshold
        count(i)=1;  %可以替换
    elseif count(i)<1-threshold
        count(i)=-1;  %取反
    else
        count(i)=0;%不能处理
    end
end

s=0;
for i=1:a
    if count(i)==0
        s=s+1;        
    end    
end
disp(s);