专利名称：一种基于栅格地图特征的小波域零水印方法
技术领域：
本发明属于多媒体信息安全技术领域，特别是针对栅格数字地图数据，提出一种 基于栅格数字地图特征的零水印方法。
背景技术：
目前已经提出了许多数字水印方法，这些方法一般可以分为两类频域水印法和 空域水印法。这些嵌入数字水印的方法都对数据的空域信息或其频域信息做了一定的修改 来嵌入水印信息，为了不让人眼发觉人为修改的痕迹，很多方法采用了基于HVS (人类视觉 系统)的视觉掩膜。用加视觉掩膜的方法在一定程度上解决了水印可感知性和鲁棒性之间 的矛盾。但是加视觉掩膜使得水印的嵌入过程复杂化，消耗计算时间太长，不利于实际应 用，而且如果别有用心者了解加视觉掩膜的方法，不能排除有能力的人用方法篡改含有水 印信息的地图，同样使人眼难以发觉的细小的改动，从而导致水印信息检测失败。因此，这 就使得数字水印的安全性受到了限制。所以最好的办法是运用零水印来解决。所谓的零水 印就是不对原始图像或数据进行修改的一种新的数字水印技术。它是利用图像或者数据的 重要特征来构造水印信息，而不是来修改这些特征。所以零水印技术很好地解决了不可见 数字水印的可感知性和鲁棒性之间的矛盾。特别是针对栅格地图的数据要求，由于栅格的 地理空间数据越来越被社会运用，它是描述关于人类赖以生存的地球的重要信息，是国家 基础性、公益性设施，是国家基础设施建设和地球科学研究的支撑性成果，是国家经济、国 防建设中不可缺少的资源，对社会持续发展起着重要作用，因此零水印的应用会对于栅格 地图数据的版权保护起到至关重要的作用。因此零水印非常适用于栅格地图的安全保护， 即不修改原始地图的内容要能满足版权的保护。
发明内容
本发明针对现有技术的不足，提供了一种基于二维栅格数字地图特征的零水印方法。本发明方法包括水印信息嵌入和水印信息检测。水印信息嵌入
所述的水印信息嵌入是指提取原始栅格地图的零水印信息，保存入库，具体包括以下 步骤
步骤(1)调制水印信息，二值水印图结合密钥和运用Arnold变换把水印信息置乱，得 到二进制比特流w，具体是设定水印图像为NXN的二值图，则水印信息长度为/=NXN，采 用离散化的Arnold变换并结合密钥得到最终水印信息比特流W=(Wl，w2, . . . , Wi, . . . , W1), 其中Wi e {0，1}，0彡Kl ；同时生成两个长度为1的二进制比特流X=(Xl，X2,... ,Xi,..., X1)和y=(yi，y2，...⑶，……)，用来保存特征数据，其中Xi=Yi=O, 0 ^ i<7 ；
步骤( 对栅格地图数据进行分块处理根据水印信息的长度Λ把栅格地图数据分成 1 土夬，且每块大小为MXM方形小块；然后进行以下步骤步骤a、对每个MXM的小块数据进行二维sym2小波分解(sym2为Symlets小波系中的 一种)，分解级数为2;
步骤b、对第i块小波系数中的LIA和Hffii系数(LIA代表第2级小波分解中的逼近 系数，Hffii代表第2级小波分解中的对角线细节系数)分别进行统计取LI^中的最小值与 最大值，均分成十个等间隔的区间，把LI^中的元素分到此十个区间之内，并返回每个范围 内元素的个数作为一个一维向量，记为V1，再取Hffii中的最小值与最大值，均分成十个等间 隔的区间，把Hffii中的元素分到此十个区间之内，并返回每个范围内元素的个数作为一个 一维向量，记为V2;
步骤c、分别对向量V1和V2的元素进行从大到小的排序，记为V1 ’和V2'，获得V1'排在向 量前面的两个元素ApB1',以及他们在向量V1中的位置HA1、HB1，获得V2’排在向量前面的两 个元素A2’、B2’，以及他们在向量V2中的位置HA2、HB2，然后做如下处理 如果 abs (Hai-Hbi) >T，则令 Xi=I ； 如果 abs (Hai-Hbi) < T，则令 Xi=O ； 如果 abs (Ha2-Hb2) >T，则令 yi=l ； 如果 abs (Ha2-Hb2) < T，则令 Yi=O ； 其中，abs ()代表取绝对值，T代表阈值。步骤( 生成零水印信息把置乱后的水印信息二进制比特流w和X进行异或操 作得到零水印信息I1, 二进制比特流w和y进行异或操作得到零水印信息I2 ；
步骤(4)保存零水印信息把生成的零水印信息I1和I2保存到水印信息库中，以备检 测时需要。水印信息检测
所述的水印信息检测是指将待检测栅格地图中提取的零水印信息与水印信息库中的 水印信息操作得到最后水印信息，具体步骤如下
步骤⑴根据待检测栅格地图数据，取出水印信息库中的水印信息IJPI2，并生成两个 长度为 7 的二进制比特流 α=(αι，α2，···，Qi,..., CI1)禾口 3=(3i，β2，···，^i,..., β D，用来保存特征数据，其中α严1=0，0彡iO;；
步骤( 对待检测栅格地图数据进行分块处理根据水印信息的长度，把栅格地图数 据分成每块大小为MXM的小块；然后进行一下步骤
步骤a、对每个MXM的小块数据进行二维sym2小波分解(sym2为Symlets小波系中的 一种)，分解级数为2;
步骤b、对第i块小波系数中的LI^和Hffii系数(LIA代表第2级小波分解中的逼近 系数，Hffii代表第2级小波分解中的对角线细节系数)分别进行统计取LI^中的最小值与 最大值，均分成十个等间隔的区间，把LI^中的元素分到此十个区间之内，并返回每个范围 内元素的个数作为一个一维向量，记为V3，再取Hffii中的最小值与最大值，均分成十个等间 隔的区间，把Hffii中的元素分到此十个区间之内，并返回每个范围内元素的个数作为一个 一维向量，记为V4;
步骤c、分别对向量V3和V4的元素进行从大到小的排序，记为V3’和V4',获得V3'排在向 量前面的两个元素C3'、D3'，以及他们在向量V3中的位置H。3、HD3，获得V4'排在向量前面的两 个元素C4’、D4',以及他们在向量V4中的位置Hc4、Hd4，然后做如下处理
5如果 abs (Hg3-Hd3) >T，则令 Qi=I ； 如果 abs (Hc3-Hd3) < T，则令 α i=0 ； 如果 abs (Hc4-HJ >T,则令 ^i=I ； 如果 abs (Hc4-H 4) < T 则令 β 尸0 ； 其中，abs()代农収绝对值，T代表阈值。步骤( 得到检测信息把α和I1进行异或操作，得到W1 ’，把β和I2进行异或 操作，得到W2'；
步骤⑷反置乱水印信息把二进制比特流Wp W；结合密钥进行Arnold反变换提取 出最终的水印图像I1 ’和I2'。本发明的有益效果为
该方法不需要对原始栅格地图数据做任何修改，根据栅格地图数据图像块颜色差别比 较明显的特征，对栅格地图数据做小波变换，并在小波域系数上提取零水印信息。该水印方 法，一方面避免把水印作为冗余信息嵌入在地图数据中引起数据改变的缺点，具有较好的 隐蔽性，可以很好地解决数字水印的不可感知性和鲁棒性之间的矛盾，也可以克服可逆数 字水印中存在的安全漏洞，是一种天然的盲水印系统，有很大的实用价值；另一方面能够抵 抗多种常见的攻击，比如缩放、平移、旋转、裁剪，具有较强的抵抗力和稳健性。
具体实施例方式本发明方法包括水印信息嵌入和水印信息检测。水印信息嵌入
所述的水印信息嵌入是指提取原始矢量地图的零水印信息，保存入库，具体包括以下 步骤
步骤(1)调制水印信息，二值水印图结合密钥和运用Arnold变换把水印信息置乱，得 到二进制比特流w，具体是设定水印图像为NXN的二值图，则水印信息长度为/=NXN，采 用离散化的Arnold变换并结合密钥得到最终水印信息比特流W=(Wl，w2, . . . , Wi, . . . , W1), 其中Wi e {0，1}，0彡Kl ；同时生成两个长度为1的二进制比特流X=(Xl，X2,... ,Xi,..., X1)和Y=G^y2,...…，…，Y1)，用来保存特征数据，其中Xi = Yi =0, 0 ^ i<7 ；
步骤( 对栅格地图数据进行分块处理根据水印信息的长度Λ把栅格地图数据分成 1块,且每块大小为MXM方形小块；然后进行以下步骤
步骤a、对每个MXM的小块数据进行二维sym2小波分解(sym2为Symlets小波系中的 一种)，分解级数为2;
步骤b、对第i块小波系数中的LIA和Hffii系数(LIA代表第2级小波分解中的逼近 系数，HH2,代表第2级小波分解中的对角线细节系数)分别进行统计将LI^和Hffii中的 元素分别分到10个等间隔范围内(取最小值与最大值，均分成十个等间隔的区间)，并返回 每个范围内元素的个数作为一个一维向量，分别记为V1和V2 ；
步骤c、分别对向量V1和V2的元素进行从大到小的排序，记为V1'和V2’，分别获得排在 向量前面的两个元素ApB1'和A2’、4’，以及他们在向量V1和V2中的位置Hai、Hbi和Ha2、Hb2， 然后做如下处理
如果 abs (Hai-Hbi) >T，则令 Xi=I ；如果 abs (Hai-Hbi) < T，则令 Xi=O ； 如果 abs (Ha2-Hb2) >T，则令 yi=l ； 如果 abs (Ha2-Hb2) < T,则令 Yi=O ；
其中，abs ()代表取绝对值，T代表阈值，对T取不同值，经不断实验，表明阈值T=2时效 果最好。步骤( 生成零水印信息把置乱后的水印信息二进制比特流w和X进行异或操 作得到零水印信息I1, 二进制比特流w和y进行异或操作得到零水印信息I2 ；
步骤(4)保存零水印信息把生成的零水印信息I1和I2保存到水印信息库中，以备检 测时需要。水印信息检测
所述的水印信息检测是指将待检测矢量地图中提取的零水印信息与水印信息库中的 水印信息操作得到最后水印信息，具体步骤如下
步骤⑴根据待检测栅格地图数据，取出水印信息库中的水印信息IdPI2，并生成两个 长度为 7 的二进制比特流 α β=(α” α2，··.，％，··.，CI1)禾口 β=(β” β2，··.，βρ···， β i)，用来保存特征数据，其中％ = β , =0, 0 ^ Kl ；；
步骤( 对待检测栅格地图数据进行分块处理根据水印信息的长度，把栅格地图数 据分成每块大小为MXM的小块；然后进行以下步骤
步骤a、对每个MXM的小块数据进行二维sym2小波分解(sym2为Symlets小波系中的 一种)，分解级数为2;
步骤b、对第i块小波系数中的LIA和Hffii系数(LIA代表第2级小波分解中的逼近 系数，HH2,代表第2级小波分解中的对角线细节系数)分别进行统计将LI^和Hffii中的 元素分别分到10个等间隔范围内(取最小值与最大值，均分成十个等间隔的区间)，并返回 每个范围内元素的个数作为一个一维向量，分别记为V3和V4 ；
步骤c、分别对向量V3和V4的元素进行从大到小的排序，记为V3'和V4’，分别并获得排 在前面的两个元素C3’、D3’和C4’、D4’，以及他们在向量V3和V4中的位置HC3、HD3和HC4、HD4，然 后做如下处理
如果 abs (Hc3-HdJ >T,则令 α i=l ； 如果 abs (Hc3-HdJ < T，则令 α i=0 ； 如果 abs (Hc4-Hd4) >T,则令 β 尸1 ； 如果 abs (Hc4-Hd4) < T，则令 β 尸0 ； 其中，abs ()代表取绝对值，T代表阈值。步骤( 得到检测信息把α和I1进行异或操作，得到W1 ’，把β和I2进行异或 操作，得到W2'；
步骤⑷反置乱水印信息把二进制比特流Wl’、W2'结合密钥进行Arnold反变换提取 出最终的水印图像。对于现在的某些针对零水印的攻击中，有可能是破坏LI^系数中特征， 也有可能是破坏Hffii系数中的特征，最后检测结果两张水印图像中只要有一张水印图像能 够识别，即可说明零水印是存在的。
权利要求
1. 一种基于栅格地图特征的小波域零水印方法，其特征在于该方法包括水印信息嵌入 和水印信息检测；所述的水印信息嵌入是指提取原始栅格地图的零水印信息，保存入库，具体包括以下 步骤步骤(1)调制水印信息，二值水印图结合密钥和运用Arnold变换把水印信息置乱，得 到二进制比特流w，具体是设定水印图像为NXN的二值图，则水印信息长度为/=NXN，采 用离散化的Arnold变换并结合密钥得到最终水印信息比特流W=(Wl，w2, . . . , Wi, . . . , W1), 其中Wi e {0，1}，0彡Kl ；同时生成两个长度为1的二进制比特流X=(Xl，X2,... ,Xi,..., X1)和y=(yi，y2，...⑶，……)，用来保存特征数据，其中Xi=Yi=O, 0 ^ i<7 ；步骤( 对栅格地图数据进行分块处理根据水印信息的长度Λ把栅格地图数据分成 1 土夬，且每块大小为MXM方形小块；然后进行以下步骤步骤a、对每个MXM的小块数据进行二维sym2小波分解，分解级数为2 ； 步骤b、对第i块小波系数中的LI^和Hffii系数分别进行统计取LI^中的最小值与 最大值，均分成十个等间隔的区间，把LI^中的元素分到此十个区间之内，并返回每个范围 内元素的个数作为一个一维向量，记为V1，再取Hffii中的最小值与最大值，均分成十个等间 隔的区间，把Hffii中的元素分到此十个区间之内，并返回每个范围内元素的个数作为一个 一维向量，记为V2;步骤c、分别对向量V1和V2的元素进行从大到小的排序，记为V1 ’和V2'，获得V1'排在向 量前面的两个元素ApB1',以及他们在向量V1中的位置HA1、HB1，获得V2’排在向量前面的两 个元素A2’、B2’，以及他们在向量V2中的位置HA2、HB2，然后做如下处理 如果 abs (Hai-Hbi) >T，则令 Xi=I ； 如果 abs (Hai-Hbi) < T，则令 Xi=O ； 如果 abs (Ha2-Hb2) >T，则令 yi=l ； 如果 abs (Ha2-Hb2) < T，则令 Yi=O ； 其中，abs()代表取绝对值，T代表阈值；步骤⑶生成零水印信息把置乱后的水印信息二进制比特流w和χ进行异或操作得 到零水印信息I1, 二进制比特流w和y进行异或操作得到零水印信息I2 ；步骤(4)保存零水印信息把生成的零水印信息I1和I2保存到水印信息库中，以备检 测时需要；所述的水印信息检测是指将待检测栅格地图中提取的零水印信息与水印信息库中的 水印信息操作得到最后水印信息，具体步骤如下步骤(5)根据待检测栅格地图数据，取出水印信息库中的水印信息IJPI2，并生成两个 长度为 7 的二进制比特流 α=(αι，α2，···，Qi,..., Ci1)禾口 ￠=(3^ β2，···，^i,..., β D，用来保存特征数据，其中α严1=0，0彡iO;；步骤(6)对待检测栅格地图数据进行分块处理根据水印信息的长度，把栅格地图数 据分成每块大小为MXM的小块；然后进行以下步骤步骤d、对每个MXM的小块数据进行二维sym2小波分解，分解级数为2 ；步骤e、对第i块小波系数中的LI^和Hffii系数分别进行统计取LI^中的最小值与最大值，均分成十个等间隔的区间，把LI^中的元素分到此十个区间之内，并返回每个范围 内元素的个数作为一个一维向量，记为V3，再取Hffii中的最小值与最大值，均分成十个等间 隔的区间，把Hffii中的元素分到此十个区间之内，并返回每个范围内元素的个数作为一个 一维向量，记为V4;步骤f、分别对向量V3和V4的元素进行从大到小的排序，记为V3’和V4'，获得V3'排在向 量前面的两个元素C3’、D3’，以及他们在向量V3中的位置H。3、HD3，获得V4'排在向量前面的两 个元素C4’、D4',以及他们在向量V4中的位置Hc4、Hd4，然后做如下处理 如果 abs (Hg3-Hd3) >T，则令 Qi=I ； 如果 abs (Hc3-Hd3) < T，则令 α i=0 ； 如果 abs (Hc4-Hd4) >T,则令 β 尸1 ； 如果 abs (Hc4-Hd4) < T，则令 β 尸0 ；步骤(7)得到检测信息把α和I1进行异或操作，得到W1 ’，把β和12进行异或操作， 得到W2'；步骤⑶反置乱水印信息把二进制比特流Wl’、W2'结合密钥进行Arnold反变换提取 出最终的水印图像I1 ’和I2'。
全文摘要
本发明公开了一种基于栅格地图特征的小波域零水印方法。现有的水印的嵌入过程复杂化，消耗计算时间太长，不利于实际应用。本发明包括水印信息嵌入和水印信息检测。水印信息嵌入是指提取原始栅格地图的零水印信息，保存入库。水印信息检测是指将待检测栅格地图中提取的零水印信息与水印信息库中的水印信息操作得到最后水印信息。本发明具有较好的隐蔽性，可以克服可逆数字水印中存在的安全漏洞。
文档编号G06T1/00GK102129660SQ20111007227
公开日2011年7月20日 申请日期2011年3月24日 优先权日2011年3月24日
发明者凌云, 林传峰, 王勋, 章志勇, 黄定军 申请人:浙江工商大学