Question

我有一个带有3600字节标头的二进制文件，13483条跟踪（列），每条标头有一个240字节的标头。我想跳过标题并将数据值读入矩阵。

我能够从文件中获取一些值，但是地震数据[50] [40]到seismicDataNH [50] [54]应该是[13,17,12,5,19,51,29， - 118，-127，-127，-50,126,126,126，-32]这不是我得到的。

我不确定我是否正确理解fread（），它是将文件读取为一行长行还是多行？我假设有一个长行，也许这就是为什么它不起作用。

这是我为阅读文件而写的代码：

#include <iostream>

using namespace std;

#define N_SAMP 1990
#define M_TR 13483

char tempArray [N_SAMP*M_TR];
char seismicData[1990][13483];
char seismicDataNH[1750][13483];

int main()
{
    FILE*seismicFile;
    seismicFile = fopen("NVGT-88-06.sgy","rb");

    if (seismicFile!=NULL)
    {
        fseek(seismicFile, 3600*sizeof(char), SEEK_CUR);
        fread(tempArray, sizeof(char), N_SAMP*M_TR, seismicFile);
        puts("\n\nRead File successfully");

        int c = 0;
        for (int in=0; in<N_SAMP; in++)
        {
            for (int im=0; im<M_TR; im++)
            {
                seismicData[in][im] = tempArray[c];
                c++;
            }
        }

        puts("\nStored in matrix");


        // Make matrix values without header values
        for (int in=240; in < N_SAMP; in++)
        {
            for(int im=0; im < M_TR; im++)
            {
                seismicDataNH[in-240][im] = seismicData[in][im];
            }
        }
        puts("Removed header");



        puts("Test values: \n");
        for (int it = 40; it<55; it++)
        {
            printf("%d\n", seismicDataNH[50][it]);
        }

        fclose(seismicFile);

    }
    return 0;
}

这里是数据文件（.sgy），如果有人想看看它： https://www.dropbox.com/s/y8aa99yqhfyacc8/NVGT-88-06.sgy?dl=0

Answer 1

根据您的描述，共有13483个连续的1990字节块（包括标题为240）。

这意味着你有一个for循环嵌套和数组索引错误的方法。

将数组定义更改为：

char seismicData[M_TR][N_SAMP];
char seismicDataNH[M_TR][N_SAMP-240];

两个嵌套for循环：

for (int im=0; im<M_TR; im++)
{
    for (int in=0; in<N_SAMP; in++)
    {
        seismicData[im][in] = tempArray[c];
        c++;
    }
}

RESP。：

for(int im=0; im < M_TR; im++)
{
    for (int in=240; in < N_SAMP; in++)
    {
        seismicDataNH[im][in-240] = seismicData[im][in];
    }
}

保持最终的for循环（打印数据）。

这应该给你预期的输出（它对我有用）。

如何在C ++中读取未知形状的二进制文件？

1 个答案: