[发明专利]基于分数傅立叶变换的图像压缩方法

说明书

(一)技术领域

本发明涉及的是一种数据压缩方法，特别是涉及一种数字图像压缩方法。 

(二)背景技术

随着计算机与数字通信技术的迅速发展，图像压缩受到了人们越来越多的关注。从本质上讲，图像压缩就是用尽可能少的代码(符号)来表示尽可能多的数据信息。现有的数字图像压缩技术可分为两大类：有损压缩和无损压缩。其中有损压缩主要包括行程长度编码、熵编码以及LZW方法；而有损压缩方法则包括色度抽样法、变换编码方法以及分形编码方法等。变换编码方法是最常用的图像压缩方法，其主要思想是使用如离散余弦变换、傅立叶变换、或者小波变换来表达图像，然后再采用一种无损编码方法实现压缩。 

分数傅立叶变换是傅立叶变换的推广形式，它将信号表示为chirp信号的组合，克服了傅立叶变换将时域信息完全摒弃的缺点，通过一种简单的方式实现了信号的从纯粹时间域到纯频率域的全过程的综合描述，能够展示出信号从纯时间域到纯频率域的所有变化特征。研究表明，通过分数傅立叶变换，可以将信号表示为紧致形式，即信号的大部分能量都集中在少数几个分数傅立叶系数上，例如chirp率为2β的chirp信号，其-π/2atan(0.5×1/β)阶分数傅立叶变换为脉冲信号。 

与本发明相关的公开报道有：1、C.Vijaya，Bhat，J.S.Signal compressionusing discrete fractional Fourier transform and set partitioning inhierarchical tree.Signal Processing.2006，86：1976-1983；2、I.S.Yetik，Kutay，M.A.，Ozaktas，H.M.Image representation and compression with thefractional Fourier transform.Optics Communications，2001，197：275-278等。 

(三)发明内容

本发明的目的在于提供一种可以将信号表示为紧致形式的特性，实现快速数字图像压缩的基于分数傅立叶变换的图像压缩方法。 

本发明的目的是这样实现的： 

第一步，计算图像的(0.5，0)、(0，0.5)、(1，0)、(0，1)阶分数傅立叶变换；第二步，根据Wigner分布与分数傅立叶变换的关系计算所有变换阶次的分数傅立叶二阶矩；第三步，根据各二阶矩的值找到最优的变换阶数，并将图像变换到最优的分数傅立叶域，从而使图像的能量集中在尽可能少的分数傅立叶系数上；第四步，忽略所有幅值小于某一阈值的系数，并采用熵编码方法对其余的分数傅立叶系数进行编码，实现图像压缩。 

本发明的有益效果是，1)计算速度快，采用快速离散分数傅立叶变换算法，本方法的时间复杂度仅为O(N×NlogN)；2)对平稳信号和非平稳的二维信号都具有较好的压缩效果。 

(四)具体实施方式

下面举例对本发明做更详细地描述： 

设I(m，n)为N×N输入图像，a，b为行、列两个方向的分数傅立叶变换阶数，并且0≤a，b≤1，Δa、Δb分别为a、b的增量。I_a，b(m，n)为I(m，n)的(a，b)阶分数傅立叶变换，w^r_a，b和w^c_a，b为I_a，b(m，n)的行向和列向二阶矩。 

本发明所提出的图像压缩算法具体实现步骤如下： 

步骤1，利用快速离散分数傅立叶变换算法计算I_0.5，0(m，n)、I_0，0.5(m，n)、I_1，0(m，n)、I_0，1(m，n)。 

步骤2，根据式(1)计算w^r_0，0，w^r_0.5，0和w^r_1，0，根据式(2)计算w^c_0，0，w^c_0，0.5和w^c_0，1