馮帆，王建華，張政保，王惠萍，戚紅軍，張海龍

1.白求恩軍醫(yī)學(xué)院基礎(chǔ)部，石家莊 050081

2.軍械工程學(xué)院計算機工程系，石家莊 050003

基于CSR-ICA模型的隱寫信息盲提取算法

馮帆1，2，王建華1，張政保2，王惠萍1，戚紅軍1，張海龍1

1.白求恩軍醫(yī)學(xué)院基礎(chǔ)部，石家莊 050081

2.軍械工程學(xué)院計算機工程系，石家莊 050003

1 引言

隨著網(wǎng)絡(luò)技術(shù)和多媒體技術(shù)的迅速發(fā)展，隱寫術(shù)已成為實現(xiàn)網(wǎng)絡(luò)隱蔽通信的全新技術(shù)手段。隱寫分析技術(shù)作為隱寫術(shù)的對抗技術(shù)，其目的是揭示出隱秘載體中秘密信息的存在，進而進行提取、篡改或破壞。鑒于隱秘信息提取的困難性，目前關(guān)于隱秘信息的提取技術(shù)研究較少，并且主要集中于特定算法的研究上[1-4]。ICA是一種統(tǒng)計計算技術(shù)，對于通常情況下，以大量樣本數(shù)據(jù)庫形式給出的多元觀測信號，它定義了一個問題生成模型，模型假設(shè)觀測數(shù)據(jù)變量是內(nèi)在變量的線性組合，且內(nèi)在變量和混合系統(tǒng)均是未知的，如果內(nèi)在變量滿足統(tǒng)計上相互獨立，且至多只有一個非高斯信號，則該變量可通過ICA的方法實現(xiàn)分離提取。ICA作為當(dāng)今信號處理領(lǐng)域的一個研究熱點，主要用于解決信號的分離提取、增強、降噪、分類等問題，近幾年來，研究者們開始將ICA模型技術(shù)用于水印及主動隱寫分析技術(shù)的研究[5-9]。

加性隱寫模型和乘性模型是數(shù)據(jù)隱寫嵌入研究應(yīng)用最多的模型，而乘性隱寫模型從數(shù)學(xué)角度上可轉(zhuǎn)化為加性模型。因此，針對加性隱寫模型，進行隱秘信息盲提取技術(shù)研究更具有通用性。隱寫信息盲提取實際上是從混疊信號中分離提取出隱寫信息的過程。對于圖像載體加性隱寫嵌入模型，隱秘信息可以看成載體和秘密信息的線性混疊信號。因絕大多數(shù)自然圖像，其概率密度函數(shù)為非高斯分布，且隱秘信息和載體信息一般為非高斯信號，唯一可能成為高斯信號的只有秘密信息，滿足ICA模型的約束條件。但是，對于隱寫分析者而言，往往不能夠得到滿足條件的兩個觀測信號。所以，隱寫信息盲提取屬于BSS問題中“欠定”ICA模型問題。研究表明[10]：源信號的非高斯性以及ICA學(xué)習(xí)算法本身的性能對源信號分離提取精度影響較大。而信號的非高斯性與其稀疏性又有著直接的關(guān)系：信號的稀疏性越好，非高斯越強。自然圖像基本上沒有稀疏性可談，但是通過相應(yīng)變換可以實現(xiàn)圖像稀疏性表示。輪廓波（Contourlet）是一個空間多分辨率分析工具，能夠根據(jù)實際應(yīng)用的具體情況，選擇所需要的空間和方向進行分析，其變換較小波變換具有更好的圖像稀疏表示，是二維圖像稀疏表示最有效的工具，圖像通過Contourlet稀疏性表示，能夠增強其非高斯性[11]。

鑒于上述分析，為便于解決問題，同時提高隱寫信息盲提取的精度，筆者將圖像稀疏性表示（CSR）與ICA有機結(jié)合，提出一種基于CSR-ICA的隱寫信息盲提取算法。算法能夠有效實現(xiàn)隱寫信息的盲提取，隱寫信息的平均提取精度達90%，運行效率高于典型FastICA算法，綜合性能優(yōu)于典型的Chandramouli算法。

2 ICA模型理論及分析

