I'm trying to read a paper called "Graph Embedding Discriminant Analysis on Grassmannian Manifolds for Improved Image Set Matching" and I came across a sentence that confused me (the last one): 



> 
> A manifold is a topological space that is locally similar to Euclidean space. At an intuitive level, manifolds can be thought of as smooth, curved surfaces embedded in higher dimensional Euclidean spaces. Riemannian manifolds are endowed with a distance measure which allows us to measure how similar two points are. We are interested in a particular class of Riemannian manifolds, known as Grassmannian manifolds. [...] **Grassmannian kernels allow us to treat the Grassmannian space as if it were a Euclidean vector space**. As a result, learning algorithms in vector spaces can be extended to their counterparts on Grassmannian manifolds.
> 
> 
> 


In my understanding, kernels are similarity functions $k(x,y)$ that can be used to map data points $x$ and $y$ to a new more "useful" space (where classes become linearly separable in the case of support vector machines for example), and the mapped points would be compared in that new space instead of the original one (please correct me if I'm wrong).


Now, I dont understand what the concept of kernels has to do with making it possible to treat points on a Grassmannian manifold as if they were in an Euclidean vector space. I see a Grassmannian manifold as a smooth (probably curvy and very complex) high dimensional surface that can be seen locally as an Euclidian space (just like a tiny bit on a 3D sphere where an ant can percieve its surroundings as a flat surface) and a geodesic would be a natural choice for distance computation on this space but I dont understand the use of Kernels in this context.


Edit: Grassmann kernels are defined in section 4 of [this](http://machinelearning.wustl.edu/mlpapers/papers/icml2008_HamL08) paper.


Thanks ! 



How can Kernel functions make a Grassmann manifold into an Euclidean vector space?
If you have you have a set $X$ and a kernel function $k: X \times X \to \mathbb R$, then you can associate to $k$ a Hilbert space $H$ of functions on $X$ in the following manner: set 
$$H\_0 = \operatorname{span} \{ k(\cdot,x) \,:\, x \in X \}\,,$$
and define on $H\_0$ the inner product
$$ \langle \sum\_i \alpha\_i k(\cdot,x\_i), \sum\_j \beta\_j k(\cdot,x\_j) \rangle
= \sum\_{i,j} \alpha\_i \beta\_j k(x\_i,x\_j)\,.$$
Then $H$ is the completion of $H\_0$ with respect to $\langle \cdot, \cdot \rangle$ and is called the *reproducing kernel Hilbert space* associated to $k$. One can identify $H$ with a subspace of the vector space of all functions on $X$.


Now on to Grassmannians: we can identify each element $x \in X$ with the function $k(\cdot, x)$, which is an element of the vector space $H$. In a way you "smear out" the point $x$ via the kernel $k(\cdot,x)$. Thus we can treat elements of the set $X$ as vectors in a vector space.


