# Patent

## Title
Ensemble sparse models for image analysis and restoration

## Abstract

Methods and systems for recovering corrupted/degraded images using approximations obtained from an ensemble of multiple sparse models are disclosed. Sparse models may represent images parsimoniously using elementary patterns from a “dictionary” matrix. Various embodiments of the present disclosure involve simple and computationally efficient dictionary design approach along with low-complexity reconstruction procedure that may use a parallel-friendly table-lookup process. Multiple dictionaries in an ensemble model may be inferred sequentially using greedy forward-selection approach and can incorporate bagging/boosting strategies, taking into account application-specific degradation. Recovery performance obtained using the proposed approaches with image super resolution and compressive recovery can be comparable to or better than existing sparse modeling based approaches, at reduced computational complexity. By including ensemble models in hierarchical multilevel learning, where multiple dictionaries are inferred in each level, further performance improvements can be obtained in image recovery, without significant increase in computational complexity.

## Background
CROSS REFERENCE TO RELATED APPLICATION
The present non-provisional utility application claims priority under 35 U.S.C. § 119(e) to provisional application No. 61/788,124 entitled “Ensemble Sparse Models For Image Analysis And Restoration,” filed on Mar. 15, 2013 and which is hereby incorporated by reference in its entirety.


FIELD
Aspects of the present disclosure relate to systems and method for the analysis and restoration of degraded images.
BACKGROUND
The elementary set of patterns available in natural images is finite, and can typically consist of edges, wedges, curves, smooth regions and textures. These elementary patterns can be stacked together as columns of a “dictionary” matrix. Any local region commonly referred to as a “patch” in a natural image can be represented as a linear combination of a small number of columns from the dictionary. This operation is referred to as sparse coding, and computing sparse codes for images is a nonlinear and complex procedure. The dictionaries for sparse coding may be pre-defined but they can also be adapted to a set of training examples. Dictionary learning is an iterative procedure that is often times is computationally expensive and usually requires a large training set is necessary to obtain good dictionaries.
Sparse approximation using learned dictionaries have been shown to provide improved performance for restoring degraded data in applications such as image denoising, inpainting, deblurring, single image super resolution, and compressive sensing.



BRIEF DESCRIPTION OF THE DRAWINGS
The foregoing and other objects, features, and advantages of the present disclosure set forth herein will be apparent from the following description of particular embodiments of those inventive concepts, as illustrated in the accompanying drawings. It should be noted that the drawings are not necessarily to scale; however, the emphasis instead is being placed on illustrating the principles of the inventive concepts. Also, in the drawings the like reference characters refer to the same parts throughout the different views. The drawings depict only typical embodiments of the present disclosure and, therefore, are not to be considered limiting in scope.
FIGS. 1-2 are illustrations of ensemble model performance tables, according to aspects of the present disclosure;
FIG. 3 is an illustration of a convergence table, according to aspects of the present disclosure;
FIG. 4 a flow chart illustrating a boosted dictionary learning process, according to aspects of the present disclosure;
FIG. 5 is an illustration of images, according to aspects of the present disclosure;
FIG. 6 is an illustration of a table, according to aspects of the present disclosure;
FIG. 7 is an example images, according to aspects of the present disclosure;
FIG. 8 illustrates an example computing device, according to aspects of the present disclosure;
FIG. 9 illustrates an example computing system, according to aspects of the present disclosure; and
FIG. 10 shows one table for compressed recovery of standard images and another table for superresolution of standard images.



DETAILED DESCRIPTION
Aspects of the present disclosure involve methods and systems for automatically learning an ensemble of sparse models that may be executed to perform image recovery. In various aspects, approximations of an image may be obtained from the multiple sparse models included within the ensemble. The approximations may be based, in part, on learned dictionaries. Two different methodologies for constructing the dictionaries may be implemented. First, a dictionary may be generated using different random subsets from the data as the dictionary (referred to herein as “bagging”). Alternatively and/or in addition, a dictionary may be generated sequentially such that the current dictionary prioritizes the poorly represented samples in the previous iteration (referred to herein as “boosting”). In some embodiments, boosted approaches can learn one dictionary per round such that the dictionary learned in a particular round is optimized for the training examples having high reconstruction error in the previous round. The approximations obtained from the individual sparse models may be suitably aggregated in order to obtain a final reconstruction of the recovered image.
As will be described in further detail below, an ensemble model approximates the data of an image as a linear combination of approximations obtained from multiple weak sparse models. Images reveal statistics that allow them to be efficiently represented using a sparse linear combination of elementary patterns. The local regions of natural images, referred to as patches, can be represented using a sparse linear combination of columns from a dictionary matrix. Given a data sample XεM, and a dictionary matrix DεMxK, the data approximated using the linear generative model as Da, where aεM is the sparse coefficient vector. This generative model that incorporates sparsity constraints in the coefficient vector, will be referred to as the “sparse model”. The dictionary can be either pre-defined or learned from the training examples themselves. Learned dictionaries have been shown to provide improved performance for restoring degraded data in applications such as denoising, inpainting, deblurring, superresolution, and compressive sensing and also in machine learning applications such as classification and clustering.
Using the proposed ensemble models may provide comparable or better performance at a much reduced computational cost. Moreover, by including ensemble models in hierarchical multilevel learning, where multiple dictionaries are inferred in each level, further performance improvements are obtained, without significant increase in computational complexity.
Sparse Coding and Dictionary Learning
Using the linear generative model, the sparse code of a data sample x can be obtained by optimizing,




  
    
      
        
          h
          ⁡
          
            (
            
              x
              ,
              D
            
            )
          
        
        =
        
          
            
              min
              a
            
            ⁢
            
              
                
                
                  x
                  -
                  Da
                
                
              
              2
              2
            
          
          +
          
            λ
            ⁢
            
              
                
                a
                
              
              1
            
          
        
      
    
    
      
        (
        1
        )
      
    
  




Here ∥a∥ is the l1 penalty that promotes sparsity of the coefficients, and the equivalence of (1) to l0 minimization has been analyzed under some strong conditions on the dictionary. Some methods to obtain sparse representations used include the Matching Pursuit (MP), Orthogonal Matching Pursuit (OMP), Order-Recursive Matching Pursuit, Basis Pursuit (BP), FOCUSS and iterated shrinkage algorithms.
In several image processing and sample machine learning applications, it may be appropriate to learn a dictionary, such that the set of training examples obtained from a probability space have a small approximation error with sparse coding. This problem can be expressed as minimizing the objective

