Python中文
首页
教程
问答
标签
搜索
登录
注册
Python计算两个大圆的交点
回答此问题可获得
20
贡献值,回答如果被采纳可获得
50
分。
<p>我试图计算两个大圆的交点(纬度和经度),这两个大圆都是由圆上的两点定义的。我一直在尝试遵循概述的方法<a href="https://stackoverflow.com/questions/2954337/great-circle-rhumb-line-intersectionp://">here</a>。 但是我得到的答案是不正确的,我的代码在下面有人知道我哪里出错了吗?在</p> <pre><code>################################################ #### Intersection of two great circles. # Points on great circle 1. glat1 = 54.8639587 glon1 = -8.177818 glat2 = 52.65297082 glon2 = -10.78064876 # Points on great circle 2. cglat1 = 51.5641564 cglon1 = -9.2754284 cglat2 = 53.35422063 cglon2 = -12.5767799 # 1. Put in polar coords. x1 = cos(glat1) * sin(glon1) y1 = cos(glat1) * cos(glon1) z1 = sin(glat1) x2 = cos(glat2) * sin(glon2) y2 = cos(glat2) * cos(glon2) z2 = sin(glat2) cx1 = cos(cglat1) * sin(cglon1) cy1 = cos(cglat1) * cos(cglon1) cz1 = sin(cglat1) cx2 = cos(cglat2) * sin(cglon2) cy2 = cos(cglat2) * cos(cglon2) cz2 = sin(cglat2) # 2. Get normal to planes containing great circles. # It's the cross product of vector to each point from the origin. N1 = cross([x1, y1, z1], [x2, y2, z2]) N2 = cross([cx1, cy1, cz1], [cx2, cy2, cz2]) # 3. Find line of intersection between two planes. # It is normal to the poles of each plane. L = cross(N1, N2) # 4. Find intersection points. X1 = L / abs(L) X2 = -X1 ilat = asin(X1[2]) * 180./np.pi ilon = atan2(X1[1], X1[0]) * 180./np.pi </code></pre> <p>我还应该提到这是在地球表面(假设是一个球体)。在</p>
0 条评论
分类:
Python问答
请先
登录
后评论
默认排序
时间排序
1 个回答
匿名
1天前
擅长:python、mysql、java
<p>DSM的解决方案在上面的注释中,您的角度是以度为单位,而sin和cos则期望弧度。 还有那条线</p> <pre><code>X1 = L / abs(L) </code></pre> <p>应该是</p> ^{pr2}$
请先
登录
后评论
针对此问题:
更多的回答
关注
89
关注
收藏
1
收藏,
216
浏览
网友 提问于 2天前
相关Python问题
如何添加虚拟方法
9 回答
如何添加表示整数的擦边字符串?
1 回答
如何添加要在Bokeh中使用的新font.ttf文件?
4 回答
如何添加要显示的矩阵XY轴编号和XY轴
2 回答
如何添加计数?
5 回答
如何添加计数器函数?
8 回答
如何添加计数器列来计算数据帧中另一列中的特定值?
3 回答
如何添加计数器来跟踪while循环中的月份和年份?
7 回答
如何添加计数并删除countplot的顶部和右侧脊椎?
1 回答
如何添加计时器wx.应用程序更新窗口对象的主循环?
7 回答
如何添加评论到帖子?PostDetailVew,Django 2.1.5
7 回答
如何添加评论拉梅尔亚姆
2 回答
如何添加诸如矩阵Python/Pandas之类的数据帧?
10 回答
如何添加谷歌地点自动完成到Flask?
4 回答
如何添加超时、python discord bot
6 回答
如何添加超过1dp的检查
8 回答
如何添加距离方法
5 回答
如何添加跟随游戏的敌人精灵
2 回答
如何添加路径以便python可以找到程序?
6 回答
如何添加身份验证/安全性以使用happybase访问HBase?
4 回答