How to read an N-dimensioned matrix from a binary file
16 次查看(过去 30 天)
显示 更早的评论
I would like to write a Matlab script to read an N-dimensioned matrix created in Fortran. Specifically, the Fortran code wrote to a binary file using the following:
write(30) a, b, c
with the following declarations for a,b,c
real a(16,20,22,6,3, 2)
real b(16,20,22,6,3,13)
real c(16,20,22,6,3,13)
In Matlab, it seems that fread handles only a vector or 2-dimesional matrix fread(fileID, sizeA) where sizeA is a vector of size n or 2-dimensional matrix [m,n].
How can I read the N-dimensioned matrix I described?
0 个评论
采纳的回答
James Tursa
2013-2-13
编辑:James Tursa
2013-2-13
Read the variable as a 1D using the total number of elements (e.g., 16*20*22*6*3*2) and then reshape the result. Since Fortran stores nD elements in memory in the same order that MATLAB does (column major) this will work. HOWEVER, unless you opened the file in Fortran as a stream, you will likely get record headers on the Fortran side that you will need to skip on the MATLAB side. Probably 4 bytes, but could be something else depending on the source machine and compiler used. E.g., something like this might work:
fid = fopen('yourfilename');
h = fread(fid,1,'*uint32'); % May need adjusting
a = fread(fid, 16*20*22*6*3*2, '*single');
b = fread(fid, 16*20*22*6*3*13, '*single');
c = fread(fid, 16*20*22*6*3*13, '*single');
a = reshape(a,[16 20 22 6 3 2]);
b = reshape(b,[16 20 22 6 3 13]);
c = reshape(c,[16 20 22 6 3 13]);
3 个评论
James Tursa
2013-2-13
编辑:James Tursa
2013-2-13
Your open statement is not stream, it is simply a sequential unformatted file. Thus I would expect to see the record header as I indicated. Since this header is machine and compiler specific (not part of the Fortran Standard) there is no guaranteed way to know for an arbitrary file what (if anything) is there. The only thing to do in a general sense is trial-and-error until you discover what works. But installations I have used have 4 bytes up front (e.g., Intel Fortran on a PC) so I would start with that as a first guess. Give my code a shot and see if it works.
更多回答(1 个)
Aniruddh Murali
2018-9-27
If I am opening a bin in Fortran using the following command
open(unit=1,file='24092018.bin',
1 form='unformatted',status="old",access='sequential')
is there a way to open the same file in Matlab?
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Fortran with MATLAB 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!