CN110298789B - Single image super-resolution reconstruction method based on TV prior - Google Patents

Single image super-resolution reconstruction method based on TV prior Download PDF

Info

Publication number
CN110298789B
CN110298789B CN201910475134.1A CN201910475134A CN110298789B CN 110298789 B CN110298789 B CN 110298789B CN 201910475134 A CN201910475134 A CN 201910475134A CN 110298789 B CN110298789 B CN 110298789B
Authority
CN
China
Prior art keywords
image
resolution
template
super
oms3
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.)
Active
Application number
CN201910475134.1A
Other languages
Chinese (zh)
Other versions
CN110298789A (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.)
Dalian Houren Technology Co ltd
Original Assignee
Liaoning Normal University
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 Liaoning Normal University filed Critical Liaoning Normal University
Priority to CN201910475134.1A priority Critical patent/CN110298789B/en
Publication of CN110298789A publication Critical patent/CN110298789A/en
Application granted granted Critical
Publication of CN110298789B publication Critical patent/CN110298789B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR 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; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/77Retouching; Inpainting; Scratch removal
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)

Abstract

The invention discloses a super-resolution reconstruction method for a single image based on TV prior, which comprises the steps of firstly utilizing a bicubic interpolation algorithm oms3 to perform up-sampling pretreatment on an image with resolution to be improved, extracting TV prior information from the pretreated image according to 28 TV direction templates, and finally introducing the TV prior information into a non-local regression frame, so that texture and edge information of a super-resolution image are better reserved, the phenomenon that the edge information of the image generated by the traditional interpolation algorithm is insufficient is overcome, the super-resolution reconstructed image is obtained, and the super-resolution effect is greatly improved.

Description

