以逗号为分割的数据如何读取

Zhuming · 发表于 2020-5-13 15:54:56

数据：2004/08/10,20:24:38.300966, 30.1426, 104.7854, 13.9, 52004/08/10,20:24:39.380542, 10.4838, 29.3264, 18.7, ...
TXT格式，在txt中没有对齐。我需要从很多这样的数据中提取出一些符合我需求的数据。数据描述：日期，时间，纬度，经度，误差值。数据比较乱，有的也有空格。
我尝试先读取一条数据，失败了，之前做的都是对齐的数据。代码也很不好意思的附上吧....

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

01program wwlln_shiyan
02implicit none
03 
04!声明变量
05character(len=20)::filename
06character(len=20)::wwlln_date,wwlln_time
07real(kind=8)::lat,lon
08real::timingerror
09integer::i
10 
11!读取文件
12filename='A20040810.txt'
13OPEN(11,file=filename,status='old',form='formatted')
14OPEN(12,file='shiyan.txt',status='replace')
15READ(11,*) wwlln_date,wwlln_time,lat,lon,timingerror
16write(*,*)wwlln_date,wwlln_time,lat,lon,timingerror
17write(12,*)wwlln_date,wwlln_time,lat,lon,timingerror
18CLOSE(11)
19CLOSE(12)
20 
21STOP
22end

Zhuming · 发表于 2020-5-13 16:00:31

补充：我以上代码出来的结果是这样的： 2004 戢三 2.1138223167426894E-317 0.0000000000000000 0.00000000
Note: The following floating-point exceptions are signalling: IEEE_DENORMAL

necrohan · 发表于 2020-5-13 16:17:04

用格式读入，把
READ(11,*) wwlln_date,wwlln_time,lat,lon,timingerror
改成
READ(11,'(a27,f,f,f)') wwlln_time,lat,lon,timingerror ! wwlln_date,

Zhuming · 发表于 2020-5-13 17:33:02

necrohan 发表于 2020-5-13 16:17
用格式读入，把
READ(11,*) wwlln_date,wwlln_time,lat,lon,timingerror
改成

你好，感谢您的回复。我按照您的提示，修改为以下代码，输出结果已经很接近了，但还是有一些问题。我先附上代码及输出结果。

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

01program wwlln_shiyan
02implicit none
03 
04!声明变量
05character(len=20)::filename
06character(len=20)::wwlln_date,wwlln_time
07real::lat,lon
08real::timingerror
09integer::i
10 
11!读取文件
12filename='A20040810.txt'
13OPEN(11,file=filename,status='old',form='formatted')
14OPEN(12,file='shiyan.txt',status='replace')
15READ(11,"(A11,A16,F7.4,F9.4,F4.1)")wwlln_date,wwlln_time,lat,lon,timingerror
16write(*,*)wwlln_date,wwlln_time,lat,lon,timingerror
17write(12,*)wwlln_date,wwlln_time,lat,lon,timingerror
18CLOSE(11)
19CLOSE(12)
20 
21STOP
22end

结果： 2004/08/10 20:24:38.300966 30.1425991 104.785400 13.0000000
对比原文：2004/08/10 20:24:38.300966 30.1426 104.7854 13.9
可以看出，还是有一些误差，是我读取格式设置不对的原因导致的。我本人正在翻书中，也希望各位大佬不吝赐教！谢谢！

Zhuming · 发表于 2020-5-13 17:38:05

另外数据原文是这样的：2004/08/10,20:24:38.300966, 30.1426, 104.7854, 13.9, 52004/08/10,20:24:39.380542, 10.4838, 29.3264, 18.7, ......
我在txt中去除了所有的空格，然后将“，”替换为空格，操作结束后的数据为下：
2004/08/10 20:24:38.300966 30.1426 104.7854 13.9 52004/08/10 20:24:39.380542 10.4838 29.3264 18.7 ......
而实际上，这样的文本我一共要处理5400个，甚至更多，肯定不可能一个接一个的处理的，我可能还需要尝试其他方法。

fcode · 发表于 2020-5-14 08:09:52

这样做，就不需要替换逗号为空格了。

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

01program wwlln_shiyan
02  implicit none
03 
04  !声明变量
05  character(len=512)::cRead
06  character(len=20)::filename
07  character(len=20)::wwlln_date,wwlln_time
08  real::lat,lon
09  real::timingerror
10  integer::i
11 
12  !读取文件
13  filename='A20040810.txt'
14  OPEN(11,file=filename,status='old',form='formatted')
15  OPEN(12,file='shiyan.txt',status='replace')
16  READ(11,"(A512)") cRead
17  call parserRead(cRead)
18  write(*,*)wwlln_date,wwlln_time,lat,lon,timingerror
19  write(12,*)wwlln_date,wwlln_time,lat,lon,timingerror
20  CLOSE(11)
21  CLOSE(12)
22   
23contains
24 
25  Subroutine parserRead(c)
26    character(len=*) :: c
27    integer :: i
28    do i = 1 , len_trim(c)
29      if( c(i:i) == '/' ) c(i:i) = '|'
30    end do
31    read(c,*) wwlln_date,wwlln_time,lat,lon,timingerror
32    do i = 1 , len_trim(wwlln_date)
33      if( wwlln_date(i:i) == '|' ) wwlln_date(i:i) = '/'
34    end do
35  End Subroutine parserRead
36   
37end program wwlln_shiyan

Zhuming · 发表于 2020-5-14 09:51:10

fcode 发表于 2020-5-14 08:09
这样做，就不需要替换逗号为空格了。
[mw_shl_code=fortran,true]program wwlln_shiyan
implicit none

学到了！！！感谢！

necrohan · 发表于 2020-5-14 11:06:40

Zhuming 发表于 2020-5-13 17:33
你好，感谢您的回复。我按照您的提示，修改为以下代码，输出结果已经很接近了，但还是有一些问题。我先附 ...

浮点数的格式读入不需要指定宽度
READ(11,"(A11,A16,F7.4,F9.4,F4.1)")wwlln_date,wwlln_time,lat,lon,timingerror
改为
READ(11,"(A11,A16,F,F,F)")wwlln_date,wwlln_time,lat,lon,timingerror
可以直接读入数据
原文件格式不需要处理空格和逗号，可以直接用
我用Inter Fortran9.0读入没问题

Zhuming · 发表于 2020-5-14 19:22:16

necrohan 发表于 2020-5-14 11:06
浮点数的格式读入不需要指定宽度
READ(11,"(A11,A16,F7.4,F9.4,F4.1)")wwlln_date,wwlln_time,lat,lon,ti ...

那应该是编译器的不同，我用的simply Fortran，这么写会报错，所以才指定宽度去写的。
不过 @fcode 给出的解决方法真的很好用。整个项目我已经做完了！
再次谢谢你！

Rookie · 发表于 2022-4-14 14:51:10

fcode 发表于 2020-5-14 08:09
这样做，就不需要替换逗号为空格了。
[mw_shl_code=fortran,true]program wwlln_shiyan
implicit none

请问这个修改后的代码，是先将数据读取到字符串“cRead”中，然后通过自定义的子程序“parserRead”将字符串“cRead”视作一个内部文件进行读取吗？

[文件读写] 以逗号为分割的数据如何读取

评分

规矩勋章

美女勋章

热心勋章

星光勋章

新人勋章

贡献勋章

管理勋章

帅哥勋章

爱心勋章

元老勋章

水王勋章