# Patent

## Title
Unsupervised representation learning with contrastive prototypes

## Abstract

The system and method are directed to a prototypical contrastive learning (PCL). The PCL explicitly encodes the hierarchical semantic structure of the dataset into the learned embedding space and prevents the network from exploiting low-level cues for solving the unsupervised learning task. The PCL includes prototypes as the latent variables to help find the maximum-likelihood estimation of the network parameters in an expectation-maximization framework. The PCL iteratively performs an E-step for finding prototypes with clustering and M-step for optimizing the network on a contrastive loss.

## Background
RELATED APPLICATIONS
The present application is a continuation of U.S. application Ser. No. 16/870,621, filed May 8, 2020, now allowed, and also claims priority to U.S. Provisional Patent Application No. 62/992,004, filed Mar. 19, 2020, which is incorporated by reference herein in its entirety.


TECHNICAL FIELD
The embodiments relate generally to representation learning, and more specifically to a prototypical contrastive learning using a prototypical contrastive learning framework.
BACKGROUND
Most existing techniques for unsupervised or self-supervised learning fall into three categories: instance-wise contrastive learning, unsupervised clustering, and learning from pre-tasks.
The instance wise contrastive learning may learn an embedding space where samples (e.g. crops) from the same instance (e.g. an image) are pulled closer and samples from different instances are pushed apart. To construct the contrastive loss for a mini-batch of samples, instance-wise contrastive learning may generate positive instance features and negative instance features for each sample. To generate the instance features, the instance-wise contrastive learning may use a memory bank approach that stores the features of all samples calculated in the previous step and selects features from the memory bank to form positive and negative pairs. To generate instance features, the instance-wise contrastive learning may also use an end-to-end approach that generates instance features using all samples within the current mini-batch. In the end-to-end approach, two encoders may be applied to the mini-batch, one encoder on the original samples, and the other encoder on the augmented version of the samples.
The instance-wise contrastive learning has various drawbacks. The task of instance discrimination may be solved by exploiting low-level differences between images. This means that the learned embeddings may not capture higher-level semantic information. Also, instance-wise learning samples a sufficiently large number of negative instances, which inevitably yields negative pairs that share similar semantic meaning and should be closer in the embedding space. However, the negative pairs are undesirably pushed apart by the contrastive loss. In other words, the instance wise discrimination learns an embedding space that only preserves the local smoothness around each instance but ignores the global semantic structure of the dataset.
The techniques based on clustering jointly learn image embeddings and cluster assignments. The clustering techniques, however may not learn transferable representations from a large scale of images. The DeepCluster clustering technique may learn from millions of images by performing iterative clustering and unsupervised representation learning. In DeepCluster, the cluster assignments may be pseudo-labels and a classification objective may be optimized. The deficiencies of this approach are that the high-dimensional features from the penultimate layer of a ConvNet are not optimal for clustering and need to be PCA-reduced. Another deficiency is that an additional linear classification layer is frequently re-initialized which interferes with representation learning.
Other techniques are self-supervised learning methods that focus on training deep neural networks to solve pretext tasks that lead to good image representations being learned. These tasks usually involve hiding certain information about the input and training the network to recover the missing information. Example tasks may include image inpainting, colorization, prediction of patch orderings and image transformations. However, these pretext tasks exploit specific structures of visual data and may be difficult to generalize to other domains.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a simplified diagram of a computing device for implementing a prototypical contrastive learning framework, according to some embodiments.
FIG. 2 is block diagram of a prototypical contrastive learning framework, according to some embodiments.
FIG. 3 is a flowchart of a method for training the prototypical contrastive learning framework, according to some embodiments.
FIG. 4 is a diagram of a histogram illustrating cluster sizes in a prototypical contrastive learning framework with fixed and estimated concentrations, according to some embodiments.
FIG. 5 is a diagram of an Expectation-Maximization algorithm for training the prototypical contrastive learning framework, according to some embodiments.
FIG. 6 is a diagram that illustrates clustering in the E-step of the Expectation-Maximization algorithm, according to some embodiments.



