### Description

Diffusion smoothing based on the finite element method (FEM) is introduced in Chung et al (2001). The following self-contained package (diffusion_smooth.zip) perform diffusion smoothing. Simply run DEMO.diffusion.m line by line. You should be able to generate Figure 1. The diffusion_smooth.m routine requires FEM.m that computes A and C matrices associated with the finite element discritization of the given surface mesh (Chung & Taylor, 2004).Instead of solving heat diffusion using FEM, it is possible to solve it as a series expansion involving a heat kernel (Seo et al. 2010, Seo & Chung, 2011). Heat kernel smoothing approach is more stable and robust. This is the approach I strongly recommend to use to smooth surface data and surface itself.

### References

[1] Chung, M.K. 2001. Statistical
Morphometry in Neuroanatomy, PhD Thesis,
McGill University.[2] Chung, M.K., Taylor, J. 2004. Diffusion Smoothing on Brain Surface via Finite Element Method, IEEE International Symposium on Biomedical Imaging (ISBI). 562.

[3] Seo, S., Chung, M.K., Vorperian, H. K. 2010. Heat kernel smoothing using Laplace-Beltrami eigenfunctions. 13th International Conference on Medical Image Computing and Computer Assisted Intervention (MICCAI). Lecture Notes in Computer Science (LNCS). 6363:505-512.

[4] Seo, S., Chung, M.K. 2011. Laplace-Beltrami eigenfunction expansion of cortical manifolds. IEEE International Symposium on Biomedical Imaging (ISBI).