Spectral grouping using the Nystrom method

Charless C. FowlkesSerge J. BelongieFan ChungJitendra Malik

article2004TPAMI1,527 citations

Proposes using the Nyström method to extrapolate spectral clustering solutions from a small subset of sample points, scaling image and video segmentation linearly with resolution while drastically cutting computational and memory costs.

Listen

Grouping visual elements into meaningful regions and objects is a foundational challenge in computer vision. While spectral graph-partitioning techniques such as Normalized Cut offer major advantages over prototype-based clustering by capturing complex cluster shapes and flexible similarity measures, their quadratic memory and cubic computational demands prevent practical deployment on high-resolution images and multi-frame video sequences.

The article demonstrates an approximation technique based on the classical Nyström method to substantially reduce the computational requirements of spectral grouping. The objective is to evaluate whether computing pairwise affinities for only a small subset of randomly sampled elements allows accurate extrapolation to full-resolution image and video segmentations.

The authors evaluate the approach through synthetic clustering benchmarks, comparative solver runtime analyses, and cross-validation across 300 natural scene images from the Corel dataset. They also test the method on practical image and video segmentation tasks incorporating color, texture, and temporal motion cues.

The key findings show that sampling fewer than one percent of total pixels—roughly 100 randomly chosen samples—reliably captures the leading graph eigenvectors and segments complex natural images. For a fixed sample size, the computational complexity scales linearly with the total number of pixels, transforming an intractable calculation into an efficient operation. In contrast to sparse eigensolvers whose convergence times degrade sharply on difficult problem instances, the Nyström approximation delivers robust and predictable performance, segmenting a multi-frame video volume in under one minute on standard hardware.

These results establish that spectral partitioning can be applied to large-scale vision tasks without arbitrary distance cutoffs or loss of long-range affinities. Practitioners can choose between a fast, one-step formulation when using positive definite similarity measures and a stable two-step approach when using general measures.

Organizations seeking to implement large-scale clustering or video segmentation should adopt the sampling-based extrapolation scheme to cut processing time and infrastructure costs. Future operational work should focus on automating the selection of the number of target clusters, as this parameter was set manually in the study.

The primary operational limitation is the requirement for appropriate sample sizing relative to scene complexity; undersampling complex scenes could cause small distinct objects to be missed. However, given the empirical stability across hundreds of natural images, confidence in the methodology's practical efficiency and grouping accuracy remains high.

  • Paper: A tutorial on spectral clustering, Ulrike von Luxburg (2007). This tutorial offers a comprehensive theoretical and practical survey of graph Laplacians and spectral clustering, consolidating the foundational concepts accelerated by Nyström-based approximations.
  • Paper: Self-Tuning Spectral Clustering, Lihi Zelnik-Manor et al. (2004). It extends spectral clustering methodology by introducing local scale estimation and automatic cluster count determination, addressing key operational limitations highlighted in the source paper.
  • Paper: Random Walks for Image Segmentation, Leo Grady (2006). It develops an alternative graph-based segmentation framework using discrete random walks and Dirichlet problems, offering a complementary linear-system perspective to spectral graph cuts.
  • Paper: Robust Subspace Segmentation by Low-Rank Representation, Guangcan Liu et al. (2010). It advances spectral segmentation to corrupted data by constructing affinity graphs via low-rank representation before applying spectral clustering.
  • Paper: Sparse Subspace Clustering: Algorithm, Theory, and Applications, Ehsan Elhamifar et al. (2012). It builds upon spectral clustering pipelines by applying l1-minimization self-expressiveness to form affinity graphs for high-dimensional subspace clustering.
  • Paper: Graph Regularized Nonnegative Matrix Factorization for Data Representation, Deng Cai et al. (2011). It combines graph Laplacian regularization with matrix factorization, generalizing manifold-preserving spectral concepts to parts-based representations.
Cover for Spectral grouping using the Nystrom method

Abstract