2.1 ICA數(shù)學(xué)模型

ICA模型分為超“欠定”（Overunderminded）ICA(m＞n)、“適定”（derminded）ICA(m=n)和“欠定”（undermined）ICA (m＜n)。關(guān)于超欠定和適定ICA問題的研究較為成熟，而對于“欠定”ICA問題研究才剛剛起步。因此本算法是在適定ICA模型基礎(chǔ)上進行研究。模型約束條件：各個源信號之間統(tǒng)計上相互獨立，若隨機變量y1,y2,…,yi,…,ym統(tǒng)計上獨立，當(dāng)且僅當(dāng)它們的聯(lián)合概率密度是其邊緣概率密度之積，見式（2）；源信號至多有一個高斯信號，對于模型y=Wz= WVx=WVAs=WA~s，如果y具有較強的非高斯型，則WA~=I，此時的分離信號則等于源信號；源信號混疊矩陣是一個常矩陣，且列滿秩。此時的ICA便能實現(xiàn)源信號的分離和提取。

ICA模型實際上是通過對“隱變量”生成過程的描述，在源信號和傳輸通道參數(shù)先驗知識未知的情況下，對源信號作一個非常弱的假設(shè)，然后根據(jù)其分布特性，確定目標(biāo)函數(shù)、選擇學(xué)習(xí)優(yōu)化算法，估計源信號混疊矩陣，進而實現(xiàn)對源信號的分離提取。假設(shè)存在n個獨立源信號s= [s1,s2,…,si,…,sn]T和m個觀測信號x=[x1,x2,…,xi,…,xm]T，且觀測信號x的n個隨機變量由信號s的m個隨機變量線性組合得到。則隨機變量xi表示為：

2.2 盲信號估計方法

常用ICA估計方法有極大化非高斯性估計算法、極小化互信息估計算法等。研究表明源信號非高斯強弱影響ICA分離精度?？紤]到算法的研究對象，在此僅對簡單而直觀的極大化非高斯性的ICA估計方法進行研究分析。在ICA算法估計中，峭度是一種經(jīng)典的非高斯度量指標(biāo)，是隨機變量的四階累計量，用于度量源信號偏離高斯性的程度。若隨機變量的均值為0，則其峭度：

若峭度kurt(x)=0，則隨機變量的分布為高斯分布；若峭度kurt(x)＜0，則隨機變量為次高斯分布；若峭度kurt(x)＞0，則隨機變量為超高斯分布，并且超高斯性越強，其峭度絕對值則越大。

3 基于CSR-ICA的隱寫信息盲提取算法

3.1 欠定ICA問題的適定化轉(zhuǎn)換

利用隱寫分析者僅有的一幅隱秘圖像z(k)，生成滿足適定ICA問題模型約束條件的另一幅觀測信號z1(k)。則圖像載體加性隱寫模型：

由上式可知，z(k)實際上是s(k)和w(k)的加權(quán)線性組合，見式（6），式（5）在β=1情況下的特例。

將隱秘載體z(k)作為新的載體，按照加性隱寫模型公式，嵌入一獨立同分布的隨機信號n(k)，得一新的含秘載體信號z1(k)，原信號的分布特性不受影響，并且z1(k)仍然是連續(xù)穩(wěn)定的非高斯信號，可通過s(k)和w(k)的加權(quán)線性組合得到，則：

3.2 源信號CSR前置處理

在混合信號輸入ICA前，首先采用CSR進行前置處理得到稀疏性較好，即非高斯性最大的輸入信號。自然圖像并不是簡單的一維分段光滑線段的簡單堆積，而是包含著內(nèi)在的幾何結(jié)構(gòu)，其突變信息通常為圖像的邊緣信息，并且位于光滑的輪廓曲線上。傳統(tǒng)的二維可分離小波是由一維小波的張量積構(gòu)成，在表示一個光滑輪廓時，其不同分辨率下的支撐空間為不同尺寸大小的正方形，對奇異曲線的逼近表現(xiàn)為奇異點的逼近。而輪廓波（Contourlet）是一種多尺度、多方向的圖像表示方法，其支撐區(qū)間具有隨尺度而長、寬比變化的“長條形結(jié)構(gòu)”，類似于線段的基結(jié)構(gòu)來逼近原始圖像。因此增強源信號的非高斯性，有利于源信號分離精度的提高。峭度可以作為信號極大化非高斯性的無量綱度量標(biāo)準(zhǔn)。Contourlet變換是圖像稀疏性表示的優(yōu)秀工具，在數(shù)字水印、隱寫分析等方面得到了成功應(yīng)用[12-15]。圖像CSR表示能夠獲得非高斯性較小波變換更強的高頻子圖。表1為不同變換方式所得到的高頻子圖峭度的大小對比。

