在R中读取湖泊深度的二进制文件

Question

我正在尝试在R中打开一个文件，它是二进制文件并用Fortran编写。该文件名为GlobalLakeDepth.dat，位于：http://www.flake.igb-berlin.de/gldbv2.tar.gz

说明指定要打开GlobalLakeDepth.dat（在Fortran中），需要执行以下操作：

An example of opening the binary file in FORTRAN90:
-- open(1, file = 'GlobalLakeDepth.dat', form='unformatted', access='direct', recl=2)

An example of reading the binary file in FORTRAN90:
-- read(1,rec=n) LakeDepth
-- where: n - record number, INTEGER(8);
          LakeDepth - mean lake depth in decimeters, INTEGER(2).

我的问题是：在Fortran中给出这些说明，如何在R中打开此文件？也就是说，这样做有“R方式”吗？

我一直在关注http://www.ats.ucla.edu/stat/r/faq/read_binary.htm的说明，但是，我仍然没有从数据文件中获取任何内容。我只需要提供有关36个大型湖泊的测量湖泊测量数据的信息。

Answer 1

您可以使用readBin来读取二进制文件。对于这个文件，我认为正确的命令是

lk <- readBin("GlobalLakeDepth.dat", n = 43200 * 21600, what = "integer", endian = "little", size = 2)

这使得非常长向量可以制作成43200 * 21600矩阵。