Spectral graph theoretic methods have recently shown great promise for the problem of image segmentation. However, due to the computational demands of these approaches, applications to large problems such as spatiotemporal data and high resolution imagery have been slow to appear. The contribution of this paper is a method that substantially reduces the computational requirements of grouping algorithms based on spectral partitioning making it feasible to apply them to very large grouping problems. Our approach is based on a technique for the numerical solution of eigenfunction problems known as the Nyström method. This method allows one to extrapolate the complete grouping solution using only a small number of samples. In doing so, we leverage the fact that there are far fewer coherent groups in a scene than pixels.

Table of Contents

  • 1 INTRODUCTION
  • 2 SPECTRAL METHODS FOR PAIRWISE CLUSTERING
  • 3 THE NYSTRÖM EXTENSION
  • 3.1 Background
  • 3.2 Matrix Completion
  • 3.3 Methods of Solution
  • 3.4 Application to Normalized Cut
  • 4 PERFORMANCE CONSIDERATIONS
  • 4.1 Approximation Properties
  • 4.2 Sampling
  • 5 SEGMENTATION RESULTS
  • 5.1 Color and Texture Segmentation
  • 5.2 Spatio-Temporal Segmentation
  • 6 CONCLUSION
  • APPENDIX A PROOF OF ONE-SHOT METHOD
  • APPENDIX B PROOF OF POSITIVE DEFINITENESS OF ε²ij
  • ACKNOWLEDGMENTS
  • REFERENCES

