Projection Bank: From High-dimensional Data to Medium-length Binary Codes

09/16/2015 ∙ by Li Liu, et al. ∙ Northumbria University IEEE 0

Recently, very high-dimensional feature representations, e.g., Fisher Vector, have achieved excellent performance for visual recognition and retrieval. However, these lengthy representations always cause extremely heavy computational and storage costs and even become unfeasible in some large-scale applications. A few existing techniques can transfer very high-dimensional data into binary codes, but they still require the reduced code length to be relatively long to maintain acceptable accuracies. To target a better balance between computational efficiency and accuracies, in this paper, we propose a novel embedding method called Binary Projection Bank (BPB), which can effectively reduce the very high-dimensional representations to medium-dimensional binary codes without sacrificing accuracies. Instead of using conventional single linear or bilinear projections, the proposed method learns a bank of small projections via the max-margin constraint to optimally preserve the intrinsic data similarity. We have systematically evaluated the proposed method on three datasets: Flickr 1M, ILSVR2010 and UCF101, showing competitive retrieval and recognition accuracies compared with state-of-the-art approaches, but with a significantly smaller memory footprint and lower coding complexity.



There are no comments yet.


page 7

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

Recent research shows very high-dimensional feature representations, e.g., Fisher Vector (FV) [23, 27, 22] and VLAD [11]

, can achieve state-of-the-art performance in many visual classification, retrieval and recognition tasks. Although these very high-dimensional representations lead to better results, with the emergence of massive-scale datasets, e.g., ImageNet

[4] with around 15M images, the computational and storage costs of these long data have become very expensive and even unfeasible. For instance, if we represent 15M samples using 51200-dimensional FVs, the storage requirement of these data is approximately 5.6TB and it will need about

arithmetic operations measuring the Euclidean distance for image retrieval on these data. Considering the trade-off between computational efficiency and performance, it is desirable to embed the high-dimensional data into a reduced feature space. However, traditional dimensionality reduction methods such as PCA


are not suitable for large-scale/high-dimensional cases. The main reasons are: (1) Most dimensionality reduction methods are based on full-matrix linear projections, which need massive computational complexity and memory storage in high-dimensional reduction circumstances; (2) The reduced representations are usually real-valued vectors. When both dimensionality and the number of samples are large, real-valued codes severely limit the efficiency for retrieval and classification tasks compared with the binary codes. Thus, recent binarization approaches

[7, 8, 25, 26, 34, 19, 14, 3, 2, 36, 23, 27, 6, 17] have been proposed to embed the original data into binary codes with a reduced dimension. The codes generated by these methods can be roughly divided into two groups, i.e., the short binary codes and the long binary codes.

Hashing short codes: Most hashing-based approaches designed for fast searching always embed relatively low-dimensional representations like GIST [21] into short binary codes (usually under 500 bits) without too much loss of information. However, recent sophisticated and state-of-the-art representations are always over ten thousand dimensions. Therefore, these hashing methods become not so effective and appropriate for the embedding of very high-dimensional data, since they cannot preserve the sufficient discriminative properties to maintain high performance if the length of the obtained binary codes is short according to [23, 27]. Although some of hashing methods can theoretically generate long binary codes for high-dimensional data, the enormous computational load and memory usage make them unpractical. For instance, one of the state-of-the-art hashing methods, Iterative Quantization (ITQ) [7], leads to unacceptable loss of retrieval accuracy compared with long codes over 1000 bits, meanwhile its computational cost becomes extremely high when the number of bits increases as shown in Fig. 1(a). Long binary codes: Opposite to hashing short codes, a few methods [6, 36] have been specially introduced to convert the high-dimensional data to long binary codes. Among them, one representative binary coding method is Bilinear Projection-based Binary Codes (BPBC) [6], which can learn two rotate matrices for efficient binary coding. However, to minimize the loss of accuracies and achieve state-of-the-art performance, the left of Fig. 1(a) shows the length of binary codes generated by BPBC have to be long enough (i.e., over 10000 bits). Such long binary codes are still not fast enough for large-scale applications. Thus, how to learn medium-length binary codes (i.e., between 1000 bits and 10000 bits) and still maintain the high accuracies becomes a challenging research topic.

Many existing binary code learning methods are based on a single linear projection matrix (e.g., random projection) to map the data from the high-dimensional space to a reduced space. No matter the training or the coding phase, the storage requirement for the single linear projection matrix remains a burden. Taking ITQ (PCA is involved as its first step) for an example, to reduce the 170400-dimensional FVs to 10000-dimension, the size of the given single projection matrix should be about 12.7GB and the number of multiplications for coding a new data sample is . Apparently, this kind of memory requirements and coding complexity of the linear projection is unrealistic for large-scale applications. To alleviate this weakness, a bilinear projection method [6] has been proposed to effectively reduce the complexity of code learning compared with the linear one. The middle and right of Fig. 1(a) illustrate the memory usage and multiplications of coding respectively for the linear projection based ITQ and the bilinear projection based BPBC. Although the size of the projection matrix for the bilinear method is dramatically reduced, the coding complexity is still relatively high, especially when the dimensionality of original data goes high. Therefore, our target is to further reduce the coding complexity and produce medium-length codes without sacrificing the accuracy.

