问题标签 [astronomy]

我有一组用三脚架上的相机拍摄的许多天文照片。使用气泡水平仪确保框架的长边与地平线平行，并且我知道每张照片中心的 alt/az（和赤道）坐标。

现在我正在编写一些 python 代码来在每个图像上覆盖一个指示器以标记北方向。对于给定的 alt/az 坐标，我可以使用 pyephem 获取北天极方向和水平方向之间的角度吗？有什么线索吗？

这段代码基本上是遍历数组 ( ra_new) 的长度，并创建图像图（由下面编码的方程式给出）。因为数组真的很长（因此有很多很多的图），我试图让 Python 将它们保存为易于阅读的通用名称。换句话说，如果我输入的文件名（为简单起见，我没有显示该代码）是“googlypants”，我希望代码循环遍历该文件并将文件保存为“GooglypantsIMAGE1”。我还需要将其保存为 png 或 jpeg。

有人可以告诉我我的错误在哪里吗？目前，代码保存为“GooglypantsIMAGEi”（其中 i 是 arra 中的元素），这当然是一个并不真正存在的文件扩展名。

I'm basically trying to plot some images based on a given set of parameters of a .fits file. However, this made me curious: what IS a .fits array? When I type in img[2400,3456] or some random values in the array, I get some output.

I guess my question is more conceptual than code-based, but, it boils down to this: what IS a .fits file, and what do the arrays and the outputs represent?

从昨天开始，我一直在寻找答案，但没有运气。所以我有一个一维光谱（.fits）文件，每个波长都有通量值。我已将它们转换为二维数组（x，y）=（波长，通量）并想编写一个程序，该程序将在某些指定的波长（x）处返回通量（y）。我试过这个：

但即使波长不同，程序也总是返回相同的通量值。请帮忙...

我的大部分代码都采用 .fits 文件并创建基于某些参数的小缩略图（它们是星系的图像，所有这些都是无关信息......）

无论如何，我设法找到一种将图像保存为 .pdf 的方法，但我不知道如何将它们保存为 .fits 文件。解决方案需要在“for”循环中，以便它可以批量保存文件，因为缩略图太多，无法逐个迭代。

最后两行是最相关的。

有什么想法吗？

我目前正在与一个业余火箭手团队（好吧，“业余”是一个延伸，我们所有人都已经完成、注册或申请了相关的大学学位课程）开展了一个将卫星送入轨道的项目。到目前为止，我主要完成了高级编程（VB.NET、Python、PHP ），对FIG-FORTH有一点涉猎，我想知道学习和使用什么语言来编写姿态控制/航空电子系统的程序运载飞船。

当然，由于速度和存储方面的效率在这里都是一个大问题，所以想到的一件事就是组装。我的问题：这样做是否可行，即在装配中对自动姿态控制系统进行编程？

用这样一种极低级的编程语言来做这件事似乎非常复杂，所以我想知道在该领域有更多经验的人是否可以权衡一下汇编是否是一个好的选择，或者仍然是低级的东西但与FORTH、FORTRAN或C相比，它们在程序大小和速度方面所招致的损失相比，处理起来稍微容易一些。

I have few tens of full sky maps, in binary format (FITS) of about 600MB each.

For each sky map I already have a catalog of the position of few thousand sources, i.e. stars, galaxies, radio sources.

For each source I would like to:

• open the full sky map
• extract the relevant section, typically 20MB or less
• run some statistics on them
• aggregate the outputs to a catalog

I would like to run hadoop, possibly using python via the streaming interface, to process them in parallel.

I think the input to the mapper should be each record of the catalogs, then the python mapper can open the full sky map, do the processing and print the output to stdout.

1. Is this a reasonable approach?
2. If so, I need to be able to configure hadoop so that a full sky map is copied locally to the nodes that are processing one of its sources. How can I achieve that?
3. Also, what is the best way to feed the input data to hadoop? for each source I have a reference to the full sky map, latitude and longitude

我最近下载并安装了 PySolar 模块，并想用它来计算太阳位置（方位角和高度）。

我遇到的问题是计算的结果与 NOAA 太阳能计算器计算的结果不匹配http://www.esrl.noaa.gov/gmd/grad/solcalc/azel.html

例如。使用 PySol

使用 NOAA 模型高度计算为 0.7，方位角计算为 237.64。

因此，这两个值之间存在差异。总结一下：

• Pysolar 计算海拔高度为 -63.03，NOAA 计算海拔高度为 0.7
• Pysolar 计算方位角为 -235.44，NOAA 计算方位角为 237.64

如果您对我做错了什么或如何计算给定位置的正确高度和方位角有任何建议，我们将不胜感激。

我到处寻找，但没有找到解决这个问题的方法。我觉得它应该相当简单，但我们会看到。

我有一个 .FITS 格式的数据立方体，我需要将其折叠成 2D FITS 图像。数据立方体有两个空间维度和一个光谱/速度维度。

只是寻找一个简单的python 例程来加载立方体并展平所有这些层（即沿光谱/速度轴整合它们）。谢谢你的帮助。

我一直在疯狂地搜索文档，但找不到这个答案。

我在 python 中生成 FITS 图像，需要为图像分配 WCS 坐标。我知道有很多方法可以通过将点源与已知目录进行匹配来做到这一点，但在这种情况下，我正在生成尘埃图，因此点源匹配不起作用（据我所知）。

所以图像是一个形状为 (240,240) 的 2D Numpy 数组。它是这样写的（x 和 y 坐标分配有点奇怪，它以某种方式工作）：

这一切都很好。分配银河坐标似乎你需要做的就是：

但是我不明白FITS图像是如何存储这些坐标的，所以不知道怎么写。我也尝试在 SAO DS9 中分配它们，但没有运气。我只需要一种将这些坐标分配给图像的简单方法。

感谢您的任何帮助，您可以提供。