%參數說明：
%cover是載體圖像的地址
%message是秘密圖像的地址
%stegocover是隱藏后圖像的地址
%level是作為載體的機體層，R為1，G為2，B為3
%data是隱藏后圖像的矩陣
function data=imagehide(cover,message,stegocover,permission,level)
%提取圖像信息并分層
cover=imread(cover,permission);
data=cover;
msg=imread(message,permission);
[row,col]=size(cover);
cover1=cover(:,:,level);
%置載體圖像R層的低四位為0
for i=1:row
    for j=1:col/3
        cover1(i,j)=bitand(cover1(i,j),240);
    end
end

%置秘密圖像的低四位為0
takemsg4=bitand(msg,240);
%將秘密圖像的高四位右移四位
shiftmsg4=bitshift(takemsg4,-4);
%圖像隱藏
for i=1:row
    for j=1:col/3
        cover1(i,j)=bitor(cover1(i,j),shiftmsg4(i,j));
    end
end

%寫回并保存
data(:,:,level)=cover1;
imwrite(data,stegocover,permission);

%提取秘密圖像信息，檢測隱藏效果
data=imread(stegocover,permission);
[row,col]=size(data);
A=data(:,:,level);
for i=1:row
    for j=1:col/3
        A(i,j)=bitand(A(i,j),15);
    end
end
A=bitshift(A,4);
%顯示結果
subplot(221);imshow(cover);title('兩幅圖像中所選的載體圖像');
subplot(222);imshow(message);title('兩幅圖像中所選的秘密圖像');
subplot(223);imshow(data);title('將秘密圖像隱藏在載體圖像后所得到的的圖像');
subplot(224);imshow(A);title('從處理后的載體圖像中提取出的秘密圖像');
% figure;
% t=double(cover(:,:,level))-double(data(:,:,level));
% imshow(t,[]);





%參數說明：
%cover是載體圖像的地址
%message是秘密圖像的地址
%stegocover是隱藏后圖像的地址
%level是作為載體的機體層，R為1，G為2，B為3
%count是替換表
%thershold是閾值
function count=fourthbitcmp(cover,message,permission,level,threshold)
%提取圖像信息并分層
cover=imread(cover,permission);
data=cover;
msg=imread(message,permission);
cover1=cover(:,:,level);
%對cover和msg的第四位進行處理
tempc=cover1;
tempm=msg;
tempc=bitand(tempc,8);
tempm=bitand(tempm,8);
temp=bitxor(tempm,tempc);
[row,col]=size(temp);

%記錄圖像每個分塊的n值
k1=0;
k2=0;
a=row*col/64;
count=zeros([1 a]);
for i=1:a
    for m=1:8
        for n=1:8
            if temp(8*k1+m,8*k2+n)==0
                count(1,i)=count(1,i)+1;
            end
        end
    end
    k2=k2+1;
    if k2*8==col
        k2=0;
        k1=k1+1;
    end
end

%計算每塊的u值與閾值進行比較
count=count/64;
for i=1:a
    if count(i)>=threshold
        count(i)=1;  %可以替換
    elseif count(i)<1-threshold
        count(i)=-1;  %取反
    else
        count(i)=0;%不能處理
    end
end

s=0;
for i=1:a
    if count(i)==0
        s=s+1;        
    end    
end
disp(s);