EMD算法原理:
MATLAB代码如下:
clear;
clc;
pout=imread('pout.tif');
[m,n]= size(pout);%读取出图片的列像素个数和行像素个数
S = floor(rand(1,(m*n)/2)*5);%产生一组二分之一像素个数的0—4的随机数,作为待隐藏的信息
carrier = reshape(pout,1,m*n);%将载体图片转化为一维数组
figure(1);
subplot(1,2,1),imshow(pout), title('原始图像');%显示原始图像;
%加密算法
for i=1:2:m*n
g1=double(carrier(i));
g2=double(carrier(i+1));
f=mod((g1+g2*2),5);
r=mod((S((i+1)/2)-f),5);
if(r==0) %r=0时,g1,g2保持不变
g1=g1;
g2=g2;
elseif(r<=2) %r小于等于二时,gr+1
if(r==1)
g1=g1+1;
else
g2=g2+1;
end
elseif(r>2) %r大于2时,g(5-r)—-5
if(5-r==1)
g1=g1-1;
else
g2=g2-1;
end
end
newcarrier(i)=uint8(g1);
newcarrier(i+1)=uint8(g2);
end
newimg=reshape(newcarrier,m,n);
subplot(1,2,2),imshow(newimg),title('嵌入信息后');%打印嵌入信息后的图片
%提取算法
for k=1:2:m*n
b1=double(newcarrier(k));
b2=double(newcarrier(k+1));
information((k+1)/2)=mod((b1+b2*2),5);
end
执行效果:
图像的像素发生了变化:
提取信息结果:
实验中踩的坑:
图像中均采用uint8的类型,f函数运算溢出的话,模将为0。所以在运算f之前,需要把g1,g2变量转换为其他类型。到显示新图像之前还要将其转换为uint8类型。