fastPCA 의 python 구현

```python
'''
fastPCA python  
'''
import numpy as np
import heapq
def fastPCA(A, k):
    r, c = A.shape
    meanVec = np.mat(np.mean(A, axis=0))
    Z = (A - np.tile(meanVec, (r, 1)))
    covMatT = np.matmul(Z, Z.T)
    D, V = np.linalg.eig(covMatT)
    # D = np.mat(D)
    temp = np.unique(heapq.nlargest(k, D))
    indx = np.argwhere(D >= temp[0]).T.flatten()
    D = D[indx]
    D = np.diag(D)
    V = V[:, indx]
    V = Z.T * V
    for i in range(k):
        V[:, i] = V[:, i]/np.linalg.norm(V[:, i])
    pcaA = Z * V
    return pcaA, V

'디지털 이미지 처리 와 기계 시각' 의 'PCA 기반 얼굴 특징 추출' 을 참고 한다.