Add simple PCA flat field util.

Here we directly call ARPACK SVD solver (through scipy) and hence to not need to rely on scikit-learn.