在实现提取水印的过程中，遇到了一些问题
首先还是根据论文中的思路来梳理一下整个提取流程
- 读入两个模型，一个原始模型ori_mesh， 一个水印模型wm_mesh。
- 将两个模型对齐（即放在同一个坐标系下）
- 分别计算两个模型的频谱系数，记ori_mesh频谱系数为Rs, Rt, Ru, 记wm_mesh的频谱系数为wRs, wRt, wRu.
- 根据下式计算Qj
- 再根据下式计算出嵌入水印的符号，这里相当于得到的是嵌入的水印b‘，只是记为了aj
- 最后根据下式得出原始水印序列b
?
根据上述步骤，我一开始的想法是在网格平台下读入两个模型，然后按照上述步骤进行计算会遇到一些麻烦，比如获取网格的时候，虽然可以用一个链表将打开的文件存放起来，但是在初始化网格的时候，由于参数都是单独的，如果全部都改写成数组或者是链表的话，那改动起来会相当麻烦。于是我决定在嵌入水印的时候，将R矩阵写入文件中，这样提取水印的时候就不必再读入原始网格，只需计算出读入的水印网格中的wR矩阵，然后从文件中读取出R矩阵，两者做差值，再按照原来的步骤进行判断即可。
?
写入和读取分别如下：
//将原始网格的频谱系数写入文件中Rs????????
ofstream
Rsfile;
????????Rsfile.open("D:\\Rs.txt",ios_base::out );
????????if ( Rsfile)
????????{
????????????for ( int
i = 0; i < m_vertexNum; i++)
????????????{
????????????Rsfile<<Rs[i]<<" ";
????????????}
????????}
????????Rsfile.close();
?
????//将原始网格的频谱系数从文件中读取出来
????????double *pRs = new
double[m_vertexNum];
????????double *pRt = new
double[m_vertexNum];
????????double *pRu = new
double[m_vertexNum];
?
????????//Rs
????????ifstream
Rsfile;
????????Rsfile.open("D:\\Rs.txt",ios_base::in );
????????if (Rsfile)
????????{
????????????// read into memory
????????????Rsfile.seekg (0, Rsfile.end);
????????????int
length = Rsfile.tellg();
????????????Rsfile.seekg (0, Rsfile.beg);
?
????????????char *buffer = new
char[length];
????????????Rsfile.read(buffer, length);
????????????Rsfile.close();
?
????????????// parse into array
????????????std::istringstream
iss(buffer);
????????????int
i = 0;
????????????while (iss >> pRs[i++]);
????????????delete [] buffer;
????????????// print or use it.
????????}
????????Rsfile.close();
在提取的过程中，主要的问题出在sign函数上，由于没有查到sign函数所属的头文件，因此我自己实现了一个sign函数，代码如下：
double
EigenDeformation::sign( double
x)
????{
????????cout <<"x = "<< x <<endl;
????????if( x > 0)
????????????return 1;
????????else
if(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
?