如何使用Astropy/Astroquery获取Sloan震级

Question

我想要创建一个脚本，从一组对象中下载量值。例如，出现在这里的：

http://mirror.sdss3.org/spectrumDetail?plateid=556&mjd=51991&fiber=312

作为输入，我有对象坐标(这是从天体查询get_spec方法获得的，因为我知道对象的mjd、板块和光纤)。我试图在天体查询站点上使用这个例子：

from astropy import coordinates as coords
from astroquery.sdss import SDSS

co = coords.SkyCoord(143.50993, 55.239775, unit="deg")
result = SDSS.query_region(co)
imgs = SDSS.get_images(co, band = ['g', 'r'])

但是，从下载的图像中我找不到震级。我怎样才能找到我的物体的大小？

Answer 1

相反，从事物的声音来看，您希望对象的属性从SDSS管道中抛出，并由坐标进行交叉匹配。

要做到这一点，请看一下SDSS.query_crossid方法。这允许您指定所需的测光列中的哪一列。

由于我不知道你正在做的科学，你必须找出什么对你的目标是最好的。请看一下这里所提供的各种量值。

举个例子，如果我想要g和r波段的模型震级作为对象，我会这样做：

In [31]: result = SDSS.query_crossid(co, photoobj_fields=['modelMag_g', 'modelMag_i'])

In [32]: print result
Out[32]: 
obj_id        objID        modelMag_g modelMag_i       obj_id1        type 
 str5         int64         float64    float64          int64         str6 
------ ------------------- ---------- ---------- ------------------- ------
 obj_0 1237654382516699210   17.57231   18.13992 1237654382516699210 GALAXY