【正文】 )))/row*col第 35 頁 共 35 頁。err = sum(sum(xor(m139。m1=mark(:)。)。181。244。216。224。181。247。210。F:\209。)。181。244。216。224。181。247。210。F:\209。202。235。243。227。198。 fm=sum((x0x2).*(x0x2))。200。235。197。196。197。197。196。161。174。203。235。182。227。198。)。181。244。216。224。181。247。210。F:\209。)。181。244。216。224。181。247。210。F:\209。219。192。220。212。)。,39。imwrite(w3,39。blue39。color39。Fontsize39。39。207。205。211。203。181。200。204。imshow(w3)。endw3=reshape(w2,row,col)。 w23=(car3(i+2)/recordx(j,3)1)/a。 w21=(car3(i)/recordx(j,1)1)/a。200。204。211。 %5203。cdr1=detcoef(cr,lr,1)。cdr3=detcoef(cr,lr,3)。db439。)。[cr,lr]=wavedec(x2,3,39。\39。198。\210。185。207。210。187。208。[x2,fs,nbits]=wavread(39。200。204。211。 %203。\39。198。\210。185。207。210。187。208。wavwrite(s1,fs,nbits,39。%marked=[y。186。208。198。210。181。186。211。203。200。199。181。181。178。%6186。db439。]。,cd239。 endendc1=[x39。 x1(t)=x1(t)*(1+a*w(j))。 j=fix(i/k)+1。a=。)。, 39。, 16, 39。, 39。241。188。161。174。title(39。subplot(2,1,1)。)。w=mark(:)。)。181。244。216。224。181。247。210。F:\209。162。197。161。174。235。193。]。 。 。%y=[。186。208。178。205。206。lenx=length(x)。202。207。198。181。200。 %204。cd1=detcoef(c,l,1)。cd3=detcoef(c,l,3)。db439。226。214。168。161。182。253。)。 [c,l]=wavedec(x,3,39。\39。198。\210。185。207。210。187。208。[x,fs,nbits]=wavread(39。 close all。)。161。174。196。161。225。208。241。188。196。245。185。253。173。211。title(39。subplot(3,2,6)。211。203。181。200。204。202。endendend%207。if a0W(p,q)=0。BLOCK2=idct2(BLOCK2)。BLOCK2=J(x:x+K1,y:y+K1)。y=(q1)*K+1。161。174。161。225。208。241。188。196。245。185。211。)。,39。I=imread(39。39。207。205。181。186。203。209。imshow(J)。)。,39。J=imread(39。Quality39。jpeg39。39。)。,39。L=imread(39。203。JPEG209。208。189。207。205。181。186。211。203。200。199。%182。BLOCK1=zeros(K,K)。J=zeros(M,M)。K=8。M=256。202。178。187。185。211。 %203。39。211。203。181。200。204。214。207。205。211。203。186。180。imshow(W)。161。174。196。161。225。190。212。elseW(p,q)=1。a=BLOCK2(1,1)/BLOCK1(1,1)1。BLOCK1=idct2(BLOCK1)。BLOCK1=I(x:x+K1,y:y+K1)。for p=1:Nfor q=1:Nx=(p1)*K+1。bmp39。39。)。,39。I=imread(39。211。203。200。204。214。207。205。181。211。203。200。199。 %180。bmp39。39。)。241。188。196。243。161。174。235。182。title(39。subplot(3,2,3)。207。205。181。186。211。203。200。199。202。 endend%207。 BLOCK=idct2(BLOCK)。 else a=1。 BLOCK=dct2(BLOCK)。 y=(q1)*K+1。161。174。235。182。)。241。188。161。174。title(39。)。,39。J=imread(39。241。188。161。174。190。212。)。241。188。170。171。188。173。title(39。)。,39。I=imread(39。241。188。173。190。212。BLOCK = zeros(K,K)。I=zeros(M,M)。200。164。241。188。161。174。N=32。182。179。207。205。 %212。235。242。204。196。161。174。235。182。, mm139。mm1=mark1(:)。[row,col]=size(mark)。39。)。%%%%%%%%%%%%%%%%%%%%%%mark=imread(39。194。194。206。203。SNR=10*log(fm/fz)%%%%%%%%%%%%%%%%%%%%%188。%%%%fz=sum(x0.*x0)。177。212。208。181。186。208。181。211。203。193。200。199。203。%%%%%%%%%%%%%%%%%%%%%%188。)。imaged=imread(39。image1=imread(num2str(filename1))。188。198。196。 %%%%%%%%%%%%208。bmp39。39。)。161。174。196。161。225。208。241。188。161。174。172。211。title(39。colormap(gray(256))。161。174。196。161。225。190。212。 else W(p,q)=255。for p=1:N for q=1:N a=caas(p,q)cas(p,q)。db139。)。[caa,chh,cvv,cdd]=dwt2(imaged,39。db139。)。[ca,ch,cv,cd] = dwt2(image1,39。39。image1=double(image1)。161。174。161。225。)。,39。imwrite(markedimage,gray(256),39。marked image39。image(markedimage)。subplot(2,2,3)。207。205。211。203。186。200。199。202。%207。db139。) 。 endendIM= idwt2(Ca,chs,cvs,cds,39。 else a=1。)。[cas,chs,cvs,cds] = dwt2(ca,39。db139。161。174。235。182。)。title(39。subplot(2,2,2)。39。 % orginal image for watermarkingimage1=double(image1)。original image39。imshow(image1)。 image1=imread(num2str(filename1))。select the image39。*.*39。200。164。161。174。N=64。182。179。207。205。%212。 clc?；谛〔ㄗ儞Q語音端點（2）clear all。 b=1。 b=b+1。))。b=1。)。)。% Damga uretme% % imgew: Damga imgesi% imge1: Orjinal imge% function W=dmg(imgew,imge1)。db139。 % k*abs(double(c1(i,j)))*W(i,j)。 else a=。size(W)。)。[caa chh cvv cdd]=dwt2(W,39。)。)。c1 = [ch cv cd]。db139。 k=。)))/row*colfunction [D]=blockdwt2(A,W)。err = sum(sum(xor(m139。m1=mark(:)。)。mark1= imread(39。39。202。235。243。227。198。 fm=sum((x0x2).*(x0x2))。200。235。197。196。197。197。196。161。174。203。235。182。227。198。x2=imaged。39。x0=image1。219。192。220。212。)。 % shows the image adaptive watermarkimwrite(mark,gray(256),39。watermark39。image(mark)。subplot(2,2,4)。 q=1。 q=q+1。mark=zeros(size(imagew))。p=1。 end endend[w2 h2]=size(watermark)。 % image adaptive watermark%watermark=watermark*255/max(watermark(:))for i = 1:row % thresholding for j = 1:col if watermark(i,j)0 watermark(i,j) = 255。204。253。161。225。161。174。)。 % shows the watermarked imageimwrite(marked,gray(256),39。Watermarked image39。image(marked)。subplot(2,2,3)。 % generates the watermarked image%markedmax = max(marked(:))。195。247。253。175。235。182。161。174。)。title(39。 subplot(2,2,2)。39。[row,col]=size(image1)。)。title(39。subplot(2,2,1)。)。,39。 [filename1,pathname]=uigetfile(39。 close all。一般研究的是不可見水印。本文主要研究魯棒水印。魯棒水印能夠抵抗各種有意或無意的攻擊。另外，此類算法相當于視頻壓縮編碼系統(tǒng)加入的后處理功能，因此幾乎不會影響現(xiàn)有視頻壓縮編碼系統(tǒng)的正常運作，可最大限度的利用現(xiàn)有硬件資源。前者要與視頻壓縮編碼算法整合，在水印嵌入過程中，此類算法具有能夠主動利用編碼過程中產(chǎn)生的各種信息對水印信息嵌入進行實時動態(tài)調(diào)控的優(yōu)勢。基于原始域視頻水印技術大多來自較為成熟的靜止語音端點技術，按照水印嵌入過程可分為兩類方案，一種是利用空間域圖像自身特性直接在空間域內(nèi)完成水印信息的嵌入:另一種是在變換域內(nèi)進行水印信息嵌入。5視頻水印概述及算法分析 視頻水印的分類