2016-04-26 55 views
1

我想创建一个脚本来从一组对象中下载幅度。例如的那些,其在这里显示:如何使用Astropy/Astroquery获取斯隆幅度

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']) 

但是,从下载的图像中找不到大小。我怎样才能找到我的物体的大小?

回答

2

从事物的声音,而不是你想要的对象的属性,SDT管道吐出来,坐标交叉匹配。

要做到这一点看看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 
+0

非常感谢你正是我需要的。如果我可以利用你的时间......你知道在SDSS中是否可以用astropy/astroquery完成这样的事情? http://www.astroml.org/sklearn_tutorial/auto_examples/plot_sdss_images.html – Delosari