scipy.io.

FortranFile#

类 scipy.io.FortranFile(filename, mode='r', header_dtype=<class 'numpy.uint32'>)[源代码]#

用于处理Fortran代码生成的非格式化顺序文件的文件对象。

参数:
filename文件对象或字符串

已打开的文件对象或文件名。

mode{‘r’, ‘w’}，可选

读写模式，默认为'r'。

header_dtype数据类型，可选

头部的数据类型。大小和字节序必须与输入/输出文件匹配。

方法

close()

关闭文件。

read_ints([dtype])

从文件中读取指定类型的一条记录，默认为整数类型（Fortran中的INTEGER*4）。

read_reals([dtype])

从文件中读取指定类型的一条记录，默认为浮点数（Fortran中的real*8）。

read_record(*dtypes, **kwargs)

从文件中读取指定类型的一条记录。

write_record(*items)

向文件写入一条记录（包括大小）。

备注

这些文件被分解为未指定类型的记录。每条记录的大小在开头给出（尽管此头部的大小并非标准），数据以无格式写入磁盘。支持BACKSPACE语句的Fortran编译器将写入大小的第二个副本，以方便向后查找。

此类别仅支持记录大小被写入两次的文件。它也不支持Intel和gfortran编译器中用于处理大于2GB且带有4字节头部的记录的子记录。

一个Fortran非格式化顺序文件的例子会这样写入：

OPEN(1, FILE=myfilename, FORM='unformatted')

WRITE(1) myvariable

由于这是一种非标准的文件格式，其内容取决于编译器和机器的字节序，因此建议谨慎使用。已知gfortran 4.8.0 和 gfortran 4.1.2 在 x86_64 上的文件可以正常工作。

建议使用Fortran直接访问文件或新式流I/O文件，这些文件可以很容易地通过numpy.fromfile读取。

示例

要创建一个非格式化顺序Fortran文件：

>>> from scipy.io import FortranFile
>>> import numpy as np
>>> f = FortranFile('test.unf', 'w')
>>> f.write_record(np.array([1,2,3,4,5], dtype=np.int32))
>>> f.write_record(np.linspace(0,1,20).reshape((5,4)).T)
>>> f.close()

要读取此文件：

>>> f = FortranFile('test.unf', 'r')
>>> print(f.read_ints(np.int32))
[1 2 3 4 5]
>>> print(f.read_reals(float).reshape((5,4), order="F"))
[[0.         0.05263158 0.10526316 0.15789474]
 [0.21052632 0.26315789 0.31578947 0.36842105]
 [0.42105263 0.47368421 0.52631579 0.57894737]
 [0.63157895 0.68421053 0.73684211 0.78947368]
 [0.84210526 0.89473684 0.94736842 1.        ]]
>>> f.close()

或者，在Fortran中：

integer :: a(5), i
double precision :: b(5,4)
open(1, file='test.unf', form='unformatted')
read(1) a
read(1) b
close(1)
write(*,*) a
do i = 1, 5
    write(*,*) b(i,:)
end do
在本页