符浩軍，劉靜禎，楊成松，張海勃

（1.西安測繪總站，陜西 西安 710054；2.地理信息工程國家重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710054；3.解放軍理工大學(xué) 野戰(zhàn)工程學(xué)院，江蘇 南京 210007）

3S 技術(shù)應(yīng)用

基于精度特征的遙感影像數(shù)據(jù)盲水印算法

符浩軍1,2，劉靜禎1，楊成松3，張海勃1

（1.西安測繪總站，陜西 西安 710054；2.地理信息工程國家重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710054；3.解放軍理工大學(xué) 野戰(zhàn)工程學(xué)院，江蘇 南京 210007）

提出一種基于數(shù)據(jù)精度特征的遙感影像盲水印算法。首先，分析了遙感影像數(shù)據(jù)特征及其對水印算法產(chǎn)生的影響，然后構(gòu)建了水印信息位、精度位集、非精度位集三者之間穩(wěn)定的同步關(guān)系，再依據(jù)待嵌入數(shù)據(jù)點(diǎn)的精度位集生成水印信息位，最后采用量化的思想將水印信息嵌入對應(yīng)的非精度位集。實(shí)驗(yàn)表明，所提出的水印算法不僅保持了遙感影像的數(shù)據(jù)精度，而且具有較好的效率，尤為重要的是，算法對于任意強(qiáng)度的數(shù)據(jù)裁剪、放大、縮小、幾何旋轉(zhuǎn)、剪切等水印攻擊都表現(xiàn)了良好的魯棒性，為遙感影像數(shù)據(jù)提供了較好的版權(quán)保護(hù)。

遙感影像；盲水??；精度位集；非精度位集；魯棒性；版權(quán)保護(hù)

隨著遙感影像數(shù)據(jù)應(yīng)用需求的增加，其版權(quán)保護(hù)問題受到眾多研究者的關(guān)注，由于數(shù)字水印技術(shù)在該領(lǐng)域有著巨大的潛力[1-2]，針對遙感影像數(shù)字水印算法的研究熱度也在日益增加。參考文獻(xiàn)[3]～[5]分別利用余弦變換、小波變換提出了相應(yīng)的水印方案，不僅保證了水印信息的不可見性，而且對常規(guī)水印攻擊具有較好的魯棒性；李麗麗[6-7]等相繼提出了基于DFT和水印分割、基于Contourlet變換的遙感影像魯棒水印算法，但算法對幾何攻擊的魯棒性較弱；陳晨[8]等提出了一種基于內(nèi)容的自適應(yīng)傅立葉變換域水印算法，但算法并未很好地考慮遙感影像精度特性；陳輝[9]等指出普通圖像僅僅是“Picture”，以滿足人類視覺的需要，而遙感圖像除了是“Picture”之外，還包含有“Data”；范承嘯[10]等結(jié)合遙感影像數(shù)據(jù)精度特點(diǎn)，提出了一種精度可控的數(shù)字水印算法，算法很好地考慮了遙感影像的精度特點(diǎn)，且對常規(guī)攻擊有較好的抵抗力，但對幾何攻擊的魯棒性不強(qiáng)；任娜[11-12]等分別采用數(shù)學(xué)映射思想和模板匹配思想，提出了基于映射機(jī)制的遙感影像水印算法和抗幾何變換的半盲水印算法，兩種方案對幾何攻擊都有較好的抵抗能力，但對旋轉(zhuǎn)、縮放等水印攻擊抵抗較弱，且半盲水印算法的實(shí)用性不夠理想；在此基礎(chǔ)上，任娜[13]等提出了一種改進(jìn)方案，雖然解決了水印抗亮度調(diào)整攻擊的問題，但仍沒有從本質(zhì)上解決原型算法[11]存在的問題。

由上述研究成果可知，遙感影像水印算法的研究取得了較為豐碩的成果，且在水印誤差控制、算法魯棒性等方面都有涉及和考慮。但是現(xiàn)有水印算法在數(shù)據(jù)特征、算法魯棒性、算法實(shí)用性等方面的考慮還不夠充分，沒有結(jié)合遙感影像數(shù)據(jù)本身特征進(jìn)行設(shè)計(jì)，導(dǎo)致許多算法在遙感影像數(shù)據(jù)版權(quán)保護(hù)應(yīng)用方面的不適用。本文針對上述問題，結(jié)合遙感影像數(shù)據(jù)特征及其水印特性，提出了一種基于精度特征的數(shù)字水印算法。

