kappa系数的计算公式

两幅栅格图的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 calculat