function yasuobi=hide_2(A1,info) %每隔两个DCT系数置换一个
A=double(A1);
%分块
[p q]=size(A);
a=ceil(p/8);
b=ceil(q/8);
block=zeros(8,8,a*b);
for m=1:a
for n=1:b
%每一小块
for s=1:8
for t=1:8
block(s,t,m*b-b+n)=A(8*m-8+s,8*n-8+t);
end
end
end
end
DCT=zeros(8,8,a*b);
%DCT变换
for m=1:a*b
DCT(:,:,m)=dct2(block(:,:,m));
end
%量化
load JpegCoeff;
NUM=zeros(8,8,a*b);
for m=1:a*b
NUM(:,:,m)=round(DCT(:,:,m)./QTAB);
end
dec=info_dec(info);
dec_double=double(dec)-48;
length_dec=length(dec);
num_block=ceil(length_dec*3/64);
count=1;x=1;
for m=1:num_block+1
for n=1:8
for t=1:8
if(count<=3*length_dec+1)
if(mod(count,3)==1)
sign=1;
if(NUM(n,t,m)<0)
sign=-1;
end
if(count~=3*length_dec+1)
NUM(n,t,m)=sign*(floor(abs(NUM(n,t,m))/10)*10+dec_double(x));
x=x+1;
end
if(count==3*length_dec+1)
NUM(n,t,m)=sign*(floor(abs(NUM(n,t,m))/10)*10+5);
end
end
end
count=count+1;
end
end
end
%zig_zag
Zig=zeros(64,a*b);
for m=1:a*b
Zig(:,m)=zig_zag_2(NUM(:,:,m));
end
%对DC 系数先做差分编码
for m=1:a*b
if(m==1)
CD(m)=Zig(1,m);
else
CD(m)=Zig(1,m-1)-Zig(1,m);
end
end
%DC编码
DC_code='';
for m=1:length(CD)
DC_code=strcat(DC_code,dc_code(CD(m)));
end
%AC编码
AC_code='';
for m=1:a*b
AC_code=strcat(AC_code,ac_code(Zig(2:64,m)));
end
yasuobi=yasuo(AC_code,DC_code,A1);
save hide_2.mat p q DC_code AC_code;
%end
matlab实现数字图像处理
需积分: 18 108 浏览量
2011-01-15
20:01:33
上传
评论
收藏 3.97MB RAR 举报
matlab实现数字图像处理 （147个子文件） 
hide_2.asv 2KB iAC_code.asv 2KB hide_1.asv 2KB hide_3.asv 2KB ac_code.asv 2KB de_hide_2.asv 1KB hide_gray.asv 1KB de_hide_1.asv 1KB pic_4_4_2.asv 1KB de_hide_3.asv 1KB jpeg_code.asv 1KB iDC_code.asv 979B dec_info.asv 873B ijpegcode.asv 865B color_L.asv 715B picture.asv 644B pic_4_2_8.asv 638B iHuffman_AC.asv 607B check.asv 600B de_hide_gray.asv 577B pic_4_2_4.asv 519B my_DCT.asv 464B pic_4_3_2.asv 426B zig_zag_1.asv 414B pic_4_2_3.asv 374B pic_4_2_1.asv 308B is_huffman.asv 277B pic_4_2_7.asv 216B pic_4_2_11.asv 169B code.asv 147B pic_4_3_1.asv 144B distance.asv 114B huffman.asv 113B pic_4_2_12.asv 89B pic_4_2_9.asv 63B pic_4_2_10.asv 63B Category.asv 61B pic_4_2_2.asv 0B
pic_4_4_3_imresize.bmp 1.72MB pic_4_4_3_imrotate.bmp 884KB pic_4_4_1.bmp 883KB pic_4_4_3_imadjust.bmp 883KB 16.bmp 558KB 28.bmp 179KB 9.bmp 153KB 24.bmp 83KB hall_black_white.bmp 59KB hall_red_circle.bmp 59KB 26.bmp 53KB 11.bmp 47KB 17.bmp 38KB 14.bmp 36KB 18.bmp 30KB 25.bmp 27KB 19.bmp 27KB hall.bmp 21KB hide_1.bmp 21KB hall_transpose.bmp 21KB hall_90.bmp 21KB hall_rebuilt.bmp 21KB hall_180.bmp 21KB 20.bmp 16KB 15.bmp 13KB 27.bmp 12KB 23.bmp 10KB 29.bmp 10KB 12.bmp 9KB 31.bmp 9KB 10.bmp 9KB 30.bmp 8KB 13.bmp 8KB 21.bmp 6KB 4.bmp 6KB 22.bmp 5KB 3.bmp 5KB 6.bmp 4KB 5.bmp 4KB 1.bmp 4KB 8.bmp 3KB 2.bmp 3KB 7.bmp 3KB hs_err_pid4040.log 16KB hide_2.m 2KB hide_1.m 2KB iAC_code.m 2KB hide_3.m 2KB ac_code.m 2KB pic_4_4_3.m 1KB de_hide_2.m 1KB de_hide_1.m 1KB jpeg_code_2.m 1KB jpeg_code.m 1KB de_hide_3.m 1KB iDC_code.m 983B pic_4_1_2.m 950B dec_info.m 945B ijpegcode_2.m 876B ijpegcode.m 870B pic_4_2_4.m 840B color_L.m 689B共 147 条
- 1
- 2
资源评论
