%---------水印信息的制造------------------
I_water = imread('sy.bmp');
Img_gray = imresize(I_water,[512,512]);
Img_bw = im2bw(Img_gray,0.9);
%----------載體圖像的讀入------------------
Img_cover = double(imread('yt.bmp'));
[nrow,ncol] = size(Img_cover);
%----------水印信息的嵌入-------------------
N=8; %嵌入位置（位平面）的選擇
Img_imbeded = zeros(nrow,ncol);%產生nr行nc列的零矩陣
for i = 1:nrow
    for j = 1:ncol
       Img_bin = dec2bin(Img_cover(i,j),8);%把水印圖像轉換為8位二進制圖像
         if Img_bw(i,j)==0;%判斷水印第（i，j）點是否為0
          Img_bin(N)='0';%源圖像第8位平面置為‘0’
         else
          Img_bin(N)='1';%源圖像第8位平面置為‘1’
         end;%源圖像第(i,j)像素點嵌入完畢
         Img_imbeded(i,j) = bin2dec(Img_bin);
    end
end
%結果顯示
subplot(121);
imshow(uint8(Img_cover));
title('載體圖像')
subplot(122);
imshow(uint8(Img_imbeded));
title('偽裝圖像')