python - 如何使用 Astropy/Astroquery 获取 Sloan 量级

Question

我想创建一个脚本来从一组对象下载震级。例如出现在这里的那些:

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

作为输入，我有对象坐标(这是我从 astroquery get_spec 方法获得的，因为我知道对象的 mjd、板和纤维)。我正在尝试使用 astroquery 网站上的示例:

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方法。这允许您指定所需的光度列。

由于我不了解您正在做的科学，您必须找出最适合您的目标的方法。看看here对于可用的各种量级。

例如，如果我想要 g 和 r 波段 model magnitudes对于您的对象，我会这样做:

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