We present two methods based on functional principal component analysis (FPCA) for the estimation of smooth derivatives of a sample of random functions, which are observed in a more than one-dimensional domain.We apply eigenvalue decomposition to a) the dual covariance matrix of the derivatives,...