基于分数傅立叶变换的图像压缩方法.pdf

上传人:Y94****206 文档编号:971022 上传时间:2018-03-22 格式:PDF 页数:4 大小:318.48KB
返回 下载 相关 举报
摘要
申请专利号:

CN200910072925.6

申请日:

2009.09.18

公开号:

CN101944235A

公开日:

2011.01.12

当前法律状态:

授权

有效性:

有权

法律详情:

授权|||实质审查的生效IPC(主分类):G06T 9/00申请日:20090918|||公开

IPC分类号:

G06T9/00; G06F17/14

主分类号:

G06T9/00

申请人:

哈尔滨工程大学

发明人:

郑丽颖; 石大明; 田凯

地址:

150001 黑龙江省哈尔滨市南岗区南通大街145号1号楼哈尔滨工程大学科技处知识产权办公室

优先权:

专利代理机构:

代理人:

PDF下载: PDF下载
内容摘要

本发明提供的是一种基于分数傅立叶变换的图像压缩方法。第一步,计算图像的(0.5,0)、(0,0.5)、(1,0)、(0,1)阶分数傅立叶变换;第二步,根据Wigner分布与分数傅立叶变换的关系计算所有变换阶次的分数傅立叶二阶矩;第三步,根据各二阶矩的值找到最优的变换阶数,并将图像变换到最优的分数傅立叶域,从而使图像的能量集中在尽可能少的分数傅立叶系数上;第四步,忽略所有幅值小于某一阈值的系数,并采用熵编码方法对其余的分数傅立叶系数进行编码,实现图像压缩。本发明计算速度快,采用快速离散分数傅立叶变换算法,本方法的时间复杂度仅为O(N×NlogN);对平稳信号和非平稳的二维信号都具有较好的压缩效果。

权利要求书