g(D)=EX[h(x,D)]  (2)

Where the columns of D, referred to as dictionary atoms, are constrained to have unit l2 norm, i.e., ∥dj∥2≦1, ∀j. If the distribution in the probability space is unknown and we only have T training examples {xi}i=1T, each with probability mass p(xl), (2) can be modified as the empirical cost function




  
    
      
        
          g
          ⁡
          
            (
            D
            )
          
        
        =
        
          
            ∑
            
              i
              =
              1
            
            T
          
          ⁢
          
            
              h
              ⁡
              
                (
                
                  
                    x
                    1
                  
                  ,
                  D
                
                )
              
            
            ⁢
            
              p
              ⁡
              
                (
                
                  x
                  i
                
                )
              
            
          
        
      
    
    
      
        (
        3
        )
      
    
  




Typically dictionary learning algorithms solve for the sparse codes using (1), and obtain the dictionary by minimizing ĝ(D), repeating the steps until convergence. We refer to this baseline algorithm as “Alt-Opt”. Since this is an alternating minimization process, it may be prudent to provide a good initial dictionary and this is performed by setting the atoms to normalized cluster centers of the data. Since the number of examples T is usually much larger than the number of dictionary atoms K, it is much more computationally intensive to obtain sparse representations with examples. Nevertheless, both learned and example-based dictionaries have found applications in inverse problems and also in machine learning applications such as clustering and classification.
Ensemble Sparse Models
According to one embodiment, data, such as images, image data, and/or the like, may be represented using a linear combination of L different sparse approximations, instead of being represented using an approximation obtained from a single sparse model. The approximation to x can be obtained by optimizing




  
    
      
        
          min
          
            
              {
              
                β
                l
              
              }
            
            
              l
              =
              1
            
            L
          
        
        ⁢
        
          
            
            
              x
              -
              
                
                  ∑
                  
                    l
                    =
                    1
                  
                  L
                
                ⁢
                
                  
                    β
                    l
                  
                  ⁢
                  
                    D
                    l
                  
                  ⁢
                  
                    a
                    l
                  
                
              
            
            
          
          2
          2
        
      
    
    
      
        (
        4
        )
      
    
  




Here each coefficient vector al is assumed to be sparse, and is obtained by solving for the optimization (1) with Dl as the dictionary. The weights {βl}l=1L control the contribution of each base model to the ensemble.
Since the ensemble combines the contributions of multiple models, it is sufficient that the dictionary for model is obtained using a “weak” training procedure. Accordingly, “weak dictionaries,” {Dl}l=1L, may be learned sequentially, using a greedy forward selection procedure, such that training examples that incurred a high approximation error with the dictionary Dl are given more importance while learning Dl+1. Furthermore, in various embodiments, an ensemble model where each individual dictionary is designed as a random subset of training samples, may be learned. In supervised learning, boosting is used to improve the accuracy of learning algorithms, using multiple weak hypotheses instead of a single strong hypothesis. The proposed ensemble sparse models are geared towards image restoration.
Two approaches may be used for learning the ensemble: (a) using a random selection and averaging (RandExAv) approach, where each dictionary is chosen as a random subset of the training examples, and (b) using a boosted approach to learn dictionaries sequentially by modifying the probability masses of the training examples in each round of learning. In the boosted approach, a method to learn the weak dictionaries by example selection using the probability distribution on the training set (BoostEx) is provided. For all cases of ensemble learning, we also provide methods to obtain the ensemble weights {βl}l=1L, from the training examples.
According to one embodiment, the learned ensemble models may be applied during image restoration. The image restoration problems that we consider here are compressive sensing using random projections and single image superresolution. When boosted ensemble models are learned for image restoration problems, the form of degradation operator specific to the application is also considered, thereby optimizing the ensemble for the application. For compressive recovery, we compare the performance of the proposed ensemble models to the single sparse model, whose dictionary is obtained using the Alt-Opt approach. Note that, the base sparse model for the example-based approaches is designed as a random subset of examples, and hence it requires minimal training.
Need for the Ensemble Model
In several scenarios, a single sparse model may be insufficient for representing the data, and using an ensemble model instead may result in a good performance. For example, statistical issues may arise, whereby several sparse models may have a similar training error when learned from a limited number of training samples. However, the performance of each of these models individually with test data can be poor. By averaging representations obtained from an ensemble, we may obtain an approximation closer to the true test data. The second reason is computational, which can occur with the case of large training sets also. The inherent issue in this case is that sparse modeling is a problem with locally optimal solution. Therefore, we may never be able to reach the global optimal solution with a single model and hence using an ensemble model may result in a lesser representation error. Note that this case is quite common in dictionary learning, since many dictionary learning algorithms only seek a local optimal solution. The third reason for using an ensemble model is representational, wherein the hypothesis space assumed cannot represent the test data sample. In the case of sparse models, this corresponds to the case where the dictionary may not provide a high-fidelity sparse approximation for a novel test data sample. This also happens in the case where the test observation is a corrupted version of the underlying test data, and there is ambiguity in obtaining a sparse representation. In this case also, it may be necessary to combine multiple sparse models to improve the estimate of the test data.
In order to simplify notation in the following analysis, let us denote the ith approximation in the ensemble model as cl=Dlal. CεRMxL, where C=[c1 . . . cL] and the weight vector is denoted as β=[β1 . . . βL]T. The individual residuals are denoted as rl=x−cl, for i=1, . . . , L, and the total residual of the approximation is given as

r=x−Cβ  (5)

We characterize the behavior of the ensemble sparse model by considering four different cases for the weights {βl}l=1L.
Unconstrained Weights:
In this case, the ensemble weights {βl}l=1L are assumed to be unconstrained and computed using the unconstrained least squares approximation




  
    
      
        
          min
          β
        
        ⁢
        
          
            
            
              x
              -
              
                C
                ⁢
                
                  
                
                ⁢
                β
              
            
            
          
          2
          2
        
      
    
    
      
        (
        6
        )
      
    
  




