Una herramienta fundamental en el análisis de formas es la incrustación virtual de la variedad riemanniana que describe la geometría de una forma en el espacio euclidiano. Se han propuesto varios métodos para incrustar formas isométricas en dominios planos, manteniendo las distancias medidas en la variedad. Recientemente, se ha prestado atención a la incrustación de formas en el espacio propio del operador Laplace-Beltrami. El espacio propio de Laplace-Beltrami preserva la distancia de difusión y es invariante bajo transformaciones isométricas. Sin embargo, las funciones propias de Laplace-Beltrami calculadas de forma independiente para diferentes formas suelen ser incompatibles entre sí. Las aplicaciones que involucran múltiples formas, como la correspondencia punto a punto, se beneficiarían enormemente si sus respectivas funciones propias se emparejaran de alguna manera. Aquí, presentamos un enfoque estadístico para emparejar las funciones propias. Consideramos los valores de las funciones propias sobre la variedad como la muestra de variables aleatorias e intentamos emparejar sus distribuciones multivariadas. La comparación de distribuciones se realiza indirectamente, utilizando estadísticas de alto orden. Mostramos que las ambigüedades de permutación y signo de las funciones propias de bajo orden pueden inferirse minimizando la diferencia de sus momentos de tercer orden. Las ambigüedades de signo de las funciones propias antisimétricas pueden resolverse explotando las relaciones invariantes isométricas entre los gradientes de las funciones propias y la normal de la superficie. Presentamos experimentos que demuestran el éxito del método propuesto aplicado a la correspondencia de puntos característicos.