修改最大像素值的可逆数字算法的实现
论文
政研论文下载论文大学下载论文大学下载关于长拳的论文浙大论文封面下载
的算法，继续与大家分享 修改最大像素值的可逆数字算法的实现 %designed by 李立宗
%lilizong@gmail.com
%t4311.m
c=imread('e:\lena.bmp');
subplot(1,3,1),imshow(c),title('嵌入水印前图像');
[m,n]=size(c);
dix=zeros(255);
for i=1:512
for j=1:512
for k=1:255
if c(i,j)==k
dix(k)=dix(k)+1;
end
end
end
end
% subplot(3,3,2),plot(dix);
max=dix(1);
for i=1:255
if max dix(i)
max=dix(i);
maxi=i;
end
end
max
c1=c;
for i=1:512
for j=1:512
if c(i,j) maxi
c1(i,j)=c(i,j)-1;
else
c1(i,j)=c(i,j);
end
end
end
dix2=zeros(255);
for i=1:512
for j=1:512
for k=1:255
if c1(i,j)==k
dix2(k)=dix2(k)+1;
end
end
end
end
% subplot(3,3,3),plot(dix2);
% subplot(3,3,4),plot(dix,'r'),hold on,plot(dix2,'g:');
mark=imread('e:\logo1.bmp');
mark=im2bw(mark);
subplot(1,3,2),imshow(mark),title('水印图像');
%读入水印图像，并将之二值化处理
% figure,imshow(mark),title('水印图像');
mark1=zeros(1,1600);
n=1;
for i=1:40
for j=1:40
mark1(n)=mark(i,j);
n=n+1;
end
end
% mark1=reshape(mark,1,1600);
%测试mark1的正确与否
% n=1;
% for i=1:40
% for j=1:40
% mark2(i,j)=mark1(n);
% n=n+1;
% end
% end
% subplot(1,3,3),imshow(mark2),title('水印图像');
c2=zeros(512,512);
c2=c1;
n=0;
for i=1:512
for j=1:512
if c2(i,j)==maxi
n=n+1;
if(n==1600)
n=1;
end
if mark1(n)==1
c2(i,j)=c2(i,j)-1;
end
end
end
end
% n=1;
% for i=1:512
% for j=1:512
% ifc2(i,j)==maxi mark1(n)==1
% c2(i,j)=c2(i,j)-1;
% n=n+1;
%
% end
% end
% end
subplot(1,3,3),imshow(c2),title('嵌入水印后图像');
markshow=zeros(1,5000);
%分配一个大数组，这样可以存储冗余量
n=1;
for i=1:512
for j=1:512
if c2(i,j)==maxi
markshow(n)=0;
n=n+1;
elseif c2(i,j)==(maxi-1)
markshow(n)=1;
n=n+1;
end
end
end
markshow;
markshow1=zeros(40,40);
%定义一个数组，用来存储恢复的水印
n=1;
for i=1:40
for j=1:40
markshow1(i,j)=markshow(n);
n=n+1;
end
end
% subplot(1,3,3),imshow(markshow1),title('恢复水印图像');
c3=c2;
for i=1:512
for j=1:512
if c2(i,j) maxi
c3(i,j)=c2(i,j)+1;
else
c3(i,j)=c2(i,j);
end
end
end
% subplot(1,3,3),imshow(c3),title('恢复原始载体图像');
浙公网安备 33021202002483