View source for One-shot learning (computer vision)

{{Short description|Object categorization problem}}
{{research paper|date=April 2016}}
'''One-shot learning''' is an [[object categorization problem]], found mostly in [[computer vision]]. Whereas most [[machine learning]]-based object categorization algorithms require training on hundreds or thousands of examples, one-shot learning aims to classify objects from one, or only a few, examples. The term '''few-shot learning''' is also used for these problems, especially when more than one example is needed.

{{TOC limit|3}}

==Motivation==
The ability to learn object categories from few examples, and at a rapid pace, has been demonstrated in humans.{{sfn|Li |Fergus|Perona|2002}}{{sfn|Thorpe|Fize|Marlot|1996}} It is estimated that a child learns almost all of the 10 ~ 30 thousand object categories in the world by age six.{{sfn|Biederman|1987}} This is due not only to the human mind's computational power, but also to its ability to synthesize and learn new object categories from existing information about different, previously learned categories. Given two examples from two object categories: one, an unknown object composed of familiar shapes, the second, an unknown, amorphous shape; it is much easier for humans to recognize the former than the latter, suggesting that humans make use of previously learned categories when learning new ones. The key motivation for solving one-shot learning is that systems, like humans, can use knowledge about object categories to classify new objects.{{sfn|Li |Fergus|Perona|2006|loc=Section 1}}{{sfn|Li|2006|loc=Section 1}}

==Background==
As with most [[classifier system|classification scheme]]s, 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?{{sfn|Li|Fergus|Perona|2006|loc=Section 2}}

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.{{sfn|Fink|2004}}
* 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.{{sfn|Bart|Ullman|2005}}
* 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.{{sfn|Murphy|et al|2004}} Alternatively context can consider camera height and scene geometry.{{sfn|Hoiem|Efros|Herbert|2005}} 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.{{sfn|Li|2006|loc=Section 2}}

==Theory==
The [[Bayesian statistics|Bayesian]] one-shot learning algorithm represents the foreground and background of images as parametrized by a mixture of constellation models.{{sfn|Burl |et al| 1996}} During the learning phase, the parameters of these models are learned using a [[conjugate prior|conjugate]] density parameter [[posterior probability|posterior]] and Variational Bayesian [[expectation–maximization algorithm|Expectation–Maximization]] (VBEM).{{sfn|Attias|1999}} 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.{{sfn|Li |et al|2006}}

===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 <math> I</math> be the query image, which contains either an example of the foreground category <math> O_{fg} </math> or only background clutter of a generic background category <math> O_{bg} </math>. Also let <math> I_t </math> be the set of training images used as the foreground category. The decision of whether <math> I </math> contains an object from the foreground category, or only clutter from the background category is:

: <math> 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})}, </math>

where the class posteriors <math> p(O_{fg} |I, I_t) </math> and <math>p(O_{bg}|I, I_t) </math> have been expanded by [[Bayes' Theorem]], yielding a ratio of [[likelihood function|likelihoods]] and a ratio of object category [[prior distribution|priors]]. We decide that the image <math>I</math> contains an object from the foreground class if <math>R</math> exceeds a certain threshold <math>T</math>. We next introduce parametric models for the foreground and background categories with parameters <math> \theta </math> and <math> \theta_{bg} </math> respectively. This foreground parametric model is learned during the learning stage from <math> I_t </math>, 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, <math> \frac{p(O_{fg})}{p(O_{bg})} </math>, and parametrizing over <math> \theta </math> and <math> \theta_{bg} </math> yields

: <math> 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}}</math>, having simplified <math>p(I | \theta, O_{fg}) </math> and <math>p(I | \theta, O_{bg})</math> to <math>p(I | \theta_{fg}) </math> and <math>p(I | \theta_{bg}).</math>

The posterior distribution of model parameters given the training images, <math> p(\theta | I_t, O_{fg}) </math> 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 <math> \delta(\theta^{ML}) </math>. 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.{{sfn|Li|Fergus|Perona|2006|loc=Section 3.1}}

