技術(shù)特征：

1.一種基于柯西統(tǒng)計(jì)建模的數(shù)字圖像水印方法，其特征按照以下步驟進(jìn)行：

約定：I指宿主圖像；表示非下采樣剪切波變換第K尺度下第J個(gè)方向的子帶；w表示長(zhǎng)度為L(zhǎng)的二值水印序列；表示從最大能量子帶取出的L個(gè)高熵塊；S表示每個(gè)高熵塊中用于水印嵌入的系數(shù)集合；

a. 初始設(shè)置

獲取宿主圖像I并初始化變量；

b. 選取宿主圖像I的NSST域能量最大子帶

宿主圖像I進(jìn)行二級(jí)NSST，選取能量最大的子帶d：

；

c. 利用乘性嵌入方式嵌入水印并得到含水印圖像

c.1 構(gòu)造出兩種嵌入強(qiáng)度函數(shù)用于嵌入水印位“1”或“0”：

，；

c.2 將子帶d進(jìn)行不重疊分塊，選取L個(gè)高熵塊；

c.3 將高熵塊中所有系數(shù)劃分成兩個(gè)部分，選取一部分系數(shù)用于水印嵌入，記為，該部分系數(shù)位于高熵塊中的位置需要滿(mǎn)足條件；

c.4 根據(jù)下式修改每個(gè)高熵塊中對(duì)應(yīng)系數(shù)得到：

c.5 將所有子帶進(jìn)行逆NSST得到含水印圖像；

d. 利用離散傅里葉變換對(duì)含水印圖像進(jìn)行旋轉(zhuǎn)攻擊校正

d.1 含水印圖像進(jìn)行離散傅里葉變換得到頻譜幅值圖像，并進(jìn)行0-255歸一化后得到圖像A；

d.2 求得圖像A中兩條能量軸上遠(yuǎn)離中心的四個(gè)極大值位置，分別與圖像A中心連線得到四條線段；

d.3 求得四個(gè)線段與過(guò)中心水平和垂直直線的夾角的平均角度，對(duì)含水印圖像逆向旋轉(zhuǎn)度，可以得到校正后含水印圖像；

e. 利用柯西分布對(duì)校正后含水印圖像的NSST域系數(shù)進(jìn)行建模

e.1 通過(guò)卡丹公式分別求得兩種嵌入強(qiáng)度函數(shù)的反函數(shù)；

e.2 利用中值法估計(jì)出子帶的噪聲偏差：

；

e.3 根據(jù)下式分別計(jì)算每一個(gè)高熵塊中在嵌入水印“1”和嵌入水印“0”兩種假設(shè)下的柯西分布概率和：

；

f. 構(gòu)造局部最優(yōu)檢測(cè)器進(jìn)行水印提取

f.1 對(duì)于每個(gè)高熵塊，利用未嵌入水印部分的NSST系數(shù)通過(guò)樣本分位數(shù)方法估計(jì)柯西分布的位置參數(shù)和形狀參數(shù)；

f.2 每個(gè)高熵塊中用于嵌入水印的系數(shù)表示為，根據(jù)下式，利用局部最大勢(shì)能檢驗(yàn)構(gòu)造局部最優(yōu)檢測(cè)器，提取每個(gè)高熵塊的具體水印信息：

；

f.3 將每個(gè)高熵塊提取的水印按順序排列，得到提取的水印序列。