void srcn(short int** pd,int* nse,int il,int jl) { int i,j,k,l; double std,avg,temp; for(i=rn;i<il-rn;i++){ for(j=rn;j<jl-rn;j++){ avg=0; temp=0; std=0; for(k=-rn;k<=rn;k++){ for(l=-rn;l<=rn;l++){ avg+=(double)pd[i+k][j+l]; } } avg=avg/((double)sn*sn); for(k=-rn;k<=rn;k++){ for(l=-rn;l<=rn;l++){ temp+=((double)pd[i+k][j+l]-avg)*((double)pd[i+k][j+l]-avg); } } std=pow(temp/((double)sn*sn),0.5); if(std<c*drt){ nse[0]++; } else{ nse[1]++; } } } }
var
This content, along with any associated source code and files, is licensed under The Code Project Open License (CPOL)