遥感图像解译 kappa系数计算程序（C++）

//kappa系数的计算程序 # include <iostream> using namespace std; void main() { double Kappa[100][100],Arr[100],iArr[100],jArr[100]; ////////矩阵数组初始化//////////////// for(int i=1;i<100;i++) { Arr[i]=0; iArr[i]=0; jArr[i]=0; for( int j=1;j<100;j++) { Kappa[i][j]=0; } } /////////////////获得混淆矩阵////////////// int n; cout<<"请输入混淆矩阵的行数（或列数）："; cin>>n; cout<<"请输入混淆矩阵中的元素："<<endl; for( i=1;i<=n;i++) { for(int j=1;j<=n;j++) { cout<<"Kappa["<<i<<"]["<<j<<"]="; cin>>Kappa[i][j]; } } ///////////////////输出混淆矩阵/////////////////////// cout<<"输入的混淆矩阵为："<<endl; for( i=1;i<=n;i++) { for(int j=1;j<=n;j++) { cout<<" "<<Kappa[i][j]<<" "; } cout<<endl; } //////////////////Kappa系数的计算////////////////////////// ///////////像元总数N的计算///////////////// double N=0.0;//像元总数 for( i=1;i<=n;i++) { for( int j=1;j<=n;j++) { N+=Kappa[i][j]; } } //cout<<N<<endl; //////////////对角线上像元之和m1////////////// double m1=0.0; for( i=1;i<=n;i++) { m1+=Kappa[i][i]; } //cout<<m1; /////////////求每一行、每一列的像元和/////////////// for( i=1;i<=n;i++) { for( int j=1;j<=n;j++) { iArr[i]+=Kappa[i][j];//一行和 jArr[i]+=Kappa[j][i];//一列和 } } /*cout<<iArr[1]<<endl; cout<<jArr[1]<<endl; cout<<iArr[2]<<endl; cout<<jArr[2]<<endl;*/ ////////////////////////每行每列对应相乘/////////////////// for( i=1;i<=n;i++) { Arr[i]= iArr[i]*jArr[i]; //cout<<Arr[i]<<endl; } ////////////////////计算乘积结果之和///////////////// double ALL=0.0; for( i=1;i<=n;i++) { ALL+=Arr[i]; } ////////////计算像元总数的平方////////////// double dN; dN=N*N; ///////////最终结果//////////////////// double kappa; kappa=(N*m1-ALL)/(dN-ALL); cout<<"总体精度："<<(m1/N)<<endl; cout<<"生产者精度："<<endl; for( i=1;i<=n;i++) { cout<<"p"<<i<<"="<<Kappa[i][i]/jArr[i]<<endl; } cout<<"使用者精度："<<endl; for( i=1;i<=n;i++) { cout<<"u"<<i<<"="<<Kappa[i][i]/iArr[i]<<endl; } cout<<"kappa="<<kappa<<endl; }