Single image super-resolution reconstruction method based on TV prior
Technical Field
The invention relates to a single image super-resolution reconstruction method, in particular to a single image super-resolution reconstruction method based on TV prior, which can enhance the image edge.
Background
At present, image super-resolution reconstruction is a quite active research field, and provides a solution to the low resolution limit from mobile phone imaging, remote sensing imaging and medical imaging. Low quality images are generally caused by low cost sensor capture, narrow band transmission and poor light interference, and when low quality images are widely used for high definition display, the visual analysis and recognition tasks tend to be displayed in higher resolution versions, so that it is necessary to estimate the high resolution image from the content of the low resolution image, i.e. image super-resolution reconstruction.
At present, people only can excavate prior knowledge as much as possible, establish a model, observe the mapping relation from a low-resolution image to a high-resolution image, and design an ultra-resolution reconstruction algorithm on the basis of the generated model. The super-resolution reconstruction method uses a simple linear function to perform up-sampling on an image, such as bilinear interpolation, bicubic interpolation, s-spline and the like, when a larger amplification factor is needed, the reconstruction-based super-resolution reconstruction algorithm cannot achieve a good effect, and has some defects in the aspects of image edge preservation and the like. There are also methods that add gradient prior information to the local image structure, which usually introduce more complex structures (patches, etc.) and auxiliary judgment information (gradients, edge detection information, etc.), but most methods fix a set of invariant model parameters for different image data contents. Therefore, when these methods process natural images having different texture contents, an over-smoothing or burring phenomenon tends to occur.
Total Variation (Total Variation) models were first used to denoise noise contaminated images, and Chan et al, 2002, promoted the TV model to image inpainting. The TV model can effectively decompose structural information and textures in an image without particularly specifying whether the textures are regular or symmetrical. In other words, the method is generic and arbitrary, and is applicable to non-uniform or anisotropic textures. However, there is no report about the application of the TV model to super-resolution reconstruction of single images.
Disclosure of Invention
The invention aims to solve the technical problems in the prior art and provides a single image super-resolution reconstruction method based on TV prior, which can enhance the image edge.
The technical solution of the invention is as follows: a single image super-resolution reconstruction method based on TV prior is sequentially carried out according to the following steps:
step 011: acquiring an Image with resolution to be improved, and recording the Image as Image _ ori;
step 012: performing up-sampling on the Image Image _ ori by a bicubic interpolation oms3 algorithm, and recording an obtained Image as Image _ oms3;
step 021: setting the size of the Image _ oms3 as M × N, moving the Image _ oms3 point by point from the upper left corner of the Image to obtain M × N blocks; starting from the first block, 28 TV direction templates with the size of 3 × 3 are established by taking the pixel to be repaired as the center, and are marked as o (n), n =1,2,3 \8230, 28, wherein n is the serial number of the direction template, 9 pixels in the template are respectively marked as { d1, d2, d3, d4, d5, d6, d7, d8, d9} from top to bottom and from left to right, and then 28 direction templates are respectively marked as: o1{ d1, d5, d6}, o2{ d4, d5, d3},3 medium d3, d5, d6}, o4{ d4, d5, d1}
o5{d6,d5,d7},o6{d4,d5,d9},o7{d6,d5,d9},o8{d4,d5,d7},o9{d2,d5,d7},o10{d2,d5,d9},o11{d2,d5,d3},o12{d2,d5,d1},o13{d8,d5,d1},o14{d8,d5,d3}o15{d8,d5,d9},o16{d8,d5,d7},o17{d4,d5,d2},o18{d6,d5,d2},o19{d4,d5,d8}o20{d6,d5,d8},o21{d1,d5,d3},o22{d1,d5,d7},o23{d3,d5,d9},o24{d7,d5,d9}o25{d4,d5,d6},o26{d2,d5,d8},o27{d1,d5,d9},o28{d7,d5,d3};
Step 022: respectively calculating the mean value of three pixels in each direction template, and recording as s (n);
step 023: calculating the difference between all s (n) and the central pixel value d5, taking the absolute value of the difference as t (n), and selecting the template with the sequence number n corresponding to the minimum value in t (n) as the optimal direction template;
step 031: starting from a first block, establishing a target block with the scale of 3 x 3 by taking a pixel to be repaired as a center, searching 25 search blocks p (m) around the target block, wherein m =1,2,3 \8230 \ 823025, calculating Euclidean distances between the 25 search blocks and the target block to be recorded as D (m), and obtaining the weight of each search block to be recorded as w (m) according to the similarity between the target block and the search blocks;
step 032: in 25 search blocks, calculating the median value of three pixel values in the direction as M (M) according to the optimal direction template;
step 033: and correspondingly multiplying the obtained w (M) and M (M) and then summing to obtain a value which is the pixel value after the current pixel point is repaired.
Step 034: and circularly traversing the whole Image _ images 3, recording the repaired Image as Image _ result and storing the Image _ result.
According to the method, firstly, images to be improved in resolution ratio are subjected to up-sampling pretreatment by utilizing a bicubic interpolation algorithm omoms3, TV prior information is extracted from the pretreated images according to 28 TV direction templates, and finally the TV prior information is introduced into a non-local regression frame, so that texture and edge information of super-resolution images are well reserved, the phenomenon that image edge information generated by a traditional interpolation algorithm is insufficient is overcome, super-resolution reconstructed images are obtained, and the super-resolution effect is greatly improved.
Drawings
Fig. 1 is a schematic diagram of performing bicubic upsampling on an original image according to an embodiment of the present invention.
FIG. 2 is a schematic diagram of an embodiment of the present invention introducing TV prior information into a non-local regression framework.
FIG. 3 is an image to be improved in resolution according to an embodiment of the present invention.
Fig. 4 is an image of the bicubic algorithm after resolution enhancement.
Fig. 5 is an image with improved resolution according to an embodiment of the invention.
Detailed Description
The invention discloses a single image super-resolution reconstruction method based on TV prior, which is sequentially carried out according to the following steps:
a step 011: acquiring an Image with resolution to be improved, and recording the Image as Image _ ori;
step 012: as shown in fig. 1, the Image _ ori is up-sampled by a bicubic interpolation oms3 algorithm, and the obtained Image is marked as Image _ oms3;
step 021: setting the size of the Image _ oms3 as M × N, moving the Image _ oms3 point by point from the upper left corner of the Image to obtain M × N blocks; starting from the first block, as shown in fig. 2, 28 TV direction templates with a scale of 3 × 3 are established with the pixel to be repaired as the center, and are denoted as o (n), n =1,2,3 \8230 \ 823028, where n is the serial number of the direction template, 9 pixels in the template are denoted as { d1, d2, d3, d4, d5, d6, d7, d8, d9} from top to bottom and from left to right, respectively, and then 28 direction templates are denoted as: o1{ d1, d5, d6}, o2{ d4, d5, d3},3 retaining ports d3, d5, d6}, o4{ d4, d5, d1}
o5{d6,d5,d7},o6{d4,d5,d9},o7{d6,d5,d9},o8{d4,d5,d7},o9{d2,d5,d7},o10{d2,d5,d9},o11{d2,d5,d3},o12{d2,d5,d1},o13{d8,d5,d1},o14{d8,d5,d3}o15{d8,d5,d9},o16{d8,d5,d7},o17{d4,d5,d2},o18{d6,d5,d2},o19{d4,d5,d8}o20{d6,d5,d8},o21{d1,d5,d3},o22{d1,d5,d7},o23{d3,d5,d9},o24{d7,d5,d9}o25{d4,d5,d6},o26{d2,d5,d8},o27{d1,d5,d9},o28{d7,d5,d3};
Step 022: respectively calculating the mean value of three pixels in each direction template, and recording as s (n);
step 023: calculating the difference between all s (n) and the central pixel value d5, taking the absolute value of the difference as t (n), and selecting the template with the sequence number n corresponding to the minimum value in t (n) as the optimal direction template;
step 031: starting from a first block, establishing a target block p0 with the scale of 3-3 by taking a pixel to be repaired as the center, searching 25 search blocks p (m) around the target block p0, wherein m =1,2,3 \8230, 25, calculating the Euclidean distance between the 25 search blocks and the target block to be recorded as D (m), and obtaining the weight of each search block to be recorded as w (m) according to the similarity between the target block and the search block;
for the target block p0, the calculation formula of the weight of the search block p1 is:
Figure DEST_PATH_IMAGE001
wherein, the formula of x (p 0, p 1) is: />
Figure 143414DEST_PATH_IMAGE002
In the formula:
Figure DEST_PATH_IMAGE003
is the square of the weighted euclidean distance of block p0 to block p1, device for selecting or keeping>
Figure 388450DEST_PATH_IMAGE004
Used for controlling the corrosion degree of the filter;
step 032: in 25 search blocks, calculating the median value of three pixel values in the direction as M (M) according to the optimal direction template;
step 033: and correspondingly multiplying the obtained w (M) and M (M) and then summing to obtain a value which is the pixel value after the current pixel point is repaired.
Step 034: and circularly traversing the whole Image _ elements 3, and recording the repaired Image as Image _ result and storing the Image.
Fig. 3 is an image to be improved in resolution according to an embodiment of the present invention.
Fig. 4 is an image of the bicubic algorithm after resolution enhancement.
FIG. 5 is an image restored by a single image super-resolution reconstruction method based on TV prior in an embodiment of the present invention.
The following table shows PSNR values after different image reconstructions by the embodiment of the present invention (NT _ oms 3) and bicubic algorithm.
Figure DEST_PATH_IMAGE005
/>

