二值信息隐藏（分块和游程编码实现）

使用分块进行信息隐藏，因为在对角线上的分块上进行的隐藏，所以

可以明显看到在对角线上有一条线，

200*200的二值图像

 

512*512的二值图像

 

（二）使用游程编码，书上的代码逻辑上有错误，还有一些函数错误，根本不能运行出结果

自己修改了得到以下结果

 

200*200的隐藏160位，可以看到微小的变化

 

512*512的隐藏160位，基本看不到变化

 

二值对角线分块隐藏代码

clc;
clear;

msgfid = fopen('hidden.txt','r'); % 隐藏信息
[msg, count] = fread(msgfid);
fclose(msgfid);

msg = dec2bin(msg,8);   %将读入的每个字节转换成二进制
msg = str2num(msg(:));
watermarklen = count*8;    % 需要隐藏的数目

[fn,pn] = uigetfile({'*.bmp','bmpfile(*.bmp)';},'chose file');  
filename = strcat(pn,fn);
I = imread(filename);  % 将载体分成 count*count 个小块，在对角线隐藏信息
[row,col] =size(I);

block(1) = floor(row/watermarklen);  % block的row
block(2) = floor(col/watermarklen);  % block的col
pixelcount = block(1)*block(2); % 一块的像素数

% 阈值,1的个数必须大于等于此值才能说明隐藏的是1
% 0 的个数必须大于等于此值才能说明隐藏的是 0
if mod(pixelcount, 2) == 0 
    threshold = pixelcount/2 + 1;
else
    threshold = ceil(pixelcount/2);
end
    
carrier = I;  % 修改的操作都在carrier上

zeroCount = zeros(1,watermarklen) ;    % 用来记录每小块上0的个数
oneCount = zeros(1,watermarklen);    % 用来记录每小块上1的个数

% 统计对角线上的小块的0和1
for n=1:watermarklen
    %使用求和来算出每个小块有多少个·1
    oneCount(1,n) = sum(sum(I(block(1)*(n-1)+1:block(1)*n, block(2)*(n-1)+1:block(2)*n)));
    zeroCount(1,n) = pixelcount - oneCount(1,n);
end
%隐藏信息
for n=1:watermarklen
    
    if msg(n,1) == 1; % 嵌入1
        
        if oneCount(1,n) < threshold; % 但是1小于阈值
            modifycount = threshold - oneCount(1,n);  % 需要修改的数目
            k =1;
            for i =block(1)*(n-1)+1:block(1)*n
                for j =block(2)*(n-1)+1:block(2)*n
                    if(carrier(i,j) == 0 && k<=modifycount)
                        carrier(i,j) = 1;
                        k = k+1;
                    end
                end
            end
        end
    else % 嵌入0
        if zeroCount(1,n) < threshold;  % 0的个数小于阈值时需要改
            modifycount = threshold - zeroCount(1,n);
            k =1;
            for i =block(1)*(n-1)+1:block(1)*n
                for j =block(2)*(n-1)+1:block(2)*n
                     if(carrier(i,j) ==1 && k<=modifycount)
                        carrier(i,j) = 0;
                        k = k+1;
                     end
                end
            end
        end
    end
end

%显示结果
imwrite(carrier,'bwmarked.bmp');  
figure;
imshow(filename);
title('原图');
figure;
imshow('bwmarked.bmp');
title('载体');

提取代码

clc;
clear;

[fn,pn] = uigetfile({'*.bmp','bmpfile(*.bmp)';},'chose file');  
filename = strcat(pn,fn);
I = imread(filename);  % 将载体分成 count*count 个小块，在对角线隐藏信息
watermarklen = 160;  % 隐藏的数目
[row,col] =size(I);

block(1) = floor(row/watermarklen);  % block的row
block(2) = floor(col/watermarklen);  % block的col
pixelcount = block(1)*block(2); % 一块的像素数

% 阈值,1的个数必须大于等于此值才能说明隐藏的是1
% 0 的个数必须大于等于此值才能说明隐藏的是 0
if mod(pixelcount, 2) == 0 
    threshold = pixelcount/2 + 1;
else
    threshold = ceil(pixelcount/2);
end

%提取出的隐藏信息
hideInfo = zeros(1,watermarklen);

for n=1:watermarklen
    
    %使用求和来算出每个小块有多少个·1
    oneCount = sum(sum(I(block(1)*(n-1)+1:block(1)*n, block(2)*(n-1)+1:block(2)*n)));
    if oneCount >= threshold % 1的数目大于阈值，则隐藏的信息为1
        hideInfo(1, n) = 1;
    else                      % 0的数目大于阈值，则隐藏的信息为0
        hideInfo(1, n) = 0;
    end
end


