问题描述:
segy地震数据读取过程中采用如下语句读取地震道:
FILE *fp = fopen(segyfilename, "r");
fseek(fp, 3600,SEEK_SET);
int ns = 999;
size_t segysize = 240 + sizeof(float)* ns;
size_t act_size_of_read = fread((char*)&segy, 1, segysize, fp);
assert( act_size_of_read == segysize);
......
程序运行失败!
调试:
发现 ns = 999 , segysize=4236 , act_size_of_read = 1683 < segysize
反复检查原始地震数据道,没有异常。
再次检查代码,发现文件打开方式有问题,更改"r"为"rb"后,运行成功!
小结:
细节决定成败。不用盲目相信自己的代码!