1 遙感影像數(shù)據(jù)水印特點(diǎn)分析

與普通圖像相比，遙感影像具有其自身的數(shù)據(jù)特征，這些特征對數(shù)字水印算法的設(shè)計(jì)提出了更高的要求：

1）數(shù)據(jù)精度特征。遙感影像的數(shù)據(jù)精度直接影響數(shù)據(jù)應(yīng)用（如空間分析）的成敗。因此，遙感影像數(shù)字水印技術(shù)需要結(jié)合其數(shù)據(jù)精度特點(diǎn)進(jìn)行設(shè)計(jì)，才能滿足其后續(xù)應(yīng)用的要求。

2）數(shù)據(jù)海量特征。與普通圖像相比，遙感影像具有海量的數(shù)據(jù)特性，這就對水印算法提出了嚴(yán)格的要求：所設(shè)計(jì)的水印算法在滿足魯棒性的前提下，還需考慮算法的時效性，所設(shè)計(jì)的水印嵌入規(guī)則復(fù)雜度要低，以滿足實(shí)用性要求。

3）數(shù)據(jù)信息量的多變性。通常遙感影像數(shù)據(jù)的尺寸較大，所包含的豐富地理信息，較多的數(shù)據(jù)信息足以支撐水印信息的嵌入，但某些滿足特定應(yīng)用的遙感影像，如瓦片地圖等，其所包含的數(shù)據(jù)信息較少，如何在較少的數(shù)據(jù)信息中嵌入足夠多的水印信息，也是算法設(shè)計(jì)者需要考慮的問題。

遙感影像數(shù)據(jù)與位平面相對應(yīng)，將其精度特征反映到位平面，可將位平面分為精度位集和非精度位集。為保證遙感影像數(shù)據(jù)的精度特征，其精度位集是不可變的，而非精度位集的變化并不會給數(shù)據(jù)精度造成任何影響。假設(shè)第2位是精度位與非精度位的臨界位，則位平面中的3～7位屬于精度位集，為保持遙感影像數(shù)據(jù)的精度特性，精度位集在任何數(shù)據(jù)處理過程中均需保持一致，因此，可利用精度位集生成水印信息；位平面中的0～2位屬于非精度位集，可將水印信息嵌入其中，并不會影響遙感影像數(shù)據(jù)的精度特征。

2 基于精度特征的數(shù)字水印算法

算法的基本思路為每個遙感影像數(shù)據(jù)點(diǎn)都對應(yīng)一個精度位集和一個非精度位集，依據(jù)精度位集生成單個水印信息，然后將其嵌入到對應(yīng)的非精度位集。通過上述水印生成/嵌入方式，可以使在水印攻擊過程中，數(shù)據(jù)點(diǎn)的水印信息位、精度位集、非精度位集保持同步，有效解決了幾何攻擊下水印信息位與數(shù)據(jù)載體同步難的問題；同時，由于單個數(shù)據(jù)點(diǎn)對應(yīng)著單個水印信息位，而這些水印信息位將組合成一個完整的水印信息集，即使在數(shù)據(jù)信息很少的情況下，也能生成足夠多的水印信息并嵌入到數(shù)據(jù)載體中，從而保證了水印信息嵌入的有效性；進(jìn)一步地，由于水印嵌入/檢測都是在空間域完成的，計(jì)算復(fù)雜度低，有效保證了算法的效率；最重要的是，水印算法是基于數(shù)據(jù)精度特征來設(shè)計(jì)和考慮的，因此算法完全滿足了遙感影像數(shù)據(jù)的精度要求。

2.1 水印信息生成