When the data x lies in the span of C, the residual will be zero, i.e., r=0. The residual that has minimum energy in the L approximations is denoted as rmin. This residual can be obtained by setting the corresponding weight in the vector β to be 1, whereas (6) computes β that achieves the best possible residual r for the total approximation. Clearly this implies

∥r∥2≦∥rmin∥2  (7)

Therefore, at worst, the total approximation will be as good as the best individual approximation. βl≧0: The ensemble weights {βl}l=1L are assumed to be non-negative in this case. The least squares approximation (4), with the constraint β≧0 will now result in a zero residual if the data x lies in the simplical cone generated by the columns of C. The simplical cone is defined as the set




  
    {
    
      
        b
        ⁢
        
          :
        
        ⁢
        
          
        
        ⁢
        b
      
      =
      
        
          ∑
          
            l
            =
            1
          
          L
        
        ⁢
        
          
            c
            l
          
          ⁢
          
            β
            l
          
        
      
    
    }
  
  .




Otherwise, the bound on the total residual given by (7) holds in this case, since rmin can be obtained by setting the appropriate weight in β to 1 in (5), and the rest to 0 in this case also





  
    
      ∑
      
        l
        =
        1
      
      L
    
    ⁢
    
      β
      l
    
  
  =
  
    1
    ⁢
    
      :
    
  




When the ensemble weights are constrained to sum to 1, the total residual can be expressed as





  
    
      
        r
        =
        
          
            ∑
            
              l
              =
              1
            
            L
          
          ⁢
          
            
              β
              l
            
            ⁢
            
              r
              l
            
          
        
      
    
    
      
        (
        8
        )
      
    
  




This can be easily obtained by replacing x as




  
    ∑
    
      l
      =
      1
    
    L
  
  ⁢
  
    
      β
      l
    
    ⁢
    x
  




in (5). Denoting the residual matrix R=[r1 . . . rL], the optimization (4) to compute the weights can also be posed as minβ∥Rβ∥2. Incorporating the constraint





  
    
      
        ∑
        
          l
          =
          1
        
        L
      
      ⁢
      
        β
        l
      
    
    =
    1
  
  ,




it can be seen that the final approximation Cβ lies in the affine hull generated by the columns of C, and the final residual, Rβ, lies in the affine hull generated by the columns of R. Clearly the final residual will be zero, only if the data x lies in the affine hull of C, or equivalently the zero vector lies in the affine hull of R. When this does not hold, the worst case bound on r given by (7) holds in this case as well.





  
    
      β
      l
    
    ≥
    0
  
  ,
  
    
      
        ∑
        
          l
          =
          1
        
        L
      
      ⁢
      
        β
        l
      
    
    =
    
      1
      :
    
  




Similar to the previous case, the total residual can be expressed as (8). As a result, the final representation Cβ lies in the convex hull generated by the columns of C, and the final residual, Rβ, lies in the convex hull generated by the columns of R. Furthermore, the final residual will be zero only if the zero vector lies in the convex hull of R. Clearly, the worst case bound on r given by (7) holds in this case.

The first case is unconstrained and it should result in the least error. The second case constrains that the solution should lie in the simplical cone spanned by the columns of C, and this should lead to higher residual energy than Case 1. Case 3 constrains the solution to lie in an affine hull, which is of L−1 dimensions compared to simplical cone in L dimensions, so it could lead to a higher error compared to Case 2. Case 4 is the subset of constraint spaces for Cases 1 to 3 and hence it will lead to the highest residual error.
Demonstration of Ensemble Representations
In order to demonstrate the performance of ensemble representations with real data, we obtain a random set of 100,000 patches, each of size 8×8, from a set of natural images. The training images consisted of a wide variety of patterns and textures. The chosen patches are then processed to remove the mean, followed by the removal of low-variance patches. Since image recovery is the important application of the proposed models, considering high-variance patches alone is beneficial. Each dictionary in the ensemble D, is obtained as a random set of K vectorized, and normalized patches. We fix the number of models in the ensemble as L=20. The test data is a random set of 1000 grayscale patches different from the training set. For each test sample, we compute the set of L approximations using the sparse model given in (1), with λ=0.2. The individual approximations are combined into an ensemble, under the four conditions on the weights, {βι}, described above. The optimal weights are computed and the mean squared norm of the residuals for all the test samples are compared in the table 100 of FIG. 1, for the dictionary sizes K={256, 1025, 2048}. compared in the table 100 of FIG. 1, for the dictionary sizes K={256, 1025, 2048}. Referring now to the table 100, we observe that the performance of the ensembles generally improve as the size of the dictionaries 102 used in the base models increase. The variation in performance across all the four cases of weights follows our reasoning in the previous section. We refer to these as “oracle” ensemble models, since the weights are optimally computed with perfect knowledge of all the individual approximations and the actual data. In practice, the weights will be precomputed from the training data.
Proposed Ensemble Sparse Representation Algorithms
The ensemble model proposed in (4) results is a good approximation for any known data. However, in order to use ensemble models in analysis and recovery of images, that are possibly corrupted or degraded, both the weights {βl}l=1L and the dictionaries {Dl}l=1L must be inferred from uncorrupted training data. Let us denote the set of training samples as X=[x1 x2 . . . xT], and the set of coefficients in base model l as Al=[al,1 al,2 . . . , al,T], where al,i is the coefficient vector of the ith sample for base model l. In the proposed ensemble learning procedures, we consider both simple averaging and boosting approaches.
Random Example Averaging Approach
The first approach chooses L random subsets of K samples from the training data itself and normalizes them to form the dictionaries {Dl}l=1L. The weights {βl} are chosen to be equal for all base models as 1/L. Note that the selection of dictionaries follows the same procedure as given in the previous demonstration described above. We refer to this ensemble approach as Random Example Averaging (RandExAv).
Boosting Approaches
The next two approaches use boosting and obtain the dictionaries and weights sequentially, such that the training examples that resulted in a poor performance with the l−1th base model are given more importance when learning the lth base model. We use a greedy forward selection procedure for obtaining the dictionaries and the weights. In each round l, the model is augmented with one dictionary Dl, and the weight αl corresponding to the dictionary is obtained. The cumulative representation for round l is given by

