One-shot learning (computer vision)

Background

As with most classification schemes, one-shot learning involves three main challenges:

• Representation: How should objects and categories be described?
• Learning: How can such descriptions be created?
• Recognition: How can a known object be filtered from enveloping clutter, irrespective of occlusion, viewpoint, and lighting?

One-shot learning differs from single object recognition and standard category recognition algorithms in its emphasis on knowledge transfer, which makes use of previously learned categories.

• Model parameters: Reuses model parameters, based on the similarity between old and new categories. Categories are first learned on numerous training examples, then new categories are learned using transformations of model parameters from those initial categories or selecting relevant parameters for a classifier.
• Feature sharing: Shares parts or features of objects across categories. One algorithm extracts "diagnostic information" in patches from already learned categories by maximizing the patches' mutual information, and then applies these features to the learning of a new category. A dog category, for example, may be learned in one shot from previous knowledge of horse and cow categories, because dog objects may contain similar distinguishing patches.
• Contextual information: Appeals to global knowledge of the scene in which the object appears. Such global information can be used as frequency distributions in a conditional random field framework to recognize objects. Alternatively context can consider camera height and scene geometry. Algorithms of this type have two advantages. First, they learn object categories that are relatively dissimilar; and second, they perform well in ad hoc situations where an image has not been hand-cropped and aligned.

Theory

The Bayesian one-shot learning algorithm represents the foreground and background of images as parametrized by a mixture of constellation models. During the learning phase, the parameters of these models are learned using a conjugate density parameter posterior and variational Bayesian expectation–maximization (VBEM). In this stage the previously learned object categories inform the choice of model parameters via transfer by contextual information. For object recognition on new images, the posterior obtained during the learning phase is used in a Bayesian decision framework to estimate the ratio of p(object | test, train) to p(background clutter | test, train) where p is the probability of the outcome.

Bayesian framework

Given the task of finding a particular object in a query image, the overall objective of the Bayesian one-shot learning algorithm is to compare the probability that object is present vs the probability that only background clutter is present. If the former probability is higher, the algorithm reports the object's presence, otherwise the algorithm reports its absence. To compute these probabilities, the object class must be modeled from a set of (1 ~ 5) training images containing examples.

To formalize these ideas, let ${\displaystyle I}$ be the query image, which contains either an example of the foreground category ${\displaystyle O_{fg}}$ or only background clutter of a generic background category ${\displaystyle O_{bg}}$. Also let ${\displaystyle I_t}$ be the set of training images used as the foreground category. The decision of whether ${\displaystyle I}$ contains an object from the foreground category, or only clutter from the background category is:

${\displaystyle R=\frac{p(O_{fg}|I,I_t)}{p(O_{bg}|I,I_t)}=\frac{p(I|I_t,O_{fg})p(O_{fg})}{p(I|I_t,O_{bg})p(O_{bg})},}$

where the class posteriors ${\displaystyle p(O_{fg}|I,I_t)}$ and ${\displaystyle p(O_{bg}|I,I_t)}$ have been expanded by Bayes' theorem, yielding a ratio of likelihoods and a ratio of object category priors. We decide that the image ${\displaystyle I}$ contains an object from the foreground class if ${\displaystyle R}$ exceeds a certain threshold ${\displaystyle T}$. We next introduce parametric models for the foreground and background categories with parameters ${\displaystyle \theta }$ and ${\displaystyle {\theta }_{bg}}$ respectively. This foreground parametric model is learned during the learning stage from ${\displaystyle I_t}$, as well as prior information of learned categories. The background model we assume to be uniform across images. Omitting the constant ratio of category priors, ${\displaystyle \frac{p(O_{fg})}{p(O_{bg})}}$, and parametrizing over ${\displaystyle \theta }$ and ${\displaystyle {\theta }_{bg}}$ yields

${\displaystyle R\propto \frac{\int p(I|\theta ,O_{fg})p(\theta |I_t,O_{fg})d\theta }{\int p(I|{\theta }_{bg},O_{bg})p({\theta }_{bg}|I_t,O_{bg})d{\theta }_{bg}}=\frac{\int p(I|\theta )p(\theta |I_t,O_{fg})d\theta }{\int p(I|{\theta }_{bg})p({\theta }_{bg}|I_t,O_{bg})d{\theta }_{bg}}}$, having simplified ${\displaystyle p(I|\theta ,O_{fg})}$ and ${\displaystyle p(I|\theta ,O_{bg})}$ to ${\displaystyle p(I|{\theta }_{fg})}$ and ${\displaystyle p(I|{\theta }_{bg}).}$

