SVD
•Scatter matrix can be big, so computation non-trivial.
•Stack data into matrix X, each row an image.  SVD gives X = UDVT
–D is diagonal with non-increasing values.
–U and V have orthonormal rows. 
•VT(:,1:k) gives first k principal components.
•matlab
–
•