前言
本文将继续讲解K-近邻算法的项目实例 - 手写识别系统。
该系统在获取用户的手写输入后,判断用户写的是什么。
为了突出核心,简化细节,本示例系统中的输入为32x32矩阵,分类结果也均为数字。但对于汉字或者别的分类情形原理都是一样的。
有了前面学习的基础,下面直接进入项目开发步骤。
第一步:收集并准备数据
在用户主目录的trainingDigits子目录中,存放的是2000个样本数据。
每个样本一个文件,其中一部分如下所示:
文件命名格式为:
分类标签_标签内序号
如 0_20.txt 就表示该样本是分类标签为0的第20个特征集。20就是个序号以区分标签内不同文件而已,没其他意义。
样本数据都是32x32矩阵:
对于这样的二维数据,如何判断样本和目标对象的距离呢?首先想到的是可以将二维降到一维。
当然也可以考虑去找找二维的距离求解方法。
下面给出降维函数:
# ==============================================
# 输入:
# 训练集文件名(含路径)
# 输出:
# 降维后的样本数据(这里一个文件一份样本数据)
# ==============================================
def img2vector(filename):
'将32x32的矩阵转换为1024一维向量' # 初始化返回向量
returnVect = numpy.zeros((1,1024)) # 打开样本数据文件
fr = open(filename) # 降维处理
for i in range(32):
lineStr = fr.readline()
for j in range(32):
returnVect[0,32*i+j] = int(lineStr[j]) return returnVect
第二步:测试算法
K临近的分类函数代码在之前的文章K-近邻分类算法原理分析与代码实现中给出了,这里直接调用:
# =================================================
# 输入:
# 空
# 输出:
# 对指定的测试集文件,指定的训练集数据进行K近邻分类
# 并打印结果信息
# =================================================
def handwritingClassTest():
'手写数字识别系统测试代码' # 分类列表
hwLabels = [] # 获取所有训练集文件名
trainingFileList = os.listdir('/home/fangmeng/trainingDigits') # 定义训练集结构体
m = len(trainingFileList)
trainingMat = numpy.zeros((m, 1024)) for i in range(m):
# 当前训练集文件名
filenameStr = trainingFileList[i]
# 文件名(filenameStr去掉.txt后缀)
fileStr = filenameStr.split('.')[0]
# 分类标签
classNumStr = int(fileStr.split('_')[0])
# 将分类标签加入分类列表
hwLabels.append(classNumStr)
# 将当前训练集文件降维后加入到训练集结构体
trainingMat[i] = img2vector('/home/fangmeng/trainingDigits/%s' % filenameStr) # 获取所有测试集文件名
testFileList = os.listdir('/home/fangmeng/testDigits')
# 错误分类记数
errorCount = 0
# 测试集文件个数
mTest = len(testFileList) print "错误的分类结果如下:"
for i in range(mTest):
# 当前测试集文件名
fileNameStr = testFileList[i]
# 文件名(filenameStr去掉.txt后缀)
fileStr = fileNameStr.split('.')[0]
# 分类标签
classNumStr = int(fileStr.split('_')[0])
# 将当前测试集文件降维
vectorUnderTest = img2vector('/home/fangmeng/testDigits/%s' % fileNameStr)
# 对当前测试文件进行分类
classifierResult = classify0(vectorUnderTest, trainingMat, hwLabels, 3) if (classifierResult != classNumStr):
print "分类结果: %d, 实际结果: %d" % (classifierResult, classNumStr)
errorCount += 1.0 print "\n总错误数: %d" % errorCount
print "\n总错误数: %f" % (errorCount/float(mTest))
运行结果:
小结
1. K-邻近算法的本质是用来分类的,要从分类的思想去思考这个算法的运用。
2. 再强调一次K-邻近算法是没有训练过程的,这点和以后学习的其他分类方法,比如决策树对比后就更清楚了。
3. K-邻近算法的效率很低,不论是从时间还是空间上看(单就这个简单项目都跑得很慢)。因此需要学习更多更优化的算法。
4. 有兴趣有时间可以考虑在hadoop/spark集群下实现这个项目或使用该算法的其他类似项目,定能大幅度提升性能。