1

到目前为止,这是我的代码:

#include <stdio.h>
#include <stdint.h>

#define N_DET 432
#define N_PROJ 500

int readSinogram(uint16_t mess[],char filename[], int sample){
    //init
    FILE *fp;
    int i;
    //open file
    fp=fopen(filename,"r");
    //read sinogram
    fseek(fp,sizeof(uint16_t)*N_DET*N_PROJ*sample*2,SEEK_SET); //move pointer
    fread(mess,sizeof(uint16_t),N_DET*N_PROJ,fp); //read sinogram
    return 0;
};


int main()
{
    uint16_t mess[N_DET*N_PROJ];
    char filename[]="C:\\Path\\MyFile.fxc";
    double curr;
    int i,DET,PROJ;

    readSinogram(mess,filename,0); //read the data

    printf("\nDET?"); //ask for input
    scanf("%u", &DET);
    printf("\nPROJ?");
    scanf("%u", &PROJ);

    curr=mess[DET+N_DET*PROJ]; //get the data
    printf("Data: %f",curr); //print the data

    return 0;
}

这会读入 .fxc 文件,该文件只是一个包含 uint16_t 格式数字的二进制文件。"readSinogram" 读取一个包含 N_DET*N_PROJ 个数的数据集。(指针移动了两倍的块大小,因为文件中有两个测量值的交替块。)读完后,可以放入 DET 和 PROJ 并查看此时的数据。

到目前为止,一切正常。数据对于某个数据范围是正确的,但是当要求太大的 DET 和/或 PROJ 时,数据是不正确的。我确实在 Matlab 中读入了相同的文件,并且可以确认其中的数据没有问题。

更准确地说:[DET+N_DET*PROJ] > 248835 以上的每个索引都将返回 52428 而不是正确的值(范围从 0 到 4088)。达到这些值的值可以正常工作。所以我猜该索引上方的“混乱”数组一定有问题。

我究竟做错了什么?提前致谢!

4

1 回答 1

1

您需要分配一个更大的数组来获得大于 248835 的索引值目前您有以下定义

#define N_DET 432
#define N_PROJ 500

这导致数组大小为432*500 = 216000

现在,如果您在值248835 248835 > 216000处索引数组,您将导致从为数组分配的内存中访问内存,这将导致未定义的行为。您需要的是一个更大的数组,可以容纳超过 248835 个条目。一个简单的#define N_DET 500就可以做到这一点,但你必须确保这是要求。基本的索引检查将有助于避免出现超出范围的问题

#define MAX_RANGE 500*500
if((DET+N_DET*PROJ) < MAX_RANGE)
    curr=mess[DET+N_DET*PROJ]; //get the data
else
   //error handling

您还需要错误处理fopen fseek并检查返回值,fread您可以选择关闭打开的文件描述符close(fd)

于 2015-03-11T23:35:59.523 回答