表1 兩種變換方式的高頻子圖概率密度函數(shù)的峭度

表1中數(shù)據(jù)顯示：圖像經(jīng)Contourlet DBF三層分解所得16個方向高頻子圖的概率密度函數(shù)的峭度最大為67.47，32個方向高頻子圖的概率密度函數(shù)的峭度最大為90.09，而圖像經(jīng)DWT三級小波分解，所得HL、LH和HH三個方向9個高頻子圖的概率密度函數(shù)的峭度最大只有36.66。因此，圖像的Contourlet變換能夠得到較DWT變換稀疏性更好的高頻子圖。所以，根據(jù)Contourlet變換高頻子圖的歸一化峭度，選擇非高斯性較強的高頻子圖作為ICA模型的輸入數(shù)據(jù)，以提高ICA估計矩陣精度。

3.3 ICA學(xué)習(xí)方法確定

算法選擇簡單直觀的極大化非高斯性ICA估計方法，通常情況下，選擇峭度作為度量非高斯性的無量綱標(biāo)準(zhǔn)，但因峭度對異常值比較敏感，所以，將歸一化峭度作為代價函數(shù)的自然梯度學(xué)習(xí)法，實現(xiàn)分離矩陣的估計，以消除異常值的影響。代價函數(shù)：

3.4 隱寫信息提取

首先對ICA模型所涉及到的觀測信號x（含密圖像）、白化信號z和源信號s（載體和隱寫信息）的估計信號y（提取的秘密信息）分別進行Contourlet變換：

式中，C(.)表示Contourlet變換，A、V和W分別為源信號的混疊矩陣、觀測信號的白化矩陣和源信號分離矩陣。

因Contourlet變換是一個線性變換，源信號混疊矩陣在變換前后保持不變。所以，基于CSR-ICA模型得到的估計矩陣與時域矩陣相同，可直接用于隱寫信息的分離和提取。此外，因圖像稀疏性表示子圖與原始圖像相比，尺寸要小很多，所以，基于CSR-ICA的源信號分離提取，不僅提高分離矩陣的估計精度，而且還縮短了算法運行的時間。

4 實驗結(jié)果與分析

為與典型的Chandramouli算法性能進行對比，實驗分別采用“三觀測信號”和“兩觀測信號”算法模型實現(xiàn)隱寫信息的盲提取。

4.1 仿真實驗

實驗采用大小為256×256的lena.bmp灰度圖像作為原始載體圖像c(k)，隱秘載體圖像z(k)通過擴頻隱寫，在DCT和DWT高頻系數(shù)中，選取嵌入長度為L的隱寫信息得到最大系數(shù)。在此隱寫信息w(k)為均值為0、方差為1的隨機序列，若隨機序列的長度為L，則L分別取10、50、100、500、1 000、3 000、7 000。

（1）“兩觀測信號”。將待測圖像z(k)作為新的載體對象，采用加性隱寫嵌入模型嵌入一獨立同分布的隨機序列，生成另一觀測信號x1(k)，并分別對x(k)、x1(k)進行j=3的LP分解，lj=2j的多方向DBF分解，優(yōu)化選擇峭度較大的高頻子圖C(x(k))、C(x1(k))作為CSR-ICA模型觀測信號[C(x(k)),C(x1(k))]T，然后采用歸一化峭度自然梯度學(xué)習(xí)法，實現(xiàn)源混疊信號的分離和提取。其中C(·)為Contourlet變換。

