求助，利用fortran的数据选取问题

river

言深深 发表于 2012-10-16 18:40

                               
登录/注册后可看大图

顶个板凳
注：x(i,j,1)是经度x(i,j,2)是纬度,x(i,j,3)，ab分别是定点经纬度,r是半径
do i=1,m

很牛叉啊，这都想的出来

大鹏

言深深 发表于 2012-10-16 23:24

                               
登录/注册后可看大图

只是做一个粗略的意思，需要进行转换的，有公式的

能不能把公式给下啊

言深深

大鹏 发表于 2012-10-17 18:19

                               
登录/注册后可看大图

能不能把公式给下啊

额，网上到处都是这样的公式啊
http://zhidao.baidu.com/question ... bc9da890fcc3126f647

“
任意两点距离计算公式为

d＝111.12cos{1/[sinΦAsinΦB十cosΦAcosΦBcos(λB—λA)]}

其中A点经度，纬度分别为λA和ΦA，B点的经度、纬度分别为λB和ΦB，d为距离。
”

SEAMAN

言深深 发表于 2012-10-16 23:24

                               
登录/注册后可看大图

只是做一个粗略的意思，需要进行转换的，有公式的

这个跟我的做法很有相似性呢，科学都是想通的

言深深

SEAMAN 发表于 2013-3-28 15:57

                               
登录/注册后可看大图

这个跟我的做法很有相似性呢，科学都是想通的

触类旁通却是是科研需要素质啊···

溫不暖

言深深 发表于 2012-10-17 19:50

                               
登录/注册后可看大图

额，网上到处都是这样的公式啊
http://zhidao.baidu.com/question/32534955.html&__bd_tkn__=70f8432979 ...

那个网页上的公式不对.....
我查了一下这个公式比较说的通
http://zhidao.baidu.com/question ... l=relate_question_2
R*{arccos[cosb*cosy*cos(a-x)+sinb*siny]}
特别注意要换算弧度和角度
fortran写出来大概是这样计算的
p=3.1417/180
r=6370(近似的地球半径)
d=r*acos(sin(lat1*p)*sin(lat2*p)+cos(lat1*p)*cos(lat2*p)*cos(lon2*p-lon1*p))