DETAILED DESCRIPTION
Unsupervised visual representation learning aims to learn image representations from pixels themselves and without relying on semantic annotations. Recent developments in unsupervised representation learning are largely driven by the task called instance discrimination. Methods based on the instance discrimination usually consist of two key components: image transformation and contrastive loss. Image transformation aims to generate multiple embeddings that represent the same image, by data augmentation, patch perturbation, or by using momentum features. The contrastive loss, which may be in the form of a noise contrastive estimator, aims to bring closer samples from the same instance and separate samples from different instances. Instance-wise contrastive learning leads to an embedding space where all instances are well-separated, and each instance is locally smooth (i.e. input perturbation leads to similar representations).
Despite their improved performance, methods based on instance discrimination share a common fundamental weakness: the semantic structure of data is not encoded by the learned representations. This problem arises because instance-wise contrastive learning considers two samples to be a negative pair as long as they are from different instances, regardless of the semantic similarity between instances. The problem is magnified by the fact that thousands of negative samples are generated to form the contrastive loss, leading to many negative pairs that share similar semantic meaning but are undesirably pushed apart in the embedding space.
The embodiments are directed to a prototypical contrastive learning (PCL) framework for unsupervised representation learning. The PCL framework explicitly encodes the semantic structure into the embedding space. A prototype in the PCL framework may be defined as a representative embedding for a group of semantically similar instances. Each instance may be assigned to several prototypes of different granularity. The PCL framework may also construct a contrastive loss which may enforce the embedding of a sample to be more similar to its assigned prototypes compared to other prototypes. In practice, the PCL framework may find prototypes by performing standard clustering on the embeddings.
In some embodiments, the PCL framework may use a bilevel Expectation-Maximization (E-M) algorithm. The E-M algorithm may find parameters of a deep neural network (DNN) that best describe the data by iteratively approximating and maximizing the likelihood function. The E-M algorithm may include additional latent variables, such as prototypes and instance assignments. The E-M algorithm may estimate the latent variables in the E-step of the E-M algorithm by performing k-means clustering. In the M-step of the E-M algorithm, the E-M algorithm may update the network parameters by minimizing the proposed contrastive loss. A proposed contrastive loss may be determined using a ProtoNCE function, described below. The E-M algorithm may determine that minimizing the ProtoNCE function is equivalent to maximizing the approximate likelihood function under the assumption that the data distribution around each prototype is an isotropic Gaussian. By using the E-M algorithm, the widely used instance discrimination task can be explained as a special case of the PCL framework, where the prototype for each instance is its augmented feature, and the Gaussian distribution around each prototype has the same fixed variance.
The embodiments of the disclosure are directed to the PCL framework for unsupervised representation learning. The learned representation not only preserves the local smoothness of each image instance, but also captures the hierarchical semantic structure of the global dataset. Further, although described with respect to images, the PCL framework may also apply to any type of unstructured data such as video, text, speech, etc.
The embodiments of the disclosure are directed to the PCL framework that includes an Expectation-Maximization (E-M) algorithm. In the E-M algorithm the iterative steps of clustering and representation learning can be interpreted as approximating and maximizing the log-likelihood function.
The embodiments of the disclosure are also directed to using the ProtoNCE function for determining the contrastive loss. Notably, the ProtoNCE function dynamically estimates the concentration for the feature distribution around each prototype. The learned prototypes contain more information about the image classes.
FIG. 1 is a simplified diagram of a computing device 100 for implementing a prototypical contrastive learning framework, according to some embodiments. As shown in FIG. 1, computing device 100 includes a processor 110 coupled to memory 120. Operation of computing device 100 is controlled by processor 110. And although computing device 100 is shown with only one processor 110, it is understood that processor 110 may be representative of one or more central processing units, multi-core processors, microprocessors, microcontrollers, digital signal processors, field programmable gate arrays (FPGAs), application specific integrated circuits (ASICs), graphics processing units (GPUs) and/or the like in computing device 100. Computing device 100 may be implemented as a stand-alone subsystem, as a board added to a computing device, and/or as a virtual machine.
Memory 120 may be used to store software executed by computing device 100 and/or one or more data structures used during operation of computing device 100. Memory 120 may include one or more types of machine readable media. Some common forms of machine readable media may include floppy disk, flexible disk, hard disk, magnetic tape, any other magnetic medium, CD-ROM, any other optical medium, punch cards, paper tape, any other physical medium with patterns of holes, RAM, PROM, EPROM, FLASH-EPROM, any other memory chip or cartridge, and/or any other medium from which a processor or computer is adapted to read.
Processor 110 and/or memory 120 may be arranged in any suitable physical arrangement. In some embodiments, processor 110 and/or memory 120 may be implemented on a same board, in a same package (e.g., system-in-package), on a same chip (e.g., system-on-chip), and/or the like. In some embodiments, processor 110 and/or memory 120 may include distributed, virtualized, and/or containerized computing resources. Consistent with such embodiments, processor 110 and/or memory 120 may be located in one or more data centers and/or cloud computing facilities.
In some embodiments, memory 120 may include non-transitory, tangible, machine readable media that includes executable code that when run by one or more processors (e.g., processor 110) may cause the one or more processors to perform the methods described in further detail herein. For example, as shown, memory 120 includes instructions for prototypical contrastive learning (PCL) framework 130 that may be used to implement and/or emulate the systems and models, and/or to implement any of the methods described further herein. The PCL framework 130 may be a “network” that may comprise any hardware or software-based framework that includes any artificial intelligence network or system, neural network or system and/or any training or learning models implemented thereon or therewith. The PCL framework 130 may include encoders that may be trained using images 140 (or other unstructured data such as video, speech, text, etc.) that the PCL framework 130 receives as input. The images 140 may comprise multiple pixels. Unlike conventional frameworks, the PCL framework 130 may be trained using images that do not include labels or tags that may identify different features of the images. After the PCL framework 130 is trained, an encoder in the PCL framework 130 may perform one or more tasks, e.g. identify a picture that is included in an image, determine image colorization, predict patch orderings, etc. In a non-structured data example, once trained an encoder in the PCL framework 130 may perform one or more tasks, namely generate clusters that include similar subsets of data in unstructured data.
In some embodiments, the PCL framework 130 may be a neural network that includes one or more components. FIG. 2 is a block diagram 200 of a PCL framework according to some embodiments. The PCL framework 130 includes a momentum encoder 205, a clustering module 235, an encoder 210, and a ProtoNCE module 215. Some or all of momentum encoder 205, clustering module 235, encoder 210, and ProtoNCE module 215 may be structured as a neural network.
An expectation-maximization (E-M) algorithm 220 may act on the momentum encoder 205, clustering module 235, encoder 210, and ProtoNCE module 215. The E-M algorithm 220 may execute during multiple iterations that occur during a configurable time period called an epoch. Further, there may be multiple epochs during which the E-M algorithm 220 may execute and train encoder 210. The E-M algorithm 220 may be divided into an E-step 225 and a M-step 230, which are both performed at each iteration. In the E-step 225, the momentum encoder 205 may identify different features in images 140 and the clustering module 235 may generate a configurable number of clusters based on the identified features. Each cluster may include a prototype, assignment and concentration of the similar features in the images 140 and may have different levels of granularity. The M-step 230 may include the encoder 210 and ProtoNCE module 215 and may train the PCL framework 130. For example, the ProtoNCE module 215 may approximate and maximize a likelihood function that is back propagated to the encoder 210. Encoder 210 is trained using the likelihood function. The momentum encoder 205 is then updated with weights of the trained encoder 210. Both E-step 225 and M-step 230 of the E-M algorithm 220 are discussed in detail below.
Each of the momentum encoder 205 and encoder 210 may be a neural network, such as a convolutional neural network. Momentum encoder 205 and encoder 210 may be structural copies of each other. In other words, momentum encoder 205 and encoder 210 may be two instances of the same neural network, but may have different weights that are assigned to the nodes of the neural network.
As discussed above, during each iteration in the epoch, the E-M algorithm 220 performs the E-step 225 and the M-step 230. At the beginning of E-step 225, momentum encoder 205 receives one or more images 140 as input. In some instances, prior to the momentum encoder 205 receiving images 140, images 140 may be augmented, e.g. cropped, color changed, etc. Momentum encoder 205 passes the received images 140 though the neural network to determine the features of images 140. The features may be the output of the last layer of the convolutional neural network that makes up the momentum encoder 205. The features may be embeddings of the momentum encoder 205.
In some embodiments, clustering module 235 may receive the features that are the output of momentum encoder 205. Clustering module 235 may cluster the features into one or more clusters, such as clusters 240A-C. Each of the clusters 240A, 240B, and 240C may be associated with a prototype. Prototypes C are shown in FIG. 2 in the center of each cluster 240A-C. Each of the clusters 240A, 240B, and 240C includes features that are close or similar to each other. The prototypes C may have different granularity. For example, cluster 240B may have a coarse-grained prototype and clusters 240A and 240C may have a fine-grained prototype. The output 245 of clustering module 235 may be prototype C, assignment S, and concentration M for each cluster 240A-C. The prototype C may be a vector that identifies a center of a cluster, such as cluster 240A, 240B or 240C. Assignment S may identify which features are assigned to each cluster 240A, 240B and/or 240C. Concentration M estimates how close or tight the features are within each cluster 240A, 240B, and 240C. For example, concentration M in cluster 240A may indicate a tight concentration, while concentration M in cluster 240B may indicate a loose concentration. Notably, feature of the same image in images 140 may belong to different clusters 240A-C.
During the M-step 230, the ProtoNCE module 215 may receive output 245 that includes prototypes C, concentrations M, and assignments S for clusters 240A-C determined during the E-step 225. The ProtoNCE module 215 may use the prototypes C, concentrations M, and assignments S to optimize the ProtoNCE function shown in Equation 6 below. The ProtoNCE module 215 may determine that minimizing the ProtoNCE function is equivalent to maximizing the approximate likelihood function under the assumption that the data distribution around each prototype is an isotropic Gaussian. The optimized ProtoNCE function may be back propagated from ProtoNCE module 215 to encoder 210.
Encoder 210 may be trained using the ProtoNCE function determined by the ProtoNCE module 215 and the images 140. For example, encoder 210 may receive and process images 140 while applying the ProtoNCE function to the weights. As discussed above, the data in the images may be augmented. During training, the weights of the encoder 210 are trained. The trained weights of encoder 210 may then update the weights of the momentum encoder 205. Because the structures of encoder 210 and momentum encoder 205 are the same, the weights from encoder 210 may update the weights of momentum encoder 205 in the same position in the structure.
In some embodiments, the weights of the momentum encoder 205 may be updated by determining an average or a moving average of the weights of the momentum encoder 205 and the weights of encoder 210. Once the weights of momentum encoder 205 are updated, the PCL framework 130 may begin the next iteration in the epoch. During the next iteration, the PCL framework 130 may be trained using the same or different images 140 to determine clusters 240A-C, the prototypes C, concentrations M, and assignments S, that are then used to further optimize the ProtoNCE function.
The iterative process may continue for a preconfigured number of epochs. Once the iterative process discussed above is concluded, the momentum encoder 205 is trained and may be applied to perform various tasks.
FIG. 3 is a simplified diagram of a method 300 for training the PCL framework, according to some embodiments. One or more of the processes 302-312 of method 300 may be implemented, at least in part, in the form of executable code stored on non-transitory, tangible, machine-readable media that when run by one or more processors may cause the one or more processors to perform one or more of the processes 302-312. The processes 302-312 may be performed at each iteration of an epoch that may include a configurable number of iterations. Further, processes 302-306 may be performed as part of the E-step 225 of the E-M algorithm 220 and processes 308-312 may be performed as part of the M-step 230.
At process 302, data from the images is received. For example, momentum encoder 205 may receive data from one or more images 140. The data from images 140 may be cropped or otherwise augmented. The data may include pixels from images 140.
At process 304, features from the images are determined. For example, momentum encoder 205, which may be structured as a convolutional neural network, may generate embeddings which are features of images 140.
At process 306, prototypes, assignments and concentrations are determined from the features. For example, clustering module 235 may receive the features determined in process 304 and generate clusters, such as clusters 240A-C, using the features. From the clusters 240A-C, clustering module 235 may determine prototypes C, assignments S, and concentrations M of each cluster 240A, 240B, and 240C. In some embodiments, the number of clusters that clustering module 235 may determine is preconfigured, and clustering module 235 determines which features are included in which one or more clusters 240A-240C.
At process 308, a ProtoNCE function is determined. For example, the ProtoNCE module 215 receives the prototypes C, assignments S, and concentration M and determines the ProtoNCE function that minimizes a proposed contrastive loss. For example, ProtoNCE module 215 may determine the ProtoNCE function by maximizing the approximate likelihood function under the assumption that the data distribution around each prototype is an isotropic Gaussian.
At process 310, an encoder is trained. For example, encoder 210 is trained using the embedding function determined in step 308 and all or a subset of images 140. Like momentum encoder 205, encoder 210 may also receive images 140. During training, the contrastive loss function may be applied to the one or more weights of the nodes in the convolution neural network included in the encoder 210 as the encoder determines features of images 140.
At process 312, the momentum encoder is updated. For example, momentum encoder 205 may be updated with the weights of encoder 210. In some embodiments, momentum encoder 205 may be updated with an average of the weights of momentum encoder 205 and weights of encoder 210. In other embodiments, momentum encoder 205 may be updated with a moving average of the weights by taking an arithmetic mean of the weights of momentum encoder 205 at previous iterations and the weights received from encoder 210.
After process 312 completes, method 300 may repeat another iteration of processes 302-312 until the iterations complete the epoch. At completion of an epoch, the method 300 may repeat for another epoch or for a configurable number of epochs. Once the PCL framework 130 completes training, encoder 210 may be included in other frameworks, including other image processing frameworks to perform different tasks.
Going back to FIG. 2, the embodiments below, describe the mathematical framework included in various components of the PCL framework 130. As discussed above, the PCL framework 130 may be trained using images 140. Given a training set X={x1, x2, . . . , xN} of N images, which may be images 140, the unsupervised visual representation learning aims to learn an embedding function ƒθ (realized via a DNN) that maps X to V={v1, v2, . . . , vN} with vi=ƒθ(xi), such that vi best describes xi.
The instance-wise contrastive learning may achieve this objective by optimizing a contrastive loss function, such as an InfoNCE function. The InfoNCE function may be defined as:



 
  
   
    
     
      L
      InfoNCE
     
     =
     
      
       ∑
       
        i
        =
        1
       
       N
      
       
      
       
        -
        log
       
       ⁢
       
        
         exp
         ⁡
         (
         
          
           
            v
            i
           
           ·
           
            v
            i
            ′
           
          
          /
          τ
         
         )
        
        
         
          ∑
          
           j
           =
           0
          
          r
         
          
         
          
           exp
           ⁡
           (
           
            
             
              v
              i
             
             ·
             
              v
              j
              ′
             
            
            /
            τ
           
           )
          
          ′
         
        
       
      
     
    
   
   
    
     Equation
     ⁢
         
     1
    
   
  
 



