C++实现网格水印之调试笔记（五）—— 提取出错

在实现提取水印的过程中，遇到了一些问题

首先还是根据论文中的思路来梳理一下整个提取流程

1. 读入两个模型，一个原始模型ori_mesh， 一个水印模型wm_mesh。
2. 将两个模型对齐（即放在同一个坐标系下）
3. 分别计算两个模型的频谱系数，记ori_mesh频谱系数为Rs, Rt, Ru, 记wm_mesh的频谱系数为wRs, wRt, wRu.
4. 根据下式计算Qj

   

5. 再根据下式计算出嵌入水印的符号，这里相当于得到的是嵌入的水印b‘，只是记为了aj

   

6. 最后根据下式得出原始水印序列b

   

?

根据上述步骤，我一开始的想法是在网格平台下读入两个模型，然后按照上述步骤进行计算会遇到一些麻烦，比如获取网格的时候，虽然可以用一个链表将打开的文件存放起来，但是在初始化网格的时候，由于参数都是单独的，如果全部都改写成数组或者是链表的话，那改动起来会相当麻烦。于是我决定在嵌入水印的时候，将R矩阵写入文件中，这样提取水印的时候就不必再读入原始网格，只需计算出读入的水印网格中的wR矩阵，然后从文件中读取出R矩阵，两者做差值，再按照原来的步骤进行判断即可。

?

写入和读取分别如下：

//将原始网格的频谱系数写入文件中Rs????????

ofstream Rsfile;

????????Rsfile.open("D:\\Rs.txt",ios_base::out );

????????if ( Rsfile)

????????{

????????????for ( inti = 0; i < m_vertexNum; i++)

????????????{

????????????Rsfile<<Rs[i]<<" ";

????????????}

????????}

????????Rsfile.close();

?

????//将原始网格的频谱系数从文件中读取出来

????????double *pRs = newdouble[m_vertexNum];

????????double *pRt = newdouble[m_vertexNum];

????????double *pRu = newdouble[m_vertexNum];

?

????????//Rs

????????ifstreamRsfile;

????????Rsfile.open("D:\\Rs.txt",ios_base::in );

????????if (Rsfile)

????????{

????????????// read into memory

????????????Rsfile.seekg (0, Rsfile.end);

????????????intlength = Rsfile.tellg();

????????????Rsfile.seekg (0, Rsfile.beg);

?

????????????char *buffer = newchar[length];

????????????Rsfile.read(buffer, length);

????????????Rsfile.close();

?

????????????// parse into array

????????????std::istringstreamiss(buffer);

????????????inti = 0;

????????????while (iss >> pRs[i++]);

????????????delete [] buffer;

????????????// print or use it.

????????}

????????Rsfile.close();

在提取的过程中，主要的问题出在sign函数上，由于没有查到sign函数所属的头文件，因此我自己实现了一个sign函数，代码如下：

double EigenDeformation::sign( doublex)

????{

????????cout <<"x = "<< x <<endl;

????????if( x > 0)

????????????return 1;

????????elseif(fabs(x) < 0.000000000001)

????????????return 0;

????????else

????????????return -1;

????}

由于精度的关系，sign函数判断的符号总是不正确，调试查看输入sign函数的值和其输出的只如下：

其中x是输入的值，P是判断后得到的值。

然后常识将输入扩大1kw倍也失败。

改用另外一个函数 ? double _copysign( double x, double y );

这个函数是C中math.h中的函数，函数功能： 以第二个参数y的符号（正或负）返回第一个参数x，于是在代码中调用形式如下：

wB[i] = _copysign(1, Q[i]) ;

结果如下：

返回符号成功，对比嵌入的水印序列

怎么两者差别这么大。。。

每次运行得到的P的结果都一样。

?

?

错误1：将Rt读入文件的时候，写成Rs了

错误2：顶点赋值方式错误

更正：

原始网格的R 水印网格的R

?

原文：http://www.cnblogs.com/scut-linmaojiang/p/4479081.html