Astropy matplotlib 和绘制星系坐标
我正在尝试使用 python 绘制星系坐标图.假设我有这些数据: data = [(0.261, -7.123, 13.03, 'Unidentified'), (-0.326, 77, 13.03, '星系')] 其中每个元组的形式为(ra、dec、flux、type). 我被要求使用 astropy + matplotlib，所以: c = SkyCoord(ra = ra
..
astropy相关内容
从 SFTP 服务器打开 Astropy FITS 文件
我有一个 Python 脚本，可以使用 Paramiko 模块通过 ssh 连接到远程服务器. 下面是我的脚本 import paramikossh = paramiko.SSHClient()ssh.set_missing_host_key_policy(paramiko.AutoAddPolicy())ssh.connect("host", username="McMissile")
..
如何使用多个嵌套的for循环加快python 2程序的速度
这段代码有多个 for 循环，我读入的列表每个都有 999 点.我想将此迭代最多10,000次.然而，即使只迭代 2 次也需要将近 10 分钟. 即使我发布了这个特定的代码，我认为对我的问题的回答可以帮助其他人更快地运行包含大量数据的代码. 您的任何建议都将受到赞赏.非常感谢. 此代码的作用:基本上，我从文本文件的列表中读取数组.每个列表(例如x1，y1，z1 ...等)每个都有
..
为什么WCS投影子图上的要素在matplotlib中放在错误的位置?
我有一个关于天文物体的拟合文件.我可以这样绘制: 来自astropy.io的 导入适合从 astropy.wcs 导入 WCShdul = fits.open(fitsfilename)[0]wcs = WCS(hdul.header)fig = plt.figure(figsize=(12,12))fig.add_subplot(111, 投影=wcs)plt.imshow(hdul.dat
..
使用hdf5文件中的astropy在skyplot上绘制均值和标准dev值
我正在尝试使用hdf5文件创建一个包含均值和标准dev值的Skyplot(使用astropy).数据链接为 https://wwwmpa.mpa-garching.mpg.de/~ensslin/research/data/faraday2020.html (法拉第天空2020年).到目前为止，我已经编写了以下代码，其中将数据从hdf5文件读取到ggl和ggb，然后将值转换为gb和gl中的银河坐标
..
AstroPy SkyCoord非常慢，该如何解决?
我正在使用AstroPy SkyCoord进行具有数百万个数据的赤道坐标到银河坐标的转换，这非常慢.任何人都有加速它的想法，否则将需要花费永远的时间来运行整个数据集.代码如下: 来自astropy导入单位的 作为u从astropy.coordinates导入SkyCoord将numpy导入为npra1 = np.loadtxt('data.txt'，usecols = {0，))dec1 =
..
如何从头开始生成拟合文件
在这篇文章中，他们解释了如何从ascii文件生成fits文件.但是，我也想知道如何将头文件和数据定义为fits文件.(将ASCII表转换为FITS图像) 例如，当我用astropy调用光谱拟合文件(从望远镜下载)时，可以分别调用数据和标头. IE 在[1]中:hdu = fits.open('observation.fits'，memmap = True)在[2]:header =
..
为什么说"ValueError:函数返回的数组在两次调用之间改变了大小"?
我正在尝试计算积分，但我不理解错误输出，它是什么意思，有什么问题以及如何解决? 输出为"ValueError:函数返回的数组在调用之间改变了大小"，我添加了完整的输出. 请帮助， 谢谢 import numpy as np import matplotlib.pyplot as plt import astropy.cosmology as cosmo from astropy import
..
在Casa中安装astropy
我正在尝试在Mac上将astropy安装到CASA 4.3.1中.我正在此处进行操作.但是运行pip.main(['install', 'astropy', '--user'])之后，我得到: AttributeError: 'module' object has no attribute 'main' 还有其他人遇到过这个问题吗? 编辑: 使用后 import pip
..
Astropy.model 2DGaussian问题
我正在尝试使2D高斯拟合图像以找到图像中最亮的对象.我正在尝试遵循此处的一维示例 http://docs.astropy .org/en/v0.3.2/modeling/index.html .我的代码看起来像 import numpy as np import astropy.io.fits as fits import os from astropy.modeling import mo
..
如何对可观察坐标列表进行排序？
我正在努力寻找从坐标列表中删除不需要的目标的最佳方法。我的坐标（Ra，12月）是使用 astropy.coordinates.SkyCoord 形成的，但是我有很多不可观测的目标，它们的磁偏角都太低，所以我想要要做的就是对我的列表进行排序，并删除所有偏角小于-10度的目标（例如，我的望远镜位于北半球）。 生成列表的代码行，为简化起见称为 radecs 并获得Ra& radecs =
..
如何在天文表中索引可观察的坐标
我问了这个问题（如何昨天通过排序可观察到的坐标列表进行排序吗？）关于排序坐标列表以删除阈值以下的某些值的问题，我从@MSeifert那里得到了一个很好的答案，但是我有一个表，其中这些坐标值与目标的其他属性（例如视在大小和Alt / Az坐标），所以我现在要问的是一种在 astropy.table 中进行掩蔽技术的方法比我上一个问题中的 astropy.coordinates.SkyCoord 列表
..
如何将以地球为中心的惯性（ECI）坐标转换为以地球为中心的地球固定（ECEF）AstroPy？其他？
我在卫星轨道的以地球为中心的惯性坐标（ECI）中具有位置（x，y，z）和速度（Vx，Vy，Vz）向量，最终希望以大地坐标（纬度，经度， & Altitude）。 根据此其他堆栈溢出问题，看来我需要转换为以地心固定地球（ECEF）坐标作为中间步骤（因此ECI-> ECEF- -> Lat / Lon / Alt）。 我知道ECI和ECEF共享相同的原点（地球质量中心），并且指向相同的z轴
..
APLpy/matplotlib:协调EPS质量图的网格Alpha级别
在普通的matplotlib轴类中，可以将网格线设置为具有一定的透明度(alpha级别).我正在尝试通过以下方法将其与APLpy软件包一起使用: fig = pyplot.figure(figsize=(18,8)) gridspec_layout = gridspec.GridSpec(1,2) gridspec_layout1 = gridspec.GridSpecFromSubpl
..
将标头复制到新文件中
我有这个脚本，可以使文件具有多个扩展名，但是我想将旧文件的标头添加到扩展名中. new_hdul = fits.HDUList() new_hdul.append(fits.PrimaryHDU(header=headermain)) new_hdul.append(fits.ImageHDU(nod1, header=header1, name='Chop1')) new_hdul.ap
..
如何将两个表与pyfits合并?
我正在使用Python 2.7.10和pyfits 3.3.之前，我使用以下代码合并两个表.但是，现在我遇到了一些错误 t1 = pyfits.open(table1)[1].columns t2 = pyfits.open(table2)[1].columns new_columns = t1 + t2 hdu = pyfits.BinTableHDU.from_columns(new_c
..
如何停止适合文件从颠倒的python导入
在python中导入拟合图像时，它会上下颠倒导入. 这是我正在使用的代码: import numpy as np from astropy.io get fits import matplotlib.pyplot as plt im = fits.getdata('myimage.fits') plt.imshow(im, cmap='gray') plt.colorbar()
..
在图像中绘制几个字段
如何在单个图像中绘制几个FITS字段?每个FITS文件都覆盖天空的相邻部分. HDU的data字段仅包含图像.因此，要在正确的坐标处绘制每个图像，我需要从header字段获取信息.但是我该如何在pyplot上传递这些信息? 我尝试使用astropy.WCS进行以下操作，以便将标头中的信息用于图表中 import matplotlib.pyplot as plt from astr
..
使用Astropy保存已编辑的.fits文件时，如何保存标题?
我正在编辑python中具有的.fits文件，但我希望标题保持完全相同.这是代码: import numpy as np from astropy.io import fits import matplotlib.pyplot as plt # read in the fits file im = fits.getdata('myfile.fits') header = fits.geth
..
固定符合FITS标准的卡意味着什么?
我正在尝试使用FITS文件.我有以下代码: from astropy.io import fits from astropy.wcs import WCS hdul = fits.open(fitsfilename)[0] wcs = WCS(hdul.header) 它给了我这些警告: 警告:Verify警告:验证报告的错误: [astropy.io.fits.verify]
..