where v′i is a positive embedding for instance i, v′j includes one positive embedding, r is negative embeddings for other instances, and r is a temperature hyper-parameter. These embeddings are obtained by feeding xi to momentum encoder 205 parametrized by θ′, v′i=ƒθ′(xi), where θ′ is a moving average of θ.

In the PCL framework 130, the prototypes c may replace v′ and a concentration estimation μ (shown as concentration M in FIG. 2) may replace the fixed temperature hyper-parameter τ. Further, the clustering and representation learning may be performed iteratively at each epoch.
The PCL framework 130 may find the network parameters θ that maximizes the likelihood function of the observed n samples:



 
  
   
    
     
      θ
      *
     
     =
     
      
       
        arg
        ⁢
          
        max
       
       θ
      
      ⁢
      
       
        ∑
        
         i
         =
         1
        
        n
       
       
        log
        ⁢
          
        
         p
         ⁡
         (
         
          
           x
           i
          
          ;
          θ
         
         )
        
       
      
     
    
   
   
    
     Equation
     ⁢
         
     2
    
   
  
 



Further, the observed data {xi}i=1n is related to the latent variable C={ci}i=1K which denotes the prototypes C of the data. In this way, the likelihood function may be re-written as:



 
  
   
    
     
      θ
      *
     
     =
     
      
       
        
         arg
         ⁢
           
         max
        
        θ
       
       ⁢
       
        
         ∑
         
          i
          =
          1
         
         n
        
        
         log
         ⁢
           
         
          p
          ⁡
          (
          
           
            x
            i
           
           ;
           θ
          
          )
         
        
       
      
      =
      
       
        
         arg
         ⁢
           
         max
        
        θ
       
       ⁢
       
        
         ∑
         
          i
          =
          1
         
         n
        
        
         log
         ⁢
         
          
           ∑
           
            
             c
             i
            
            ∈
            C
           
          
          
           p
           ⁡
           (
           
            
             x
             i
            
            ,
            
             
              c
              i
             
             ;
             θ
            
           
           )
          
         
        
       
      
     
    
   
   
    
     Equation
     ⁢
         
     3
    
   
  
 



