Fortran Coder

查看: 11598|回复: 7

[数值问题] 求两坐标点距离

[复制链接]

67

帖子

16

主题

0

精华

专家

F 币
275 元
贡献
201 点
发表于 2017-3-18 21:43:30 | 显示全部楼层 |阅读模式
本帖最后由 wxy 于 2017-3-19 12:41 编辑

我想实现输入一个坐标点跟文件中所有坐标点依次进行距离计算然后输出距离最近点,但是现在的代码出现图中的问题,计算出来全是是第一个坐标的值, 麻烦帮看一下应该怎么改
[Fortran] 纯文本查看 复制代码
real function distance(lat1, lon1, lat2, lon2)
real , parameter ::  R=6378.137
real lon1 , lat1 , lon2 , lat2
distance = (R*cos(lon1)*cos(lat1)+R*cos(lon2)*cos(lat2))/sqrt(sin(lon1)**2*sin(lon2)**2+sin(lat1)**2*sin(lat2)**2)
end function distance
program none
    implicit none
    integer i,m
    character(len=12):: site
    real :: lat , lon , lat1 , lon1 , distance
    real , allocatable ::  lat2(:) , lon2(:)

    open(7,file="find_close.txt")
    open(10,file="find.txt")
    m=46154 !行数
    allocate(lat2(m),lon2(m))

    Write (*,*) 'lon:'!输入经度
    Read (*, '(f7.3)') lon

    Write (*,*) 'lat:'!输入纬度
    Read (*, '(f6.3)') lat

    do i=1,m
        lon1=lon
        lat1=lat
        read(7,'(a9,2x,f6.3,2x,f7.3)') site , lat2(i) , lon2(i)
        write(10,'(a9,2x,f6.2,2x,f7.2,2x,f9.3)') site , lat2(i) , lon2(i) , distance(lat1, lon1, lat2, lon2)
    end do
end program
QQ截图20170319120522.png

新建文本文档.txt

8.12 KB, 下载次数: 2

63

帖子

9

主题

0

精华

专家

超凡脱俗

F 币
474 元
贡献
237 点
发表于 2017-3-19 12:51:13 | 显示全部楼层
[Fortran] 纯文本查看 复制代码
do i=1,m
        lon1=lon
        lat1=lat
        read(7,'(a9,2x,f6.3,2x,f7.3)') site , lat2(i) , lon2(i)
        write(10,'(a9,2x,f6.2,2x,f7.2,2x,f9.2)') site , lat2(i) , lon2(i) , distance(lat1, lon1, lat2(i), lon2(i))
    end do

QQ图片20170319124945.png
天下英雄出我辈,一入江湖岁月催。

鸿图霸业谈笑间,不胜人生一场醉。

67

帖子

16

主题

0

精华

专家

F 币
275 元
贡献
201 点
 楼主| 发表于 2017-3-19 14:12:28 | 显示全部楼层
本帖最后由 wxy 于 2017-3-19 14:14 编辑


再请教一下,我把计算出来的距离放入动态数组near,用minval找到距离最小值,如果要把最小值所在的那一行输出,应该怎么改呢,minloc和minval第一次用,
[Fortran] 纯文本查看 复制代码
real function distance(lat1, lon1, lat2, lon2)
real , parameter ::  R=6378.137
real lon1 , lat1 , lon2 , lat2
distance = abs(R*cos(lon1)*cos(lat1)+R*cos(lon2)*cos(lat2))/sqrt(sin(lon1)**2*sin(lon2)**2+sin(lat1)**2*sin(lat2)**2)
end function distance
program none
    implicit none
    integer i,m
    character(len=12):: site
    real :: lat , lon , lat1 , lon1 , distance , loc , nearmin
    real , allocatable ::  lat2(:) , lon2(:) , near(:)

    open(7,file="find_close.txt")
    open(10,file="near.txt")
    m=46154 !行数
    allocate(lat2(m),lon2(m),near(m))

    Write (*,*) 'lon:'!输入经度
    Read (*, '(f7.3)') lon

    Write (*,*) 'lat:'!输入纬度
    Read (*, '(f6.3)') lat

    do i=1,m
        lon1=lon
        lat1=lat
        read(7,*) site , lat2(i) , lon2(i)
        write(10,'(a9,2x,f6.2,2x,f7.2,2x,f11.2)') site , lat2(i) , lon2(i) , distance(lat1, lon1, lat2(i), lon2(i))
    end do
    close(10)