Xl=(1−+αl)Xl−1+αlDlAl  (9)

Note that the weights of the greedy forward selection algorithm, αl, and the weights of the ensemble model, βl, are related as




  
    
      
        
          β
          l
        
        =
        
          
            α
            l
          
          ⁢
          
            
              ∏
              
                t
                =
                
                  l
                  +
                  1
                
              
              L
            
            ⁢
            
              
            
            ⁢
            
              (
              
                1
                -
                
                  α
                  t
                
              
              )
            
          
        
      
    
    
      
        (
        10
        )
      
    
  




From (9), it can be seen that Xl lies in the affine hull of Xl−1 and DlAl. Furthermore, from the relationship between the weights {αl} and {βl} given in (10), it is clear that




  
    
      ∑
      
        l
        =
        1
      
      L
    
    ⁢
    
      β
      l
    
  
  =
  1




and hence the ensemble model uses the constraints given in Case 3. Only the Cases 3 and 4 lead to an efficient greedy forward selection approach for the ensemble model in (4), and we use Case 3 since it leads to a better approximation performance (FIG. 1).

In boosted ensemble learning, the importance of the training samples in a particular round is controlled by modifying their probability masses. Each round consists of (a) learning a dictionary Dl corresponding to the round, (b) computing the approximation for the current round l, (c) estimating the weight αl, (d) computing the residual energy for the training samples, and (e) updating the probability masses of the training samples for the next round. Since the goal of ensemble approaches is to have only weak individual models, Dl is obtained using naive dictionary learning procedures as described later in this section. The dictionaries for the first round are obtained by fixing uniform probability masses for each training example in the first round, (i.e.), p1(xl)=1/T for i={1, 2, . . . , T}. Assuming that Dl is known, the approximation for the current round is computed by coding the training samples X with the dictionary using (1). The weight αι is computed such that the error between the training samples and the cumulative approximation Xl is minimized. Using (9), this optimization can be expressed a




  
    
      
        
          min
          
            α
            l
          
        
        ⁢
        
          
            
            
              
                X
                l
              
              -
              
                [
                
                  
                    
                      (
                      
                        1
                        -
                        
                          α
                          
                            l
                            ⁢
                            
                              
                            
                          
                        
                      
                      )
                    
                    ⁢
                    
                      X
                      
                        l
                        -
                        1
                      
                    
                  
                  +
                  
                    
                      α
                      l
                    
                    ⁢
                    
                      D
                      l
                    
                    ⁢
                    
                      A
                      l
                    
                  
                
                ]
              
            
            
          
          F
          2
        
      
    
    
      
        (
        11
        )
      
    
  




and can be solved in closed form with the optimal value given as,




  
    
      
        
          α
          l
        
        =
        
          
            
              T
              r
            
            ⁡
            
              [
              
                
                  
                    (
                    
                      X
                      -
                      
                        X
                        
                          l
                          -
                          1
                        
                      
                    
                    )
                  
                  T
                
                ⁢
                
                  (
                  
                    
                      
                        D
                        l
                      
                      ⁢
                      
                        A
                        l
                      
                    
                    -
                    
                      X
                      
                        l
                        -
                        1
                      
                    
                  
                  )
                
              
              ]
            
          
          
            
              
                
                
                  
                    
                      D
                      l
                    
                    ⁢
                    
                      A
                      l
                    
                  
                  -
                  
                    X
                    
                      l
                      -
                      1
                    
                  
                
                
              
              F
              2
            
            ⁢
            
              
            
          
        
      
    
    
      
        (
        12
        )
      
    
  




wherein Tr denotes the trace of the matrix. The residual matrix for all the training samples in round l is given by Rl=X−DlAl. The energy of the residual for the ith training sample is given as el(i)=∥rl,i∥22. If the dictionary in round l provides a large approximation error for sample i, then that sample will be given more importance in round l+1. This will ensure that the residual error for sample i in round l+1 will be small. The simple scheme of updating the probability masses as pl+1(xi)=el(i), upweights the badly represented samples and downweights the well-represented ones for the next round.
Given a training set {xi}i=1L, and its probability masses {pl(xi)}i=1L, we will propose two simple approaches for learning the dictionaries corresponding to the individual sparse models.
BoostEx: From (3), it is clear that the learned dictionary atoms are close to training samples that have higher probabilities. Therefore, in the BoostEx method, the dictionary for round l is updated by choosing K data samples based on the non-uniform weight distribution, {pl(xi)}i=1T, and normalizing them. This scheme will ensure that those samples with high approximation errors in the previous round, will be better represented in the current round.
Demonstration of the Proposed Approaches
The performance of the proposed ensemble schemes for dictionaries of three different sizes K={256, 1024, 2048} are compared. The same training set is used with the RandExAv and BoostEx schemes. The dictionaries {Dι}ι=1L and the weights {βι}ι=1L are obtained with the above schemes for L=20. The individual approximations in the training set are obtained using (1) with the sparsity penalty set as λ=0.2. For each sample in the test set, the individual representations are computed using (1) with λ=0.2. The ensemble approximation the ith test sample is obtained as Σι=1LβιDιaι,i. FIG. 2 provides an illustration of a table 200 that compares the performances of the proposed schemes 202, 204, and 206 for different dictionary sizes 208, 210, and 212. The minimum error obtained across all individual approximations 214, 216, and 218 is also shown for comparison, with both methods and the different dictionary sizes (i.e. dictionaries 208, 210, and 212). It can be seen that the proposed schemes satisfy the basic property of the ensemble, i.e., the ensemble approximation performs better than the best constituent individual approximation of table 100 if FIG. 1. As the number of number of approximations in the ensemble increase, the average mean squared error (MSE) 220 for the proposed methods reduce, as shown in table 300 of FIG. 3 for a dictionary size of 1024. Clearly, increasing the number of models in the ensemble results in a better approximation, but the MSE flattens out as the number of rounds increase.
Application: Image Restoration
In restoration applications, it is necessary to solve an inverse problem, in order to estimate the test data y from

z=Φ(y)+n  (14)