In order to optimize the function in Equation 3, the PCL framework 130 may use a surrogate function to lower-bound Equation 3, as follows:



 
  
   
    
     
      
       ∑
       
        i
        =
        1
       
       n
      
      
       log
       ⁢
       
        
         ∑
         
          
           c
           i
          
          ∈
          C
         
        
        
         p
         ⁡
         (
         
          
           x
           i
          
          ,
          
           
            c
            i
           
           ;
           θ
          
         
         )
        
       
      
     
     =
     
      
       
        ∑
        
         i
         =
         1
        
        n
       
       
        log
        ⁢
        
         
          ∑
          
           
            c
            i
           
           ∈
           C
          
         
         
          
           Q
           ⁡
           (
           
            c
            i
           
           )
          
          ⁢
          
           
            p
            ⁡
            (
            
             
              x
              i
             
             ,
             
              
               c
               i
              
              ;
              θ
             
            
            )
           
           
            Q
            ⁡
            (
            
             c
             i
            
            )
           
          
         
        
       
      
      ≥
      
       
        ∑
        
         i
         =
         1
        
        n
       
       
        
         ∑
         
          
           c
           i
          
          ∈
          C
         
        
        
         
          Q
          ⁡
          (
          
           c
           i
          
          )
         
         ⁢
         log
         ⁢
         
          
           p
           ⁡
           (
           
            
             x
             i
            
            ,
            
             
              c
              i
             
             ;
             θ
            
           
           )
          
          
           Q
           ⁡
           (
           
            c
            i
           
           )
          
         
        
       
      
     
    
   
   
    
     Equation
     ⁢
        
     4
    
   
  
 



