基于三维矩阵wdct变换的彩色半调图像压缩方法
【专利摘要】本发明公开了一种基于三维矩阵WDCT变换的彩色半调图像压缩方法。其步骤：（1）对输入的彩色半调图像进行颜色空间转换；（2）三维矩阵WDCT正变换；（3）统计并建立长游程系数位置扫描规律；（4）系数重排；（5）游程Huffman（RLH）编码；（6）输出码流；（7）输入码流；（8）RLH解码；（9）系数恢复；（10）三维矩阵WDCT反变换；（11）颜色空间反变换；（12）输出解压图。本发明考虑到彩色半调图像在空间结构和通道间的信息冗余和相关性，把彩色图像的多个通道处理统一起来，对其进行WDCT变换，既去除了各通道间的相关性，也使得图像能量分布相对集中；本发明中的长游程系数扫描规律不受图像内容，及其颜色表示模型、彩色半调技术限制，具有广泛适用性。
【专利说明】基于三维矩阵WDCT变换的彩色半调图像压缩方法
【技术领域】
[0001]本发明属于图像处理【技术领域】的图像压缩方法，具体涉及一种彩色半调图像经颜色空间转换及WDCT变换后，统计系数长游程扫描规律再进行压缩的方法，可支持远程打印、传真过程中彩色半调图像的压缩。
技术背景
[0002]在彩色半调图像远程打印、传真过程中，为减少图像存档所需的存储空间、缩短传输时间，有效解决传输效率和有限带宽之间的矛盾，需要对彩色半调图像进行压缩。
[0003]目前彩色半调图像的压缩技术主要有分颜色通道基于模板的压缩算法和矢量压缩算法。2008年，朴成范等人在专利“对半色调图像编码和解码的方法和设备” [P]中提出一种针对CMYK颜色空间下的点聚簇类有序抖动彩色半调图像的压缩方法。该方法针对彩色半调图像各颜色通道的特性设计了相应的模板，并采用分通道方式进行算术编码，解压时需附带模板信息。2009年，朴成范等人又在文献“Sung-Bum Park(朴成范)，Woo-SungShim, et al.Novel context template design scheme for lossless color halftoneimage compression[C].Proc.0f SPIE-1S&T Electronic Imaging,，，中改进了之前的模板。算法主要思想是对每个像素的相邻区域进行线性采样，并用线性采样特征和两个约束条件构成模板，然后分通道进行算术编码，编码效率较JBIG标准有所提高。以上这两种方法计算复杂度高，没有去除彩色半调图像各通道间的相关性，且只适用于CMYK颜色空间下的点聚簇类有序抖动彩色半调图像。2012年，曹慧敏在文献“曹慧敏.灰度和彩色半调图像无损压缩算法研究[D].西安:西安建筑科技大学，2012”中提出将彩色半调图像作为一个整体进行压缩的方法。其基本思想是将彩色半调图像的RGB三分量作为一个矢量，每八个连续的矢量作为一个编码单元，根据R、G、B三个分量中的元素取值情况将待编码矢量块划分为纯色块和混色块，最终采用Huffman编码法进行压缩编码。这种压缩方法与图像内容和半调类型关系密切，图像内容和半调类型不同，则压缩比差异较大，尤其是它应用于误差分散类半调图像时，甚至会出现数据膨胀现象。
[0004]综上所述，朴成范等人提出的分通道彩色半调图像压缩方法虽然能实现数据压缩，但未考虑通道间的相关性，其本质还是灰度半调图像压缩方法的扩展，且仅适用于点聚簇类有序抖动彩色半调图像，适用范围窄。曹慧敏提出的矢量化彩色半调图像压缩算法复杂度低，考虑了彩色通道间的相关性，但未考虑消减像素间的相关性，且只能对有序抖动类彩色半调图像进行有限的压缩，对误差分散类彩色半调图像压缩时，还可能出现数据膨胀现象，达不到图像压缩的目的。
【发明内容】
[0005]本发明针对现有彩色半调图像压缩方法的不足，提出一种基于三维矩阵WDCT变换的彩色半调图像压缩新方法，它通过三维矩阵WDCT变换，既消除了通道间的相关性，又消减了像素间的冗余，在此基础上统计并生成了变换后系数的长游程位置扫描规律，将系数相同的像素排在一起，增大了游程长度，使RLH编码方法可以扬其所长，数据压缩比大大提高，且算法复杂度低。此外，本发明不受图像内容及其颜色表示空间、半调方法的限制，适用范围广。
[0006]为了实现上述技术任务，本发明采用如下技术方案予以实现:
[0007](I)对输入的彩色半调图像进行颜色空间转换:将输入的彩色半调图像转换到YCbCr颜色空间，其中，亮度分量Y、色差分量Cb、Cr的取值范围均为[0，255]；
[0008](2 )对YCbCr颜色空间的图像进行WDCT正变换:
[0009]①对YCbCr颜色空间的图像做分块；
[0010]②对分块子图集做WDCT正变换，得到变换系数并采用JPEG量化表进行量化；
[0011](3)建立各通道的WDCT系数长游程位置扫描规律:对量化后的WDCT系数进行分通道统计，建立各通道量化系数的长游程位置扫描规律；
[0012](4)对量化后的WDCT系数重新排序:根据建立的各通道长游程系数扫描规律，分通道重排各子图集的WDCT系数；
[0013](5)对重排后的系数做RLH编码:
[0014]①对重排后的系数做游程统计；
[0015]②以游程信息作为Huffman编码的信源，进行RLH编码，形成图像的编码数据；
[0016](6)构造压缩数据形成输出码流:
[0017]①将图像分块的尺度参数作为压缩信息的第一部分数据，为其分配一个字节的码流；
[0018]②将块内的各通道长游程位置扫描规律描述参数作为压缩信息的第二部分数据，为其分配的存储空间大小由分块的尺度参数、颜色通道数确定；
[0019]③最后将图像的RLH编码数据作为压缩信息的第三部分数据，合成输出码流；
[0020](7)接收输入码流，压缩数据拆分:
[0021]①取出首字节码流信息，作为图像分块的尺度参数；
[0022]②根据分块的尺度参数、颜色通道数取出相应字节的数据，作为码流信息的第二部分数据，它是块内的各通道长游程位置扫描规律描述参数；
[0023]③剩余的码流信息为图像的编码数据；
[0024](8)对图像编码数据解码:
[0025]①对图像编码数据进行Huffman解码；
[0026]②对Huffman解码数据进行游程解码；
[0027](9)依据各通道的系数位置扫描规律对解码后的各个系数进行位置恢复:依据解码出的块内各通道长游程位置扫描规律描述参数，分通道恢复解码后系数的位置。
[0028](10)对恢复位置后的系数做WDCT反变换:
[0029]①对恢复位置后的系数做反量化；
[0030]②对反量化系数做WDCT反变换；
[0031]③将反变换后的子图集进行重组；
[0032](11)颜色空间反变换:
[0033]①将超出亮度分量Y和色差分量Cb、Cr取值范围的数据进行相应的映射；
[0034]②将YCbCr颜色空间的图像转换到输入图像颜色空间；[0035]( 12)输出解压图:将颜色空间反变换结果进行二值化处理并输出。
[0036]上述步骤(2)进行三维矩阵WDCT正变换的步骤如下:
[0037]第一，将YCbCr颜色空间的图像按块不重叠的原则分割成长方体子图块，将长方体子图块依次向后排列，形成三维长矩阵，将三维长矩阵分割成正方体子图集；
[0038]第二，对正方体子图集做WDCT横向、竖向、纵向正变换；
[0039]第三，采用JPEG量化表对变换后的系数集做量化处理；
[0040]上述步骤(3)建立各通道的WDCT系数长游程位置扫描规律的步骤如下:
[0041]第一，将量化后系数集按颜色通道抽取，形成分通道系数子集；
[0042]第二，分通道统计系数子集中各个位置上出现次数最多的系数，计算其出现概率，记录其位置信息；
[0043]第三，分通道将系数、概率、位置信息按照系数优先、概率次优的顺序做降序排列，排序后不同通道的对应位置信息就是所在通道中WDCT系数的长游程位置扫描规律。
[0044]与现有技术相比，本发明具有以下优点:
[0045]第一，本发明充分考虑通道间的相关性和像素间的冗余性，利用WDCT变换消除通道间的相关性、消减同一 通道内像素间的冗余，使大量分散的图像相关信息集中起来，在变换域对彩色半调图像进行压缩，新方法的压缩性能明显超过现有方法，也丰富了彩色半调图像压缩方法的种类。
[0046]第二，本发明构造了三维矩阵WDCT系数的长游程位置分布规律统计方法，它可以自适应地总结出三维矩阵WDCT系数长游程扫描规律，不受图像内容及其颜色表示模型、半调方法的影响，具有良好的普适性。
【专利附图】
【附图说明】
[0047]图1为本发明流程图。
[0048]图2为本发明处理过程中以512X512图像为例的变化情况示意图。
[0049]图3为本发明对几种典型彩色半调图像的实验结果图。
【具体实施方式】
[0050]参照图1，以RGB彩色半调图像为例，用本发明压缩的【具体实施方式】如下。
[0051]步骤(1).对输入的彩色半调图像进行颜色空间转换:
[0052]①将输入的RGB彩色半调图像按照JPEG标准转换到YCbCr颜色空间，其中，亮度分量Y、色差分量Cb、Cr的取值范围均为[0，255]，转换公式如下；
【权利要求】
1.一种基于三维矩阵WDCT变换的彩色半调图像压缩方法，包括以下步骤: (O对输入的彩色半调图像进行颜色空间转换:将输入的彩色半调图像转换到YCbCr颜色空间，其中，亮度分量Y、色差分量Cb、Cr的取值范围均为[O，255]； (2)对YCbCr颜色空间的图像进行WDCT正变换: ①对YCbCr颜色空间的图像做分块； ②对分块子图集做WDCT正变换，得到变换系数并采用JPEG量化表进行量化； (3)建立各通道的WDCT系数长游程位置扫描规律:对量化后的WDCT系数进行分通道统计，建立各通道量化系数的长游程位置扫描规律； (4)对量化后的WDCT系数重新排序:根据建立的各通道长游程系数扫描规律，分通道重排各子图集的WDCT系数； (5)对重排后的系数做RLH编码: ①对重排后的系数做游程统计； ②以游程信息作为Huffman编码的信源，进行RLH编码，形成图像的编码数据； (6)构造压缩数据形成输出码流: ①将图像分块的尺度参数作为压缩信息的第一部分数据，为其分配一个字节的码流； ②将块内的各通道长游程位置扫描规律描述参数作为压缩信息的第二部分数据，为其分配的存储空间大小由分块的尺度参数、颜色通道数确定； ③最后将图像的RLH编码数据作为压缩信息的第三部分数据，合成输出码流； (7)接收输入码流，压缩数据拆分: ①取出首字节码流信息，作为图像分块的尺度参数； ②根据分块的尺度参数、颜色通道数取出相应字节的数据作为码流信息的第二部分数据，它是块内的各通道长游程位置扫描规律描述参数； ③剩余的压缩信息为图像的编码数据； (8)对图像编码数据解码: ①对图像编码数据进行Huffman解码； ②对Huffman解码数据进行游程解码； (9)依据各通道的系数位置扫描规律对解码后的各个系数进行位置恢复:依据解码出的块内各通道长游程位置扫描规律描述参数，分通道恢复解码后系数的位置。 (10)对恢复位置后的系数做WDCT反变换: ①对恢复位置后的系数做反量化； ②对反量化系数做WDCT反变换； ③将反变换后的子图集进行重组； (11)颜色空间反变换: ①将超出亮度分量Y和色差分量Cb、Cr取值范围的数据进行相应的映射； ②将YCbCr颜色空间的图像转换到输入图像颜色空间； (12)输出解压图:将颜色空间反变换结果进行二值化处理并输出。
2.根据权利要求1所述的基于WDCT变换的彩色半调图像压缩方法，其特征在于:所述步骤(2)进行WDCT正变换的步骤是: 第一，将YCbCr颜色空间的图像按块不重叠的原则分割成长方体子图块，将长方体子图块依次向后排列，形成三维长矩阵，将三维长矩阵分割成正方体子图集； 第二，对正方体子图集做WDCT横向、竖向、纵向变换； 第三，采用JPEG量化表对变换后的系数集做量化处理。
3.根据权利要求1所述的基于三维矩阵WDCT变换的彩色半调图像压缩方法，其特征在于:所述步骤(3)建立各通道的WDCT系数长游程位置扫描规律的步骤是: 第一，将量化后系数集按通道排列； 第二，分通道统计系数集中每个位置上出现次数最多的系数，计算其出现概率，记录其位置信息； 第三，分通道将系数、概率、位置信息按照系数优先、概率次优的顺序做降序排列，排序后不同通道的对应位置 信息就是所在通道中WDCT系数的长游程位置扫描规律。
【文档编号】H04N19/93GK103763566SQ201410007576
【公开日】2014年4月30日 申请日期:2014年1月7日 优先权日:2014年1月7日
【发明者】孔月萍, 高黄玮, 路婷婷, 张晓杰, 代冰辉, 雷红涛 申请人:西安建筑科技大学