In this paper, we propose a novel binarization method for high-dimensional data. The proposed method first de-aggregates the original very high-dimensional representations into several groups of short representations according to their intrinsic data properties along the dimensions. After that, for each group of short representations, a small projection will be learned via the max-margin constraint to optimally preserve the data similarity. We denote our method as Binary Projection Bank (BPB), since a bank of small projections will be finally generated in our method instead of learning conventional linear or bilinear projections as illustrated in Fig. 1(b). The contributions of this paper include: (1) We propose a medium-length binary code learning method, which outperforms state-of-the-art linear and bilinear methods; (2) In spite of the reduced code length, our method only requires low and constant memory usage and coding complexity; (3) A kernelized version (KBPB) has also been proposed for better performance.

2 Related Work

There are a few works specifically focusing on high-dimensional data reduction. One of popular methods is Product Quantization (PQ) [10]

. Prior to PQ, however, a random rotation is always needed to balance the variance of high-dimensional data according to 

[11]. As we discussed before, such rotation requires high computational complexity. Recently, an efficient high-dimensional reduction method based on feature merging [5, 9, 16], termed Pseudo-supervised Kernel Alignment (PKA) [15], has achieved good performance but with cheaper computation. Besides, aiming for large-scale tasks, some binary reduction techniques for high-dimensional data have also been introduced. Perronnin et al. [23] proposed the “” binarization scheme and compared with Locality Sensitive Hashing (LSH) [3] and Spectral Hashing (SpH) [34] on the compressed FVs. Hashing Kernel (HK) [28] is utilized for high-dimensional signature compression as well in [27]. Most recently, the Bilinear Projection-based Binary Codes (BPBC) [6] is proposed to achieve more efficient binary coding. However, experiments in [23, 27, 6] manifest these methods require very long codes to yield acceptable performance.

3 Binary Projection Bank

3.1 Notation and Motivation

We are given training data in a -dimensional space: . The goal of this paper is to generate binary codes for each data point, for which the similarity and structure in the original data is preserved.

Traditional algorithms take -dimensional data as input and use a single projection to form the linear prediction function . Actually it can be regarded that projection matrix consists of

linear classifiers (binary output:

) over the original feature space. However, for realistic high-dimensional data with noise and redundancy of dimension, learning single projections across the entire high-dimensional feature space is unwise and costs very high computational complexity. To tackle this problem, we aim to split the original high-dimensional feature space into subspaces by merging the similar dimensions together similar to [15]. In this way, linear classifiers (i.e., projection vectors) will be explored since a subspace spanned by the dimensions with the similar property should only require one linear classifier. Hence, for each of subspace, only one small projection vector will be learned in our architecture. Table 1 summarizes the resource requirements for different projection schemes.

Storage Coding complexity
Table 1: Storage and coding complexity of different projection schemes. The sizes of two matrices in the bilinear projection are and .

In this paper, we propose a Binary Projection Bank (BPB) algorithm. To effectively decompose the original data space into subspaces, we first employ a K-means clustering scheme along dimensions. Particularly, we take

as K-means input, which is the rows of matrix , and divide them into clusters consisting of as illustrated:


As mentioned in [15], the decomposition via K-means can successfully preserve the intrinsic data structure and simultaneously group the dimension with the similar property together. We denote the number of the dimensions in the -th cluster by , where . Then the subspace spanned by the dimensions in the -th cluster is . In each subspace, a linear classifier is learned to generate one bit for the data. Hence, and is the reduced dimension. In the following, we give the detailed formulations and learning steps of BPB.

3.2 Formulation of BPB

To preserve the data similarity, we first construct a pseudo label for each data pair according to their k-nearest neighbors in original data space as follows:


where is the set of k-nearest neighbors of . Besides, we also define for . In BPB learning phase, our goal is to minimize the distances of positive pairs and maximize the distances of the negative pairs in each subspace generated by K-means clustering.

For the -th subspace, we use

to represent the data in this subspace. We tactfully transfer learning projections to learning linear classifiers via pair-wise labels generated from unlabeled data. By adopting linear classifier

similar to the SVM framework, positive pairs are positioned in the same side of the hyperplane while negative pairs are expected to be placed at different sides of the hyperplane.

In fact, we can denote , then the classifier becomes . Therefore, it is equivalent to the linear classifier without the bias . In the following computation, we omit and the binary code for each data in the -th subspace can be acquired as follows:


where is the coefficient of the classifier for the -th subspace. With the above requirement and the maximum margin criterion for the positive and negative pairs, we have the following optimization problem:


where is for the margin regularization. It is noticeable that if , the constraint becomes , which strictly constrains that every point is out of the margin. Using the hinge loss term, we can rewrite the optimization problem in (4) as the following objective function:


where is the balance parameter to control the importance of the two terms. Since we cannot directly obtain the optimal in our objective function, a gradient descent scheme has to be applied here. Let us denote


Taking the derivative of with respect to , we can easily obtain the gradient of :


Note that in our implementation, if , we can set , where is a small nonzero random vector. The same scheme has also been used in [13, 32].

Therefore, we utilize the gradient descent method and have the following update rule to optimize :


where is the step length.

We repeat the optimization problem in (4) for all the subspaces and concatenate the binary bits together to form the final binary code. The final binary code for high-dimensional data can be illustrated as:

Adaptive gradient descent (AGD):

Furthermore, for fast convergence, we also associate the optimization procedure with an adaptive step length. We first initialize . For the -th iteration, if , we enlarge the step length in the next iteration to accelerate the convergence, otherwise, we decrease to its half size: . In the experiments, we also set an upper bound for the number of iteration for the gradient descent. Thus, we stop the iteration when the number of iteration reaches a maximum or the difference is less than a small threshold.

Figure 2: Comparison of the average pairwise error on subspaces and pairwise error on original data space with respect to reduced code length on the ILSVR2010 dataset with 64000-d FV. All the results are the means of 50 runs with 70-iteration of AGD.

Average pairwise error of subspaces vs. pairwise error of uncompressed data:

We also analyze the empirical error rate of pairwise data in the projected space. Suppose is the solution of minimizing in Eq.(5) acquired by AGD for the -th subspace, . Then the average error rate on pair-wise data for all the subspaces is defined as , where represents the cardinality of the set. On the other hand, we also compute the solution of minimizing , for the original -dimensional data . Finally, we compare with the error rate in Fig. 2. It is observed that is lower than at the medium code length (under 10000 bits), which indicates that the data distribution in subspaces has much better separability than the one in the original space. However, when , the number of dimensions in each subspace will shrink to a very small value. In this case, data in subspaces are difficult to be linearly separated by classifiers for the current BPB. Thus, we will extend BPB to the version with non-linear kernels for better performance.

3.3 Kernel BPB

In this section, we introduce our algorithm with kernel functions, i.e., kernel BPB (KBPB), since the kernel method can theoretically and empirically be able to solve the linear inseparability problem mentioned in above. Although the kernel method would cost high computational complexity for high-dimensional data, in our method, the kernel function will only be performed in small subspaces which are spanned by the dimensions in each cluster.

In the -th subspace, suppose data are mapped to a Hilbert space by a mapping function and the kernel function is the inner product function in the Hilbert space. As defined in the Kernelized Locality-Sensitive Hashing (KLSH) [12] and Kernel-Based Supervised Hashing (KSH) [18], we uniformly select samples from the training data (we call them basis samples) to reduce the coding complexity of the kernel trick from to (the effect on selection of will be discussed in experiments). Then we establish the prediction function with the kernel as follows (without loss of generality, suppose the choices of basis samples are the first samples for all the subspaces):


where , are the coefficients and is the bias. It is actually the linear classifier for the data in the Hilbert space. The binary codes with sufficient information should be zero-centered [18, 12, 7], which renders that . To satisfy this condition, we set . Introducing in Eq. (10), the prediction function becomes:


where , . It is easy to observe that is the coefficient vector of the hyperplane for the data in the Hilbert space. With the similar constraints, we have the following optimization problem for KBPB:


where . Naturally, the corresponding objective function will be:

where is the balance parameter as in BPB. Let us denote and Then their derivatives with respect to can be computed as:


where and , Similar to BPB, if , we can set , where is a small nonzero random vector. Therefore, we have the update rule for KBPB as follows:


where is the step length, which is also adaptively tuned by AGD.

Finally, having calculated the coefficients of the kernelized prediction function for all the subspaces, the binary codes for the original data can be expressed as:

4 Experiments

4.1 Large-scale image retrieval

The proposed BPB and KBPB algorithms are first evaluated for the image similarity search task. Two realistic large-scale image datasets are used in our experiments: Flickr 1M and ILSVR2010. For Flickr 1M, we downloaded close to one million web images with 55 groups from Flickr inspired by [30, 27]. For each image in Flickr 1M, we extract -d SIFT features in patches of around interest points detected by [20]. The ILSVR2010111 dataset is a subset of the ImageNet [4] dataset and contains 1.2 million images from categories. The publicly available dense -d SIFT features [4] are used.