!======================================================================================================================
    open(10,file="near.txt")
    do i=1,m
        read(10,*)site , lat2(i) , lon2(i) , near(i)
        nearmin= minval(near(:))
        loc=minloc(near(:),1)
        print *, nearmin, loc
    end do

end program


QQ截图20170319141418.png

63

帖子

9

主题

0

精华

专家

超凡脱俗

F 币
474 元
贡献
237 点
发表于 2017-3-19 14:28:27 | 显示全部楼层
wxy 发表于 2017-3-19 14:12
再请教一下,我把计算出来的距离放入动态数组near,用minval找到距离最小值,如果要把最小值所在的那一行输出 ...

参考这个示例
QQ图片20170319142811.png
天下英雄出我辈,一入江湖岁月催。

鸿图霸业谈笑间,不胜人生一场醉。

67

帖子

16

主题

0

精华

专家

F 币
275 元
贡献
201 点
 楼主| 发表于 2017-3-19 14:57:54 | 显示全部楼层

谢谢 现在已经成功输出最小值的位置了,如果我想直接把最小值那一行输出,在知道行号的情况下,用什麽语句呢,以前只写过从第几行开始读,但没写过指定第几行输出,print *,怎样才能把指定行输出
[Fortran] 纯文本查看 复制代码
    open(10,file="near.txt")
    do i=1,m
        read(10,*)site , lat2(i) , lon2(i) , near(i)
    end do
        write(*,*)"最小值",minval(near(1:m))
        write(*,*)"最小值位置",minloc(near(1:m))
end program

QQ截图20170319144850.png

63

帖子

9

主题

0

精华

专家

超凡脱俗

F 币
474 元
贡献
237 点
发表于 2017-3-19 15:34:13 | 显示全部楼层
wxy 发表于 2017-3-19 14:57
谢谢 现在已经成功输出最小值的位置了,如果我想直接把最小值那一行输出,在知道行号的情况下,用什麽语句呢 ...

额,这个都找到了,你写到文件里的数据不都是按行,一行一行的排列好的吗,既然找到了想找的行数,直接把相应的数据输出就行
天下英雄出我辈,一入江湖岁月催。

鸿图霸业谈笑间,不胜人生一场醉。

67

帖子

16

主题

0

精华

专家

F 币
275 元
贡献
201 点
 楼主| 发表于 2017-3-19 15:48:23 | 显示全部楼层
Jackdaw 发表于 2017-3-19 15:34
额,这个都找到了,你写到文件里的数据不都是按行,一行一行的排列好的吗,既然找到了想找的行数,直接把 ...

嗯嗯 谢谢

67

帖子

16

主题

0

精华

专家

F 币
275 元
贡献
201 点
 楼主| 发表于 2017-3-19 22:34:40 | 显示全部楼层
本帖最后由 wxy 于 2017-3-19 22:49 编辑
Jackdaw 发表于 2017-3-19 15:34
额,这个都找到了,你写到文件里的数据不都是按行,一行一行的排列好的吗,既然找到了想找的行数,直接把 ...

额 最后一步还是做不出来,还是得麻烦您帮看一下第五十行是 k=minloc(near(1:m))
[Fortran] 纯文本查看 复制代码
    open(10,file="near.txt")
    do i=1,m
        read(10,*)site(i) , lat2(i) , lon2(i) , near(i)
    end do
        write(*,*)"最小值",minval(near(1:m))
        write(*,*)"最小值位置",minloc(near(1:m))
!======================================================================================================================
!输出最近坐标点所在行
!======================================================================================================================
    k=minloc(near(1:m))
    do i=1,m
        if (i==k)then
            read(10,'(a9,2x,f6.2,2x,f7.2,2x,f11.2)') site(i) , lat2(i) , lon2(i) , near(i)
            print *, site(i) , lat2(i) , lon2(i) , near(i)
            exit
        end if
    end do
end program

QQ截图20170319223400.png
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-6-24 18:32

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表