數(shù)字水印matlab程序(編輯修改稿)
2025-08-03 14:04 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 oise0=10*rand(size(WImage2))。 WImage2=WImage2+noise0。 figure。imshow(WImage2,[ ])。title(39。adding uniform noise 添加噪音39。)。 watemark_detect(WImage2,Y1,I0,waterMark1)。case 3 WImage3=idctF1。 H=fspcial(39。gaussian高斯39。,[10,10],5)。 WImage3=imfilter(WImage3,H)。 figure。 imshow(WImage3,[ ])。 title(through filter [10,10] 濾波39。)。 watemark_detect(WImage3,Y1,I0,waterMark1)。case 4 WImage4=idctF1。 WImage4(1:128,1。128)=256。 figure。 imshow(WImage4)。 title(39。cutting part of the image 剪切39。)。 watemark_detect(WImage4,Y1,I0,waterMark1)。 case 5 WImage5=idctF1。 WImage5=im2double(WImage5)。 um=10。 dctm=dctmtx(8)。 p1=dctm。 p2=dctm.39。 imageDCT=blkproc(WImage5,[8,8],39。p1*p2*x39。,dctm,dctm.39。)。 DCTvar=im2col(imageDCT,[8,8],39。distinct39。).39。 n=size(DCTvar,1)。 DCTvar=(sum(DCTvar.*DCTvar)(sum(DCTvar)/n).^2)/n。 [dum,order]=sort(DCTvar)。 um=64um。 mask=ones(8,8)。 mask(order(1:um))=zeros(1,um)。 im88=zeros(9,9)。 im88(1:8,1:8)=mask。 im128128=kron(im88(1:8,1:8),ones(16))。 dctm=dctmtx(8)。 p1=dctm.39。 p2=mask(1。8,1:8)。 p3=dctm。Wimage5=bikproc(imageDCT,[8,8],39。p1*()*p339。,dctm.39。,mask(1:8,1:8),dctm)。 figure。 imshow(Wimage5)。 title(39。JPEG Image 壓縮39。)。 watemark_detect(WImage5,Y1,I0,waterMark1)。case 6 WImage6=idctF1。 WImage6=imrotate(WImage6,45,39。bilinear39。,39。corp39。)。 figure。 imshow(Wimage6)。 title(39。rotate 45 旋轉(zhuǎn)39。)。 watemark_detect(WImage6,Y1,I0,waterMark1)。case 0 break。otherwise error(39。you have a valid value(您的輸入錯(cuò)誤)39。)。endd=input(39。please input you choice(請(qǐng)輸入您的選擇):39。)。end%結(jié)束f=imread(39。39。)。%將含水印圖像f歸一化,以便于攻擊處理。m=max(max(I_W)。I_W=double(I_W )./double(m)。%攻擊attack=0。switch attack case 0, attackf= I_W。 att=39。未攻擊39。 case 1,%%1. JPEG 壓縮imwrite(I_W,39。39。,39。jpg39。,39。quality39。,30)。attackf=imread(39。39。)。attackf=double(attackf)/255。att=39。JPEG壓縮39。 case 2,% %2. 高斯低通濾波h=fspecial(39。gaussian39。,3,1)。attackf=filter2(h, I_W)。att=39。高斯低通濾波39。 case 3,%%3. 直方圖均衡化attackf=histeq(I_W)。att=39。直方圖均衡化39。 case 4,%%4. 圖像增亮attackf=imadjust(I_W,[],[,1])。att=39。圖像增亮39。 case 5,%%5. 圖像變暗attackf=imadjust(I_W,[],[0,])。att=39。圖像變暗39。 case 6,%%6. 增加對(duì)比度attackf=imadjust(I_W,[,],[])。att=39。增加對(duì)比度39。 case 7,%%7. 降低對(duì)比度attackf=imadjust(I_W,[],[,])。att=39。降低對(duì)比度39。 case 8,%%8. 添加高斯噪聲attackf=imnoise(I_W,39。gaussian39。,0,)。att=39。添加高斯噪聲39。 case 9,%%9. 椒鹽噪聲attackf=imnoise(I_W,39。salt amp。 pepper39。,)。att=39。椒鹽噪聲39。 case 10,%%10. 添加乘積性噪聲attackf=imnoise(I_W,39。speckle39。,)。att=39。添加乘積性噪聲39。end。%攻擊后處理I_W=attackf.*double(m)。figure(2)。imshow(uint8(I_W)。%顯示水印嵌入圖攻擊后效果title(att)。imwrite(uint8(I_W),39。39。)。%提取水印% clear。a=imread(39。39。)。t=sdwt_ex(double(a),39。db239。,tkey)。%根據(jù)密鑰樹(shù)分解[w,map]=extract(t,tkey)。%抽取水印[r,c]=size(w)。figure(3)。for i=1:r subplot(ceil(r/3),3,i) imshow(255100*abs(uint8(reshape(w(i,:),map(1),map(2)))))。 title(strcat(39。抽取水印圖39。,num2str(i)))。end。V=imread(39。D:\lena\39。)。Alfa=。LENGTH=2500。 [i u]=size(V)。 %計(jì)算V的規(guī)格r=100。 %設(shè)置分解矩陣的秩W=rand(i,r) %初始化WH，為非負(fù)數(shù)H=rand(r,u)maviter=100。 %最大迭代次數(shù)for iter=1:maviter W=W.*((V./(W*H))*H39。)。 %注意這里的三個(gè)公式和文中的是對(duì)應(yīng)的 W=W./(ones(i,1)*sum(W))。 H=H.*(W39。*(V./(W*H)))。end[Y1,10]=sort(H)L=imread(39。D:\lena\39。)。[Y0,10]=sort(L)x=m*n。k=LENGTH。M=zeros(x,1)。%%修改幅值較大的n個(gè)頻域成分的幅值，嵌入水?。ㄒ?yàn)閮蓚€(gè)問(wèn)題不同，所以有兩個(gè)注釋符）for i=1:x if k=1 M(x)=Y1(x)*(1+alfa*Y0(k))。 k=k1。 else M(x)=Y1(x)。 end x=x1。endN=zeros(x,1)。x=m*n。for i=1:x N(I1(i))=M(i)。enda=1。for j=1:nfor i=1:m N(i,j)=N(a)。 a=a+1。endend%%重構(gòu)圖像，得到嵌入水印的圖象B=W.*H。figure(1)。i
語(yǔ)文相關(guān)推薦
鄂ICP備17016276號(hào)-1