Knowls

  1. Knowl 1 — Nyström Matrix Completion for Large Affinity Matrices

    model/method

    Let W∈RN×NW \in \mathbb{R}^{N \times N} be a symmetric pairwise affinity matrix defined over NN data points (such as pixels or voxels). Let n≪Nn \ll N sample points be chosen, and denote m=N−nm = N - n. When the rows and columns of WW are permuted so that the nn samples appear first, WW is partitioned into subblocks as:

    W=[ABBTC]W = \begin{bmatrix} A & B \\ B^T & C \end{bmatrix}

    where A∈Rn×nA \in \mathbb{R}^{n \times n} represents affinities among the sample points, B∈Rn×mB \in \mathbb{R}^{n \times m} represents affinities between sample points and the remaining mm points, and C∈Rm×mC \in \mathbb{R}^{m \times m} represents affinities among the unselected points.

    The Nyström extension approximates the unobserved block CC by BTA−1BB^T A^{-1} B, yielding the low-rank approximated matrix W^\hat{W}:

    W^=[ABBTBTA−1B]=[ABT]A−1[AB]\hat{W} = \begin{bmatrix} A & B \\ B^T & B^T A^{-1} B \end{bmatrix} = \begin{bmatrix} A \\ B^T \end{bmatrix} A^{-1} \begin{bmatrix} A & B \end{bmatrix}

    When AA is rank-deficient or ill-conditioned due to sample redundancies, the Moore-Penrose pseudoinverse A+A^+ replaces A−1A^{-1}. The approximation error on the unobserved block is quantified by the Schur complement norm ∥C−BTA−1B∥\|C - B^T A^{-1} B\|, which measures how well the rows of BB span the pairwise interactions in CC.

  2. Knowl 2 — One-Shot Orthonormal Eigendecomposition for Positive Definite Submatrices

    theoretical result

    When the sample affinity submatrix A∈Rn×nA \in \mathbb{R}^{n \times n} of the partitioned affinity approximation W^=[ABT]A−1[AB]\hat{W} = \begin{bmatrix} A \\ B^T \end{bmatrix} A^{-1} \begin{bmatrix} A & B \end{bmatrix} is positive definite, the exact orthonormal eigendecomposition of W^∈RN×N\hat{W} \in \mathbb{R}^{N \times N} satisfying W^=VΛSVT\hat{W} = V \Lambda_S V^T and VTV=InV^T V = I_n is computed in a single n×nn \times n diagonalization step.

    Let A1/2A^{1/2} be the symmetric positive definite square root of AA, and define the symmetric matrix S∈Rn×nS \in \mathbb{R}^{n \times n} by:

    S=A+A−1/2BBTA−1/2S = A + A^{-1/2} B B^T A^{-1/2}

    Let S=USΛSUSTS = U_S \Lambda_S U_S^T be the eigendecomposition of SS, where US∈Rn×nU_S \in \mathbb{R}^{n \times n} is orthogonal (USTUS=InU_S^T U_S = I_n) and ΛS∈Rn×n\Lambda_S \in \mathbb{R}^{n \times n} is the diagonal matrix of eigenvalues. The N×nN \times n matrix VV containing the leading orthonormal eigenvectors of W^\hat{W} is:

    V=[ABT]A−1/2USΛS−1/2V = \begin{bmatrix} A \\ B^T \end{bmatrix} A^{-1/2} U_S \Lambda_S^{-1/2}

    This construction satisfies both VTV=InV^T V = I_n and VΛSVT=W^V \Lambda_S V^T = \hat{W} without ever explicitly computing or materializing the (N−n)×(N−n)(N-n) \times (N-n) block BTA−1BB^T A^{-1} B.

  3. Knowl 3 — Implicit Degree Vector Computation for Nyström Normalized Cut

    equation

    In Normalized Cut spectral clustering, constructing the normalized affinity matrix D−1/2WD−1/2D^{-1/2} W D^{-1/2} requires the row-sum degree vector d=W1Nd = W \mathbf{1}_N, where D=diag(d)D = \text{diag}(d) and 1k\mathbf{1}_k denotes a kk-dimensional column vector of all ones. For the Nyström-approximated affinity matrix W^=[ABBTBTA−1B]\hat{W} = \begin{bmatrix} A & B \\ B^T & B^T A^{-1} B \end{bmatrix} with A∈Rn×nA \in \mathbb{R}^{n \times n} and B∈Rn×mB \in \mathbb{R}^{n \times m} (m=N−nm = N - n), the approximated degree vector d^=W^1N∈RN\hat{d} = \hat{W} \mathbf{1}_N \in \mathbb{R}^N is computed directly without calculating BTA−1BB^T A^{-1} B as:

    d^=[A1n+B1mBT1n+BTA−1B1m]=[ar+brbc+BTA−1br]\hat{d} = \begin{bmatrix} A \mathbf{1}_n + B \mathbf{1}_m \\ B^T \mathbf{1}_n + B^T A^{-1} B \mathbf{1}_m \end{bmatrix} = \begin{bmatrix} a_r + b_r \\ b_c + B^T A^{-1} b_r \end{bmatrix}

    where ar=A1n∈Rna_r = A \mathbf{1}_n \in \mathbb{R}^n and br=B1m∈Rnb_r = B \mathbf{1}_m \in \mathbb{R}^n are the row sums of AA and BB, respectively, and bc=BT1n∈Rmb_c = B^T \mathbf{1}_n \in \mathbb{R}^m is the column sum vector of BB.

    Given d^\hat{d}, the entries of the normalized subblocks Aˉ∈Rn×n\bar{A} \in \mathbb{R}^{n \times n} and Bˉ∈Rn×m\bar{B} \in \mathbb{R}^{n \times m} of D^−1/2W^D^−1/2\hat{D}^{-1/2} \hat{W} \hat{D}^{-1/2} (where D^=diag(d^)\hat{D} = \text{diag}(\hat{d})) are computed as:

    Aˉij=Aijd^id^j,i,j∈{1,…,n}\bar{A}_{ij} = \frac{A_{ij}}{\sqrt{\hat{d}_i \hat{d}_j}}, \quad i, j \in \{1, \dots, n\}

    Bˉij=Bijd^id^j+n,i∈{1,…,n},  j∈{1,…,m}\bar{B}_{ij} = \frac{B_{ij}}{\sqrt{\hat{d}_i \hat{d}_{j+n}}}, \quad i \in \{1, \dots, n\}, \; j \in \{1, \dots, m\}

  4. Knowl 4 — Nyström Normalized Cut Spectral Clustering Algorithm

    algorithm

    The Nyström Normalized Cut algorithm computes approximate leading spectral embedding coordinates for NN data points using affinities evaluated on nn randomly sampled points (n≪Nn \ll N) and partitions the dataset via kk-means.

    Input: Unnormalized affinity submatrices A∈Rn×nA \in \mathbb{R}^{n \times n} and B∈Rn×mB \in \mathbb{R}^{n \times m} (m=N−nm = N - n), number of desired eigenvectors NEN_E, number of clusters kk
    Output: Cluster assignment vector c∈{1,…,k}Nc \in \{1, \dots, k\}^N
    ar←A1na_r \leftarrow A \mathbf{1}_n
    br←B1mb_r \leftarrow B \mathbf{1}_m
    bc←BT1nb_c \leftarrow B^T \mathbf{1}_n
    d1←ar+brd_1 \leftarrow a_r + b_r
    d2←bc+BTA−1brd_2 \leftarrow b_c + B^T A^{-1} b_r
    d^←[d1d2]\hat{d} \leftarrow \begin{bmatrix} d_1 \\ d_2 \end{bmatrix}
    For i=1i = 1 to nn, j=1j = 1 to nn:
        Aˉij←Aij/d^id^j\bar{A}_{ij} \leftarrow A_{ij} / \sqrt{\hat{d}_i \hat{d}_j}
    For i=1i = 1 to nn, j=1j = 1 to mm:
        Bˉij←Bij/d^id^j+n\bar{B}_{ij} \leftarrow B_{ij} / \sqrt{\hat{d}_i \hat{d}_{j+n}}
    Asi←(Aˉ)−1/2A_{\text{si}} \leftarrow (\bar{A})^{-1/2}
    Q←Aˉ+AsiBˉBˉTAsiQ \leftarrow \bar{A} + A_{\text{si}} \bar{B} \bar{B}^T A_{\text{si}}
    Compute eigendecomposition Q=USΛSUSTQ = U_S \Lambda_S U_S^T
    V←[AˉBˉT]AsiUSΛS−1/2V \leftarrow \begin{bmatrix} \bar{A} \\ \bar{B}^T \end{bmatrix} A_{\text{si}} U_S \Lambda_S^{-1/2}
    Sort the columns of VV in ascending order according to the eigenvalues of I−ΛSI - \Lambda_S
    For i=1i = 1 to NEN_E, j=1j = 1 to NN:
        Eij←Vj,i+1/Vj,1E_{ij} \leftarrow V_{j, i+1} / V_{j, 1}
    Cluster the NN column vectors of E∈RNE×NE \in \mathbb{R}^{N_E \times N} into kk groups using kk-means
    return cluster assignments cc

    The algorithm operates in O(n3+n2N)O(n^3 + n^2 N) time, allowing pairwise spectral grouping to scale linearly with the total element count NN.

  5. Knowl 5 — Positive Definiteness of the Exponential Chi-Square Histogram Kernel

    theoretical result

    Let hi,hj∈RKh_i, h_j \in \mathbb{R}^K be discrete probability histograms with positive bins (hi(k)>0h_i(k) > 0 for all ii and k∈{1,…,K}k \in \{1, \dots, K\}). The χ2\chi^2 distance between hih_i and hjh_j is given by:

    χij2=12∑k=1K(hi(k)−hj(k))2hi(k)+hj(k)=1−2∑k=1Khi(k)hj(k)hi(k)+hj(k)\chi_{ij}^2 = \frac{1}{2} \sum_{k=1}^K \frac{(h_i(k) - h_j(k))^2}{h_i(k) + h_j(k)} = 1 - 2 \sum_{k=1}^K \frac{h_i(k) h_j(k)}{h_i(k) + h_j(k)}

    The affinity kernel Wij=exp⁡(−χij2/σ)W_{ij} = \exp(-\chi_{ij}^2 / \sigma) with scale σ>0\sigma > 0 is strictly positive definite. This holds because the matrix QQ with entries Qij=2∑k=1Khi(k)hj(k)hi(k)+hj(k)Q_{ij} = 2 \sum_{k=1}^K \frac{h_i(k) h_j(k)}{h_i(k) + h_j(k)} is positive definite, as demonstrated by the integral identity 1a+b=∫01xa+b−1dx\frac{1}{a+b} = \int_0^1 x^{a+b-1} dx, which yields the quadratic form:

    cTQc=2∑k=1K∫01(∑i=1ncihi(k)xhi(k)−1/2)2dx>0c^T Q c = 2 \sum_{k=1}^K \int_0^1 \left( \sum_{i=1}^n c_i h_i(k) x^{h_i(k) - 1/2} \right)^2 dx > 0

    for any non-zero real vector cc. Since Wij=exp⁡(−1/σ)exp⁡(Qij/σ)W_{ij} = \exp(-1/\sigma) \exp(Q_{ij}/\sigma) and the exponential of a positive definite matrix is positive definite, WW is positive definite. This guarantees that affinity submatrices AA generated by Gaussian-weighted χ2\chi^2 histogram comparisons are positive definite, ensuring the one-shot Nyström solution is applicable.

  6. Knowl 6 — Two-Step Orthonormalization for Indefinite Sample Affinity Submatrices

    model/method

    When the sample affinity submatrix A∈Rn×nA \in \mathbb{R}^{n \times n} is indefinite, the symmetric square root A1/2A^{1/2} is undefined, preventing the one-shot method from being used. Instead, orthonormal eigenvectors of W^=[ABBTBTA−1B]\hat{W} = \begin{bmatrix} A & B \\ B^T & B^T A^{-1} B \end{bmatrix} are computed in two steps:

    1. Diagonalize the sample submatrix A=UAΛAUATA = U_A \Lambda_A U_A^T, and evaluate the standard Nyström extension to form the non-orthogonal approximate eigenvector matrix U∈RN×nU \in \mathbb{R}^{N \times n}:

    U=[UABTUAΛA−1]U = \begin{bmatrix} U_A \\ B^T U_A \Lambda_A^{-1} \end{bmatrix}

    1. Define Z=UΛA1/2∈RN×nZ = U \Lambda_A^{1/2} \in \mathbb{R}^{N \times n} such that W^=ZZT\hat{W} = Z Z^T. Compute the n×nn \times n inner product matrix ZTZZ^T Z and its eigendecomposition ZTZ=FΣFTZ^T Z = F \Sigma F^T, where F∈Rn×nF \in \mathbb{R}^{n \times n} is orthogonal and Σ∈Rn×n\Sigma \in \mathbb{R}^{n \times n} contains the eigenvalues.

    2. The leading orthonormal eigenvectors V∈RN×nV \in \mathbb{R}^{N \times n} and eigenvalues Λ∈Rn×n\Lambda \in \mathbb{R}^{n \times n} of W^\hat{W} are:

    V=ZFΣ−1/2,Λ=ΣV = Z F \Sigma^{-1/2}, \quad \Lambda = \Sigma

    which satisfy VTV=InV^T V = I_n and W^=VΣVT\hat{W} = V \Sigma V^T. This two-step process requires an additional O(n3)O(n^3) operations and introduces greater loss of numerical precision than the one-shot method.

  7. Knowl 7 — Computational Complexity of Nyström Spectral Partitioning

    theoretical result

    For an image or dataset of NN pixels and nn randomly selected sample points (n≪Nn \ll N):

    1. Affinity computations: Evaluates only nN−n(n−1)/2n N - n(n-1)/2 entries (blocks AA and BB) instead of all N(N−1)/2N(N-1)/2 entries of the dense affinity matrix, reducing storage requirements from O(N2)O(N^2) to O(nN)O(n N).
    2. Implicit degree vector computation: Requires O(n2+n(N−n))=O(nN)O(n^2 + n(N-n)) = O(n N) operations.
    3. Inner matrix diagonalization (SS or ZTZZ^T Z): Requires O(n3)O(n^3) operations.
    4. Eigenvector extrapolation: Forming the full N×nN \times n matrix VV takes O(n2N)O(n^2 N) operations.

    The overall computational complexity is O(n3+n2N)O(n^3 + n^2 N), which scales linearly with image size NN for a fixed sample size nn. In comparison, exact dense eigensolvers require O(N3)O(N^3) operations and O(N2)O(N^2) memory, while sparse iterative solvers (e.g., Lanczos) require O(niter⋅n⋅N)O(n_{\text{iter}} \cdot n \cdot N) operations, where the iteration count nitern_{\text{iter}} to convergence varies widely and can degrade significantly on ill-conditioned affinity matrices.

  8. Knowl 8 — Eigenvector Subspace Repeatability and Sample Efficiency on Natural Images

    empirical result

    The stability and sample efficiency of the Nyström Normalized Cut approximation was evaluated on 300 natural scene images from the Corel dataset (image size 240×160=38,400240 \times 160 = 38,400 pixels). Repeatability between two sets of leading NE=4N_E = 4 eigenvectors U,V∈RN×NEU, V \in \mathbb{R}^{N \times N_E} computed from different, independent random sample subsets was evaluated using the normalized subspace projection metric:

    1NE∥UTV∥F2=1−12NE∥UUT−VVT∥F2\frac{1}{N_E} \|U^T V\|_F^2 = 1 - \frac{1}{2 N_E} \|U U^T - V V^T\|_F^2

    which is invariant to internal orthogonal rotations of the eigensubspaces and yields 1.01.0 for identical spans.

    Across 190 unique pairwise comparisons over 20 random trials per image, the agreement metric increased monotonically from ≈0.72\approx 0.72 at n=20n = 20 samples to >0.93> 0.93 at n=400n = 400 samples. This demonstrates that drawing n=100n = 100 to 400400 random samples—less than 1%1\% of total image pixels—is sufficient to capture the leading eigensubspace of natural image affinity matrices with high numerical repeatability.

  9. Knowl 9 — Spatiotemporal Video Segmentation via Space-Time Nyström Embedding

    model/method

    Spatiotemporal video segmentation represents an entire video sequence of dimensions X×Y×TX \times Y \times T as a single three-dimensional space-time volume graph of N=X⋅Y⋅TN = X \cdot Y \cdot T voxels. Pairwise affinities combine local appearance cues (brightness, color, texture) and motion cues across space and time.

    By drawing nn samples across the sequence (e.g., n=100n = 100 voxels drawn uniformly from the first, middle, and last frames), the Nyström approximation computes only an n×Nn \times N affinity matrix and obtains the leading kk eigenvectors of the normalized Laplacian in O(n3+n2N)O(n^3 + n^2 N) time.

    Clustering the leading eigenvector coordinates with kk-means segments the full 3D volume simultaneously. This automatically establishes segment correspondence across video frames without requiring separate frame-to-frame optical flow tracking, motion layer models, or explicit boundary initialization. For a video sequence of 120×120×5120 \times 120 \times 5 voxels (N=72,000N = 72,000), segmentation with n=100n = 100 completes in under one minute.

