更新:Skyfield刚刚进行了重大修订,包括扩展文档和角度分离方法 – 请参阅接受的答案.

我正在使用Skyfield来计算两个物体之间的明显角度间隔.我没有在包装内找到方法,所以我通过计算两个明显位置矢量之间的点积来“发明”了一种方法.

这是目前最好的方法吗?在Skyfield的范围内它是否基本正确?

def separation(seconds, lat, lon):

lat, lon, seconds = float(lat), float(lon), float(seconds) # necessary it seems

place = earth.topos(lat, lon)

jd = JulianDate(utc=(2016, 3, 9, 0, 0, seconds))

mpos = place.at(jd).observe(moon).apparent().position.km

spos = place.at(jd).observe(sun).apparent().position.km

mlen = np.sqrt((mpos**2).sum())

slen = np.sqrt((spos**2).sum())

sepa = ((3600.*180./np.pi) *

np.arccos(np.dot(mpos, spos)/(mlen*slen)))

return sepa

from skyfield.api import load, now, JulianDate

import numpy as np

from scipy.optimize import minimize

data = load('de421.bsp')

sun = data['sun']

earth = data['earth']

moon = data['moon']

sep = separation(12000, 32.5, 215.1)

print "sun-moon aparent separation: ", sep, " arcsec"

解决方法:

Skyfield现在支持一种方法,该方法返回两个位置之间的角度间隔:

标签:python,skyfield

python中arcsec_python – 更好的方法来计算Skyfield中两个物体的明显角度分离?相关推荐

  1. python中偶数怎么表示_Python程序可计算列表中的偶数和奇数

    在本文中,我们将学习下面给出的问题陈述的解决方案. 问题陈述-我们得到一个列表,我们需要计算列表中的偶数和奇数. 共有三种方法,如下 - 方法1:使用蛮力方法 示例list1 = [21,3,4,6, ...

  2. rms c语言 函数,C中任何更快的RMS值计算?

    我正在用C编写一个小型8位微控制器的软件.部分代码是读取电流互感器(ZCT)的ADC值,然后计算RMS值.流过ZCT的电流是正弦曲线但可能会失真.我的代码如下: float adc_value,ins ...

  3. java中如何使用反射调用方法以及获得类中的属性

    使用反射获取类中的方法: 1):获取方法所在类的字节码对象. 2):获取方法. ------------------------------------------------------------ ...

  4. 8种在JavaScript数组中查找指定元素的方法(用于开发中数据的处理)

    1.Array.prototype.includes() includes() 方法用来判断一个数组是否包含一个指定的值,如果包含则返回 true,否则返回 false. 该方法支持两个参数value ...

  5. php中创建函数的正确方法,如何在PHP中创建一个函数

    我正在尝试创建一个名为saveorder的函数.到目前为止这是我的代码: //function foo($arg_1, $arg_2, /* ..., */ $arg_n) function save ...

  6. 有限元方法编程计算电场中各个区域分布的电势

    用有限元方法解二维的拉普拉斯方程,求出有边缘二维面上的电势分布 这里以计算10x10区域为演示,当差异小于设置值后停止演算 Code: #include <QCoreApplication> ...

  7. c语言中除号和百分号的用法,注册会计师计算题中能打百分号吗

    注册会计师财务管理机考想要输入x的平方,就可以输入x^2 ^是shift+,输入^的时候鼠标的光标会跳到上标位置,别用鼠标点上标,输入完成之后,可以按Tab键回到正常的位置. 扩展资料: 其他常见快捷 ...

  8. Python精讲:在Python中遍历字典的三大方法详解

    欢迎你来到站长在线的站长学堂学习Python知识,本文学习的是<在Python中遍历字典的三大方法详解>.本知识点主要内容有:使用字典对象的items()方法可以遍历字典的项和字典的&qu ...

  9. 6种在 Python 中从 List 中删除重复项的方法

    来源 | https://medium.com/@cookbug/six-ways-to-remove-duplicates-from-list-in-python-970d998b1384 翻译 | ...

  10. python string length_如何使用python获取字符串长度?哪些方法?

    掌握多种python技巧,对于我们更好的灵活应用python是非常重要的,比如接下来给大家介绍的获取字节长度,那大家脑海里就该有印象了,有几种方法呢?一起来看下吧~ 1.使用len()函数 这是最直接 ...

最新文章

  1. python中json dumps_python中json.loads,dumps,jsonify使用
  2. LHC大神问的矩阵转置问题
  3. highcharts图表组件常见问题:highcharts图表组件错误集合分析大放送
  4. lambdafiltermapreduce函数的基本使用
  5. 简单01背包 POJ3211 Washing Clothes 多种衣服分别dp
  6. NET 2.0中WinForm自定义的程序配置存放到哪里去了
  7. 机房收费--主界面(控件随窗体成比例改变)
  8. dd命令测试linux磁盘io情况,【LINUX】正确的使用dd进行磁盘读写速度测试
  9. linux主从库配置文件,linux 数据库主从同步配置
  10. 为什么学Python
  11. 栈-----括号匹配+表达式计算
  12. SCCM 2007 R7使用手记
  13. bios sgx需要开启吗_Win10改Win7,UEFI改Legacy 启动,修改BIOS大全
  14. 记一次云服务器ssh爆破维护实例
  15. 至简设计系列_7段数码管显示
  16. python网络数据采集 第二版_Python网络数据采集 (影印版)第2版
  17. ps如何把自己的图与样机结合_如何用PS制作样机,并且看起来毫无痕迹
  18. 基于Android Studio开发的笔记APP
  19. 自制万能xp镜像让重做系统变得简单
  20. android雪花飘落效果,【OpenGL】Shader实例分析(七)- 雪花飘落效果

热门文章

  1. [CSS备忘] css3零散
  2. 查询工商单位注册信息 平台
  3. 类的加载过程详解之过程五:类的Unloading(卸载)
  4. MySQL中exists和in的区别
  5. 面试题之请描述一下Java类加载过程
  6. 15张图带你彻底明白spring循环依赖,再也不用怕了
  7. hashCode()和equals()的若干问题
  8. SpringBoot中的定时任务和异步任务
  9. springcloud整合sentinel
  10. 呼吁的词汇用过“人联网晚一天应用就是人类的重大损失”