[go: up one dir, main page]

CN108364258A - A kind of method and system improving image resolution ratio - Google Patents

A kind of method and system improving image resolution ratio Download PDF

Info

Publication number
CN108364258A
CN108364258A CN201810120907.XA CN201810120907A CN108364258A CN 108364258 A CN108364258 A CN 108364258A CN 201810120907 A CN201810120907 A CN 201810120907A CN 108364258 A CN108364258 A CN 108364258A
Authority
CN
China
Prior art keywords
image
subgraph
amplified
cartoon
texture
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
CN201810120907.XA
Other languages
Chinese (zh)
Other versions
CN108364258B (en
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.)
Zhejiang Normal University CJNU
Original Assignee
Zhejiang Normal University CJNU
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 Zhejiang Normal University CJNU filed Critical Zhejiang Normal University CJNU
Priority to CN201810120907.XA priority Critical patent/CN108364258B/en
Publication of CN108364258A publication Critical patent/CN108364258A/en
Application granted granted Critical
Publication of CN108364258B publication Critical patent/CN108364258B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4053Scaling of whole images or parts thereof, e.g. expanding or contracting based on super-resolution, i.e. the output image resolution being higher than the sensor resolution
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/23Clustering techniques
    • G06F18/232Non-hierarchical techniques
    • G06F18/2321Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
    • G06F18/23213Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4007Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/40Image enhancement or restoration using histogram techniques

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Probability & Statistics with Applications (AREA)
  • Image Processing (AREA)

Abstract

The invention discloses a kind of method and system improving image resolution ratio.Method includes:Original image is decomposed using opposite total variance method, obtains cartoon subgraph and texture subgraph;The cartoon subgraph is amplified using interpolation method, obtains amplified cartoon subgraph;The texture subgraph is amplified using homotopy Method, obtains amplified texture subgraph;The amplified cartoon subgraph and amplified texture subgraph are synthesized, high-definition picture is obtained.The present invention is amplified, then synthesized using different methods by the cartoon subgraph and texture subgraph different to amplification demand, obtains clearer high-definition picture, while improving high-definition picture quality, reduce processing time respectively.

Description

一种提高图像分辨率的方法及系统A method and system for improving image resolution

技术领域technical field

本发明涉及图像放大领域,特别涉及一种提高图像分辨率的方法及系统。The invention relates to the field of image enlargement, in particular to a method and system for improving image resolution.

背景技术Background technique

单幅图像的超分辨率重建技术是指利用一幅低分辨率图像来放大图像,获得一幅高分辨率图像的技术,以提高图像的分辨率和清晰度。目前该技术已经广泛应用于图像和视频检测和显示的各种仪器中,如医学影像、视频监控、遥感图像、高清电视等领域。The super-resolution reconstruction technology of a single image refers to the technology of using a low-resolution image to enlarge the image to obtain a high-resolution image, so as to improve the resolution and clarity of the image. At present, this technology has been widely used in various instruments for image and video detection and display, such as medical imaging, video surveillance, remote sensing images, high-definition television and other fields.

为此,专家和学者们已经提出很多单幅图像超分辨率的方法。这些单幅图像的超分辨率重建方法主要分为基于像素点插值、基于图像投影和基于学习的方法。其中,在基于插值的方法中主要包括双线性插值法和双三次插值法。这类方法计算量小,可满足图像实时放大的要求,但是通过这类方法重建的高分辨率图像的质量较低,细节部分不丰富,图像比较模糊。基于图像投影的方法有迭代反投影法、和凸集投影法POCS(projection onconvex set)等,这类方法一般需要图像多次的迭代投影后才能获得效果比较好的高分辨率的图像,这类图像超分辨率方法容易受到图像噪声的干扰,所获得图像不稳定,同时也有图像细节部分丢失问题。基于学习的超分辨率方法主要是对大量的高低分辨率图像样本进行训练最终获得图像的先验知识,再对图像进行超分辨率重建,经典的方法有最大后验概率估计法MAP(maximum a posterior),领域嵌入方法NE(NeighborEmbedding),基于样例学习的方法;最近,Yang等利用线性规划求解低分辨率图像块的稀疏表示,并将此表示系数与高分辨率字典进行线性组合得到高分辨率图像块;这类方法的缺点在于训练需要很长的时间和过程,同时训练要求的存储空间很大,其优点在于所得到的高分辨率图像的质量一般来说要好于以上两类方法。因此如何在提高高分辨率图像的质量的同时减少处理时间,成为一个亟待解决的技术问题。For this reason, experts and scholars have proposed many single image super-resolution methods. These single-image super-resolution reconstruction methods are mainly divided into pixel-point interpolation-based, image-projection-based, and learning-based methods. Among them, interpolation-based methods mainly include bilinear interpolation and bicubic interpolation. This kind of method has a small amount of calculation and can meet the requirements of real-time image enlargement, but the quality of the high-resolution image reconstructed by this kind of method is low, the details are not rich, and the image is relatively blurred. Image projection-based methods include iterative back-projection method and convex set projection method POCS (projection onconvex set), etc. This type of method generally requires multiple iterative projections of the image to obtain a better high-resolution image. The image super-resolution method is easily disturbed by image noise, the obtained image is unstable, and there is also the problem of partial loss of image details. The super-resolution method based on learning is mainly to train a large number of high- and low-resolution image samples to finally obtain the prior knowledge of the image, and then perform super-resolution reconstruction on the image. The classic method has the maximum a posteriori probability estimation method MAP (maximum a posterior), domain embedding method NE (NeighborEmbedding), a method based on example learning; recently, Yang et al. used linear programming to solve the sparse representation of low-resolution image blocks, and linearly combined the representation coefficients with high-resolution dictionaries to obtain high-resolution Resolution image blocks; the disadvantage of this type of method is that training takes a long time and process, and the storage space required for training is large, and its advantage is that the quality of the obtained high-resolution image is generally better than the above two types of methods . Therefore, how to reduce the processing time while improving the quality of high-resolution images has become an urgent technical problem to be solved.

发明内容Contents of the invention

本发明的目的是提供一种提高图像分辨率的方法及系统,以在提高高分辨率图像的质量的同时减少处理时间。The object of the present invention is to provide a method and system for increasing image resolution, so as to reduce processing time while improving the quality of high-resolution images.

为实现上述目的,本发明提供了如下方案:To achieve the above object, the present invention provides the following scheme:

一种提高图像分辨率的方法,所述方法包括如下步骤:A method for improving image resolution, said method comprising the steps of:

采用相对总变差法对原图像进行分解,获得卡通子图像和纹理子图像;The original image is decomposed by relative total variation method to obtain cartoon sub-image and texture sub-image;

采用插值法对所述卡通子图像进行放大,获得放大后的卡通子图像;Using an interpolation method to enlarge the cartoon sub-image to obtain the enlarged cartoon sub-image;

采用同伦方法对所述纹理子图像进行放大,获得放大后的纹理子图像;Using a homotopy method to amplify the texture sub-image to obtain the amplified texture sub-image;

对所述放大后的卡通子图像和放大后的纹理子图像进行合成,获得高分辨率图像。Combining the enlarged cartoon sub-image and the enlarged texture sub-image to obtain a high-resolution image.

可选的,所述采用相对总变差法对图像进行分解,获得卡通子图像和纹理子图像;具体包括:Optionally, the relative total variation method is used to decompose the image to obtain cartoon sub-images and texture sub-images; specifically including:

对原图像进行分割,获得分割后的图像;Segment the original image to obtain the segmented image;