where Q(c1) denotes some distribution over c's (Σci∈C Q(ci)=1), and the last step of derivation uses Jensen's inequality. Further, to make the inequality hold with equality, the PCL framework 130 may require for




 
  
   
    
     
      p
      ⁡
      (
      
       
        x
        i
       
       ,
       
        
         c
         i
        
        ;
        θ
       
      
      )
     
     
      Q
      ⁡
      (
      
       c
       i
      
      )
     
    
   
   
     
   
  
 



to be a constant. In this case:




 
  
   
    
     
      Q
      ⁡
      (
      
       c
       i
      
      )
     
     =
     
      
       
        p
        ⁡
        (
        
         
          x
          i
         
         ,
         
          
           c
           i
          
          ;
          θ
         
        
        )
       
       
        
         ∑
         
          
           c
           i
          
          ∈
          C
         
        
        
         p
         ⁡
         (
         
          
           x
           i
          
          ,
          
           
            c
            i
           
           ;
           θ
          
         
         )
        
       
      
      =
      
       
        
         p
         ⁡
         (
         
          
           x
           i
          
          ,
          
           
            c
            i
           
           ;
           θ
          
         
         )
        
        
         p
         ⁡
         (
         
          
           x
           i
          
          ;
          θ
         
         )
        
       
       =
       
        p
        ⁡
        (
        
         
          
           c
           i
          
          ;
          
           x
           i
          
         
         ,
          
         θ
        
        )
       
      
     
    
   
   
    
     Equation
     ⁢
        
     5
    
   
  
 



Further, by ignoring the constant −Σi=1nΣci∈C Q(ci)log Q(ci) in Equation 4, the PCL framework 130 may maximize:



 
  
   
    
     
      ∑
      
       i
       =
       1
      
      n
     
     
      
       ∑
       
        
         c
         i
        
        ∈
        C
       
      
      
       
        Q
        ⁡
        (
        
         c
         i
        
        )
       
       ⁢
       log
       ⁢
       
        p
        ⁡
        (
        
         
          x
          i
         
         ,
         
          
           c
           i
          
          ;
          θ
         
        
        )
       
      
     
    
   
   
    
     Equation
     ⁢
        
     6
    
   
  
 



During the E-step 225 of the E-M algorithm 220, the PCL framework 130 aims to estimate p(ci; xi, θ). To achieve this, the clustering module 235 may perform k-means clustering on the features v′i=ƒθ′(xi) of images 140 identified by momentum encoder 205 to obtain k clusters. Prototype ci may be defined as cluster centroid for the i-the cluster. Then the clustering module 235 computes p(ci; xi, θ)= (xi∈ci), where  (xi∈ci)=1 if xi belongs to the cluster represented by ci. Otherwise  (xi∈ci)=0.
During the M-step 230, the ProtoNCE module 215 maximizes the lower-bound of Equation 6 as follows:



 
  
   
    
     
      
       ∑
       
        i
        =
        1
       
       n
      
      
       
        ∑
        
         
          c
          i
         
         ∈
         C
        
       
       
        
         Q
         ⁡
         (
         
          c
          i
         
         )
        
        ⁢
        log
        ⁢
          
        
         p
         ⁡
         (
         
          
           x
           i
          
          ,
          
           
            c
            i
           
           ;
           θ
          
         
         )
        
       
      
     
     =
     
      
       
        ∑
        
         i
         =
         1
        
        n
       
       
        
         ∑
         
          
           c
           i
          
          ∈
          C
         
        
        
         
          p
          ⁡
          (
          
           
            c
            i
           
           ;
           
            
             x
             
              i
              ⁢
              ′
             
            
            ⁢
            θ
           
          
          )
         
         ⁢
         log
         ⁢
           
         
          p
          ⁡
          (
          
           
            x
            i
           
           ,
           
            
             c
             i
            
            ;
            θ
           
          
          )
         
        
       
      
      =
      
       
        ∑
        
         i
         =
         1
        
        n
       
       
        
         ∑
         
          
           c
           i
          
          ∈
          C
         
        
        
         1
         ⁢
         
          (
          
           
            x
            i
           
           ∈
           
            c
            i
           
          
          )
         
         ⁢
         log
         ⁢
           
         
          p
          ⁡
          (
          
           
            x
            i
           
           ,
           
            
             c
             i
            
            ;
            θ
           
          
          )
         
        
       
      
     
    
   
   
    
     Equation
     ⁢
        
     7
    
   
  
 



Under the assumption of a uniform prior over cluster centroids, the p(xi, ci; θ) may be represented as follows:



 
  
   
    
     
      p
      ⁡
      (
      
       
        x
        i
       
       ,
       
        
         c
         i
        
        ;
        θ
       
      
      )
     
     =
     
      
       
        p
        ⁡
        (
        
         
          
           x
           i
          
          ;
          
           c
           i
          
         
         ,
         θ
        
        )
       
       ⁢
       
        p
        ⁡
        (
        
         
          c
          i
         
         ;
         θ
        
        )
       
      
      =
      
       
        1
        k
       
       ·
       
        p
        ⁡
        (
        
         
          
           x
           i
          
          ;
          
           c
           i
          
         
         ,
          
         θ
        
        )
       
      
     
    
   
   
    
     Equation
     ⁢
        
     8
    
   
  
 



where the prior probability p(ci; θ) for each ci is set to




 
  
   1
   k
  
  .
 



