fortran大型数组的二进制存储-matlab读取

fortran大型数组的二进制存储-matlab读取-可靠方法

转载▼

标签： 文化

fortran语言是一门高效的数值计算语言，用其来计算大型数组时非常方便的，而一些科学工作的后期处理往往要涉及到可视化的问题，而这个问题并非是fortran的专长，我个人比较偏好matlab，因此写小文来阐述一个从fortran中生成的数组或矩阵如何移植到matlab中，以便进行后期数据处理分析。
过程大概可以分化为(1)生成数组--->(2)存为二进制格式的文件--->(3)matlab的读取
为了使得这个过程简洁明了，假定我们要生成的数组为a(100,121,21),即为100*121*21的矩阵或数组，里面的数据为float型(kind=4)，即这个数组是一个较大的数组,为了生成和检验我们的方法是否有效，我们假定a(31,28,20)=1.，而其他的数据均为0.，这样我们可以最后仅仅检验点a(31,28,20)上的数据是否为1即可。

第一步：生成数组
program main
implicit none
real(kind=4)::a(100,121,21)
integer::k
a=0.
a(31,28,20)=1. ！！！！注释：这里将点a(31,28,20)=1.，其他点都为零
第二步：存为二进制格式的文件
open(12,file='test.bin',access='direct',form='unformatted',recl=100*121)
！！！！ 这里我们把100*121即数组的前二维当做一个记录，而每个数据为4字节，故recl=100*121*4/4，
！！！！若为双精度recl=100*121*8/4
do k=1,21,1
write(12,rec=k) a(:,:,k) ！！！！存储
enddo
end program main

第三步：matlab的读取
！！！！由于matlab中的数组是按列存的，而fortran是按行存的，故不能直接进行读取，需要用到一个matlab程序进行转换，即我给出的Newshape程序，这里L=100,M=121,N=21，调用方式为Aout=Newshape(100,121,21)
function[Aout]=Newshape(L,M,N)
fid=fopen('test.bin','r')
[A,count]=fread(fid,inf,'float')
for k=1:N
for j=1:M
for i=1:L
Aout(i,j,k)=A(M*L*(k-1)+L*(j-1)+i);
end
end
end

结束语：matlab中输入Aout(31,28,20),就可以得出1，说明了我们的方法真实可靠。