储存数据后再读取，为何读取的数据值与原始数据不一样？

13833943155

b = f_ssta(:);
fid=fopen('E:\WMN\M\obs_ssta_de_trending\data\1993-dec2021.det.00.dat','a');  %自动生成.dat文件
fwrite(fid, b, 'float32');   % 将数据写入文件
fclose(fid);


以上是储存代码；

fnm = 'E:\WMN\M\obs_ssta_de_trending\data\1993-dec2021.det.00.dat'
fid=fopen(fnm,'r','b');
data=fread(fid,'*float32');
fclose(fid);
data1 = zeros(134, 61, 348);
idx = 1;  % 存储数据的索引
for t = 1:348
    for j = 1:61
        for i = 1:134
            % 读取数据
            data1(i, j, t) =data(idx);
            % 更新索引
            idx = idx + 1;
        end
    end
end

wjy_ecnu

感觉应该是一样的  唯一看着像出问题的地方是fopen 建议以r的形式打开 a的形式如果源文件存在 会以追加的形式在后面写入 那自然两个文件不一样了

另外后面三层for循环效率太低了 可以用reshape提高效率

13833943155

wjy_ecnu 发表于 2023-8-15 11:45
感觉应该是一样的  唯一看着像出问题的地方是fopen 建议以r的形式打开 a的形式如果源文件存在 会以追加的形 ...

a这个地方，后来修改了。但是原始数据和储存再读取数据，只有数值为0的都还是0，其他的数值没有一样的

13833943155

原始数据数值和存储再读取的数值对比：