专利名称：基于空间预测与变换的遥感图像压缩方法
技术领域：
本发明涉及一种遥感图像压縮的方法。这是一种基于空间预测与变换的针对细节丰富的 遥感图像的压縮方法。
背景技术：
随着航天和遥感技术的飞速发展，遥感图像在日常生活、工农业生产和国防等领域受到 越来越广泛的应用。由于星载传感器的分辨率不断提高，其图像数据采集速率己经远远超过 了传输信道的吞吐速率，所以必须对遥感图像进行高保真的压縮。近年来，图像数据的压縮 取得了很大的发展，涌现了很多非常优秀的压縮算法。目前受到广泛应用，并己纳入新的国 际标准的基于小波的压縮方法是一种成熟的、具有高保真压縮效果的压縮方法，但针对遥感 图像的压縮，该方法有一些严重的缺陷，比如，图像中经小波分解压縮后重建的"锐边"，不 可避免的被"钝化"(即损失了高频信息)，并且随着压縮倍率的提高，"钝化"将越严重，直接 造成图像中纹理细节的损失；该方法必须以大面积图像为基础(随着小波分解所用基础面积 的减小，恢复图像质量将急剧降低)，并在多次二维小波分解基础上再做无失真编码后，才能 得到高保真的压縮效果，但这种以大面积作为处理基础的算法是难于实现实时压縮的；基于 小波的压缩方法是通过牺牲高频信息、保存低频信息的方法来换取压缩比的，这对于遥感图 像来说，因为其特殊的应用背景，往往小目标才是它关心的重要对象。而图像中的小目标与 其所在的局部环境的对比度较小，经过小波变换后，其信息大多数集中在第一级小波分解的 三个高频分量中，很容易被量化掉。所以对遥感图像是不利的。
基于空间重采样的遥感图像压缩算法在压缩比为4倍左右的应用情况下，其压縮质量不 低于基于小波的压縮算法。但是因为基于空间重采样的遥感图像压縮算法在对图像的结构描 述上过于简单，随着压縮比的上升，其压縮质量迅速下降；基于结构描述的图像压縮算法以 4x4的像素块作为基本的图像子区，用作压縮处理的基本单元，每个图像子区根据其纹理结 构的复杂程度，分为平坦区、粗纹理区和细纹理区三类，以子区内图像的平均方差作为分类 的准则，其压縮质量较基于空间重采样的遥感图像压縮算法有显著提高。
虽然基于结构描述的图像压縮算法考虑了图像中纹理的分布，进行分区压縮，但是该算 法只是用了4x4大小的块，并且只进行了三种分类模式，显然，这种方法并不能覆盖图像中
的所有纹理。
发明内容
本发明要解决的技术问题是针对现有技术的上述不足，延续了基于结构描述的图像压 缩算法的思想，借鉴了 &264视频编码标准中的帧内编码技术，提供一种对细节纹理丰富的
遥感图像的压縮方法，提高压縮效率。
为实现这样的方法，本发明的技术方案中，以块为单位进行处理。对每一块而言，先利 用相邻的已编码解码块对当前块进行预测，得到预测块后，用当前块减去预测块得到残差块，
这样处理能够降低相邻块的相关性，提高编码效率。接下来对残差块进行整数DCT变换和量 化。之所以采用DCT变换，这是因为DCT被认为是对K-L变换的最理想的逼近，它能够使 图像能量经变换后尽可能集中于尽可能少的几个系数上，从而得到较高的压縮比。根据有关 研究，DCT变换在高信噪比的情况下比小波变换具有更好的压縮性能；最后，对经过变换编 码的数据进行熵编码。
本发明提供一种基于空间预测与变换的遥感图像压縮方法，其特征在于包括如下具体步
骤
(1) 对图像进行分块把图像分成若干个16x16大小的宏块，如果图像的行或列不能 被16整除，则对图像进行补齐，把图像扩展，使其行或列被16整除；补齐的方法是如果 行不能被16整除，则对最后一行进行复制，然后填充扩展的行；如果列不能被16整除，则 对最后一列进行复制，然后填充扩展的列；
(2) 预测分析对一个16x16的块进行预测分析，先利用相邻的已编码解码块对当前 块进行预测，得到预测块后，用当前块减去预测块得到残差块；接下来对残差块进行整数DCT 变换和量化；经过上面的分析，得到最终进行编码的块大小16x16或者4x4，以及最佳预测 模式。
(3) 预测编码根据B得到的块大小及最佳预测模式对该块进行变换、量化；此外， 还要进行变换和量化的块进行反量化和反变换的处理，然后与预测值一起生成重建块；
(4) 对C中得到的数据进行zigzag扫描，对预测编码后的数据进行重新排列；然后进 行熵编码；
(5) 重复B、 C、 D，遍历图像中所有的16x16宏块。
图1是基于空间预测和变换的遥感图像压縮方法流程图； 图2是以绝对误差和为代价函数的模式选择流程图。
具体实施例方式
为了更好的理解本发明的技术方案，以下结合附图对本发明的实施方式作进一步描述。 本发明基于空间预测与变换的遥感图像压缩方法流程如图1所示，各部分具体实施细节 如下
1.对图像进行分块；
把图像分成若干个16x16大小的宏块，如果图像的行或列不能被16整除，则对图像进行 补齐，把图像扩展，使其行或列被16整除；补齐的方法为如果行不能被16整除，则对最 后一行进行复制，然后填充扩展的行；如果列不能被16整除，则对最后一列进行复制，然后 填充扩展的列；
2.对待编码块进行预测分析；
遥感图像不同于一般的自然图像，它所含的高频信息和纹理区域较多，所以本实施例采 用16x16和4x4两种预测类型对待编码块进行预测分析，在选择编码预测方式时以绝对误差 和SAD (Sum of Absolute Difference)作为判断准则。
绝对误差和是指当前待编码块与预测块之间的差值绝对值之和，即残差之和。假设当前 块用S(z'，/)表示，预测值为pm/(/J)， i，j'-0，l，…,iV，则
<formula>formula see original document page 5</formula>
绝对误差和(SAD)值越小，说明图像的残差值越小，从而说明当前待编码块与预测块 越接近，相应预测模式的预测的精度越高；反之，绝对误差和(SAD)值越大，说明图像的 残差值越大，从而预测块不能很好的反应当前待编码块，相应预测模式的预测精度越低。
为了综合考虑各种预测模式的不同偏重方向性，当对宏块进行编码时会遍历所有的预测 模式，并借助模式判决准则从中选出最优的模式作为最终的编码模式。
图2是以SAD代价函数模型进行模式选择的流程-
(1) 根据任一种预测模式构造出相应块的预测值P;
(2) 计算原始块与预测块P之间的绝对误差和&4Z);
(3) 对于16xl6的块来说，每种预测模式下的绝对误差和SAD值即是该模式下的代价 cost一16,即cost一16-SAD;对于4x4的块来说，每种预测模式下的代价cost_4的定义如下
cosr_4 = &4D + A*i a&
其中，；L = 0.85*2(w—12)",《p表示量化参数；Rate表示采用这种预测模式需要的比特数。
(4) 对于16x16的块来说，按照步骤(1)、步骤(2)、步骤(3)计算每种可能预测模
式的cost_16值，选取cost_16值最小的作为最佳的16x16预测模式；
(5) 对于4x4的块来说，按照步骤(1)、步骤(2)、步骤(3)计算每种可能预测模式 的cost—4值，选取cost_4值最小的作为该块最佳的预测模式；
(6) 按照步骤(1)、步骤(2)、步骤(3)、步骤(4)分别找出16个4x4块的最佳模 式及其最小的代价值cost_4，并把16个cost—4相加得到cost4—sum;
(7) 对于两种预测类型，根据上面的模式选择结果，选取costJ6和cost4一sum中最小 的那种作为最佳预测类型，这种预测类型下的预测模式就是最佳预测模式。
3. 对预测残差值进行变换编码和量化；
对残差值进行4x4的DCT变换和量化，所采用的变换和量化技术是H.264标准中的技术。
整数DCT对图像信号作用后可以获得与DCT类似的频率分解；同时，正交变换的目的 是减少变换系数的相关性，使得变换系数的能量集中在较少的几个分量上，变换前后的能量 并没有损失，在去除空间相关性上，整数DCT与DCT具有相同的效果。
在本实施例中还用到Hadamard变换，Hadamard变换是针对16x16的块经过整数变换后 的16个直流DC系数的变换，经过整数4x4DCT变换后，由于能量大多数集中在每个块的 DC系数中，因此我们利用四阶的Hadamard变换来去除宏块中4x4直流系数的+&关性。
在本实施例中为了进行待编码块的预测，要对已作过整数DCT和量化的块进行反量化和 整数DCT的反变换。
结合附图2来看，预测值加上反变换之后的残差值就是重建块，即重建块-尸+ D:，其
中Z)二表示经过反量化和反变换后的残差值。重建编码块的目的是为了对相邻未编码块进行预 测。
4. 对经过变换编码和量化处理后得图像进行扫描和熵编码；
对经过预测编码后的数据进行之字形(zigzag)扫描，然后用基于上下文的自适应二进制 算术编码(CABAC)进行熵编码。
CABAC的编码过程由下面三个基本步骤组成
(1) 二进制转换；
(2) 上下文建模，编码的自适应性就在这一步体现；
(3) 二进制算术编码。
5. 最后对所有的16xl6宏块进行2、 3、 4的操作步骤。
权利要求
1、一种基于空间预测与变换的遥感图像压缩方法，其特征在于包括如下具体步骤(1)对图像进行分块把图像分成若干个16×16大小的宏块，如果图像的行或列不能被16整除，则对图像进行补齐，把图像扩展，使其行或列被16整除；(2)预测分析对一个16×16的块进行预测分析，先利用相邻的已编码解码块对当前块进行预测，得到预测块后，用当前块减去预测块得到残差块；接下来对残差块进行整数DCT变换和量化；经过上面的分析，得到最终进行编码的块大小16×16或者4×4，以及最佳预测模式；(3)预测编码根据步骤(2)得到的块大小及最佳预测模式对该块进行变换、量化，此外，还要进行变换和量化的块进行反量化和反变换的处理，然后与预测值一起生成重建块；(4)对步骤(3)中得到的数据进行之字形扫描，对预测编码后的数据进行重新排列，然后进行熵编码；(5)重复步骤(2)、步骤(3)、步骤(4)，遍历图像中所有的16×16宏块。
2、 一种基于空间预测与变换的遥感图像压缩方法，其特征在于所述步骤(I)中的补 齐的方法是如果行不能被16整除，则对最后一行进行复制，然后填充扩展的行；如果列不 能被16整除，则对最后一列进行复制，然后填充扩展的列。
全文摘要
本发明提出一种基于空间预测与变换的遥感图像压缩方法，该方法充分利用当前块与相邻块的空间信息，根据相邻已编码重建的块来预测当前的像素值，对预测残差进行变换编码进一步消除空间冗余，在变换方面，根据不同子块的特点自适应调整编码子区的大小，使用整数4×4离散余弦变换；同时，对16×16宏块经过整数4×4离散余弦变换的直流系数进行Hadamard变换，进一步改善大面积平坦图像区的编码性能；本发明在空间域进行预测，避免了频率域预测时丢失细小细节的缺点，提高了算法的压缩性能。
文档编号H04N7/26GK101383972SQ20081022350
公开日2009年3月11日 申请日期2008年9月28日 优先权日2008年9月28日
发明者周付根, 王艳芹 申请人:北京航空航天大学