0
0

[HUST]多媒体安全实验：图像隐写JSTEG，F4,F5的实现

引言：在变换域进行隐写相比于直接隐写在图像的像素值信息中的方法，具有更好的鲁棒性，也更难从肉眼上被察觉。

本次实验使用MATLAB进行编写，由于在此之前我基本没有使用过MATLAB进行变成，因而可能会有很多不雅观的地方，希望各位能多多包容。而且因为本人水平有限，所以还望大家包容指正。

1.JSTEG

  1.     相对于直接嵌入灰度而言，选择使用DCT隐写的一大原因有是因为DCT系数表示的是各种图像变化的程度，因而，从理论上而言，一张自然图像的DCT图像应该呈现一种山峰形的变化趋势。以0对称，左右逐渐下降的形状。

JSTEG的隐写方法是直接更改DCT绝对值大于1的系数的最低位，将信息通过这种方式嵌入图像之中。在解析时跳过绝对值小于等于1的DCT系数，便可以直接读出隐写的信息。

这种方法显然足够简单，而缺点是JSTEG由于直接更改了最低位，因而会产生值对现象，打个比方，从理论上来说，1出现的频次应该会比2高上很多，但是如果选择了JSTEG方法的话，就会发现，1和2的出现频次会随着嵌入的信息变多，而越来越接近，因而如果分析者发现了这种值对现象，便不难发现了这些图像被处理过。

原图像DCT系数直方图：

而嵌入之后的图像可以看到明显的值对现象：

不难看出之间的差距。

F4嵌入：

  1.     对于F4而言，他将用正奇数和负偶数代表秘密消息1、负奇数和正偶数代表秘密消息0。一旦不一致则将绝对值-1.
  3.     这种做法的优点是相较于JSTEG而言，不会发生值对现象（毕竟他相当于全方位的平移）。

F5嵌入:

  1.     F5是F4的一种直接改进，对于F5而言，他可以有很多有很多形式，具体而言，他可以用2^(n+1)-1位来表示n个比特信息，对于n=1而言，第一个比特信息是b1=a1⊕a2,b2=a2⊕a3,那么，无论何种情况，你只需要在3个DCT系数中修改一个，即可以表示两个字符，这样子的好处是大大增加了DCT嵌入效率，坏处是对于同样要嵌入2个比特信息，正常只需要2个DCT系数，而在这里则需要3个，这就是所谓的有得必有失吧。

JSTEG代码实现:

  1.     首先去这个网站下载：http://dde.binghamton.edu/download/nsf5simulator/

