[发明专利]一种基于DCT变换的矢量地图数据盲水印方法

摘要

本发明公开了一种基于DCT变换的矢量地图数据盲水印方法。传统水印算法还很难直接应用在矢量地图数据上。本发明方法的包括水印信息嵌入和水印信息提取。水印信息嵌入是把用户选择的水印信息嵌入到原始矢量地图中得到一个新的矢量地图，具体包括调制水印信息、对矢量地图数据进行分块处理、和嵌入水印信息。水印信息提取是从待检测的矢量地图中提取水印信息内容，具体包括对待检测的矢量地图数据进行分块处理、提取水印信息和生成水印信息图。本发明方法具有较好的鲁棒性，并且是一种盲检算法，具有较好的实用价值。

主权项

一种基于DCT变换的矢量地图数据盲水印方法，其特征在于该方法的包括水印信息嵌入和水印信息提取；所述的水印信息嵌入是把用户选择的水印信息嵌入到原始矢量地图中得到一个新的矢量地图，具体步骤包括：步骤(1)调制水印信息，二值水印图结合密钥和运用Arnold变换把水印信息置乱，得到二进制比特流w：设定水印图像为N×N的二值图，则水印信息长度为l＝N×N，采用离散化的Arnold变换并结合密钥得到最终水印信息比特流w＝(w1，w2，...，wi，...，wl)，其中wi∈{0，1}，0≤i＜l；步骤(2)对矢量地图数据进行分块处理：运用四叉树算法把地图划分成矩形网格，使得每个矩形网格内所含的顶点数不小于水印信息长度l×α，并对划分所得的树状结构矩形网格按深度优先排序，其中α表示嵌入每位水印信息所需要的载体数据个数；步骤(3)嵌入水印信息：按照深度优先顺序遍历每一块网格，采用步骤a至步骤e完成水印嵌入，具体步骤是：步骤a.把网格内的所有顶点进行排序，排序规则为：从网格的左上角开始，按照先从左到右、再从上到下的顺序对顶点进行排序；步骤b.截取顶点坐标：对排好序的所有顶点坐标值减去该网格内最左下顶点坐标值，得到一组顶点坐标差值，并对这组顶点坐标差值放大m倍，m取10、100、1000或10000；步骤c.DCT变换：把连续α个顶点坐标差值进行一维DCT变换，得到一组长度为α的DCT系数值d，保留d0和da-1取剩下系数值的最小值dmin，即dmin＝min{d1，...，da-2}，然后水印信息位wi和系数值关系做以下判断操作：如果wi等于1，而且fabs(dmin)＞fabs(da-1)，则da～1＝fabs(dmin)×β×sign(da-1)；如果wi等于1，而且fabs(dmin)≤fabs(da-1)，则da-1保持不变；如果wi等于0，而且fabs(dmin)＜fabs(da-1)，则da-1＝fabs(dmin)×1/β×sign(da-1)；如果wi等于0，而且fabs(dmin)≥fabs(da-1)，则da-1保持不变；其中，函数fabs表示取绝对值，函数sign表示取数值的正负符号，β表示放大倍数，β＞1；步骤d.重复步骤c直到所有水印信息嵌入到相应的DCT系数中去；步骤e.DCT反变换：依次把嵌入水印信息的DCT系数做反变换，顶点坐标差值缩小m倍，将缩小后的顶点坐标差值加上网格内最左下顶点坐标值，得到新的坐标值数据；所述的水印信息提取是从待检测的矢量地图中提取水印信息内容，包括以下步骤：步骤(I)对待检测的矢量地图数据进行分块处理：运用四叉树算法把地图划分成矩形网格，使得每个矩形网格内所含的顶点数不小于水印信息长度l×α，并对划分所得的树状结构矩形网格按深度优先排序；步骤(II)提取水印信息：遍历每一块网格块，采用步骤f至步骤i提取水印，具体步骤是：步骤f.把网格内的所有顶点进行排序；排序规则为：从网格的左上角开始，按照先从左到右、再从上到下的顺序对顶点进行排序；步骤g.截取顶点坐标：对排好序的所有顶点坐标值减去该网格内最左下顶点坐标值，得到一组顶点坐标差值，并对这组顶点坐标差值放大m倍，m取10、100、1000或10000；步骤h.DCT变换：把连续α个顶点坐标差值进行一维DCT变换，得到一组长度为α的DCT系数值d，保留d0和da-1取剩下系数值的最小值dmin，即dmin＝min{d1，...，da-2}，通过判断得到水印信息位wi，判断规则为：如果fabs(dmin)＞fabs(da-1)，则wi＝0；fabs(dmin)≤fabs(da-1)，则wi＝1；步骤i.重复步骤h直到从相应的DCT系数中提取出所有水印信息；步骤(III)生成水印信息图：把提取出的水印信息，通过密钥和Arnold反置乱变换操作，重新生成得到二值水印图。