watermark插件的github地址：https://github.com/saucxs/watermark
一、簡單
阿里巴巴內網的不可見水印用的是什么算法？
這樣就可以加水印了。我選擇了百度的首頁，好像百度對這個有防護。
哇，嚇到我了，原來可以這么玩。
View Code
html代碼
View Code
寫了插件，這個是測試地址
http://www.chengxinsong.cn/watermark_test/
包括，測試，重置，顯示，隨機，四個部分。
特性：1、測試對水印參數屬性，重置水印屬性參數，顯示此時的水印屬性參數，隨機產生水印屬性參數；
2、水印按鈕組是position值fixed，可以浮現在頁面之上，不占字節。
3、對系統的各個部分頁面進行水印的測試。
二、復雜
頻域制定數字盲水印
信號是有頻率的，一個信號可以看做是無數個不同階的正弦信號的的疊加。
上式為傅里葉變換公式，是指時域信號（對於信號我們說時域，因為是與時間有關的，而圖像我們往往說空域，與空間有關），是指頻率。
1、編碼的目的有二，一是對水印加密，二控制水印能量的分布。以下是疊加數字盲水印的實驗。
（1）原圖像。尺寸300*240 ，漢子一枚，
（2）水印照片。
（3）水印編碼。編碼方式采用隨機序列編碼，通過編碼，水印分布到隨機分布到各個頻率，並且對水印進行了加密。
（4）原圖像頻域。經歷的是傅里葉變換，下圖變換后的頻域圖像
（5）水印圖像頻域。經歷的是傅里葉變換，下圖變換后的頻域圖像
（6）合並水印和原圖。之后，將疊加水印的頻譜進行傅里葉逆變換，得到疊加數字水印后的圖像，，將圖像頻域和水印編碼進行合並。看不出來已經加了水印吧，哈哈哈
實際上，我們是把水印以噪聲的形式添加到原圖像中。
（7）水印圖與原圖的殘差（看不出來殘差區別，需要調整對比度才能看得出來）
（8）最終的均方差（MSE）和信噪比（PSNR）
（9）下圖是原圖頻譜豎過來的樣子，其能量主要集中在低頻。
那么，為什么頻譜發生了巨大的變化，而在空域卻變化如此小呢？這是因為我們避開了圖像的主要頻率。
合並之后
（10）水印提取是水印疊加的逆過程，
(11)提取后，得到水印。
代碼
clc;clear;close all;
alpha = 1;
%% read data
im = double(imread('G:\2017學習\Work\圖片水印\test.jpg'))/255;
mark = double(imread('G:\2017學習\Work\圖片水印\watermark.png'))/255;
figure, imshow(im),title('original image');
figure, imshow(mark),title('watermark');
%% encode mark
imsize = size(im);
%random
TH=zeros(imsize(1)*0.5,imsize(2),imsize(3));
TH1 = TH;
TH1(1:size(mark,1),1:size(mark,2),:) = mark;
M=randperm(0.5*imsize(1));
N=randperm(imsize(2));
save('G:\2017學習\Work\圖片水印\encode.mat','M','N');
for i=1:imsize(1)*0.5
for j=1:imsize(2)
TH(i,j,:)=TH1(M(i),N(j),:);
end
end
% symmetric
mark_ = zeros(imsize(1),imsize(2),imsize(3));
mark_(1:imsize(1)*0.5,1:imsize(2),:)=TH;
for i=1:imsize(1)*0.5
for j=1:imsize(2)
mark_(imsize(1)+1-i,imsize(2)+1-j,:)=TH(i,j,:);
end
end
figure,imshow(mark_),title('encoded watermark');
%% add watermark
FA=fft2(im);
figure,imshow(FA);title('spectrum of original image');
FB=FA+alpha*double(mark_);
figure,imshow(FB); title('spectrum of watermarked image');
FAO=ifft2(FB);
figure,imshow(FAO); title('watermarked image');
%imwrite(uint8(FAO),'watermarked image.jpg');
RI = FAO-double(im);
figure,imshow(uint8(RI)); title('residual');
%imwrite(uint8(RI),'residual.jpg');
xl = 1:imsize(2);
yl = 1:imsize(1);
[xx,yy] = meshgrid(xl,yl);
figure, plot3(xx,yy,FA(:,:,1).^2+FA(:,:,2).^2+FA(:,:,3).^2),title('spectrum of original image');
figure, plot3(xx,yy,FB(:,:,1).^2+FB(:,:,2).^2+FB(:,:,3).^2),title('spectrum of watermarked image');
figure, plot3(xx,yy,FB(:,:,1).^2+FB(:,:,2).^2+FB(:,:,3).^2-FA(:,:,1).^2+FA(:,:,2).^2+FA(:,:,3).^2),title('spectrum of watermark');
%% extract watermark
FA2=fft2(FAO);
G=(FA2-FA)/alpha;
GG=G;
for i=1:imsize(1)*0.5
for j=1:imsize(2)
GG(M(i),N(j),:)=G(i,j,:);
end
end
for i=1:imsize(1)*0.5
for j=1:imsize(2)
GG(imsize(1)+1-i,imsize(2)+1-j,:)=GG(i,j,:);
end
end
figure,imshow(GG);title('extracted watermark');
%imwrite(uint8(GG),'extracted watermark.jpg');
%% MSE and PSNR
C=double(im);
RC=double(FAO);
MSE=0; PSNR=0;
for i=1:imsize(1)
for j=1:imsize(2)
MSE=MSE+(C(i,j)-RC(i,j)).^2;
end
end
MSE=MSE/360.^2;
PSNR=20*log10(255/sqrt(MSE));
MSE
PSNR