We represent each image in both datasets using two high-dimensional representations: Fisher Vector (FV) and VLAD. In respect to FV, the Gaussian Mixture Model is implemented on SIFT features with

Gaussians for both datasets. In this way, the dimension of the final FV for each image is . While, for VLAD representations, the K-means clustering has been used to cluster the SIFT features into centers and aggregate them into VLAD vectors of dimensions. These VLAD vectors are also power and normalized [24]. In terms of both datasets, we randomly select images as the query and the remaining images are regarded as the gallery database. For evaluation, we first report the semantic precision at and retrieved images (according to the ground-truth) for both Flickr 1M and ILSVR2010, and then the precision-recall curves are illustrated as well. Additionally, we report the size of projection storage and the coding time (the average time used for each data) for some state-of-the-art methods. Our experiments are completed using Matlab 2014a on a server configured with a 6-core processor and 64GB of RAM running the Linux OS.

Compared methods and settings:

In our experiments, we compare the proposed method with nine coding methods including four real-valued dimensionality reduction methods: Principal Component Analysis (PCA), the projection via Gaussian random rotation (RR), Product Quantization (PQ)

[10] and Pseudo-supervised Kernel Alignment (PKA) [15], and five binary coding methods: the sign function binarization, “” binarization [23], Locality Sensitive Hashing (LSH) [3], Spectral Hashing (SpH) [34], Bilinear Projection-based Binary Codes (BPBC) [6] and Circulant Binary Embedding (CBE) [36]. We use the publicly available codes of LSH, SpH, PQ , CBE and PCA, and implement RR, PKA and BPBC ourselves. Additionally, two natural baselines: randomly sampling the dimensions to form subspaces without replacement (RandST+BPB) and learning multiple bits with ITQ in each subspace (Kmeans+ITQ) are also included in our experiments. All of the above methods are then evaluated for compressing FV and VLAD representations into three different medium-lengthed codes: . Considering the feasibility on the training phase of all the methods, in this experiment, 150K data are randomly selected from the gallery database of Flickr 1M and ILSVR2010 respectively to form the training set. Besides, we also randomly choose another 50K data samples from each of the datasets as a cross-validation set for parameter tuning. Under the same experimental setting, all the parameters used in the compared methods have been strictly chosen according to their original papers.

For the proposed BPB/KBPB, the pairwise label of each data pair is determined by their nearest neighbors. The balance parameter for each dataset is selected from one of the values in the range of , which yields the best performance on the cross-validation set. The maximum number of the iteration of AGD is fixed at , which has been proved to converge well for the objective function. For KBPB, we adopt as the number of basis samples for both Flickr 1M and ILSVR2010. We use the polynomial kernel and the RBF kernel to implement KBPB and KBPB, respectively. The best value of for KBPB is selected via cross-validation and the value of for KBPB is determined adaptively based on the method in [1]. In fact, any kernel function satisfying the Mercer’s condition can be used in KBPB. In BPB/KBPB, since the coding procedure in each subspace is independent, we implement the parallel computation scheme to speed up the training time. Considering the uncertainty of the K-means clustering, all the reported results by our methods are the averages of 50 runs.

Methods Fisher Vector (64000-d) VLAD (32000-d)
Precision@top 50 Precision@top 100 Precision@top 50 Precision@top 100
8000 bit 6400 bit 4000 bit 8000 bit 6400 bit 4000 bit 4000 bit 3200 bit 2000 bit 4000 bit 3200 bit 2000 bit
Original 0.383 0.383 0.383 0.355 0.355 0.355 0.370 0.370 0.370 0.339 0.339 0.339
PCA 0.371 0.352 0.314 0.354 0.316 0.284 0.365 0.341 0.302 0.337 0.295 0.261
PQ 0.160 0.132 0.121 0.154 0.128 0.111 0.142 0.114 0.103 0.138 0.109 0.090
PKA 0.380 0.374 0.320 0.352 0.337 0.301 0.376 0.366 0.307 0.335 0.317 0.276
RR+PQ 0.292 0.279 0.251 0.283 0.268 0.234 0.291 0.267 0.233 0.265 0.249 0.213
Sign 0.281 0.281 0.281 0.277 0.277 0.277 0.271 0.271 0.271 0.262 0.262 0.262
0.273 0.273 0.273 0.262 0.262 0.262 - - - - - -
LSH 0.319 0.294 0.270 0.304 0.289 0.267 0.315 0.287 0.250 0.287 0.272 0.244
SpH 0.259 0.288 0.301 0.244 0.273 0.296 0.254 0.267 0.277 0.225 0.240 0.265
BPBC 0.343 0.338 0.314 0.328 0.303 0.289 0.328 0.312 0.294 0.311 0.281 0.267
CBE 0.382 0.379 0.377 0.356 0.351 0.342 0.360 0.351 0.342 0.332 0.327 0.321
Kmeans+ITQ(20bits) 0.365 0.338 0.313 0.338 0.320 0.310 0.363 0.350 0.332 0.324 0.311 0.299
RandST+BPB 0.352 0.335 0.321 0.331 0.316 0.303 0.351 0.345 0.331 0.319 0.305 0.296
BPB 0.385 0.381 0.376 0.356 0.349 0.345 0.365 0.357 0.350 0.338 0.333 0.329
KBPB 0.391 0.388 0.375 0.358 0.353 0.347 0.370 0.362 0.355 0.342 0.340 0.332
KBPB 0.398 0.391 0.379 0.367 0.359 0.350 0.378 0.376 0.363 0.349 0.345 0.336