Claims (1)

1. A single image super-resolution reconstruction method based on TV prior is characterized by comprising the following steps in sequence:
step 011: acquiring an Image with resolution to be improved, and recording the Image as Image _ ori;
step 012: performing up-sampling on the Image Image _ ori by a bicubic interpolation oms3 algorithm, and recording an obtained Image as Image _ oms3;
step 021: setting the size of the Image _ oms3 as M x N, moving the Image _ oms3 point by point from the upper left corner of the Image to obtain M x N blocks; starting from the first block, 28 TV direction templates with the size of 3 × 3 are established by taking the pixel to be repaired as the center, and are marked as o (n), n =1,2,3 \8230, 28, wherein n is the serial number of the direction template, 9 pixels in the template are respectively marked as { d1, d2, d3, d4, d5, d6, d7, d8, d9} from top to bottom and from left to right, and then 28 direction templates are respectively marked as: o1{ d1, d5, d6}, o2{ d4, d5, d3},3 retaining ports d3, d5, d6}, o4{ d4, d5, d1}
o5{d6,d5,d7},o6{d4,d5,d9},o7{d6,d5,d9},o8{d4,d5,d7},o9{d2,d5,d7},o10{d2,d5,d9},o11{d2,d5,d3},o12{d2,d5,d1},o13{d8,d5,d1},o14{d8,d5,d3}o15{d8,d5,d9},o16{d8,d5,d7},o17{d4,d5,d2},o18{d6,d5,d2},o19{d4,d5,d8}o20{d6,d5,d8},o21{d1,d5,d3},o22{d1,d5,d7},o23{d3,d5,d9},o24{d7,d5,d9}o25{d4,d5,d6},o26{d2,d5,d8},o27{d1,d5,d9},o28{d7,d5,d3};
Step 022: respectively calculating the mean value of three pixels in each direction template, and recording as s (n);
step 023: calculating the difference between all s (n) and the central pixel value d5, taking the absolute value of the difference as t (n), and selecting the template with the sequence number n corresponding to the minimum value in t (n) as the optimal direction template;
step 031: starting from a first block, establishing a target block with the scale of 3 x 3 by taking a pixel to be repaired as a center, searching 25 search blocks p (m) around the target block, wherein m =1,2,3 \8230 \ 823025, calculating Euclidean distances between the 25 search blocks and the target block to be recorded as D (m), and obtaining the weight of each search block to be recorded as w (m) according to the similarity between the target block and the search blocks;
step 032: in 25 search blocks, calculating the median of three pixel values in the direction as M (M) according to the optimal direction template;
step 033: correspondingly multiplying the obtained w (M) and M (M) and summing to obtain a value which is the pixel value of the current pixel point after restoration;
step 034: and circularly traversing the whole Image _ elements 3, and recording the repaired Image as Image _ result and storing the Image.
CN201910475134.1A 2019-06-03 2019-06-03 Single image super-resolution reconstruction method based on TV prior Active CN110298789B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910475134.1A CN110298789B (en) 2019-06-03 2019-06-03 Single image super-resolution reconstruction method based on TV prior

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910475134.1A CN110298789B (en) 2019-06-03 2019-06-03 Single image super-resolution reconstruction method based on TV prior