In some embodiments, the distribution around each prototype is an isotropic Gaussian, which leads to:




 
  
   
    
     
      p
      ⁡
      (
      
       
        
         x
         i
        
        ;
        
         c
         i
        
       
       ,
       θ
      
      )
     
     =
     
      
       exp
       ⁡
       (
       
        
         -
         
          
           (
           
            
             v
             i
            
            -
            
             c
             s
            
           
           )
          
          2
         
        
        
         2
         ⁢
         
          σ
          s
          2
         
        
       
       )
      
      /
      
       
        ∑
        
         j
         =
         1
        
        k
       
        
       
        exp
        ⁡
        (
        
         
          -
          
           
            (
            
             
              v
              i
             
             -
             
              c
              j
             
            
            )
           
           2
          
         
         
          2
          ⁢
          
           σ
           j
           2
          
         
        
        )
       
      
     
    
   
   
    
     Equation
     ⁢
        
     9
    
   
  
 



where vi=ƒθ(x1) and xi∈cs. If the ProtoNCE module 215 applies the l2-normalization to both v and c, then (v−c)2=2−2v·c. Combining Equations 3, 4, 6, 7, 8, and 9, the maximum log-likelihood estimation may be written as:




 
  
   
    
     
      θ
      *
     
     =
     
      
       
        arg
        ⁢
          
        max
       
       θ
      
      ⁢
      
       
        ∑
        
         i
         =
         1
        
        N
       
       
        
         -
         log
        
        ⁢
        
         
          exp
          ⁡
          (
          
           
            v
            i
           
           ·
           
            
             c
             s
            
            
             u
             s
            
           
          
          )
         
         
          
           ∑
           
            j
            =
            1
           
           k
          
          
           exp
           ⁡
           (
           
            
             v
             i
            
            ·
            
             
              c
              j
             
             
              u
              j
             
            
           
           )
          
         
        
       
      
     
    
   
   
    
     Equation
     ⁢
        
     10
    
   
  
 



which is in the same form as the InfoNCE loss in Equation 1. Here μ∝σ2 denotes the concentration level of the feature distribution around a prototype (smaller μ means more concentrated distribution). Therefore, instance-wise contrastive learning can be interpreted as a special case of prototypical contrastive learning, where the prototypes are instance features (i.e. C=V′), and the concentration of the distribution around each instance is the same (i.e. μ=τ).

The ProtoNCE module 215 may sample r negative prototypes to calculate the normalization term. ProtoNCE module 215 may cluster samples M times with different number of clusters K={km}m=1M, which has a more robust probability estimation of prototypes that encode the hierarchical structure. Furthermore, a loss to instance discrimination may be added to retain the property of local smoothness. The ProtoNCE function used by the ProtoNCE module 215 to determine ProtoNCE loss may be defined as:



 
  
   
    
     
      L
      
       P
       ⁢
       r
       ⁢
       o
       ⁢
       t
       ⁢
       o
       ⁢
       N
       ⁢
       C
       ⁢
       E
      
     
     =
     
      
       ∑
       
        i
        =
        1
       
       N
      
      
       -
       
        (
        
         
          log
          ⁢
          
           
            exp
            ⁡
            (
            
             
              
               v
               i
              
              ·
              
               v
               i
               ′
              
             
             /
             τ
            
            )
           
           
            
             ∑
             
              j
              =
              0
             
             r
            
            
             exp
             ⁡
             (
             
              
               
                v
                i
               
               ·
               
                v
                j
                ′
               
              
              /
              τ
             
             )
            
           
          
         
         +
         
          
           1
           M
          
          ⁢
          
           
            ∑
            
             m
             =
             1
            
            M
           
            
           
            log
            ⁢
            
             
              exp
              ⁡
              (
              
               
                
                 v
                 i
                
                ·
                
                 c
                 s
                 m
                
               
               /
               
                μ
                s
                m
               
              
              )
             
             
              
               ∑
               
                j
                =
                0
               
               r
              
              
               exp
               ⁡
               (
               
                
                 
                  v
                  i
                 
                 ·
                 
                  c
                  j
                  m
                 
                
                /
                
                 μ
                 j
                 m
                
               
               )
              
             
            
           
          
         
        
        )
       
      
     
    
   
   
    
     Equation
     ⁢
        
     11
    
   
  
 



As illustrated in FIG. 2, the distribution of embeddings around each prototype c has different level of concentration. μ may denote the concentration estimation (concentration M), where a smaller μ indicates larger concentration. Here μ may be calculated using the momentum features {v′z}z=1Z that are within the same cluster as a prototype c. The desired μ should be small (high concentration) if (1) the average distance between v′z and c is small, and (2) the cluster contains more feature points (i.e. Z is large). Therefore, μ may be defined as:



 
  
   
    
     μ
     =
     
      
       
        ∑
        
         z
         =
         1
        
        Z
       
        
       
        
         
         
          
           v
           z
           ′
          
          -
          c
         
         
        
        2
       
      
      
       Z
       ⁢
         
       
        log
        ⁡
        (
        
         Z
         +
         α
        
        )
       
      
     
    
   
   
    
     Equation
     ⁢
        
     12
    
   
  
 



where α is a smooth parameter to ensure that small clusters do not have an overly large μ. Also μ may be normalized for each set of prototypes Cm such that they have a mean of τ.

