如何在 Python 中围绕特定像素旋转图像？我正在尝试对一组夜空图像进行反旋转。由于星星围绕北极星旋转，我可以将北极星定义为旋转中心并旋转每个图像以排列星星。

我已经查看了使用 ehem 计算日出和日落的示例，并且效果很好。

当我尝试计算这两次之间的中点时，我遇到了麻烦。这是我所拥有的：

我得到：

日出：2010/11/2 12:47:40 UTC
日落：2010/11/2 23:24:25 UTC
中午：5:18:22.679044

这就是我卡住的地方。我是一个 python 初学者，坦率地说，我不是一个程序员。

任何建议都将受到欢迎！

是否有一个 Objective-C（或 C）库（与核心位置兼容）可以告诉我任何给定日历日的日出和日落时间？

我想计算给定日期和时间的地球位置（相对于太阳）和轴旋转。可以假设太阳在 0,0,0 坐标处是静止的。例如，由于月球引力引起的非常小的偏转也可以忽略不计。任何准确度都在一个度左右就足够了。

是否有任何库/源/数据可以帮助我完成此任务？

问候，

我正在进行研究，这将有助于确定观察到的空间的大小以及自大爆炸以来经过的时间。希望您能提供帮助！

我有双线性数据，我想对其执行两个线性回归。有一个斜率和截距发生变化的点，我需要（编写一个程序来）找到这个点。

想法？

~ 大约太阳正午

作为更新，部分混淆是太阳正午是自公元前 4713 年 1 月 1 日格林威治中午以来所有计算开始的地方。

Date.ordinal_to_jd 的正确使用并没有弥补这一事实。因此，通过像这样添加或减去 12 小时：

我们应该得到更少的错误。既然我们的计算是从昨天中午开始的，那么我们使用哪个呢？

该代码来自此页面Sunrise_equation的两个方程式。

我从这里的一位用户那里得到的第一个答案是，我们不了解 0.0009 和 lw / 360 的使用。lw / 360 似乎是本初子午线的小数日。至于 0.0009，它必须是自公元前 4713 年 1 月 1 日格林威治中午以来以秒为单位的少量方差。有关更多信息，请参阅 IAU 标准

根据此页面，我计算为 0.007776 秒。

我有一些来自 Date 类的信息，不包括方法细节。

附带说明一下，Ruby 有一种计算朱利安日期的方法真是太好了。我正在研究NOAA的 Javascript 代码。

这是我受到链接启发而写的一堂课。

现在这让我到了那里，但我仍在研究如何返回一个数字，例如由最上面的方程计算的数字。试试这个，我们可以看到我们在 JDN 中确实得到了 0.5。谁是对的？红宝石还是 NOAA？

NOAA 使用从 jd 中减去的 2000 年 1 月 1 日值 2451545.0 来获得像这样的分数世纪时间

我的任务相对简单：对于输入文件中的每一行，测试该行是否满足给定的一组条件，如果满足，则将该行的特定列写入新文件。我已经编写了一个执行此操作的 python 脚本，但我需要一些帮助：1）提高速度，2）在列名方面的最佳工作方式（因为列号可能因文件而异），以及 3 ) 指定我的过滤条件和所需输出列的最佳方式。

1）我使用的文件包含天文图像的光度测量。每个文件大约 1e6 行乘 150 列浮点数，通常大小超过 1GB。我有一个旧的 AWK 脚本，它将在大约 1 分钟内处理这样的文件；我的 python 脚本需要 5 到 7 分钟。我经常需要调整过滤条件并重新运行几次，直到输出文件是我想要的，所以速度绝对是可取的。我发现 for 循环非常快。这就是我在循环中做事的方式，这会减慢它的速度。与将整行读入内存相比，使用 itemgetter 只选择我想要的列是一个很大的改进，但我不确定我能做些什么来进一步提高速度。这能和 AWK 一样快吗？

2）我想根据列名而不是列号来工作，因为特定数量（光子计数、背景、信噪比等）的列号可以在文件之间更改。在我的 AWK 脚本中，我总是需要检查指定条件和输出列的列号是否正确，即使过滤和输出适用于相同的数量。我在 python 中的解决方案是创建一个字典，为每个数量分配一个列号。当一个文件有不同的列时，我只需要指定一个新的字典。也许有更好的方法来做到这一点？

3) 理想情况下，我只需要指定输入和输出文件的名称、过滤条件和要输出的所需列，它们会在我的脚本顶部找到，所以我不需要去搜索代码只是为了调整一些东西。我的主要问题是未定义的变量。例如，典型条件是“SNR > 4”，但在开始从测光文件中读取行之前，实际上并未为“SNR”（信噪比）分配值。我的解决方案是使用字符串和 eval/exec 的组合。再次，也许有更好的方法？

我根本没有接受过计算机科学的培训（我是天文学的研究生）——我通常只是把一些东西拼凑在一起并调试直到它起作用。但是，针对我以上三点的优化对我的研究来说变得非常重要。我为这篇冗长的帖子道歉，但我觉得这些细节会有所帮助。除了清理/编码风格之外，您对我的任何和所有建议都将不胜感激。

非常感谢，杰克

我正在尝试求解开普勒方程，作为找到给定时间的轨道物体真正异常的一步。但事实证明，开普勒方程很难求解，维基百科页面描述了使用微积分的过程。好吧，我不知道微积分，但我知道求解方程涉及无限数量的集合，这些集合产生越来越接近正确答案的近似值。

从数学上看我看不出如何通过计算来做到这一点，所以我希望有更好的数学背景的人可以帮助我。我怎样才能通过计算解决这个野兽？

FWIW，我正在使用 F#——我可以计算这个方程所需的其他元素，这只是我遇到问题的部分。

我也对在给定时间、近点距离和偏心率的情况下近似真实异常的方法持开放态度

我正在尝试构建一个天文学应用程序，我应该在给定的纬度、经度和时间显示天顶处和附近的星星。

我知道位置天文学的基础知识，令我困惑的是如何从大量数据中创建数据库！

我的查询是：

1. 更喜欢哪个目录？喜欢 Hipparcos 目录或 Henry Draper？它的优点和缺点是什么？
2. 需要哪些字段？我已经确定了一些，但还应该添加什么，为什么？我识别的字段：星名、一些唯一的星号、RA、Dec、星等、星的颜色。
3. 我试图找出从 Hipparcos 目录中获取Common Star 名称的表，但没有任何运气：http: //vizier.u-strasbg.fr/cgi-bin/VizieR-2 ?-source=I/239 。我在任何一张桌子上都找不到。

谢谢，拉吉

我正在使用 PyEphem 并想计算阴影的长度（假设在地下种植了一根单位长度的棍子）。长度将由 cot(phi) 给出，其中 phi 是太阳仰角（如果我错了，请纠正我）。我不确定在太阳上使用什么场？在下面的示例中，我使用角度 alt：

请在下面查看我的解释：

1. 如果切线是无限的，则表明太阳直接在头顶，没有阴影。
2. 如果切线为零，则表明太阳在地平线上，阴影无限长。
3. 我不知道如何解释 cot(phi) 的负面结果。有人能帮我吗？

最后，在给定 ephem.Observer() 的情况下，我对如何使用 PyEphem 从阴影长度向后工作到下一次太阳将投射该长度的阴影感到困惑。

我会很感激这方面的帮助。