PCA OCR Identification

PCA_OCR.py

u, sigma, test = train(read(), 30)
right = 0.0
total = 0.0
for img in test:
    matrix = np.subtract(img[1], np.mean(img[1])) # Center
    matrix = np.reshape(matrix, (-1, 1))
    covImg = np.dot(matrix, np.reshape(matrix, (1, -1))).T

    index = 0;
    best = 0.0

    for i in range(10):
      if best < np.mean(covImg * sigma[i]):
          index = i
          best = np.mean(covImg * sigma[i])
    if img[0] == index:
        right += 1.0
    total += 1.0
    print (right/total), total, img[0], index