To create this morphing sequence, I first defined the correspondence points between the two images, then I computed the delaunay triangulation of the intermediate face and solving a 6x6 linear system for each triangle I found the affine matrices that transform the triangulation of the final face to the triangulation of the face to be morphed. Finally I multiplied each pixel of the target image by the correspondent affine matrix and intepolated the result with interp2.
To compute the mean face I first calculated the average of the markers position, then I morphed all the faces into it and finally I took the average color.
