最近修改论文已经到了最后的关口，结果被质疑分类精度评估的方法。
我原来的方法是利用ArcGIS建立格网，格网中心点作为我的目视解译点，然后利用ArcGIS收集该点的影像分类结果，再把数据属性导出用python计算混淆矩阵。
被质疑之后我使用了分层抽样法，原理我就不赘述了，通过分层抽样法可以计算出一幅影像每个类别需要的抽样点为多少。
之后在分类影像上每个类别的点随机撒在该类别的位置，再人工目视解译，与影像分类结果一起计算混淆矩阵即可。
我首先尝试了Erdas，首先有着ArcGIS的电脑安装erdas后会出现一定的问题，之后我会再写一篇文章概述。并且最后发现erdas抽样的点很难导入ArcGIS或者有个excel表格这样方便计算混淆矩阵，于是我放弃了。
之后网上搜索到ArcMap 10.4以上版本可以较为方便的，进行分层抽样撒点，但是按照高版本的有风险，我目前用10.2
ArcGIS 10.2没有办法让点能够撒到分类影像上那一类的位置上去，但是我最后发现ENVI可以使用这个功能。
正文如下：----------------------------------------不想看我叙述的可以看这里------------------------------
1.分类
首先我的分类结果是tif影像，1、2、3这些值来表示各个类别，tif不能直接在ENVI中进行抽样撒点，要在ENVI中分类，采用决策树，也就是classification-decision tree，首先先建立一个决策树，注意这是让你选b1的波段，先不要选，之前我选了后来换成别的影像存在错误。保存之后，使用classification-decision tree-execute exiting decision tree，选择之前保存的决策树txt。选择你的分类结果影像，然后就可以得到分类影像。
2.抽样撒点
使用envi的classification-post classification-generate random sample-using ground truth image
之后选择你需要的影像，点确定
之后可以选择类（注意这里我给了一个类别给影像上没有分类的，就是none，可以不选它）
点OK后，选择分层抽样，和diproportionnate，然后点击set class sample sizes，就可以设定每个类别撒多少点了
其他选择不变保存roi
之后打开影像
选择overlay-region of interest，就可以看到生成的点（其实就是roi)，注意这里的ROI可能存在之前生成的，不是你这次的，删掉。以下窗口选择，选择File-output roi to shapefile这个，
如下，选择你需要的类别保存即可，基本是全选
3.shapefile获取影像分类点+目视解译
将shp文件在ArcGIS打开，之后选择spatial analysis tools-extraction-extract values to points,选择你的shp文件和影像，最后可以得到一个矢量文件，里面有一栏RASTERVALU属性就是影像的值。
之后只要打开属性表，add field增加新的一列，作为目视解译的结果，然后右键使用栅格计算器，使它等于RASTERVALU这一列即可（这一步是为了减小工作量，当然也可以自己赋值），然后就可以目视解译改这个值啦。
4.计算混淆矩阵
得到目视解译后的点的矢量文件，可以使用python打开shp文件自带的dbf文件（其实就是属性表）计算混淆矩阵和kappa系数即可
代码如下
from sklearn.metrics import confusion_matrix
from sklearn.metrics import accuracy_score
from sklearn.metrics import cohen_kappa_score
from sklearn.metrics import precision_recall_fscore_support
from sklearn.metrics import classification_report
from sklearn.metrics import precision_score
from sklearn.metrics import recall_score
import pandas
import os
from dbfread import DBF
rootdir = "H:\\result_2_type\\points"
list = os.listdir(rootdir)
colnames = ['OID_', 'RASTERVALU', 'true']
true=[]
result=[]
for i in range(0,len(list)):
if not list[i].endswith(".dbf"):
continue
folder = os.path.join(rootdir,list[i])
print(folder)
table = DBF(folder, encoding='GBK')
df = pandas.DataFrame(iter(table))
result_temp = df['RASTERVALU'].tolist()
true_temp = df['true'].tolist()
score=accuracy_score(true_temp, result_temp)
print('score=',score)
kappa=cohen_kappa_score(true_temp, result_temp)
print('kappa=',kappa)
report=classification_report(true_temp, result_temp,digits=4)
print('report=',report)
true.extend(true_temp)
result.extend(result_temp)