hamming Distance Metric Learning

Introduction

This paper tries to propose a method to learn mappings from high dimensional data to binary codes. One of the main advantages of using binary space is that one can do exact KNN classification in sublinear time. Like other metric learning method this paper also tries to optimize some cost function which is based one a similarity measure between data points. One choice of similarity measure in binary space is Euclidean distance which produces unsatisfactory results. Another choice is Hamming distance, which is the total number of positions at which the corresponding bits are different.

The task is to learn a mapping from b(x) that project p-dimensional real valued input x onto a q dimensional binary code while preserving some notion of similarity. This mapping, which is called hash function is parameterized by a matrix w such that:
$b(x,w)=sign(f(w,x))$

In a previous paper, the authors tried to used a loss function which bears some similarity to the hinge function used in SVM. It includes a hyper-parameter which is a threshold in Hamming space that differentiates neighbors from non-neighbors. such that similar points are mapped to binary codes that do differ in more than P bits and disimilar points should map to points closer no more than P bits. For two binary codes $h$ and $g$ with hamming distance $||h-g||_H$ and a similarity label $s \in {0,1}$ the pairwise hinge loss function is defined as:

However in practice finding value of P is not easy. Moreover in some datasets the relative pairwise distance is important not the precise numerical value. As a results in this paper authors define the loss function in terms of the relative similarity. To define relative similarity it is assumed that dataset include triplet of items $(x,x^+,x^-)$ such that $x$ is more similar to $x^+$ than $x^-$. With this assumption the ranking loss on triplet of binary codes $(h,h^+,h^-)$ is:

$L_{triple}(h,h^+,h^-)=[||h-h^+||-||h-h^-||+1]_+$

Optimization

Given a training set of triplets, $D={(x_i,x_i^+,x_i^-)}_{i=1}^n$, the objective is to minimize the sum of ranking loss for all training samples and a simple regularizer on the vector of unknown parameters $w$:

$L(w)=\sum_{(x,x^+,x^-) \in D}l_{triple}(b(x,w),b(x^+,w),b(x^-,w))+\frac{\lambda}{2}||w||$
which is a discontinuous and non-convex function and optimization is not trivial. The discontinuity is because of the sign function and can be mitigated through construction of upper bound on the empirical loss. To do that one can rewrite function b as following
$b(x,w)=sign(f(x,w) = \underset{h\in H} {arg\,max} h^Tf(x,w)$

where $H={-1,+1}^q$.

Upper bound on empirical loss

The upper bound on loss has the following form:

$l_{triple}((b(x,w),b(x^+,w),b(x^-,w)) \leq max_{g,g^+,g^-}\{l_{triple}(g,g^+,g^-)+\g^Tf(x,w)+g^^+^Tf(x^+,w)=g^-^Tf(x^-,w)\}$