（2）“三觀測信號”。將（1）中的x1(k)替換Chandramouli算法中原始信號的估計信號x?(k)，并與已知的非高斯隨機序列f(k)進行線性組合，生成另外兩個觀測信號，按照（1）中相同的方法生成[C(x(k)),C(x1(k)),C(x2(k))]T，實現(xiàn)隱寫信息的分離提取。

初始值隨機產(chǎn)生，算法運行20次所提取隱寫信息長度平均錯誤率見表2，DCT域隱寫嵌入長度分別為50和100的秘密信息的提取見圖1。

表2 基于SCR-ICA的信息長度平均估計錯誤率

圖1 秘密信息嵌入長度提取

此外，算法因?qū)τ^測信號進行Controulet變換，使圖像大小縮小了很多，降低了ICA學(xué)習(xí)算法的計算復(fù)雜度，提高了算法學(xué)習(xí)效率。表3是LP=3，DBF分解方向數(shù)為16的高頻子圖混疊信號與原始圖像混疊信號采用典型的固定點學(xué)習(xí)算法的耗時情況。

表3 算法耗時

4.2 結(jié)果分析

表2中的數(shù)據(jù)顯示：算法隱寫信息提取的正確率最大為95%、最小為86%，明顯高于Chandramouli文獻[5]提出的“三觀測信號”模型算法的提取正確率（最大70%、最小62%），雖然與Chandramouli文獻[6]中的“兩觀測信號”模型算法的實驗結(jié)果（最大為100%、最小為92%）相比，低了幾個百分點，但因其算法要求兩觀測信號必須是不同嵌入強度下的相同載體和秘密信息線性組合生成，限制條件過于苛刻，對于隱寫分析者來講，實現(xiàn)起來難度較大，甚至根本無法實現(xiàn)，如果隱寫分析者具備Chandramouli算法中觀測信號生成條件，可將其引入Chandramouli算法中，進一步提高算法的性能；如果不具備其條件，采用本文提出的算法，效果則會更好。

表3的數(shù)據(jù)顯示，F(xiàn)astICA算法的運行時間為0.361 6，而SCR-FastICA算法的運行時間為0.004 1，因此，采用圖像CSR對ICA模型輸入數(shù)據(jù)進行前置處理能夠提高算法運行效率。

5 結(jié)束語

本文對通用隱寫嵌入模型、隱寫信息盲提取原理以及ICA模型技術(shù)進行了綜合分析和研究，給出了隱寫信息盲提取“欠定”ICA問題“適定”化解決的思路，在此基礎(chǔ)上，將圖像Contourlet稀疏性表示與ICA模型技術(shù)進行有機的結(jié)合，提出了一種基于CSR-ICA模型實現(xiàn)隱寫信息盲提取算法。算法僅需一幅隱寫圖像，即可實現(xiàn)隱寫信息的提取，克服了典型Chandramouli算法約束條件的限制，降低了提取信息的錯誤率，同時提高了算法的運行效率，其綜合性能優(yōu)于典型Chandramouli算法。

[1]Fridrich J，Soukal D，Goljan M.Maximum likelihood estimation of secret message length embedded using±K steganography in spatial domain[C]//Proceedings of the SPIE-Security Steganography and Watermarking of Multimedia Contents VII，San Jose，CA，2005：595-606.

[2]張衛(wèi)明，李世取，劉九芬.對空域圖像LSB隱寫術(shù)的提取攻擊[J].計算機學(xué)報，2007，30（9）：1625-1630.

[3]Lou Der-Chyuan，Liu Chiang-Lung，Lin Chih-Linn.Message estimation for universal steganalysis using multi-classification support vector machine[J].Computer Standards&Interfaces，2009，31（2）：420-427.

[4]Kodovsky J，F(xiàn)ridrich J.Quantitative structural steganalysis of Jsteg[J].IEEE Transactions on Information Forensics and Security，2010，15（4）：681-693.

[5]Chandramoouli R.A mathematical framework for active stenanalysis[J].Multimedia Systems，2003，9（3）：301-311.

[6]Chandramouli R，Subbalakshmi K P.Active steganalysis of spread spectrum image steganography[C]//2003 IEEE International Symposium on Circuits and Systems（ISCAS 2003）：General&Nonlinear Circuits and Systems，2003：830-833.

