Fortran Coder

标题: 关于读取指定行 [打印本页]

作者: wxy 时间: 2017-3-4 16:42
标题: 关于读取指定行
本帖最后由 wxy 于 2017-3-4 17:41 编辑

请问我想提取文件中的指定行数据，应该怎么改

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

program find_close

  implicit none

  integer m,i

  integer ::  sol , country

  integer , allocatable :: site(:)

  real , allocatable :: lat(:) , long(:)

  open(8,file="soil.sol")

  open(10,file="find_close.txt")

  !read(8,*)m  !//m无法再文件中获取，需手动输入或获取文件行数减去3，得到。

    m=414503 !文件行数

  allocate(site(m),lat(m),long(m))



  do i = 1,m

    if(mod(i-3,9)==0)then

    read(8,*) site(i) , country , lat(i) , long(i) , sol

    end if

  enddo

  write(10,*) site(i) , country , lat(i), long(i) , sol

end program find_close

QQ截图20170304164156.png (40.32 KB, 下载次数: 536)

QQ截图20170304171959.png (36.88 KB, 下载次数: 529)

作者: vvt 时间: 2017-3-4 18:23
以下代码，会搜索文件中以 * 开头的行，并且读取它后面的第三行。（而其他无用的行有多少，无关精要）

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

program find_close

  implicit none

  integer m,i,k

  integer ::  sol , country

  integer , allocatable :: site(:)

  real , allocatable :: lat(:) , long(:)

  character(len=32) :: c

  open(8,file="soil.sol")

  open(10,file="find_close.txt")

  !read(8,*)m  !//m无法再文件中获取，需手动输入或获取文件行数减去3，得到。

  m=414503 !//m 此时的含义为，有多少个数据（而不是有多少行）

  allocate(site(m),lat(m),long(m))

O:Do i = 1 , m

    do 

      Read(8,'(a)',iostat=k) c

      if ( k/=0 ) exit O

      if ( c(1:1) == '*' ) exit

    end do

    Read(8,*)

    read(8,*) site(i) , country , lat(i) , long(i) , sol

    write(10,*) site(i) , country , lat(i), long(i) , sol

  End Do O

end program find_close

作者: wxy 时间: 2017-3-4 19:49
本帖最后由 wxy 于 2017-3-4 20:25 编辑

