-3

我对具有 C 和 C++ 基本知识的编程非常陌生。我一直在尝试在用于图像分析的 .dat 文件上运行代码。

基本上，这个 .dat 文件具有簇的 x 和 y 位置以及它们的半径。使用此代码，我试图将我的图像（缩放像素后尺寸为 200 毫米 x 300 毫米）划分为 5 毫米 x 5 毫米的网格，并找到半径 = 1 的点数和半径 = 100 的点数。

这是我的代码：

#include<stdio.h>
int main()
{ 
    int i=0,g,h,ng[350][250],ns[350][250],ntot[350][250];
    float a[6001],b[6001],c[6001];
    FILE *fres1, *fres2;

    for(g=5;g<=350;g=g+5)
    { 
        for(h=5;h<=250;h=h+5)
        {
            ng[g][h]=0;
            ns[g][h]=0;
            ntot[g][h]=0;   
        } 
    }

    fres1 = fopen("trial_final_1.dat","r+"); 
    fres2 = fopen("result_boxes.dat","w");

    if(fres1 == NULL)
    {
        printf("Error!! Cannot open file \n" );
        return 1;
    } 
    else
    {
        printf("File opened successfully, kindly check your folder for result_boxes file :) \n");
        fprintf(fres2,"x(mm)        y(mm)       no. of glass beads        no. of stain-less steel balls     total no. of balls \n");
        while(!feof(fres1))
        {
            fscanf(fres1,"%f %f %f",&a[i],&b[i],&c[i]);           
            g = ((a[i]/5) + 1)*5;
            h = ((b[i]/5) + 1)*5;
            if(c[i] == 100)
            {
                ng[g][h] = ng[g][h]+1;  
            }
            else
            {
                ns[g][h] = ns[g][h] + 1;
            }
            ntot[g][h] = ntot[g][h] + 1;
            ++i;
        }  
        g = 5;
        h = 5;
        for(g=5;g<=350;g=g+5)
        { 
            for(h=5;h<=250; h=h+5)
            {
                fprintf(fres2,"%d       %d            %d                        %d                            %d\n",g,h,ng[g][h],ns[g][h],ntot[g][h]);  

            }
        }
    }

    fclose(fres1);
    fclose(fres2);
    return 0;    

}

我的输入 .dat 文件（即“trial_final_1.dat”）包含以下数据：

150.505951 0.797619 100.000000
172.327438 5.976130 100.000000
137.538651 11.089217 100.000000
151.304276 10.139803 100.000000
137.008926 13.175000 100.000000
120.400734 13.943015 1.000000
136.759262 14.199074 100.000000

在运行我的代码时，最初我得到输出“文件已成功打开，请检查您的文件夹中的 result_boxes 文件:)”，但在此之后我收到一条消息，指出 file.exe 已停止工作。请帮我解决这个问题。我也尝试过使用while(fscanf(fres1,"%f %f %f",&a[i],&b[i],&c[i])!=EOF )而不是while(!feof(fres1))，但两种情况下的输出保持不变。

4

1 回答 1

1

您没有对数组进行任何边界检查。快速浏览您的代码表明您超出了 的范围ng，ns并且ntot使用此循环：

for(g=5;g<=350;g=g+5)
{ 
    for(h=5;h<=250; h=h+5)
    {
        ...
    }
}

上面的问题是您使用<=而不是<，或者您的数组尺寸是一个元素太小。

fres2在写入之前，您也不会检查是否是有效的文件句柄，也不会检查您的调用是否fscanf成功。例如：

fscanf(fres1,"%f %f %f",&a[i],&b[i],&c[i]);           
g = ((a[i]/5) + 1)*5;
h = ((b[i]/5) + 1)*5;

如果上述调用失败，您可能正在使用a[i]or的未初始化值b[i]，然后您继续将其用作数组索引而无需任何边界检查。这是灾难的根源，而且它可能总是会发生，因为您的循环会在您尝试阅读之前检查 EOF。

至少你应该这样做：

if( 3 != fscanf(fres1,"%f %f %f",&a[i],&b[i],&c[i]) ) break;

还有一件事......你认为当i达到 6001 时会发生什么？

于 2013-06-11T05:30:25.413 回答