The posterior distribution of model parameters given the training images, ${\displaystyle p(\theta |I_t,O_{fg})}$ is estimated in the learning phase. In this estimation, one-shot learning differs sharply from more traditional Bayesian estimation models that approximate the integral as ${\displaystyle \delta ({\theta }^{ML})}$. Instead, it uses a variational approach using prior information from previously learned categories. However, the traditional maximum likelihood estimation of the model parameters is used for the background model and the categories learned in advance through training.

Object category model

For each query image ${\displaystyle I}$ and training images ${\displaystyle I_t}$, a constellation model is used for representation. To obtain this model for a given image ${\displaystyle I}$, first a set of N interesting regions is detected in the image using the Kadir–Brady saliency detector. Each region selected is represented by a location in the image, ${\displaystyle X_i}$ and a description of its appearance, ${\displaystyle A_i}$. Letting ${\displaystyle X=\sum _{i=1}^N{X}_i,A=\sum _{i=1}^N{A}_i}$ and ${\displaystyle X_t}$and ${\displaystyle A_t}$ the analogous representations for training images, the expression for R becomes:

${\displaystyle R\propto \frac{\int p(X,A|\theta ,O_{fg})p(\theta |X_t,A_t,O_{fg})d\theta }{\int p(X,A|{\theta }_{bg},O_{bg})p({\theta }_{bg}|X_t,A_t,O_{bg})d{\theta }_{bg}}=\frac{\int p(X,A|\theta )p(\theta |X_t,A_t,O_{fg})d\theta }{\int p(X,A|{\theta }_{bg})p({\theta }_{bg}|X_t,A_t,O_{bg})d{\theta }_{bg}}}$

The likelihoods ${\displaystyle p(X,A|\theta )}$ and ${\displaystyle p(X,A|{\theta }_{bg})}$ are represented as mixtures of constellation models. A typical constellation model has P(3 ~ 7) parts, with N(~100) interest regions. Thus a P-dimensional vector h assigns one region of interest (out of N regions) to each model part (for P parts). Thus h denotes a hypothesis (an assignment of interest regions to model parts) for the model and a full constellation model is represented by summing over all possible hypotheses h in the hypothesis space ${\displaystyle H}$. Finally the likelihood is written

${\displaystyle p(X,A|\theta )=\sum _{\omega =1}^{\mathrm{\Omega }}\sum _{\mathbf{\text{h}}\in H}p(X,A,\mathbf{\text{h}},\omega |\theta ).}$

The different ${\displaystyle \omega }$'s represent different configurations of parts, whereas the different hypotheses h represent different assignations of regions to parts, given a part model ${\displaystyle \omega }$. The assumption that the shape of the model (as represented by ${\displaystyle X}$, the collection of part locations) and appearance are independent allows one to consider the likelihood expression ${\displaystyle p(X,A,\mathbf{\text{h}},\omega |\theta )}$ as two separate likelihoods of appearance and shape.

Appearance

The appearance of each feature is represented by a point in appearance space (discussed below in implementation). "Each part ${\displaystyle p}$ in the constellation model has a Gaussian density within this space with mean and precision parameters ${\displaystyle {\theta }_{p,\omega }^A={\mu }_{p,\omega }^A,{\mathrm{\Gamma }}_{p,\omega }^A}$." From these the appearance likelihood described above is computed as a product of Gaussians over the model parts for a give hypothesis h and mixture component ${\displaystyle \omega }$.

Shape

The shape of the model for a given mixture component ${\displaystyle \omega }$ and hypothesis h is represented as a joint Gaussian density of the locations of features. These features are transformed into a scale and translation-invariant space before modelling the relative location of the parts by a 2(P - 1)-dimensional Gaussian. From this, we obtain the shape likelihood, completing our representation of ${\displaystyle p(X,A,\mathbf{\text{h}},\omega |\theta )}$. To reduce the number of hypotheses in the hypothesis space ${\displaystyle H}$, only those hypotheses that satisfy the ordering constraint that the x-coordinate of each part is monotonically increasing are considered. This eliminates ${\displaystyle P!}$ hypotheses from ${\displaystyle H}$.

