dct數字水印演算法
① DCT域圖像數字水印演算法,新手,這個水印嵌入程序看不明白,誰能幫我把每條注釋下Q282092728,謝謝了
M=256;%原圖像長度
N=32;%水印圖像長度
K=8; %圖像分塊大小
I=zeros(M,M);%創建一個零矩陣,用於存放載體圖像
J=zeros(N,N); %創建一個零矩陣,用於存放水印圖像
BLOCK=zeros(K,K);%創建一個零矩陣,用於存放圖像分塊
%顯示水印圖像
subplot(1,8,2);
J=imread('14','bmp'); %讀入水印圖像
imshow(J); %顯示水印圖像
title('水印圖像');
%顯示原圖像
subplot(1,3,2);
I=imread('11','bmp'); %讀入原始圖像
imshow(I); %顯示原始圖像
title('原始公開圖像');
%嵌入水印
tem=1;%創建變數tem,沒用
for p=1:N%水印圖像行循環
for q=1:N%水印圖像列循環
x=(p-1)*K+1; %x為載體圖像行坐標
y=(q-1)*K+1; %y為載體圖像列坐標
BLOCK=I(x:x+K-1,y:y+K-1); %BLOCK為載體圖像I的分塊,分塊大小為K*K,
%初始值為I(0:K-1,0:K-1)
BLOCK=dct2(BLOCK);%對BLOCK進行二維DCT變換,得到新的BLOCK即%DCT系數矩陣BLOCK
if J(p,q)==0%如果水印圖像的第(p,q)個像素為0
a=-1;%嵌入參數為-1
else
a=1; %若如果水印圖像的第(p,q)個像素為1嵌入參數為1
end
BLOCK(2,1)=BLOCK(2,1)*(1+a*0.01); %對載體圖像的分塊的DCT系數矩陣
%BLOCK
BLOCK=idct2(BLOCK);%對DCT系數矩陣進行反變換,得到嵌入水印後的載體
%圖像分塊BLOCK
I(x:x+K-1,y:y+K-1)=BLOCK;%用嵌入水印後的圖像分塊BLOCK代替載體圖像
%的對應分塊
End%水印圖像列循環結束
End%水印圖像行循環結束
%顯示嵌入水印後的圖像
subplot(1,3,3);
imshow(I);%顯示嵌入水印後圖像
title('嵌入水印後的圖像');
imwrite(I,'embedded.bmp','bmp');%將嵌入水印後圖像寫成bmp文件
夠詳細吧,希望對你有用。
② 基於DCT數字水印演算法的Matlab實現源代碼
讀取圖片
讀取水印圖片(二維碼)
原圖片做成灰度
水印圖做成二值圖
原圖做8*8DCT處理
選擇DC成分以外的元素(每個8*8塊里選擇1個元素)
生成偽隨機亂數
根據亂數做成密鑰
將水印圖根據密鑰與演算法嵌入到第6步驟
嵌入後將每個元素替換原來的第6步元素
做反變換
圖片保存,嵌入完畢
提取
讀取嵌入水印後的圖片
生成嵌入時的偽隨機亂數(生成亂數之前可以加一些幾何攻擊)
根據亂數做成密鑰
嵌入水印後的圖片做8*8DCT處理
選擇DC成分以外的元素(選取跟原來嵌入時相同)
根據演算法與密鑰提取水印圖
提取的水印圖與原水印圖做比較
保存提取的水印圖,提取完畢