数字盲水印工具dct.exe

happy886rr

http://bcn.bathome.net/s/tool/index.html?key=dct

使用了DCT离散余弦变换，可以在图片中嵌入全盲数字水印，利用模拟量子态，成功解决信息的稳定性。图片即使被调色，被ps处理，其嵌入的数字水印依然可读取。本工具采取非开源策略。

使用说明:

1. 
   
2. 加数字盲水印:
   
3.     dct [输入图片] [输出图片] [水印字符串]
   
4. 解数字盲水印:
   
5.     dct [输入图片]
   

复制代码

核心函数

1. // 离散余弦变换核心
   
2. int DCT2D(double* cMatrix, double* pBmpMatrix, double* cMatrixT, int Nx, int Ny, int N, DCT2Mode mode, byte* insertData, int insertLen)
   
3. {
   
4.         // 如果在编码模式下
   
5.         if(mode == DCT2_ENCODE_MODE)
   
6.         {
   
7.                 // 检测插入长度是否合法
   
8.                 if
   
9.                 (
   
10.                     (insertLen <= 0)          ||
    
11.                     (insertLen > Nx * Ny / 8) ||
    
12.                     (Nx == 0)                 ||
    
13.                     (Ny == 0)
    
14.                 )
    
15.                 {
    
16.                         // 插入数据量过大，则退出
    
17.                         fprintf(stdout, "You can't inserted bytes to the picture\n");
    
18.                         exit(1);
    
19.                 }
    
20.         }
    
21. 
    
22.         // 辅助数组
    
23.         double* pTmp = (double*)malloc(N * N * sizeof(double));
    
24.         int pBmpMatrixWith = Nx * N;
    
25. 
    
26.         // 插入索引位
    
27.         int insertIndex = 0;
    
28. 
    
29.         // 分块DCT变换
    
30.         for(int j = 0; j < Ny; j ++)
    
31.         {
    
32.                 for(int i = 0; i < Nx; i ++)
    
33.                 {
    
34.                         int xOffset = i * N;
    
35.                         int yOffset = j * N * pBmpMatrixWith;
    
36. 
    
37.                         // 正交化乘以系数矩阵
    
38.                         for(int dy = 0; dy < N; dy ++)
    
39.                         {
    
40.                                 for(int dx = 0; dx < N; dx ++)
    
41.                                 {
    
42.                                         double sum = 0.0f;
    
43.                                         for(int sn = 0; sn < N; sn ++)
    
44.                                         {
    
45.                                                 sum += (*(cMatrix + (dy * N + sn))) * (*(pBmpMatrix + (yOffset + sn * pBmpMatrixWith) + (xOffset+ dx)));
    
46.                                         }
    
47. 
    
48.                                         *(pTmp + (dy * N + dx)) = sum;
    
49.                                 }
    
50.                         }
    
51. 
    
52.                         // 正交化乘以转置系数矩阵
    
53.                         for(int dy = 0; dy < N; dy ++)
    
54.                         {
    
55.                                 for(int dx = 0; dx < N; dx ++)
    
56.                                 {
    
57.                                         double sum = 0.0f;
    
58.                                         for(int sn = 0; sn < N; sn ++)
    
59.                                         {
    
60.                                                 sum += (*(pTmp + (dy * N + sn))) * (*(cMatrixT + (sn * N) + dx));
    
61.                                         }
    
62. 
    
63.                                         *(pBmpMatrix + (yOffset + dy * pBmpMatrixWith) + (xOffset+ dx)) = sum;
    
64.                                 }
    
65.                         }
    
66. 
    
67.                         // 逆变换模式
    
68.                         if(mode == DCT2_INVERSE_TRANSFORM)
    
69.                         {
    
70.                                 continue;
    
71.                         }
    
72. 
    
73.                         // 编码模式
    
74.                         if(mode == DCT2_ENCODE_MODE)
    
75.                         {
    
76.                                 double* pA = pBmpMatrix + (yOffset + 5 * pBmpMatrixWith) + (xOffset+ 5);
    
77.                                 double* pB = pBmpMatrix + (yOffset + 6 * pBmpMatrixWith) + (xOffset+ 6);
    
78. 
    
79.                                 // 采集密文bit位
    
80.                                 int insertV = ((int)(*(insertData + insertIndex / 8))) & ((int)(0x80 >> (insertIndex & 7)));
    
81. 
    
82.                                 if(insertV == 0)
    
83.                                 {
    
84.                                         if(*pA > *pB)
    
85.                                         {
    
86.                                                 SWAP_DOUBLE(*pA, *pB);
    
87.                                         }
    
88. 
    
89.                                         *pB += 1.0f;
    
90.                                 }
    
91.                                 else
    
92.                                 {
    
93.                                         if(*pA < *pB)
    
94.                                         {
    
95.                                                 SWAP_DOUBLE(*pA, *pB);
    
96.                                         }
    
97.                                         *pB -= 1.0f;
    
98.                                 }
    
99. 
    
100.                                 // 插入计数器
     
101.                                 insertIndex ++;
     
102.                                 continue;
     
103.                         }
     
104. 
     
105.                         // 解码模式
     
106.                         if(mode == DCT2_DECODE_MODE)
     
107.                         {
     
108.                                 double* pA = pBmpMatrix + (yOffset + 5 * pBmpMatrixWith) + (xOffset+ 5);
     
109.                                 double* pB = pBmpMatrix + (yOffset + 6 * pBmpMatrixWith) + (xOffset+ 6);
     
110.                                 byte* pInsertV = (byte*)(insertData + (insertIndex >> 3));
     
111. 
     
112.                                 if(*pA > *pB)
     
113.                                 {
     
114.                                         *pInsertV |= (byte)(0x80 >> (insertIndex & 7));
     
115.                                 }
     
116. 
     
117.                                 // 插入计数器
     
118.                                 insertIndex ++;
     
119.                                 continue;
     
120.                         }
     
121.                 }
     
122.         }
     
123. 
     
124.         // 释放辅助数组
     
125.         free(pTmp);
     
126.         return 0;
     
127. }

复制代码

CrLf

研究了半天代码，看不懂...放弃

happy886rr

回复 2# CrLf
离散余弦变换，就是把密文转化为图片频域能量添加到图片上，再逆变换成图片，参见这个帖子https://www.cnblogs.com/zkwarrior/p/5980191.html
说白了就是防止盗图用的，我的代码中只用了1.0f倍的能量，如果改成1000被，那么别人把图片压缩，调色，密文依然能还原，也就是密文不是嵌入到二进制数据里，而是嵌入到图片的灰度频率里。你通过ps软件对图片压缩、修改、调色、甚至直接拍照， 其解密的密文依然不变。 用于网站防盗图、信息隐藏等方面。

CrLf

66666，其实看过这篇（好像是知乎），但是还是很难理解
我原以为会在核心代码里看到 sin，然而并没有

happy886rr

回复 4# CrLf
调用数学库函数慢，都提前求好，转化为整数。因为它是分块处理的，就是每次处理一个8*8的小格子，那些cos系数提前算出来并量化为整数比值，这样可以用整形替代浮点型。
说白了就是cos之类的，提前就算好，存在一个static const int 系数[8][8]的数组里。

CrLf

回复 5# happy886rr


    soga

zgzxp

加入水印的文件再截图后，解密不了原水印