===Object category model===
For each query image <math>I</math> and training images <math> I_t</math>, a [[constellation model]] is used for representation.{{sfn|Burl |et al| 1996}}{{sfn|Weber|Welling|Perona|2000}}{{sfn|Fergus|Perona |Zisserman|2003}} To obtain this model for a given image <math>I</math>, first a set of N interesting regions is detected in the image using the [[Kadir–Brady saliency detector]].{{sfn|Kadir |Brady| 2001}} Each region selected is represented by a location in the image, <math>X_i</math> and a description of its appearance, <math> A_i </math>. Letting <math> X = \sum_{i = 1}^N X_i, A = \sum_{i = 1}^N A_i </math> and <math>X_t </math>and <math>A_t</math> the analogous representations for training images, the expression for R becomes:

: <math> 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}}</math>

The likelihoods <math> p(X, A|\theta) </math> and <math> p(X, A|\theta_{bg}) </math> are represented as [[mixture model|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 <math>H</math>. Finally the likelihood is written

: <math>p(X,A|\theta) = \sum_{\omega=1}^{\Omega} \sum_{\textbf{h} \in H} p(X,A,\textbf{h}, \omega | \theta). </math>

The different <math>\omega</math>'s represent different configurations of parts, whereas the different hypotheses '''h''' represent different assignations of regions to parts, given a part model <math>\omega</math>. The assumption that the shape of the model (as represented by <math>X</math>, the collection of part locations) and appearance are independent allows one to consider the likelihood expression <math>p(X,A,\textbf{h}, \omega | \theta) </math> as two separate likelihoods of appearance and shape.{{sfn|Li|Fergus|Perona|2006|loc=Section 3.2}}

====Appearance====
The appearance of each feature is represented by a point in appearance space (discussed below in implementation). "Each part <math> p </math> in the constellation model has a [[Normal distribution|Gaussian density]] within this space with mean and precision parameters <math> \theta_{p,\omega}^{A} = { \mu_{p,\omega}^{A}, \Gamma_{p,\omega}^{A} } </math>." 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 <math> \omega </math>.{{sfn|Li|Fergus|Perona|2006|loc=Section 3.2.1}}

====Shape====
The shape of the model for a given mixture component <math> \omega </math> 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 <math>p(X,A, \textbf{h}, \omega | \theta)</math> . In order to reduce the number of hypotheses in the hypothesis space <math> H</math>, only those hypotheses that satisfy the ordering constraint that the x-coordinate of each part is monotonically increasing are considered. This eliminates <math>P! </math> hypotheses from <math> H </math>.{{sfn|Li|Fergus|Perona|2006|loc=Section 3.2.1}}

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

==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 <math> X </math>, 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 <math> A </math>, the appearance parameter, is then formed from the first 10 principal components of each patch.{{sfn|Li|Fergus|Perona|2006|loc=Section 5.1}}

===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 <math> p(\theta| X_t, A_t, O_{fg})</math> 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&nbsp;GHz machine with a 4-part model and < 10 training images.{{sfn|Li|Fergus|Perona|2006|loc=Sections 4, 5.2}}

==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 <math> X_t </math> and through [[principal component analysis|PCA]], <math> A_t </math>.
* Next, the prior model parameters are computed from 30 models <math> \theta_t </math>, 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 <math> p(\theta| X_t, A_t, O_{fg})</math> towards parts of the parameter space corresponding to coherent models. Only one mixture component is used, letting <math> \Omega = 1 </math>. 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 <math> I = T(I_L) </math>.

A [[Siamese neural network]] works in tandem on two different input vectors to compute comparable output vectors.<ref>{{cite AV media| url-status = live| archive-url = https://ghostarchive.org/varchive/youtube/20211210/4S-XDefSjTM| archive-date = 2021-12-10| url = https://www.youtube.com/watch?v=4S-XDefSjTM| title = Few-Shot Learning (2/3): Siamese Networks | website=[[YouTube]]}}{{cbignore}}</ref>

===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

: <math> E = \sum_{p=1}^{P}H(\nu(p)), </math>

"where <math>\nu(p)</math> is the binary random variable defined by the values of a particular pixel p across all of the images, <math>H( )</math> is the discrete entropy function of that variable, and <math> 1\leq p \leq P </math> is the set of pixel indices for the image."

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

Sketch of algorithm:

* Initialize <math>U_I</math>'s to the identity.
* Compute the joint pixelwise entropies of the current set of images.
* For each image <math>I_i</math>, iterate through all possible affine transformations <math>A</math> (rotation, x-translation, y-translation, x-scale, y-scale, x-shear, y-shear) and test if <math>AU_i</math> decreases the joint pixelwise entropies. If so, set <math>U_i = AU_i</math>.
* Repeat previous step until convergence.

At the end of the algorithm, <math>U_i(I) = I_{L_i}</math>, and <math>T = U_i^{-1}</math> transforms the latent image back into the originally observed image.{{sfn|Miller|Matsakis|Viola|2000|Section 3}}

===Classification===
To use this model for classification, it must be estimated with the maximum posterior probability given an observed image <math>I</math>. Applying Bayes' rule to <math>P(c_j|I)</math> and parametrization by the transformation <math>T</math> gives a difficult integral that must be approximated, and then the best transform <math>T</math> (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 algorithm|nearest neighbor classifier]] based on [[Hausdorff distance]] between images can classify the latent (and thus the test image) as belonging to a particular class <math>c_j</math>.

To find <math>T</math>, 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 <math>c_j</math>, congealing provides a corresponding <math>T_{\text{test}} = U_{\text{test}}^{-1}</math> that maps I to its latent. The latent can then be classified.{{sfn|Miller|Matsakis|Viola|2000|loc=Section 4}}

===Single-example classification===
Given a set of transformations <math>B_i</math> obtained from congealing many images of a certain category, the classifier can be extended to the case where only one training <math>I_t</math> example of a new category <math>c</math> is allowed. Applying all the transformations <math>B_i</math> sequentially to <math>I_t</math> creates an artificial training set for <math>c</math>. This artificial data set can be made larger by borrowing transformations from many already known categories. Once this data set is obtained, <math>I</math>, a test instance of <math>c</math>, 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.{{sfn|Miller|Matsakis|Viola|2000|loc=Section 7}}

==See also==
* [[Variational Bayesian methods]]
* [[Variational message passing]]
* [[Expectation–maximization algorithm]]
* [[Bayesian inference]]
* [[Feature detection (computer vision)|Feature detection]]
* [[Association rule learning]]
* [[Hopfield network]]
* [[Zero-shot learning]]

==Citations==
{{reflist}}

==References==
*{{cite journal|first=Fei Fei|last= Li|title=Knowledge transfer in learning to recognize visual object classes|journal=International Conference on Development and Learning (ICDL)|year=2006|url=http://vision.stanford.edu/documents/Li_ICDL2006.pdf}}
*{{cite journal|first1=Fei Fei|last1=Li |first2=R. |last2=Fergus |first3=P. |last3=Perona|title=One-Shot learning of object categories|journal=IEEE Transactions on Pattern Analysis and Machine Intelligence |volume=28|issue=4 |pages=594–611|year=2006 |doi=10.1109/TPAMI.2006.79 |pmid=16566508 |s2cid=6953475 |url=http://vision.stanford.edu/documents/LiFergusPerona2006.pdf}}
*{{cite journal|last1=Miller|last2=Matsakis|last3=Viola|url=https://people.cs.umass.edu/~elm/papers/cvpr2000.pdf|title=Learning from One Example through Shared Densities on Transforms|journal=Proc. Computer Vision and Pattern Recognition|year=2000}}
*{{cite journal|first1=F.F. |last1=Li |first2=R. |last2=VanRullen |first3=C.|last3=Coch|first4= P. |last4=Perona|title=Rapid natural scene categorization in the near absence of attention|journal=PNAS|volume=99 |issue=14|pages=9596–9601|year= 2002|doi=10.1073/pnas.092277599 |pmid=12077298 |pmc=123186 |bibcode=2002PNAS...99.9596L |doi-access=free }}
*{{cite journal|first1=S. |last1=Thorpe |first2=D. |last2=Fize|first3= C.|last3= Marlot|url=https://www.allpsych.uni-giessen.de/rauisch12/readings/ThorpeEtal.Nature.1996.pdf|title=Speed of processing in the human visual system|journal=Nature|year=1996 |volume=381|issue=6582 |pages=520–522, 1996|doi=10.1038/381520a0 |pmid=8632824 |bibcode=1996Natur.381..520T |s2cid=4303570 }}
*{{cite journal|first=I. |last=Biederman |url=http://courses.nus.edu.sg/course/eleclf/ee6901/Recognition%20by%20Components.pdf|title=Recognition-by-Components: a theory of human understanding|journal=Psychological Review|volume=94|pages=115–147|year=1987|issue=2 |doi=10.1037/0033-295X.94.2.115 |pmid=3575582 }}
*{{cite journal|first=M. |last=Fink|title=Object classification from a single example utilizing class relevance pseudo-metrics|journal=NIPS|year=2004|citeseerx=10.1.1.91.7461}}
*{{cite journal|last1=Bart |last2=Ullman |url=https://www.inf.tu-dresden.de/content/institutes/ki/is/HS_SS08_Papers/BartUllmanCVPR05.pdf|title=Cross-generalization: learning novel classes from a single example by feature replacement|journal=CVPR|year=2005}}
*{{cite journal|first1=K. |last1=Murphy|first2=A. |last2=Torralba|first3= W.T. |last3=Freeman|url=https://www.cs.ubc.ca/~murphyk/Papers/nips2003_long.pdf|title=Using the forest to see the trees: a graphical model relating features, objects, and scenes|journal=NIPS|year=2004}}
*{{cite journal|first1=D. |last1=Hoiem|first2= A.A. |last2=Efros|first3= M. 
|last3=Herbert|url=https://www.cs.utexas.edu/~grauman/courses/spring2009/slides/luis_pres.pdf|title=Geometric context from a single image|journal=ICCV|year=2005}}
*{{cite journal|first=H. |last=Attias|title=Inferring Parameters and Structure of Latent Variable Models by Variational Bayes|journal=Proc. Of the 15th Conf. In Uncertainty in Artificial Intelligence|pages=21–30|year=1999|arxiv=1301.6676}}
*{{cite journal|first1=M. |last1=Burl|first2= M. |last2=Weber|first3= P. |last3=Perona|url=https://link.springer.com/content/pdf/10.1007/BFb0054769.pdf|title=A Probabilistic Approach to Object Recognition Using Local Photometry and Global Geometry|journal=Proc. European Conf. Computer Vision|series=Lecture Notes in Computer Science|pages=628–641|year=1996|volume=1407|doi=10.1007/BFb0054769|isbn=978-3-540-64613-6}}
*{{cite journal|first1=R. |last1=Fergus |first2=P. |last2=Perona |first3= A. |last3=Zisserman|url=http://host.robots.ox.ac.uk/pascal/VOC/pubs/fergus03.pdf|title=Object Class Recognition by Unsupervised Scale-Invariant Learning|journal=Proc. Computer Vision and Pattern Recognition|pages=264–271|year= 2003}}
*{{cite journal|first1=M. |last1=Weber |first2=M. |last2=Welling|first3= P. |last3=Perona|url=https://link.springer.com/content/pdf/10.1007/3-540-45054-8_2.pdf|title=Unsupervised Learning of Models for Recognition|journal=Proc. European Conf. Computer Vision|series=Lecture Notes in Computer Science |pages=101–108|year= 2000|volume=1842 |doi=10.1007/3-540-45054-8_2 |isbn=978-3-540-67685-0 }}
*{{cite journal|first1=T.|last1=Kadir|first2=M.|last2=Brady|title=Scale, Saliency, and Image Description|url=https://www.researchgate.net/publication/220660282|journal=International Journal of Computer Vision|volume=45|issue=2|pages=83–105|year=2001|doi=10.1023/A:1012460413855|s2cid=825395}}

[[Category:Learning in computer vision]]