where Φ(.) is the corruption operator and n is the additive noise. If the operator Φ(.) is linear, we can represent it using the matrix Φ. With the prior knowledge that y is sparsely representable in a dictionary D according to (1), (14) can be expressed as z=ΦDa+n. Restoring x now reduces to computing the sparse codes a by solving




  
    
      
        
          
            min
            a
          
          ⁢
          
            
              
              
                z
                -
                
                  Φ
                  ⁢
                  
                    
                  
                  ⁢
                  Da
                
              
              
            
            2
            2
          
        
        +
        
          λ
          ⁢
          
            
              
              a
              
            
            1
          
        
      
    
    
      
        (
        15
        )
      
    
  




and finally estimating y=Da. In the proposed ensemble methods, the final estimate of x is obtained as a weighted average of the individual approximations. Furthermore, in the boosting approach, the degradation operation can be included when learning the ensemble. This is achieved by degrading the training data as ΦX, and obtaining the approximation with the coefficients computed using (15) instead of (1). The procedure to obtain boosted dictionaries using degraded data and computing the final approximation is illustrated in FIG. 4. As illustrated, multiple rounds of boosting 402, 404, and 406 are calculated or otherwise executed from training data 400 to generate the final approximation 410. Stated differently, the final approximation 410 is estimated sequentially using the weights {αι}ι=1L, but it is equivalent to computing {βι}ι=1L using (10) and computing the ensemble estimate Σι=1LβιDιAι.
Compressive Recovery
In compressed sensing (CS), the N-dimensional observation z is obtained by projecting the M-dimensional data y onto a random linear subspace, where N<<M. In this case, the entries of the degradation matrix φεNxM are obtained as i.i.d. realizations of a Gaussian or Bernoulli random variable. Compressive recovery can be effectively performed using conventional dictionaries or ensemble dictionaries. In addition, the proposed idea of ensemble learning can be incorporated in existing learning schemes to achieve improved recovery performance. In particular, the multilevel dictionary learning algorithm can be very easily adapted to compute ensemble representations. Before discussing the experimental setup, and the results of the proposed methods, we will describe the modification to multilevel dictionary learning for improving the compressed recovery performance with learned dictionaries.
Improved Multilevel Dictionaries:
The multilevel dictionary (MLD) learning algorithm is a hierarchical procedure where the dictionary atoms in each level are obtained using a one dimensional subspace clustering procedure. Multilevel dictionaries have been shown to generalize well to novel test data, and have resulted in high performance in compressive recovery. We propose to employ the RandExAv procedure in each level of multilevel learning to reduce overfitting and thereby improve the accuracy of the dictionaries in representing novel test samples. In each level, L different dictionaries are drawn as random subsets of normalized training samples. For each training sample, a 1-sparse representation is computed with each individual dictionary, and the approximations are averaged to obtain the ensemble representation for that level. Using the residual vectors as the training data, this process is repeated for multiple levels. The sparse approximation for a test sample is computed in a similar fashion. Since the sparse code computation in each level is performed using simple correlation operations, the computation complexity is not increased significantly by employing ensemble learning. In our simulations, we will refer to this approach as Example-based Multilevel Dictionary learning (Ex-MLD).
Results:
The training set is the same as that described above. For the baseline Alt-Opt approach, we train a single dictionary with K=256 using 100 iterations with the sparsity penalty λtr set to 0.1. The ensemble learning procedures BoostEx and RandExAv are trained with L=50 and K=256 for sparsity penalty λtr=0.1. The boosted ensembles are trained by taking the random projection operator into consideration, as discussed in Section IV for the reduced measurements, N={8, 16, 32}. For the Ex-MLD method, both the number of levels and the number of atoms in each level were fixed at 16. In each level, we obtained L=50 dictionaries to compute the ensemble representation.
The recovery performance of the proposed ensemble models is evaluated using the set of standard images shown in Table I 1000 of FIG. 10. Each image 1001 is divided into non-overlapping patches of size 8×8, and random projection is performed with the number of measurements set at N={8, 16, 32}. For the Alt-Opt procedure, the individual patches are recovered using (15), and for the ensemble methods, the approximations computed using the L individual dictionaries are combined. The penalty λte is set to 0.1 for sparse coding in all cases. For each method, the PSNR values were obtained by averaging the results over 10 iterations with different random measurement matrices, and the results are reported in Table I. It was observed that the proposed ensemble methods outperform the Alt-Opt methods in all cases. In particular, we note that the simple RandExAv performs better than the boosting approach, although it was shown previously that boosting approaches show a superior performance. The reason for this discrepancy is that boosting aggressively reduces error with training data, and hence may lead to overfitting with degraded test data. Whereas, the RandExAv method provides the same importance to all individual approximations both during the training and the testing phases. As a result, it provides a better generalization in the presence of degradation. We also note that similar behavior has been observed with ensemble classification methods. Random sampling methods such as bagging perform better than boosting with noisy examples, since bagging exploits classification noise to produce more diverse classifiers. Furthermore, we observed that the proposed Ex-MLD method performed significantly better than all approaches, particularly for lower number of measurements. FIG. 5 shows images 502-510 recovered using the different approaches, when N was fixed at 8. As it can be observed, the Ex-MLD and RandExAv methods (i.e. 510 and 508) provide PSNR gains of 2.38 dB and 1.34 dB respectively, when compared to the Alt-Opt approach.
Single Image Superresolution
Single image superresolution (SISR) attempts to reconstruct a high-resolution image using just a single low-resolution image. It is a severely ill-posed problem and in sparse representation based approaches, the prior knowledge that natural image patches can be represented as a sparse linear combination of elementary patches, is used. The degraded test image is represented as Z=ΦY, where the operator Φ the blurs the high-resolution image Y and then downsamples it. Note that Y and Z denote vectorized high- and low-resolution images respectively. Each overlapping patch obtained from the degraded image is denoted as z. The paired dictionary learning procedure (PairDict) has been very effective in recovering the high-resolution patches. This method initially creates degraded counterparts of the high-resolution training images, following which gradient-based features are extracted from the low-resolution patches and the features are appended to the corresponding vectorized high-resolution patches. These augmented features are used to train a paired dictionary




   
  
    (
    
      
        
          
            D
            lo
          
        
      
      
        
          
            D
            hi
          
        
      
    
    )
  




