For the spatial kernel k Y , we use a Gaussian RBF, which is a widely used choice for vector data.
we define k Y (y i , y j ; γ) = exp(−γd 2 i,j ), where d 2 ij is the squared Euclidean distance between y i and y j , and γ is a parameter of the kernel function.