Publications (2)

Publication Number Publication Date
CN110298789A CN110298789A (en) 2019-10-01
CN110298789B true CN110298789B (en) 2023-04-14

Family

ID=68027535

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910475134.1A Active CN110298789B (en) 2019-06-03 2019-06-03 Single image super-resolution reconstruction method based on TV prior

Country Status (1)

Country Link
CN (1) CN110298789B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109064402A (en) * 2018-07-27 2018-12-21 四川大学 Based on the single image super resolution ratio reconstruction method for enhancing non local total variation model priori

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103700062B (en) * 2013-12-18 2017-06-06 华为技术有限公司 Image processing method and device

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109064402A (en) * 2018-07-27 2018-12-21 四川大学 Based on the single image super resolution ratio reconstruction method for enhancing non local total variation model priori

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
全局和局部结构内容自适应正则化的单幅图像超分辨模型;唐松泽;肖亮;黄伟;刘鹏飞;中国图象图形学报;20(1);全文 *

Also Published As

Publication number Publication date
CN110298789A (en) 2019-10-01

Similar Documents

Publication Publication Date Title
CN111080724B (en) Fusion method of infrared light and visible light
Lu et al. Underwater image super-resolution by descattering and fusion
Kiku et al. Beyond color difference: Residual interpolation for color image demosaicking
CN109146788B (en) Super-resolution image reconstruction method and device based on deep learning
WO2022110638A1 (en) Human image restoration method and apparatus, electronic device, storage medium and program product
CN108154474B (en) Super-resolution image reconstruction method, device, medium and equipment
CN110163237B (en) Model training and image processing method, device, medium and electronic equipment
CN112767290B (en) Image fusion method, image fusion device, storage medium and terminal device
CN106920214B (en) Super-resolution reconstruction method for space target image
Mousavi et al. Sparsity-based color image super resolution via exploiting cross channel constraints
CN107680043B (en) Single image super-resolution output method based on graph model
CN110557584A (en) image processing method and device, and computer readable storage medium
CN109785233B (en) Image super-resolution reconstruction method
CN108416736B (en) Image super-resolution reconstruction method based on secondary anchor point neighborhood regression
WO2022198381A1 (en) Imaging processing method and device
Zhang et al. Deep spatial adaptive network for real image demosaicing
Liu et al. Single image super-resolution using a deep encoder–decoder symmetrical network with iterative back projection
WO2020078102A1 (en) Image enhancement method and apparatus, and computer-readable storage medium
Chen et al. Dynamic degradation intensity estimation for adaptive blind super-resolution: A novel approach and benchmark dataset
CN109064394B (en) Image super-resolution reconstruction method based on convolutional neural network
CN110298789B (en) Single image super-resolution reconstruction method based on TV prior
CN106558021B (en) Video enhancement method based on super-resolution technology
CN116403226A (en) Unconstrained fold document image correction method, system, equipment and storage medium
Wang et al. Super resolution for compressed screen content video
Cao et al. Super-resolution algorithm through neighbor embedding with new feature selection and example training

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
TR01 Transfer of patent right

Effective date of registration: 20231120

Address after: Room 2401, 2402, 2403, 2407, 24/F, Building B, No. 6, Yiyang Road, Qixianling, High-tech Industrial Park, Dalian, Liaoning, 116000

Patentee after: Dalian Houren Technology Co.,Ltd.

Address before: No. 116500, Shahekou Road, Dalian City, Liaoning Province

Patentee before: LIAONING NORMAL University

TR01 Transfer of patent right