Как использовать Astropy/Astroquery для получения значений Слоана

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 метод. Это позволяет вам указать, какой из столбцов фотометрии вы хотите.

Поскольку я не знаю науки, которой вы занимаетесь, вам придется решить, что лучше всего подходит для ваших целей. Посмотрите здесь на различные меры величины, которые доступны.

В качестве примера, если бы я хотел, чтобы значения 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