chapter4图像信息伪装技术.docVIP

chapter 04 图像信息伪装技术 listing 4.1 幻方的求取 %文件名：magicsquares.m %程序员：祝素雯 %编写时间：2011.11.26 %函数功能：本函数将完成n阶二维幻方的求取，要求n为奇数 %输入格式举例：result = magicsquares(5) %参数说明： %n为阶数 %result为求得的二维幻方矩阵 function result = magicsquares(n) if mod(n, 2) == 0 error(n要求为奇数); end result = zeros(n); %根据P128算法，先将1放在第1行的中间一列 j = floor(n / 2) + 1; %中间1列 result(1, j) = 1; %从2到n^2每个的数依次存放在前一个书的上一行、后一列的位置上... %但下列特殊情况除外 i = n + 1; %便于以后从第n行开始考虑起 for k = 2 : n * n i = i - 1; j = j + 1; %特殊情况：如果前一个数存放在第1行第n列的位置上， %则将当前数字存放在前一个数的下面（即第2行第n列） if i 1 j n i = i + 2; j = j - 1; else %特殊情况：如果前一个数存放在第1行，则将当前数字存放在第n行 if i 1 i = n; end %特殊情况：如果前一个数存放在第n列，则将当前数字存放在第1行 if j n j = 1; end end if result(i,j) == 0 result(i,j) = k; else %特殊情况：如果按上述规则确定的位置上已经有数存放了(result(i,j)!=0) %则直接将当前数字存放在前一个数的下面（即行数+1，列数不变） i = i + 2; j = j - 1; result(i,j) = k; end end 测试结果 listing 4.2 幻方置乱主函数 %文件名：magicreplace.m %程序员：祝素雯 %编写时间：2011.11.26 %函数功能：本函数将完成对输入的图像信号按幻方置换策略进行置乱 %输入格式举例：result = magicreplace(secretimage, 1, 1983) %参数说明： %matrix为输入图像矩阵 %key为控制密钥 %cord为1表示置乱变换，为0表示复原变换 %result为置乱后的结果 function result = magicreplace(matrix, eord, key) %----------要求输入为图像矩阵而非图像文件---------- % matrix = imread(image); % %matrix = double(matrix) / 255; %----------这里以11阶幻方为例，可以利用magicsquares函数---------- %-----------------生成高阶幻方，带入进行置乱操作----------------- margin = 11; %分析原图像尺寸并补遗，使其能被分成若干完整的11*11子块 [m,n] = size(matrix); rowadd = margin - mod(m, margin); coladd = margin - mod(n, margin); if rowadd == margin rowadd = 0; end if coladd == margin coladd = 0; end input = uint8(zeros([m + rowadd n + coladd])); input(1:m, 1:n) = matrix(1:m, 1:n); %密钥生成随机数 rand(seed, key); control = randint(1,1,[1 margin * margin]); %生成11阶标准幻方 %magic = magicsquares(margin); %这里取了一个已生成的11阶幻方为例 magic = [ 38 52 66 69 83 97 100 114 7 21 24 61 7