fortran大型数组的二进制存储-matlab读取
来源:互联网 发布:c标签和js函数 编辑:程序博客网 时间:2024/05/17 23:38
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,说明了我们的方法真实可靠。
过程大概可以分化为(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即可。
第一步:生成数组
第二步:存为二进制格式的文件
!!!! 这里我们把100*121即数组的前二维当做一个记录,而每个数据为4字节,故recl=100*121*4/4,
!!!!若为双精度recl=100*121*8/4
第三步: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
end
结束语:matlab中输入Aout(31,28,20),就可以得出1,说明了我们的方法真实可靠。
- fortran大型数组的二进制存储-matlab读取
- (转)Fortran、matlab 和 C语言中数组的存储方式
- 存储、读取二进制的图像
- 字符串的二进制存储读取
- Fortran和 C语言中数组的存储方式
- matlab读取二进制数据文件的方法
- FORTRAN的数组-笔记
- Fortran 数组的声明
- 二进制流的存储和读取
- 二、Matlab图像的读取和存储
- Matlab图像的读取和存储
- Python读取MATLAB存储的mat文件
- python读取存储matlab的 .mat文件
- FORTRAN 数组的下标问题
- Fortran 判断数组的 维度。
- matlab与fortran的数据传递
- NET 二进制形式的文件的存储与读取
- C++与MATLAB数组的存储结构
- 最简化的DirectX 11开发环境的配置
- 虚拟摄像头驱动原理及开发
- struts2 中同一个action的实现中对应多个input的处理方法
- 进步了!
- 在GitHub上分享和展示你的代码
- fortran大型数组的二进制存储-matlab读取
- hibernate3中出现 could not initialize proxy - no Session 错误的解决办法以及lazy load的介绍
- EasyUI Tree+Asp.net实现权限树或目录树导航
- 关于platform_device和platform_driver的匹配
- 联合查询与多表查询基础详解
- hdu1355 The Peanuts
- 资源内存管理
- hibernate前2天的简单总结
- struts2 iterator标签使用