In some embodiments, in the ProtoNCE loss (Equation 11), μsm acts as a scaling factor on the similarity between an embedding vi and its prototype csm. With the proposed the similarity for embeddings in a loose cluster (lager μ) are down-scaled, pulling them closer to the prototype. On the contrary, embeddings in a tight cluster (smaller μ) have an up-scaled similarity, thus less encouraged to approach the prototype. Therefore, representation learning with the ProtoNCE function yields more balanced clusters with similar concentration. This prevents a trivial solution where most embeddings collapse to a single cluster, a problem that could only be heuristically addressed by data-resampling in DeepCluster. FIG. 4 is a histogram 400 of a cluster size for prototypical contrastive learning for 50000 clusters with a fixed and estimated concertation Histogram 400 illustrates that using a different concertation μ for each prototype yields more balanced clusters with similar sizes, which leads to better representation learning.
In some embodiments, minimizing the proposed ProtoNCE loss may be considered as simultaneously maximizing the mutual information between V and all the prototypes {V′, C1, . . . , CM}. This leads to better representation learning, for two reasons. First, the encoder 210 may learn the shared information among prototypes and ignore the individual noise that exists in each prototype. The shared information is more likely to capture higher-level semantic knowledge. Second, when compared to instance features, prototypes have a larger mutual information (MI) with the class labels. Furthermore, training the encoder 210 using the ProtoNCE loss function may increase the MI between the instance features (or their assigned prototypes) and the ground-truth class labels for all images in a training dataset.
In some embodiments, the PCL framework 130 can provide more insights into the nature of the learned prototypes. The optimization in Equation 10 is similar to optimizing the cluster-assignment probability p(s; xi, θ) using the cross-entropy loss, where the prototypes c represent weights for a linear classifier. With k-means clustering, the linear classifier has a fixed set of weights as the mean vectors for the representations in each cluster,



 
  c
  =
  
   
    1
    Z
   
   ⁢
   
    
     ∑
     
      z
      =
      1
     
     Z
    
    
     
      v
      z
      ′
     
     .
    
   
  
 



A similar idea has been used for few-shot learning, where a non-parametric prototypical classifier performs better than a parametrized linear classifier.

In some embodiments, the PCL framework 130 may be trained using the ImageNet-1M dataset, which contains approximately 1.28 million images in 1000 classes. Momentum encoder 205 or encoder 210 may be a ResNet-50, whose last fully-connected layer outputs a 128-D and L2-normalized feature. PCL framework 130 may perform data augmentation on images 140 with random crop, random color jittering, random horizontal flip, and random grayscale conversion. The PCL framework 130 may use a SGD as an optimizer, with a weight decay of 0.0001, a momentum of 0.9, and a batch size of 256. The PCL framework 130 may train for 200 epochs, where the PCL framework 130 may warm-up the network in the first 20 epochs by only using the InfoNCE loss. The initial learning rate is 0.03 and may be multiplied by 0.1 at 120 and 160 epochs. In terms of the hyper-parameters, we set τ=0.1, α=10, and number of clusters K={25000, 50000, 100000}. We use the GPU k-means implementation in faiss which takes approximately 10 seconds. The clustering is performed every epoch, which introduces 219 1⁄3 computational overhead due to a forward pass through the dataset. The number of negatives for ProtoNCE module 213 is set as k=16000.
FIG. 5 is a diagram of an E-M algorithm for training the prototypical contrastive learning framework, according to some embodiments. In FIG. 5, algorithm 1 (E-M algorithm 220) includes pseudo-code for prototypical contrastive learning.
At line 1, algorithm 1 receives input which includes an encoder function ƒθ, the training dataset X which could be images 140 or other unstructured data, and a number of clusters K={km}m=1M.
At line 2, a momentum encoder 205 is initialized to 0, which may be the weights of encoder 210.
At line 3, a number of epoch are initialized using the MaxEpoch variable.
At line 4, the momentum features V′ from the training dataset X are generated using the momentum encoder 210.
At lines 4-8 the E-step 225 is performed, where the clustering module 235 clusters V′ features into km clusters, returns prototypes Cm (lines 6) and estimates the concentration distribution um around each prototype using Equation 12 (line 7).
At lines 9-14 and M-step 230 is performed. In the M-step 230, the images 140 (or other unstructured data) in the training dataset X may be loaded in minibatches and passed through encoder 210 and momentum encoder 205 at lines 10. The ProtoNCE module 215 determines a loss function using the features from the encoder 210 and momentum encoder 205 at line 11 and as shown in Equation 11. At line 12, the encoder 210 is trained using the loss function which updates the weights of the encoder 210. At line 13, the weights of the momentum encoder 205 are updated with the weights of the encoder 210.
FIG. 6 is a diagram that illustrates clustering in the E-step of the E-M algorithm, according to some embodiments. FIG. 6 illustrates three clusters. Cluster 602 is a fine-grained cluster that includes images of a horse and a man. Cluster 604 is a coarse-grained cluster that includes images of a horse. Clusters 602 and 604 overlap as both include images of a horse ad different levels of granularity. Cluster 606 is a coarse-grained cluster that includes images of an eagle. Notably, cluster 606 does not overlap with clusters 602 and 606 as the images do not have common features.
Some examples of computing devices, such as computing device 100 may include non-transitory, tangible, machine readable media that include executable code that when run by one or more processors (e.g., processor 110) may cause the one or more processors to perform the processes of method 300. Some common forms of machine readable media that may include the processes of method 300 are, for example, floppy disk, flexible disk, hard disk, magnetic tape, any other magnetic medium, CD-ROM, any other optical medium, punch cards, paper tape, any other physical medium with patterns of holes, RAM, PROM, EPROM, FLASH-EPROM, any other memory chip or cartridge, and/or any other medium from which a processor or computer is adapted to read.
This description and the accompanying drawings that illustrate inventive aspects, embodiments, implementations, or applications should not be taken as limiting. Various mechanical, compositional, structural, electrical, and operational changes may be made without departing from the spirit and scope of this description and the claims. In some instances, well-known circuits, structures, or techniques have not been shown or described in detail in order not to obscure the embodiments of this disclosure. Like numbers in two or more figures represent the same or similar elements.
In this description, specific details are set forth describing some embodiments consistent with the present disclosure. Numerous specific details are set forth in order to provide a thorough understanding of the embodiments. It will be apparent, however, to one skilled in the art that some embodiments may be practiced without some or all of these specific details. The specific embodiments disclosed herein are meant to be illustrative but not limiting. One skilled in the art may realize other elements that, although not specifically described here, are within the scope and the spirit of this disclosure. In addition, to avoid unnecessary repetition, one or more features shown and described in association with one embodiment may be incorporated into other embodiments unless specifically described otherwise or if the one or more features would make an embodiment non-functional.
This application is further described with respect to the attached document in Appendix I., entitled “Unsupervised Representation Learning with Contrastive Prototypes,” 14 pages, which is considered part of this disclosure and the entirety of which is incorporated by reference.
Although illustrative embodiments have been shown and described, a wide range of modification, change and substitution is contemplated in the foregoing disclosure and in some instances, some features of the embodiments may be employed without a corresponding use of other features. One of ordinary skill in the art would recognize many variations, alternatives, and modifications. Thus, the scope of the invention should be limited only by the following claims, and it is appropriate that the claims be construed broadly and in a manner consistent with the scope of the embodiments disclosed herein.

