數(shù)字水印matlab程序-文庫吧資料
2025-07-13 14:04本頁面
【正文】 體圖像39。 psnr_cover=double(a0)。39。 k2=randn(1,8)。 %%%%%%%%%%%%%%%5 以下生成水印信息 %% mark=I。 rm=dimI(1)。水印圖像39。 subplot(2,3,1)。 I=ceil(I)。 I=rgb2gray(I)。39。 start_time=cputime。endclear all。嵌入水印后的圖象39。imshow(B)。endend%%重構(gòu)圖像，得到嵌入水印的圖象B=W.*H。for j=1:nfor i=1:m N(i,j)=N(a)。for i=1:x N(I1(i))=M(i)。endN=zeros(x,1)。 else M(x)=Y1(x)。%%修改幅值較大的n個(gè)頻域成分的幅值，嵌入水?。ㄒ?yàn)閮蓚€(gè)問題不同，所以有兩個(gè)注釋符）for i=1:x if k=1 M(x)=Y1(x)*(1+alfa*Y0(k))。k=LENGTH。)。end[Y1,10]=sort(H)L=imread(39。 H=H.*(W39。)。 %設(shè)置分解矩陣的秩W=rand(i,r) %初始化WH，為非負(fù)數(shù)H=rand(r,u)maviter=100。 [i u]=size(V)。Alfa=。D:\lena\39。end。抽取水印圖39。 subplot(ceil(r/3),3,i)figure(3)。%根據(jù)密鑰樹分解[w,map]=extract(t,tkey)。db239。)。a=imread(39。)。imwrite(uint8(I_W),39。imshow(uint8(I_W)。%攻擊后處理I_W=attackf.*double(m)。添加乘積性噪聲39。,)。 case 10,%%10. 添加乘積性噪聲attackf=imnoise(I_W,39。椒鹽噪聲39。,)。salt amp。att=39。gaussian39。att=39。att=39。att=39。att=39。att=39。att=39。,3,1)。 case 2,% %2. 高斯低通濾波h=fspecial(39。JPEG壓縮39。39。quality39。jpg39。39。%%1. JPEG 壓縮 case 1,未攻擊39。 case 0,%攻擊attack=0。m=max(max(I_W)。39。)。endd=input(39。you have a valid value(您的輸入錯(cuò)誤)39。case 0 break。)。 title(39。 figure。corp39。bilinear39。case 6 WImage6=idctF1。)。 title(39。 figure。,dctm.39。Wimage5=bikproc(imageDCT,[8,8],39。8,1:8)。 p1=dctm.39。 im128128=kron(im88(1:8,1:8),ones(16))。 im88=zeros(9,9)。 mask=ones(8,8)。 [dum,order]=sort(DCTvar)。 n=size(DCTvar,1)。distinct39。)。p1*p2*x39。 p2=dctm.39。 dctm=dctmtx(8)。 WImage5=im2double(WImage5)。 watemark_detect(WImage4,Y1,I0,waterMark1)。cutting part of the image 剪切39。 imshow(WImage4)。128)=256。case 4 WImage4=idctF1。)。 imshow(WImage3,[ ])。 WImage3=imfilter(WImage3,H)。gaussian高斯39。case 3 WImage3=idctF1。)。title(39。 figure。 noise0=10*rand(size(WImage2))。while d~=0switch d case 1watermark_detect(idctF1,Y1,I0,waterMark1)。please input you choice(1,2,3,4,5,6):39。)。 %旋轉(zhuǎn)45度disp(39。510 quality JPEG pressing 壓縮39。)。%剪切disp(39。3adding filter [10 10] 濾波39。)。%添加噪音disp(39。1smoothing patterns39。)。disp(39。input you choice according to the following image processing operation:39。)。 title(39。endendendfigure(5)。if (I_W(x,y)P(x,y))0F(p,q)=1。%提取水印for p=1:size/Bfor q=1:size/Bx=(p1)*B+1。 end end w0=w1。 end if j1= =0 j1=colum。 end if j1colum j1=mod(j1,colum)。 j1=i+2*j。plot(array2)。endresultStr=char(extractedInfo)。 infochar=infochar+bitset(0,n,1)。if array2(8*(m1)+n)==1for n=1:8extractedInfo=zeros(InfoStrSize,1)。提取水印39。plot(array2)。 k=k+1。 end Alpha= Alpha1。elseend array2(l)=double(tmp2)。 tmp=tmp/Alpha。 %block_dct1(1,1)=block_dct1(1,1)*(1+Alpha*mark(k))。 T=sum(T)。 T=0。if m=1|n=1BW_8_8=BW(x:x+block1,y:y+block1)。dct是有損壓縮如jpeg使用的技術(shù)。Block_dct1=dct2(block_dct1)。 %取原始圖像小格中的像素點(diǎn)到block_dct1矩陣中。 y=(n1)*block+1。 l=1。array2=zeros(1,InfoStrSize*8)。)。I=imread(39。 )。 D=imread(39。%提取水印2：39。)%保存該圖像D=uint8(D)。title(39。endend%顯示嵌入水印后的圖像subplot(3,2,4)。 D(x:x+block1,y:y+block1)=Block_dct。 % block_dct1(1,1)=block_dct1(1,1)*(1+Alpha*mark(k))。 l=l+1。 %block_dct1(1,1)=block_dct1(1,1)*(1+Alpha*mark(k))。 T=sum(T)。 T=0。if m=1|n=1BW_8_8=BW(x:x+block1,y:y+block1)。Dct是可逆的運(yùn)算%對(duì)二維數(shù)組進(jìn)行離散余弦變換。%取原始圖像小格中的像素點(diǎn)到block_dct1矩陣中。 y=(n1)*block+1。 k=1。)。Title(39。subplot(3,2,3)。%BW=edge(I,’Sobel’)。)。)%顯示prewitt為算子的邊緣圖BW=edge(i,39。title(39。)。i=imread(39。)。title(39。 array(8*(m1)+n)=bitget(Infochar,n)。為99for n=1:8%% 39。%將字符串轉(zhuǎn)換為位數(shù)組array=zeros(1,InfoStrSize*8)。dcf39。%產(chǎn)生水印，并顯示水印信息；subplot(3,2,1)。BW
語文相關(guān)推薦
鄂ICP備17016276號(hào)-1