The “Original” indicates uncompressed FV/VLAD. The “Sign” refers to directly using the sign function on original vectors. “[23] scheme is specifically designed for FV and the dimension of reduced codes via “” is fixed at =. KBPB indicates KBPB with the polynomial kernel and KBPB indicates KBPB with the RBF kernel. The results of BPB and KBPB are mean accuracies of 50 runs. For Original, PCA, PKA and RR, the Euclidean distance is used to measure the retrieval. For RR+PQ, the asymmetric distance (ASD) [10] is adopted and Hamming distance is used for the rest of compared methods. Kmeans+ITQ(20bits) indicates using Kmeans to split the dimensions into subspaces and then apply ITQ to learning 20 bits codes for each subspace. RandST+BPB denotes randomly split the dimensions into subspaces without replacement and adopt BPB optimization scheme to learn codes.

Table 2: Retrieval results (semantic precision) comparison on Flickr 1M with 64000-dimensional FV and 32000-dimensional VLAD.
Methods Fisher Vector (64000-d) VLAD (32000-d)
Precision@top 50 Precision@top 100 Precision@top 50 Precision@top 100
8000 bit 6400 bit 4000 bit 8000 bit 6400 bit 4000 bit 4000 bit 3200 bit 2000 bit 4000 bit 3200 bit 2000 bit
Original 0.214 0.214 0.214 0.175 0.175 0.175 0.186 0.186 0.186 0.149 0.149 0.149
PCA 0.185 0.177 0.159 0.154 0.138 0.125 0.176 0.160 0.143 0.139 0.116 0.101
PQ 0.157 0.153 0.138 0.132 0.125 0.110 0.154 0.142 0.127 0.115 0.106 0.089
PKA 0.214 0.205 0.183 0.179 0.168 0.155 0.188 0.172 0.157 0.152 0.143 0.124
RR+PQ 0.182 0.163 0.157 0.151 0.137 0.128 0.174 0.166 0.151 0.132 0.117 0.108
Sign 0.152 0.152 0.152 0.131 0.131 0.131 0.141 0.141 0.141 0.109 0.109 0.109
0.175 0.175 0.175 0.144 0.144 0.144 - - - - - -
LSH 0.185 0.171 0.160 0.162 0.155 0.143 0.175 0.159 0.143 0.144 0.128 0.117
SpH 0.160 0.169 0.153 0.139 0.151 0.123 0.136 0.152 0.149 0.108 0.127 0.120
BPBC 0.187 0.183 0.176 0.165 0.153 0.146 0.180 0.173 0.154 0.139 0.127 0.120
CBE 0.214 0.207 0.197 0.178 0.164 0.155 0.184 0.174 0.167 0.145 0.138 0.131
Kmeans+ITQ(20bits) 0.201 0.197 0.184 0.167 0.158 0.150 0.172 0.165 0.155 0.137 0.131 0.125
RandST+BPB 0.202 0.193 0.182 0.161 0.149 0.140 0.175 0.164 0.155 0.138 0.129 0.121
BPB 0.210 0.205 0.197 0.175 0.166 0.158 0.186 0.177 0.165 0.148 0.142 0.134
KBPB 0.218 0.208 0.201 0.177 0.172 0.163 0.190 0.181 0.171 0.151 0.145 0.137
KBPB 0.226 0.221 0.208 0.181 0.176 0.168 0.198 0.190 0.180 0.157 0.150 0.143
Table 3: Retrieval results comparison (semantic precision) on ILSVR2010 with 64000-dimensional FV and 32000-dimensional VLAD.

Results comparison:

We list the retrieval results comparison of different methods at top 50 and 100 retrieval results on the Flickr 1M and ILSVR2010 datasets in Table 2 and Table 3, respectively. Generally, FV gains slightly better results than VLAD on both datasets. Meanwhile, the accuracies on the ILSVR2010 dataset are lower than those on the Flickr 1M dataset, since there are more categories and larger intra-class variations in ILSVR2010. It is noticeable that PQ achieves the low precision on Flickr 1M, while RR+PQ can lead to more reasonable results. The reason is that for high-dimensional representations, there may exist unbalanced variance that influences the performance. Thus, randomly rotating the high-dimensional data prior to PQ222In [10], PQ can achieve competitive results without random rotation. However, they focus on relatively low-dimensional SIFT/GIST features whose variance already tends to be roughly balanced. is recommended in [11]. Nevertheless, due to that the images in ILSVR2010 are textured with the dominant object which leads to relatively balanced variance, the basic PQ can achieve modest results on ILSVR2010. PCA and PKA have remarkable accuracies as real-valued compression techniques on both datasets and CBE is regarded as the strongest baseline of binary coding methods according to its performance. LSH, SpH and the “” scheme can obtain similar results on both datasets and using sign function directly on uncompressed FV/VLAD is proved to be the worst binarization method. Additionally, Kmeans+ITQ(20bits) can achieve slightly better performance than RandST+BPB, but both significantly lower than BPB.

From Table 2 and Table 3, our BPB algorithm consistently outperforms all the compared methods at every code length and leads to competitive accuracies with CBE and original FV/VLAD. Moreover, KBPB can achieve better performance than BPB since the kernel method can theoretically and empirically solve the problem of linear inseparability of subspaces with relatively lower dimensions (average dimension of each subspace is ). Thus, KBPB gives significantly better performance when is large, i.e., on relatively long binary codes. The best performance on both datasets has been achieved by KBPB with the RBF kernel. Especially, when the code length decreases, the retrieval accuracies from all compared methods (expect SpH) dramatically drop, but the accuracies of our methods only slightly change showing the robustness of the proposed methods on medium-dimensional binary coding. Currently, we use hard-assignment K-means for our work. In Fig. 3, we have also evaluated the possibility to use soft-assignment clustering for our methods. The results illustrate that for the medium-dimensional codes (i.e., between 1000 bits and 10000 bits), soft-assignment clustering based BPB and RBF-KBPB can achieve competitive results with hard-assignment BPB and KBPB. However, in the extreme condition (i.e., code dimfeature dim), the soft-clustering based methods can still produce the reasonable results without lossing much of accuracy, while the current hard-clustering methods fail, since the feature dimensions can be re-used during soft-clustering. Thus, from Fig. 3, Table 2 and Table 3, we can observe our current version of projection bank can indeed achieve better performance for the medium-dimensional codes compared with other methods. Besides, Fig. 4 presents the precision-recall curves of all compared methods on both datasets with 8000 bits for FV and 4000 bits for VLAD, respectively. From all these figures, we can further discover that, for both datasets, BPB/KBPB outperform other high-dimensional compression methods with the medium-lengthed codes by comparing the retrieval precision and the Area Under the Curve (AUC).

Figure 3: Comparison with soft-kmeans on BPB
3200 Projections 1562.53 1576.37 0.34 0.49 0.49 24.40
4000 storage 1953.19 1970.49 0.42 0.49 0.49 30.51
6400 (MB) 3125.04 3152.50 0.65 0.49 0.49 48.83
8000 3906.38 3940.38 0.80 0.49 0.49 61.17
3200 Training 2.21 1102.72 848.21 304.15 314.27 603.52
4000 time 3.30 1563.48 1021.32 340.24 389.56 718.33
6400 (s) 4.52 1918.52 1336.10 422.44 472.33 802.51
8000 8.64 2204.30 1559.08 464.23 505.16 890.24
3200 Coding 246.71 599.47 5.46 21.64 0.24 2.41
4000 time 370.24 791.43 13.45 22.32 0.23 6.64
6400 (ms) 534.10 1182.57 30.57 21.03 0.24 19.10
8000 721.51 1603.10 80.02 23.15 0.24 55.91

The training and coding time in this table are both 50-run averaged runtime. Parallel computation is adopted to speed up our training phase.

Table 4: Comparison of computational cost on different code lengths for 64000-d FV reduction on ILSVR2010 (Data stored in double precision).
(a) Flickr 1M (FV) (b) Flickr 1M (VLAD) (c) ILSVR2010 (FV) (d) ILSVR2010 (VLAD)
Figure 4: Comparison of precision-recall curves on Flickr 1M and ILSVR2010 datasets with 8000 bits FV and 4000 bits VLAD.
(a) Flickr 1M (FV) (b) Flickr 1M (FV) (c) ILSVR2010 (VLAD) (d) ILSVR2010 (VLAD)
Figure 5: (a) and (c) show the mean of 50 runs of retrieval accuracies of KBPB (with the RBF kernel) vs. parameter on Flickr 1M and ILSVR2010. (b) and (d) show the parameter sensitivity analysis of on Flickr 1M and ILSVR2010 at 6400 bits and 3200 bits, respectively.

Complexity and parameter sensitivity analysis:

Table 4 illustrates the comparison of the memory usage for projections, the training time and the coding time on ILSVR2010. RR+PQ costs the largest memory space and more time for coding since the full-matrix projection (i.e., RR) is involved. Compared with RR+PQ, BPBC and CBE need much lower memory costs and time complexity for training and coding. Our BPB is slightly time-consuming than CBE in the training phase but the most efficient one for coding. Meanwhile, KBPB costs more memory space than BPBC and BPB but is still more efficient for coding than BPBC. In addition, Fig. 5 reports the effect of performance by varying two essential parameters and . In terms of the number of basis samples used in KBPB with the RBF kernel, when , the retrieval accuracy curves become approximately stable on both datasets with FV and VLAD, respectively. It indicates that our KBPB can lead to relatively robust results with . As we can see, for balance parameter , our methods (both BPB and KBPB) can achieve the good performance when and on Flick 1M and ILSVR2010, respectively.

4.2 Large-scale action recognition

Finally, we evaluate our methods for action recognition on the UCF101 dataset [29] which contains 13320 videos from 101 action categories. We strictly follow the 3-split train/test setting in [29] and report the average accuracies as the overall results. The 426-dimensional default Dense Trajectory Features (DTF) [31] are extracted from each video, and GMM and K-means are used to cluster them into 200 visual words for FV and VLAD respectively. Thus, the length of FV is and the length of VLAD is . For our methods, we fix and , which are both selected via cross-validation set, and other parameters are the same as the previous retrieval experiments. In this experiment, we apply the linear SVM333According to [27, 6], hashing kernel [28, 33]

renders to an unbiased estimation of the dot-product in the original space. Thus, binary codes can also be directly fed into a linear SVM.

for action recognition. From the relevant results shown in Table 5, it can be observed that the recognition accuracies computed by all methods have generally smaller differences compared with the diversity of performance in retrieval tasks. The reason is that the supervised SVM training can compensate the discriminative power between different methods, whereas the unsupervised retrieval cannot. Our BPB and KBPB can not only achieve competitive results with original features, but also perform better than other compression methods on medium-lengthed codes with FV and VLAD. Moreover, KBPB consistently gives the best performance.

Methods Fisher Vector (170400-d) VLAD (85200-d)
17040 bit 11360 bit 8520 bit 8520 bit 5680 bit 4260 bit
Original 80.33 80.33 80.33 77.95 77.95 77.95
PCA 78.62 78.31 75.4 77.03 76.28 74.1
RR+PQ 77.25 77.67 75.50 75.38 75.21 74.03
PKA 80.30 78.88 76.54 77.21 77.00 76.4
PQ 75.90 74.84 74.31 72.85 72.01 70.99
sign 75.26 75.26 75.26 74.41 74.41 74.41
76.78 75.20 74.56 - - -
LSH 74.19 73.02 71.88 72.40 71.11 70.4
SpH 71.36 73.04 75.28 69.35 72.97 74.83
BPBC 77.21 76.40 75.89 75.91 74.73 73.22
CBE 80.65 78.23 76.47 77.91 75.34 74.03
BPB 80.02 79.26 78.30 77.53 76.38 75.52
KBPB 80.74 80.35 79.37 78.28 77.31 76.54
KBPB 82.18 81.55 80.71 78.69 77.52 76.90
Table 5: Comparison of action recognition performance (%) on the UCF 101 dataset.

5 Conclusion and Future Work

