Сравнение изображений

    
    int xc = W/2;
    int yc = H/2;
 
    int r = 500;
    int d = 20;
 
    float corr_max = -1.0;
 
    int i_max = 0; 
    int j_max = 0;
 
    int progress = 0;
 
    float img1_avar = 0.0;
    for(int y = yc-r; y <= yc+r; ++y){
        for(int x = xc-r; x <= xc+r; ++x)
            if(x>=0 && y>=0 && x<W && y<H){
                img1_avar += arr1[x+y*W];
            }
    }
    img1_avar /= pow(2*r+1, 2);
 
    for(int j = yc-d; j <= yc+d; ++j){
        for(int i = xc-d; i <= xc+d; ++i){
            if(i>=0 && j>=0 && i<W && j<H){
                float img2_avar = 0.0;
                for(int y = j-r; y <= j+r; ++y){
                    for(int x = i-r; x <= i+r; ++x){
                        if(x>=0 && y>=0 && x<W && y<H){
                            img1_avar += arr2[x+y*W];  
                        }
                    }
                }
 
                img1_avar /= pow(2*r+1, 2);
 
                float s1 = 0;
                float s2 = 0;
                float s3 = 0;
 
                for(int y = j-r; y <= j+r; ++y){
                    for(int x = i-r; x <= i+r; ++x){
                        if(x>=0 && y>=0 && x<W && y<H &&
                          (x-i+xc)>=0 && (y-j+yc)>=0 && (x-i+xc)<W && (y-j+yc)<H){
                            float a1 = arr1[x-i+xc + (y-j+yc)*W] - img1_avar;
                            float a2 = arr2[x+y*W] - img2_avar;
 
                            s1 += a1*a2;
                            s2 += a1*a1;
                            s3 += a2*a2;
                        }
                    }
                }
 
                float corr = s1/(sqrt(s2)*sqrt(s3));
                if(corr > corr_max){
                    corr_max = corr;
                    i_max = i;
                    j_max = j;
                }
            }
            progress ++;
        }
 
        cout << 100*progress/((2*d+1)*(2*d+1)) << endl;
    }
 
    int dx = xc - i_max;
    int dy = yc - j_max;