利用最小化窗函数的图像全变分的方法,对分割后的图像进行分解,获得卡通子图像和纹理子图像。The segmented image is decomposed to obtain cartoon sub-image and texture sub-image by using the method of image total variation with minimum window function.

可选的,所述对原图像进行分割,获得分割后的图像,具体包括:Optionally, the segmenting the original image to obtain the segmented image specifically includes:

分别利用模板和原图像进行卷积,得到x方向的梯度图像Ix和y方向的梯度图像Iytemplate and Carry out convolution with original image, obtain the gradient image I x of x direction and the gradient image I y of y direction;

计算所述x方向的梯度图像和所述y方向的梯度图像的梯度图像模值IgCalculating the gradient image modulus I g of the gradient image in the x direction and the gradient image in the y direction;

其中,(i,j)为像素点坐标;Among them, (i, j) is the pixel coordinates;

将梯度图像模值大于第一阈值的像素点作为边缘像素点,将梯度图像模值大于第二阈值的像素点作为候选边缘像素点;Taking the pixels whose gradient image modulus is greater than the first threshold as edge pixels, and the pixels whose gradient image modulus is greater than the second threshold as candidate edge pixels;

选取多个边缘像素点,分别以每个边缘像素点为起始点,连接与所述边缘像素点在八个方向上直接相邻的边缘像素点或候选边缘像素点,获得若干封闭区域,并将位于大的封闭区域内的小的封闭区域删除,获得若干没有包含关系的封闭区域;Select a plurality of edge pixels, take each edge pixel as a starting point, connect the edge pixels or candidate edge pixels directly adjacent to the edge pixel in eight directions, obtain several closed areas, and The small closed area located in the large closed area is deleted, and several closed areas without inclusion relationship are obtained;

对每个没有包含关系的封闭区域做灰度直方图,并对所述灰度直方图采用K-means方法来进行分割,获得分割后的图像。Make a grayscale histogram for each closed area that has no inclusion relationship, and use the K-means method to segment the grayscale histogram to obtain a segmented image.

可选的,所述窗函数为:Optionally, the window function is:

点(xp,yp)是当前中心像素点,(xq,yq)是(xp,yp)的全变分的像素点,σ2 p,q是(xq,yq)和(xp,yp)所在物体内的像素值的方差,cp,q为相乘因子,Tg为第三阈值。Point (x p , y p ) is the current center pixel point, (x q , y q ) is the total variation pixel point of (x p , y p ), σ 2 p,q is the variance of the pixel values in the object where (x q ,y q ) and (x p ,y p ) are located, c p,q is the multiplication factor, and T g is the third threshold.

可选的,所述采用插值法对卡通子图像进行放大,获得放大后的卡通子图像,具体包括:Optionally, the interpolation method is used to enlarge the cartoon sub-image to obtain the enlarged cartoon sub-image, which specifically includes:

对所述卡通子图像的像素点采用插值模板函数进行插值,获得放大后的卡通子图像;所述插值模板函数为:The pixels of the cartoon sub-image are interpolated using an interpolation template function to obtain the enlarged cartoon sub-image; the interpolation template function is:

L(x,y)=N1*o(x,y)*d(x,y)*sinc(α*x/σ2)L(x,y)=N 1 *o(x,y)*d(x,y)*sinc(α*x/σ 2 )

*sinc(α*y/σ2)*exp(β*(x2+y2)/2σ2) (3)*sinc(α*y/σ 2 ) * exp(β * (x 2 +y 2 )/2σ 2 ) (3)

其中,o(x,y)表示是否利用当前像素点的值来估计未知像素点的值,若当前像素点和未知像素点在同一个物体中,则o(x,y)=1,即利用当前像素点;否则,o(x,y)=0,即不利用当前像素点,d(x,y)表示当前像素点和插值模板中心点的距离;函数sinc(x)=(sin(x)/x);σ2是参与对未知像素点插值的所有像素点的均方差;α和β是常数。Among them, o(x, y) indicates whether to use the value of the current pixel to estimate the value of the unknown pixel. If the current pixel and the unknown pixel are in the same object, then o(x, y)=1, that is, use The current pixel; otherwise, o(x,y)=0, that is, the current pixel is not used, and d(x,y) represents the distance between the current pixel and the center point of the interpolation template; the function sinc(x)=(sin(x )/x); σ 2 is the mean square error of all pixels involved in the interpolation of unknown pixels; α and β are constants.

可选的,所述采用同伦方法对所述纹理子图像进行放大,获得放大后的纹理子图像,具体包括:Optionally, the amplifying the texture sub-image by using a homotopy method to obtain the amplified texture sub-image specifically includes:

采用K-SVD字典学习法(K-SVD字典学习法为K均值分类法和奇异值分解svd法相结合的方法)离线地获得所述纹理子图像的图像块字典,所述图像块字典包括高分辨率图像块字典和低分辨率图像块字典;Using the K-SVD dictionary learning method (the K-SVD dictionary learning method is a method combining the K-means classification method and the singular value decomposition svd method) to obtain the image block dictionary of the texture sub-image offline, the image block dictionary includes high-resolution High-resolution image block dictionary and low-resolution image block dictionary;

对所述纹理子图像,利用所述图像块字典和正交匹配跟踪方法进行在线放大,获得初始高分辨率图像;For the texture sub-image, use the image block dictionary and the orthogonal matching tracking method to perform online amplification to obtain an initial high-resolution image;

对所述初始高分辨率图像进行最近邻的加边处理;获得加边处理后的高分辨图像;performing nearest-neighbor edge processing on the initial high-resolution image; obtaining a high-resolution image after edge processing;

对所述加边处理后的高分辨率图像进行第一次同伦处理,获得第一高分辨率图像;performing the first homotopy processing on the high-resolution image after the edge processing to obtain the first high-resolution image;

对所述第一高分辨率图像进行第二次同伦处理,获得放大后的高分辨率纹理子图像。The second homotopy processing is performed on the first high-resolution image to obtain an enlarged high-resolution texture sub-image.

可选的,所述对所述放大后的卡通子图像和放大后的纹理子图像进行合成,获得高分辨率图像,具体包括:Optionally, said synthesizing the enlarged cartoon sub-image and the enlarged texture sub-image to obtain a high-resolution image specifically includes:

采用公式(4)对所述放大后的卡通子图像和放大后的纹理子图像进行合成,获得高分辨率图像;Adopt formula (4) to synthesize the cartoon sub-image after the enlargement and the texture sub-image after the enlargement, obtain a high-resolution image;

fH=fc+ft1*G(ft) (4)f H =f c +f t1 *G(f t ) (4)

其中,fH为高分辨率的图像,ft放大后的纹理图像,fc为放大后的卡通图像,G(ft)为对图像ft求Robert梯度的模值,λ1为常数。Among them, f H is a high-resolution image, f t is an enlarged texture image, f c is an enlarged cartoon image, G( ft ) is the modulus value of Robert's gradient for image f t , and λ 1 is a constant.

一种提高图像分辨率的系统,所述系统包括:A system for increasing image resolution, the system comprising:

图像分解模块,用于采用相对总变差法对原图像进行分解,获得卡通子图像和纹理子图像;The image decomposition module is used to decompose the original image using the relative total variation method to obtain cartoon sub-images and texture sub-images;

卡通子图像放大模块,用于采用插值法对所述卡通子图像进行放大,获得放大后的卡通子图像;The cartoon sub-image enlargement module is used to amplify the cartoon sub-image by interpolation to obtain the enlarged cartoon sub-image;

纹理子图像放大模块,用于采用同伦方法对所述纹理子图像进行放大,获得放大后的纹理子图像;A texture sub-image enlarging module, configured to amplify the texture sub-image using a homotopy method to obtain an enlarged texture sub-image;

图像合成模块,用于对所述放大后的卡通子图像和放大后的纹理子图像进行合成,获得高分辨率图像。The image synthesis module is used to synthesize the enlarged cartoon sub-image and the enlarged texture sub-image to obtain a high-resolution image.

可选的,所述图像分解模块具体包括:Optionally, the image decomposition module specifically includes:

图像分割子模块,用于对原图像进行分割,获得分割后的图像;The image segmentation sub-module is used to segment the original image to obtain the segmented image;

图像分解子模块,用于利用最小化窗函数的图像全变分的方法,对分割后的图像进行分解,获得卡通子图像和纹理子图像。The image decomposition sub-module is used to decompose the segmented image to obtain cartoon sub-images and texture sub-images by using the image total variation method of minimizing the window function.

可选的,所述图像分割子模块具体包括:Optionally, the image segmentation submodule specifically includes:

卷积单元,用于分别利用模板对原图像进行卷积,得到x方向的梯度图像Ix和y方向的梯度图像IyConvolutional units for exploiting templates separately and Carry out convolution to original image, obtain the gradient image I x of x direction and the gradient image I y of y direction;

梯度图像模值计算单元,用于计算所述x方向的梯度图像和所述y方向的梯度图像的梯度图像模值Ig A gradient image modulus calculation unit, configured to calculate the gradient image modulus I g of the gradient image in the x direction and the gradient image in the y direction;

其中,(i,j)为像素点坐标;Among them, (i, j) is the pixel coordinates;

像素点选取单元,用于将梯度图像模值大于第一阈值的像素点作为边缘像素点,将梯度图像模值大于第二阈值的像素点作为候选边缘像素点;A pixel point selection unit, configured to use pixels whose gradient image modulus is greater than the first threshold as edge pixels, and use pixels whose gradient image modulus is greater than the second threshold as candidate edge pixels;

封闭区域获取单元,用于选取多个边缘像素点,分别以每个边缘像素点为起始点,连接与所述边缘像素点在八个方向上直接相邻的边缘像素点或候选边缘像素点,获得若干封闭区域,并将位于大的封闭区域内的小的封闭区域删除,获得若干没有包含关系的封闭区域;The closed area acquisition unit is used to select a plurality of edge pixel points, and respectively use each edge pixel point as a starting point to connect edge pixel points or candidate edge pixel points directly adjacent to the edge pixel point in eight directions, Obtain several closed areas, delete the small closed areas located in the large closed area, and obtain several closed areas without inclusion relationship;

图像分割单元,用于对每个没有包含关系的封闭区域做灰度直方图,并对所述灰度直方图采用K-means方法来进行分割,获得分割后的图像。The image segmentation unit is configured to make a grayscale histogram for each closed area that has no inclusion relationship, and use the K-means method to segment the grayscale histogram to obtain a segmented image.

根据本发明提供的具体实施例,本发明公开了以下技术效果:According to the specific embodiments provided by the invention, the invention discloses the following technical effects:

本发明公开了一种提高图像分辨率的方法及系统。方法包括:采用相对总变差法对原图像进行分解,获得卡通子图像和纹理子图像;采用插值法对所述卡通子图像进行放大,获得放大后的卡通子图像;采用同伦方法对所述纹理子图像进行放大,获得放大后的纹理子图像;对所述放大后的卡通子图像和放大后的纹理子图像进行合成,获得高分辨率图像。本发明通过对放大需求不同的卡通子图像和纹理子图像,采用不同的方法分别放大,再进行合成,得到更清晰的高分辨率图像,在提高高分辨率图像质量的同时,减少了处理时间。The invention discloses a method and system for improving image resolution. The method includes: using the relative total variation method to decompose the original image to obtain a cartoon sub-image and a texture sub-image; using an interpolation method to amplify the cartoon sub-image to obtain an enlarged cartoon sub-image; using a homotopy method to obtain a cartoon sub-image The texture sub-image is enlarged to obtain an enlarged texture sub-image; the enlarged cartoon sub-image and the enlarged texture sub-image are synthesized to obtain a high-resolution image. The present invention uses different methods to amplify cartoon sub-images and texture sub-images with different amplifying requirements, and then synthesizes them to obtain clearer high-resolution images, which reduces the processing time while improving the quality of high-resolution images .

附图说明Description of drawings

为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the accompanying drawings required in the embodiments. Obviously, the accompanying drawings in the following description are only some of the present invention. Embodiments, for those of ordinary skill in the art, other drawings can also be obtained according to these drawings without paying creative labor.

图1为本发明提供的一种提高图像分辨率的方法的流程图;Fig. 1 is a flow chart of a method for improving image resolution provided by the present invention;

图2为本发明提供的一种提高图像分辨率的方法在进行第一次同伦处理的初始图像块;Fig. 2 is that a kind of method for improving image resolution provided by the present invention is carrying out the initial image block of homotopy processing for the first time;

图3为本发明提供的一种提高图像分辨率的方法在进行第二次同伦处理的初始图像块;Fig. 3 is that a kind of method for improving image resolution provided by the present invention is carrying out the initial image block of the second homotopy processing;

图4为本发明提供的一种提高图像分辨率的系统的结构框图;Fig. 4 is a structural block diagram of a system for improving image resolution provided by the present invention;

图5为本发明提供的人物图像的各种不同处理方法的效果对比图;Fig. 5 is a comparison diagram of the effects of various different processing methods for character images provided by the present invention;

图6为本发明提供的蝴蝶图像的各种不同处理方法的效果对比图;Fig. 6 is a comparison diagram of effects of various different processing methods of butterfly images provided by the present invention;

其中,图5中的(a)为原始图像;(b)为双三次插值图像;(c)为Yang提出的方法的图像;(d)为Zedye提出的方法的图像;(e)为本发明的方法的图像;图6中的(a)为原始图像;(b)为双三次插值图像;(c)为Yang提出的方法的图像;(d)为Zedye提出的方法的图像;(e)为本发明的方法的图像。Wherein, (a) in Fig. 5 is the original image; (b) is the bicubic interpolation image; (c) is the image of the method proposed by Yang; (d) is the image of the method proposed by Zedye; (e) is the present invention The image of the method; (a) in Figure 6 is the original image; (b) is the bicubic interpolation image; (c) is the image of the method proposed by Yang; (d) is the image of the method proposed by Zedye; (e) is an image of the method of the present invention.

具体实施方式Detailed ways

本发明的目的是提供一种提高图像分辨率的方法及系统,以在提高高分辨率图像的质量的同时减少处理时间。The object of the present invention is to provide a method and system for increasing image resolution, so as to reduce processing time while improving the quality of high-resolution images.

为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对发明作进一步详细的说明。In order to make the above objects, features and advantages of the present invention more comprehensible, the invention will be further described in detail below in conjunction with the accompanying drawings and specific embodiments.

基于插值的方法对平滑区域的放大具有较好的效果,基于学习的方法对图像的细节部分的放大具有较好的效果。在本发明所提出的方法中将两者有机地结合起来,以获得更好的性能。首先,利用改进的图像分解方法把图像分解为卡通(cartoon)和纹理(texture)部分。其中,卡通部分提取的是图像的结构信息,像素点值只在物体交界处有较大变化,在物体内部的像素点值变化小,图像平滑。纹理部分提取的是图像的细节部分,其中的像素点值的变化较大。然后,对卡通部分,利用插值的方法进行放大;对纹理部分,利用同伦方法进行放大。最后,把以上两个部分的处理结果进行较好的合成,得到更清晰的高分辨率图像。以下,将对所提出方法的每个细节加以更详细地说明。The method based on interpolation has a better effect on the enlargement of the smooth area, and the method based on learning has a better effect on the enlargement of the detailed part of the image. In the method proposed by the present invention, the two are organically combined to obtain better performance. First, an improved image decomposition method is used to decompose the image into cartoon and texture parts. Among them, the cartoon part extracts the structural information of the image, and the pixel value only has a large change at the junction of the object, and the change of the pixel value inside the object is small, and the image is smooth. The texture part extracts the detail part of the image, where the pixel value changes greatly. Then, for the cartoon part, use the method of interpolation to enlarge; for the texture part, use the method of homotopy to enlarge. Finally, the processing results of the above two parts are better synthesized to obtain a clearer high-resolution image. In the following, each detail of the proposed method will be explained in more detail.

如图1所示,本发明提供了一种提高图像分辨率的方法,所述方法包括如下步骤:As shown in Figure 1, the present invention provides a kind of method that improves image resolution, and described method comprises the following steps:

步骤101,采用相对总变差法对原图像进行分解,获得卡通子图像和纹理子图像。在原有的卡通和纹理分解方法RTV中,是利用最小化加窗函数的图像全变分的方法来获得卡通部分的。这个方法有如下缺点:1.图像的物体边缘处的全变分不应统计到需要最小化的目标函数中,这种像素值的变化是不同物体一般具有不同像素值的正常变化;2.窗函数值的大小和范围应当随着物体的像素值而自适应地变化,以提高卡通部分提取的效果。例如,对于像素值在物体内较快变化时,窗函数应当取较大范围和较大值。为此,本发明提出的采用相对总变差法对原图像进行分解,获得卡通子图像和纹理子图,具体包括:Step 101, using relative total variation method to decompose the original image to obtain cartoon sub-image and texture sub-image. In the original cartoon and texture decomposition method RTV, the cartoon part is obtained by minimizing the total variation of the image with the window function. This method has the following disadvantages: 1. The total variation at the edge of the object of the image should not be counted into the objective function that needs to be minimized. The change of this pixel value is a normal change that different objects generally have different pixel values; 2. The window The size and range of the function value should change adaptively with the pixel value of the object, so as to improve the effect of cartoon part extraction. For example, when the pixel value changes quickly in the object, the window function should take a larger range and a larger value. For this reason, the relative total variation method proposed by the present invention decomposes the original image to obtain cartoon sub-images and texture sub-images, specifically including:

首先,对原图像进行分割,获得分割后的图像,具体包括:First, segment the original image to obtain the segmented image, including:

分别利用模板和原图像进行卷积,得到x方向的梯度图像Ix和y方向的梯度图像Iytemplate and Carry out convolution with original image, obtain the gradient image I x of x direction and the gradient image I y of y direction;

计算所述x方向的梯度图像和所述y方向的梯度图像的梯度图像模值IgCalculating the gradient image modulus I g of the gradient image in the x direction and the gradient image in the y direction;

其中,(i,j)为像素点坐标;Among them, (i, j) is the pixel coordinates;

将梯度图像模值大于第一阈值的像素点作为边缘像素点,将梯度图像模值大于第二阈值的像素点作为候选边缘像素点;Taking the pixels whose gradient image modulus is greater than the first threshold as edge pixels, and the pixels whose gradient image modulus is greater than the second threshold as candidate edge pixels;

选取多个边缘像素点,分别以每个边缘像素点为起始点,连接与所述边缘像素点在八个方向上直接相邻的边缘像素点或候选边缘像素点,获得若干封闭区域,并将位于大的封闭区域内的小的封闭区域删除,获得若干没有包含关系的封闭区域;Select a plurality of edge pixels, take each edge pixel as a starting point, connect the edge pixels or candidate edge pixels directly adjacent to the edge pixel in eight directions, obtain several closed areas, and The small closed area located in the large closed area is deleted, and several closed areas without inclusion relationship are obtained;

对每个没有包含关系的封闭区域做灰度直方图,并对所述灰度直方图采用K-means方法来进行分割,获得分割后的图像。Make a grayscale histogram for each closed area that has no inclusion relationship, and use the K-means method to segment the grayscale histogram to obtain a segmented image.

其次,利用最小化窗函数的图像全变分的方法,对分割后的图像进行分解,获得卡通子图像和纹理子图像。所述窗函数为:Second, the segmented image is decomposed to obtain the cartoon sub-image and the texture sub-image by using the method of image total variation that minimizes the window function. The window function is:

点(xp,yp)是当前中心像素点,(xq,yq)是(xp,yp)的全变分的像素点,σ2 p,q是(xq,yq)和(xp,yp)所在物体内的像素值的方差,cp,q为相乘因子,其和所在物体内的检测出的边缘像素点数成正比,Tg为第三阈值。Point (x p , y p ) is the current center pixel point, (x q , y q ) is the total variation pixel point of (x p , y p ), σ 2 p,q is the variance of the pixel values in the object where (x q ,y q ) and (x p ,y p ) are located, c p,q is the multiplication factor, which is equal to the detected edge pixels in the object Points are proportional, T g is the third threshold.

步骤102,采用插值法对所述卡通子图像进行放大,获得放大后的卡通子图像。本发明提出的使用插值的方法来完成对卡通子图像的高分辨率图像中未知像素的估计和计算和以往的方法不同,本发明的插值法在插值过程中将遵循如下两个原则:1、只利用和未知像素点在同一物体中的像素点来进行插值,以尽量避免插值后像素点模糊的情况。2、自适应地调整插值模板的大小,对于图像中平滑的部分将取范围比较大的模板,以让更多一些的像素参与插值,对于图像中变化较大的部分,采用较小范围的模板,只利用和当前像素相关性大的像素来进行插值。这样,也可以避免插值后像素点模糊的情况。Step 102, using an interpolation method to enlarge the cartoon sub-image to obtain an enlarged cartoon sub-image. The method of using interpolation proposed by the present invention to complete the estimation and calculation of unknown pixels in the high-resolution image of the cartoon sub-image is different from previous methods. The interpolation method of the present invention will follow the following two principles in the interpolation process: 1. Only use pixels in the same object as unknown pixels for interpolation to avoid blurring of pixels after interpolation. 2. Adaptively adjust the size of the interpolation template. For the smooth part of the image, a template with a relatively large range will be used to allow more pixels to participate in the interpolation. For the part with large changes in the image, a template with a smaller range will be used , and only use pixels with a large correlation with the current pixel for interpolation. In this way, the blurring of pixels after interpolation can also be avoided.

具体包括,对所述卡通子图像的像素点采用插值模板函数进行插值,获得放大后的卡通子图像;所述插值模板函数为:Specifically comprising, interpolating the pixels of the cartoon sub-image using an interpolation template function to obtain an enlarged cartoon sub-image; the interpolation template function is:

L(x,y)=N1*o(x,y)*d(x,y)*sinc(α*x/σ2)L(x,y)=N 1 *o(x,y)*d(x,y)*sinc(α*x/σ 2 )

*sinc(α*y/σ2)*exp(β*(x2+y2)/2σ2) (3) * sinc(α * y/σ 2 ) * exp(β * (x 2 +y 2 )/2σ 2 ) (3)

其中,o(x,y)表示是否利用当前像素点的值来估计未知像素点的值,若当前像素点和未知像素点在同一个物体中,则o(x,y)=1,即利用当前像素点;否则,o(x,y)=0,即不利用当前像素点,d(x,y)表示当前像素点和插值模板中心点的距离;函数sinc(x)=(sin(x)/x);σ2是参与对未知像素点插值的所有像素点的均方差;α和β是常数。由于在实际应用中,需要对插值模板函数进行截断来使用,这里只取L(x,y)>T1(T1为需要在实验中确定的函数)的像素参与插值。N1为一归一化常数,使所有的参与插值的系数的累加值为1。Among them, o(x, y) indicates whether to use the value of the current pixel to estimate the value of the unknown pixel. If the current pixel and the unknown pixel are in the same object, then o(x, y)=1, that is, use The current pixel; otherwise, o(x,y)=0, that is, the current pixel is not used, and d(x,y) represents the distance between the current pixel and the center point of the interpolation template; the function sinc(x)=(sin(x )/x); σ 2 is the mean square error of all pixels involved in the interpolation of unknown pixels; α and β are constants. Since the interpolation template function needs to be truncated for use in practical applications, only pixels with L(x,y)>T 1 (T 1 is a function that needs to be determined in experiments) are selected to participate in the interpolation. N 1 is a normalization constant, so that the cumulative value of all the coefficients participating in the interpolation is 1.

步骤103,采用同伦方法对所述纹理子图像进行放大,获得放大后的纹理子图像。本发明所述提出的对纹理子图像处理的方法中,离线的高分辨率字典和低分辨率的字典的学习过程采用K-SVD字典学习方法先得出关于纹理部分高分辨率字典Dh和低分辨率字典Dl。对于所述纹理子图像的在线放大过程,首先利用得到的图像块字典和正交匹配跟踪OMP(Orthogonal Matching Pursuit)方法进行放大,得到初始的高分辨率的图像。接着,本发明对此初始的高分辨率图像进行最近邻的加边处理,然后运行两次同伦方法进行在线恢复。在对图像加边处理中,未知像素值由离其最近的像素点值确定,所加边的高度为M1,宽度为N1,M1×N1为进行每次高分辨率纹理图像块处理的图像块的大小。在第一次同伦处理中,采用图2所示的初始高分辨率图像块进行更新。在第二次同伦处理中,采用图3所示的初始高分辨率图像块进行更新。其中,在图2中显示了一个高分辨率图像中待放大的图像块的示意图,其中带阴影部分的是已经放大的区域,白色部分是需要进行像素值估计的部分。在图3中,显示了一个在第二次同伦处理中的图像块,其中带阴影的部分认为是已知像素值的部分,空白部分认为是需要再次更新的部分。Step 103, using a homotopy method to enlarge the texture sub-image to obtain an enlarged texture sub-image. In the method for texture sub-image processing proposed by the present invention, the learning process of off-line high-resolution dictionary and low-resolution dictionary adopts the K-SVD dictionary learning method to first obtain the high-resolution dictionary D h and Low resolution dictionary D l . For the online enlargement process of the texture sub-image, firstly, the obtained image block dictionary and the Orthogonal Matching Pursuit (OMP) method are used for enlargement to obtain an initial high-resolution image. Next, the present invention performs nearest neighbor edge processing on the initial high-resolution image, and then runs the homotopy method twice for online recovery. In the process of adding edges to the image, the unknown pixel value is determined by the value of the nearest pixel point, the height of the added edge is M 1 , the width is N 1 , and M 1 ×N 1 is for each high-resolution texture image block The size of the processed image blocks. In the first homotopy process, the initial high-resolution image block shown in Fig. 2 is used for updating. In the second homotopy processing, the initial high-resolution image block shown in Fig. 3 is used for updating. Wherein, FIG. 2 shows a schematic diagram of an image block to be enlarged in a high-resolution image, where the shaded part is the enlarged region, and the white part is the part where pixel value estimation is required. In Fig. 3, an image block in the second homotopy processing is shown, where the shaded part is considered to be a part with known pixel values, and the blank part is considered to be a part that needs to be updated again.

具体包括:采用K-SVD字典学习法离线地获得所述纹理子图像的图像块字典,所述图像块字典包括高分辨率图像块字典和低分辨率图像块字典;对所述纹理子图像,利用所述图像块字典和正交匹配跟踪方法进行在线放大,获得初始高分辨率图像;对所述初始高分辨率图像进行最近邻的加边处理;获得加边处理后的高分辨图像;对所述加边处理后的高分辨率图像进行第一次同伦处理,获得第一高分辨率图像;对所述第一高分辨率图像进行第二次同伦处理,获得放大后的高分辨率纹理子图像。Specifically include: using the K-SVD dictionary learning method to obtain the image block dictionary of the texture sub-image offline, the image block dictionary includes a high-resolution image block dictionary and a low-resolution image block dictionary; for the texture sub-image, Using the image block dictionary and the orthogonal matching tracking method to perform online amplification to obtain an initial high-resolution image; perform nearest-neighbor edge processing on the initial high-resolution image; obtain a high-resolution image after edge processing; The high-resolution image after the edge processing is subjected to the first homotopy processing to obtain the first high-resolution image; the first high-resolution image is subjected to the second homotopy processing to obtain the enlarged high-resolution image Rate texture subimage.

步骤104,对所述放大后的卡通子图像和放大后的纹理子图像进行合成,获得高分辨率图像,具体包括:采用公式(4)对所述放大后的卡通子图像和放大后的纹理子图像进行合成,获得高分辨率图像;Step 104, combining the enlarged cartoon sub-image and the enlarged texture sub-image to obtain a high-resolution image, specifically including: using formula (4) to synthesize the enlarged cartoon sub-image and the enlarged texture The sub-images are synthesized to obtain a high-resolution image;

fH=fc+ft1*G(ft)(4)f H =f c +f t1 *G(f t )(4)

其中,fH为高分辨率的图像,ft放大后的纹理图像,fc为放大后的卡通图像,G(ft)为对图像ft求Robert梯度的模值,λ1为常数。Among them, f H is a high-resolution image, f t is an enlarged texture image, f c is an enlarged cartoon image, G( ft ) is the modulus value of Robert's gradient for image f t , and λ 1 is a constant.

图4为本发明提供的一种提高图像分辨率的系统的结构框图。参见图4,一种提高图像分辨率的系统,所述系统包括:FIG. 4 is a structural block diagram of a system for improving image resolution provided by the present invention. Referring to Fig. 4, a system for improving image resolution, the system includes:

图像分解模块401,用于采用相对总变差法对原图像进行分解,获得卡通子图像和纹理子图像;所述图像分解模块具体包括:图像分割子模块,用于对原图像进行分割,获得分割后的图像;图像分解子模块,用于利用最小化窗函数的图像全变分的方法,对分割后的图像进行分解,获得卡通子图像和纹理子图像。The image decomposition module 401 is used to decompose the original image by using the relative total variation method to obtain the cartoon sub-image and the texture sub-image; the image decomposition module specifically includes: an image segmentation sub-module, which is used to segment the original image to obtain The segmented image; the image decomposition sub-module, which is used to decompose the segmented image to obtain cartoon sub-images and texture sub-images by using the image total variation method of minimizing the window function.

所述图像分割子模块具体包括:卷积单元,用于分别利用模板对原图像进行卷积,得到x方向的梯度图像Ix和y方向的梯度图像Iy;梯度图像模值计算单元,用于计算所述x方向的梯度图像和所述y方向的梯度图像的梯度图像模值Ig The image segmentation sub-module specifically includes: a convolution unit for using templates respectively and The original image is convolved to obtain the gradient image I x in the x direction and the gradient image I y in the y direction; the gradient image modulus calculation unit is used to calculate the gradient image in the x direction and the gradient image in the y direction Gradient image modulus I g ;

其中,(i,j)为像素点坐标;像素点选取单元,用于将梯度图像模值大于第一阈值的像素点作为边缘像素点,将梯度图像模值大于第二阈值的像素点作为候选边缘像素点;封闭区域获取单元,用于选取多个边缘像素点,分别以每个边缘像素点为起始点,连接与所述边缘像素点在八个方向上直接相邻的边缘像素点或候选边缘像素点,获得若干封闭区域,并将位于大的封闭区域内的小的封闭区域删除,获得若干没有包含关系的封闭区域;图像分割单元,用于对每个没有包含关系的封闭区域做灰度直方图,并对所述灰度直方图采用K-means方法来进行分割,获得分割后的图像。Wherein, (i, j) is pixel point coordinates; Pixel point selection unit is used for taking the pixel point whose gradient image modulus value is greater than the first threshold as an edge pixel point, and using the pixel point whose gradient image modulus value is greater than the second threshold value as a candidate Edge pixel points; a closed area acquisition unit, used to select a plurality of edge pixel points, and each edge pixel point is used as a starting point to connect edge pixel points or candidate points directly adjacent to the edge pixel point in eight directions Edge pixels, to obtain several closed areas, and delete the small closed areas located in the large closed area, to obtain several closed areas without containment relationship; the image segmentation unit is used to gray each closed area without containment relationship Intensity histogram, and the K-means method is used to segment the gray histogram to obtain a segmented image.

卡通子图像放大模块402,用于采用插值法对所述卡通子图像进行放大,获得放大后的卡通子图像;The cartoon sub-image enlargement module 402 is used to enlarge the cartoon sub-image by interpolation to obtain the enlarged cartoon sub-image;

纹理子图像放大模块403,用于采用同伦方法对所述纹理子图像进行放大,获得放大后的纹理子图像;A texture sub-image enlarging module 403, configured to amplify the texture sub-image using a homotopy method to obtain an enlarged texture sub-image;

图像合成模块404,用于对所述放大后的卡通子图像和放大后的纹理子图像进行合成,获得高分辨率图像。The image synthesis module 404 is configured to synthesize the enlarged cartoon sub-image and the enlarged texture sub-image to obtain a high-resolution image.

本发明还通过实验的方式验证本发明所提出的方法和系统的性能,具体包括:The present invention also verifies the performance of the method and system proposed by the present invention through experiments, specifically including:

选用windows 8操作系统,Matlab2013b软件、8GB内存,2.9GHz处理器作为实验平台,选用实验样本库的91幅图像作为训练样本,对待测测试图像从平均峰值信噪比PSNR(peak signal to noise ratio)和RMSE(Root Mean Square Error)。如公式(5)和(6)所示:Windows 8 operating system, Matlab2013b software, 8GB memory, and 2.9GHz processor were selected as the experimental platform, and 91 images from the experimental sample library were selected as training samples. The average peak signal-to-noise ratio PSNR (peak signal to noise ratio) and RMSE (Root Mean Square Error). As shown in formulas (5) and (6):

其中,M×N是图像的大小,x(i,j)和分别是原始图像和重构图像在坐标(i,j)处的图像值。where M×N is the size of the image, x(i,j) and are the image values at coordinates (i, j) of the original image and the reconstructed image, respectively.

在本发明的实验中,由于人眼对图像的亮度分量很敏感,本发明先把彩色图像分为YCbCr三个通道,然后只取Y(亮度)通道分量进行如下的超分辨率处理。低分辨率图像被重建放大三倍,对于图像的纹理部分,选取训练的图像块为50000,K-SVD的迭代次数设置为40,字典的原子数设置为1000。重叠矩阵大小为3x3,重叠像素为1。实验效果如图5-6所示。In the experiment of the present invention, since the human eye is very sensitive to the brightness component of the image, the present invention first divides the color image into three channels of YCbCr, and then only takes the Y (brightness) channel component to carry out the following super-resolution processing. The low-resolution image is reconstructed and enlarged three times. For the texture part of the image, 50000 image blocks are selected for training, the number of iterations of K-SVD is set to 40, and the number of atoms in the dictionary is set to 1000. The overlap matrix size is 3x3 and the overlap pixel is 1. The experimental results are shown in Figure 5-6.

从图5、6可以看出,利用双三次插值法处理的结果过于平滑,同时从人物和蝴蝶的细节来看,伴有锯齿效应。Yang的方法恢复的已经较好,但是脸上的雀斑等没有恢复,丢失的细节较多。本发明提出的方法注重于细节的恢复,所以在细节方面更胜一筹。It can be seen from Figures 5 and 6 that the results processed by the bicubic interpolation method are too smooth, and at the same time, from the details of the characters and butterflies, there is a jagged effect. Yang's method has recovered better, but the freckles on the face have not been recovered, and many details have been lost. The method proposed by the present invention pays attention to the recovery of details, so it is superior in details.

本发明采用客观评价方法对三种不同的重建方法的PSNR,RMSE进行计算。选用标准测试图像进行比较,从表1和表2可以看出,本发明的方法与双三次插值方法、Yang的方法、Zeyde方法相比,具有更好的性能,由此可见,本发明的方法能够获得较好的重建图像的质量。本发明使用的方法比Zeyde的方法提高了约0.2db。The present invention uses an objective evaluation method to calculate PSNR and RMSE of three different reconstruction methods. Select the standard test image for comparison, as can be seen from Table 1 and Table 2, the method of the present invention has better performance compared with the bicubic interpolation method, the method of Yang, the Zeyde method, thus it can be seen that the method of the present invention A better reconstruction image quality can be obtained. The method used in the present invention has an improvement of about 0.2db over Zeyde's method.

表1不同方法处理后图像的PSNR值Table 1 PSNR values of images processed by different methods

图像image 双三次插值bicubic interpolation Yang方法Yang method Zeyde方法Zeyde method 本发明方法The method of the invention zebrazebra 26.634026.6340 27.953127.9531 28.496828.4968 28.814728.8147 ppt3ppt3 23.706623.7066 24.977724.9777 25.178925.1789 25.533325.5333 monarchmonarch 29.425929.4259 30.713330.7133 31.092131.0921 31.374831.3748 lennalenna 31.677831.6778 32.638432.6384 32.995132.9951 33.225533.2255 flowersflowers 27.230127.2301 28.248828.2488 28.423628.4236 28.630328.6303 comiccomic 23.115623.1156 23.902323.9023 23.972623.9726 24.143324.1433

表2不同方法处理后图像的RMSE值Table 2 RMSE values of images processed by different methods

图像image 双三次插值bicubic interpolation Yang方法Yang method Zeyde方法Zeyde method 本发明方法The method of the invention zebrazebra 11.880611.8806 10.206710.2067 9.58749.5874 9.24289.2428 ppt3ppt3 16.642216.6422 14.376514.3765 14.047314.0473 13.485713.4857 monarchmonarch 8.61488.6148 7.42817.4281 7.11117.1111 6.88336.8833 lennalenna 11.092611.0926 9.86519.8651 9.66859.6685 9.44129.4412 flowersflowers 11.092611.0926 9.86519.8651 9.66859.6685 9.44129.4412 comiccomic 17.813917.8139 16.271416.2714 16.140216.1402 15.826115.8261

