关于二进制文件读取
本帖最后由 wxy 于 2017-9-1 11:38 编辑图一是二进制文件的基本属性，这是一个491*411的格点数据，我在读取时只读到1*411就end of file 了，我让recl=mrow*mcol,结果只输出到1*1，我知道是RECL这里设置的有问题，但不清楚这种二维的数据该怎么设置recl,
program main
implicit none
integer*4,parameter ::Mrow=411, Mcol=491,mxp=201801
real*4 ::grid1(Mcol,Mrow)
integer*4 irow,icol,status, kk0
open(7,file='xd10kmlr1.bil',form="unformatted",access='direct',recl=mcol*4,iostat=status)
open(8,file="bil.txt")
kk0 = 0
do icol=1,mcol
do irow=1,mrow
kk0=kk0+1
read(7,rec=kk0)grid1(icol,irow)
print*,icol ,irow,grid1(icol,irow)
!read(20, rec=kk0)(grid2(icol,irow,i),i=1,MXP)
write(8,"(i3,2x,i3,2x,f7.2)") icol ,irow,grid1(icol,irow)
end do
end do
close(7)
end 可以考虑用 stream 方式读写，而不用 direct fcode 发表于 2017-9-1 11:40
可以考虑用 stream 方式读写，而不用 direct
嗯嗯 成功了 谢谢 本帖最后由 wxy 于 2017-9-2 12:23 编辑
fcode 发表于 2017-9-1 11:40
可以考虑用 stream 方式读写，而不用 direct
还想请教一下，我在读取32位数据时可以491*411全部读出来，但换成16位格式数据时
只读取出一半，而且读取出的值是很长一串数值，并不是想要的结果，我怎样可以设定成16位呢，我用的gfortranprogram main
implicit none
integer,parameter ::Mrow=411, Mcol=491
integer ::grid(Mcol,Mrow)
integer ::irow,icol
open(7,file='xdem10km1.bil',access='stream',form="unformatted")
open(8,file="bil.txt")
do icol=1,mcol
do irow=1,mrow
read(7)grid(icol,irow)
print*,icol ,irow,grid(icol,irow)
!read(20, rec=kk0)(grid2(icol,irow,i),i=1,MXP)
write(8,*) icol ,irow,grid(icol,irow)
end do
end do
close(7)
end
wxy 发表于 2017-9-2 12:21
还想请教一下，我在读取32位数据时可以491*411全部读出来，但换成16位格式数据时
只读取出一半，而且读取 ...
1. 16位格式数据, 2-Byte INTEGER? 請試試
INTEGER, PARAMETER :: ibyte2 = SELECTED_INT_KIND(r= 4) !2-Byte INTEGER
INTEGER(KIND=ibyte2) :: grid(mcol,mrow)
2. DIRECT/STREAM access, 不保證完全相容, 也可能compiler dependent
3.图一是二进制文件的基本属性，这是一个491*411的格点数据，我在读取时只读到1*411就end of file
open(7,file='xd10kmlr1.bil',form="unformatted",access='direct',recl=mcol*4,iostat=status)
-> open(7,file='xd10kmlr1.bil',form="unformatted",access='direct',REAL=4,iostat=status) ! 4-byte
-> open(7,file='xd10kmlr1.bil',form="unformatted",access='direct',RECL=1,iostat=status) ! one "4-byte" data
RECL=?, compiler dependent
chiangtp 发表于 2017-9-2 15:44
1. 16位格式数据, 2-Byte INTEGER? 請試試
INTEGER, PARAMETER :: ibyte2 = SELECTED_INT_KIND(r= 4) ...
嗯嗯 明白了 谢谢 wxy 发表于 2017-9-2 12:21
还想请教一下，我在读取32位数据时可以491*411全部读出来，但换成16位格式数据时
只读取出一半，而且读取 ...
我用你的代码读到一半就end of file了 请问该怎么改呢
页:
[1]