During morphogenesis, epithelia (cell sheets) undergo complex deformations as they stretch, bend, and twist to form the embryo. Often these changes in shape create multivalued surfaces that can be problematic for strain measurements. This paper presents a method for quantifying deformation of such surfaces. The method requires four-dimensional spatiotemporal coordinates of a finite number of surface markers, acquired using standard imaging techniques. From the coordinates of the markers, various deformation measures are computed as functions of time and space using straightforward matrix algebra. This method accommodates sparse randomly scattered marker arrays, with reasonable errors in marker locations. The accuracy of the method is examined for some sample problems with exact solutions. Then, the utility of the method is illustrated by using it to measure surface stretch ratios and shear in the looping heart and developing brain of the early chick embryo. In these examples, microspheres are tracked using optical coherence tomography. This technique provides a new tool that can be used in studies of the mechanics of morphogenesis.