Matlab和Fortran之间的数据解释

Question

我有一些解释二进制文件的Fortran代码。 由于另一个问题，我问我了解了fortran代码的工作原理，但是我需要将其更改为MATLAB代码的帮助。 这是Fortran代码：

  IMPLICIT NONE

  DOUBLE PRECISION      SCALE

  CHARACTER*72          BLINE
  CHARACTER*72          DATA_FILE_FULL_NAME
  CHARACTER*6           DATA_FILE_NAME

  CHARACTER*4           CH4, CH4F
  REAL*4                RL4
  EQUIVALENCE          (RL4,CH4)

  CHARACTER*2           C2
  LOGICAL               LFLAG

  INTEGER*2             I2

  if(i2.eq.13881) LFLAG=.FALSE.

  c2='69'
  LFLAG=.TRUE.

  DATA_FILE_FULL_NAME='./'//DATA_FILE_NAME//'.DAT'

  OPEN(UNIT=20, FILE=DATA_FILE_FULL_NAME, ACCESS='DIRECT',
 .     RECL=72, status='OLD')
    READ(20,REC=1) BLINE
    CH4f=BLINE(7:10)
    call flip(4,lflag,ch4f,ch4)
    SCALE=RL4

  RETURN
  END  

c   ..................................................
    subroutine flip(n,lflag,ch1,ch2)
c   ..................................................

  integer*4        n, i
  character*(*)    ch1, ch2
  logical          lflag

  if(lflag) then
    do i=1,n
      ch2(i:i)=ch1(n-i+1:n-i+1)
    enddo
  else
    ch2=ch1
  endif

  return
  end   

因此，基本上（我相信我理解正确），Fortran代码正在检查字节序，并在机器和数据不匹配时翻转字节序。 此外，数据存储在内存中，同时为CH4和RL4保留，因此在定义CH4之后通过调用RL4，可以将数据简单地转换为REAL * 4类型而不是CHARACTER * 4类型。 现在，我需要弄清楚如何将其移植到Matlab中。 我已经可以读取原始数据，但是当我尝试各种形式的解释时，总是会得到错误的答案。 到目前为止，我已经尝试过：

fid=fopen(LMK_file,'r');
BLINE=fread(fid, 72,'char=>char');
CH4=BLINE(7:10);
SCALE=single(CH4(1)+CH4(2)+CH4(3)+CH4(4));

SCALE=int8(CH4(1)+256*int8(CH4(2))+256^2*int8(CH4(3))+256^3*int8(CH4(4));

在MATLAB中，但两者都给了我相同的错误答案（这很有意义，因为它的作用几乎相同）。

有什么建议么？

Answer 1

对于正在寻找该问题答案的其他人，这是我想出的一切，这要归功于评论中的帮助。

fid=fopen(LMK_file,'r');  %open the file for reading
IGNORE=fread(fid,6,'*char');  %ignore the first 6 bytes.  This could also be performed using fseek
SCALE=fread(fid,1,'*float32','b');  %reading in the scale value as a floating point 32 bit number. 

其中b表示文件的大结尾。 通过使用'* float32'选项，fread自动将接下来的4个字符解释为32位浮点数。

这将从给定文件返回正确的比例值。