课设lsb算法.docxVIP

信息隐藏课设LSB算法实验步骤隐藏提取及测试（1）隐藏算法，保存在hide_lsb.m文件中。function o = hide_lsb(block, data, I)%function o = hide_lsb(block, data, I)%block:the minimal ceil to hide%data:the information%I:source imagesi = size(I);lend = length(data);%将图像划分为M*N个小块N = floor(si(2) / block(2));M = min(floor(si(1) / block(1)), ceil(lend / N));o = I;fori = 0 : M-1%计算每小块垂直方向起止位置rst = i * block(1) + 1;red = (i + 1) * block(1);for j = 0 : N-1%计算每小块隐藏的秘密信息的序号idx = i * N + j + 1;ifidx lendbreak;end;%取每小块隐藏的秘密信息bit = data(idx);%计算每小块水平方向起止位置cst = j * block(2) + 1;ced = (j + 1) * block(2);%将每小块最低位平面替换为秘密信息 o(rst:red, cst:ced) = bitset(o(rst:red, cst:ced), 1, bit);end;end;（2）提取算法，保存在dh_lsb.m文件中。function out = dh_lsb(block, I)%function out = dh_lsb(block, I)%block:the minimal ceil to dehide%I:stegoed imagesi = size(I);%将图像划分为M*N个小块N = floor(si(2) / block(2));M = floor(si(1) / block(1));out = [];%计算比特1判决阈值：即每小块半数以上元素隐藏的是比特1时，判决该小块嵌入的信息为1thr = ceil((block(1) * block(2) + 1) / 2);idx = 0;fori = 0 : M-1%计算每小块垂直方向起止位置rst = i * block(1) + 1;red = (i + 1) * block(1);for j = 0 : N-1%计算每小块将要数据的秘密信息的序号idx = i * N + j + 1;%计算每小块水平方向起止位置cst = j * block(2) + 1;ced = (j + 1) * block(2);%提取小块最低位平面，统计1比特个数，判决输出秘密信息tmp = sum(sum(bitget(I(rst:red, cst:ced), 1)));if(tmp = thr)out(idx) = 1;elseout(idx) = 0;end;end;end;（3）测试脚本，保存在test.m文件中。fid = 1;len = 10;%随机生成要隐藏的秘密信息d = randsrc(1, len, [0 1]);block = [3, 3];[fn, pn] = uigetfile({*.bmp, bmp file(*.bmp);}, 选择载体);s = imread(strcat(pn, fn));ss = size(s);if(length(ss) = 3)I = rgb2gray(s);else I = s;end;si = size(I);sN = floor(si(1) / block(1)) * floor(si(2) / block(2));tN = length(d);%如果载体图像尺寸不足以隐藏秘密信息，则在垂直方向上复制填充图像ifsN tNmultiple = ceil(tN / sN);tmp = [];fori = 1:multipletmp = [tmp; I];end; I = tmp;end;%调用隐藏算法，把隐蔽载体写至硬盘。stegoed = hide_lsb(block, d, I);imwrite(stegoed, hide.bmp, bmp);[fn, pn] = uigetfile({*.bmp, bmp file(*.bmp);}, 选择隐蔽载体);y = imread(strcat(pn, fn));sy = size(y);if(length(sy) = 3) I = rgb2gray(y);else I = y;end; %调用提取算法，获得秘密信息。out = dh_lsb(block,