kappa系数,如何计算每一分类的kappa和总体kappa

kappa系数,如何计算每一分类的kappa和总体kappa

Kappa系数：是一种计算分类精度的方法。它是通过把所有地表真实分类中的像元总数（N）乘以混淆矩阵对角线（Xkk）的和，再减去某一类地表真实像元总数与被误分成该类像元总数之积对所有类别求和的结果，再除以总像元数的平方差减去某一类中地表真实像元总数与该类中被分类像元总数之积对所有类别求和的结果所得到的。
两幅栅格图的kappa计算公式为 k = (Po-Pc)/(1-Pc)
设栅格总象元数为n，真实栅格为1的象元数为a1，为0的象元数为a0，模拟栅格为1的象元数为b1，为0的象元数为b0，两个栅格对应象元值相等的象元数为s，则
Po = s/n， Pc = (a1*b1+a0*b0)/(n*n)
在arcmap里，a1、a0、b1、b0从属性表可读出，s用raster calculator配合con()函数不难求出。
kappa计算结果为-1~1，但通常kappa是落在 0~1 间，可分为五组来表示不同级别的一致性：0.0~0.20极低的一致性(slight)、0.21~0.40一般的一致性(fair)、0.41~0.60 中等的一致性(moderate)、0.61~0.80 高度的一致性(substantial)和0.81~1几乎完全一致(almost perfect)。追问请不要复制黏贴，我要的是分项kappa