hideInfo = reshape(hideInfo,length(hideInfo)/8,8);  %将二进制字符串分成n行，8列
bin = num2str(hideInfo); % 二进制数值转变成二进制字符串
dec = bin2dec(bin);%二进制字符串转换成十进制数组
recoverdata = native2unicode(dec);% 本机编码转换成unicode
fprintf(1, '恢复出的信息： %s
',recoverdata);

游程编码隐藏代码

clc;
clear;

threshold = 10;
[fn,pn] = uigetfile({'*.bmp','bmpfile(*.bmp)';},'chose file');  
filename = strcat(pn,fn);
oi = imread(filename);
si = size(oi); % 保存读入图像的行列,以便恢复图像
oi = oi(:);  % 变成一列, 按列变
[len, col] = size(oi); %len= row,1
carrier = oi;   % 以后数据的修改都在carrier上

j=1;
i=1;
%统计游程长度
while(i <= len)
    last = oi(i); % 依次取出每个数
    count = 1;
    while( i+1 <= len && oi(i+1)== last) % 下一个与上一个相同
        i = i + 1;
        count = count + 1;
    end
    RLE(j) = count;
    j = j + 1;
    i = i + 1;
end

msgfid = fopen('hidden.txt', 'r');
[msg, msgcount] = fread(msgfid);
fclose(msgfid);

msg = dec2bin(msg,8);   %将读入的每个字节转换成二进制
msg = str2num(msg(:));  % 
msgcount = msgcount*8;    % 需要隐藏的数目

i = 1; % msg 的索引
count = 1; % RLE的索引,只为奇数 1 3 5...
sum = 0;
while count < length(RLE)-1
    addSum = RLE(count) + RLE(count+1);   %相邻两个数相加的和，从1开始
    
    if count >= 3  
        sum = sum + RLE(count - 1) + RLE(count - 2); % 当前count的在数组中的下标-1
    end
    
    if(addSum >= threshold)  % 只有当相邻游程和大于给定数的时候,才隐藏信息
        
        if RLE(count) >= RLE(count+1)
            modifyPoint = sum + RLE(count); % 修改两个游程的临界值,长的那边的
        else
            modifyPoint = sum + RLE(count) + 1;
        end
        
        if(msg(i) == 0)  % 嵌入的信息为0，但是游程为奇数，需要修改
            if(mod(RLE(count), 2) == 1) %RLE(count) == '奇数'   
                carrier(modifyPoint) = mod(carrier(modifyPoint)+1 ,2); % 修改点的值取反 mod(carrier(modifyPoint)+1 ,2)
            end
        else  % 嵌入的信息为1，但是游程为偶数，需要修改
            if mod(RLE(count), 2) == 0 %RLE(count) == '偶数' 
                carrier(modifyPoint) = mod(carrier(modifyPoint) + 1 ,2); % 修改点的值取反 mod(carrier(modifyPoint)+1 ,2)
            end
        end
        
        i = i + 1; % 隐藏一个位后， msg索引加1
    end
    
    count = count + 2; %不管有没有隐藏信息，count都要增加
    
    if i == msgcount  % 所有的信息都已隐藏就退出循环
        break;
    end
end

if i < msgcount
    error('不能隐藏全部信息！');
end

imwrite(reshape(carrier, si(1),si(2)), 'hide.bmp');
figure;
subplot(121);
imshow(reshape(oi, si(1),si(2)));
title('原图');
subplot(122);
imshow('hide.bmp');
title('载体');
    

提取代码

clc;
clear;


%也可以增大隐藏信息的间隔 count 1 5 9 13 以后实现
threshold = 3;
[fn,pn] = uigetfile({'*.bmp','bmpfile(*.bmp)';},'chose file');  
filename = strcat(pn,fn);
oi = imread(filename);
oi = oi(:);  % 变成一列, 按列变
[len, col] = size(oi); %len= row,1

j=1;
i=1;
%统计游程长度
while(i <= len)
    last = oi(i); % 依次取出每个数
    count = 1;
    while( i+1 <= len && oi(i+1)== last) % 下一个与上一个相同
        i = i + 1;
        count = count + 1;
    end
    RLE(j) = count;
    j = j + 1;
    i = i + 1;
end

i = 1; 
msgcount = 160;
msg = zeros(msgcount,1);
count = 1; % RLE的索引,只为奇数 1 3 5...
while i <= msgcount
   
    addSum = RLE(count) + RLE(count+1);   %相邻两个数相加的和，从1开始   
    if(addSum >= threshold)  % 只有当相邻游程和大于给定数的时候,才隐藏信息
        if mod(RLE(count),2) == 0
            msg(i) = 0;
        else
            msg(i) = 1;
        end
        i = i + 1;
    end
    count = count + 2;     
end

hideInfo = reshape(msg,length(msg)/8,8);  %将二进制字符串分成n行，8列
bin = num2str(hideInfo); % 二进制数值转变成二进制字符串
dec = bin2dec(bin);%二进制字符串转换成十进制数组
recoverdata = native2unicode(dec);% 本机编码转换成unicode
fprintf(1, '恢复出的信息： %s
',recoverdata);