function LSB
%LSB
%圖像水印的LSB算法
clear all;
close all;
c=zeros(200*200/8,8);
%行為文本的字符數目，列為8，每個字符的ASCILL值。
A=zeros(256,256);
B=zeros(200,200);
C=zeros(200,200);
%初始化 
 
subplot(2,2,1);
grid on
[A,map]=imread('lena.BMP','bmp'); 
imshow(A,[]);
title('原始圖像');
%顯示原始圖象 
 
for i=1:200 
  for j=1:200 
    temp=A(i,j);
    B(i,j)=bitset(temp,1,0); 
  end 
end 
subplot(2,2,2);
imshow(B,[]);
title('預處理后的圖象'); 
%顯示LSB置位0后的圖象 
 
fid = fopen('temp.txt','r');
F = fread(fid,'uchar');
maxI=length(F); 

for i=1:maxI 
  a=F(i);  
  for j=1:8 
    c(i,9-j)=bitget(a,j); 
  end 
end 
fclose(fid)
%讀入并存儲待隱藏文件 
 
 
c=reshape(c,200,200); 
for i=1:200
  for j=1:200 
    temp=B(i,j); 
    C(i,j)=bitset(temp,1,c(i,j)); 
  end 
end 
subplot(2,2,3); 
imshow(C,[]); 
title('嵌入文本后的圖象'); 
%顯示嵌入文本文件后的圖象 

S=C; 
%含有文本信息的圖象S 
for i=1:200 
  for j=1:200 
    t(i,j)=bitget(c(i,j),8); 
  end 
end 
t=reshape(t,200*200/8,8); 
T=setstr(F') %提取出來的文本放在T