将该文件拖入你的工作目录中即可正常运行（也就是当前脚本的文件夹)

  1.     你需要一个原始的图像，并且将它命名为cover.jpg（当然你也可以跟随代码进行改动）
  3.     你修改后的图像是stego.jpg
  5.   代码实现:
  1. %%% setup
  2. COVER = 'cover.jpg'; % cover image (grayscale JPEG image)
  3. STEGO = 'stego.jpg'; % resulting stego image that will be created
  4. SEED=99;
  5. myInformation='Hello,Matlab'
  6. msgBin = de2bi(int8(myInformation),8,'left-msb');
  7. len = size(msgBin,1).*size(msgBin,2);
  8. myInfo = reshape(double(msgBin).',len,1).';
  9. tic;
  10. [nzAC] = nsf5_simulation(COVER,STEGO,myInfo);
  11. T = toc;
  13. fprintf('-----\n');
  14. fprintf('nsF5 simulation finished\n');
  15. fprintf('cover image: %s\n',COVER);
  16. fprintf('stego image: %s\n',STEGO);
  17. fprintf('PRNG seed: %i\n',SEED);
  18. %fprintf('relative payload: %.4f bpac\n',ALPHA);
  19. fprintf('number of nzACs in cover: %i\n',nzAC);
  21. fprintf('elapsed time: %.4f seconds\n',T);
  23. function [AC]=nsf5_simulation(COVER,STEGO,message)
  24. try
  25.     jobj=jpeg_read(COVER);
  26.     DCT=jobj.coef_arrays{1}
  27. catch
  28.     error('ERROR (problem with the cover image)');
  29. end
  31. AC=numel(DCT)-numel(DCT(1:8:end,1:8:end));%计算AC系数数量，每个DCT分块是8*8的
  32. if(length(message)>AC)
  33.     error('ERROR (too long message)');
  34. end
  35. %打乱Changeable的排列顺序
  36. idD=1;
  37. len=length(message);
  38. for id=1:len
  39.     while(abs(DCT(idD))<=1)
  40.         message(id);
  41.         DCT(idD)=0;
  42.         idD=idD+1;
  43.         if(idD>=AC)
  44.             break;
  45.         end
  46.     end
  47.     if(DCT(idD)<0)
  48.         message(id)
  49.         DCT(idD)
  50.     if(message(id)==mod(DCT(idD),2))
  51.         
  52.         if(mod(DCT(idD),2)==1)
  53.         DCT(idD)=DCT(idD)-sign(DCT(idD));
  54.         else
  55.         DCT(idD)=DCT(idD)+sign(DCT(idD));
  56.         end    
  57.     end
  58.     else
  59.                message(id)
  60.         DCT(idD)
  61.     if(message(id)~=mod(DCT(idD),2))
  62.     
  63.         if(mod(DCT(idD),2)==1)
  64.             DCT(idD)=DCT(idD)-sign(DCT(idD));
  65.         else
  66.             DCT(idD)=DCT(idD)+sign(DCT(idD));
  67.         end
  68.     end
  69.   
  70.         
  71.     end
  72.       idD=idD+1;
  73. end
  75. try
  76.     jobj.coef_arrays{1}=DCT;
  77.     jobj.optimize_coding=1;
  78.     jpeg_write(jobj,STEGO);
  79. catch
  80.     error('ERROR (problem with saving the stego image)')
  81. end
  82. end
  83. function res=invH(y)
  84. to_minimize=@(x) (H(x)-y)^2;
  85. res=fminbnd(to_minimize,eps,0.5-eps);
  86. end
  87. function res=H(x)
  88. res=-x*log2(x)-(1-x)*log2(1-x);
  89. end

F4代码实现：

  1. %%% setup
  2. COVER = 'cover.jpg'; % cover image (grayscale JPEG image)
  3. STEGO = 'stego.jpg'; % resulting stego image that will be created
  4. SEED=99;
  5. output=fopen("test.txt","w");
  6. myInformation='SEEYOUINHE'
  7. msgBin = de2bi(int8(myInformation),8,'left-msb');
  8. len = size(msgBin,1).*size(msgBin,2);
  9. myInfo = reshape(double(msgBin).',len,1).';
  10. tic;
  11. [nzAC] = nsf5_simulationF4(COVER,STEGO,myInfo);
  12. T = toc;
  14. fprintf('-----\n');
  15. fprintf('nsF5 simulation finished\n');
  16. fprintf('cover image: %s\n',COVER);
  17. fprintf('stego image: %s\n',STEGO);
  18. fprintf('PRNG seed: %i\n',SEED);
  19. %fprintf('relative payload: %.4f bpac\n',ALPHA);
  20. fprintf('number of nzACs in cover: %i\n',nzAC);
  22. fprintf('elapsed time: %.4f seconds\n',T);
  24. function [AC]=nsf5_simulationF4(COVER,STEGO,message)
  25. output=fopen("test.txt","w");
  26. try
  27.     jobj=jpeg_read(COVER);
  28.     DCT=jobj.coef_arrays{1}
  29. catch
  30.     error('ERROR (problem with the cover image)');
  31. end
  33. AC=numel(DCT)-numel(DCT(1:8:end,1:8:end));%计算AC系数数量，每个DCT分块是8*8的
  34. if(length(message)>AC)
  35.     error('ERROR (too long message)');
  36. end
  37. %打乱Changeable的排列顺序
  38. idD=1;
  39. len=length(message)
  40. for id=1:len
  41.     while(abs(DCT(idD))<=1)
  42.         if(DCT(idD)<0 && message(id)==1)
  43.             DCT(idD)=0
  44.         elseif(DCT(idD)>0 && message(id)==0)
  45.             DCT(idD)=0
  46.         elseif(DCT(idD)~=0)
  47.             message(id)
  48.             DCT(idD)
  49.             break
  50.         else
  51.         end
  52.         
  53.         idD=idD+1;
  54.         if(idD>=AC)
  55.             break;
  56.         end
  57.     end
  58.     if(DCT(idD)<=-1)
  59.     if(message(id)==mod(DCT(idD),2))
  60.         fprintf("%d",id);
  61.         
  62.         message(id)
  64.         sign(DCT(idD));
  65.         DCT(idD)=DCT(idD)+1;
  66.                 DCT(idD)
  67.         idD=idD+1;
  68.         continue
  69.     end
  70.             message(id)
  71.             DCT(idD)
  72.     fprintf(output,"%d\n",DCT(idD));
  73.     
  74.     idD=idD+1;
  75.     continue
  76.     end
  77.     if(DCT(idD)>=1)
  78.     if(message(id)~=mod(DCT(idD),2))
  79.         message(id)
  80.     
  81.         sign(DCT(idD));
  82.         DCT(idD)=DCT(idD)-1;
  83.             DCT(idD)
  84.         idD=idD+1;
  85.         continue;
  86.     end
  87.         message(id)
  88.         DCT(idD)
  89.     fprintf(output,"%d\n",DCT(idD));
  90.     idD=idD+1;
  91.     continue
  92.     end
  93.     
  94. end
  96. try
  97.     jobj.coef_arrays{1}=DCT;
  98.     jobj.optimize_coding=1;
  99.     jpeg_write(jobj,STEGO);
  100. catch
  101.     error('ERROR (problem with saving the stego image)')
  102. end
  103. end
  104. function res=invH(y)
  105. to_minimize=@(x) (H(x)-y)^2;
  106. res=fminbnd(to_minimize,eps,0.5-eps);
  107. end
  108. function res=H(x)
  109. res=-x*log2(x)-(1-x)*log2(1-x);
  110. end

F5代码实现：

  1. %%% setup
  2. COVER = 'cover.jpg'; % cover image (grayscale JPEG image)
  3. STEGO = 'stego.jpg'; % resulting stego image that will be created
  4. SEED=99;
  5. output=fopen("test.txt","w");
  6. myInformation='HELLOWOLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824HELLOWORLDU201911824';
  7. msgBin = de2bi(int8(myInformation),8,'left-msb');
  8. len = size(msgBin,1).*size(msgBin,2);
  9. myInfo = reshape(double(msgBin).',len,1).';
  10. Extend=2-mod(len,2)
  11. for extend=1:Extend
  12.     myInfo(len+extend)=0
  13. end
  14. [nzAC] = nsf5_simulationF5(COVER,STEGO,myInfo);
  15. T = toc;
  17. fprintf('-----\n');
  18. fprintf('nsF5 simulation finished\n');
  19. fprintf('cover image: %s\n',COVER);
  20. fprintf('stego image: %s\n',STEGO);
  21. fprintf('PRNG seed: %i\n',SEED);
  22. %fprintf('relative payload: %.4f bpac\n',ALPHA);
  23. fprintf('number of nzACs in cover: %i\n',nzAC);
  25. fprintf('elapsed time: %.4f seconds\n',T);
  26. function [mba1,mba2,mba3]=my_AND(a1,a2,a3,b1,b2)
  27. if(a1>0)
  28.     m1=mod(a1,2);
  29. else
  30.     ma1=a1*-1+1;
  31.     m1=mod(ma1,2);
  32. end
  33. if(a2>0)
  34.     m2=mod(a2,2);
  35. else
  36.     ma2=a2*-1+1;
  37.     m2=mod(ma2,2);
  38. end
  39. if(a3>0)
  40.     m3=mod(a3,2);
  41. else
  42.     ma3=a3*-1+1;
  43.     m3=mod(ma3,2);
  44. end
  45. mb1=mod(m1+m2,2);
  46. mb2=mod(m2+m3,2);
  47. if(mb1~=b1&&mb2==b2)
  48.     mba2=a2;
  49.     mba3=a3;
  50.     mba1=a1-sign(a1);
  51. elseif(mb1==b1&&mb2~=b2)
  52.     mba3=a3-sign(a3);
  53.     mba1=a1;
  54.     mba2=a2;
  55. elseif(mb1~=b1&&mb2~=b2)
  56.     mba2=a2-sign(a2);
  57.     mba1=a1;
  58.     mba3=a3;
  59. else
  60.     mba1=a1;
  61.     mba2=a2;
  62.     mba3=a3;
  63. end
  64. end
  65. function [AC]=nsf5_simulationF5(COVER,STEGO,message)
  66. output=fopen("test.txt","w");
  67. try
  68.     jobj=jpeg_read(COVER);
  69.     DCT=jobj.coef_arrays{1}
  70. catch
  71.     error('ERROR (problem with the cover image)');
  72. end
  74. AC=numel(DCT)-numel(DCT(1:8:end,1:8:end));%计算AC系数数量，每个DCT分块是8*8的
  75. AC=AC/(3/2);
  76. if(length(message)>AC)
  77.     error('ERROR (too long message)');
  78. end
  79. idD=1;
  80. len=length(message)
  81. for id=1:2:len
  82.     place1=0;
  83.     place2=0;
  84.     place3=0;
  85.     value1=0;
  86.     value2=0;
  87.     value3=0;
  88.     Value1=0;
  89.     Value2=0;
  90.     Value3=0;
  91.     a1=message(id);
  92.     a2=message(id+1);
  93.     while(abs(DCT(idD))<=0)
  94.         DCT(idD)=0;
  95.         idD=idD+1;
  96.         if(idD>=AC)
  97.             break;
  98.         end
  99.     end
  100.     place1=idD;
  101.     value1=DCT(idD);
  102.     idD=idD+1;
  103.     while(abs(DCT(idD))<=0)
  104.         DCT(idD)=0;
  105.         idD=idD+1;
  106.         if(idD>=AC)
  107.             break;
  108.         end
  109.     end
  110.     place2=idD;
  111.     value2=DCT(idD);
  112.     idD=idD+1;
  113.     while(abs(DCT(idD))<=0)
  114.         DCT(idD)=0;
  115.         idD=idD+1;
  116.         if(idD>=AC)
  117.             break;
  118.         end
  119.     end
  120.     
  121.     place3=idD;
  122.     value3=DCT(idD);
  123.     idD=idD+1;
  124.     
  125.     [Value1,Value2,Value3]=my_AND(value1,value2,value3,a1,a2);
  126.     while(Value1==0||Value2==0||Value3==0)
  127.        if(Value1==0&&Value2~=0&&Value3~=0)
  128.            DCT(place1)=0;
  129.            value1=value2;
  130.            place1=place2;
  131.            place2=place3;
  132.            value2=value3;
  133.            
  134.          while(abs(DCT(idD))<=0)
  135.             DCT(idD)=0;
  136.             idD=idD+1;
  137.             if(idD>=AC)
  138.                 break;
  139.             end
  140.          end
  141.             value3=DCT(idD)
  142.             place3=idD
  143.             idD=idD+1;
  144.        [Value1,Value2,Value3]=my_AND(value1,value2,value3,a1,a2); 
  145.        elseif(Value1~=0&&Value2==0&&Value3~=0)
  146.            DCT(place2)=0;
  147.            place2=place3;
  148.            value2=value3;
  149.            
  150.          while(abs(DCT(idD))<=0)
  151.             DCT(idD)=0;
  152.             idD=idD+1;
  153.             if(idD>=AC)
  154.                 break;
  155.             end
  156.          end
  157.             value3=DCT(idD);
  158.             place3=idD;
  159.             idD=idD+1;
  160.        [Value1,Value2,Value3]=my_AND(value1,value2,value3,a1,a2);
  161.        elseif(Value1~=0&&Value2~=0&&Value3==0)
  162.            DCT(place3)=0;
  163.          while(abs(DCT(idD))<=0)
  164.             DCT(idD)=0;
  165.             idD=idD+1;
  166.             if(idD>=AC)
  167.                 break;
  168.             end
  169.          end
  170.             value3=DCT(idD)
  171.             place3=idD
  172.             idD=idD+1;
  173.        [Value1,Value2,Value3]=my_AND(value1,value2,value3,a1,a2); 
  174.        elseif(Value1==0&&Value2==0&&Value3~=0)
  175.            DCT(place1)=0;
  176.            DCT(place2)=0;
  177.            value1=value3
  178.            place1=place3
  179.          while(abs(DCT(idD))<=0)
  180.             DCT(idD)=0;
  181.             idD=idD+1;
  182.             if(idD>=AC)
  183.                 break;
  184.             end
  185.          end
  186.             value2=DCT(idD)
  187.             place2=idD
  188.             idD=idD+1;
  189.            while(abs(DCT(idD))<=0)
  190.             DCT(idD)=0;
  191.             idD=idD+1;
  192.             if(idD>=AC)
  193.                 break;
  194.             end
  195.          end
  196.             value3=DCT(idD);
  197.             place3=idD;
  198.             idD=idD+1;
  199.        [Value1,Value2,Value3]=my_AND(value1,value2,value3,a1,a2); 
  200.       elseif(Value1==0&&Value2~=0&&Value3==0)
  201.            DCT(place1)=0;
  202.            DCT(place3)=0;
  203.            value1=value2
  204.            place1=place2
  205.          while(abs(DCT(idD))<=0)
  206.             DCT(idD)=0;
  207.             idD=idD+1;
  208.             if(idD>=AC)
  209.                 break;
  210.             end
  211.          end
  212.             value2=DCT(idD)
  213.             place2=idD
  214.             
  215.             idD=idD+1;
  216.            while(abs(DCT(idD))<=0)
  217.             DCT(idD)=0;
  218.             idD=idD+1;
  219.             if(idD>=AC)
  220.                 break;
  221.             end
  222.          end
  223.             value3=DCT(idD)
  224.             place3=idD
  225.             idD=idD+1;
  226.        [Value1,Value2,Value3]=my_AND(value1,value2,value3,a1,a2);
  227.        
  228.         elseif(Value1~=0&&Value2==0&&Value3==0)
  229.            DCT(place2)=0;
  230.            DCT(place3)=0;
  231.          while(abs(DCT(idD))<=0)
  232.             DCT(idD)=0;
  233.             idD=idD+1;
  234.             if(idD>=AC)
  235.                 break;
  236.             end
  237.          end
  238.             value2=DCT(idD)
  239.             place2=idD
  240.             
  241.             idD=idD+1;
  242.            while(abs(DCT(idD))<=0)
  243.             DCT(idD)=0;
  244.             idD=idD+1;
  245.             if(idD>=AC)
  246.                 break;
  247.             end
  248.          end
  249.             value3=DCT(idD)
  250.             place3=idD
  251.             
  252.             idD=idD+1;
  253.        [Value1,Value2,Value3]=my_AND(value1,value2,value3,a1,a2); 
  254.        else
  255.            DCT(place1)=0;
  256.            DCT(place2)=0;
  257.            DCT(place3)=0;
  258.          while(abs(DCT(idD))<=0)
  259.             DCT(idD)=0;
  260.             idD=idD+1;
  261.             if(idD>=AC)
  262.                 break;
  263.             end
  264.          end
  265.             value1=DCT(idD)
  266.             place1=idD
  267.             
  268.             idD=idD+1;
  269.            while(abs(DCT(idD))<=0)
  270.             DCT(idD)=0;
  271.             idD=idD+1;
  272.             if(idD>=AC)
  273.                 break;
  274.             end
  275.          end
  276.             value2=DCT(idD)
  277.             place2=idD
  278.             
  279.             idD=idD+1;
  280.             while(abs(DCT(idD))<=0)
  281.             DCT(idD)=0;
  282.             idD=idD+1;
  283.             if(idD>=AC)
  284.                 break;
  285.             end
  286.          end
  287.             value3=DCT(idD)
  288.             place3=idD
  289.             
  290.             idD=idD+1;
  291.        [Value1,Value2,Value3]=my_AND(value1,value2,value3,a1,a2); 
  292.            
  293.        end
  294.        
  295.     end
  296.     temp1=message(id);
  297.     temp2=message(id+1);
  298.     DCT(place1)=Value1;
  299.     DCT(place2)=Value2;
  300.     DCT(place3)=Value3;
  301. end
  303. try
  304.     jobj.coef_arrays{1}=DCT;
  305.     jobj.optimize_coding=1;
  306.     jpeg_write(jobj,STEGO);
  307. catch
  308.     error('ERROR (problem with saving the stego image)')
  309. end
  310. end
  311. function res=invH(y)
  312. to_minimize=@(x) (H(x)-y)^2;
  313. res=fminbnd(to_minimize,eps,0.5-eps);
  314. end
  315. function res=H(x)
  316. res=-x*log2(x)-(1-x)*log2(1-x);
  317. end

JSTEG内容提取:

  1. STEGO='stego.jpg';
  2. output=fopen("secret.txt","w");
  3. SEED=99;
  4. mlen=80;
  5. tic;
  6. messageste=nsf5_extract(STEGO,mlen);
  7. length(messageste)
  8. for j=1:length(messageste)
  9.     if(rem(j,8)==0)
  10.         fprintf(output,"%d",messageste(j));
  11.         %fprintf(output,'\r\n');
  12.         continue
  13.     end
  14.     fprintf(output,"%d",messageste(j));
  15. end
  16. %fscanf(output,'%s',messageste);
  17. T=toc;
  18. fprintf('-------\n');
  19. fprintf('nsF5 extract finished\n');
  20. fprintf('elapsed time:%.4f seconds\n',T);
  21. fclose(output);
  22. function message=nsf5_extract(STEGO,mlen)
  23. try
  24.     jobj=jpeg_read(STEGO);
  25.     DCT=jobj.coef_arrays{1};
  26. catch
  27.     error('ERROR (problem with the STEGO image)');
  28. end
  30. AC=numel(DCT)-numel(DCT(1:8:end,1:8:end));
  31. idD=1;
  32. for id=1:mlen
  33.     while(abs(DCT(idD))<=1)
  34.         idD=idD+1;
  35.     end
  36.     if(DCT(idD)>0)
  37.     if(mod(DCT(idD),2)==1)
  38.         message(1,id)=1;
  39.     else
  40.         message(1,id)=0;
  41.     end
  42.     else
  43.           if(mod(DCT(idD),2)~=1)
  44.             message(1,id)=1;
  45.           else
  46.             message(1,id)=0;
  47.           end 
  48.     end
  49.     idD=idD+1;
  50. end
  51. end

F4内容提取

  1. STEGO='stego.jpg';
  2. output=fopen("secret.txt","w");
  3. SEED=99;
  4. mlen=80;
  5. tic;
  6. messageste=nsf5_extractF4(STEGO,mlen);
  7. length(messageste)
  8. for j=1:length(messageste)
  9.     if(rem(j,8)==0)
  10.         fprintf(output,"%d",messageste(j));
  11.         %fprintf(output,'\r\n');
  12.         continue
  13.     end
  14.     fprintf(output,"%d",messageste(j));
  15. end
  16. %fscanf(output,'%s',messageste);
  17. T=toc;
  18. fprintf('-------\n');
  19. fprintf('nsF5 extract finished\n');
  20. fprintf('elapsed time:%.4f seconds\n',T);
  21. fclose(output);
  22. function message=nsf5_extractF4(STEGO,mlen)
  23. try
  24.     jobj=jpeg_read(STEGO);
  25.     DCT=jobj.coef_arrays{1};
  26. catch
  27.     error('ERROR (problem with the STEGO image)');
  28. end
  30. AC=numel(DCT)-numel(DCT(1:8:end,1:8:end));
  31. idD=1;
  32. for id=1:mlen
  33.     while(DCT(idD)==0)
  34.         idD=idD+1;
  35.     end
  36.     DCT(idD)
  37.     if(DCT(idD)<0)
  38.     if(mod(DCT(idD),2)==1)
  39.         message(1,id)=0;
  40.     end
  41.     if(mod(DCT(idD),2)==0)
  42.         message(1,id)=1;
  43.     end
  44.     
  45.     end
  46.     if(DCT(idD)>0)
  47.     if(mod(DCT(idD),2)==1)
  48.         message(1,id)=1;
  49.     end
  50.     if(mod(DCT(idD),2)==0)
  51.         message(1,id)=0;
  52.     end
  53.      end
  54.     
  55.     idD=idD+1;
  56. end
  57. end

F5内容提取

  1. STEGO='stego.jpg';
  2. output=fopen("secret.txt","w");
  3. SEED=99;
  4. mlen=80;
  5. tic;
  6. messageste=nsf5_extractF5(STEGO,mlen);
  7. length(messageste)
  8. for j=1:length(messageste)
  9.     if(rem(j,8)==0)
  10.         fprintf(output,"%d",messageste(j));
  11.         %fprintf(output,'\r\n');
  12.         continue
  13.     end
  14.     fprintf(output,"%d",messageste(j));
  15. end
  16. %fscanf(output,'%s',messageste);
  17. T=toc;
  18. fprintf('-------\n');
  19. fprintf('nsF5 extract finished\n');
  20. fprintf('elapsed time:%.4f seconds\n',T);
  21. fclose(output);
  22. function message=nsf5_extractF5(STEGO,mlen)
  23. mlen=mlen+3-mod(mlen,3);
  24. try
  25.     jobj=jpeg_read(STEGO);
  26.     DCT=jobj.coef_arrays{1};
  27. catch
  28.     error('ERROR (problem with the STEGO image)');
  29. end
  31. AC=numel(DCT)-numel(DCT(1:8:end,1:8:end));
  32. idD=1;
  33. for id=1:2:mlen
  34.     while(DCT(idD)==0)
  35.         idD=idD+1;
  36.     end
  37.     value1=DCT(idD)
  38.     idD=idD+1;
  39.     while(DCT(idD)==0)
  40.         idD=idD+1;
  41.     end
  42.     value2=DCT(idD)
  43.     idD=idD+1;
  44.     while(DCT(idD)==0)
  45.         idD=idD+1;
  46.     end
  47.     value3=DCT(idD)
  48.     idD=idD+1;
  49.     if(value1>0)
  50.         m1=mod(value1,2);
  51.     else
  52.         m1=mod(value1+1,2);
  53.     end    
  54.     if(value2>0)
  55.         m2=mod(value2,2);
  56.     else
  57.         m2=mod(value2+1,2);
  58.     end
  59.     if(value3>0)
  60.         m3=mod(value3,2);
  61.     else
  62.         m3=mod(value3+1,2);
  63.     end
  64.     message(id)=mod((m1+m2),2);
  65.     message(id+1)=mod((m2+m3),2);
  66.     
  67. end
  68. end
标签： matlab
本文转载自: https://blog.csdn.net/weixin_45695828/article/details/124211945
版权归原作者 弗里曼博士 所有， 如有侵权，请联系我们删除。
登录后发布评论

“[HUST]多媒体安全实验：图像隐写JSTEG，F4,F5的实现”的评论:

还没有评论
关于作者
...
overfit同步小助手
文章同步
相关阅读

网络安全法-网络运行安全

使用selenium/drissionpage时如何阻止chrome自动跳转http到https

docker desktop 里部署的Open WebUI 管理员密码忘记了的处理方法

在ubuntu20.04中搭建onsite比赛运行环境

利用开源的低代码表单设计器FcDesigner高效管理和渲染复杂表单结构

Kafka学习笔记

【前端】浏览器输入url到页面呈现发生了什么？