[7]劉曉芹，王嘉禎，徐波.基于HMT和ICA的主動隱寫分析算法[J].武漢大學(xué)學(xué)報：信息版，2008，33（10）：1060-1066.

[8]Malik H.Blind watermark estimation attack for spread spectrum watermarking[J].Informatica，2009，33（1）：49-68.

[9]Umaamaheshvari A，Thanushkodi K.Digital image watermarking based on independent component analysis and ridgelet transform[J].International Journal of Computer Science and Network Security，2011，11（4）：14-17.

[10]史習(xí)智.盲信號處理——理論與實踐[M].上海：上海交通大學(xué)出版社，2008.

[11]焦李成，侯彪，王爽，等.圖像多尺度幾何分析理論與應(yīng)用——后小波分析理論與應(yīng)用[M].上海：上海交通大學(xué)出版社，2008.

[12]樓偶俊，王鉦旋.基于特征點模板的Contourlet域抗幾何攻擊水印算法研究[J].計算機學(xué)報，2009，32（2）：308-317.

[13]Sajedi H，Jamzad M.Using contourlet transform and cover selection for secure steganography[J].International Journal of Information Security，2010，9（5）：337-352.

[14]石慧，林闖，李明楚，等.基于Contourlet變換的抗幾何攻擊數(shù)字水印算法[J].光電子·激光，2011，22（10）：1575-1581.

[15]Sajedi H，Jamzad M.CBS：Contourlet-Based Steganalysis method[J].Journal of Signal Processing Systems for Signal，Image，and Video Technology，2010，61（3）：367-373.

FENG Fan1，2,WANG Jianhua1,ZHANG Zhengbao2,WANG Huiping1,QI Hongjun1,ZHANG Hailong1

1.Department of Basic Courses,Bethune Military Medical College,Shijiazhuang 050081,China
2.Department of Computer Engineering,Ordnance Engineering College,Shijiazhuang 050003,China

Aiming at steganography adding model,a blind extraction algorithm of embedded messages is proposed based on CSR-ICA.It attains estimated signal of original cover signals using linear combination just one digital image,which consists with constrain of ICA model.Then signals with optimal normalization kurtosis are selected and inputted to ICA model by optimization algorithm.In order to eliminate effect of outlier,normalization kurtosis is taken as objective function of separation algorithm on precision.The algorithm has higher comprehensive performance and conquers the localization of proposed algorithm by Chandramouli.Average extraction correct rate is 90%.The results of simulation experiments prove its validity further.

Contourlet Sparse Representation（CSR）;Independent ComponentAnalysis（ICA）model;steganography;blind extraction

針對加性隱寫模型，提出一種基于CSR-ICA的隱寫信息盲提取算法。算法僅需一幅隱寫圖像，在滿足ICA模型線性約束條件下得到載體信號的估計信號，通過Contourlet稀疏性表示（CSR）對模型輸入信號進行前置處理，優(yōu)化選取歸一化峭度性較大的信號作為模型輸入信號，將歸一化峭度作為分離算法學(xué)習(xí)的目標(biāo)函數(shù)，避免異常值給分離算法帶來的誤差。算法具有較好的綜合性能,并且克服了Chandramouli算法的局限性，提取正確率平均為90%。仿真實驗結(jié)果給出了算法的有效性驗證。

Contourlet稀疏性表示（CSR）；獨立成分分析（ICA）模型；隱寫；盲提取

A

TP393

10.3778/j.issn.1002-8331.1203-0628

FENG Fan,WANG Jianhua,ZHANG Zhengbao,et al.Blind extraction algorithm of embedded messages based on CSR-ICA. Computer Engineering and Applications,2013,49（24）：70-73.

河北省基金項目（No.052135，No.11213527D）。

馮帆（1971—），女，博士，主要研究方向為信息安全；王建華（1977—），女，博士，主要研究方向為圖像處理；張政保（1964—），男，博士，主要研究方向為信息安全。E-mail：fengfan_2357@163.com

2012-03-27

2012-06-14

1002-8331（2013）24-0070-04

CNKI出版日期：2012-08-08http://www.cnki.net/kcms/detail/11.2127.TP.20120808.1548.019.html