In this paper, we have presented a novel binarization approach called Binary Projection Bank (BPB) for high-dimensional data, which exploits a group of small projections via the max-margin constraint to optimally preserve the intrinsic data similarity. Different from the conventional linear or bilinear projections, the proposed method can effectively map very high-dimensional representations to medium-dimensional binary codes with a low memory requirement and a more efficient coding procedure. BPB and the kernelized version KBPB have achieved better results compared with state-of-the-art methods for image retrieval and action recognition applications. In the future, we will focus more on using soft-assignment clustering based projection bank methods.


  • [1] X. Bai, X. Yang, L. J. Latecki, W. Liu, and Z. Tu. Learning context-sensitive shape similarity by graph transduction. T-PAMI, 32(5):861–874, 2010.
  • [2] Z. Cai, L. Liu, M. Yu, and L. Shao. Latent structure preserving hashing. In BMVC, 2015.
  • [3] M. S. Charikar. Similarity estimation techniques from rounding algorithms. In STOC, 2002.
  • [4] J. Deng, W. Dong, R. Socher, L.-J. Li, K. Li, and L. Fei-Fei. Imagenet: A large-scale hierarchical image database. In CVPR, 2009.
  • [5] B. Fulkerson, A. Vedaldi, and S. Soatto. Localizing objects with smart dictionaries. In ECCV. 2008.
  • [6] Y. Gong, S. Kumar, H. A. Rowley, and S. Lazebnik. Learning binary codes for high-dimensional data using bilinear projections. In CVPR, 2013.
  • [7] Y. Gong and S. Lazebnik. Iterative quantization: A procrustean approach to learning binary codes. In CVPR, 2011.
  • [8] J.-P. Heo, Y. Lee, J. He, S.-F. Chang, and S.-E. Yoon. Spherical hashing. In CVPR, 2012.
  • [9] H. Jegou, M. Douze, and C. Schmid. Hamming embedding and weak geometric consistency for large scale image search. In ECCV. 2008.
  • [10] H. Jegou, M. Douze, and C. Schmid. Product quantization for nearest neighbor search. PAMI, 33(1):117–128, 2011.
  • [11] H. Jégou, M. Douze, C. Schmid, and P. Pérez. Aggregating local descriptors into a compact image representation. In CVPR, 2010.
  • [12] B. Kulis and K. Grauman. Kernelized locality-sensitive hashing. T-PAMI, 34(6):1092–1104, 2012.
  • [13] N. Kwak. Principal component analysis based on l1-norm maximization. T-PAMI, 30(9):1672–1680, 2008.
  • [14] Y. Lin, R. Jin, D. Cai, S. Yan, and X. Li. Compressed hashing. In CVPR, 2013.
  • [15] L. Liu and L. Wang. A scalable unsupervised feature merging approach to efficient dimensionality reduction of high-dimensional visual data. In ICCV, 2013.
  • [16] L. Liu, L. Wang, and C. Shen. A generalized probabilistic framework for compact codebook creation. In CVPR, 2011.
  • [17] L. Liu, M. Yu, and L. Shao. Multiview alignment hashing for efficient image search. IEEE Transactions on Image Processing, 24(3):956–966, 2015.
  • [18] W. Liu, J. Wang, R. Ji, Y.-G. Jiang, and S.-F. Chang. Supervised hashing with kernels. In CVPR, 2012.
  • [19] W. Liu, J. Wang, S. Kumar, and S.-F. Chang. Hashing with graphs. In ICML, 2011.
  • [20] D. G. Lowe. Object recognition from local scale-invariant features. In ICCV, 1999.
  • [21] A. Oliva and A. Torralba. Modeling the shape of the scene: A holistic representation of the spatial envelope.

    International journal of computer vision

    , 42(3):145–175, 2001.
  • [22] F. Perronnin and C. Dance. Fisher kernels on visual vocabularies for image categorization. In CVPR, 2007.
  • [23] F. Perronnin, Y. Liu, J. Sánchez, and H. Poirier. Large-scale image retrieval with compressed fisher vectors. In CVPR, 2010.
  • [24] F. Perronnin, J. Sánchez, and T. Mensink. Improving the fisher kernel for large-scale image classification. In ECCV. 2010.
  • [25] J. Qin, L. Liu, M. Yu, Y. Wang, and L. Shao. Fast action retrieval from videos via feature disaggregation. In BMVC, 2015.
  • [26] R. Salakhutdinov and G. Hinton. Semantic hashing. International Journal of Approximate Reasoning, 50(7):969–978, 2009.
  • [27] J. Sánchez and F. Perronnin. High-dimensional signature compression for large-scale image classification. In CVPR, 2011.
  • [28] Q. Shi, J. Petterson, G. Dror, J. Langford, A. Smola, and S. Vishwanathan. Hash kernels for structured data. JMLR, 10:2615–2637, 2009.
  • [29] K. Soomro, A. R. Zamir, and M. Shah. Ucf101: A dataset of 101 human actions classes from videos in the wild. arXiv preprint arXiv:1212.0402, 2012.
  • [30] G. Wang, D. Hoiem, and D. Forsyth. Learning image similarity from flickr groups using fast kernel machines. T-PAMI, 34(11):2177–2188, 2012.
  • [31] H. Wang, A. Klaser, C. Schmid, and C.-L. Liu. Action recognition by dense trajectories. In CVPR, 2011.
  • [32] H. Wang, X. Lu, Z. Hu, and W. Zheng. Fisher discriminant analysis with l1-norm. IEEE Transactions on Cybernetics, 44(6):828–842, 2014.
  • [33] K. Weinberger, A. Dasgupta, J. Langford, A. Smola, and J. Attenberg. Feature hashing for large scale multitask learning. In ICML, pages 1113–1120, 2009.
  • [34] Y. Weiss, A. Torralba, and R. Fergus. Spectral hashing. In NIPS, 2008.
  • [35] S. Wold, K. Esbensen, and P. Geladi. Principal component analysis. Chemometrics and intelligent laboratory systems, 2(1):37–52, 1987.
  • [36] F. X. Yu, S. Kumar, Y. Gong, and S.-F. Chang. Circulant binary embedding. ICML, 2014.