vvt 发表于 2017-3-4 18:23
以下代码，会搜索文件中以 * 开头的行，并且读取它后面的第三行。（而其他无用的行有多少，无关精要）
[mw_ ...

谢谢我试了一下，出现整型错误，请问是不是integer,allocatable :: site(:)这里有问题因为文件中对应的是CHN+一串数字，我其实主要想提取出坐标经纬度lat、long和编号site ，这种情况应该改呢

QQ截图20170304202207.png (38.31 KB, 下载次数: 435)

作者: vvt 时间: 2017-3-4 20:18
请上传一段数据文件，以便调试

作者: wxy 时间: 2017-3-4 20:24
本帖最后由 wxy 于 2017-3-4 21:28 编辑

vvt 发表于 2017-3-4 20:18
请上传一段数据文件，以便调试

这是其中一段，谢谢。

soil.txt

15.69 KB, 下载次数: 12

作者: wxy 时间: 2017-3-4 22:28

vvt 发表于 2017-3-4 20:18
请上传一段数据文件，以便调试

上传了一小段麻烦帮看一下

作者: vvt 时间: 2017-3-4 22:31
是的，的确是因为类型不对。
"CHN003830" 不能直接被转换成 integer
"LO" 也是如此，如果你不需要它，可以不读它。

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

program find_close

  implicit none

  integer m,i,k

  character(len=12) ::  country

  character(len=30) , allocatable :: site(:)

  real , allocatable :: lat(:) , long(:)

  character(len=32) :: c

  open(8,file="soil.txt")

  open(10,file="find_close.txt")

  !read(8,*)m  !//m无法再文件中获取，需手动输入或获取文件行数减去3，得到。

  m=414503 !//m 此时的含义为，有多少个数据（而不是有多少行）

  allocate(site(m),lat(m),long(m))

O:Do i = 1 , m

    do 

      Read(8,'(a)',iostat=k) c

      if ( k/=0 ) exit O

      if ( c(1:1) == '*' ) exit

    end do

    Read(8,*)

    read(8,*) site(i) , country , lat(i) , long(i)

    write(10,*) trim(site(i)) , trim(country) , lat(i), long(i)

  End Do O

end program find_close

作者: wxy 时间: 2017-3-4 23:15
本帖最后由 wxy 于 2017-3-4 23:43 编辑

vvt 发表于 2017-3-4 22:31
是的，的确是因为类型不对。
"CHN003830" 不能直接被转换成 integer
"LO" 也是如此，如果你不需要它，可以 ...

嗯嗯谢谢成功了

作者: wxy 时间: 2017-3-5 01:05

vvt 发表于 2017-3-4 22:31
是的，的确是因为类型不对。
"CHN003830" 不能直接被转换成 integer
"LO" 也是如此，如果你不需要它，可以 ...

还想请教一下另一个问题我想实现输入一个经纬度跟提取到的坐标相匹配(如果刚好数值完全一样就匹配,如果没有就寻找最近的坐标)，并把坐标相对应的编号输出也就是 “CHN003830”但我想不出怎么实现寻找最近的点

QQ截图20170305010259.png (10.06 KB, 下载次数: 514)

作者: vvt 时间: 2017-3-5 09:28
这个问题就复杂了。
如果数据量不大，可以用数组存储，然后查找（关于距离，因为有乘方的关系，所以计算速度比较慢，可以只计算 X 的差和 Y 的差。效果也差不多）
如果数据量大，可以先对数组进行一次排序，然后二分法查找。

作者: li913 时间: 2017-3-5 12:07
如果满足下面的要求，可以多费点功夫做数据前期处理：
1、坐标数据量大，且短时间内不会改变；
2、频繁调用该功能；
3、对计算速度要求很高。

QQ截图20170305115617.png (24.6 KB, 下载次数: 354)

作者: wxy 时间: 2017-3-5 13:08
本帖最后由 wxy 于 2017-3-5 14:03 编辑

li913 发表于 2017-3-5 12:07
如果满足下面的要求，可以多费点功夫做数据前期处理：
1、坐标数据量大，且短时间内不会改变；
2、频繁调用 ...

我现在正学着排序，然后我在想文件中经纬度是固定对应的，我的代码里把经度纬度分别放进两个动态数组里，如果这样排序是不是对应关系就打乱了，是不是可以把经纬度两个一维数组放进一个二维数组里再排序，直接把经纬度用二维数组进行读取，不知道我的想法对不对

作者: wxy 时间: 2017-3-5 13:58
本帖最后由 wxy 于 2017-3-5 14:01 编辑

vvt 发表于 2017-3-5 09:28
这个问题就复杂了。
如果数据量不大，可以用数组存储，然后查找（关于距离，因为有乘方的关系，所以计算速 ...

之前的代码是把经纬度分别放入两个数组，我想把经纬度放入一个二维动态数组再排序，但二维动态数组还不太会用出现了图片中的错误，麻烦帮我看一下怎么改

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

program near

  implicit none

  integer m,i,k

  character(len=12) ::  country

  character(len=30) , allocatable :: site(:)

  real , allocatable ::  array(:,:)

  character(len=32) :: c

  open(8,file="soil.sol")

  open(10,file="find_close.txt")

  !read(8,*)m  !//m无法再文件中获取，需手动输入或获取文件行数减去3，得到。

  m=414503 !//m 此时的含义为，有多少个数据（而不是有多少行）

  allocate(array(m,m),site(m))

O:Do i = 1 , m

    do

      Read(8,'(a)',iostat=k) c

      if ( k/=0 ) exit O

      if ( c(1:1) == '*' ) exit

    end do

    Read(8,*)

    read(8,*) site(i) , country ,array(i,i)





    write(10,*) trim(site(i)) , array(i,i)

  End Do O

end program near

QQ截图20170305140026.png (35.69 KB, 下载次数: 389)

作者: li913 时间: 2017-3-5 16:49
array数组太大，600多G，内存不够用。allocate(array(2,m),site(m))

作者: wxy 时间: 2017-3-5 17:33
本帖最后由 wxy 于 2017-3-5 22:16 编辑

li913 发表于 2017-3-5 16:49
array数组太大，600多G，内存不够用。allocate(array(2,m),site(m))

嗯嗯成功了谢谢开始学着弄下一步

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