我需要使用xmatch从星象查询包交叉匹配一个大的本地目录与2 2MASS。像往常一样,我把我的本地适配表装上了不稳定的东西:
from astropy.io import fits
hdu = fits.open(root+'mycat.fits')然后,按照天体查询文档中描述的语法,尝试在该表(表为天体查询文档)中使用xmatch
from astroquery.xmatch import XMatch
table = XMatch.query(cat1=hdu[2],
cat2='vizier:II/246/out',
max_distance=1 * u.arcsec, colRA1='RA',
colDec1='Dec')但是,得到以下错误:
AttributeError: 'BinTableHDU' object has no attribute 'read'天体查询文档上的示例只展示了如何提供本地CSV文件。但是我的目录有大约700万个条目,所以把它作为一个ASCII CSV文件传递是不方便的。
我应该如何传递我的FITS表作为输入?谢谢!
发布于 2016-05-27 05:23:17
虽然xmatch可以接受一个文件对象作为输入,但该文件对象必须是一个Vizier样式的.csv表。你需要首先把你的FITS表转换成一个不稳定的表。
from astropy.table import Table
myTable = Table(data=hdu[2].data)https://stackoverflow.com/questions/37437692
复制相似问题