根据本发明提供的具体实施例,本发明公开了以下技术效果:According to the specific embodiments provided by the invention, the invention discloses the following technical effects:

1、本发明在图像放大之前,把图像首先分离出纹理和卡通部分,其中纹理部分主要包含图像中的细节部分,卡通部分包含图像中的结构部分。因此,在图像放大过程中,需要对纹理部分进行更精细的放大,以克服传统方法中在图像放大后丢失图像细节的问题。同时,本发明不同于传统的TV(total variation)分解方法,拥有比TV方法更好的分解方式,并在实验结果中证明了它的效果,所以本发明先对这个方法进行了改进,然后将图像分解为基于物体的结构部分和纹理部分。1. Before the image is enlarged, the present invention first separates the image into texture and cartoon parts, wherein the texture part mainly includes the detail parts in the image, and the cartoon part includes the structure parts in the image. Therefore, in the process of image enlargement, finer enlargement of the texture part is required to overcome the problem of losing image details after image enlargement in traditional methods. At the same time, the present invention is different from the traditional TV (total variation) decomposition method, and has a better decomposition method than the TV method, and its effect has been proved in the experimental results, so the present invention first improves this method, and then The image is decomposed into an object-based structure part and a texture part.

2、对于图像中的卡通部分,其高频分量较少,使用传统的bi-cubic方法对此进行放大,就已经能够取得较好的性能,同时可以节约图像放大所需的时间。2. For the cartoon part in the image, its high-frequency components are less, and using the traditional bi-cubic method to amplify this part can already achieve better performance and save the time required for image enlargement.

3、对于图像中的纹理部分,需要在放大时保持其细节。因此,需要采用效果好的基于学习的方法。在图像训练阶段,为了克服Zeyde方法不能保证高低分辨率字典具有相同的稀疏表示的缺点,加入了联合训练的思想以确保高低图像块在过完备字典上具有相同的表示系数。对纹理部分的在线放大阶段,可以采用传统的OMP(正交匹配跟踪)方法,其优点是运行速度快且放大效果较好。而且本发明使用L1-homotopy(基于1范数的同伦)方法,和OMP方法比较,此方法的重构效果更好。3. For the texture part in the image, it is necessary to maintain its details when zooming in. Therefore, effective learning-based methods are required. In the image training phase, in order to overcome the shortcomings of the Zeyde method that the high and low resolution dictionaries have the same sparse representation, the idea of joint training is added to ensure that the high and low image blocks have the same representation coefficients on the overcomplete dictionary. The traditional OMP (Orthogonal Matching Pursuit) method can be used for the online enlargement stage of the texture part, and its advantage is that the operation speed is fast and the enlargement effect is good. Moreover, the present invention uses the L1-homotopy (homotopy based on 1 norm) method, and compared with the OMP method, the reconstruction effect of this method is better.

本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。Each embodiment in this specification is described in a progressive manner, each embodiment focuses on the difference from other embodiments, and the same and similar parts of each embodiment can be referred to each other. As for the system disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and for the related information, please refer to the description of the method part.

本文中应用了具体个例对发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例,基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。In this paper, specific examples are used to illustrate the principle and implementation of the invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea. The described embodiments are only part of the embodiments of the present invention. , not all of the embodiments, based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative work, all belong to the protection scope of the present invention.

Claims (10)

1. a kind of method improving image resolution ratio, which is characterized in that described method includes following steps:
Original image is decomposed using opposite total variance method, obtains cartoon subgraph and texture subgraph;
The cartoon subgraph is amplified using interpolation method, obtains amplified cartoon subgraph;
The texture subgraph is amplified using homotopy Method, obtains amplified texture subgraph;
The amplified cartoon subgraph and the amplified texture subgraph are synthesized, high resolution graphics is obtained Picture.
2. a kind of method improving image resolution ratio according to claim 1, which is characterized in that described to be become using relatively total Poor method decomposes image, obtains cartoon subgraph and texture subgraph;It specifically includes:
Original image is split, the image after being divided;
Using the full variational approach of image for minimizing window function, the image after segmentation is decomposed, obtains cartoon subgraph With texture subgraph.
3. a kind of method improving image resolution ratio according to claim 2, which is characterized in that described to be carried out to original image Segmentation, the image after being divided specifically include:
It is utilized respectively templateWithConvolution is carried out with original image, obtains the gradient map in the directions x As IxWith the gradient image I in the directions yy
Calculate the gradient image modulus value I of the gradient image in the directions x and the gradient image in the directions yg
Wherein, (i, j) is pixel point coordinates;
Pixel using the gradient image modulus value more than first threshold is big by the gradient image modulus value as edge pixel point In second threshold pixel as candidate edge pixel;
Multiple edge pixel points are chosen, respectively using each edge pixel point as starting point, connection and the edge picture Vegetarian refreshments the edge pixel point of direct neighbor or candidate edge pixel on eight directions, obtain several enclosed areas Domain, and the small closed area in the big closed area is deleted, obtain the enclosed area of several not inclusion relations Domain;
Grey level histogram is done to the closed area of each not inclusion relation, and K- is used to the grey level histogram Means methods are split, image after being divided.
4. a kind of method improving image resolution ratio according to claim 2, which is characterized in that
The window function is:
Point (xp,yp) it is Current central pixel point, (xq,yq) it is (xp,yp) full variation pixel,σ2 p,qIt is (xq,yq) and (xp,yp) where pixel value in object variance, cp,qFor multiplier, TgFor third threshold value.
5. a kind of method improving image resolution ratio according to claim 1, which is characterized in that described to use interpolation method pair The cartoon subgraph is amplified, and is obtained amplified cartoon subgraph, is specifically included:
It uses interpolation stencil function into row interpolation the pixel of the cartoon subgraph, obtains the amplified cartoon subgraph Picture;The interpolation stencil function is:
L (x, y)=N1*o(x,y)*d(x,y)*sinc(α*x/σ2)
*sinc(α*y/σ2)*exp(β*(x2+y2)/2σ2) (3)
Wherein, o (x, y) indicates whether the value using current pixel point to estimate the value of unknown pixel point, if current pixel point and Unknown pixel point is in the same object, then o (x, y)=1, that is, utilize current pixel point;Otherwise, o (x, y)=0, i.e., do not utilize Current pixel point, d (x, y) indicate the distance of current pixel point and interpolation template center point;Function sinc (x)=(sin (x)/ x);σ2It is the mean square deviation participated in all pixels point of unknown pixel point interpolation;α and β is constant.
6. a kind of method improving image resolution ratio according to claim 1, which is characterized in that described to use homotopy Method The texture subgraph is amplified, amplified texture subgraph is obtained, specifically includes:
Obtain the image block dictionary of the texture subgraph, described image block dictionary packet offline using K-SVD dictionary learnings method Include high-definition picture block dictionary and low-resolution image block dictionary;
To the texture subgraph, amplified online using described image block dictionary and orthogonal matching tracking method, is obtained just Beginning high-definition picture;
The edged processing of arest neighbors is carried out to the initial high-resolution image;Obtain edged treated full resolution pricture;
Homotopy processing for the first time is carried out to the edged treated high-definition picture, obtains the first high-definition picture;
Second of homotopy processing is carried out to first high-definition picture, obtains amplified high-resolution texture subgraph.
7. it is according to claim 1 it is a kind of improve image resolution ratio method, which is characterized in that it is described to the amplification after Cartoon subgraph and the amplified texture subgraph synthesized, obtain high-definition picture, specifically include:
The amplified cartoon subgraph and the amplified texture subgraph are synthesized using formula (4), obtained High-definition picture;
fH=fc+ft1*G(ft) (4)
Wherein, fHFor high-resolution image, ftAmplified texture image, fcFor amplified cartoon image, G (ft) it is to figure As ftAsk the modulus value of Robert gradients, λ1For constant.
8. a kind of system improving image resolution ratio, which is characterized in that the system comprises:
Picture breakdown module obtains cartoon subgraph and texture for being decomposed to original image using opposite total variance method Image;
Cartoon subgraph amplification module obtains amplified card for being amplified to the cartoon subgraph using interpolation method Logical subgraph;
Texture subgraph amplification module is obtained amplified for being amplified to the texture subgraph using homotopy Method Texture subgraph;
Image synthesis unit, for being closed to the amplified cartoon subgraph and the amplified texture subgraph At acquisition high-definition picture.
9. a kind of system improving image resolution ratio according to claim 8, which is characterized in that described image decomposing module It specifically includes:
Image segmentation submodule, for being split to original image, the image after being divided;
Picture breakdown submodule, for using the full variational approach of image for minimizing window function, being carried out to the image after segmentation It decomposes, obtains cartoon subgraph and texture subgraph.
10. a kind of system improving image resolution ratio according to claim 9, which is characterized in that described image segmentation Module specifically includes:
Convolution unit, for being utilized respectively templateWithConvolution is carried out to original image, is obtained The gradient image I in the directions xxWith the gradient image I in the directions yy
Gradient image magnitude calculation unit, the ladder of the gradient image of gradient image and the directions y for calculating the directions x Spend image modulus value Ig
Wherein, (i, j) is pixel point coordinates;
Pixel selection unit, for the gradient image modulus value to be more than to the pixel of first threshold as edge pixel point, The gradient image modulus value is more than the pixel of second threshold as candidate edge pixel;
Closed area acquiring unit is respectively with each edge pixel point for choosing multiple edge pixel points Initial point, connection and the edge pixel o'clock the edge pixel point of direct neighbor or candidate edge picture on eight directions Vegetarian refreshments obtains several closed areas, and the small closed area in the big closed area is deleted, and acquisition is several not to be had There is the closed area of inclusion relation;
Image segmentation unit does grey level histogram for the closed area to each not inclusion relation, and to the ash Degree histogram is split using K-means methods, the image after being divided.
CN201810120907.XA 2018-02-07 2018-02-07 A method and system for improving image resolution Active CN108364258B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810120907.XA CN108364258B (en) 2018-02-07 2018-02-07 A method and system for improving image resolution

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810120907.XA CN108364258B (en) 2018-02-07 2018-02-07 A method and system for improving image resolution