Conjugate densities

To compute ${\displaystyle R}$, the integral ${\displaystyle \int p(X,A|\theta )p(\theta |X_t,A_t,O_{fg})d\theta }$ must be evaluated, but is analytically intractable. The object category model above gives information about ${\displaystyle p(X,A|\theta )}$, so what remains is to examine ${\displaystyle p(\theta |X_t,A_t,O)}$, the posterior of ${\displaystyle \theta }$, and find a sufficient approximation to render the integral tractable. Previous work approximates the posterior by a ${\displaystyle \delta }$function centered at ${\displaystyle {\theta }^{\ast }}$, collapsing the integral in question into ${\displaystyle p(X,A|{\theta }^{\ast })}$. This ${\displaystyle {\theta }^{\ast }}$ is normally estimated using a maximum likelihood (${\displaystyle {\theta }^{\ast }={\theta }^{ML}}$) or maximum a posteriori (${\displaystyle {\theta }^{\ast }={\theta }^{MAP}}$) procedure. However, because in one-shot learning, few training examples are used, the distribution will not be well-peaked, as is assumed in a ${\displaystyle \delta }$function approximation. Thus instead of this traditional approximation, the Bayesian one-shot learning algorithm seeks to "find a parametric form of ${\displaystyle p(\theta )}$ such that the learning of ${\displaystyle p(\theta |X_t,A_t,O_{fg})}$ is feasible". The algorithm employs a Normal-Wishart distribution as the conjugate prior of ${\displaystyle p(\theta |X_t,A_t,O_{fg})}$, and in the learning phase, variational Bayesian methods with the same computational complexity as maximum likelihood methods are used to learn the hyperparameters of the distribution. Then, since ${\displaystyle p(X,A|\theta )}$ is a product of Gaussians, as chosen in the object category model, the integral reduces to a multivariate Student's T distribution, which can be evaluated.

Implementation

Feature detection and representation

To detect features in an image so that it can be represented by a constellation model, the Kadir–Brady saliency detector is used on grey-scale images, finding salient regions of the image. These regions are then clustered, yielding a number of features (the clusters) and the shape parameter ${\displaystyle X}$, composed of the cluster centers. The Kadir–Brady detector was chosen because it produces fewer, more salient regions, as opposed to feature detectors like multiscale Harris, which produces numerous, less significant regions.

The regions are then taken from the image and rescaled to a small patch of 11 × 11 pixels, allowing each patch to be represented in 121-dimensional space. This dimensionality is reduced using principal component analysis, and ${\displaystyle A}$, the appearance parameter, is then formed from the first 10 principal components of each patch.

Learning

To obtain shape and appearance priors, three categories (spotted cats, faces, and airplanes) are learned using maximum likelihood estimation. These object category model parameters are then used to estimate the hyper-parameters of the desired priors.

Given a set of training examples, the algorithm runs the feature detector on these images, and determines model parameters from the salient regions. The hypothesis index h assigning features to parts prevents a closed-form solution of the linear model, so the posterior ${\displaystyle p(\theta |X_t,A_t,O_{fg})}$ is estimated by variational Bayesian expectation–maximization algorithm, which is run until parameter convergence after ~ 100 iterations. Learning a category in this fashion takes under a minute on a 2.8 GHz machine with a 4-part model and < 10 training images.

Experimental results

Motorbike example

To learn the motorbike category:

• Six training images are selected from the motorbike category of the Caltech 4 Data Set and the Kadir–Brady detector is applied, giving ${\displaystyle X_t}$ and through PCA, ${\displaystyle A_t}$.
• Next, the prior model parameters are computed from 30 models ${\displaystyle {\theta }_t}$, 10 from each of the three learned categories: spotted cats, faces, and airplanes. This prior encodes the knowledge that "models lacking visual consistency [ie background clutter] occupy a different part of the parameter space [from] coherent models."
• In learning, which is performed next, the prior biases the posterior ${\displaystyle p(\theta |X_t,A_t,O_{fg})}$ towards parts of the parameter space corresponding to coherent models. Only one mixture component is used, letting ${\displaystyle \mathrm{\Omega }=1}$. The estimation of the posterior is shown below.
• Finally, the figures below show the learned motorbike model with shape and appearance of parts, and the corresponding features.
• For recognition tests, the model above is applied to 50 images which contain motorbikes and 50 which do not. The image below shows an ROC curve, measuring the probability of detection over the probability of false detection, as well as some recognized examples.

Shared densities on transforms

Another algorithm uses knowledge transfer by model parameters to learn a new object category that is similar in appearance to previously learned categories. An image is represented as either a texture and shape, or as a latent image that has been transformed, denoted by ${\displaystyle I=T(I_L)}$.

A Siamese neural network works in tandem on two different input vectors to compute comparable output vectors.

Congealing

In this context, congealing is "the simultaneous vectorization of each of a set of images to each other". For a set of training images of a certain category, congealing iteratively transforms each image to minimize the images' joint pixelwise entropies E, where

${\displaystyle E=\sum _{p=1}^{P}H(\nu (p)),}$

"where ${\displaystyle \nu (p)}$ is the binary random variable defined by the values of a particular pixel p across all of the images, ${\displaystyle H()}$ is the discrete entropy function of that variable, and ${\displaystyle 1\le p\le P}$ is the set of pixel indices for the image."

The congealing algorithm begins with a set of images ${\displaystyle I_i}$ and a corresponding transform matrix ${\displaystyle U_i}$, which at the end of the algorithm will represent the transformation of ${\displaystyle I_i}$ into its latent ${\displaystyle I_{L_i}}$. These latents ${\displaystyle I_{L_i}}$ minimize the joint pixel-wise entropies. Thus the task of the congealing algorithm is to estimate the transformations ${\displaystyle U_i}$.

Sketch of algorithm:

• Initialize ${\displaystyle U_I}$'s to the identity.
• Compute the joint pixelwise entropies of the current set of images.
• For each image ${\displaystyle I_i}$, iterate through all possible affine transformations ${\displaystyle A}$ (rotation, x-translation, y-translation, x-scale, y-scale, x-shear, y-shear) and test if ${\displaystyle A{U}_i}$ decreases the joint pixelwise entropies. If so, set ${\displaystyle U_i=A{U}_i}$.
• Repeat previous step until convergence.

At the end of the algorithm, ${\displaystyle U_i(I)=I_{L_i}}$, and ${\displaystyle T=U_i^{-1}}$ transforms the latent image back into the originally observed image.

Classification

To use this model for classification, it must be estimated with the maximum posterior probability given an observed image ${\displaystyle I}$. Applying Bayes' rule to ${\displaystyle P(c_j|I)}$ and parametrization by the transformation ${\displaystyle T}$ gives a difficult integral that must be approximated, and then the best transform ${\displaystyle T}$ (that which maps the test image to its latent image) must be found. Once this transformation is found, the test image can be transformed into its latent, and a nearest neighbor classifier based on Hausdorff distance between images can classify the latent (and thus the test image) as belonging to a particular class ${\displaystyle c_j}$.

To find ${\displaystyle T}$, the test image I is inserted into the training ensemble for the congealing process. Since the test image is drawn from one of the categories ${\displaystyle c_j}$, congealing provides a corresponding ${\displaystyle T_{\text{test}}=U_{\text{test}}^{-1}}$ that maps I to its latent. The latent can then be classified.

Single-example classification

Given a set of transformations ${\displaystyle B_i}$ obtained from congealing many images of a certain category, the classifier can be extended to the case where only one training ${\displaystyle I_t}$ example of a new category ${\displaystyle c}$ is allowed. Applying all the transformations ${\displaystyle B_i}$ sequentially to ${\displaystyle I_t}$ creates an artificial training set for ${\displaystyle c}$. This artificial data set can be made larger by borrowing transformations from many already known categories. Once this data set is obtained, ${\displaystyle I}$, a test instance of ${\displaystyle c}$, can be classified as in the normal classification procedure. The key assumption is that categories are similar enough that the transforms from one can be applied to another.