EMD信息隱藏算法的MATLAB實現（基礎）

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類型。