请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 118|回复: 0

[混合编程] IDL之由方向和距离求坐标(二分逼近法)

[复制链接]

新浪微博达人勋

发表于 2020-6-30 14:58:49 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
问题由来:从已知坐标往X方向L米处的坐标是多少?

目前,我只知道距离怎么求,求坐标没有现成函数,于是我利用二分法来逼近精确坐标。
思路是,我先选坐标原点和大约1万米外的点(近似为经纬度0.1)为直线两端,
然后选择端点的中间点为新端点,
然后新端点和旧端点继续折中出新端点,......,
直到其中一个端点与原点的距离在L±1m的精度范围内。

000.png
代码与解释:

;vec是单位矢量,代表方向,笛卡尔坐标,从正东逆时针,[0,1]表示正北
function neighbor,core,vec,distance

  ;不确保为双精度,返回值很可能是6位有效数字
  x1=double(core[0]);经度
  y1=double(core[1]);纬度
  x2=x1+0.1*vec[0];以大约1万米外的位置为初始“远端点”
  y2=y1+0.1*vec[1]
  mx=(x1+x2)/2;折中点
  my=(y1+y2)/2
  
  constraint=1;精确度,单位米,不达到这个精度循环不停止

  repeat begin
    d1=map_2points(x1,y1,core[0],core[1],/METERS);近端点与原点的距离
    d2=map_2points(x2,y2,core[0],core[1],/METERS);远端点与原点的距离
    ;如果折中点较远,则折中点为新的远端点:
    if(map_2points(mx,my,core[0],core[1],/METERS) gt distance+constraint)then begin
      x2=mx
      y2=my
      mx=(x1+x2)/2
      my=(y1+y2)/2
    endif
    ;如果折中点较近,则折中点为新的近端点:
    if(map_2points(mx,my,core[0],core[1],/METERS) lt distance+constraint)then begin
      x1=mx
      y1=my
      mx=(x1+x2)/2
      my=(y1+y2)/2
    endif
    ;如果远近2个端点其中之一与原点的距离在误差范围内就结束循环
  endrep until (abs(d1-distance) le constraint or abs(d2-distance) le constraint)


  d1=map_2points(x1,y1,core[0],core[1],/METERS)
  d2=map_2points(x2,y2,core[0],core[1],/METERS)
  x=d1 gt d2?x2:x1;二端点近者为返回值
  y=d1 gt d2?y2:y1
  return,[x,y]
end

测试:
IDL> neighbor([120.0,30.0],[0,1],500)
% Compiled module: NEIGHBOR.
% Compiled module: MAP_2POINTS.
       120.00000000000000       30.004492187566939

该算法的python版见《python之由方向和距离求坐标(二分逼近法)》:
http://bbs.06climate.com/forum.php?mod=viewthread&tid=96382&extra=
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

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