such that each low-resolution and its corresponding high-resolution patches share the same sparse code. For a low-resolution test patch z, the sparse code a is obtained using Dlo, and the corresponding high-resolution counterpart is recovered as y=Dhia. An initial estimate Y0 of the high-resolution image is obtained by appropriately averaging the overlapped high-resolution patches. Finally, a global reconstruction constraint is enforced by projecting Y0 on to the solution space of ΦY=Z,





  
    
      
        
          
            min
            Y
          
          ⁢
          
            
              
              
                Z
                -
                
                  Φ
                  ⁢
                  
                    
                  
                  ⁢
                  Y
                
              
              
            
            2
            2
          
        
        +
        
          c
          ⁢
          
            
              
              
                Y
                -
                
                  Y
                  0
                
              
              
            
            2
            2
          
        
      
    
    
      
        (
        16
        )
      
    
  




to obtain the final reconstruction. As an alternative, in the example-based procedure (ExDict) the dictionaries Dlo and Dhi are directly fixed as the features extracted from low-resolution patches and vectorized high resolution patches respectively. Similar to the PairDict method, the global reconstruction constraint in (16) is imposed for the final reconstruction.
In our simulations, standard grayscale images 1004 (Table II 1002 of FIG. 10) are magnified by a factor of 2, using the proposed approaches. In addition to the PairDict and ExDict methods, simple bicubic interpolation is also used as a baseline method. We also obtained paired dictionaries with 1024 atoms using 100,000 randomly chosen patches of size 5×5 from the grayscale natural images in the training set. The sparsity penalty used in training was λtr=0.15. The training set was reduced to the size of 20,000 samples and used as the dictionary for the ExDict method. For ensemble learning, L was fixed at 50 and the approximation for each data sample was obtained using just a 1-sparse representation.
For different number of training samples, we compared the training times for PairDict (40 iterations) and BoostEx (L=50) algorithms in table 800 of FIG. 6. The computation times were obtained using a single core of a 2.8 GHz Intel i7 Linux machine with 8 GB RAM. The BoostEx approach 814 has lesser computational complexity for training compared to the PairDict approach 812. The ExDict procedure requires no training and for RandExAv, training time is just the time for randomly selecting K samples from the training set of T samples, for L rounds. Clearly, the complexity incurred for this is extremely low.
For the test images, SISR is performed using the baseline PairDict 812 and ExDict approaches using a sparsity penalty of λte=0.2. The recovery performance of the proposed algorithms are reported in Table II. For PairDict 812, as well the proposed ensemble methods, the dictionary size is fixed at 1024, whereas all the examples are used for training with the ExDict approach. We observed from our results that an ensemble representation with a simplified sparse coding scheme (1-sparse) matched the performance of the baseline methods, as illustrated in images 902, 904, 906, 908, and 910 of FIG. 7.
Several approaches exist for recovering images from their degraded observations. Recently, methods based on sparse representations have gained prominence and provide state-of-the-art performance. The use of sparse models is motivated by its ability to effectively represent natural images as a sparse linear combination of elementary patterns. These patterns are typically stacked together as columns of a “dictionary” matrix, and the dictionary can be inferred from training examples using a variety of learning algorithms.
The important challenges associated with the use of sparse models include: (i) computational complexity of obtaining sparse codes, (ii) scalability of dictionary learning algorithms, and (iii) non-availability of a large training dataset. Various embodiments of the current invention address these challenges by proposing to use an ensemble of multiple sparse models that are learned using low-complexity approaches. The approximation with each sparse model is obtained using a parallel-friendly “Correlate-and-Maximum” procedure. The final reconstruction is obtained by suitably aggregating the approximations from the individual models. The dictionary for each sparse model can be designed either using different random subsets from the training data (bagging), or sequentially such that the current dictionary emphasizes those data samples that were not well represented in the previous iteration (boosting).
Our results with two image recovery problems, compressive sensing and image superresolution, demonstrate that the simple recovery procedure with the proposed ensemble method matches or improve over the performance of several existing methods, at reduced computational complexity. Incorporating the ensemble dictionaries in hierarchical multilevel learning and obtaining multiple dictionaries in each level leads to further performance improvements.
Various embodiments of the present disclosure include some or all of the following:
1. The proposed scheme uses an ensemble of sparse models for recovering degraded images, where each model is obtained from uncorrupted training data and combined using bagging/boosting strategies.
2. Learning each model in the ensemble is a low-complex and simple procedure. The dictionary for each model is a random subset selected from the training examples and individual approximations with each model are obtained using simple correlation operations.
3. For a degraded observation, the computation of individual approximations can be parallelized as they involve only a simple “Correlate-and-Maximum” procedure. The final approximation is just an aggregation of the individual approximations with weights that are pre-computed in the training stage.
4. Depending on the application, the knowledge of degradation and the error from the previous approximation can be incorporated to “boost” the approximations and obtain better models. Recovery performance is further improved by incorporating the proposed ensemble models in the hierarchical multilevel dictionary learning algorithm.
5. The performance of existing image recovery algorithms that use sparse models is dependent on the availability of a rich set of training examples. Due to training complexity and data availability, a typical learning algorithm may only use a small set of samples to infer sparse models. By using an ensemble of weak sparse models, a number of embodiments address the statistical and computational challenges in learning, and thereby provide high fidelity image recovery.
6. When the form of degradation is known, it can be incorporated during the training process itself.
7. In contrast to other sparse modeling methods, the complex dictionary design procedure during training is replaced by a simple table-lookup process.
8. The proposed system may be highly scalable as it is parallel-friendly. Computation of the sparse code for a test sample using each model in the ensemble can be performed on a GPU. Sparse coding is simplified into a “Correlate-and-Maximum” procedure that can be easily fit in a single CUDA-GPU kernel. By storing the look-up tables in the GPU memory and performing all computations in the GPU, the full potential of GPUs can be harnessed enabling real time applications.
Advancement in imaging technologies is beneficial in several modern-day computing paradigms. Recovering images from different forms of degradation is a fundamental problem in a variety of imaging modalities. Efficient and low-complexity algorithms are important for commercial applications in mobile image processing and medical imaging, as examples. Availability of large-scale data has added another new dimension to computational challenges. The limited size and cost constraints of mobile devices prohibits the use of bulky and expensive optics and sensors for good imaging. Image recovery algorithms can bridge the gap between low-cost sensing and high-quality image rendering.
Although several image processing solutions exist for these problems, the more recent sparse modeling methods have had a significant impact. The improved performance of sparse models comes at the price of increase in complexity and dependency on a rich set of training examples. Aspects of the present disclosure propose the use of an ensemble of multiple “weak” sparse models to obtain improved estimates even under limited training conditions. Furthermore, by simplifying the sparse coding and dictionary inference processes to “Correlate-and-Maximum” procedures, certain embodiments make it viable and economical to use sparse models in commercial systems. In various embodiments, the proposed method can enhance low resolution/degraded images to high resolution images using a parallel implementation with multi-core CPUs and GPUs on mobile devices, desktops, and servers, for instance. This new technology can make a strong impact in the areas of HDTV, medical/satellite imaging, mobile image processing, and remote sensing, as examples.
FIG. 8 illustrates an exemplary embodiment of computer system 600, all of which or a portion of which can be suitable for implementing the techniques described above. As an example, a different or separate one of chassis 602 (and its internal components) can be suitable for implementing the techniques described above. Furthermore, one or more elements of computer system 600 (e.g., refreshing monitor 606, keyboard 604, and/or mouse 610, etc.) can also be appropriate for implementing the techniques described above. Computer system 600 comprises chassis 602 containing one or more circuit boards (not shown), Universal Serial Bus (USB) port 612, Compact Disc Read-Only Memory (CD-ROM) and/or Digital Video Disc (DVD) drive 616, and hard drive 614. A representative block diagram of the elements included on the circuit boards inside chassis 602 is shown in FIG. 9. Central processing unit (CPU) 710 in FIG. 9 is coupled to system bus 714 in FIG. 9. In various embodiments, the architecture of CPU 710 can be compliant with any of a variety of commercially distributed architecture families.
Continuing with FIG. 9, system bus 714 also is coupled to memory storage unit 708, where memory storage unit 708 comprises both read only memory (ROM) and random access memory (RAM). Non-volatile portions of memory storage unit 708 or the ROM can be encoded with a boot code sequence suitable for restoring computer system 600 (FIG. 8) to a functional state after a system reset. In addition, memory storage unit 708 can comprise microcode such as a Basic Input-Output System (BIOS). In some examples, the one or more memory storage units of the various embodiments disclosed herein can comprise memory storage unit 708, a USB-equipped electronic device, such as, an external memory storage unit (not shown) coupled to universal serial bus (USB) port 612 (FIGS. 6-7), hard drive 614 (FIGS. 6-7), and/or CD-ROM or DVD drive 616 (FIGS. 6-7). In the same or different examples, the one or more memory storage units of the various embodiments disclosed herein can comprise an operating system, which can be a software program that manages the hardware and software resources of a computer and/or a computer network. The operating system can perform basic tasks such as, for example, controlling and allocating memory, prioritizing the processing of instructions, controlling input and output devices, facilitating networking, and managing files. Some examples of common operating systems can comprise Microsoft® Windows® operating system (OS), Mac® OS, UNIX® OS, and Linux® OS.
As used herein, “processor” and/or “processing module” means any type of computational circuit, such as but not limited to a microprocessor, a microcontroller, a controller, a complex instruction set computing (CISC) microprocessor, a reduced instruction set computing (RISC) microprocessor, a very long instruction word (VLIW) microprocessor, a graphics processor, a digital signal processor, or any other type of processor or processing circuit capable of performing the desired functions. In some examples, the one or more processors of the various embodiments disclosed herein can comprise CPU 710.
In the depicted embodiment of FIG. 9, various I/O devices such as disk controller 704, graphics adapter 724, video controller 702, keyboard adapter 726, mouse adapter 706, network adapter 720, and other I/O devices 722 can be coupled to system bus 714. Keyboard adapter 726 and mouse adapter 706 are coupled to keyboard 604 (FIGS. 8 and 9) and mouse 610 (FIGS. 8 and 9), respectively, of computer system 600 (FIG. 8). While graphics adapter 724 and video controller 702 are indicated as distinct units in FIG. 9, video controller 702 can be integrated into graphics adapter 724, or vice versa in other embodiments. Video controller 702 is suitable for refreshing monitor 606 (FIGS. 8 and 9) to display images on a screen 608 (FIG. 8) of computer system 600 (FIG. 8). Disk controller 704 can control hard drive 614 (FIGS. 8 and 9), USB port 612 (FIGS. 8 and 9), and CD-ROM drive 616 (FIGS. 8 and 9). In other embodiments, distinct units can be used to control each of these devices separately.
In some embodiments, network adapter 720 can comprise and/or be implemented as a WNIC (wireless network interface controller) card (not shown) plugged or coupled to an expansion port (not shown) in computer system 600 (FIG. 8). In other embodiments, the WNIC card can be a wireless network card built into computer system 600 (FIG. 8). A wireless network adapter can be built into computer system 600 by having wireless communication capabilities integrated into the motherboard chipset (not shown), or implemented via one or more dedicated wireless communication chips (not shown), connected through a PCI (peripheral component interconnector) or a PCI express bus of computer system 600 (FIG. 8) or USB port 612 (FIG. 8). In other embodiments, network adapter 720 can comprise and/or be implemented as a wired network interface controller card (not shown).
Although many other components of computer system 600 (FIG. 8) are not shown, such components and their interconnection are well known to those of ordinary skill in the art. Accordingly, further details concerning the construction and composition of computer system 600 and the circuit boards inside chassis 602 (FIG. 8) are not discussed herein.
When computer system 600 in FIG. 8 is running, program instructions stored on a USB-equipped electronic device connected to USB port 612, on a CD-ROM or DVD in CD-ROM and/or DVD drive 616, on hard drive 614, or in memory storage unit 708 (FIG. 9) are executed by CPU 710 (FIG. 9). A portion of the program instructions, stored on these devices, can be suitable for carrying out at least part of the techniques described above.
Although computer system 600 is illustrated as a desktop computer in FIG. 8, there can be examples where computer system 600 may take a different form factor while still having functional elements similar to those described for computer system 600. In some embodiments, computer system 600 may comprise a single computer, a single server, or a cluster or collection of computers or servers, or a cloud of computers or servers. Typically, a cluster or collection of servers can be used when the demand on computer system 600 exceeds the reasonable capability of a single server or computer.
Various embodiments of the subject matter described herein include various combinations of the acts, structure, components, and features described herein, shown in the drawings, or that are known in the art. Moreover, certain procedures can include acts such as manufacturing, obtaining, or providing components that perform functions described herein or in the documents that are incorporated by reference. Furthermore, various embodiments include advertising and selling products that perform functions described herein, that contain structure described herein, or that include instructions to perform functions described herein, as examples. Such products may be obtained or provided through distributors, dealers, or over the Internet, for instance. The subject matter described herein also includes various means for accomplishing the various functions or acts described herein, in the documents that are incorporated by reference, or that are apparent from the structure and acts described.
Further, as used herein, the word “or”, except where indicated otherwise, does not imply that the alternatives listed are mutually exclusive. Even further, where alternatives are listed herein, it should be understood that in some embodiments, fewer alternatives may be available, or in particular embodiments, just one alternative may be available, as examples.