1: 一种基于分数傅立叶变换的图像压缩方法, 其特征是 : 第一步, 计算图像的 (0.5, 0)、 (0, 0.5)、 (1, 0)、 (0, 1) 阶分数傅立叶变换 ; 第二步, 根据 Wigner 分布与分数傅立叶变换 的关系计算所有变换阶次的分数傅立叶二阶矩 ; 第三步, 根据各二阶矩的值找到最优的变 换阶数, 并将图像变换到最优的分数傅立叶域, 从而使图像的能量集中在尽可能少的分数 傅立叶系数上 ; 第四步, 忽略所有幅值小于某一阈值的系数, 并采用熵编码方法对其余的分 数傅立叶系数进行编码, 实现图像压缩。
2: 根据权利要求 1 基于分数傅立叶变换的图像压缩方法, 其特征是其具体步骤为 : 步骤 1, 利用快速离散分数傅立叶变换算法计算 I0.5, n)、 I0, n)、 I1, n)、 I0, 0(m, 0.5(m, 0(m, n) ; 1(m, 步 骤 2, 根据式 r 计 算 wr0, wr0.5, 根据式 0, 0 和 w 1, 0, c 计算 wc0, wc0, 0, 0.5 和 w 0, 1, 其中, 为图像的总能量 ; 步 骤 3, 分 别 令 a = 0, Δa, 2Δa, …, 1, b = 0, Δb, 2Δb, …, 1, 根 据 式 c 分别计算 wra, b 和 w a, b, r r 其 中,系 数 c1 = wr0, c2 = wr0.5, c3 = wr1, c4 = wc0, c5 = wc0, 0, 0-0.5(w 0, 0+w 1, 0), 0, 0, c c c6 = wc0, 0+w 0, 1), 1 ; 0.5-0.5(w 0, 步骤 4, 令目标函数为 r c J = w a, b×w a, b 计算使 J 取最小值的变换阶数 (a′, b′ ) ; 步骤 5, 利用快速离散分数傅立叶变换算法计算 Ia′, n) ; b′ (m, 步骤 6, 若 |Ia′, n)| < 0.2max(|Ia′, n)|) 则令 Ia′, n) = 0, 其中 max b′ (m, b′ (m, b′ (m, 为取最大值函数 ; 步骤 7, 将步骤 6 计算得到的 Ia′, n) 用最大熵方法编码, 实现图像压缩 ; b′ (m, 上述步骤中 I(m, n) 为 N×N 输入图像, a, b 为行、 列两个方向的分数傅立叶变换阶数, 并且 0 ≤ a, b ≤ 1, Δa、 Δb 分别为 a、 b 的增量 ; Ia, n) 为 I(m, n) 的 (a, b) 阶分数傅立 b(m, r c 叶变换, w a, n) 的行向和列向二阶矩。 b 为 Ia, b 和 w a, b(m,

说明书


基于分数傅立叶变换的图像压缩方法

    ( 一 ) 技术领域
     本发明涉及的是一种数据压缩方法, 特别是涉及一种数字图像压缩方法。 ( 二 ) 背景技术
     随着计算机与数字通信技术的迅速发展, 图像压缩受到了人们越来越多的关注。 从本质上讲, 图像压缩就是用尽可能少的代码 ( 符号 ) 来表示尽可能多的数据信息。现有 的数字图像压缩技术可分为两大类 : 有损压缩和无损压缩。其中有损压缩主要包括行程长 度编码、 熵编码以及 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 的增量。Ia,b(m, n) 为 I(m, n) 的 (a, b) 阶分数傅立叶 r c 变换, w a, n) 的行向和列向二阶矩。 b 和 w a, b 为 Ia, b(m,
     本发明所提出的图像压缩算法具体实现步骤如下 :
     步骤 1, 利用快速离散分数傅立叶变换算法计算 I0.5, n)、 I0, n)、 I1, n)、 0(m, 0.5(m, 0(m, I0, n)。 1(m, r c
     步骤 2, 根据式 (1) 计算 wr0, wr0.5, 根据式 (2) 计算 wc0, wc0, 0, 0 和 w 1, 0, 0, 0.5 和 w 0, 1
     其中,为图像的总能量。步骤 3, 分别令 a = 0, Δa, 2Δa,…, 1, b = 0, Δb, 2Δb,…, 1, 根据式 (3)、 (4) r c 分别计算 w a, b 和 w a, b。
     r r 其中, 系数 c1 = wr0, c2 = wr0.5, c3 = wr1, c4 = wc0, c5 = wc0, 0, 0-0.5(w 0, 0+w 1, 0), 0, 0, c c c6 = wc0, 0+w 0, 1), 1 0.5-0.5(w 0,
     步骤 4, 令目标函数为
     计算使 J 取最小值的变换阶数 (a′, b′ )。
     步骤 5, 利用快速离散分数傅立叶变换算法计算 Ia′, n)。 b′ (m,
     步骤 6, 若 |Ia′, n)| < 0.2max(|Ia′, n)|) 则令 Ia′, n) = 0, 其中 b′ (m, b′ (m, b′ (m, max 为取最大值函数。
     步骤 7, 将步骤 6 计算得到的 Ia′, n) 用最大熵方法编码, 实现图像压缩。 b′ (m,
     4

基于分数傅立叶变换的图像压缩方法.pdf_第1页
第1页 / 共4页
基于分数傅立叶变换的图像压缩方法.pdf_第2页
第2页 / 共4页
基于分数傅立叶变换的图像压缩方法.pdf_第3页
第3页 / 共4页
点击查看更多>>
资源描述

《基于分数傅立叶变换的图像压缩方法.pdf》由会员分享,可在线阅读,更多相关《基于分数傅立叶变换的图像压缩方法.pdf(4页珍藏版)》请在专利查询网上搜索。

1、10申请公布号CN101944235A43申请公布日20110112CN101944235ACN101944235A21申请号200910072925622申请日20090918G06T9/00200601G06F17/1420060171申请人哈尔滨工程大学地址150001黑龙江省哈尔滨市南岗区南通大街145号1号楼哈尔滨工程大学科技处知识产权办公室72发明人郑丽颖石大明田凯54发明名称基于分数傅立叶变换的图像压缩方法57摘要本发明提供的是一种基于分数傅立叶变换的图像压缩方法。第一步,计算图像的05,0、0,05、1,0、0,1阶分数傅立叶变换;第二步,根据WIGNER分布与分数傅立叶变换的。

2、关系计算所有变换阶次的分数傅立叶二阶矩;第三步,根据各二阶矩的值找到最优的变换阶数,并将图像变换到最优的分数傅立叶域,从而使图像的能量集中在尽可能少的分数傅立叶系数上;第四步,忽略所有幅值小于某一阈值的系数,并采用熵编码方法对其余的分数傅立叶系数进行编码,实现图像压缩。本发明计算速度快,采用快速离散分数傅立叶变换算法,本方法的时间复杂度仅为ONNLOGN;对平稳信号和非平稳的二维信号都具有较好的压缩效果。51INTCL19中华人民共和国国家知识产权局12发明专利申请权利要求书1页说明书2页CN101944242A1/1页21一种基于分数傅立叶变换的图像压缩方法,其特征是第一步,计算图像的05,。

3、0、0,05、1,0、0,1阶分数傅立叶变换;第二步,根据WIGNER分布与分数傅立叶变换的关系计算所有变换阶次的分数傅立叶二阶矩;第三步,根据各二阶矩的值找到最优的变换阶数,并将图像变换到最优的分数傅立叶域,从而使图像的能量集中在尽可能少的分数傅立叶系数上;第四步,忽略所有幅值小于某一阈值的系数,并采用熵编码方法对其余的分数傅立叶系数进行编码,实现图像压缩。2根据权利要求1基于分数傅立叶变换的图像压缩方法,其特征是其具体步骤为步骤1,利用快速离散分数傅立叶变换算法计算I05,0M,N、I0,05M,N、I1,0M,N、I0,1M,N;步骤2,根据式计算WR0,0,WR05,0和WR1,0,根。

4、据式计算WC0,0,WC0,05和WC0,1,其中,为图像的总能量;步骤3,分别令A0,A,2A,1,B0,B,2B,1,根据式分别计算WRA,B和WCA,B,其中,系数C1WR0,0,C2WR05,005WR0,0WR1,0,C3WR1,0,C4WC0,0,C5WC0,0505WC0,0WC0,1,C6WC0,1;步骤4,令目标函数为JWRA,BWCA,B计算使J取最小值的变换阶数A,B;步骤5,利用快速离散分数傅立叶变换算法计算IA,BM,N;步骤6,若|IA,BM,N|02MAX|IA,BM,N|则令IA,BM,N0,其中MAX为取最大值函数;步骤7,将步骤6计算得到的IA,BM,N用最。

5、大熵方法编码,实现图像压缩;上述步骤中IM,N为NN输入图像,A,B为行、列两个方向的分数傅立叶变换阶数,并且0A,B1,A、B分别为A、B的增量;IA,BM,N为IM,N的A,B阶分数傅立叶变换,WRA,B和WCA,B为IA,BM,N的行向和列向二阶矩。权利要求书CN101944235ACN101944242A1/2页3基于分数傅立叶变换的图像压缩方法一技术领域0001本发明涉及的是一种数据压缩方法,特别是涉及一种数字图像压缩方法。二背景技术0002随着计算机与数字通信技术的迅速发展,图像压缩受到了人们越来越多的关注。从本质上讲,图像压缩就是用尽可能少的代码符号来表示尽可能多的数据信息。现有。

6、的数字图像压缩技术可分为两大类有损压缩和无损压缩。其中有损压缩主要包括行程长度编码、熵编码以及LZW方法;而有损压缩方法则包括色度抽样法、变换编码方法以及分形编码方法等。变换编码方法是最常用的图像压缩方法,其主要思想是使用如离散余弦变换、傅立叶变换、或者小波变换来表达图像,然后再采用一种无损编码方法实现压缩。0003分数傅立叶变换是傅立叶变换的推广形式,它将信号表示为CHIRP信号的组合,克服了傅立叶变换将时域信息完全摒弃的缺点,通过一种简单的方式实现了信号的从纯粹时间域到纯频率域的全过程的综合描述,能够展示出信号从纯时间域到纯频率域的所有变化特征。研究表明,通过分数傅立叶变换,可以将信号表示。

7、为紧致形式,即信号的大部分能量都集中在少数几个分数傅立叶系数上,例如CHIRP率为2的CHIRP信号,其/2ATAN051/阶分数傅立叶变换为脉冲信号。0004与本发明相关的公开报道有1、CVIJAYA,BHAT,JSSIGNALCOMPRESSIONUSINGDISCRETEFRACTIONALFOURIERTRANSFORMANDSETPARTITIONINGINHIERARCHICALTREESIGNALPROCESSING2006,8619761983;2、ISYETIK,KUTAY,MA,OZAKTAS,HMIMAGEREPRESENTATIONANDCOMPRESSIONWITHT。

8、HEFRACTIONALFOURIERTRANSFORMOPTICSCOMMUNICATIONS,2001,197275278等。三发明内容0005本发明的目的在于提供一种可以将信号表示为紧致形式的特性,实现快速数字图像压缩的基于分数傅立叶变换的图像压缩方法。0006本发明的目的是这样实现的0007第一步,计算图像的05,0、0,05、1,0、0,1阶分数傅立叶变换;第二步,根据WIGNER分布与分数傅立叶变换的关系计算所有变换阶次的分数傅立叶二阶矩;第三步,根据各二阶矩的值找到最优的变换阶数,并将图像变换到最优的分数傅立叶域,从而使图像的能量集中在尽可能少的分数傅立叶系数上;第四步,忽略所有。

9、幅值小于某一阈值的系数,并采用熵编码方法对其余的分数傅立叶系数进行编码,实现图像压缩。0008本发明的有益效果是,1计算速度快,采用快速离散分数傅立叶变换算法,本方法的时间复杂度仅为ONNLOGN;2对平稳信号和非平稳的二维信号都具有较好的压缩效果。四具体实施方式说明书CN101944235ACN101944242A2/2页40009下面举例对本发明做更详细地描述0010设IM,N为NN输入图像,A,B为行、列两个方向的分数傅立叶变换阶数,并且0A,B1,A、B分别为A、B的增量。IA,BM,N为IM,N的A,B阶分数傅立叶变换,WRA,B和WCA,B为IA,BM,N的行向和列向二阶矩。001。

10、1本发明所提出的图像压缩算法具体实现步骤如下0012步骤1,利用快速离散分数傅立叶变换算法计算I05,0M,N、I0,05M,N、I1,0M,N、I0,1M,N。0013步骤2,根据式1计算WR0,0,WR05,0和WR1,0,根据式2计算WC0,0,WC0,05和WC0,1001400150016其中,为图像的总能量。0017步骤3,分别令A0,A,2A,1,B0,B,2B,1,根据式3、4分别计算WRA,B和WCA,B。001800190020其中,系数C1WR0,0,C2WR05,005WR0,0WR1,0,C3WR1,0,C4WC0,0,C5WC0,0505WC0,0WC0,1,C6WC0,10021步骤4,令目标函数为00220023计算使J取最小值的变换阶数A,B。0024步骤5,利用快速离散分数傅立叶变换算法计算IA,BM,N。0025步骤6,若|IA,BM,N|02MAX|IA,BM,N|则令IA,BM,N0,其中MAX为取最大值函数。0026步骤7,将步骤6计算得到的IA,BM,N用最大熵方法编码,实现图像压缩。说明书CN101944235A。

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 物理 > 计算;推算;计数


copyright@ 2017-2020 zhuanlichaxun.net网站版权所有
经营许可证编号:粤ICP备2021068784号-1