自己用matlab写的jpeg压缩，解压程序

自己用matlab写的灰度图像压缩解压程序 jpegencoder可以打开图像压缩后存储为c.jpg（只限于灰度图像） jpegdecoder只适用于把jpegencoder的c.jpg解压，还不可以对任意图像解压，有待以后完善 作者 ar_dong

clear all
close all
clc
bmp1=imread（‘c.jpg‘）;
filename=‘c.jpg‘;
fid = fopen（filename‘r‘）;
a=fread（fid2）;
if a（1）~=255
    error（‘%s不是图片‘‘a.jpg‘）
end
if a（2）~=216
     error（‘%s不是图片‘‘a.jpg‘）
end
a=double（fread（fid4））;
if a（1）~=255
    error（‘%s不是图片‘‘a.jpg‘）
end
if a（2）~=224
    error（‘%s图片无法解析‘‘a.jpg‘）
end
count=（a（3）*256+a（4）-2）;
a=double（fread（fidcount））;
a=double（fread（fid4））;
if a（1）~=255
    error（‘%s不是图片‘‘a.jpg‘）
end
if a（2）~=219
    error（‘%s图片无法解析‘‘a.jpg‘）
end
count=（a（3）*256+a（4）-3）;
a=fread（fid1）;
if a（1）~=0
    error（‘%s图片无法解析‘‘a.jpg‘）
end

if count~=64
    error（‘%s图片无法解析‘‘a.jpg‘）
end

a=double（fread（fidcount））;
table=zeros（8）;
huffman00=zeros（655363）;
huffman10=zeros（655363）;
table=[ a（01）a（02）a（06）a（07）a（15）a（16）a（28）a（29）
        a（03）a（05）a（08）a（14）a（17）a（27）a（30）a（43）
        a（04）a（09）a（13）a（18）a（26）a（31）a（42）a（44）
        a（10）a（12）a（19）a（25）a（32）a（41）a（45）a（54）
        a（11）a（20）a（24）a（33）a（40）a（46）a（53）a（55）
        a（21）a（23）a（34）a（39）a（47）a（52）a（56）a（61）
        a（22）a（35）a（38）a（48）a（51）a（57）a（60）a（62）
        a（36）a（37）a（49）a（50）a（58）a（59）a（63）a（64）];
a=fread（fid4）;
if a（1）~=255
    error（‘%s不是图片‘‘a.jpg‘）
end
if a（2）~=192
    error（‘%s图片无法解析‘‘a.jpg‘）
end
count=（a（3）*256+a（4）-2）;
a=double（fread（fidcount））;
if a（1）~=8
    error（‘%s图片无法解析‘‘a.jpg‘）
end
hight=a（2）*256+a（3）
whide=a（4）*256+a（5）
bmp=zeros（hightwhide）;
b=[11170]‘;
if isequal（a（6:9）b）
    a=0;
else
    error（‘%s图片无法解析‘‘a.jpg‘）
end
a=fread（fid4）;
if a（1）~=255
    error（‘%s不是图片‘‘a.jpg‘）
end
if a（2）~=196
    error（‘%s图片无法解析‘‘a.jpg‘）
end
count=（a（3）*256+a（4）-2-1-16）;
a=fread（fid1）;
b=fread（fid16）;
b1=b;
c=fread（fidcount）;
count=0;
clear huffman;
j=1;
for i=1:16
    while b（i）>0
        huffman（j1）=count;
        huffman（j2）=i;
        huffman（j3）=c（j）;
        j=j+1;
        count=count+1;
        b（i）=b（i）-1;
    end
    count=count*2;
end
%%huffman
fanhuffman=zeros（655363）;
i=1;
e=size（huffman）;
e=e（1）; 
while huffman（i2）<17
    b=huffman（i1）;
    c=huffman（i2）;%%长度
    d=huffman（i3）;%%类
    addr_start=b*2^（16-c）;
    addr_end=（b+1）*2^（16-c）-1;
    huffman（i4）=addr_start;
    huffman（i5）=addr_end;
    fanhuffman（addr_start+1:addr_end+11）=floor（d/16）;%%段
    fanhuffman（addr_start+1:addr_end+12）=rem（d16）;%%类
    fanhuffman（addr_start+1:addr_end+13）=c;%%基础码长
    i=i+1;
    if i>e
        break
    end  
end

if a==0
    huffman00=fanhuffman;
elseif a==1
    huffman01=fanhuffman;
elseif a==16
    huffman10=fanhuffman;
elseif a==17
    huffman11=fanhuffman;
else 
    error（‘%s图片存在错误‘‘a.jpg‘）
end
a=fread（fid4）;
if a（1）~=255
    error（‘%s不是图片‘‘a.jpg‘）
end
if a（2）~=196
    error（‘%s图片无法解析‘‘a.jpg‘）
end
count=（a（3）*256+a（4）-2-1-16）;
a=fread（fid1）;
b=fread（fid16）;
b1=b;
c=fread（fidcount）;
count=0;
clear huffman;
j=1;
for i=1:16

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      17534  2011-07-27 08:20  jpg压缩解压\a.jpg

     文件      44571  2011-07-21 09:28  jpg压缩解压\b.jpg

     文件    1440054  2010-02-28 02:09  jpg压缩解压\Bliss.bmp

     文件      31740  2011-08-07 15:18  jpg压缩解压\c.jpg

     文件      17406  2011-07-27 16:02  jpg压缩解压\c.jpg.bak

     文件     481078  2011-08-07 15:17  jpg压缩解压\gray.bmp

     文件        264  2011-07-27 08:59  jpg压缩解压\huffman1.mat

     文件       1168  2011-07-27 09:00  jpg压缩解压\huffman2.mat

     文件       7023  2011-07-27 16:11  jpg压缩解压\jpegdecoder_dong.m

     文件       6632  2011-08-07 15:03  jpg压缩解压\jpegencoder_dong.asv

     文件       6710  2011-08-07 15:18  jpg压缩解压\jpegencoder_dong.m

     文件        256  2011-07-27 08:50  jpg压缩解压\table.mat

    ..A.SH.     24064  2011-07-29 17:26  jpg压缩解压\Thumbs.db

     文件         64  2011-07-27 15:56  jpg压缩解压\Untitled.m

     文件        233  2011-07-27 08:19  jpg压缩解压\Untitled2.m

     文件    1624502  2011-07-27 10:00  jpg压缩解压\未命名.bmp

     文件       2113  2011-07-27 10:00  jpg压缩解压\未命名.GIF

     文件       9211  2011-07-27 10:00  jpg压缩解压\未命名.JPG

     文件       4525  2011-07-27 10:00  jpg压缩解压\未命名.PNG

     文件      10176  2011-07-27 10:00  jpg压缩解压\未命名.TIF

     文件      68734  2011-07-27 10:00  jpg压缩解压\未命名1.bmp

     文件        193  2011-08-07 15:22  jpg压缩解压\说明.txt

     目录          0  2011-08-07 15:19  jpg压缩解压

----------- ---------  ---------- -----  ----

              3798251                    23