## Claims

1. A method for performing image recovery comprising:
generating, using at least one processor, an ensemble of sparse models from uncorrupted training data, comprising:
learning a dictionary from the uncorrupted training data by randomly selecting a subset of the uncorrupted training data to form the dictionary, and
learning at least one weight from the uncorrupted training data using an unconstrained least squares approximation, the at least one weight corresponding to the dictionary;

using respective sparse models of the ensemble of sparse models to generate a plurality of individual approximations of a first image, wherein each individual approximation of the plurality of individual approximations corresponds to a particular sparse model of the respective sparse models; and
aggregating, using the at least one processor, the individual approximations to generate a second image, by
calculating a linear combination of the plurality of individual approximations with pre-determined weights applied to each individual approximation of the plurality of individual approximations.

2. The method of claim 1, wherein the first image is a degraded image and the second image is an improved image, wherein the respective sparse models are individual weak sparse models, and wherein aggregating the individual approximations comprises calculating a weighted average of the individual approximations.

3. The method of claim 1, wherein generating an ensemble of sparse models comprises learning at least one dictionary and at least one weight corresponding to the at least one dictionary from the uncorrupted training data sequentially, using a greedy forward selection procedure.

4. The method of claim 1, wherein generating an ensemble of sparse models comprises:
degrading uncorrupted training data; and
learning at least one dictionary from the degraded uncorrupted training data.