根據(jù)每個遙感影像數(shù)據(jù)點(diǎn)的屬性值獲取對應(yīng)的位平面，設(shè)第i位是精度位和非精度位的臨界位，其中i∈[0,1,…,6,7]，那么其對應(yīng)的精度位集為Bp∈[i+1, i+2,…,7]，若將精度位集對應(yīng)的位平面取出，形成一 個新數(shù)據(jù)，其對應(yīng)的值域?yàn)镠p∈[0,27-i-1]；與之類似，對應(yīng)的非精度位集為Bu∈[0,1,…,i]，其對應(yīng)的值域?yàn)镠u∈[0,2i+1-1]。水印信息的生成由Hp決定，本文采用偽隨機(jī)m序列為水印信息，生成方法見參考文獻(xiàn) [14]，不再詳述。需要強(qiáng)調(diào)的是w取值為[-1,1]，且在水印信息生成過程中，為有效降低水印檢測的誤判率，需進(jìn)行如下控制：

1）對生成的水印信息集進(jìn)行取值分布控制，使其取值概率滿足式（1）要求：

2）每個遙感影像數(shù)據(jù)點(diǎn)都有一個對應(yīng)的水印信息位，眾多水印信息位組成一個完整的水印信息集，以保證能生成足夠多的水印信息，且所有水印信息能有效嵌入數(shù)據(jù)載體中。

3）有效降低精度位集所生成的水印與原始非精度位集所提取水印之間的聯(lián)系，保證二者的獨(dú)立性。

2.2 水印信息嵌入

對于每個遙感影像數(shù)據(jù)點(diǎn)，采用量化的方式將生成的水印信息嵌入對應(yīng)的非精度位集，水印信息嵌入方式為：

當(dāng)w=+1時，水印嵌入方式為：

當(dāng)w=-1時，水印嵌入方式為：

2.3 水印信息檢測

水印信息檢測流程與水印嵌入流程基本相逆，主要步驟為：

1）對于每個待檢測的遙感影像數(shù)據(jù)點(diǎn)，根據(jù)其精度特征，獲取其對應(yīng)的精度位集和非精度位集。

2）根據(jù)精度位集生成原始水印信息位w，再按照式（4）從非精度位集中提取相應(yīng)的水印信息位w'。

3）遍歷待檢測遙感影像的每個數(shù)據(jù)點(diǎn)，對每個數(shù)據(jù)點(diǎn)進(jìn)行上述處理，由此得到原始水印信息位集W和提取的水印信息位集W '。

4）計(jì)算W和W '的相關(guān)檢測系數(shù)，將其與水印閾值比較，由此判定待檢測遙感影像數(shù)據(jù)的版權(quán)歸屬，這里相關(guān)檢測系數(shù)NC的計(jì)算公式為[15]：