## Claims

1. A method for training a prototypical contrastive learning (PCL) framework to cluster images, comprising:
receiving images at a momentum encoder of the PCL framework;
determining, using the momentum encoder, features from the images;
clustering, using the momentum encoder, the images into clusters according to the features;
determining at least one prototype and at least one concentration of the clusters;
determining, a contrastive loss function from the at least one prototype and the at least one concentration of the clusters;
training an encoder of the PCL framework using the contrastive loss function and a subset of images in the images, wherein the training updates the weights of the encoder based on the contrastive loss function;
updating weights of the momentum encoder using weights of the encoder; and
determining second clusters from second images using the updated momentum encoder.

2. The method of claim 1, wherein the encoder is a convolutional neural network.

3. The method of claim 1, wherein a neural network structure of the momentum encoder is the same as a neural network structure of the encoder.

4. The method of claim 1, wherein the updated weights of the momentum encoder are a moving average of the weights of the momentum encoder and the weights of the encoder.

5. The method of claim 1, further comprising:
dividing the images into multiple subsets of images; and
training the encoder of the PCL framework on the subsets of images until the contrastive loss function is minimized.

6. The method of claim 1, wherein the updated momentum encoder is trained to perform at least one task associated with processing the images.

7. The method of claim 1, further comprising:
training the encoder and updating the momentum encoder over a configurable number of iterations.

8. A system for training a prototypical contrastive learning (PCL) framework to cluster images, comprising:
at least one memory configured to store the PCL framework and images; and
a processor coupled to the memory and configured to execute instructions for training the PCL framework, the instructions comprising:
receiving the images at a momentum encoder of the PCL framework;
determining, using the momentum encoder, features from the images;
clustering, using the momentum encoder, the images into clusters according to the features;
determining at least one prototype and at least one concentration of the clusters;
determining, a contrastive loss function from the at least one prototype and the at least one concentration of the clusters;
training an encoder of the PCL framework using the contrastive loss function and a subset of images in the images;
updating weights of the momentum encoder using weights of the encoder; and
determining second clusters from second images using the updated momentum encoder.

9. The system of claim 8, wherein the encoder is a convolutional neural network.

10. The system of claim 8, wherein a neural network structure of the momentum encoder is the same as a neural network structure of the encoder.

11. The system of claim 8, wherein the training updates the weights of the encoder based on the contrastive loss function.

12. The system of claim 8, wherein the updated weights of the momentum encoder are a moving average of the weights of the momentum encoder and the weights of the encoder.

13. The system of claim 8, wherein the processor is further configured to perform instructions comprising:
dividing the images into multiple subsets of images; and
training the encoder of the PCL framework on the subsets of images until the contrastive loss function is minimized.

14. The system of claim 8, wherein the updated momentum encoder is trained to perform at least one task associated with processing the images.

15. The system of claim 8, wherein the processor is further configured to perform instructions comprising:
training the encoder and updating the momentum encoder over a configurable number of iterations.

16. A non-transitory computer-readable medium storing instructions thereon, that when executed by a processor, cause the processor to perform operations that train a prototypical contrastive learning (PCL) framework to cluster images, the operations comprising:
receiving images at a momentum encoder of the PCL framework;
determining, using the momentum encoder, features from the images;
clustering, using the momentum encoder, the images into clusters according to the features;
determining at least one prototype and at least one concentration of the clusters;
determining, a contrastive loss function from the at least one prototype and the at least one concentration of the clusters;
training an encoder of the PCL framework using the contrastive loss function and a subset of images in the images, wherein the training updates the weights of the encoder based on the contrastive loss function;
updating weights of the momentum encoder using weights of the encoder; and
determining second clusters from second images using the updated momentum encoder.

17. The non-transitory computer-readable medium of claim 16, wherein a neural network structure of the momentum encoder is the same as a neural network structure of the encoder.

18. The non-transitory computer-readable medium of claim 16, wherein the updated weights of the momentum encoder are a moving average of the weights of the momentum encoder and the weights of the encoder.

19. The non-transitory computer-readable medium of claim 16, wherein the operations further comprise:
dividing the images into multiple subsets of images; and
training the encoder of the PCL framework on the subsets of images until the contrastive loss function is minimized.

20. The non-transitory computer-readable medium of claim 16, wherein the updated momentum encoder is trained to perform at least one task associated with processing the images.