Coverage note — Specific qualitative image segmentation outputs (tiger image, flowerbed/garden video, jumping leap video) and standard texton/halftoning filterbank construction details from cited prior work were omitted in favor of the core mathematical principles, algorithms, complexity analyses, and quantitative findings.

References

  1. 1.D. Achlioptas, F. McSherry, and B. Schölkopf, ‘Sampling Techniques for Kernel Methods,’ Proc. Neural Information Processing Systems Conf., pp. 335-342, 2002.
  2. 2.E.H. Adelson and J.R. Bergen, ‘Spatiotemporal Energy Models for the Perception of Motion,’ J. Optical Society Am. A, vol. 2, no. 2, pp. 284-299, 1985.
  3. 3.C.T.H. Baker, The Numerical Treatment of Integral Equations. Oxford: Clarendon Press, 1977.
  4. 4.S. Belongie, C. Fowlkes, F. Chung, and J. Malik, ‘Spectral Partitioning with Indefinite Kernels Using the Nyström Extension,’ Proc. European Conf. Computer Vision, 2002.
  5. 5.C. Berg, J.P.R. Christensen, and P. Ressel, Harmonic Analysis on Semigroups. Springer-Verlag, 1984.
  6. 6.C.M. Bishop, Neural Networks for Pattern Recognition. Oxford Univ. Press, 1995.
  7. 7.R.C. Bolles, H.H. Baker, and D.H. Marimont, ‘Epipolar-Plane Image Analysis: An Approach to Determining Structure from Motion,’ Int’l. J. Computer Vision, vol. 1, pp. 7-55, 1987.
  8. 8.K. Boyer, D. Fagerström, M. Kubovy, P. Johansen, and S. Sarkar, ‘POCV99 Breakout Session Report: Spatiotemporal Grouping,’ Perceptual Organization for Artificial Vision Systems, S. Sarkar and K.L. Boyer, eds., 2000.
  9. 9.J.M. Buhmann, ‘Data Clustering and Learning,’ The Handbook of Brain Theory and Neural Networks, M.A. Arbib, ed., pp. 278-281, 1995.
  10. 10.O. Chapelle, P. Haffner, and V. Vapnik, ‘SVMs for Histogram Based Image Classification,’ IEEE Trans. Neural Networks, vol. 10, no. 5, pp. 1055-1064, Sept. 1999.
  11. 11.F.R.K. Chung, Spectral Graph Theory. Am. Math. Soc., 1997.
  12. 12.C. Fowlkes, S. Belongie, and J. Malik, ‘Efficient Spatiotemporal Grouping Using the Nyström Method,’ Proc. IEEE Conf. Computer Vision and Pattern Recognition, Dec. 2001.
  13. 13.A. Frieze and R. Kannan, ‘Quick Approximation to Matricies and Applications,’ Combinatorica, vol. 19, pp. 175-220, 1999.
  14. 14.A. Frieze, R. Kannan, and S. Vempala, ‘Fast Monte-Carlo Algorithms for Finding Low-Rank Approximations,’ Proc. IEEE Symp. Foundations of Computer Science, pp. 370-378, 1998.
  15. 15.Y. Gdalyahu, D. Weinshall, and M. Werman, ‘Stochastic Image Segmentation by Typical Cuts,’ Proc. Conf. Computer Vision and Pattern Recognition, 1999.
  16. 16.S. Gepshtein and M. Kubovy, ‘The emergence of visual objects in space-time,’ Nat’l Academy of Science USA, vol. 97, no. 14, pp. 8186-8191, 2000.
  17. 17.D. Haussler, ‘Convolution Kernels on Discrete Structure,’ technical report, Univ. of California at Santa Cruz, 1999.
  18. 18.J. Malik, S. Belongie, T. Leung, and J. Shi, ‘Contour and Texture Analysis for Image Segmentation,’ Int’l. J. Computer Vision, vol. 43, no. 1, pp. 7-27, June 2001.
  19. 19.S. Mika, B. Schölkopf, A.J. Smola, K.-R. Müller, M. Scholz, and G. Rätsch, ‘Kernel PCA and De-Noising in Feature Spaces,’ Advances in Neural Information Processing Systems 11, pp. 536-542, 1999.
  20. 20.A.Y. Ng, M.I. Jordan, and Y. Weiss, ‘On Spectral Clustering: Analysis and an Algorithm,’ Proc. Neural Information Processing Systems Conf., 2002.
  21. 21.E.J. Nyström, ‘Über die Praktische Auflösung von Linearen Integralgleichungen mit Anwendungen auf Randwertaufgaben der Potentialtheorie,’ Commentationes Physico-Mathematicae, vol. 4, no. 15, pp. 1-52, 1928.
  22. 22.P. Perona and W.T. Freeman, ‘A Factorization Approach to Grouping,’ Proc. Fifth European Conf. Computer Vision, 1998.
  23. 23.W.H. Press, S.A. Teukolsky, W.T. Vetterling, and B.P. Flannery, Numerical Recipies in C, second ed. Cambridge Univ. Press, 1992.
  24. 24.J. Puzicha and S. Belongie, ‘Model-Based Halftoning for Color Image Segmentation,’ Proc. Int’l Conf. Pattern Recognition, vol. 3, pp. 629-632, 2000.
  25. 25.J. Puzicha, T. Hofmann, and J. Buhmann, ‘Non-Parametric Similarity Measures for Unsupervised Texture Segmentation and Image Retrieval,’ Computer Vision and Pattern Recognition, 1997.
  26. 26.S. Sarkar and K.L. Boyer, ‘Quantitative Measures of Change Based on Feature Organization: Eigenvalues and Eigenvectors,’ Proc. Conf. Computer Vision and Pattern Recognition, 1996.
  27. 27.B. Schölkopf, A. Smola, and K.-R. Müller, ‘Nonlinear Component Analysis as a Kernel Eigenvalue Problem,’ Neural Computation, vol. 10, pp. 1299-1319, 1998.
  28. 28.G.L. Scott and H.C. Longuet-Higgins, ‘An Algorithm for Associating the Features of Two Images,’ Proc. Royal Soc. London, vol. B-244, pp. 21-26, 1991.
  29. 29.J. Shi and J. Malik, ‘Motion Segmentation and Tracking Using Normalized Cuts,’ Proc. Int’l Conf. Computer Vision, Jan. 1998.
  30. 30.J. Shi and J. Malik, ‘Normalized Cuts and Image Segmentation,’ IEEE Trans. Pattern Analysis and Machine Intelligence, vol. 22, no. 8, pp. 888-905, Aug. 2000.
  31. 31.Y. Weiss, ‘Smoothness in Layers: Motion Segmentation Using Nonparametric Mixture Estimation,’ Proc. IEEE Conf. Computer Vision and Pattern Recognition, pp. 520-526, 1997.
  32. 32.Y. Weiss, ‘Segmentation Using Eigenvectors: A Unifying View,’ Proc. Seventh Int’l. Conf. Computer Vision, pp. 975-982, 1999.
  33. 33.Y. Weiss and E. Adelson, ‘A Unified Mixture Framework for Motion Segmentation: Incorporating Spatial Coherence and Estimating the Number of Models,’ Proc Conf. Computer Vision and Pattern Recognition, pp. 321-326, June 1996.
  34. 34.C. Williams and M. Seeger, ‘Using the Nyström Method to Speed Up Kernel Machines,’ Advances in Neural Information Processing Systems 13: Proc. 2000 Conf., T.K. Leen, T.G. Dietterich, and V. Tresp, eds., pp. 682-688, 2001.

