Systems and method of training networks for real-world super resolution with unknown degradations

A method and apparatus are provided. The method includes generating a dataset for real-world super resolution (SR), training a first generative adversarial network (GAN), training a second GAN, and fusing an output of the first GAN and an output of the second GAN.

FIELD

The present disclosure relates generally to image super resolution, and more particularly, to system and method for designing efficient super resolution deep convolutional neural networks by cascade network training, cascade network trimming, and dilated convolutions.

BACKGROUND

Super resolution imaging generates a high resolution (HR) image from a low resolution (LR) image. Super resolution (SR) imaging has wide applicability, from surveillance and face/iris recognition to medical image processing, as well as the straightforward improvement of the resolution of images and video. Many algorithms/systems have been proposed for performing SR, from interpolations (Li, Xin and Orchard, Michael, New edge-directed interpolation. IEEE Transactions on Image Processing (TIP), vol. 10, issue 10, pp. 1521-1527 (Oct. 2001), which is incorporated by reference in its entirety), contour features (Tai, Yu-Wing; Liu, Shuaicheng; Brown, Michael; and Lin, Stephen,Super resolution using edge prior and single image detail synthesis.2010 IEEE Int'l Conference on Computer Vision and Pattern Recognition (CVPR), pp. 2400-2407, which is incorporated by reference in its entirety), and statistical image priors (Kim, Kwang In and Kwon, Younghee.Single-image super-resolution using sparse regression and natural image prior. IEEE Transactions on Pattern Analysis and Machine Intelligence (TPAMI), vol. 32, no. 6, pp. 1127-1133 (Jan. 2010), which is incorporated by reference in its entirety), to example-based methods which learn from a dictionary of patches, such as neighbor embedding (Chang, Hong; Yeung, Dit-Yan; and Xiong, Yimin,Super-resolution through neighbor embedding,2004 CVPR, pp. 275-282, which is incorporated by reference in its entirety) and sparse coding (Yang, Jianchao; Wright, John; Huang, Thomas; and Ma, Yi,Image super-resolution via sparse representation. IEEE TIP, vol. 19, no. 11, pp. 2861-2873 (Nov. 2010), which is incorporated by reference in its entirety).

Recently, convolutional neural networks (CNNs) have provided a significant improvement in SR accuracy. See, e.g., Dong, Chao; Loy, Chen Change; He, Kaiming; and Tang, Xiaoou,Learning a deep convolutional network for image super-resolution,2014 European Conference on Computer Vision (ECCV), pp. 184-199 (hereinafter, “Dong et al. 2014”), which is incorporated by reference in its entirety. Sometimes referred to as “SRCNNs” (i.e., super-resolution convolutional neural networks), their accuracy can be limited by a small structure, e.g., 3-layers, and/or small context reception field. In response, researchers have proposed increasing the size of SRCNNs, but most proposals use a prohibitively large number of parameters, and many of the SRCNNs under discussion cannot be executed in real-time. Due to the large network sizes being proposed, it can be very difficult to even guess at the appropriate training settings, i.e., learning rate, weight initialization, and weight decay. As a result, training may not converge at all or fall into a local minimum.

SUMMARY

Accordingly, the present disclosure has been made to address at least the problems and/or disadvantages described herein and to provide at least the advantages described below.

According to one embodiment, a method includes generating a dataset for real-world SR, training a first generative adversarial network (GAN), training a second GAN, and fusing an output of the first GAN and an output of the second GAN.

According to one embodiment, an apparatus includes one or more non-transitory computer-readable media, at least one processor which, when executing instructions stored on the one or more non-transitory computer-readable media, performs the steps of generating a dataset for real-world SR, training a first GAN, training a second GAN, and fusing an output of the first GAN and an output of the second GAN

DETAILED DESCRIPTION

Hereinafter, embodiments of the present disclosure are described in detail with reference to the accompanying drawings. It should be noted that the same elements are designated by the same reference numerals although they are shown in different drawings. In the following description, specific details such as detailed configurations and components are merely provided to assist in the overall understanding of the embodiments of the present disclosure. Therefore, it should be apparent to those skilled in the art that various changes and modifications of the embodiments described herein may be made without departing from the scope of the present disclosure. In addition, descriptions of well-known functions and constructions are omitted for clarity and conciseness. The terms described below are terms defined in consideration of the functions in the present disclosure, and may be different according to users, intentions of the users, or custom. Therefore, the definitions of the terms should be determined based on the contents throughout the specification.

The present disclosure may have various modifications and various embodiments, among which embodiments are described below in detail with reference to the accompanying drawings. However, it should be understood that the present disclosure is not limited to the embodiments, but includes all modifications, equivalents, and alternatives within the scope of the present disclosure.

The terms herein are merely used to describe various embodiments of the present disclosure but are not intended to limit the present disclosure. Singular forms are intended to include plural forms unless the context clearly indicates otherwise. In the present disclosure, it should be understood that the terms “include” or “have” indicate existence of a feature, a number, a step, an operation, a structural element, parts, or a combination thereof, and do not exclude the existence or probability of addition of one or more other features, numerals, steps, operations, structural elements, parts, or combinations thereof.

Various embodiments may include one or more elements. An element may include any structure arranged to perform certain operations. Although an embodiment may be described with a limited number of elements in a certain arrangement by way of example, the embodiment may include more or less elements in alternate arrangements as desired for a given implementation. It is worthy to note that any reference to “one embodiment” or “an embodiment” means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment. The appearance of the phrase “one embodiment” (or “an embodiment”) in various places in this specification does not necessarily refer to the same embodiment.

This disclosure provides a new approach, or more accurately, several new techniques, for creating an SRCNN. Herein, the term “cascade trained super resolution convolutional neural network” (CT-SRCNN) may refer to all of the new techniques described herein together, or to one or more of the new techniques, which should be made clear by the context in which the term is used. Different from existing approaches that train all the layers from the beginning with unsupervised weight initialization, CT-SRCNN starts training with a small network (e.g., 3 layers). New layers are gradually inserted into the network when the current network cannot adequately reduce the training error.

With this “cascade training” strategy, convergence is made easier, and the accuracy is consistently increased when more layers are used. But while the depth increases, the relative complexity of the network does not, because of the nature of the new layers. More specifically, all the weights of the new layers in CT-SRCNN are randomly initialized, and the learning rate is fixed. This is a great advantage compared to approaches which need to spend a great deal of time and resources tuning the parameters. One specific example of a CT-SRCNN with 13 layers (as shown and discussed further below), the accuracy is competitive with the state-of-the-art image SR networks, while having an execution speed more than 5 times faster, and using only 115th of the parameters.

In this disclosure, “cascade network trimming” is described, which further refines the CT-SRCNN model by reducing the storage and computational complexities, as well as another method to further improve the efficiency of super resolution deep convolutional neural networks by deploying a form of “dilated convolution,” instead of performing the complete conventional convolutional calculations, which may further reduce the CT-SRCNN model complexity.

The rest of the disclosure discusses these three different schemes/features of the CT-SRCNN in order:

II. Cascade network trimming; and

Although these three methods/techniques are discussed in the context of the CT-SRCNN, each method/technique could be applied individually or severally to other SR schemes or CNN networks, as would be understood by one of ordinary skill in the art.

FIG.1illustrates an exemplary block diagram of a method for constructing a cascade trained super resolution convolutional neural network (CT-SRCNN), according to one embodiment.

At110, the training set is prepared, meaning a set of low resolution (LR) images with corresponding high resolution (HR) images by which the CT-SRCNN “learns” a model to use when attempting to create high resolution images from low resolution images. In this embodiment, at120, each LR image is bicubic upsampled and the LR/HR patches cropped in preparation for training. See, e.g., Dong et al. 2014 and Dong, Chao; Loy, Chen Change; He, Kaiming; and Tang, Xiaoou,Image super-resolution using deep convolutional networks, IEEE Transactions on Pattern Analysis and Machine Intelligence (TPAMI), vol. 38, no. 2, pp. 295-307 (Feb. 2016) (hereinafter, “Dong et al. 2016a”), which is incorporated by reference in its entirety, for more details concerning this step. As would be understood by one of ordinary skill in the art, there are a variety of pre-training preparation techniques and this disclosure is not limited to this bicubic upsampling and LR/HR patching as a pre-training preparation technique.

At130, cascade training is performed in accordance with the present disclosure. Embodiments of cascade training according to specific embodiments of the present disclosure are described below. At140, cascade network trimming is performed in accordance with the present disclosure. Embodiments of network trimming according to specific embodiments of the present disclosure are further described below. At150, the process is complete and the CT-SRCNN system is ready for real world use.

Although these different processes (i.e., cascade training and cascade network trimming) are described and shown inFIG.1as separate and distinct stages/steps, there may be overlap between these functions in actual implementations in accordance with the present disclosure.

I. Cascade Training

FIG.2illustrates an exemplary diagram of cascade training, according to one embodiment. At205, the process of training begins.

At210, the training starts at stage i=1. The fledgling network starts with b number of layers, and c number of layers are added in each stage where the training error converges (220) or remains higher than a threshold (250). Thus, at each training stage i, a CNN with c*(i−1)+b layers is trained. When stage i=1, the CNN with the first b number of layers is trained. After stage i=1, cascade training starts to add intermediate layers to the b number of layers, specifically c number of layers at a time, as necessary.

At220, it is determined whether the network has started to converge, e.g., whether the training error has stopped decreasing by a certain amount (from the previous stage). If it has (i.e., the CNN is converging), c number of intermediate layers are added at230, and the next iteration begins at240(i=i+1). During this iterative process, the new layers may be set to any arbitrary weighting, as the intermediate layers will have no influence on the weight matrix sizes of the other layers. Indeed, all existing layers inherit their previous weight matrix. This cascade training iterative process continues, making the CNN deeper and deeper, until the training error is smaller than a threshold at250, and then the CNN model is output at255.

FIGS.3A and3Billustrate some of the differences between cascade training and existing training methods.

InFIG.3A, an example of the flowchart inFIG.2is shown. InFIG.3A, b number of layers equals three, as shown at the top (310), which represents the first CNN to be trained, and the number of layers c added in each stage is one. Each new layer has its weights set randomly, while each pre-existing layer inherits its weights from the previous stage. With every newly inserted intermediate layer, the CNN becomes deeper. At each stage, the deeper CNN is trained again. Since most of the weights are inherited from the previous stage, the continuous re-training is relatively easy, even with a fixed learning rate.

Existing methods, however, as shown inFIG.3B, start with a “complete” set of layers which need to be tuned at the same time. Training all of the layers at the same time as shown inFIG.3Bis far more complex than the scheme shown inFIG.3Adue to the slow convergence, where cascade training trains shallower networks until convergence, incrementally inserts layers with random weights while keeping previously trained layers intact, and retrains the whole network until a deeper network converges. Moreover, cascade training can simply fix the learning rate and generate new layers with random weights.

FIGS.4A and4Billustrate a beginning CNN and an ending CNN, respectively, after cascade training according to one embodiment.

Let x denote an interpolated LR image and y denote its matching HR image. Given a training set {(xi, yi), i=1, . . . , N} with N samples, the goal for the CT-SRCNN is to learn a model g that predicts the HR output 9=g(x). During training, mean square error

(MSE)⁢12⁢∑i=1Nyi-yˆi
is minimized over the training set.

InFIG.4A, cascade training starts from a 3-layer model (b=3). The first layer (410) consists of 64 9×9 filters, and the second (413) and the third layer (415) consist of 32 5×5 filters. All the weights (of new layers) are randomly initialized by a Gaussian function with σ=0.001, and all convolutions have stride one. “Stride” is one of the hyperparameters of a convolutional layer, and controls how the depth columns around the spatial dimensions (width and height) are allocated-to put it another way, stride indicates how the filter convolves around the input volume, namely, “stride one” indicates that the filter convolves around the input volume one pixel at a time, “stride two” indicates the filter convolves two pixels at time, etc. See, e.g., Definition of “Convolutional neural network,” downloaded on Jun. 5, 2017 from Wikipedia at https://en.wikipedia.org/wiki/Convolutional neural network; “A Beginner's Guide to Understanding Convolutional Networks—Part 2,” downloaded on Jun. 5, 2017 from https://adeshpande3.github.io/A-Beginner %27s-Guide-To-Understanding-Convolutional-Neural-Networks-Part-2/; both of which are incorporated by reference in their entireties.

Returning toFIG.4A, when the MSE of the current stage stops decreasing significantly, e.g., the error decreases less than 3% in an epoch, the training goes to the next stage. See, e.g. step220ofFIG.2. To accelerate the training in this embodiment, two new layers are inserted into the network for each stage (i.e., c=2 in step230inFIG.2). Thus, the training starts from 3 layers, as shown atFIG.4A, and then proceeds to 5 layers, 7 layers, . . . , and finally 13 layers after five (5) stages. Each new layer consists of 32 3×3 filters. This size ensures a smaller network even when the CNN is becoming progressively deeper. The new intermediate layers are inserted immediately before the last 32 5×5 filters layer415. The weights from any layer existing in the preceding stage inherits the weights from the previous stage, and the weights of the two new layers are always randomly initialized (Gaussian distribution with σ=0.001). Since new convolutional layers will reduce the size of the feature map, 2 pixels are zero-padded in each new intermediate 3×3 layer. As a result, all the stages in cascade training have the same size of the output, so that the training samples can be shared.

As a network goes deeper, it usually becomes more difficult for the training with existing methods to converge. For example, the SRCNN in Dong et al. 2016a failed to show superior performance with more than three layers. In Kim, Jiwon; Lee, Jung Kwon; and Lee, Kyoung Mu,Accurate image super-resolution using very deep convolutional networks,2016 CVPR, pp. 1646-1654, which is incorporated by reference in its entirety (hereinafter, “VDSR”), a high initial learning rate is tuned and gradually decreased. But when using a large diverse training set (e.g., more than 30 million patches from 160,000 images), the high learning rate does not work well. A potential reason for this is that the high learning rate leads to vanishing/exploding gradients.

In CT-SRCNN, only a few weights are randomly initialized in each stage, so the convergence is relatively easy. A fixed learning rate 0.0001 for all layers in CT-SRCNN without any decay is also feasible. In order to accelerate the training, only the first stage need be changed, e.g., the learning rate of the first stage can be set to 0.001. In experiments/simulations, the 13-layer CT-SRCNN like the one inFIG.4Bhas already achieved state-of-the-art accuracy, while using many less parameters compared to other networks such as VDSR or Kim, Jiwon; Lee, Jung Kwon; and Lee, Kyoung Mu,Deeply-recursive convolutional network for image super-resolution,2016 CVPR, pp. 1637-1645, which is incorporated by reference in its entirety (hereinafter, “DRCN”). To the contrary, direct training of a randomly initialized deeper network requires a lot of effort in parameter tuning to ensure best convergence in these other networks, even though experiments have shown these networks may fail to converge with acceptable error.

As shown in Table 1 below, when two image qualities metrics, the peak signal to noise ratio (PSNR) and the structure similarity measure (SSIM), are measured, it can be seen that the CT-SRCNN achieves better quality and faster speed. Moreover, the CT-SRCNN retrieves more details compared to VDSR and DRCN.

Given an L-layer in CNN, assume the ithlayer has ηi−1input channels, a ki×kiconvolution kernel, and ηifilters. The number of parameters in the ithlayer is ηi−1×ηi×ki×ki. The bias term is ignored in this calculation. Then the overall number of parameters is Σi=1Lηi−1×ηi×ki×ki. Thus, for example, in a 3-layer CT-SRCNN with 64-32-1 filters in each layer, η0=1, η1=64, η2=32,η3=1, k1=9,k2=5,k3=5, so the overall number of parameters is 1×64×9×9+64×5×5×32+1×32×5×5×1=57,184.

PSNR/SSIM are utilized to measure the image reconstruction quality. PSNR is the ratio between the maximum possible power of an image pixel and the power of corrupting noise that affects the fidelity. It is calculated as

PSNR=2⁢0⁢log1⁢0⁢2⁢5⁢5M⁢S⁢E,
where the MSE is calculated between the ground truth and a reconstructed image (SR output). The larger the PSNR, the better the image quality. The maximum value of PSNR is infinite. See, e.g., definition of “Peak signal-to-noise ratio,” downloaded on Jun. 27, 2017 from Wikipedia at https://en.wikipedia.org/wiki/Peak signal-to-noise ratio, which is incorporated by reference in its entirety.

SSIM is a perception-based model that considers image degradation as perceived change in structural information, while also incorporating the luminance masking and contrast masking. It is more consistent with human vision than PSNR. SSIM is calculated as

SSIM=20⁢log10⁢(2⁢μx⁢μy+c1)⁢(2⁢σx⁢y+c2)(μx2+μy2+c1)⁢(σx2+σy2+c2),
where x is the reconstructed image, y is the reference image (ground truth), μ is the mean, σ is the variance, σxyis the covariance between x and y, c1=6.5025, and c2=58.5225. SSIM lays between [0,1]. If x is a perfect copy of y, the SSIM will be 1. See, e.g., definition of “Structural Similarity,” downloaded on Jun. 27, 2017 from Wikipedia at https://en.wikipedia.org/wiki/Structural similarity, which is incorporated by reference in its entirety.

Most neural networks have redundancy. Removing such redundancy clearly improves efficiency. In embodiments of the present disclosure, a major number of filters and/or weights may be removed from certain layers with a minor loss in accuracy.

This technique/approach (cascade network trimming) can be used with the cascade training described above, or can be used independent of cascade training. Given a deep convolutional neural network with acceptable accuracy or performance, techniques/approaches for reducing network size, computational complexity, and/or processing speed while keeping the network depth the same and not degrading the accuracy are always needed.

Similar to cascade training, cascade network trimming also includes an iterative process. In each stage, filters are trimmed from only d layers, which means that, for an L-layer network, the (L−(i−1)d−1)th layer to (L-id)th layer are trimmed in stage i. For example, when trimming d=2 layers from a 13-layer CT-SRCNN, the 12thand 11thlayers are trimmed in the first stage i=1, and then the network is fine-tuned. When it converges, the second stage i=2 begins with trimming the 9thand 10thlayers. This procedure is iteratively repeated until all of the layers are trimmed. Although the 13thlayer is ignored in the above procedure, the procedure may also be considered as trimming 12thand 13thlayer in the first stage, and trimming 10thand 11thlayer in the second stage, etc.

FIG.5illustrates an exemplary diagram of cascade network trimming, according to one embodiment. At505, the process of trimming begins with a trained CNN with L layers.

At510, the trimming starts at stage i=1. As mentioned above, only d layers of the total L-layer CNN are trimmed in a stage. Thus, the (L−(i−1)d−1)th layer to (L-id)th layer are trimmed in stage i at510. At520, fine tuning is performed. At530, it is determined whether the training error has stopped decreasing by a certain amount (from the previous stage). If it has, it is determined whether the total number of stages multiplied by the layers trimmed per stage is greater than or equal to the total number of layers at540(“(id>=L)?”). If the training error has not stopped decreasing at530, the method returns to fine tuning at520.

If it is determined that the total number of stages multiplied by the layers trimmed per stage is greater than or equal to the total number of layers at540(“(id>=L)?”), the process ends and the trimmed CNN model is output at565. If it is determined that the total number of stages multiplied by the layers trimmed per stage is less than the total number of layers at540(“(id>=L)?”), the method begins the next stage at550(“i=i+1”).

FIGS.6A and6Billustrate some of the differences between network trimming methods in accordance with one embodiment.

InFIG.6A, one layer of the CNN is trimmed per stage, and fine tuning is performed between each stage, in accordance with an embodiment of the present disclosure. By contrast, all of the layers of the CNN inFIG.6Bare both fine-tuned and trimmed at the same time. Tuning and trimming all of the layers at the same time as shown inFIG.6Bis far more complex than the scheme shown inFIG.6A.

Cascade network trimming is done by trimming whole filters from the layers. To recover any lost accuracy, trimming is done layer by layer, with fine-tuning till convergence after each trimmed layer or group of layers.

As shown inFIG.7, once a filter is trimmed, the adjacent layer will also be influenced. InFIG.7, a filter710(block of dotted lines) is trimmed from the ith layer, ηi=ηi−1, some weights (indicated by the dotted lines within the filters) in the (i+1)th layer will also be trimmed. So trimming the filter in the ith layer will reduce the computational cost for both the ith and the (1+1)th layer. In CNN, the number of input channels of the (i+1)th layer is equal to the number of filters (output channel) of the ith layer.

InFIG.7, assume there are ηi=4 filters and =5 input channels in the ith layer, and ηi+1=10 filters and ηi=4 input channels in the (i+1)th layer before the trimming. If filter710is trimmed from the ith layer, the trimmed ηiwill be reduced to 3, and the ηi+1 is still 10. The slices720in the (i+1)th layer are the trimmed weights, which correspond to the multiplications. As mentioned in the last section, there will be ηi−1×ki×ki×ηi×wi×himultiplications in the ith layer, and ηi×ki+1×ki+i×ηi+i×wi+1×hi+1multiplications in the (i+1)th layer. Since ηiis reduced, the number of multiplications in both the ith layer and (i+1)th layer are also reduced.

An appropriate criteria is used to decide which filters are to be trimmed. In this embodiment, a measurement of relative importance is used. More specifically, the relative importance Ri,jof the jth filter in the ith layer is defined by the square sum of all the weights in jth filter, where Wi,jis the weights matrix of the jth filter in the ith layer, as shown in Equation (1):

Accordingly, the filters with the smallest Ri,jare removed. As discussed above, when filter710is trimmed from the ith layer, some weights720in the (i+1)th layer will also be trimmed, resulting in W′i+1,j. Thus, when calculating Ri+1,j, either the non-trimmed weights Wi+1,j(also referred to as “independent trimming”) are used in Equation (3), or the trimmed weights W′i+1,jare used in Equation (2):

The algorithm below provides an exemplary high-level description of the iterative process for trimming filters from the layers.

With different rates/thresholds ∈weightsand ∈filters,i, different trimmed models may be created. Since filter trimming influences the adjacent layers, fine-tuning will be needed to retrieve the accuracy in most cases where filter trimming is used. By contrast, weight pruning has a relatively smaller influence. With an appropriate trimming rate (e.g., less than 0.2), the accuracy will not decrease much even without fine-tuning.

Dilated convolution, also known as á trous convolution, is a type of convolution which was originally developed for wavelet decomposition (see Holschneider, M.; Kronland-Martinet, R.; Morlet, J.; and Tchamitchian, Ph.,A Real-Time Algorithm for Signal Analysis with the Help of the Wavelet Transformin WAVELETS: TIME-FREQUENCY METHODS AND PHASE SPACE, J. M. Combes et al., eds., pp. 286-297 (1987), which is incorporated by reference in its entirety), but has been applied to semantic segmentation, particularly in order to get dense features (see, e.g., Yu, Fisher and Koltun, Vladlen,Multi-scale context aggregation by dilated convolutions,2016 Int'l Conference on Learning Representations (ICLR) (hereinafter, “Yu et al. 2016”), which is incorporated by reference in its entirety).

In a purely convolutional network composed of convolutions layers without pooling, the receptive field of units can only grow linearly layer by layer because the feature maps are generated based on convolving adjacent pixels from the input. A feasible way to increase the receptive field is to convolve the input pixels from a larger region. This is similar to using a ‘dilation kernel’ in dilation convolution instead of using the conventional dense kernel for conventional convolution.

Suppose F is a discrete function, K is a convolution kernel, and the dilated convolution*d is a generalized version of typical convolution, as defined by Equation (3) below, where d is the dilation factor. The conventional convolution is a simple 1-dilated convolution (Le., when d=1).

One advantage of applying dilated convolution in a CNN is that the dilated version has a larger reception field, as illustrated inFIGS.8A and8B. The dilated convolutional filter is obtained by upsampling the original filter, i.e., by inserting zeros between its elements. Hence, by design the dilated filter has a structured pattern of zero elements. Compared to weight pruning, where the zero elements have a random pattern and location, dilated filters have a structured pattern for the zero weights, and are much more useful for reducing the computational complexity in hardware and software. Hence, for super resolution, embodiments of the present disclosure deploy the dilated filters differently than their typical usage, by keeping the same receptive field and instead using it to reduce the computational complexity in comparison to the non-dilated filter with the same receptive field.

FIGS.8A and8Billustrate some of the differences between dilated convolution in accordance with one embodiment and conventional convolution, respectively. InFIG.8B, conventional convolution is performed with stride two, while, inFIG.8A, there is a 2-dilated convolution (meaning that the multiplication-and-accumulation operation in the convolution is applied every 2 pixels, instead of every pixel) with stride one according to an embodiment of the present disclosure. AlthoughFIGS.8A and8Bhave the same feature map size (with padding for the dilated version), the reception field of the 2-dilated feature map is larger compared to the convolutional one. In a CNN, the input and output are 2-D feature maps soFIGS.8A or8Bare on the x-direction or y-direction only.

FIG.8Billustrates an example of a conventional convolution with a size 3 kernel and stride 2, where the input is a 7-pixel signal (represented by 7 circles). InFIG.8B, every 3 adjacent pixels are convolved (as indicated by the connecting lines) with the kernel and then generate an output (a square) of the feature map, beginning with the 1stto 3rdpixels (the lined circles) and the first output (the lined square) of the feature map. The next convolution inFIG.8Bare the 3rdto 5thpixels because the stride is 2, and the next output (the black square) of the feature map consists of elements, with receptive field 3.

By contrast,FIG.8Aillustrates an example of a 2-dilated convolution with a size 3 kernel and stride 1. In d-dilated convolution, the convolution is applied every d pixels. So the first output (the lined square) of the feature map is generated by convolving the 1st, 3rd, and 5thpixels (lined circles) with the 3×3 kernel. Then the next output (the black square) is generated by convolving the 2nd, 4th, and 6thpixels.

In an embodiment where all the layers in the CNN are convolutions with stride one, the dilated convolution may be applied in a different way. Given a k×k convolution kernel with stride one, the reception field of the resulting feature map is k×k. If 2-dilated convolution is used, the reception field of the resulting feature map is (2k−1)×(2k−1). For example, the 9×9 1-dilated layer410and a 5×5 1-dilated layer413of the CT-SRCNN inFIGS.4A and4Bmay be replaced by a 5×5 2-dilated layer and 3×3 2-dilated layer, respectively, instead. The resulting network will have the same size reception field, but less parameters due to the smaller kernel size.

Accordingly, in one embodiment, once a CT-SRCNN with a 9×9 1-dilated layer and two 5×5 1-dilated layers is trained, those layers may be replaced by a 5×5 2-dilated layer and two 3×3 2-dilated layers before fine tuning is performed. Unlike Yu et al. 2016, a dilated CT-SRCNN according to an embodiment of the present disclosure does not need any zero padding in the dilated layer.

As mentioned above, many researchers are attempting to increase the accuracy and efficiency of SRCNNs by, for example, using more layers (e.g., VDSR) or a deeply recursive structure (e.g., DRCN). Other researchers similarly propose to use more complicated networks. Wang, Zhaowen; Liu, Ding; Yang, Jianchao; Han, Wei; and Huang, Thomas,Deep networks for image super-resolution with sparse prior,2015 IEEE Int'l Conference on Computer Vision (ICCV), pp. 370-378, which is incorporated herein by reference, integrated a sparse representation prior with feed-forward network based on the learned iterative shrinkage and thresholding algorithm. VDSR increased the number of layers to 20 and used small filters and a high learning rate with adjustable gradient clipping; the same group also designed a deep recursive CNN with recursive-supervision and skip-connection in DRCN. Dahl, Ryan; Norouzi, Mohammad; and Shlens, Jonathon,Pixel Recursive Super Resolution, arXiv 1702.00783 [22 Mar. 2017], which is incorporated herein by reference, combined the ResNet with a Pixel Recursive Super Resolution, which showed promising results on face and bed SR where super resolution is applied to bed images).

Others prefer to use perception loss instead of the mean square error (MSE) for the training error, which is closer to natural texture and human vision. Sønderby, Casper; Caballero, Jose; Theis, Lucas; Shi, Wenzhe; and Huszár, Ferenc,Amortised MAP Inference for Image Super-resolution, arXiv 1610.04490 [21 Feb. 2017], which is incorporated herein by reference, introduced a method for amortised MAP inference, which calculated the MAP estimation directly using CNN. Johnson, Justin; Alahi, Alexandre; and Fei-Fei, Li,Perceptual losses for real-time style transfer and super-resolution,2016 ECCV, pp. 694-711, which is incorporated herein by reference, proposed the use of perceptual loss functions for training feedforward networks for image transformation tasks. Ledig, Christian, et al.,Photo-realistic single image super-resolution using a generative adversarial network, arXiv 1609.04802 [13 Apr. 2017], which is incorporated herein by reference, employed a very deep residual network (ResNet), and further presented the super resolution generative adversarial network (SRGAN) to obtain the images similar to natural texture.

However, although the works listed above improved the accuracy of the SR system, the improved accuracy was at the cost of having more layers/parameters and/or more difficult hyperparameter tuning procedures. In other words, any advance in accuracy was counter-balanced by extreme increases in complexity.

Other researchers focused on improving efficiency by extracting the feature maps in LR space and training using upscaling filters. Shi, Wenzhe, et al.,Real-time Single Image and Video Super-Resolution Using an Efficient sub pixel convolutional neural network,2016 CVPR, pp. 1874-1883, which is incorporated herein by reference, introduced an efficient subpixel convolutional layer which learned an array of upscaling filters to upscale the LR feature maps into the HR output. Dong, Chao; Loy, Chen Change; and Tang, Xiaoou,Accelerating the super-resolution convolutional neural network.2016 ECCV, pp. 391-407, which is incorporated herein in its entirety (hereinafter, “Dong et al. 2016b”), re-designed the SRCNN by adding smaller filters, a deconvolution layer, and feature space shrinkage to accelerate the speed without losing the accuracy.

However, due to the use of the upscaling layer, the patch size and context reception field of these networks will be relatively small. As a result, the accuracy is relatively lower compared to extracting feature maps from upsampled LR space.

By contrast, the CT-SRCNN described herein can go deeper, thereby achieving high accuracy, without the heavy tuning of parameters. The network size of the CT-SRCNN is much smaller compared to the state-of-the-art solutions, such as those listed above. The CT-SRCNN can also process 20-25 frames/second of video with a resolution of 720×480 in a single GPU. This efficiency could be further enhanced by network trimming and dilated convolution.

In this disclosure, a cascade training method is described which trains deep CNN for super resolution with both high accuracy and efficiency. The cascade training ensures that the network might consistently go deeper with a relatively smaller size. The network trimming and dilated convolution described herein further reduce the network complexity. The experimental results on benchmark image and video datasets show that the disclosed method herein achieves competitive performance compared to other state-of-the-art solutions, but at much higher speeds.

Although described in the framework of image super-resolution, the techniques described herein can be generalized any type of CNN for any type of purpose, such as denoising, or image restoration.

FIG.9illustrates an exemplary diagram of the present apparatus, according to one embodiment. An apparatus900includes at least one processor910and one or more non-transitory computer readable media920. The at least one processor910, when executing instructions stored on the one or more non-transitory computer readable media920, performs the steps of training a CNN having three or more layers; performing cascade training on the trained CNN to add one or more intermediate layers until a training error is less than a threshold; and performing network trimming of the CNN output from the cascade training. Moreover, the one or more non-transitory computer-readable media920stores instructions for the at least one processor910to perform the steps of training a CNN having three or more layers; performing cascade training on the trained CNN to add one or more intermediate layers until a training error is less than a threshold; and performing network trimming of the CNN output from the cascade training.

FIG.10illustrates an exemplary flowchart for manufacturing and testing the present apparatus, according to one embodiment.

At1050, the apparatus (in this instance, the chipset described above) is manufactured, including at least one processor and one or more non-transitory computer-readable media. When executing instructions stored on the one or more non-transitory computer readable media, the at least one processor performs the steps of training a CNN having three or more layers; performing cascade training on the trained CNN to add one or more intermediate layers until a training error is less than a threshold; and performing network trimming of the CNN output from the cascade training. The one or more non-transitory computer-readable media store instructions for the at least one processor to perform the steps of training a CNN having three or more layers; performing cascade training on the trained CNN to add one or more intermediate layers until a training error is less than a threshold; and performing network trimming of the CNN output from the cascade training.

At1060, the apparatus (in this instance, a chipset) is tested. Testing1060includes testing whether the apparatus has at least one processor which, when executing instructions stored on one or more non-transitory computer readable media, performs the steps of training a CNN having three or more layers; performing cascade training on the trained CNN to add one or more intermediate layers until a training error is less than a threshold; and performing network trimming of the CNN output from the cascade training; and testing whether the apparatus has the one or more non-transitory computer-readable media which store instructions for the at least one processor to perform the steps of training a CNN having three or more layers; performing cascade training on the trained CNN to add one or more intermediate layers until a training error is less than a threshold; and performing network trimming of the CNN output from the cascade training.

Experimental Validation

A. Cascade Training

In Table A-I, the PSNR/SSIM of a cascade trained CNN in accordance with the present disclosure is compared to non-cascade trained CNN with unsupervised weight initialization from VDSR. It can be seen that with the same network architecture, the PSNR/SSIM of CT-SRCNN is clearly better than non-cascade training.

FIG.11is an exemplary diagram illustrating the convergence speed of cascade trained CNNs according to one embodiment vs. non-cascade trained CNNs. The CT-SRCNN is found to converge faster compared to non-CT-SRCNN. The accuracy of the CT-SRCNN consistently increases when more layers are utilized. This indicates that cascade network training also trains SRCNNs deeper and deeper. Cascade network training performs better compared to conventional training in both accuracy and convergence speed.

In Table A-II, the number of parameters, PSNR, SSIM, and time per image of a CT-SRCNN-13 in accordance with the present disclosure is compared to known SR networks in scale3.

Table A-III shows that the cascade trimmed CT-SRCNN (where 4 out of the 13 layers are trimmed) achieves similar performance to the non-cascade trimmed CT-SRCNN, but the network size is reduced 20%. Cascade network trimming according to the present disclosure is also applied to another network, namely, the fast SR-CNN (FSRCNN) (see Dong et al. 2016b). This network consists of 7 convolutional layers and one deconvolution layer. Similar to trimming the CT-SRCNN according to an embodiment above, 2 layers of the FSRCNN are also trimmed in each stage. Table A-III shows that network cascade trimming according to the present disclosure is also effective for FSRCNN.

There is a trade-off between the trimming rate and the accuracy. If only 2 layers (the 7thand 8th) are trimmed, there is almost no accuracy loss, while 30% of the parameters are removed. If all 8 layers are trimmed (Cascade trimmed FSRCNN 8 layer, trim 8 layers), the accuracy is still better compared to the official model (FSRCNN official lite version), with a smaller network size (3,400 compared to 3,900 parameters).

Table A-IV shows the experimental results of a dilated 13-layer CT-SRCNN. The dilation is applied for the first 9×9 layer, the second 5×5 layer, and the last 5×5 layer. Instead, 5×5, 3×3, and 3×3 2-dilated convolutional layers are utilized. It can be seen that the dilated version of CT-SRCNN can achieve similar PSNR/SSIM to the non-dilated version, but the network size is clearly reduced.

Image enhancement techniques include image and video super resolution for retrieving a high-resolution image(s) from low-resolution input(s), image denoising for generating a clean image from a given noisy input, and compressed image restoration that improves the image quality of a decoded compressed image. Furthermore, different network architectures may be implemented for different image enhancement tasks.

Image compression reduces irrelevance and redundancy of an image in order to store or transmit the image at low bit rates. Image compression is a basic element of image processing that has been used in a typical imaging device. Traditional image coding standards (e.g., JPEG JPEG2000, Better Portable Graphics (BPG)) attempt to distribute the available bits for every nonzero quantized transform coefficient in the whole image. While the compression ratio increases, the bits per pixel (bpp) decreases as a result of the use of bigger quantization steps, which causes the decoded image to have blocking artifacts or noises. To overcome such an issue, the quality of a decoded image may be improved using a post-processing deblocking or denoising method. Typical methods include post-filtering. However, such post-processing methods are very time-consuming because solving the optimal solutions involves computationally expensive iterative processes. Therefore, it is difficult to apply them to practical applications.

Image denoising generates a clean image X from a given noisy image Y which follows an image degradation model Y=D(X). For the additive white Gaussian noise (AWGN) model, the ith observed pixel is yi=D(xi)=xi+ηiwhere ηi˜N(0, σ2) is independent and identically distributed (i.i.d) Gaussian noise with zero mean and variance σ2. AWGN has been used to model the signal-independent thermal noise and other system imperfections. Degradation due to low light shot noise is signal dependent and has often been modeled using Poisson noise, where D(xi)=pi, such that (pi)˜Piis a Poisson random variable with mean xi. However, this noise approaches a Gaussian distribution for average light conditions as P(λ)≈N(λ, λ), for large enough λ. Hence, the noise due to capturing by an imaging device is better modeled as a Poisson noise with AWGN, referred to as Poisson-Gaussian noise, such that yi=D(xi)=αpi+ηifor some scalar α>0.

For image denoising, an input is a noisy image, and an output is a clean image. Additional systems disclosed herein may apply a cascade training network architecture similar to image super resolution as described above but removing the upsample module in the input. The cascade training network architecture may be further applicable to blind denoising, where the noise level is unknown.

System and methods disclosed herein may train a deep CNN for image denoising. The systems and methods may receive a noisy image Y as input and predicts a clean image X as output. Given a training set {(Xi, Yi), i=1, . . . , N} with N samples, the system learns a model S that predicts the clean image {circumflex over (X)}1=S(Yi). The training aims to minimize the mean square error

(MSE)⁢12⁢∑i=1NXi-X^12
over the training set.

For compressed image restoration, an input is a decoded compressed image and an output is a refined image. The systems and methods described below may apply a cascade training network architecture similar to image super resolution as described above but removing the upsample module in the input.

Furthermore, the retrieval from a decoded image to an uncompressed image may be considered as a mapping between two feature maps. The system may apply a neural network to perform retrieval from a decoded image. The system may train a deep CNN from the decoded image to uncompressed ground truth. The CNN receives a decoded image as input and predicts a clean image as output. Given a training set with N samples, the system learns a model that predicts the retrieved image. The training aims to minimize the MSE over the training set.

A residual network (ResNet) has demonstrated considerable performance in computer vision applications such as image classification or super resolution. The systems and methods may provide a denoising residual network (DN-ResNet). DN-ResNet includes residual blocks (ResBlocks) that are gradually inserted into the network stage by stage during training. Such a training strategy allows the resulting DN-ResNet to converge fast and be more computationally efficient than typical denoising networks. In one embodiment, the system modifies the ResBlocks to have a learnable weighted skip connection to provide better denoising performance. The DN-ResNet provides a deep CNN trained for blind denoising of Poisson-Gaussian corrupted images. By cascading a plurality of weighted ResBlocks (e.g., 5), the DN-ResNet achieves state-of-the-art performance on three denoising problems, Gaussian, Poisson, and Poisson-Gaussian, for both known noise-level (non-blind denoising) and unknown noise level (blind denoising). The speed of DN-ResNet is many times faster than previous denoising networks. The DN-ResNet also works well for issues related to compressed image restoration problems. Thus, the DN-ResNet can be generalized to other applications.

FIG.12is an exemplary diagram of a conventional ResBlock1200, according to one embodiment.FIG.13is an exemplary diagram of a simplified ResBlock1300, according to one embodiment.FIG.14is an exemplary diagram of a weighted ResBlock1400, according to one embodiment.

Referring toFIGS.12,13, and14, the DN-ResNet may include a basic element such as a simplified ResBlock1300or a weighted ResBlock1400. Different from the conventional ResBlock1200, the batch normalization (BN) layers1202and1204, and the rectified linear unit (ReLU) layer1206after the addition are removed, since removing such layers does not harm the performance of feature-map based ResNet. In addition, the simplified ResBlock1300may be modified as shown in the weighted ResBlock1400to have learnable weighted skip connections where the skip connection of each ResBlock1400passes through a scale layer1404having several learnable weights {wi, . . . wn}, where n is the number of feature maps at the skip connection.

As the DN-ResNet goes deeper, training and hyper-parameter tuning becomes increasingly difficult. The system may train a deep neural network by cascading the simplified ResBlocks, also referred to as cascade trained ResNet (CT-ResNet). The cascade training separates the entire training into stages and proceeds one by one. The system provides training of CT-ResNet beginning from a simple 3-layer CNN model. The first layer may include 64 9×9 filters, the second layer may include 32 5×5 filters, and the last layer may include 1 5×5 filter. The convolutions may have stride one, and the weights may be randomly initialized from a Gaussian distribution with, for example, a sigma 0.001.

FIG.15is an exemplary diagram of a cascade trained system (a CT-ResNet)1500, according to one embodiment. After the 3-layer CNN is trained, the system1500cascades the ResBlocks stage by stage. In each stage, one new ResBlock is inserted. In the example shown, the training starts from 3 layers, and proceeds to 5 layers, 7 layers, and so on. Each convolutional layer in the ResBlock may include 32 3×3 filters. This ensures a smaller network when going deeper. The new layers are inserted just before the last 5×5 layer. The weights of pre-existing layers are inherited from the previous stage, and the weights of the new ResBlocks are randomly initialized. Hence, only a few weights of CT-ResNet are randomly initialized at each stage, so the convergence is relatively easy. For example, using a fixed learning rate 0.0001 for all layers without any decay is feasible.

Since new convolutional layers reduce the size of the feature map, the system may zero pad 2 pixels in each new 3×3 layer. As a result, all the stages in the cascade training have the same size as the output, so that the training samples could be shared.

FIG.16is an exemplary diagram of color image decoding, according to one embodiment. The system may train CT-ResNets on different color channels, such as red/green/blue (RGB) channels or luma/blue-difference/red-difference (YCbCr) channels, respectively. 2,000 training images are utilized to generate the training data. In the testing, after the compressed image is decoded (e.g., by JPEG 2000 or BPG), the trained CT-ResNets are applied on each channel of the decoded image that is separated into the RGB channels. The retrieved images are further fused to derive a final output. The present system may use the JPEG 2000 (CR=159) and BPG (QF=40) to compress/decode the image, and train CT-ResNets on RGB channels and YCbCr channels.

Further denoising performance improvements are provided by using the edge-aware loss function instead of the conventional mean square error (MSE). By incorporating a depthwise separable ResBlock (DS-ResBlock) into DN-ResNet. DN-DS-ResNet may be fine-tuned from DN-ResNet by the cascade training described above. The ResBlocks in DN-ResNet may be replaced by DN-DS-ResBlocks stage by stage. As a result, further complexity cost reduction is provided, with acceptable accuracy loss.

Although the network size of the cascade trained DN-ResNet is relatively small (e.g., 150K parameters when cascading up to 13 layers), the network size may be further reduced by using depthwise separable DN-ResNet.

FIG.17is a diagram of depthwise separable convolution, according to an embodiment. Referring toFIG.17, the standard convolution layer1702is factored into a depthwise convolution and a 1×1 pointwise convolution1706. The standard convolution layer1702has M input channels and N×K×K filters. In the depthwise convolution layer1704, the M input channels from the standard convolution layer1702are replaced by M depthwise convolution layers with one K×K filter and N 1×1 convolutional layers with M input channels.

As a result, the number of the multiplications are reduced from
M×K×K×N×W×H
to
M×K×K×W×H+M×N×W×H
and the reduction in computation is

FIG.18is a diagram of ResBlocks, according to an embodiment. Referring toFIG.18, a depthwise separable ResBlock1802and a ResBlock1804in DN-ResNet, similar to the ResBlock ofFIG.13, are shown. The standard convolutional layers (Cony) in the ResBlock1804are replaced by depthwise separable convolutional layers (DW-Conv).

In DN-ResNet, the convolutional layers in the ResBlocks have 32 3×3 filters, and the input channel is also32. In ResBlock1804, as an example, the size of the feature map is 640×480, the number of multiplications is therefore
640×480×3×3×32×32×2 or 5.6×109.

In the DS-ResBlock1802, the number of multiplications is
640×480×3×3×32+640×480×32×32 or 9×108.

Thus, the computational cost of the DS-ResBlock1802is reduced 6 times compared to the ResBlock1804.

It is possible to apply the same cascade training procedure to construct the DN-DS-ResNet by cascading DS-ResBlocks. However, since the weights are randomly initialized in cascade training, there will be a large amount of training time. As disclosed herein, another way to train DN-DS-ResNet based on existing DN-ResNet, is referred to as “cascade evolving”.

FIG.19is a diagram of cascade evolving, according to an embodiment. Referring toFIG.19, given a DN-ResNet, to obtain a DN-DS-ResNet, all ResBlocks1902can be replaced by DS-ResBlocks1904, and the entire network can be fine-tuned. Performing this in a single run, the fine-tuning will not converge well. Instead, the ResBlocks1902may be replaced one by one. In each fine-tuning stage (e.g., evolving stage 1, evolving stage 2, evolving stage 3, etc.), only one ResBlock is replaced by a DS-ResBlock, followed by a fine-tuning, as shown inFIG.19.

Similar to cascade training, the weight in the new DS-ResBlock is randomly initialized, and the weight in all other layers is inherited. The replacement starts at the end of the network to make sure there is less influence to the whole network. In each evolving stage, since most of the weights are inherited, the convergence will be relatively easy.

Denoising networks typically aim to minimize the mean square error (MSE)

over the training set. Provided herein is an edge-aware MSE, where pixels in the edges are given higher weights compared to non-edge pixels. An edge-aware loss function may be given as:

where M is an edge map, N is the total number of pixels, and w is a constant. The second term adds a constraint to the loss function. An issue with image denoising occurs because the edges are more difficult to retrieve from a noisy image, especially when the noise level is high. Applying the edge-aware loss function above, the constraint makes the edges less difficult to retrieve. In addition, since high-frequency information such as an edge is more sensitive in human vision, increasing the accuracy of edge pixels with the loss function contributes to the perceptive quality.

In experimental trials, for image denoising, the PASCAL VOC2010dataset is used to generate the training samples. 1,000 testing images are used to evaluate the performance of the DN-ResNet described above, while the remaining images are used for training. Random Gaussian/Poisson/Poisson-Gaussian noisy images are generated with different noise levels. 33×33 noisy patches and corresponding 17×17 clean patches are cropped. Different noise variances σ{circumflex over ( )}2 are considered, where σ∈ {10, 25, 50, 75}. Before corrupting with Poisson or Poisson-Gaussian noise, the input image pixel values are scaled to have max a peak value from the set peak ∈ {1, 2, 4, 8}. For the Poisson-Gaussian noise, σ∈{0.1, 0.2, 0.5, 1, 2, 3, 6, 12} and peak=10σ.

For compressed image restoration, decoded images are obtained. The 33×33 decoded patches and the corresponding 17×17 restored patches are further extracted from the training set. The PSNR is utilized to evaluate the performance. The networks are trained on the Y-channel, but the networks may also be directly applied to the RGB channels without losing much quality.

DN-ResNets from 3 layers to 13 layers are tested on Gaussian, Poisson, and Poisson-Gaussian denoising using PASCAL VOC dataset. These DN-ResNets are trained by cascading the ResBlocks (e.g., ResBlock1300) on known noise level. The MSE loss is utilized for all models. The PSNR consistently increases along with using more layers. From 3 layers to 13 layers, the PSNR is increased 0.4-0.5 dB at all sigmas and peaks. Although the deepest network shown is a 13-layer DN-ResNet, the accuracy may still be further improved by cascading more layers. The cascade training is compared with one-shot training, where a 13-layer DN-ResNet is trained from unsupervised weight initialization. The one-shot training of DN-ResNet-13 results in 0.3 dB lower PSNR than cascade training for all tests. Since the cascade training can be considered as a “partial-supervised initialization,” its convergence will be easier compared to the one-shot training based on unsupervised weight initialization. In Table 2 below, the bold entries represented the best result.

The DN-ResNet was also trained by different edge-aware loss functions described above, as well as utilizing DN-ResNet for blind denoising. As opposed to non-blind denoising where multiple networks are trained for each noise level respectively, only one DN-ResNet is trained for blind denoising by mixing all Gaussian/Poisson/Poisson-Gaussian noises with different noise levels. In Table 3, utilizing DN-ResNet for blind denoising will not decrease much compared to non-blind denoising. This trade-off is valuable since blind denoising does not require a time-consuming noise level estimation. In addition, utilizing an edge-aware loss function can improve the PSNR by 0.1-0.15 dB and enhance the perceptive quality as well. The best group is directly using the gradient magnitude generated from the Sobel operator.

The Sobel operator is used in image processing and computer vision, particularly within edge detection algorithms where it creates an image emphasizing edges. The operator uses two 3×3 kernels which are convolved with the original image to calculate approximations of the derivatives—one for horizontal changes, and one for vertical. If we define A as the source image, and Gx and Gy are two images which at each point contain the horizontal and vertical derivative approximations respectively, the computations are as follows:

The final gradient map G can be obtained by

The DN-ResNet is also constructed by different types of ResBlocks for the blind denoising networks. In Table 4, the DN-DS-ResNet construction by the DS-ResBlocks only decreases the PSNR less than 0.1 dB compared to DN-ResBlock, but the computational cost (e.g., number of multiplications and accumulations (MACs)) and the network size are significantly reduced. This indicates the effectiveness of improving the network efficiency for the above described DS-ResBlock. Further, if the DN-DS-ResNet is constructed by one-shot fine-tuning DN-ResNet, the accuracy will decrease a significant amount. This indicates the effectiveness of the above described cascade evolving. Using the DS-ResBlock together with the edge-aware loss function, high accuracy and less computational cost can be achieved.

The disclosed DN-ResNet and DN-DS-ResNet achieve state-of-the-art performance on Gaussian, Poisson/Poisson-Gaussian denoising, with better efficiency and smaller model size compared to existing deep CNNs. The disclosed networks are effective for both known noise level and unknown noise level.

Besides image denoising, the disclosed DN-ResNet may also be applied on compressed image restoration. For all the compression methods including JPEG, JPEG 2000, and BPG, the DN-ResNet is able to improve the quality of the decoded images. 1-2 dB, 0.5-1.5 dB, and 0.3-0.5 dB gain can be observed for JPEG, JPEG 2000, and BPG respectively.

The DN-ResNet for image denoising achieves both high accuracy and efficiency. Cascade training is efficient and effective in training efficient deep ResNets. The denoising accuracy can be further enhanced by adding learnable weights at the skip connections.

Image SR generates an HR image from a given LR image by attempting to recover the missing information. Recently, deep CNNs have been deployed to solve the image super-resolution problem, as they demonstrate significant accuracy improvements.

Due to the lack of real-world LR-HR patches, images are bicubic-downsampled to create the LR-HR training pairs. This results in clean and noise free LR images. Unfortunately, in the real-world scenario where the images come straight from the camera, there will always be additional noise or unknown degradations. As a result, the state-of-the-art CNN methods trained only to reconstruct images artificially downsampled using bicubic downsampling may lead to dramatic artifacts when applied to real-world images. Disclosed herein is a method to train a real-world SR system, which gives SR outputs with good perceptual quality.

FIG.20illustrates a flowchart of a method for real-world super-resolution, according to an embodiment. At2002, a dataset is generated for real-world SR. The dataset may be generated by downsampling low quality images by a generic degradation model to be used as LR images and directly using the corresponding high quality images to be used as HR images, or by directly using the low quality images as LR images and super resolving the high quality images from the LR images by a generic SR network to be used as HR images.

To train a generic SR network, the SR dataset may be generated based on multiple degradations from the image processing artifacts. More specifically, an LR image x is generated from an HR image y following the degradation model formulated as the follows:
x=D(y*k)+n

where D is the downsampling operation, k is a blur kernel and n is noise. The noise is not necessarily additive.

Multiple downsampling methods are considered, such as nearest neighbor, bilinear, bicubic, and lanczos. The downsampling method is randomly selected when generating the LR patches.

The blur kernel setting of SR is usually simple. The most commonly-used isotropic Gaussian blue kernel parameterized by the standard deviation may be utilized. The standard deviation of the Gaussian kernel may be randomly sampled in the range of [0.2,3] and the kernel size may be fixed at 15×15.

Most of the real-world LR images are noisy due to some image processing artifacts. Some real-world noises include Gaussian, Poisson, or Poisson-Gaussian components. Therefore, Gaussian, Poisson, and Poisson-Gaussian noises may be randomly selected when generating the LR images. The parameters are based on artifacts where the sigma of Gaussian noise ranges from [0, 25], the peak of Poisson noise is uniformly sampled from [50, 150]. When generating the Poisson-Gaussian noise, a similar Poisson peak range may be used, but the Gaussian sigma may be reduced to the range of [0, 5].

In real-world SR, a high-quality dataset of LR-HR images which matches the domain of target imaging device may be crucial for the performance of the SR networks. The present method and system provides an effective methodology for creating a high-quality mobile SR dataset. More specifically, the present method and system creates a new mobile SR dataset using registered image pairs at same scale, without the knowledge of the exact HR-LR degradation model. Given a set of registered mobile-digital single lens reflex (DSLR) images which are relatively easy to obtain in real world, the present method and system super-resolves the DSLR images to create the HR images. The LR images in the dataset are the corresponding mobile patches.

There are multiple ways to super resolve the DLSR images. Applying a simple bicubic upsampling algorithm to super resolve the DSLR patches can generate high-quality HR patches for training the mobile SR networks. The best performance is achieved by using a generic model trained on a generic SR dataset, where the LR-HR pairs are generated by simulating multiple degradations by randomly generating downsampling/blur/noises.

At2004, GANs are trained for generation of the super resolved images, where two GANs are trained by using different discriminators and hyper-parameters in an enhanced super-resolution GAN (ESRGAN). This provides two SR networks with complementary characteristics. The SR network may be trained using a GAN, (e.g., an ESRGAN framework) while using the RCAN as the generator instead of ResNet. During the testing, the generator (RCAN) may be directly used to estimate the HR image from a given LR image. To further improve the perceptual quality, two RCANs may be trained by using different discriminators and hyperparameters in ESRGAN. This provides two SR networks with complementary characteristics. The final SR prediction is a pixel-wise ensemble of these two RCANs.

FIG.21illustrates a diagram of an SR approach, according to an embodiment. An SR dataset2102is input into two GANs for training. The first GAN2104processes LR images in an RCAN generator2108, which produces fake estimated HR images2110. Real HR images2112and the estimated HR images2110are processed by a standard discriminator2114, and the GAN2104produces a real/fake decision2116. The second GAN2120processes LR images2122in an RCAN generator2124, which produces fake estimated HR images2126.

Real HR images2128and the estimated HR images2126are processed by a relativistic discriminator2130, and the GAN2120produces a real/fake decision2132.

FIG.22illustrates a diagram of an RCAN2200, according to an embodiment. The RCAN may be based on a residual in residual (RIR) structure, which includes several residual groups with long skip connections. Each residual group includes some residual blocks (ResBlock) with short skip connections. In each ResBlock, the channel attention mechanism may be utilized to adaptively rescale channel-wise features by considering interdependencies among channels. The RCAN may be used as the generator when training the GAN networks.

The generative adversarial networks (e.g., super resolution GAN (SRGAN)) takes advantage of the strength of the GAN networks to model the space of natural images, and uses perceptual and adversarial losses to guide the SR network to favor output images that reside on the manifold of natural images. After that, several modifications related to the perceptual-driven, GAN-based approach in SRGAN are provided. In one embodiment, the ESRGAN framework is utilized due to the usage of relativistic discriminator, which is able to create sharper edges and more realistic texture details.

The present method and system may train two GANs based on the generated real-world SR dataset. Both of their generators may be implemented by RCAN. Two discriminators are utilized: (1) a relativistic discriminator, which predicts if the real image (HR image) is more realistic than the fake image (HR output); and (2) a standard discriminator, which simply predicts the real image is real, and the fake image is fake.

During the training, the loss function of the relativistic generator LGRincludes the L1image loss, perceptual loss Lpand the adversarial loss LαR, as follows:
LGR=Lp+λLαR+ηL1

L1=xi[∥G(xi)−γi∥] calculates the L1distance between the super-resolved image G(xi) from the RCAN generator network G(·) and the ground truth HR image γi.xi[∥·∥] denotes the operation of taking an average for all the images in the mini-batch. The perceptual loss Lpcalculates a feature maps distance of xiand yiwith the usage of a pre-trained 19 layer VGG network. The adversarial loss LαRis based on the relativistic GAN discriminator and is defined as follows:
LαR=−xr[log(1−DR(xr,xf)]−xf[log(DR(xf,xr)]

in which DR(xr, xf)=σ(C(xr)−xf[C(xf)]) represents the discriminator network based on the relativistic GAN, C(xf) is the non-transformed discriminator output, and σ is the sigmoid function. DR(xr, xf) predicts if the real image is more realistic than the fake image rather than deciding if an input image is absolutely real or fake. If the output image from the generator network G is denoted as xfand the corresponding real image is denoted as xr, then the corresponding discriminator loss function may be defined as follows:
LDR=−xr[log(DR(xr,xf)]−xf[log(1−DR(xf,xr)]

The hyperparameters λ and η determine the contribution of different loss components in the final loss function. The parameter η may be increased to reduce the quantitative error of the estimate while increasing the adversarial loss weight may result in an improvement in perceptual quality of the results. Another GAN can be trained with an RCAN generator but with a different generator loss LDSbased on the standard GAN, and may be shown as follows:
LGS=Lp+ΔLαS+ηL1

where LGSis the adversarial loss based on the standard GAN whose corresponding discriminator can be described as DS(x)=σ(C(x)).

At2006, the outputs from the trained GANs are fused. The outputs may be fused according to an illumination threshold. If the illumination level of the second GAN's output (the GAN using the relativistic discriminator) is lower than the illumination threshold, then the pixel-wise ensemble of the two outputs from the GANS are used as the final output. If the illumination level of the second GAN's output is higher than the illumination threshold, the output of the second GAN is used as the final output.

The final SR prediction may be a pixel-wise ensemble of these two RCANs. The SR outputs generated by the GAN with relativistic generator show good perceptual quality in high-frequency regions. In contrast, the SR outputs generated by the GAN with standard generator generate less artifacts in the smooth regions of some low-illumination images. The fusion uses the SR estimates of the above two GANs. A selective average technique may be utilized based on the median brightness of all pixels in the image in order to improve the visual quality on low-illumination images. The HR output from the GAN trained using relativistic GAN loss as γSRR, and the HR output from the RCAN model trained using standard GAN loss function is γSRS. The fused output image γSRfusedis derived as:

where Ymedianis the median of the pixel intensity values of all the pixels in the Y (luminance) component of the YCbCr color space representation, α is a constant between [0, 1]. Two GAN models trained based on different adversarial losses with different complementary effects are fused, so that the perceptual quality of the two images used for fusion are close. This ensures that artifacts may be reduced in some regions of low-illumination images while not sacrificing the overall perceptual quality.

FIG.23illustrates a diagram of fusing results from trained GANs, according to an embodiment. An image2302is input into a first GAN2304that produces a smooth output with less artifacts (e.g., using a standard decimator) and into a second GAN2306that produces a sharp output with artifacts (e.g., using a relativistic decimator). The two outputs are fused into a pixel-wise ensemble2308to produce an output2310.

The steps and/or operations described above in relation to an embodiment of the present disclosure may occur in a different order, or in parallel, or concurrently for different epochs, etc., depending on the specific embodiment and/or implementation, as would be understood by one of ordinary skill in the art. Different embodiments may perform actions in a different order or by different ways or means. As would be understood by one of ordinary skill in the art, some drawings are simplified representations of the actions performed, their descriptions herein simplified overviews, and real-world implementations would be much more complex, require more stages and/or components, and would also vary depending on the requirements of the particular implementation. Being simplified representations, these drawings do not show other required steps as these may be known and understood by one of ordinary skill in the art and may not be pertinent and/or helpful to the present description.

Similarly, some drawings are simplified block diagrams showing only pertinent components, and some of these components merely represent a function and/or operation well-known in the field, rather than an actual piece of hardware, as would be understood by one of ordinary skill in the art. In such cases, some or all of the components/modules may be implemented or provided in a variety and/or combinations of manners, such as at least partially in firmware and/or hardware, including, but not limited to one or more application-specific integrated circuits (“ASICs”), standard integrated circuits, controllers executing appropriate instructions, and including microcontrollers and/or embedded controllers, field-programmable gate arrays (“FPGAs”), complex programmable logic devices (“CPLDs”), and the like. Some or all of the system components and/or data structures may also be stored as contents (e.g., as executable or other machine-readable software instructions or structured data) on a non-transitory computer-readable medium (e.g., as a hard disk; a memory; a computer network or cellular wireless network or other data transmission medium; or a portable media article to be read by an appropriate drive or via an appropriate connection, such as a DVD or flash memory device) so as to enable or configure the computer-readable medium and/or one or more associated computing systems or devices to execute or otherwise use or provide the contents to perform at least some of the described techniques.

One or more processors, simple microcontrollers, controllers, and the like, whether alone or in a multi-processing arrangement, may be employed to execute sequences of instructions stored on non-transitory computer-readable media to implement embodiments of the present disclosure. In some embodiments, hard-wired circuitry may be used in place of or in combination with software instructions. Thus, embodiments of the present disclosure are not limited to any specific combination of hardware circuitry, firmware, and/or software.

The term “computer-readable medium” as used herein refers to any medium that stores instructions which may be provided to a processor for execution. Such a medium may take many forms, including but not limited to, non-volatile and volatile media. Common forms of non-transitory computer-readable media include, for example, a floppy disk, a flexible disk, hard disk, magnetic tape, or any other magnetic medium, a CD-ROM, any other optical medium, punch cards, paper tape, any other physical medium with patterns of holes, a RAM, a PROM, and EPROM, a FLASH-EPROM, any other memory chip or cartridge, or any other medium on which instructions which can be executed by a processor are stored.

Some embodiments of the present disclosure may be implemented, at least in part, on a portable device. “Portable device” and/or “mobile device” as used herein refers to any portable or movable electronic device having the capability of receiving wireless signals, including, but not limited to, multimedia players, communication devices, computing devices, navigating devices, etc. Thus, mobile devices include (but are not limited to) user equipment (UE), laptops, tablet computers, portable digital assistants (PDAs), mp3 players, handheld PCs, instant messaging devices (IMD), cellular telephones, global navigational satellite system (GNSS) receivers, watches, or any such device which can be worn and/or carried on one's person.

Various embodiments of the present disclosure may be implemented in an integrated circuit (IC), also called a microchip, silicon chip, computer chip, or just “a chip,” as would be understood by one of ordinary skill in the art, in view of the present disclosure. Such an IC may be, for example, a broadband and/or baseband modem chip.

While several embodiments have been described, it will be understood that various modifications can be made without departing from the scope of the present disclosure. Thus, it will be apparent to those of ordinary skill in the art that the present disclosure is not limited to any of the embodiments described herein, but rather has a coverage defined only by the appended claims and their equivalents.