问题标签 [astronomy]

在给定固定对象和观察者的情况下，如何确定 pyEphem 中目标的视差角？

我正在尝试使用 pyEphem 创建一个仪器模拟器，我需要确定仪器旋转器转动以使场旋转的速率（这是一个 alt-az 望远镜），所以我真的在尝试确定视差角的变化率。

谢谢，乔什

我正在尝试安装天文学图像软件 gaia。我正在运行 Snow Leopard 10.6.8 并安装了带有开发工具的 xcode 工具 3.2.6。我也相信我有必要的正确 gfortran 编译器。

我现在如何安装 gaia？我需要运行哪些命令。当我下载它时，它会解压缩并且只有一堆文件夹。

任何帮助是极大的赞赏。

谢谢，马克斯

我有一个 numpy 2D 直方图，它使用两个空间坐标（x_array 和 y_array）对数据点进行分类，并对每个 bin 中第三个 data_array 组件的值进行平均。

代码如下：

等等（x 和 y 在某些地方是颠倒的，我知道......这就是我让它工作的方式）。

我现在需要做的是对这些数据进行过采样（或网格化），即使用重叠的 bin 采样几次，这样整个采样就会移动，并且也向上（单独）移动半个 bin 宽度，以获得更平滑的图像。

我显然可以增加 bin 的数量来提高分辨率，但这只能在一定程度上起作用……在某个点上，有些 bin 没有任何数据点，因为它们是随机分布的。

我已经在这个上工作了一段时间，我只是被难住了。感谢您的任何帮助，您可以提供。

我有一个太阳能电池板指向某个方向（它是法线向量）。我想计算它与太阳当前位置之间的角度。我正在使用 pyephem，我有两对方位角和高度的信息。

找到面板法线向量和指向太阳的向量之间的角度的最简单方法是什么？

这是Josh r 代码的 php 实现，用于计算给定日期和时间的太阳位置：

这是 MvG help 之后的更正代码：

这已经在刚果（赤道附近）纬度/经度进行了测试：-4.77867 / 11.86364，日期为 2013 年 9 月 1 日 10:00。在这种情况下，正确答案是：仰角 = 67.77503 方位角 = 54.51532

感谢您帮助调试此 php 代码！

格雷格·法布尔。

有没有办法从分布在大多数 Linux 发行版上的时区数据库中提取历史闰秒时刻？我正在寻找 python 中的解决方案，但是在命令行上工作的任何东西都可以。

我的用例是在 gps 时间（基本上是自 1980 年第一颗 GPS 卫星开启以来的秒数）和 UTC 或本地时间之间进行转换。UTC 不时调整闰秒，而 gps-time 线性增加。这相当于在 UTC 和TAI之间转换之间进行转换。TAI 也忽略闰秒，因此 TAI 和 gps-time 应该始终以相同的偏移量演变。在工作中，我们使用 gps-time 作为全球同步天文观测的时间标准。

我有在 gps-time 和 UTC 之间转换的工作函数，但我必须硬编码一个闰秒表，我在这里得到了（该文件tzdata2013xx.tar.gz包含一个名为leapseconds）。每隔几年，当宣布新的闰秒时，我必须手动更新此文件。我更愿意从标准 tzdata 中获取此信息，该信息每年通过系统更新数次自动更新。

我很确定这些信息隐藏在 /usr/share/zoneinfo/. 我已经能够使用struct.unpack（man tzfile给出一些关于格式的信息）来提取其中的一些，但我从来没有让它完全工作。是否有任何标准软件包可以访问此信息？我知道pytz，它似乎从同一个数据库中获取标准 DST 信息，但它不能访问闰秒。我还找到了 tai64n，但查看它的源代码，它只包含一个硬编码表。

编辑

受到 steveha 的回答和pytz/tzfile.py中的一些代码的启发，我终于得到了一个可行的解决方案（在 py2.5 和 py2.7 上测试）：

结果

虽然这确实解决了我的问题，但我可能不会选择这个解决方案。相反，我将按照 Matt Johnson 的建议在我的代码中包含leap-seconds.list 。这似乎是用作 tzdata 来源的权威列表，并且可能每年由 NIST 更新两次。这意味着我将不得不手动进行更新，但是这个文件很容易解析并且包含一个过期日期（似乎缺少 tzdata）。

我正在使用 VSOP2000 数据文件，并且能够计算时间 T 的太阳和月亮的笛卡尔坐标......问题是如何将该坐标（在 J2000 的时间范围内）转换为当前日期（时间 T )? 我只是想知道算法。

我写了这个程序：

并得到了这些结果：

换句话说，月亮从极地（总是向上）到地平线以下，而实际上并没有升起或落下。

我意识到 89.5 度是一种极端情况，但为什么这段代码不能工作呢？

作为注释http://aa.usno.navy.mil/data/docs/RS_OneYear.php显示在此间隔内有一个月落（JD 2456623 = 2013-11-26 12:00:00）：

另一个奇怪的地方：如果我注释掉“$rst->get_transit();” 行，我得到以下结果：

上升/设置仍然很奇怪（这很好，因为状态是-1），但它们是不同的。我一直认为 get_lunar_rst() 返回一个固定的结构，但显然不是？检查运输时间的行为改变了结构？

编辑：好的，我运行了两次程序，根本没有改变它并得到：

换句话说，2456623 集无缘无故地发生了变化。

我想在星系的位置堆叠 WMAP 地图。是否可以使用healpy执行此操作？healpy 中的哪个函数可以让我得到一个点周围的天空投影？谢谢

是否可以使用 PyEphem 来计算：

1. 给定一个起始年份，找到主要/次要月球静止的下一年？
2. 使用这些日期来查找给定 X、Y 位置的月出/月落方位角？

我似乎无法在 PyEphem 文档中找到有关停顿的信息。如果 PyEphem 不能做到这一点，有没有人知道可以提供帮助的另一块或软件或库？