在 Matlab 中将未格式化的数据读入 4D 数组（与 IDL）

Question

我有一个未格式化（我认为是二进制）数据的 .img 文件，我想将其作为浮点读入 Matlab 中的 4D 数组中。我注意到Matlab中用于读取二进制数据的“fread”函数只能自动读入一维或二维数组，所以我尝试了这样的操作：

function fileArr4D= load4DFile(dim1, dim2, dim3, dim4, filename)

fileArr4D= zeros(dim1, dim2, dim3, dim4); %initialize array
fid= fopen(strcat('files/', filename));
for l= 1:dim4
    for k= 1:dim3
        fileArr4D(:, :, k, dim4)= fread(fid, [dim1, dim2], 'float');
    end
end
fclose(fid);
fileArr4D= flipdim(fileArr4D, 2); %flips along dimension 2 (of 4)

end

最后一行是翻转数据（特定于数据的东西），并且我将其包含在内的唯一原因是以防万一它是错误的根源。

不管怎样，这个函数给了我一个 4D 数组，但不是我想要的值。我实际上正在将一些代码从IDL翻译到Matlab，在IDL中，实现是：

pp=fltarr(dim1,dim2,dim3,dim4)
file='pathto/filename.img'
openr,1,file
readu,1,pp
pp=reverse(pp,2)
close,1

但是，将Matlab中fileArr4D返回的数组与IDL中的pp进行比较，数组的最小值、最大值和平均值都不同。关于这两种情况的不同做法，或者在 Matlab 中将二进制数据读入 4D 数组的更自然的方式，您有什么想法吗？

顺便说一句，如果我切换 Matlab 函数中的循环顺序（for k= 1:dim3, for l= 1:dim4, ... 而不是 for l= 1:dim4 ，对于 k= 1:dim3, ...），我得到了正确的最小值和最大值，但平均值错误（与之前的平均值相同）。

任何将 IDL 片段的行为引入 Matlab 的帮助都将不胜感激。

Answer 1

由于您想使用 fread 来访问数据，并且它只能读取一维或二维数组，因此我建议将数据作为一维数组加载。一旦您获得一维数组的数据，只要您知道维度应该是什么，您就可以对一维数组调用 reshape 函数，将其重新格式化为适当的 4D 数组尺寸。