Calculate an uncorrelated multilinear principal component basis representation for functional data on two-dimensional domains