//计算区域删减后总的保留区域的方差

         d[0]=0;

         for(k=i-m;k<i;k++)

               for(l=j+k-i;l<j;l++)

               {d[0]+=(AA[k*width+l]-u)*(AA[k*width+l]-u);}

               d[0]=d[0]/(2*m+1)/(2*n+1)*8;

          ////////////////////////////

          d[1]=0;

          for(k=i-m;k<i;k++)

              for(l=j-k+i;l<j;l++)

              {d[1]+=(AA[k*width+l]-u)*(AA[k*width+l]-u);}

              d[1]=d[1]/(2*m+1)/(2*n+1)*8;

           //////////////////////////////

          d[2]=0;

          for(k=i-m;k<i;k++)

              for(l=j-k+i;l<j+m;l++)

              {d[2]+=(AA[k*width+l]-u)*(AA[k*width+l]-u);}

              d[2]=d[2]/(2*m+1)/(2*n+1)*8;

          //////////////////////////////

        d[3]=0;

        for(k=i;k<i+m;k++)

              for(l=j+k-i;l<j+m;l++)

              {d[3]+=(AA[k*width+l]-u)*(AA[k*width+l]-u);}

              d[3]=d[3]/(2*m+1)/(2*n+1)*8;

         //////////////////////////////

        d[4]=0;

        for(k=i;k<i+m;k++)

              for(l=j;l<j+k-i;l++)

              {d[4]+=(AA[k*width+l]-u)*(AA[k*width+l]-u);}

             d[4]=d[4]/(2*m+1)/(2*n+1)*8;

        //////////////////////////////

        d[5]=0;

        for(k=i;k<i+m;k++)

              for(l=j-k+i;l<j;l++)

              {d[5]+=(AA[k*width+l]-u)*(AA[k*width+l]-u);}

               d[5]=d[5]/(2*m+1)/(2*n+1)*8;

        //////////////////////////////

        d[6]=0;

        for(k=i;k<i+m;k++)

              for(l=j-m;l<j-k+i;l++)