Publications (2)

Publication Number Publication Date
CN108364258A true CN108364258A (en) 2018-08-03
CN108364258B CN108364258B (en) 2022-09-27

Family

ID=63004989

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810120907.XA Active CN108364258B (en) 2018-02-07 2018-02-07 A method and system for improving image resolution

Country Status (1)

Country Link
CN (1) CN108364258B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI823756B (en) * 2022-06-27 2023-11-21 大陸商威視芯半導體(合肥)有限公司 Image resolution enhancement by value transfer

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020171660A1 (en) * 2001-05-02 2002-11-21 Eastman Kodak Company Block sampling based method and apparatus for texture synthesis
US20160210777A1 (en) * 2015-01-16 2016-07-21 Disney Enterprises, Inc. Image decomposition and path-space motion estimation
CN106127688A (en) * 2016-06-30 2016-11-16 北京大学 A kind of super-resolution image reconstruction method and system thereof
CN106204447A (en) * 2016-06-30 2016-12-07 北京大学 The super resolution ratio reconstruction method with convolutional neural networks is divided based on total variance
CN107341765A (en) * 2017-05-05 2017-11-10 西安邮电大学 A kind of image super-resolution rebuilding method decomposed based on cartoon texture

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020171660A1 (en) * 2001-05-02 2002-11-21 Eastman Kodak Company Block sampling based method and apparatus for texture synthesis
US20160210777A1 (en) * 2015-01-16 2016-07-21 Disney Enterprises, Inc. Image decomposition and path-space motion estimation
CN106127688A (en) * 2016-06-30 2016-11-16 北京大学 A kind of super-resolution image reconstruction method and system thereof
CN106204447A (en) * 2016-06-30 2016-12-07 北京大学 The super resolution ratio reconstruction method with convolutional neural networks is divided based on total variance
CN107341765A (en) * 2017-05-05 2017-11-10 西安邮电大学 A kind of image super-resolution rebuilding method decomposed based on cartoon texture

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
GUODONG JING 等: "Single-image Super-Resolution based on Decomposition and Sparse Representation", 《2010 INTERNATIONAL CONFERENCE ON MULTIMEDIA COMMUNICATIONS》 *
LI XU 等: "Structure Extraction from Texture via Relative Total Variation", 《ACM TRANSACTIONS ON GRAPHICS》 *
周尚波 等: "分数阶偏微分方程在图像处理中的应用", 《计算机应用》 *
闫敬文等: "《压缩感知及应用》", 31 October 2015, 国防工业出版社 *
黄攀峰 等著: "《空间绳系机器人技术》", 31 August 2014, 中国宇航出版社 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI823756B (en) * 2022-06-27 2023-11-21 大陸商威視芯半導體(合肥)有限公司 Image resolution enhancement by value transfer

Also Published As

Publication number Publication date
CN108364258B (en) 2022-09-27

Similar Documents

Publication Publication Date Title
Huang et al. Robust single-image super-resolution based on adaptive edge-preserving smoothing regularization
US9087390B2 (en) High-quality upscaling of an image sequence
US8731337B2 (en) Denoising and artifact removal in image upscaling
Russell et al. Exploiting the sparse derivative prior for super-resolution and image demosaicing
Sandeep et al. Single image super-resolution using a joint GMM method
Mousavi et al. Sparsity-based color image super resolution via exploiting cross channel constraints
WO2015180053A1 (en) Method and apparatus for rapidly reconstructing super-resolution image
Yu et al. Learning to super-resolve blurry images with events
WO2009150882A1 (en) Image registration processing device, region expansion processing device, and image quality improving device
CN106934806A (en) It is a kind of based on text structure without with reference to figure fuzzy region dividing method out of focus
Del Gallego et al. Multiple-image super-resolution on mobile devices: an image warping approach
Walha et al. Resolution enhancement of textual images: a survey of single image‐based methods
Liu et al. Single image super-resolution using a deep encoder–decoder symmetrical network with iterative back projection
Banerjee et al. Super-resolution of text images using edge-directed tangent field
CN106558021B (en) Video enhancement method based on super-resolution technology
Peng et al. Efficient image resolution enhancement using edge-directed unsharp masking sharpening for real-time ASIC applications
Zhao et al. Single depth image super-resolution with multiple residual dictionary learning and refinement
CN108364258B (en) A method and system for improving image resolution
Ning et al. Multi-frame image super-resolution reconstruction using sparse co-occurrence prior and sub-pixel registration
Zhang et al. Bilateral upsampling network for single image super-resolution with arbitrary scaling factors
Xie et al. Feature dimensionality reduction for example-based image super-resolution
CN109951666A (en) Super-resolution restoration method based on surveillance video
Bareja et al. An improved iterative back projection based single image super resolution approach
Luong et al. An image interpolation scheme for repetitive structures
Amintoosi et al. Precise image registration with structural similarity error measurement applied to superresolution

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