0

背景：整个程序设计为在参考图像和1800个目标图像之间进行2D DIC，（用于断层重建）在我的代码中，有这个for循环块

for (k=0; k<kmax; k++)
{
    K=nm12+(k*(h-n+1))/(kmax-1);
    printf("\nk=%d\nL= ", K);
    for (l=0; l<lmax; l++)
    {
        ///For each subset, calculate and store its mean and standard deviation.   
        ///Also want to know the sum and sum of squares of subset, but in two sections, stored in fm/df[k][l][0 and 1].
        L=nm12+(l*(w-n+1))/(lmax-1);
        printf("%d ", L);
        fm[k][l][0]=0;
        df[k][l][0]=0;
        fm[k][l][1]=0;
        df[k][l][1]=0;
        ///loops are j then i as it is more efficient (saves m-1 recalculations of b=j+L; 
        for (j=0; j<m; j++)
        {
            b=j+L;
            for (i=0; i<M; i++)
            {
                a=i+K;
                fm[k][l][0]+=ref[a][b];
                df[k][l][0]+=ref[a][b]*ref[a][b];
            }

            for (i=M; i<m; i++)
            {
                a=i+K;
                fm[k][l][1]+=ref[a][b];
                df[k][l][1]+=ref[a][b]*ref[a][b];
            }
        }
        fm[k][l][2]=m2r*(fm[k][l][1]+fm[k][l][0]);
        df[k][l][2]=sqrt(df[k][l][1]+df[k][l][0]-m2*fm[k][l][2]*fm[k][l][2]);
        a+=1;
    }
}

每次 l 达到 10 时，该行df[k][l][2]=sqrt(df[k][l][1]+df[k][l][0]-m2*fm[k][l][2]*fm[k][l][2]);似乎不再被执行。我的意思是调试器显示 df[k][l][2] 的值没有正确地从零更改为总和。此外，只要 l>=10，df[k][l][0 和 1] 就保持固定，而不管 k 和 l。

kmax=15, lmax=20, n=121, m=21, M=(3*m)/4=15, nm12=(n-m+1)/2=50。

数组 fm 和 df 是双数组，声明double fm[kmax][lmax][3], df[kmax][lmax][3];

此外，该行a+=1;只是用作检查 df[k][l][2] 值的断点，并且对代码功能没有影响。

任何关于为什么会发生这种情况、如何​​修复等的帮助将不胜感激！

编辑：更多信息。

数组 ref（包含参考图像像素值）是一个动态数组，使用 malloc 分配内存，在此代码块中：

    double **dark, **flat, **ref, **target, **target2, ***gm, ***dg;
    dark=(double**)malloc(h * sizeof(double*));
    flat=(double**)malloc(h * sizeof(double*));
    ref=(double**)malloc(h * sizeof(double*));
    target=(double**)malloc(h * sizeof(double*));
    target2=(double**)malloc(h * sizeof(double*));
    size_t wd=w*sizeof(double);
    for (a=0; a<h; a++)
    {
        dark[a]=(double*)malloc(wd);
        flat[a]=(double*)malloc(wd);
        ref[a]=(double*)malloc(wd);
        target[a]=(double*)malloc(wd);
        target2[a]=(double*)malloc(wd);
    }

其中 h=1040 和 w=1388 是图像的尺寸。

4

1 回答 1

2

您没有过多提及您正在使用的编译器、IDE 或框架。但隔离问题的一种方法是创建一个新的小型（控制台）项目，其中仅包含您发布的代码段。这样，您将消除大多数类型的输入/线程/堆栈/内存/编译器等问题。
如果没有，它足够小，可以将整个示例发布到 stackoverflow 上，供我们拆开思考。

因此，您应该为您的算法创建一个独立的单元测试。

于 2013-08-14T15:07:58.173 回答