式中，bi=1-XOR(Wi,Wi')；i取值范圍為1≤i≤L；L為水印長度。

3 實(shí)驗(yàn)與分析

對所提出的水印算法進(jìn)行相關(guān)的實(shí)驗(yàn)和分析，所采用的實(shí)驗(yàn)數(shù)據(jù)為兩幅不同尺寸的遙感影像地圖，圖1中左幅尺寸為256×256，右幅尺寸為7 340×5 746。實(shí)驗(yàn)內(nèi)容包括可視化分析和魯棒性分析等。在實(shí)驗(yàn)過程中，將第2位作為精度位和非精度位的臨界位。

3.1 可視化分析

依據(jù)所提出的水印算法，對兩幅實(shí)驗(yàn)數(shù)據(jù)嵌入水印信息，得到含水印信息的遙感影像地圖，具體如圖2所示。

從主觀視覺上對比原始遙感影像地圖和含水印遙感影像地圖可知，水印嵌入前后沒有對遙感影像地圖造成視覺上的影響，因此水印算法具有較好的隱蔽性；客觀上，計(jì)算原始遙感影像地圖和含水印遙感影像地圖的峰值信噪比，計(jì)算方法如式（6）所示，得到小尺寸數(shù)據(jù)的PSNR=48.352 4、大尺寸數(shù)據(jù)的PSNR=49.795 3，也驗(yàn)證了水印算法的不可感知性。

圖1 實(shí)驗(yàn)數(shù)據(jù)

圖2 含水印實(shí)驗(yàn)數(shù)據(jù)

對原始遙感影像地圖進(jìn)行水印檢測，相關(guān)檢測系數(shù)分別為0.001 6、0.001 3。因此可知，由精度位集所生成的水印與原始非精度位集所提取水印之間關(guān)聯(lián)性很低，基本可忽略不計(jì)，所設(shè)計(jì)的算法有效地降低了水印的誤判率。對含水印遙感影像地圖進(jìn)行水印檢測，其相關(guān)檢測系數(shù)均為1.00，說明算法對遙感影像地圖具有很好的適應(yīng)性。同時實(shí)驗(yàn)中假定第2位是精度位和非精度位的臨界位，那么水印嵌入所引起的數(shù)據(jù)誤差被嚴(yán)格控制在3個灰度級以內(nèi)，因此算法能較好地控制水印嵌入所引起的數(shù)據(jù)誤差，保持了遙感影像地圖的數(shù)據(jù)精度特征。

3.2 魯棒性分析

為了驗(yàn)證水印算法對任意水印攻擊方式在任意攻擊強(qiáng)度下的抵抗能力，在Photoshop平臺下對含水印遙感影像地圖進(jìn)行了一系列不同強(qiáng)度不同方式的水印攻擊，本文選用數(shù)據(jù)裁剪、放大、縮小、剪切、幾何旋轉(zhuǎn)、亮度調(diào)整、飽和度調(diào)整等水印攻擊方式，檢測攻擊后的數(shù)據(jù)相應(yīng)的水印信息，表1、2給出了兩幅實(shí)驗(yàn)數(shù)據(jù)在各種攻擊后的相關(guān)實(shí)驗(yàn)結(jié)果。

分析表1、2可知，對于數(shù)據(jù)裁剪、放大、縮小、剪切、幾何旋轉(zhuǎn)、亮度調(diào)整等水印攻擊方式，水印算法具有很好的抵抗能力，即使在水印攻擊強(qiáng)度很大的情況下，算法的魯棒性也很穩(wěn)定；尤其是針對裁剪、剪切等攻擊，即使在含水印數(shù)據(jù)只有幾百個點(diǎn)的情況下，水印檢測系數(shù)仍保持為1.00，說明水印算法對信息含量小的小尺寸數(shù)據(jù)仍然有較強(qiáng)的適應(yīng)性。

表1 小尺寸數(shù)據(jù)在各種攻擊實(shí)驗(yàn)下的檢測結(jié)果

表2 大尺寸數(shù)據(jù)在各種攻擊實(shí)驗(yàn)下的檢測結(jié)果

表3 算法魯棒性的橫向比較

將本文算法與現(xiàn)有成果進(jìn)行橫向比較，由表3可知，本文算法在抗幾何攻擊的魯棒性方面要優(yōu)于其他算法，尤其是對旋轉(zhuǎn)、縮放等攻擊的優(yōu)勢較為明顯；且本文算法是基于空間域進(jìn)行設(shè)計(jì)的，而參考文獻(xiàn)[10]是通過變換域進(jìn)行水印嵌入/檢測的，因此本文算法效率占優(yōu)；參考文獻(xiàn)[11]屬于半盲水印算法，本文算法屬于盲水印算法，因此本文算法的實(shí)用性較優(yōu)。綜上所述，本文所提出的水印算法具有良好的魯棒性和實(shí)用性。

4 結(jié) 語

根據(jù)遙感影像數(shù)據(jù)的精度特征，本文提出了一種適合于遙感影像數(shù)據(jù)版權(quán)保護(hù)的水印算法，并從可視化、魯棒性等方面對算法進(jìn)行了實(shí)驗(yàn)與分析。結(jié)果表明，算法不僅對數(shù)據(jù)放大、縮小、幾何旋轉(zhuǎn)、亮度調(diào)整、裁剪等水印攻擊具有較好的魯棒性，而且對于高強(qiáng)度數(shù)據(jù)剪切處理，也具有很好的抵抗能力，有效解決了現(xiàn)有水印算法無法適應(yīng)小數(shù)據(jù)水印嵌入/檢測的難題，算法尤為適合于小尺寸遙感影像數(shù)據(jù)（如網(wǎng)絡(luò)瓦片地圖等）的版權(quán)保護(hù)。另一方面，由于本文算法所采用的水印信息是基于數(shù)據(jù)自身特征生成的，因此具有較好的可移植性，可推廣應(yīng)用到柵格地圖版權(quán)保護(hù)等領(lǐng)域，但本算法并不適用于多版權(quán)水印嵌入的情況，下一步可結(jié)合多版權(quán)需求，研究基于內(nèi)容的影像數(shù)據(jù)自適應(yīng)分塊策略，在此基礎(chǔ)上對本算法進(jìn)行改進(jìn)，從而解決本算法不滿足多版權(quán)應(yīng)用要求的問題。

[1]孫圣和,陸哲明,牛夏牧,等.數(shù)字水印技術(shù)及應(yīng)用[M].北京:科學(xué)出版社,2004:10-14

[2]楊義先,鈕心忻.數(shù)字水印理論和技術(shù)[M]. 北京:高等教育出版社,2006:20-21

[3]耿迅,龔志輝,張春美.基于HVS和整數(shù)小波變換的遙感圖像水印算法[J].測繪通報,2007(8):20-22

[4]朱長青,任娜.一種基于偽隨機(jī)序列和DCT的遙感影像水印算法[J].武漢大學(xué)學(xué)報(信息科學(xué)版),2011,36(12):1 427-1 429

[5]王向陽,楊紅穎,鄔俊.基于內(nèi)容的離散余弦變換域自適應(yīng)遙感圖像數(shù)字水印算法[J].測繪學(xué)報,2005,34(4):324-330

[6]李麗麗,孫勁光.基于DFT和水印分割的遙感影像數(shù)字水印方案[J].計(jì)算機(jī)系統(tǒng)應(yīng)用,2011,20(9):204-206

[7]李麗麗,孫勁光.基于Contourlet變換的遙感影像數(shù)字水印方案[J].計(jì)算機(jī)應(yīng)用與軟件,2012,29(2):60-63

[8]陳晨,何建農(nóng).DFT域數(shù)字水印算法在遙感圖像中的應(yīng)用[J].計(jì)算機(jī)與現(xiàn)代化,2008(11):77-79

[9]陳輝,郭科,鄭文峰.數(shù)字水印技術(shù)應(yīng)用于遙感圖像版權(quán)保護(hù)的評測標(biāo)準(zhǔn)研究[J].物探化探計(jì)算技術(shù),2008,30(5):436-441

[10]范承嘯,符浩軍.一種精度可控的柵格地理數(shù)據(jù)數(shù)字水印算法[J].地理空間信息,2014,12(2):18-20

[11]任娜,朱長青,王志偉.抗幾何攻擊的高分辨率遙感影像半盲水印算法[J].武漢大學(xué)學(xué)報(信息科學(xué)版),2011,36(3):329-332

[12]任娜,朱長青,王志偉.基于映射機(jī)制的遙感影像盲水印算法[J].測繪學(xué)報,2011,40(5):623-627

[13]任娜,朱長青,佟德宇.改進(jìn)的基于遙感影像線性亮度調(diào)整的水印檢測算法[J].南京師范大學(xué)學(xué)報(工程技術(shù)版),2014(3):52-56

[14]楊成松.矢量地理數(shù)據(jù)數(shù)字水印模型與算法研究[D].鄭州:信息工程大學(xué),2011:43-44

[15]符浩軍,朱長青,李永練,等.抗幾何變換的二值柵格地圖盲水印算法[J].測繪科學(xué)技術(shù)學(xué)報,2013,30(1):91-94

P237

B

1672-4623（2017）07-0045-04

10.3969/j.issn.1672-4623.2017.07.014

符浩軍，博士，工程師，主要研究方向?yàn)榈乩頂?shù)據(jù)共享與安全、多源地名數(shù)據(jù)獲取與融合等。

2016-05-03。

項(xiàng)目來源：地理信息工程國家重點(diǎn)實(shí)驗(yàn)室開放研究基金資助項(xiàng)目（SKLGIE2015-M-4-5）；中國博士后科學(xué)基金資助項(xiàng)目（2016M592925）；國家自然科學(xué)基金青年基金資助項(xiàng)目（41401518）。