5. The method of claim 1, wherein each individual approximation of the individual approximations is determined using a low-complexity correlate-and-maximum procedure and wherein the sparse models represent images parsimoniously using elementary patterns from a dictionary matrix.

6. A system for performing image recovery comprising: at least one processor to:
generate an ensemble of sparse models from uncorrupted training data, comprising
learning a dictionary from the uncorrupted training data by randomly selecting a subset of the uncorrupted training data to form the dictionary; and
learning at least one weight from the uncorrupted training data using an unconstrained least squares approximation, the at least one weight corresponding to the dictionary;

use respective sparse models of the ensemble of sparse models to generate a plurality of individual approximations of a first image, wherein each individual approximation of the plurality of individual approximations corresponds to a particular sparse model of the respective sparse models; and
aggregate the individual approximations to generate a second image,
wherein each individual approximation of the individual approximations is determined using a low-complexity correlate-and-maximum procedure, and wherein the sparse models represent images parsimoniously using elementary patterns from a dictionary matrix.

7. The system of claim 6, wherein aggregating the individual approximations comprises calculating a linear combination of the plurality of individual approximations with pre-determined weights applied to each individual approximation of the plurality of individual approximations.

8. The system of claim 6, wherein the first image is a degraded image and the second image is an improved image, wherein the respective sparse models are individual weak sparse models, and wherein aggregating the individual approximations comprises calculating a weighted average of the individual approximations.

9. The system of claim 6, wherein generating an ensemble of sparse models comprises learning at least one dictionary and at least one weight corresponding to the at least one dictionary from the uncorrupted training data sequentially, using a greedy forward selection procedure.

10. The system of claim 6, wherein generating an ensemble of sparse models comprises:
degrading uncorrupted training data; and
learning at least one dictionary from the degraded uncorrupted training data.

11. A non-transitory computer-readable medium encoded with instructions for performing image recovery, the instructions executable by a processor, comprising:
generating an ensemble of sparse models from uncorrupted training data, comprising
learning a dictionary from the uncorrupted training data by randomly selecting a subset of the uncorrupted training data to form the dictionary; and learning at least one weight from the uncorrupted training data using an unconstrained least squares approximation, the at least one weight corresponding to the dictionary

using respective sparse models of the ensemble of sparse models to generate a plurality of individual approximations of a first image, wherein each individual approximation of the plurality of individual approximations corresponds to a particular sparse model of the respective sparse models; and
aggregating the individual approximations to generate a second image,
wherein the first image is a degraded image and the second image is an improved image, wherein the respective sparse models are individual weak sparse models, and wherein aggregating the individual approximations comprises calculating a weighted average of the individual approximations.

12. The non-transitory computer-readable medium of claim 11, wherein aggregating the individual approximations comprises calculating a linear combination of the plurality of individual approximations with pre-determined weights applied to each individual approximation of the plurality of individual approximations.

13. The non-transitory computer-readable medium of claim 11, wherein generating an ensemble of sparse models comprises:
learning at least one dictionary from the uncorrupted training data by randomly selecting a subset of the uncorrupted training data to form the at least one dictionary; and
learning at least one weight from the uncorrupted training data using an unconstrained least squares approximation, the at least one weight corresponding to the at least one dictionary.

14. The non-transitory computer-readable medium of claim 11, wherein generating an ensemble of sparse models comprises learning at least one dictionary and at least one weight corresponding to the at least one dictionary from the uncorrupted training data sequentially, using a greedy forward selection procedure.

15. The non-transitory computer-readable medium of claim 11, wherein each individual approximation of the individual approximations is determined using a low-complexity correlate-and-maximum procedure, and wherein the sparse models represent images parsimoniously using elementary patterns from a dictionary matrix.