數(shù)字水印技術(shù)LSB算法MATLAB源代碼

1、close allclear all;% save start timestart_time=cputime;% read in the cover objectfile_name=E:PhdMyPaperSecond paper - ExperimentExperimentCover and Secretlmagele na.bmp：cover_object,map=imread(file _n ame);% read in the secret imagefile_name=E:PhdMyPaperSecond paper - ExperimentExperimentCover and S

2、ecretImagehat.bmp;secret,map1=imread(file_ name);secret仁secret;% convert to double for n ormalizatio n, the n back aga insecret=double(secret);secret=fix(secret./2);%rounds the eleme nts of A toward zero, result ing in an array of in tegerssecret=uin t8(secret);% determ ine size of cover objectMc=si

3、ze(cover_object,1);%HeightNc=size(cover_object,2);%Width% determ ine size of secret objectMm=size(secret,1);%HeightNm=size(secret,2);%Width% title the secret object out to cover object size to gen erate watermarkfor ii = 1:Mcfor jj = 1:Ncif (secret(ii,jj)75)watermark(ii,jj)=255;elsewatermark(ii,jj)=

4、0;endendend% now we set the lsb of cover_object(ii,jj) to the value of watermark(ii,jj) watermarked_image=cover_object;for ii = 1:Mcfor jj = 1:Ncif (watermark(ii,jj)=O)atatwatermarkedm age(ii,jj)=bitset(watermarkedm age(ii,jj),1,0);%Set bitcerta in positi onelsewatermarkedm age(ii,jj)=bitset(waterma

5、rkedm age(ii,jj),1,1);%Set bitcerta in positi onendendend% write the watermarked image out to a fileimwrite(watermarkedmage,lsb_watermarked.bmp,bmp);%=Extracti on=% determ ine size of watermarked imageMw=size(watermarkedmage,1); %HeightNw=size(watermarkedmage,2); %Width% use lsb of watermarked image

6、 to recover watermarkfor ii = 1:Mwfor jj = 1:NwRecoveryWatermark(ii,jj)=bitget(watermarkedm age(ii,jj),1);endend% scale the recovered watermarkRecoveryWatermark=2*double(RecoveryWatermark);%=%=SharpAttack= SharpedWatermarkImage = watermarkedm age;H = fspecial(u nsharp,0.001);%the alpha from 0.005 to

7、 0.1 20 numberssharpImage = imfilter(SharpedWatermarklmage,H,replicate);figure(6);colormap(gray(256);imshow(sharpImage);title(sharp Image);imwrite(sharplmage, gray(256), sharplmage.bmp);% determ ine size of watermarked imageMSw=size(sharpImage,1); %HeightNSw=size(sharplmage,2); %Width% use lsb of wa

8、termarked image to recover watermarkfor ii = 1:MSwfor jj = 1:NSwRecoverySharpWatermark(ii,jj)=bitget(abs(sharplmage(ii,jj),1);endend% scale the recovered watermarkRecoverySharpWatermark=2*double(RecoverySharpWatermark);psn rSharp = psn r(RecoveryWatermark, RecoverySharpWatermark)figure(7);colormap(g

9、ray(256);imshow(RecoverySharpWatermark,)title(RecoverySharpWatermark);%=%=Gaussia n Blur=Gaussia nBlurWatermarkedlmage = watermarked_image;H = fspecial(gaussian,size(GaussianBlurWatermarkedlmage),0.035); %the alpha from 0.005 to 0.1 20 numbersGaussia nBlurlmage = imfilter(Gaussia nBlurWatermarkedlma

10、ge,H,replicate);figure(8);colormap(gray(256);imshow(Gaussia nBlurlmage)title(Gaussia n Blur Image);imwrite(Gaussia nBlurlmage, gray(256), Gaussia nBlurlmage.bmp);% determ ine size of watermarked imageMSw=size(Gaussia nBlurlmage,1); %HeightNSw=size(Gaussia nBlurlmage,2); %Width% use lsb of watermarke

11、d image to recover watermarkfor ii = 1:MSwfor jj = 1:NSwRecoveryGaussia nBlurWatermark(ii,jj)=bitget(Gaussia nBlurlmage(ii,jj),1);end end% scale the recovered watermarkRecoveryGaussia nBlurWatermark=2*double(RecoveryGaussia nBlurWatermark);psn rBlur = psn r(RecoveryWatermark, RecoveryGaussia nBlurWa

12、termark) figure(9);colormap(gray(256);imshow(RecoveryGaussia nBlurWatermark) title(RecoveryGaussia nWatermark);% % %=% %HistogramhistogramWatermarkedlmage = watermarkedm age;histogramimage = histeq(watermarked_image,10); % from 10 to 200 figure(10);colormap(gray(256);image(histogramimage);title(hist

13、ogramimage Image);imwrite(histogramimage, gray(256), histogramimage.bmp);MSw=size(histogramlmage,1);%HeightNSw=size(histogramlmage,2);%Width% use lsb of watermarked image to recover watermarkfor ii = 1:MSwfor jj = 1:NSwRecoveryHistogramWatermarked(ii,jj)=bitget(histogramlmage(ii,jj),1); endend% scal

14、e the recovered watermarkRecoveryHistogramWatermarked=2*double(RecoveryHistogramWatermarked); psn rHistorgram = psn r(RecoveryWatermark, RecoveryHistogramWatermarked) figure(11);colormap(gray(256);imshow(RecoveryHistogramWatermarked,) title(RecoveryHistogramWatermarked);% %=Noise=NoiseWatermarked =

15、watermarked_image;Noiseimage = imno ise(NoiseWatermarked,gaussia n,0.001);figure(12);colormap(gray(256);imshow(NoiseImage);title( no iselmage);imwrite(Noiselmage, gray(256), no iselmage.bmp);MSw=size(Noiselmage,1); %HeightNSw=size(Noiselmage,2); %Width% use lsb of watermarked image to recover watermarkfor ii = 1:MSwfor jj = 1:NSwRecoveryNoiseWatermarked(ii,jj)=bitget(Noiselmage(ii,jj),1);endend% scale the recovered watermarkRecoveryNoiseWatermarked=2*double(R