CN109255770A - 一种新型图像变换域降采样方法 - Google Patents

一种新型图像变换域降采样方法 Download PDF

Info

Publication number
CN109255770A
CN109255770A CN201810775713.3A CN201810775713A CN109255770A CN 109255770 A CN109255770 A CN 109255770A CN 201810775713 A CN201810775713 A CN 201810775713A CN 109255770 A CN109255770 A CN 109255770A
Authority
CN
China
Prior art keywords
image
transform domain
block
matrix
norm
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
CN201810775713.3A
Other languages
English (en)
Other versions
CN109255770B (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 Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201810775713.3A priority Critical patent/CN109255770B/zh
Publication of CN109255770A publication Critical patent/CN109255770A/zh
Application granted granted Critical
Publication of CN109255770B publication Critical patent/CN109255770B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T9/00—Image coding
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00—Image enhancement or restoration
    • G06T5/10—Image enhancement or restoration using non-spatial domain filtering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20021—Dividing image into blocks, subimages or windows
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20048—Transform domain processing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Compression Of Band Width Or Redundancy In Fax (AREA)
  • Compression Or Coding Systems Of Tv Signals (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

本发明属于图像处理领域，提供一种新型图像变换域降采样方法，用以克服传统的图像降采样方法中因高频信息损失而引起的图像质量退化问题；本发明首先将图像划分为大小为N×N图像宏块，然后在像素域建立基于图像宏块的L2‑范数和L1‑范数联合优化目标函数，通过最小化目标函数，在变换域得到具有较强稀疏特性的(N/2)×(N/2)大小的变换系数块。本发明通过最小化像素域的L2‑范数和L1‑范数联合优化，在变换域实现对图像信号的降采样；并且在保持图像原始特征的同时，增强图像在变换域的稀疏特性，有效降低变换域系数块的数据量，为进一步的数据压缩提供有力的保障。

Description

一种新型图像变换域降采样方法
技术领域
本发明属于图像处理领域，涉及数字图像的降采样技术，具体提供一种新型图像变换域降采样方法。
背景技术
图像是人类获取信息的主要途径，对与视觉紧密相关的图像进行高效的处理，具有十分重要的意义；图像处理的应用范围从宇航领域扩展到生物医学、信息科学、资源环境科学、天文学、物理学、工业、农业、国防、教育、艺术各个领域与行业，对人们的日常生活产生了重大的影响；然而，由于存储和传输原始的图像需要耗费巨大的空间和带宽，因此对图像数据进行压缩处理十分有必要。经研究发现，图像信号经过变换以后，信息主要集中在少部分的低频变换系数上，对这些低频系数进行高效的处理，能够有效节约传输带宽和存储空间，达到去除图像空间冗余的目的。
为了能够在变换域去除图像内部的空间冗余，传统的变换域降采样方法采用了直接去除部分高频系数的策略，如文献“Down-scaling for better transformcompression”，这种降采样方法实现起来比较简单，但由于丢弃了代表图像细节的高频信息，容易造成图像质量的严重退化。基于此，为了在变换域降采样后能够有效保留图像的高频信息，本发明提出将图像变换域的稀疏化和降采样相结合，通过联合优化L2-范数和L1-范数，实现高效的图像变换域降采样。
发明内容
本发明的目的在于提供一种新型图像变换域降采样方法，用以克服传统的图像降采样方法中因高频信息损失而引起的图像质量退化问题；本发明通过最小化像素域的L2-范数和L1-范数联合优化目标函数，在变换域实现对图像信号的降采样；并且在保持图像原始特征的同时，增强图像在变换域的稀疏特性，有效降低变换域系数块的数据量。
为实现上述目的，本发明采用的技术方案为：
一种新型图像变换域降采样方法，包括以下步骤：
步骤1.图像的预处理
将大小为W×H的图像，按照JPEG图像压缩标准中的图像分块方法划分为大小为N×N的图像宏块，记为bi，其中，W为图像的宽度、H为图像的高度，i为图像块的索引；
步骤2.图像宏块的列向量化
对每一个N×N的图像宏块bi，按照图像块列向量化方法，转换为N2×1的列向量，记为xi；
步骤3.构建变换矩阵
首先，构建N2×N2的变换矩阵D：其中，C为N×N的离散余弦变换矩阵、为Kronecker乘法中的乘法算子；
然后，将变换矩阵D的第1～8个列向量，第17～24个列向量，第33～40个列向量，……，第[(N-1)2+1]/2～[N(N-1)]/2个列向量依次取出，并按照从左到右的顺序组成一个矩阵，记为Λ；
步骤4.产生图像宏块的降采样系数
首先，构建目标函数为：其中，Xi为图像宏块bi的降采样系数矢量，||·||2表示L2-范数，||·||1表示L1-范数；
然后，采用快速迭代阈值收敛算法求解上述目标函数，得到的最优解
最后，将中的系数按照从上到下，从左到右的顺序构建一个(N/2)×(N/2)的系数块，记为
本发明的有益效果在于：
本发明提供一种新型图像变换域降采样方法，能够直接产生图像在变换域具有强稀疏特性的变换系数，这些变换系数在用于图像压缩编码时，能够有效降低编码的码率；同时，所产生的变换系数能够充分保留原始图像的高频信息，在利用这些系数进行图像重建时能够保证较高的重建质量；因此，将本发明所提供的图像变换域降采样方法应用于实际的图像编码中，能够有效提升图像的编码效率。
附图说明
图1为本发明一种新型图像变换域降采样方法的流程示意图。
图2为实施例中测试图像中所有16×16的系数块中左上角64个低频系数的统计分布图。
图3为实施例中基于本发明得到的离散余弦变换系数统计分布图。
具体实施方式
下面结合附图和实施例对本发明进一步详细说明。
本发明提供一种新型图像变换域降采样方法，为了方便描述本发明的内容，首先做以下术语定义：
定义1：传统的JPEG图像压缩标准中的图像分块方法
传统的JPEG图像压缩标准中的图像分块方式以图像的左上角为起点，按照从左到右、从上到下的顺序，将原始图像划分为多个互不重叠的等尺寸的图像块，具体描述过程参见文献“JPEG(Joint Photographic Experts Group):ISO/IEC IS 10918–1/ITU-TRecommendation T.81,Digital Compression and Coding of Continuous-ToneStillImage,1993”；
定义2：传统的图像块列向量化方法
传统的图像块列向量化方法是把一个大小为n×n的图像块的每个列向量按照从左到右的顺序依次去除，排列成一个(n2×1)的列向量的方法，这里n代表图像块的高度和宽度。
定义3：传统的产生离散余弦变换矩阵的方法
传统的产生离散余弦变换矩阵的方法是根据离散余弦变换的定义，计算变换矩阵中的每一个变换系数，产生任意大小的离散余弦变换矩阵的方法，具体描述过程参见文献“数字视频编码技术原理”，高文、赵德斌、马思伟著，科学出版社；
定义4：传统的Kronecker乘法
传统的Kronecker乘法可以表示为其中，表示Kronecker乘法算子，A是大小为m×n的矩阵，aij为矩阵A中的元素，i代表横坐标索引，j代表纵坐标索引。
B是大小为p×q的矩阵，bij为矩阵B中的元素，i代表横坐标索引，j代表纵坐标索引。
C是大小为mp×nq的矩阵，
具体描述过程参见文献“矩阵分析与应用(第2版)”，张贤达著，清华大学出版社；
定义5：传统的快速迭代阈值收敛算法
传统的快速迭代阈值收敛算法(A fast iterative shrinkage-thresholdingalgorithm)是一种基于梯度下降思想的优化方法，与传统梯度下降方法相比，迭代步骤中修正方式的选择更为合理，收敛速度更快，具体步骤以及默认参数参见文献“A fastiterative shrinkage-thresholding algorithm for linear inverse problems”；
本发明主要采用仿真实验的方式验证该算法的可行性，所有步骤都经过实验验证。为了实现图像变换域降采样，具体实施步骤如下：
步骤1，图像的预处理
将大小为W×H的图像，按照传统的JPEG图像压缩标准中的图像分块方法划分为大小为16×16的图像宏块，记为bi，这里，W代表图像的宽度，H代表图像的高度，i代表图像块的索引下标；
步骤2，图像宏块的列向量化
对每一个N×N的图像宏块bi，按照传统的图像块列向量化方法，转换为大小为256×1的列向量，记为xi。
步骤3，产生变换矩阵
首先，按照传统的产生离散余弦变换矩阵的方法，产生大小为16×16的离散余弦变换矩阵，记为C；
接着，用传统的Kronecker乘法和离散余弦变换矩阵C,产生一个大小为256×256的变换矩阵，记为D，这里，符号代表传统Kronecker乘法中的乘法算子；
最后，用矩阵D的第1～8个列向量、第17～24个列向量、第33～40个列向量、……、第113～120个列向量按照从左到右的顺序组成一个矩阵，记为Λ；
步骤4，产生图像宏块的降采样系数
首先，定义图像宏块bi的降采样系数矢量为Xi；
其次，定义λ＝2，用Xi、步骤2中得到的列向量xi和步骤3中的得到的矩阵Λ构建目标函数为这里λ为目标方程中的调节参数；
接着，用传统的快速迭代阈值收敛算法求解目标函数的最优解，将得到的最优解，记为
最后，将中的系数按照从上到下，从左到右的顺序产生一个8×8的系数块，记为
将实施例应用于Lena、Goldhill和Peppers三幅分辨率均为512×512的经典图像中，附图2是测试图像中所有16×16的系数块中左上角64个低频系数的统计分布图，附图3是基于本发明得到的离散余弦变换系数统计分布图。对比附图2和附图3可以明显的发现，使用本发明所提供的方法产生的低频系数更加稀疏，因此为进一步的高效压缩，提供了保障。
以上所述，仅为本发明的具体实施方式，本说明书中所公开的任一特征，除非特别叙述，均可被其他等效或具有类似目的的替代特征加以替换；所公开的所有特征、或所有方法或过程中的步骤，除了互相排斥的特征和/或步骤以外，均可以任何方式组合。

Claims (1)

1.一种新型图像变换域降采样方法，包括以下步骤：
步骤1.图像的预处理
将大小为W×H的图像，按照JPEG图像压缩标准中的图像分块方法划分为大小为N×N的图像宏块，记为bi，其中，W为图像的宽度、H为图像的高度，i为图像块的索引；
步骤2.图像宏块的列向量化
对每一个N×N的图像宏块bi，按照图像块列向量化方法，转换为N2×1的列向量，记为xi；
步骤3.构建变换矩阵
首先，构建N2×N2的变换矩阵D：其中，C为N×N的离散余弦变换矩阵、为Kronecker乘法中的乘法算子；
然后，将变换矩阵D的第1～8个列向量，第17～24个列向量，第33～40个列向量，……，第[(N-1)2+1]/2～[N(N-1)]/2个列向量依次取出，并按照从左到右的顺序组成一个矩阵，记为Λ；
步骤4.产生图像宏块的降采样系数
首先，构建目标函数为：其中，Xi为图像宏块bi的降采样系数矢量；
然后，采用快速迭代阈值收敛算法求解上述目标函数，得到的最优解
最后，将中的系数按照从上到下，从左到右的顺序构建一个(N/2)×(N/2)的系数块，记为
CN201810775713.3A 2018-07-16 2018-07-16 一种图像变换域降采样方法 Active CN109255770B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810775713.3A CN109255770B (zh) 2018-07-16 2018-07-16 一种图像变换域降采样方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810775713.3A CN109255770B (zh) 2018-07-16 2018-07-16 一种图像变换域降采样方法

Publications (2)

Publication Number Publication Date
CN109255770A true CN109255770A (zh) 2019-01-22
CN109255770B CN109255770B (zh) 2022-11-08

Family

ID=65049149

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810775713.3A Active CN109255770B (zh) 2018-07-16 2018-07-16 一种图像变换域降采样方法

Country Status (1)

Country Link
CN (1) CN109255770B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110049242A (zh) * 2019-04-18 2019-07-23 腾讯科技（深圳）有限公司 一种图像处理方法和装置

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5740284A (en) * 1993-11-30 1998-04-14 Polaroid Corporation Coding method and apparatus for resampling and filtering images using discrete cosine transforms
US6141456A (en) * 1997-12-31 2000-10-31 Hitachi America, Ltd. Methods and apparatus for combining downsampling and inverse discrete cosine transform operations
CN102542542A (zh) * 2012-01-02 2012-07-04 西安电子科技大学 基于非局部稀疏模型的图像去噪方法
CN102769745A (zh) * 2012-06-21 2012-11-07 哈尔滨工业大学 插值依赖的图像自适应下采样方法
CN103700074A (zh) * 2013-12-23 2014-04-02 电子科技大学 基于离散余弦变换系数分布的自适应压缩感知采样方法
CN103871041A (zh) * 2014-03-21 2014-06-18 上海交通大学 基于认知正则化参数构建的图像超分辨率重构方法
CN104616027A (zh) * 2015-02-06 2015-05-13 华东交通大学 一种非相邻的图结构稀疏人脸识别方法
CN105306936A (zh) * 2015-07-17 2016-02-03 福州大学 一种基于块压缩感知的图像编码方法
WO2016050870A1 (en) * 2014-09-30 2016-04-07 The Provost, Fellows, Foundation Scholars, & The Other Members Of Board, Of The College Of The Holy & Undiv. Trinity Of Queen Elizabeth Near Dublin A generalized frequency division multiplexing transceiver
CN105678697A (zh) * 2015-12-30 2016-06-15 北京工业大学 一种基于dct域本征变换的人脸图像超分辨率重建方法
CN107610049A (zh) * 2017-08-21 2018-01-19 华侨大学 基于稀疏正则化技术和加权引导滤波的图像超分辨率方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5740284A (en) * 1993-11-30 1998-04-14 Polaroid Corporation Coding method and apparatus for resampling and filtering images using discrete cosine transforms
US6141456A (en) * 1997-12-31 2000-10-31 Hitachi America, Ltd. Methods and apparatus for combining downsampling and inverse discrete cosine transform operations
CN102542542A (zh) * 2012-01-02 2012-07-04 西安电子科技大学 基于非局部稀疏模型的图像去噪方法
CN102769745A (zh) * 2012-06-21 2012-11-07 哈尔滨工业大学 插值依赖的图像自适应下采样方法
CN103700074A (zh) * 2013-12-23 2014-04-02 电子科技大学 基于离散余弦变换系数分布的自适应压缩感知采样方法
CN103871041A (zh) * 2014-03-21 2014-06-18 上海交通大学 基于认知正则化参数构建的图像超分辨率重构方法
WO2016050870A1 (en) * 2014-09-30 2016-04-07 The Provost, Fellows, Foundation Scholars, & The Other Members Of Board, Of The College Of The Holy & Undiv. Trinity Of Queen Elizabeth Near Dublin A generalized frequency division multiplexing transceiver
CN104616027A (zh) * 2015-02-06 2015-05-13 华东交通大学 一种非相邻的图结构稀疏人脸识别方法
CN105306936A (zh) * 2015-07-17 2016-02-03 福州大学 一种基于块压缩感知的图像编码方法
CN105678697A (zh) * 2015-12-30 2016-06-15 北京工业大学 一种基于dct域本征变换的人脸图像超分辨率重建方法
CN107610049A (zh) * 2017-08-21 2018-01-19 华侨大学 基于稀疏正则化技术和加权引导滤波的图像超分辨率方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
SHUYUAN ZHU等: "Compression-Dependent Transform-Domain Downward Conversion for Block-Based Image Coding", 《IEEE TRANSACTIONS ON IMAGE PROCESSING》 *
THONG T. DO等: "Fast and Efficient Compressive Sensing Using Structurally Random Matrices", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 *
何兴宇等: "基于多观测向量序列降采样恢复的稀疏矩阵重构", 《系统工程与电子技术》 *
左功玉: "基于Kronecker积的压缩感知雷达成像方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110049242A (zh) * 2019-04-18 2019-07-23 腾讯科技（深圳）有限公司 一种图像处理方法和装置

Also Published As

Publication number Publication date
CN109255770B (zh) 2022-11-08

Similar Documents

Publication Publication Date Title
Liu et al. Compressive sampling-based image coding for resource-deficient visual communication
CN106101725B (zh) 一种基于压缩感知理论和空域下采样技术的图像压缩方法
US20140072240A1 (en) Method and apparatus for converting an image, and method and apparatus for inverse converting an image
JPH09107464A (ja) 圧縮されたデータ形式で与えられる画像を変倍する方法
CN106204670B (zh) 一种基于优化下采样算法和压缩感知理论的图像压缩方法
EP2145476B1 (en) Image compression and decompression using the pixon method
CN105976409B (zh) 一种基于压缩感知理论的图像压缩方法
CN111881920B (zh) 一种大分辨率图像的网络适配方法及神经网络训练装置
AU2002259268C1 (en) Apparatus and method for encoding and computing a discrete cosine transform using a butterfly processor
Kekre et al. ImageCompression Using Real Fourier Transform, Its Wavelet Transform And Hybrid Wavelet With DCT
Hassan et al. Color image compression based on DCT, differential pulse coding modulation, and adaptive shift coding
Siddeq et al. A novel 2D image compression algorithm based on two levels DWT and DCT transforms with enhanced minimize-matrix-size algorithm for high resolution structured light 3D surface reconstruction
CN103700074A (zh) 基于离散余弦变换系数分布的自适应压缩感知采样方法
Mafijur Rahman et al. A low complexity lossless Bayer CFA image compression
Jeong et al. An overhead-free region-based JPEG framework for task-driven image compression
CN107146259A (zh) 一种基于压缩感知理论的新型图像压缩方法
CN107146260B (zh) 一种基于均方误差的图像压缩感知采样方法
CN109255770B (zh) 一种图像变换域降采样方法
CN112767503B (zh) 一种基于jpeg压缩编码的数字病理全切片图像快速分析方法
CN104683818B (zh) 基于双正交不变集多小波的图像压缩方法
Divakara et al. High speed area optimized hybrid da architecture for 2d-dtcwt
Zheng et al. Sequential convolution and runge-kutta residual architecture for image compressed sensing
Kang et al. Generating images in compressed domain using generative adversarial networks
CN108184127A (zh) 一种可配置的多尺寸dct变换硬件复用架构
CN108989812B (zh) 一种基于图像压缩的去块效应方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant