圖像數(shù)字水印的方案(doc畢業(yè)設(shè)計(jì)論文)(編輯修改稿)
2025-06-06 00:30 本頁面
【文章內(nèi)容簡介】 結(jié)果：④ 水印全過程（經(jīng)空域壓縮檢測水印）程序源代碼%水印加入程序Q=input(39。請輸入放縮因子的值(建議小于1)：Q＝39。)P=input(39。請輸入您所希望的圖像放縮系數(shù)值():P=39。)subplot(3,3,1)RGB=imread(39。浙江臺州學(xué)院39。,39。jpg39。)。imshow(RGB)。title(39。原始圖像39。)。subplot(3,3,2)imshow(RGB(:,:,3))。title(39。B分量39。)。subplot(3,3,3)N=dct2(RGB(:,:,3))。imshow(log(abs(N)),[]),colorbar。title(39。B分量能量分布39。)。subplot(3,3,4)I=imread(39。lena139。,39。bmp39。)。imshow(I)。title(39?；叶人D像39。)。subplot(3,3,5)M=dct2(I)。imshow(log(abs(M)),[]),colorbar。title(39。水印能量分布39。)。subplot(3,3,7)J=M(1:128,1:128)。J(128:364,128:400)=0。J=rot90(J)。J=rot90(J)。J(365:600,401:750)=0。J=rot90(J)。J=rot90(J)。N=N+Q*J。K=idct2(N)。RGB(:,:,3)=K。imshow(RGB)。title(39。加入水印后圖像39。)。subplot(3,3,8)I=imresize(RGB,P,39。nearest39。)。imshow(I)。title(39。壓縮P倍圖像39。)。subplot(3,3,9)J=imresize(I,1/P,39。nearest39。)。imshow(J)。title(39。再放大P倍還原圖像39。)%水印提取程序subplot(3,3,6)RGB1=imread(39。浙江臺州學(xué)院39。,39。jpg39。)。N=dct2(J(:,:,3))。M=dct2(RGB1(:,:,3))。M=(NM)/Q。B=idct2(M(236:365,350:401))。Y=mat2gray(B)。imshow(Y)。title(39。經(jīng)放縮后提取的水印圖像39。)：運(yùn)行結(jié)果： 圖像水印的dwt算法%以下是水印提取算法clear all。clc。 %保存時(shí)間start_time=cputime。figure(1)。 %讀出原始圖像subplot(1,2,1)。input=imread(39。39。)。imshow(input)。title(39。原始圖像39。)。 %讀出水印圖像subplot(1,2,2)。watermarked_image=imread(39。39。)。imshow(watermarked_image,[])。title(39。水印圖像39。)。 %三色分離input=double(input)。watermarked_image=double(watermarked_image)。inputr=input(:,:,1)。watermarked_imager=watermarked_image(:,:,1)。inputg=input(:,:,2)。watermarked_imageg=watermarked_image(:,:,2)。inputb=input(:,:,3)。watermarked_imageb=watermarked_image(:,:,3)。 %水印圖像R的分解[Cwr,Swr]=WAVEDEC2(watermarked_imager,2,39。haar39。)。 %圖像R的分解[Cr,Sr]=WAVEDEC2(inputr,2,39。haar39。)。 %水印圖像G的分解[Cwg,Swg]=WAVEDEC2(watermarked_imageg,2,39。haar39。)。 %圖像R的分解[Cg,Sg]=WAVEDEC2(inputg,2,39。haar39。)。 %水印圖像B的分解[Cwb,Swb]=WAVEDEC2(watermarked_imageb,2,39。haar39。)。 %圖像B的分解[Cb,Sb]=WAVEDEC2(inputb,2,39。haar39。)。 %提取水印小波系數(shù)%提取水印R的小波系數(shù)r=。for k=0:3 whr(k+1,:)=Cwr(1+size(Cwr,2)/4+k*size(Cwr,2)/16:... size(Cwr,2)/4+(k+1)*size(Cwr,2)/16)... Cr(1+size(Cr,2)/4+k*size(Cr,2)/16:... size(Cr,2)/4+(k+1)*size(Cr,2)/16)。 wvr(k+1,:)=Cwr(1+size(Cwr,2)/2+k*size(Cwr,2)/16:... size(Cwr,2)/2+(k+1)*size(Cwr,2)/16)... Cr(1+size(Cr,2)/2+k*size(Cr,2)/16:... size(Cr,2)/2+(k+1)*size(Cr,2)/16)。 wdr(k+1,:)=Cwr(1+3*size(Cwr,2)/4+k*size(Cwr,2)/16:... 3*size(Cwr,2)/4+(k+1)*size(Cwr,2)/16)... Cr(1+3*size(Cr,2)/4+k*size(Cr,2)/16:... 3*size(Cr,2)/4+(k+1)*size(Cr,2)/16)。endwhr=(whr(1,:)+whr(2,:)+whr(3,:)+whr(4,:))/(4*r)。wvr=(wvr(1,:)+wvr(2,:)+wvr(3,:)+wvr(4,:))/(4*r)。wdr=(wdr(1,:)+wdr(2,:)+wdr(3,:)+wdr(4,:))/(4*r)。war=(Cwr(1:size(Cwr,2)/16)Cr(1:size(Cr,2)/16))/r。 %提取水印小波系數(shù)%提取水印G的小波系數(shù)g=。for k=0:3 whg(k+1,:)=Cwg(1+size(Cwg,2)/4+k*size(Cwg,2)/16:... size(Cwg,2)/4+(k+1)*size(Cwg,2)/16)... Cg(1+size(Cg,2)/4+k*size(Cg,2)/16:... size(Cg,2)/4+(k+1)*size(Cg,2)/16)。 wvg(k+1,:)=Cwg(1+size(Cwg,2)/2+k*size(Cwg,2)/16:... size(Cwg,2)/2+(k+1)*size(Cwg,2)/16)... Cg(1+size(Cg,2)/2+k*size(Cg,2)/16:... size(Cg,2)/2+(k+1)*size(Cg,2)/16)。 wdg(k+1,:)=Cwg(1+3*size(Cwg,2)/4+k*size(Cwg,2)/16:...
公司管理相關(guān)推薦
鄂ICP備17016276號-1