Citation

MLA
Fowlkes, C., et al. “Spectral Grouping Using the Nystrom Method”. IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 26, no. 2, 2004, pp. 214–25, https://doi.org/10.1109/TPAMI.2004.1262185.
APA
Fowlkes, C., Belongie, S., Fan Chung, & Malik, J. (2004). Spectral grouping using the nystrom method. IEEE Transactions on Pattern Analysis and Machine Intelligence, 26(2), 214–225. https://doi.org/10.1109/TPAMI.2004.1262185
Chicago
Fowlkes, C., S. Belongie, Fan Chung, and J. Malik. 2004. “Spectral Grouping Using the Nystrom Method”. IEEE Transactions on Pattern Analysis and Machine Intelligence 26 (2): 214–25. https://doi.org/10.1109/TPAMI.2004.1262185.
Harvard
Fowlkes, C. et al. (2004) “Spectral grouping using the nystrom method”, IEEE Transactions on Pattern Analysis and Machine Intelligence, 26(2), pp. 214–225. Available at: https://doi.org/10.1109/TPAMI.2004.1262185.
Vancouver
1. Fowlkes C, Belongie S, Fan Chung, Malik J (2004) Spectral grouping using the nystrom method. IEEE Transactions on Pattern Analysis and Machine Intelligence 26:214–225

BibTeX

@article{Fowlkes_2004, title={Spectral grouping using the nystrom method}, volume={26}, ISSN={0162-8828}, url={http://dx.doi.org/10.1109/TPAMI.2004.1262185}, DOI={10.1109/tpami.2004.1262185}, number={2}, journal={IEEE Transactions on Pattern Analysis and Machine Intelligence}, publisher={Institute of Electrical and Electronics Engineers (IEEE)}, author={Fowlkes, C. and Belongie, S. and Fan Chung and Malik, J.}, year={2004}, month=Feb, pages={214–225} }
Metadata:Crossref

Access the Paper

This paper is available from its original source. Click below to access the PDF.

Open PDF