前言

       本文将继续讲解K-近邻算法的项目实例 - 手写识别系统。

       该系统在获取用户的手写输入后，判断用户写的是什么。

       为了突出核心，简化细节，本示例系统中的输入为32x32矩阵，分类结果也均为数字。但对于汉字或者别的分类情形原理都是一样的。

       有了前面学习的基础，下面直接进入项目开发步骤。

第一步：收集并准备数据

       在用户主目录的trainingDigits子目录中，存放的是2000个样本数据。

       每个样本一个文件，其中一部分如下所示：

      

       文件命名格式为：

              分类标签_标签内序号

       如 0_20.txt 就表示该样本是分类标签为0的第20个特征集。20就是个序号以区分标签内不同文件而已，没其他意义。

       样本数据都是32x32矩阵：

      

       对于这样的二维数据，如何判断样本和目标对象的距离呢？首先想到的是可以将二维降到一维。

       当然也可以考虑去找找二维的距离求解方法。

       下面给出降维函数：

1 # ============================================== 2 # 输入： 3 #        训练集文件名(含路径) 4 # 输出： 5 #        降维后的样本数据(这里一个文件一份样本数据) 6 # ============================================== 7 def img2vector(filename): 8     '将32x32的矩阵转换为1024一维向量' 9     10     # 初始化返回向量11     returnVect = numpy.zeros((1,1024))12     13     # 打开样本数据文件14     fr = open(filename)15     16     # 降维处理17     for i in range(32):18         lineStr = fr.readline()19         for j in range(32):20             returnVect[0,32*i+j] = int(lineStr[j])21             22     return returnVect

第二步：测试算法

       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集群下实现这个项目或使用该算法的其他类似项目，定能大幅度提升性能。