import cv2
import numpy as np
SZ = 20
bin_n = 16
affine_flags = cv2.WARP_INVERSE_MAP | cv2.INTER_LINEAR
# 使用图片的二阶矩对图片进行抗扭斜处理
def deskew(img):
m = cv2.moments(img) # 获取图片的矩
if abs(m['mu02']) < 1e-2:
return img.copy()
skew = m['mu11'] / m['mu02']
M = np.float32([[1, skew, -0.5 * SZ * skew], [0, 1, 0]])
img = cv2.warpAffine(img, M, (SZ, SZ), flags=affine_flags)
return img
# 计算图像的hog描述符
def hog(img):
gx = cv2.Sobel(img, cv2.CV_32F, 1, 0)
gy = cv2.Sobel(img, cv2.CV_32F, 0, 1)
mag, ang = cv2.cartToPolar(gx, gy)
bins = np.int32(bin_n * ang / (2 * np.pi))
bin_cells = bins[:10, :10], bins[10:, :10], bins[:10, 10:], bins[10:, 10:]
mag_cells = mag[:10, :10], mag[10:, :10], mag[:10, 10:], mag[10:, 10:]
hists = [np.bincount(b.ravel(), m.ravel(), bin_n) for b, m in zip(bin_cells, mag_cells)]
hist = np.hstack(hists)
return hist
# 获取图片,拆分图片
img = cv2.imread('digits.png', 0)
# 将数据分割成50行,100列
cells = [np.hsplit(row, 100) for row in np.vsplit(img, 50)]
train_cells = [i[:50] for i in cells] # 前50列作为训练数据
test_cells = [i[50:] for i in cells] # 后50列作为测试数据
deskewed = [list(map(deskew, row)) for row in train_cells] # 对训练数据进行抗扭斜数据
hogdata = [list(map(hog, row)) for row in deskewed] # 获取数据的hog符
trainData = np.float32(hogdata).reshape(-1, 64) # 行数不知道,转变成64列
labels = np.repeat(np.arange(10), 250)[:, np.newaxis] # 获取标签,
# 创建SVM分类器,设置各项参数
svm = cv2.ml.SVM_create()
svm.setKernel(cv2.ml.SVM_LINEAR)
svm.setType(cv2.ml.SVM_C_SVC)
svm.setC(2.67)
svm.setGamma(5.383)
svm.train(trainData, cv2.ml.ROW_SAMPLE, labels)
svm.save('svm_data.dat')
# 对测试数据进行处理
deskewed = [list(map(deskew, row)) for row in test_cells]
hogdata = [list(map(hog, row)) for row in deskewed]
testData = np.float32(hogdata).reshape(-1, bin_n * 4)
ret, result = svm.predict(testData)
mask = result == labels
correct = np.count_nonzero(mask)
print(correct * 100.0 / len(result), '%')
print(len(result))
print(correct)