C++相关系数法影像匹配代码

以下是基于C语言的影像匹配代码，使用相关系数法：

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

int main()
{
    int i, j, k, l, m, n;
    int row1 = 512, col1 = 512; //第一幅图像的大小
    int row2 = 512, col2 = 512; //第二幅图像的大小
    double **img1, **img2; //图像的二维数组
    double **corr; //相关系数矩阵
    double maxcorr = -1.0; //最大相关系数
    int maxcorr_i = 0, maxcorr_j = 0; //最大相关系数的位置

    //分配内存空间
    img1 = (double **)malloc(row1 * sizeof(double *));
    img2 = (double **)malloc(row2 * sizeof(double *));
    corr = (double **)malloc((row1 - row2 + 1) * sizeof(double *));
    for (i = 0; i < row1; i++)
    {
        img1[i] = (double *)malloc(col1 * sizeof(double));
    }
    for (i = 0; i < row2; i++)
    {
        img2[i] = (double *)malloc(col2 * sizeof(double));
    }
    for (i = 0; i < row1 - row2 + 1; i++)
    {
        corr[i] = (double *)malloc((col1 - col2 + 1) * sizeof(double));
    }

    //读取图像数据
    FILE *fp1 = fopen("img1.raw", "rb");
    FILE *fp2 = fopen("img2.raw", "rb");
    for (i = 0; i < row1; i++)
    {
        for (j = 0; j < col1; j++)
        {
            fread(&img1[i][j], sizeof(double), 1, fp1);
        }
    }
    for (i = 0; i < row2; i++)
    {
        for (j = 0; j < col2; j++)
        {
            fread(&img2[i][j], sizeof(double), 1, fp2);
        }
    }
    fclose(fp1);
    fclose(fp2);

    //计算相关系数
    for (i = 0; i < row1 - row2 + 1; i++)
    {
        for (j = 0; j < col1 - col2 + 1; j++)
        {
            double mean1 = 0.0, mean2 = 0.0, var1 = 0.0, var2 = 0.0, cov = 0.0;
            for (k = 0; k < row2; k++)
            {
                for (l = 0; l < col2; l++)
                {
                    mean1 += img1[i + k][j + l];
                    mean2 += img2[k][l];
                }
            }
            mean1 /= (double)(row2 * col2);
            mean2 /= (double)(row2 * col2);
            for (k = 0; k < row2; k++)
            {
                for (l = 0; l < col2; l++)
                {
                    var1 += pow(img1[i + k][j + l] - mean1, 2.0);
                    var2 += pow(img2[k][l] - mean2, 2.0);
                    cov += (img1[i + k][j + l] - mean1) * (img2[k][l] - mean2);
                }
            }
            var1 /= (double)(row2 * col2 - 1);
            var2 /= (double)(row2 * col2 - 1);
            cov /= (double)(row2 * col2 - 1);
            corr[i][j] = cov / sqrt(var1 * var2);
            if (corr[i][j] > maxcorr)
            {
                maxcorr = corr[i][j];
                maxcorr_i = i;
                maxcorr_j = j;
            }
        }
    }

    //输出结果
    printf("最大相关系数: %lf\n", maxcorr);
    printf("最大相关系数位置: (%d, %d)\n", maxcorr_i, maxcorr_j);

    //释放内存空间
    for (i = 0; i < row1; i++)
    {
        free(img1[i]);
    }
    for (i = 0; i < row2; i++)
    {
        free(img2[i]);
    }
    for (i = 0; i < row1 - row2 + 1; i++)
    {
        free(corr[i]);
    }
    free(img1);
    free(img2);
    free(corr);

    return 0;
}

需要注意的是，这里假设输入的是双精度浮点数类型的图像数据，存储在.raw格式的文件中。如果需要使用其他类型的图像数据，需要进行相应的修改。此外，还需要根据实际情况调整图像的大小和相关系数矩阵的大小。