Fortran Coder

标题: fortran读取有缺失的数据问题 [打印本页]

作者: daydaycode 时间: 2020-3-24 15:52
标题: fortran读取有缺失的数据问题
本帖最后由 daydaycode 于 2020-3-24 17:22 编辑

data.txt (12.64 KB, 下载次数: 3)

fortran小白请教下fortran处理txt文件问题。需要处理的数据为地面自动站气象数据，数据包含的数据类型主要有字符串和浮点及整数等类型格式，数据为62247*212的矩阵，但是在57418行到最后的前6列都没有值，这样在数据读取的时候就会提示读取文件末尾报错，请问想正确读入这类型的数据，该怎么编写fortran代码。自己参考别人的代码，写的fortran程序，不能处理以上的问题，代码及数据和报错信息如下，还请各位前辈指教下怎么修改

[Fortran] 纯文本查看 复制代码

PROGRAM predeal_dimianzidong

IMPLICIT NONE



INTEGER                                ::i,j                !读取数据的行数

INTEGER, PARAMETER      ::row=62247,col=212

CHARACTER,DIMENSION(row,col)     ::data_read,data_1



OPEN(UNIT=3, FILE='SURF_CHN_MUL_HOR_2019080100.txt', STATUS='OLD',ACTION='READ',FORM='FORMATTED',ACCESS='sequential')        !打开文件

!READ(3,*) ((data_read(i,j),j=1,col),i=1,row)

DO i = 1,row

READ(3,*)(data_read(i,j),j=1,col)

END DO

!data_1 = transpose(data_read) 

print*,SHAPE(data_read)

!PRINT*,data_1

CLOSE(UNIT=3)        !关闭文件

STOP

END PROGRAM predeal_dimianzidong

截图.png (118.2 KB, 下载次数: 541)

数据结构类型

2.png (25.6 KB, 下载次数: 524)

报错情况

4.png (59.44 KB, 下载次数: 553)

读取每行数据

3.png (36.77 KB, 下载次数: 581)

读取每行数据

微信截图_20200324172002.png (87.63 KB, 下载次数: 527)

缺失处截图

作者: li913 时间: 2020-3-24 16:07
本帖最后由 kyra 于 2020-3-24 17:19 编辑

把你的数据上传，只需要57418行前后几行的数据就行。

[Fortran] 纯文本查看 复制代码

DO i = 1,row

if(i<57418) then

READ(3,*) data_read(i,:)

else

READ(3,*) data_read(i,7:)

end if

END DO

作者: daydaycode 时间: 2020-3-24 16:19

li913 发表于 2020-3-24 16:07
把你的数据上传，只需要57418行前后几行的数据就行。

DO i = 1,row

您好，数据已上传，麻烦给看一下

作者: daydaycode 时间: 2020-3-24 16:38

li913 发表于 2020-3-24 16:07
把你的数据上传，只需要57418行前后几行的数据就行。

DO i = 1,row

按照您的代码修改，还是不行，我用论坛读取非固定格式每行的数据数那个程序，发现我的数据第一行只有184个数据，还有的行没有211个数据，只有210个，像这种情况该怎么处理呢。

作者: vvt 时间: 2020-3-24 17:15
你首先需要找到一个规则。再什么情况下，你认为第一列是缺失的，什么情况下你认为是第二列缺失的？

比如：

姓名年龄体重
张三 30 55.3
李四 40

你如何判断第二行的40，是年龄还是体重？

作者: daydaycode 时间: 2020-3-24 17:21

vvt 发表于 2020-3-24 17:15
你首先需要找到一个规则。再什么情况下，你认为第一列是缺失的，什么情况下你认为是第二列缺失的？

比如： ...

这个要怎么来定呢，我也不能确定它是哪一列会有缺失，我上传了有缺失的地方截图，你可以看下。

作者: vvt 时间: 2020-3-24 17:24
我既然回复你的帖子，肯定都认真看过你提供的所有内容了。

我没有办法判断，因为你不管说 40 是年龄，还是体重，都是不违背任何规则的。

你是问题的提出者，应该你来确定问题的规则。然后根据问题的规则来书写合适的代码。

作者: vvt 时间: 2020-3-24 17:25
你可以看看你的数据是如何得来的，是否可以转存为其他数据格式？比如 csv，用逗号间隔开。

比如
姓名  , 年龄 , 体重
张三  , 30 , 55.3
李四   ,,   40

因为第3行有两个逗号，我可以认定，40应该是体重，而年龄缺失。

作者: li913 时间: 2020-3-25 10:56
看起来你的文件中数据是以tab分隔的，你可以用excel打开数据，检查一下是否有误，然后另存为csv格式。读取csv格式就方便了。

[Fortran] 纯文本查看 复制代码

program test

  implicit none

  integer i

  INTEGER, PARAMETER      ::row=11,col=212

  character(40) dat(col-6, row)

  dat = 'undefined'

  open(10,file='1.csv')

  read(10,*)

  do i = 1, row

    read(10,*) dat(:,i)

  end do

  close(10)

end program

作者: daydaycode 时间: 2020-3-25 13:47

li913 发表于 2020-3-25 10:56
看起来你的文件中数据是以tab分隔的，你可以用excel打开数据，检查一下是否有误，然后另存为csv格式。读取c ...

十分感谢，现在问题解决了。请问，像这种情况只能把txt转成csv吗，还有别的方法吗？

作者: daydaycode 时间: 2020-3-25 14:24

vvt 发表于 2020-3-24 17:25
你可以看看你的数据是如何得来的，是否可以转存为其他数据格式？比如 csv，用逗号间隔开。

比如

好的，多谢。

作者: li913 时间: 2020-3-25 17:50

daydaycode 发表于 2020-3-25 13:47
十分感谢，现在问题解决了。请问，像这种情况只能把txt转成csv吗，还有别的方法吗？ ...

有的，更麻烦。读取整行，以tab为分隔符，将一行分为多个记录，存入数组。

作者: daydaycode 时间: 2020-3-26 10:42

li913 发表于 2020-3-25 17:50
有的，更麻烦。读取整行，以tab为分隔符，将一行分为多个记录，存入数组。 ...

好的，多谢。

作者: Rookie 时间: 2022-4-20 14:54

vvt 发表于 2020-3-24 17:25
你可以看看你的数据是如何得来的，是否可以转存为其他数据格式？比如 csv，用逗号间隔开。

比如

请问这种缺失的数据，读取的时候会保存为什么？

作者: fcode 时间: 2022-4-20 17:14

Rookie 发表于 2022-4-20 14:54
请问这种缺失的数据，读取的时候会保存为什么？

没有统一规则。看具体“保存该文件的程序”的策略。

作者: Rookie 时间: 2022-4-22 08:11

fcode 发表于 2022-4-20 17:14
没有统一规则。看具体“保存该文件的程序”的策略。

好的。关于字符串的缺失，我的电脑直接输出空白。所以我在派生数据的声明阶段，直接初始化了。
PS.我所编译的程序，在读取外部文件的TXT文件时候，到达文档结尾的时候，IOSTAT仍为0.不知道前辈们有没有遇到过类似的情况？

作者: Rookie 时间: 2022-4-22 08:18

Rookie 发表于 2022-4-22 08:11
好的。关于字符串的缺失，我的电脑直接输出空白。所以我在派生数据的声明阶段，直接初始化了。
PS.我所编 ...

刚刚发现是txt文档数据结尾有两个空白行，才会导致iostat仍为0.

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)