CN102547295A - 一种图像数据无损压缩的编码方法 - Google Patents

一种图像数据无损压缩的编码方法 Download PDF

Info

Publication number
CN102547295A
CN102547295A CN2012100407057A CN201210040705A CN102547295A CN 102547295 A CN102547295 A CN 102547295A CN 2012100407057 A CN2012100407057 A CN 2012100407057A CN 201210040705 A CN201210040705 A CN 201210040705A CN 102547295 A CN102547295 A CN 102547295A
Authority
CN
China
Prior art keywords
array
difference matrix
order difference
code stream
coding
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012100407057A
Other languages
English (en)
Other versions
CN102547295B (zh
Inventor
高健
饶珺
刘万
宋奥
王忠元
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Shanghai for Science and Technology
Original Assignee
University of Shanghai for Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Shanghai for Science and Technology filed Critical University of Shanghai for Science and Technology
Priority to CN201210040705.7A priority Critical patent/CN102547295B/zh
Publication of CN102547295A publication Critical patent/CN102547295A/zh
Application granted granted Critical
Publication of CN102547295B publication Critical patent/CN102547295B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

本发明公开了一种图像数据无损压缩的编码方法方法，它主要包括：从图像矩阵中提取一阶混合差分矩阵；依一阶混合差分矩阵元素的统计分布，确定合并次数与编码截断阶数；根据合并次数与编码截断阶数对一阶混合差分矩阵元素进行一次编码，实现一次压缩；对一次压缩编码二值码流进行二次编码，实现二次编码压缩。该方法能根据图像数据的统计特征选取压缩量最大情形下所对应的数据合并次数q与编码截断阶数p，依所选合并次数和编码截断阶数构造一种即时编码，灵活性强，算法简单，易于实现，编码效率高；该方法的一次压缩编码表的构成简单；该方法的编码前缀及编码后缀结构统一，解码过程无需匹配，解码效率高。

Description

一种图像数据无损压缩的编码方法
技术领域
本发明涉及一种图像数据无损压缩的编码方法，属于数字图像无损压缩技术领域。
背景技术
在信息环境中，信息大部分是以数字化的方式表示、存储和传输的，其中图像占了绝大部分。但是数字化的图像需要占用巨大的空间，为了对图像数据进行有效的处理、存储和传输，必须对图像数据进行压缩存储。因此，数字图像信息压缩编码也就成为多媒体信息处理中的关键技术之一。如何在保证无损的前提下尽可能提高图像的压缩效率，对于与图像有关的应用而言，具有非常重要的意义。
赫夫曼方法是前缀即时码中编码效率最优的编码方法，其对信源的统计特性没有特殊要求，对编码设备的要求也比较简单，综合性能优于香农码和费诺码。但其解码过程需通过码表匹配来实现，因此解码效率低下。一些次优编码方法，如B1码，B2码以及移位码等，虽然编解码效率得到提高，但编码形式固定、单一，对信源的统计特性有特殊要求，压缩效率较低。
发明内容
本发明的目的在于针对现有技术存在的问题提出一种图像数据无损压缩的编码方法，该方法具有灵活性强，算法简单，易于实现，编解码效率高。
为了达到上述目的，本发明采用下述技术方案：
一种图像数据无损压缩的编码方法，其特征在于，该方法具体步骤如下：
（1）、从原始图像逐行读取数据生成图像矩阵D m×n，从图像矩阵D m×n 中生成横向一阶差分矩阵X m×n 和纵向一阶差分矩阵Y m×n ；
（2）、从横向一阶差分矩阵X m×n 和纵向一阶差分矩阵Y m×n 中生成混合一阶差分矩阵XY m×n 和混合一阶差分矩阵XY m×n 子块的状态二值码流D；
（3）从混合一阶差分矩阵XY m×n 中逐行提取元素符号，0表示正号，1表示负号，生成混合一阶差分矩阵XY m×n 的符号二值码流S和非负混合一阶差分矩阵Z m×n ；
（4）、统计非负混合一阶差分矩阵Z m×n 中各元素出现的个数，生成元素个数统计数组A;
（5）、从元素个数统计数组A中生成元素映射数组T S ，依元素映射数组T S 对非负混合一阶差分矩阵Z m×n 中各元素重新分配数值；
（6）、对元素个数统计数组A分别进行q=0,1,2,3,4次数据合并，建立5组统计数组A 0,A 1,A 2,A 3,A 4和对应的5组组号数组F 0,F 1,F 2,F 3,F 4；
 （7）、依次对5组统计数组A 0,A 1,A 2,A 3,A 4和5组组号数组F 0,F 1,F 2,F 3,F 4，统计编码截断阶数p分别为1,2,3,…,32阶下的编码比特数，选取编码比特数最小时的合并次数q c 和编码截断阶数p c ；
(8）、逐行提取非负混合一阶差分矩阵Z m×n 元素，依次放入数组E中，依步骤（7）选取的2个参数，即合并次数q c 和编码截断阶数p c ，对数组E中各元素进行合并与编码，生成一次压缩编码二值码流T 1；
（9）、统计一次压缩编码二值码流T 1中连续0和连续1的长度，得到一次压缩编码二值码流T 1的长度分布数组L；
（10）、用长度分布数组L的长度商数组Q和余数二值码流R表示长度分布数组L，实现二次编码压缩，生成二次编码二值码流T 2。
本发明的一种图像数据无损压缩的编码方法具有如下特点：
（1）、该方法能根据图像数据的统计特征选取压缩量最大情形下所对应的数据合并次数q与编码截断阶数p，依所选合并次数和编码截断阶数构造一种即时编码，灵活性强，算法简单，易于实现，编码效率高；（2）、该方法的一次压缩编码表的构成简单，对于数值小于编码截断阶数的压缩元素，由编码前缀，即1元即时码表示，对于数值大于等于编码截断阶数的压缩元素，由编码前缀和编码后缀联合表示；（3）、该方法的一次压缩编码表的主要部分为其前缀编码，前缀编码可按位平面层面分解，并依位平面层面摆放，基于图像数据的连续性，这种二值码流的层面摆放方式为后续的二次压缩编码提供可能，两次编码后的整体压缩效率高；（4）、该方法的编码前缀及编码后缀结构统一，解码过程无需匹配，解码效率高。
附图说明
图1是本发明的一种图像数据无损压缩的编码方法的流程图；
图2是生成混合一阶差分矩阵XY m×n 方法示意图；
图3是数组                                                
Figure 858257DEST_PATH_IMAGE001
在合并次数q c 和编码截断阶数p c 下的编码表；
图4是序列按位平面层面的分解图。
具体实施方式
下面结合附图和256级灰度的任意图像，对本发明作进一步描述。
如图1所示，本发明的一种图像数据无损压缩的编码方法，其具体步骤如下：
步骤（1）、从原始图像逐行读取数据生成图像矩阵D m×n，从图像矩阵D m×n 中生成横向一阶差分矩阵X m×n 与纵向一阶差分矩阵Y m×n ，具体如下：
从宽为m，高为n的原始图像中逐行提取图像象素值，将提取所得象素值置入图像矩阵D m×n 中,其中D m×n =(d ij )  m×n ，
分别计算第i行、第j列处的横向一阶差分值x ij 与纵向一阶差分值y ij ，具体计算公式如下：
      x ij =d ij -d i,j-1   ，i=1,2,…,m,   j=2,3,…,n                       （1）
      y ij =d ij -d i-1,j   ， i=2,3,…,m,   j=1,2,…,n                       （2）
其中，x ij 为横向一阶差分矩阵X m×n  在(i,j)处的横向一阶差分值，y ij 为纵向一阶差分矩阵Y m×n  在(i,j)处的纵向一阶差分值，横向一阶差分矩阵X m×n 中第1列差分元素值、纵向一阶差分矩阵Y m×n 中第1行差分元素值的具体计算公式如下：
x i1=d i1 -d i-1,1     ，i=2,3,…,m                                    （3）
y 1j =d 1j -d 1,j-1    ，j=2,3,…,n                                   （4）
    横向一阶差分矩阵X m×n 和纵向一阶差分矩阵Y m×n 第1行第1列差分元素值置为0；
步骤（2）、从横向一阶差分矩阵X m×n 和纵向一阶差分矩阵Y m×n 中生成混合一阶差分矩阵XY m×n 和混合一阶差分矩阵XY m×n 子块的状态二值码流D，具体如下：
对横向一阶差分矩阵X m×n 和纵向一阶差分矩阵Y m×n 按4×4大小的子块划分，各得
（m×n）/16个子块，分别计算横向一阶差分矩阵X m×n 和纵向一阶差分矩阵Y m×n 相同位置子块元素绝对值之和，取绝对值之和最小的子块对应的各个元素置入混合一阶差分矩阵XY m×n 相同位置子块处，生成一阶混合差分矩阵XY m×n =(xy ij ) m×n ，参见附图2，对横向一阶差分矩阵X m×n 和纵向一阶差分矩阵Y m×n 相同位置子块，记
Figure 713081DEST_PATH_IMAGE002
，
Figure 492818DEST_PATH_IMAGE003
，若
Figure 888027DEST_PATH_IMAGE004
，则将横向一阶差分矩阵X m×n 子块数据置入混合一阶差分矩阵XY m×n 相同位置子块处，此时子块的状态标记为0；若，则将纵向一阶差分矩阵Y m×n 子块数据置入混合一阶差分矩阵XY m×n 相同位置子块处，此时子块的状态标记为1，将子块的状态标记0或1依次置入子块的状态二值码流D中；
步骤（3）、从混合一阶差分矩阵XY m×n 中逐行提取元素符号，0表示正号，1表示负号，生成混合一阶差分矩阵XY m×n 的符号二值码流S和非负混合一阶差分矩阵Z m×n ，具体如下：
从混合一阶差分矩阵XY m×n =(xy ij ) m×n 中逐行提取矩阵的各元素符号，正号元素记为0，负号元素记为1，依次将符号元素0或1置入一阶混合差分矩阵XY m×n 的符号二值码流S，对混合一阶差分矩阵XY m×n =(xy ij ) m×n 各元素取绝对值，生成非负混合一阶差分矩阵Z m×n =(|xy ij |) m×n ；
步骤（4）、统计非负混合一阶差分矩阵Z m×n 中各元素出现的个数，生成元素个数统计数组A，具体如下：
统计非负混合一阶差分矩阵Z m×n =(|xy ij |) m×n 中各元素出现的个数，按元素出现次数的大小置入元素个数统计数组A=(n 0,n 1,…,n e )中，其中n 0≥n 1≥…≥n e ，n 0,n 1,…,n e 分别对应非负混合一阶差分矩阵Z m×n =(|xy ij |) m×n 中元素值a 0,a 1,a 2,…,a e 出现的个数, 0≤e≤255，其中|xy ij |
Figure 462545DEST_PATH_IMAGE006
{a 0,a 1,a 2,…,a e }，表示非负混合一阶差分矩阵Z m×n 中任一点元素值|xy ij |仅在集合{ a 0,a 1,a 2,…,a e }中取值；
步骤（5）、从元素个数统计数组A中生成元素映射数组T S ，依元素映射数组T S 对非负混合一阶差分矩阵Z m×n 各元素重新分配数值，具体如下：
依次将元素个数统计数组A=(n 0,n 1, …,n i ,…,n e )中n i 所对应的非负一阶混合差分矩阵Z m×n 中元素a i ,置入元素映射数组T S 中，生成元素映射数组T S =(a 0,a 1, …,a i ,…,a e )，逐行提取非负混合一阶差分矩阵Z m×n 各元素|xy ij |，|xy ij |为非负混合一阶差分矩阵Z m×n 中第i行第j列元素，如果|xy ij |=a k  ，一阶非负混合差分矩阵Z m×n 中第i行第j列元素值|xy ij |置为k，
其中a k 是映射数组T S =(a 0,a 1,a 2,…,a e )中的第k+1个元素，0≤k≤e；
步骤（6）、对步骤（4）所得元素个数统计数组A分别进行q=0,1,2,3,4次数据合并，建立5组统计数组A 0,A 1,A 2,A 3,A 4和对应的5组组号数组F 0,F 1,F 2,F 3,F 4，具体如下：
为论述方便，假定e=255，则
0次合并统计数组A 0=A=(n 0,n 1,…,n 255)；   
0次合并组号数组F 0=(0,1,2,…,255)；
1次合并统计数组A 1=(n 0+n 1, n 2+n 3…, n 254+n 255)=
Figure 729579DEST_PATH_IMAGE007
；
1次合并组号数组F 1=(0,1,2,…,127)；
2次合并统计数组A 2=(n 0+n 1+ n 2+n 3,…, n 252+n 253+ n 254+n 255)=；
2次合并组号数组F 2=(0,1,2,…,63)；
3次合并统计数组A 3=(n 0+n 1+…+n 7,…, n 248+n 249+…+n 255)=
Figure 36505DEST_PATH_IMAGE009
；
3次合并组号数组F 3=(0,1,2,…,31)；
4次合并统计数组A 4=(n 0+n 1+…+n 15,…, n 240+n 241+…+n 255)=
Figure 29869DEST_PATH_IMAGE010
；
4次合并组号数组F 4=(0,1,2,…,15)，
其中，A为步骤（4）所述元素个数统计数组，组号数组F k 中元素a i 经k+1次合并后对应到组号数组F k+1中元素round[a i /2]，round为取整函数，0≤k≤3；
步骤（7）、依次对5组统计数组A 0,A 1,A 2,A 3,A 4和5组组号数组F 0,F 1,F 2,F 3,F 4，统计编码截断阶数p分别为1,2,3,…,32阶下的编码比特数，选取编码比特数最小时的合并次数q c 和编码截断阶数p c ，具体如下：
编码截断阶数p为可用1元即时二进制码0,10,110,…,11…10，分别表示元素0,1,2,…,p-1，其中，11…10由p-1个1和1个0组成，大于p-1的元素的二进制码，由p个1的二进制码组成其前缀，由b位二进制码组成其后缀，其中参数b由公式（6）确定，
如果编码截断阶数为p，合并次数为q，对图像矩阵D m×n 编码占用的比特数G为
Figure 784199DEST_PATH_IMAGE012
   （5）
其中， B 0(D)为步骤（2）所述子块的状态二值码流D占用比特数；B 1(S)为步骤（3）所述混合一阶差分矩阵元素的符号二值码流S占用比特数；B 2(T s )为步骤（5）所述元素映射数组T s 占用比特数，
Figure 458894DEST_PATH_IMAGE013
对应的组号元素为
Figure 982279DEST_PATH_IMAGE014
，参数b按如下公式确定：
记
Figure 412123DEST_PATH_IMAGE015
，
如果0≤f≤1,则b=1；
如果2≤f≤3,则b=2；
如果4≤f≤7,则b=3；
如果8≤f≤15,则b=4；                                              （6）
如果16≤f≤31,则b=5；
如果32≤f≤63,则b=6；
如果64≤f≤127,则b=7；
如果128≤f≤255,则b=8，
对0,1,2,3,4次合并后的5组统计列表A 0,A 1,A 2,A 3 ,A 4和5组组号数组F 0,F 1,F 2,F 3,F 4，分别统计编码截断阶数为p=1,2, …,32阶下占用的编码比特数，由5种合并，32个编码截断阶数，共计得5×32=160种组合形式，按公式（5）分别计算160种组合的编码占用的比特数G，从160种组合中选取编码占用的比特数G最小时对应的编码截断阶数p=p c 和元素合并次数q=q c ，其中1≤p c ≤32，0≤q c ≤4；
步骤（8）、逐行提取非负混合一阶差分矩阵Z m×n 元素，依次放入数组E中，依步骤（7）选取的2个参数，即合并次数q c 和编码截断阶数p c ，对数组E中各元素进行合并与编码，生成一次压缩编码二值码流T 1,具体如下：
逐行提取矩阵Z m×n =(|xy ij |) m×n 元素，依次放入数组E中，E的表达式为
E=(|xy 11|,|xy 12|,…, |xy 1n |,|xy 21|,|xy 22|,…, |xy 2n |,…, |xy  m1|,|xy  m2|,…, |xy  mn |)，
 按步骤（7）所选合并次数q c ，对数组E中各元素整除
Figure 591432DEST_PATH_IMAGE016
，商置入数组
Figure 866555DEST_PATH_IMAGE001
，余数置入数组，
Figure 782876DEST_PATH_IMAGE017
中各余数的二进制表示位数为q c  位，参照附图3，对数组
Figure 511797DEST_PATH_IMAGE001
中各元素进行编码，小于p c 的元素0,1,2,…,p c -1对应的编码为；大于p c 的元素对应的编码为，其中参数b的取值由步骤（7）所述公式（6）确定，
将各元素编码的前缀部分，即1元即时码0,10,110,…,11…10,11…11，依数组
Figure 524807DEST_PATH_IMAGE001
中元素顺序进行纵向分解，按位平面层面依次置入一次压缩编码二值码流T 1中，接着在一次压缩编码二值码流T 1中置入后缀部分的b位二进制数，然后在一次压缩编码二值码流T 1中置入余数数组
Figure 296454DEST_PATH_IMAGE017
对应的二值码流、步骤（3）所述混合一阶差分矩阵XY m×n 的符号二值码流S和步骤（2）所述混合一阶差分矩阵XY m×n 子块的状态二值码流D，生成一次压缩编码二值码流T 1。
例如某序列为{3,3,2,3,5,4,0,0,4,1,2,1,2,6,2,1,1,1,2,1}，假定数据合并次数为0，编码截断阶数为5，参照附图4，前缀部分分解为5层位平面，1层位平面为11111100111111111111，2层位平面为11111110101110001，3层位平面110111100100，4层位平面为000111，5层位平面1001，按位平面层面的摆放结果如下：11111100111111111111,11111110101110001,  110111100100, 000111, 1001；
步骤（9）、统计一次压缩编码二值码流T 1中连续0和连续1的长度，得到一次压缩编码二值码流T 1的长度分布数组L，具体如下：
统计步骤（8）一次压缩编码二值码流T 1中0段和1段的长度，得到长度分布数组L=(l 1,l 2,l 3,…,l r)，
如一次压缩编码二值码流T 1=(00010000111110011100001111111111000001010111),则长度分布数组L=(3,1,4,5,4,10,5,1,1,1,1,3)；
步骤（10）、用长度分布数组L的长度商数组Q和余数二值码流R表示长度分布数组L，实现二次编码压缩，生成二次编码二值码流T 2,具体如下：
对长度分布数组L=(l 1,l 2,…,l r )中各元素按公式（7）计算，得到长度商数组Q和余数二值码流R，其中
Q=round((l 1+1)/2,(l 2+1)/2,…,(l r +1)/2)                          （7）
round表示取整函数，余数置入余数二值码流R=(b 1,b 2,…, b i ,…,b r )中，其中，当(l i +1)/2为偶数时，b i 为0，当(l i +1)/2为奇数时，b i 为1。
长度商数组Q中元素为1时，编码为0或1，如果该原始长度由0组成，则编码为0，如果该原始长度由1组成，则编码为1；长度商数组Q中元素为2时，编码为00或11，如果该原始长度由0组成，则编码为00，如果该原始长度由1组成，则编码为11。如某二值码流串0111111000111111111111，其长度分布数组L=(1,6,3,12)，长度商数组Q编码为011100111111，余数二值码流R为0101。
将长度商数组Q的二进制编码与余数二值码流R组合后，生成二次编码二值码流T 2。
为了验证使用本发明的方法压缩效果，通过对Lena、Columbia等100多幅图像进行压缩测试，验证了本发明的方法压缩效果，实验结果得出，对于70%以上图像而言，经二次编码后，图像每象素的平均比特数低于步骤（2）所述的混合一阶差分矩阵XY m×n 所对应的信源数据的香农熵。
以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围。

Claims (9)

1.一种图像数据无损压缩的编码方法，其特征在于，包括如下步骤：
（1）、从原始图像逐行读取数据生成图像矩阵D m×n，从图像矩阵D m×n 中生成横向一阶差分矩阵X m×n 和纵向一阶差分矩阵Y m×n ；
（2）、从横向一阶差分矩阵X m×n 和纵向一阶差分矩阵Y m×n 中生成混合一阶差分矩阵XY m×n 和混合一阶差分矩阵XY m×n 子块的状态二值码流D；
（3）从混合一阶差分矩阵XY m×n 中逐行提取元素符号，0表示正号，1表示负号，生成混合一阶差分矩阵XY m×n 的符号二值码流S和非负混合一阶差分矩阵Z m×n ；
（4）、统计非负混合一阶差分矩阵Z m×n 中各元素出现的个数，生成元素个数统计数组A;
（5）、从元素个数统计数组A中生成元素映射数组T S ，依元素映射数组T S 对非负混合一阶差分矩阵Z m×n 中各元素重新分配数值；
（6）、对元素个数统计数组A分别进行q=0,1,2,3,4次数据合并，建立5组统计数组A 0,A 1,A 2,A 3,A 4和对应的5组组号数组F 0,F 1,F 2,F 3,F 4；
 （7）、依次对5组统计数组A 0,A 1,A 2,A 3,A 4和5组组号数组F 0,F 1,F 2,F 3,F 4，统计编码截断阶数p分别为1,2,3,…,32阶下的编码比特数，选取编码比特数最小时的合并次数q c 和编码截断阶数p c ；
(8）、逐行提取非负混合一阶差分矩阵Z m×n 元素，依次放入数组E中，依步骤（7）选取的2个参数，即合并次数q c 和编码截断阶数p c ，对数组E中各元素进行合并与编码，生成一次压缩编码二值码流T 1；
（9）、统计一次压缩编码二值码流T 1中连续0和连续1的长度，得到一次压缩编码二值码流T 1的长度分布数组L；
（10）、用长度分布数组L的长度商数组Q和余数二值码流R表示长度分布数组L，实现二次编码压缩，生成二次编码二值码流T 2。
2.如权利要求1所述的一种图像数据无损压缩的编码方法，其特征在于：在（2）中所述的从横向一阶差分矩阵X m×n 和纵向一阶差分矩阵Y m×n 中生成混合一阶差分矩阵XY m×n 和混合一阶差分矩阵XY m×n 子块的状态二值码流D，具体如下：
对横向一阶差分矩阵X m×n 和纵向一阶差分矩阵Y m×n 按4×4大小的子块划分，分别计算相同位置两差分矩阵子块元素绝对值之和，取绝对值之和最小的子块对应的各个元素置入混合一阶差分矩阵XY m×n 相同位置子块处，生成一阶混合差分矩阵XY m×n =(xy ij ) m×n ，对横向一阶差分矩阵X m×n 和纵向一阶差分矩阵Y m×n 相同位置子块，记 
Figure 2012100407057100001DEST_PATH_IMAGE001
，
Figure 691485DEST_PATH_IMAGE002
，若
Figure 2012100407057100001DEST_PATH_IMAGE003
，则将横向一阶差分矩阵X m×n 子块数据置入混合一阶差分矩阵XY m×n 相同位置子块处，此时子块的状态标记为0；若
Figure 644398DEST_PATH_IMAGE004
，则将纵向一阶差分矩阵Y m×n 子块数据置入混合一阶差分矩阵XY m×n 相同位置子块处，此时子块的状态标记为1，将子块的状态标记0或1依次置入子块的状态二值码流D中。
3.如权利要求1所述的一种图像数据无损压缩的编码方法，其特征在于：在（4）中所述的统计非负混合一阶差分矩阵Z m×n 中各元素出现的个数，生成元素个数统计数组A，具体如下：
统计非负混合一阶差分矩阵Z m×n =(|xy ij |) m×n 中各元素出现的个数，按元素出现次数的大小置入元素个数统计数组A=(n 0,n 1,…,n e )中，其中n 0≥n 1≥…≥n e ，n 0,n 1,…,n e 分别对应非负混合一阶差分矩阵Z m×n =(|xy ij |) m×n 中元素值a 0,a 1,a 2,…,a e 出现的个数, 0≤e≤255，其中|xy ij |
Figure 2012100407057100001DEST_PATH_IMAGE005
{a 0,a 1,a 2,…,a e }，表示非负混合一阶差分矩阵Z m×n 中任一点元素值|xy ij |仅在集合{ a 0,a 1,a 2,…,a e }中取值。
4.如权利要求1所述的一种图像数据无损压缩的编码方法，其特征在于：在（5）中所述的从元素个数统计数组A中生成元素映射数组T S ，依元素映射数组T S 对非负混合一阶差分矩阵Z m×n 各元素重新分配数值，具体如下：
依次将元素个数统计数组A=(n 0,n 1, …,n i ,…,n e )中n i 所对应的非负一阶混合差分矩阵Z m×n 中元素a i ,置入元素映射数组T S 中，生成元素映射数组T S =(a 0,a 1, …,a i ,…,a e )，逐行提取非负混合一阶差分矩阵Z m×n 各元素|xy ij |，|xy ij |为非负混合一阶差分矩阵Z m×n 中第i行第j列元素，如果|xy ij |=a k  ，一阶非负混合差分矩阵Z m×n 中第i行第j列元素值|xy ij |置为k，其中a k 是映射数组T S =(a 0,a 1,a 2,…,a e )中的第k+1个元素，0≤k≤e。
5.如权利要求1所述的一种图像数据无损压缩的编码方法，其特征在于：在步骤（6）中所述的对步骤（4）所得元素个数统计数组A分别进行q=0,1,2,3,4次数据合并，建立5组统计数组A 0,A 1,A 2,A 3,A 4和对应的5组组号数组F 0,F 1,F 2,F 3,F 4，具体如下：
为论述方便，假定e=255，则
0次合并统计数组A 0=A=(n 0,n 1,…,n 255)；   
0次合并组号数组F 0=(0,1,2,…,255)；
1次合并统计数组A 1=(n 0+n 1, n 2+n 3…, n 254+n 255)=
Figure 254502DEST_PATH_IMAGE006
；
1次合并组号数组F 1=(0,1,2,…,127)；
2次合并统计数组A 2=(n 0+n 1+ n 2+n 3,…, n 252+n 253+ n 254+n 255)=
Figure 2012100407057100001DEST_PATH_IMAGE007
；
2次合并组号数组F 2=(0,1,2,…,63)；
3次合并统计数组A 3=(n 0+n 1+…+n 7,…, n 248+n 249+…+n 255)=；
3次合并组号数组F 3=(0,1,2,…,31)；
4次合并统计数组A 4=(n 0+n 1+…+n 15,…, n 240+n 241+…+n 255)=
Figure 2012100407057100001DEST_PATH_IMAGE009
；
4次合并组号数组F 4=(0,1,2,…,15)，
其中，A为步骤（4）所述元素个数统计数组，组号数组F k 中元素a i 经k+1次合并后对应到组号数组F k+1中元素round[a i /2]，round为取整函数，0≤k≤3。
6.如权利要求1所述的一种图像数据无损压缩的编码方法，其特征在于：在上述步骤（7）中所述的依次对5组统计数组A 0,A 1,A 2,A 3,A 4和5组组号数组F 0,F 1,F 2,F 3,F 4，统计编码截断阶数p分别为1,2,3,…,32阶下的编码比特数，选取编码比特数最小时的合并次数q c 和编码截断阶数p c ，具体如下：
编码截断阶数p为可用1元即时二进制码0,10,110,…,11…10，分别表示元素0,1,2,…,p-1，其中，11…10由p-1个1和1个0组成，大于p-1的元素的二进制码，由p个1的二进制码组成其前缀，由b位二进制码组成其后缀，其中参数b由公式（6）确定，
如果编码截断阶数为p，合并次数为q，对图像矩阵D m×n 编码占用的比特数G为
  （5）
其中， B 0(D)为步骤（2）所述子块的状态二值码流D占用比特数；B 1(S)为步骤（3）所述混合一阶差分矩阵元素的符号二值码流S占用比特数；B 2(T s )为步骤（5）所述元素映射数组T s 占用比特数，
Figure 2012100407057100001DEST_PATH_IMAGE011
对应的组号元素为
Figure 109828DEST_PATH_IMAGE012
，参数b按如下公式确定：
记
Figure 2012100407057100001DEST_PATH_IMAGE013
，
如果0≤f≤1,则b=1；
如果2≤f≤3,则b=2；
如果4≤f≤7,则b=3；
如果8≤f≤15,则b=4；                                              （6）
如果16≤f≤31,则b=5；
如果32≤f≤63,则b=6；
如果64≤f≤127,则b=7；
如果128≤f≤255,则b=8，
对0,1,2,3,4次合并后的5组统计列表A 0,A 1,A 2,A 3 ,A 4和5组组号数组F 0,F 1,F 2,F 3,F 4，分别统计编码截断阶数为p=1,2, …,32阶下占用的编码比特数，由5种合并，32个编码截断阶数，共计得5×32=160种组合形式，按公式（5）分别计算160种组合的编码占用的比特数G，从160种组合中选取编码占用的比特数G最小时对应的编码截断阶数p=p c 和元素合并次数q=q c ，其中1≤p c ≤32，0≤q c ≤4。
7.如权利要求1所述的一种图像数据无损压缩的编码方法，其特征在于：在上述步骤（8）所述的逐行提取非负混合一阶差分矩阵Z m×n 元素，依次放入数组E中，依步骤（7）选取的2个参数，即合并次数q c 和编码截断阶数p c ，对数组E中各元素进行合并与编码，生成一次压缩编码二值码流T 1,具体如下：
逐行提取矩阵Z m×n =(|xy ij |) m×n 元素，依次放入数组E中，E的表达式为
E=(|xy 11|,|xy 12|,…, |xy 1n |,|xy 21|,|xy 22|,…, |xy 2n |,…, |xy  m1|,|xy  m2|,…, |xy  mn |)，
按步骤（7）所选合并次数q c ，对数组E中各元素整除
Figure 380580DEST_PATH_IMAGE014
，商置入数组，余数置入数组
Figure 539028DEST_PATH_IMAGE016
，中各余数的二进制表示位数为q c  位，对数组
Figure 2012100407057100001DEST_PATH_IMAGE017
中各元素进行编码，小于p c 的元素0,1,2,…,p c -1对应的编码为
Figure 115820DEST_PATH_IMAGE018
；大于p c 的元素对应的编码为
Figure 2012100407057100001DEST_PATH_IMAGE019
，其中参数b的取值由步骤（7）所述公式（6）确定，
将各元素编码的前缀部分，即1元即时码0,10,110,…,11…10,11…11，依数组
Figure 67727DEST_PATH_IMAGE020
中元素顺序进行纵向分解，按位平面层面依次置入一次压缩编码二值码流T 1中，接着在一次压缩编码二值码流T 1中置入后缀部分的b位二进制数，然后在一次压缩编码二值码流T 1中置入余数数组
Figure 2012100407057100001DEST_PATH_IMAGE021
对应的二值码流、步骤（3）所述混合一阶差分矩阵XY m×n 的符号二值码流S和步骤（2）所述混合一阶差分矩阵XY m×n 子块的状态二值码流D，生成一次压缩编码二值码流T 1。
8.如权利要求1所述的一种图像数据无损压缩的编码方法，其特征在于：在上述步骤（9）中所述的统计一次压缩编码二值码流T 1中连续0和连续1的长度，得到一次压缩编码二值码流T 1的长度分布数组L，具体如下：
统计步骤（8）一次压缩编码二值码流T 1中0段和1段的长度，得到长度分布数组L=(l 1,l 2,l 3,…,l r)。
9.如权利要求1所述的一种图像数据无损压缩的编码方法，其特征在于：在上述步骤（10）中所述的用长度分布数组L的长度商数组Q和余数二值码流R表示长度分布数组L，实现二次编码压缩，生成二次编码二值码流T 2,具体如下：
对长度分布数组L=(l 1,l 2,…,l r )中各元素按公式（7）计算，得到长度商数组Q和余数二值码流R，其中
Q=round((l 1+1)/2,(l 2+1)/2,…,(l r +1)/2)                            （7）
round表示取整函数，余数置入余数二值码流R=(b 1,b 2,…, b i ,…,b r )中，其中，当(l i +1)/2为偶数时，b i 为0，当(l i +1)/2为奇数时，b i 为1，
长度商数组Q中元素为1时，编码为0或1，如果该原始长度由0组成，则编码为0，如果该原始长度由1组成，则编码为1；为2时，编码为00或11，如果该原始长度由0组成，则编码为00，如果该原始长度由1组成，则编码为11；以此类推，如某二值码流串0111111000111111111111，其长度分布数组L=(1,6,3,12)，长度商数组Q编码为011100111111，余数二值码流R为0101，
将长度商数组Q的二进制编码与余数二值码流R组合后，生成二次编码二值码流T 2。
CN201210040705.7A 2012-02-22 2012-02-22 一种图像数据无损压缩的编码方法 Expired - Fee Related CN102547295B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210040705.7A CN102547295B (zh) 2012-02-22 2012-02-22 一种图像数据无损压缩的编码方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210040705.7A CN102547295B (zh) 2012-02-22 2012-02-22 一种图像数据无损压缩的编码方法

Publications (2)

Publication Number Publication Date
CN102547295A true CN102547295A (zh) 2012-07-04
CN102547295B CN102547295B (zh) 2014-07-30

Family

ID=46353093

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210040705.7A Expired - Fee Related CN102547295B (zh) 2012-02-22 2012-02-22 一种图像数据无损压缩的编码方法

Country Status (1)

Country Link
CN (1) CN102547295B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108810553A (zh) * 2018-05-02 2018-11-13 上海大学 一种基于稀疏化处理的移动节点监测数据序列压缩方法
WO2023051335A1 (zh) * 2021-09-30 2023-04-06 华为技术有限公司 数据编码方法、数据解码方法以及数据处理装置
CN115955571A (zh) * 2023-03-10 2023-04-11 深圳市启明智显科技有限公司 一种用于嵌入式设备的图像储存方法、装置及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101126807A (zh) * 2007-09-20 2008-02-20 西安电子科技大学 遥感图像二倍无损及近无损压缩编码方法
US20080292028A1 (en) * 2005-10-31 2008-11-27 Lg Electronics, Inc. Method and Apparatus for Signal Processing and Encoding and Decoding Method, and Apparatus Therefor
CN101329725A (zh) * 2008-07-30 2008-12-24 电子科技大学 基于梯度投影和形态学的指纹图像分割方法
CN101984463A (zh) * 2010-11-02 2011-03-09 中兴通讯股份有限公司 全景图合成方法及装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080292028A1 (en) * 2005-10-31 2008-11-27 Lg Electronics, Inc. Method and Apparatus for Signal Processing and Encoding and Decoding Method, and Apparatus Therefor
CN101126807A (zh) * 2007-09-20 2008-02-20 西安电子科技大学 遥感图像二倍无损及近无损压缩编码方法
CN101329725A (zh) * 2008-07-30 2008-12-24 电子科技大学 基于梯度投影和形态学的指纹图像分割方法
CN101984463A (zh) * 2010-11-02 2011-03-09 中兴通讯股份有限公司 全景图合成方法及装置

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108810553A (zh) * 2018-05-02 2018-11-13 上海大学 一种基于稀疏化处理的移动节点监测数据序列压缩方法
CN108810553B (zh) * 2018-05-02 2021-01-12 上海大学 一种基于稀疏化处理的移动节点监测数据序列压缩方法
WO2023051335A1 (zh) * 2021-09-30 2023-04-06 华为技术有限公司 数据编码方法、数据解码方法以及数据处理装置
CN115955571A (zh) * 2023-03-10 2023-04-11 深圳市启明智显科技有限公司 一种用于嵌入式设备的图像储存方法、装置及系统

Also Published As

Publication number Publication date
CN102547295B (zh) 2014-07-30

Similar Documents

Publication Publication Date Title
CN100596029C (zh) Ldpc码校验矩阵构造方法及利用该方法的编码解码装置
CN101442676B (zh) 基于cavlc编码的h.264视频加密方法
CN100517979C (zh) 一种数据压缩及解压缩方法
Wong et al. Simultaneous arithmetic coding and encryption using chaotic maps
CN102396219A (zh) 对于具有稀疏直方图的图像的嵌入式图形编码
CN107018419A (zh) 一种基于ambtc的图像压缩编码方法
CN102547295B (zh) 一种图像数据无损压缩的编码方法
CN100423582C (zh) 一种用于将待编码数据进行二进制化编码的方法和装置
EP2200182A4 (en) METHOD AND DEVICE FOR ENCODING THE LOW DENSITY GENERATOR MATRIX CODE
CN101661631A (zh) 基于局部最小生成树的点模型快速无损压缩处理方法
CN102263560B (zh) 差分编码方法和系统
CN105005806B (zh) 一种矩阵点图标识、生成方法、生成系统及防伪标签
CN102651795B (zh) 游长缩减的二元序列压缩编码方法
CN104093027A (zh) 用于彩色图像的联合标量嵌入式图形编码
Zhang et al. Generalization of the ZZW embedding construction for steganography
CN103973310A (zh) 一种基于异或逻辑运算折半划分的测试数据压缩方法
Mukesh et al. Enhancing AES algorithm with arithmetic coding
CN101657973A (zh) 具有采用位精度进行编码和解码的程序的记录介质及其装置
KR101578608B1 (ko) 일반 코드 생성을 위한 데이터 부호화 장치 및 데이터 복호화 장치
CN104702404B (zh) 一种基于纠错码与低轮aes的加密方法及系统
Sachnev et al. An improved matrix encoding scheme for JPEG steganography
Medeni et al. Steganographic algorithm based on error-correcting codes for gray scale images
CN102377996A (zh) 上下文自适应可变长编码方法及装置
Chong et al. Stepping-random code: a rateless erasure code for short-length messages
CN101436865B (zh) 对经过Turbo编码后的比特流进行打孔的方法及装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent of invention or patent application
CB03 Change of inventor or designer information

Inventor after: Gao Jian

Inventor after: Rao Jun

Inventor after: Liu Wan

Inventor after: Song Ao

Inventor after: Wang Zhongyuan

Inventor after: Shu Qiwei

Inventor before: Gao Jian

Inventor before: Rao Jun

Inventor before: Liu Wan

Inventor before: Song Ao

Inventor before: Wang Zhongyuan

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: GAO JIAN RAO JUN LIU WAN SONG AO WANG ZHONGYUAN TO: GAO JIAN RAO JUN LIU WAN SONG AO WANG ZHONGYUAN SHU QIWEI

C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140730

Termination date: 20170222