In [58]:
%autosave 20
%matplotlib inline
import matplotlib.pyplot as plt
import numpy as np
import astropy.units as u
from astropy.coordinates import SkyCoord
from astroquery.vizier import Vizier
from astroquery.simbad import Simbad
from astropy import table
Autosaving every 20 seconds
In [9]:
m67 = Simbad.query_object('M67')
display(m67)
# SkyCoord.guess_from_table(m67['RA', 'DEC'])
center = SkyCoord(ra=m67['RA'], dec=m67['DEC'],
unit=(u.hourangle, u.deg))
center
Table masked=True length=1
MAIN_ID RA DEC RA_PREC DEC_PREC COO_ERR_MAJA COO_ERR_MINA COO_ERR_ANGLE COO_QUAL COO_WAVELENGTH COO_BIBCODE "h:m:s" "d:m:s" mas mas deg object str13 str13 int16 int16 float32 float32 int16 str1 str1 object NGC 2682 08 51 18 +11 48.0 4 4 -- -- 0 E 2005ApJ...619..824X
Out[9]:
<SkyCoord (ICRS): (ra, dec) in deg
[(132.825, 11.8)]>
In [19]:
vizier = Vizier(
column_filters={},
row_limit=10000
)
twomass_table = vizier.query_region(
center,
radius=5*u.arcmin,
catalog='2MASS'
)[0]
display(twomass_table)
# twomass = SkyCoord.guess_from_table(twomass_table)
twomass = SkyCoord(
ra=twomass_table['RAJ2000'],
dec=twomass_table['DEJ2000'],
unit=(u.deg, u.deg),
)
sdss_table = vizier.query_region(
center,
radius=5*u.arcmin,
catalog='SDSS7'
)[0]
display(sdss_table)
# sdss = SkyCoord.guess_from_table(sdss_table)
sdss = SkyCoord(
ra=sdss_table['RA_ICRS'],
dec=sdss_table['DE_ICRS'],
unit=(u.deg, u.deg),
)
Table masked=True length=271
_q RAJ2000 DEJ2000 _2MASS Jmag e_Jmag Hmag e_Hmag Kmag e_Kmag Qflg Rflg Bflg Cflg Xflg Aflg deg deg mag mag mag mag mag mag int32 float64 float64 bytes17 float32 float32 float32 float32 float32 float32 bytes3 bytes3 bytes3 bytes3 uint8 uint8 1 132.810285 11.743501 08511446+1144366 13.859 0.036 13.479 0.041 13.334 0.044 AAA 222 222 ccc 0 0 1 132.810279 11.744696 08511446+1144409 14.176 0.028 13.650 0.023 13.598 0.036 AAA 222 222 ccc 0 0 1 132.793054 11.730462 08511033+1143496 16.286 0.104 15.553 0.121 15.463 0.172 ABC 222 111 000 0 0 1 132.806171 11.740373 08511348+1144253 16.506 0.123 16.310 -- 15.894 -- BUU 200 100 000 0 0 1 132.786225 11.733243 08510869+1143596 16.539 0.126 16.255 0.229 15.621 -- BDU 220 110 000 0 0 1 132.817854 11.742462 08511628+1144328 12.308 0.021 12.021 0.020 11.954 0.020 AAA 222 111 000 0 0 1 132.839180 11.733575 08512140+1144008 14.821 0.037 14.325 0.049 14.082 0.052 AAA 222 111 cc0 0 0 1 132.840989 11.721581 08512183+1143176 13.274 0.023 13.003 0.023 12.952 0.025 AAA 222 111 000 0 0 1 132.848917 11.736172 08512374+1144102 16.669 0.142 16.011 0.192 15.790 0.229 BCD 222 111 000 0 0 ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... 1 132.892817 11.828841 08513427+1149438 11.804 0.022 11.569 0.022 11.522 0.020 AAA 222 111 000 0 0 1 132.892725 11.848458 08513425+1150544 12.585 0.023 12.341 0.023 12.257 0.023 AAA 222 111 000 0 0 1 132.886145 11.841920 08513267+1150309 14.164 0.030 13.679 0.032 13.614 0.039 AAA 222 111 000 0 0 1 132.883952 11.834337 08513214+1150036 11.514 0.022 11.226 0.020 11.198 0.020 AAA 222 111 000 0 0 1 132.885801 11.814464 08513259+1148520 10.645 0.022 10.541 0.020 10.526 0.018 AAA 222 111 000 0 0 1 132.885731 11.844612 08513257+1150406 11.730 0.022 11.657 0.020 11.617 0.020 AAA 222 111 000 0 0 1 132.802904 11.878436 08511269+1152423 8.650 0.018 8.122 0.018 7.976 0.018 AAA 111 111 000 0 0 1 132.833993 11.879987 08512015+1152479 12.772 0.022 12.545 0.020 12.457 0.021 AAA 222 111 000 0 0 1 132.856334 11.877449 08512552+1152388 12.150 0.022 11.897 0.020 11.810 0.021 AAA 222 111 000 0 0 1 132.840646 11.877172 08512175+1152378 10.105 0.023 9.816 0.022 9.760 0.018 AAA 222 111 000 0 0
Table masked=True length=1401
_q mode cl SDSS m_SDSS zsp umag e_umag gmag e_gmag rmag e_rmag imag e_imag zmag e_zmag RA_ICRS DE_ICRS ObsDate Q mag mag mag mag mag mag mag mag mag mag deg deg yr int32 uint8 uint8 bytes19 bytes1 float64 float32 float32 float32 float32 float32 float32 float32 float32 float32 float32 float64 float64 float64 uint8 1 2 3 J085116.77+114829.2 -- 20.841 0.431 22.441 0.804 24.681 6.832 24.270 8.114 18.390 0.166 132.819886 11.808115 2006.0844 3 1 1 6 J085116.78+115038.8 -- 14.967 0.006 15.049 0.013 14.243 0.009 13.369 0.009 13.202 0.005 132.819920 11.844135 2005.9314 3 1 2 6 J085116.78+115038.9 -- 15.047 0.006 14.535 0.009 14.196 0.009 15.187 0.011 13.249 0.004 132.819926 11.844151 2006.0161 3 1 1 6 J085116.81+114541.6 -- 15.873 0.005 14.439 0.003 13.921 0.003 13.735 0.003 13.726 0.004 132.820053 11.761561 2005.9314 3 1 1 3 J085116.81+114807.4 -- 17.625 0.074 19.837 0.141 16.103 0.013 14.947 0.016 16.534 0.113 132.820063 11.802067 2005.9314 3 1 2 6 J085116.81+114541.6 * -- 16.432 0.008 14.911 0.005 15.822 0.014 14.156 0.005 14.214 0.007 132.820069 11.761566 2006.0844 3 1 1 6 J085116.85+114938.6 -- 19.266 0.030 16.806 0.004 15.733 0.004 15.367 0.004 15.167 0.005 132.820234 11.827405 2005.9314 3 1 2 6 J085116.85+114938.6 * -- 19.193 0.027 16.791 0.004 15.736 0.003 15.356 0.003 15.152 0.005 132.820239 11.827412 2006.0161 3 1 1 6 J085116.97+115009.3 -- 14.885 0.004 14.539 0.011 14.919 0.011 15.225 0.014 13.245 0.005 132.820726 11.835924 2005.9314 3 ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... 1 2 6 J085116.32+114711.3 -- 24.398 0.716 21.991 0.069 20.450 0.029 19.081 0.015 18.308 0.020 132.818016 11.786476 2006.0844 3 1 2 6 J085116.43+114834.3 -- 25.592 0.811 21.966 0.078 21.777 0.104 20.221 0.042 20.405 0.131 132.818471 11.809527 2006.0161 3 1 1 6 J085116.46+114607.3 -- 23.583 0.844 20.807 0.033 19.231 0.013 18.222 0.009 17.652 0.017 132.818608 11.768719 2005.9314 3 1 1 6 J085116.47+114900.1 -- 21.381 0.139 18.740 0.009 17.361 0.005 16.798 0.005 16.471 0.008 132.818628 11.816693 2005.9314 3 1 2 6 J085116.47+114900.1 * -- 21.117 0.107 18.719 0.008 17.369 0.005 16.795 0.005 16.503 0.008 132.818630 11.816701 2006.0161 3 1 2 6 J085116.50+114723.9 -- 22.393 0.218 22.069 0.077 22.271 0.137 22.261 0.189 23.057 0.419 132.818769 11.789983 2006.0844 3 1 2 3 J085116.61+114529.3 -- 24.631 9.999 17.000 0.020 16.128 0.015 16.188 0.019 14.333 0.017 132.819241 11.758138 2006.0844 3 1 1 3 J085116.65+114528.9 -- 24.358 3.652 16.021 0.015 16.623 0.024 16.473 0.024 22.823 9.999 132.819412 11.758052 2005.9314 3 1 1 6 J085116.68+114529.3 -- 14.846 0.005 14.886 0.011 14.748 0.011 14.636 0.012 12.490 0.004 132.819537 11.758144 2005.9314 3 1 2 6 J085116.69+114529.3 -- 14.571 0.006 15.403 0.012 15.798 0.014 14.651 0.013 13.414 0.014 132.819564 11.758164 2006.0844 3
In [36]:
plt.figure(figsize=(6,6))
plt.plot(twomass.ra.to_value(u.deg),
twomass.dec.to_value(u.deg),
'x', color='blue')
plt.plot(sdss.ra.to_value(u.deg),
sdss.dec.to_value(u.deg),
'.', color='green')
plt.scatter(center.ra.to_value(u.deg),
center.dec.to_value(u.deg),
marker='*', color='red')
Out[36]:
<matplotlib.collections.PathCollection at 0x7f22e4809940>
In [39]:
idx_close_to_center = center.separation(twomass) < 1 * u.arcmin
twomass_table[idx_close_to_center]
Out[39]:
Table masked=True length=13
_q RAJ2000 DEJ2000 _2MASS Jmag e_Jmag Hmag e_Hmag Kmag e_Kmag Qflg Rflg Bflg Cflg Xflg Aflg deg deg mag mag mag mag mag mag int32 float64 float64 bytes17 float32 float32 float32 float32 float32 float32 bytes3 bytes3 bytes3 bytes3 uint8 uint8 1 132.811505 11.789971 08511476+1147238 11.747 0.021 11.495 0.022 11.430 0.018 AAA 222 111 000 0 0 1 132.815406 11.804343 08511569+1148156 13.244 0.024 12.965 0.027 12.955 0.026 AAA 222 111 ddd 0 0 1 132.827869 11.784062 08511868+1147026 11.017 0.021 10.619 0.020 10.529 0.018 AAA 222 111 000 0 0 1 132.829189 11.789187 08511900+1147210 15.819 0.074 15.092 0.089 15.042 0.116 AAB 222 111 sss 0 0 1 132.829927 11.798506 08511918+1147546 12.566 0.021 12.247 0.020 12.182 0.020 AAA 222 111 000 0 0 1 132.821257 11.804460 08511710+1148160 8.140 0.027 7.526 0.018 7.385 0.021 AAA 111 111 000 0 0 1 132.830069 11.807475 08511921+1148269 15.151 0.042 14.600 0.054 14.312 0.060 AAA 222 111 000 0 0 1 132.823341 11.804540 08511760+1148163 14.616 0.023 11.767 -- 11.455 -- AUU 200 100 c00 0 0 1 132.817966 11.786548 08511631+1147115 16.762 0.164 15.665 -- 15.393 -- CUU 266 200 c00 0 0 1 132.836275 11.808737 08512070+1148314 15.468 0.048 14.843 0.066 14.562 0.074 AAA 222 111 000 0 0 1 132.814363 11.792049 08511544+1147313 11.708 0.021 11.475 0.020 11.430 0.018 AAA 222 111 000 0 0 1 132.832875 11.810451 08511989+1148376 15.084 0.043 14.508 0.048 14.445 0.067 AAA 222 111 000 0 0 1 132.817027 11.787326 08511608+1147143 15.947 0.086 15.164 -- 14.736 -- AUU 266 200 c00 0 0
In [43]:
idx_twomass, idx_sdss, d2d, d3d = sdss.search_around_sky(
twomass, 2*u.arcsec
)
twomass_table[idx_twomass]
Out[43]:
Table masked=True length=556
_q RAJ2000 DEJ2000 _2MASS Jmag e_Jmag Hmag e_Hmag Kmag e_Kmag Qflg Rflg Bflg Cflg Xflg Aflg deg deg mag mag mag mag mag mag int32 float64 float64 bytes17 float32 float32 float32 float32 float32 float32 bytes3 bytes3 bytes3 bytes3 uint8 uint8 1 132.810285 11.743501 08511446+1144366 13.859 0.036 13.479 0.041 13.334 0.044 AAA 222 222 ccc 0 0 1 132.810285 11.743501 08511446+1144366 13.859 0.036 13.479 0.041 13.334 0.044 AAA 222 222 ccc 0 0 1 132.810279 11.744696 08511446+1144409 14.176 0.028 13.650 0.023 13.598 0.036 AAA 222 222 ccc 0 0 1 132.810279 11.744696 08511446+1144409 14.176 0.028 13.650 0.023 13.598 0.036 AAA 222 222 ccc 0 0 1 132.793054 11.730462 08511033+1143496 16.286 0.104 15.553 0.121 15.463 0.172 ABC 222 111 000 0 0 1 132.793054 11.730462 08511033+1143496 16.286 0.104 15.553 0.121 15.463 0.172 ABC 222 111 000 0 0 1 132.806171 11.740373 08511348+1144253 16.506 0.123 16.310 -- 15.894 -- BUU 200 100 000 0 0 1 132.806171 11.740373 08511348+1144253 16.506 0.123 16.310 -- 15.894 -- BUU 200 100 000 0 0 1 132.786225 11.733243 08510869+1143596 16.539 0.126 16.255 0.229 15.621 -- BDU 220 110 000 0 0 ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... 1 132.885731 11.844612 08513257+1150406 11.730 0.022 11.657 0.020 11.617 0.020 AAA 222 111 000 0 0 1 132.885731 11.844612 08513257+1150406 11.730 0.022 11.657 0.020 11.617 0.020 AAA 222 111 000 0 0 1 132.802904 11.878436 08511269+1152423 8.650 0.018 8.122 0.018 7.976 0.018 AAA 111 111 000 0 0 1 132.802904 11.878436 08511269+1152423 8.650 0.018 8.122 0.018 7.976 0.018 AAA 111 111 000 0 0 1 132.833993 11.879987 08512015+1152479 12.772 0.022 12.545 0.020 12.457 0.021 AAA 222 111 000 0 0 1 132.833993 11.879987 08512015+1152479 12.772 0.022 12.545 0.020 12.457 0.021 AAA 222 111 000 0 0 1 132.856334 11.877449 08512552+1152388 12.150 0.022 11.897 0.020 11.810 0.021 AAA 222 111 000 0 0 1 132.856334 11.877449 08512552+1152388 12.150 0.022 11.897 0.020 11.810 0.021 AAA 222 111 000 0 0 1 132.840646 11.877172 08512175+1152378 10.105 0.023 9.816 0.022 9.760 0.018 AAA 222 111 000 0 0 1 132.840646 11.877172 08512175+1152378 10.105 0.023 9.816 0.022 9.760 0.018 AAA 222 111 000 0 0
In [59]:
idx_sdss, d2d, d3d = twomass.match_to_catalog_sky(sdss)
t = table.hstack((twomass_table, sdss_table[idx_sdss]))
t['sep'] = d2d
t
/opt/conda/lib/python3.6/site-packages/astropy/utils/metadata.py:360: MergeConflictWarning: Cannot merge meta key 'ID' types <class 'str'> and <class 'str'>, choosing ID='II_294_sdss7'
MergeConflictWarning)
/opt/conda/lib/python3.6/site-packages/astropy/utils/metadata.py:360: MergeConflictWarning: Cannot merge meta key 'name' types <class 'str'> and <class 'str'>, choosing name='II/294/sdss7'
MergeConflictWarning)
/opt/conda/lib/python3.6/site-packages/astropy/utils/metadata.py:360: MergeConflictWarning: Cannot merge meta key 'description' types <class 'str'> and <class 'str'>, choosing description='* output of the SDSS photometric catalog'
MergeConflictWarning)
Out[59]:
Table masked=True length=271
_q_1 RAJ2000 DEJ2000 _2MASS Jmag e_Jmag Hmag e_Hmag Kmag e_Kmag Qflg Rflg Bflg Cflg Xflg Aflg _q_2 mode cl SDSS m_SDSS zsp umag e_umag gmag e_gmag rmag e_rmag imag e_imag zmag e_zmag RA_ICRS DE_ICRS ObsDate Q sep deg deg mag mag mag mag mag mag mag mag mag mag mag mag mag mag mag mag deg deg yr deg int32 float64 float64 bytes17 float32 float32 float32 float32 float32 float32 bytes3 bytes3 bytes3 bytes3 uint8 uint8 int32 uint8 uint8 bytes19 bytes1 float64 float32 float32 float32 float32 float32 float32 float32 float32 float32 float32 float64 float64 float64 uint8 float64 1 132.810285 11.743501 08511446+1144366 13.859 0.036 13.479 0.041 13.334 0.044 AAA 222 222 ccc 0 0 1 1 6 J085114.46+114436.5 -- 17.255 0.009 15.653 0.004 15.071 0.003 14.892 0.004 14.784 0.005 132.810253 11.743496 2006.0844 3 3.17266619626919e-05 1 132.810279 11.744696 08511446+1144409 14.176 0.028 13.650 0.023 13.598 0.036 AAA 222 222 ccc 0 0 1 1 6 J085114.45+114440.9 -- 17.981 0.013 16.210 0.004 15.399 0.003 15.203 0.004 15.341 0.005 132.810246 11.744697 2006.0844 3 3.232459429986805e-05 1 132.793054 11.730462 08511033+1143496 16.286 0.104 15.553 0.121 15.463 0.172 ABC 222 111 000 0 0 1 1 6 J085110.32+114349.5 -- 22.852 0.280 20.067 0.017 18.661 0.009 18.001 0.008 17.614 0.014 132.793028 11.730441 2006.0844 3 3.300088183610788e-05 1 132.806171 11.740373 08511348+1144253 16.506 0.123 16.310 -- 15.894 -- BUU 200 100 000 0 0 1 2 6 J085113.48+114425.6 * -- 25.300 1.040 21.810 0.071 20.392 0.030 18.975 0.013 18.226 0.025 132.806196 11.740456 2005.9314 3 8.65339410341914e-05 1 132.786225 11.733243 08510869+1143596 16.539 0.126 16.255 0.229 15.621 -- BDU 220 110 000 0 0 1 1 6 J085108.61+114359.3 -- 20.928 0.065 19.075 0.010 18.351 0.008 18.136 0.009 18.019 0.017 132.785914 11.733164 2006.0844 3 0.00031458271139627445 1 132.817854 11.742462 08511628+1144328 12.308 0.021 12.021 0.020 11.954 0.020 AAA 222 111 000 0 0 1 1 6 J085116.28+114432.8 -- 16.157 0.011 16.007 0.019 15.607 0.012 15.057 0.018 14.451 0.012 132.817833 11.742445 2006.0844 3 2.6678360157593366e-05 1 132.839180 11.733575 08512140+1144008 14.821 0.037 14.325 0.049 14.082 0.052 AAA 222 111 cc0 0 0 1 1 6 J085121.39+114400.8 -- 20.172 0.040 17.701 0.006 16.579 0.004 16.210 0.005 15.980 0.007 132.839139 11.733563 2006.0844 3 4.189846167827402e-05 1 132.840989 11.721581 08512183+1143176 13.274 0.023 13.003 0.023 12.952 0.025 AAA 222 111 000 0 0 1 1 6 J085121.83+114317.6 -- 16.307 0.007 15.082 0.005 15.596 0.012 14.390 0.005 14.292 0.006 132.840975 11.721581 2006.0844 3 1.3708049052631302e-05 1 132.848917 11.736172 08512374+1144102 16.669 0.142 16.011 0.192 15.790 0.229 BCD 222 111 000 0 0 1 2 6 J085123.71+114410.6 -- 25.270 1.038 21.171 0.042 19.661 0.018 18.613 0.011 17.967 0.021 132.848824 11.736295 2005.9314 3 0.00015303644818535168 ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... 1 132.892817 11.828841 08513427+1149438 11.804 0.022 11.569 0.022 11.522 0.020 AAA 222 111 000 0 0 1 2 6 J085134.27+114943.8 * -- 14.876 0.005 13.322 0.003 14.191 0.009 12.781 0.001 13.041 0.004 132.892807 11.828832 2006.0161 3 1.3296539552404174e-05 1 132.892725 11.848458 08513425+1150544 12.585 0.023 12.341 0.023 12.257 0.023 AAA 222 111 000 0 0 1 1 6 J085134.25+115054.4 -- 15.299 0.005 14.207 0.003 14.482 0.011 13.989 0.012 13.608 0.004 132.892711 11.848453 2005.9314 3 1.458550850287712e-05 1 132.886145 11.841920 08513267+1150309 14.164 0.030 13.679 0.032 13.614 0.039 AAA 222 111 000 0 0 1 1 6 J085132.67+115030.8 -- 18.088 0.016 16.176 0.003 15.460 0.004 15.241 0.004 15.108 0.005 132.886143 11.841908 2005.9314 3 1.2158599918403357e-05 1 132.883952 11.834337 08513214+1150036 11.514 0.022 11.226 0.020 11.198 0.020 AAA 222 111 000 0 0 1 2 6 J085132.14+115003.5 -- 14.740 0.005 14.954 0.010 14.723 0.009 14.594 0.010 13.274 0.005 132.883938 11.834324 2006.0161 3 1.8887997936612012e-05 1 132.885801 11.814464 08513259+1148520 10.645 0.022 10.541 0.020 10.526 0.018 AAA 222 111 000 0 0 1 1 6 J085132.59+114851.9 -- 14.417 0.007 14.000 0.008 11.589 0.001 13.202 0.008 13.250 0.012 132.885798 11.814436 2005.9314 3 2.815355613377517e-05 1 132.885731 11.844612 08513257+1150406 11.730 0.022 11.657 0.020 11.617 0.020 AAA 222 111 000 0 0 1 1 6 J085132.57+115040.5 -- 14.558 0.006 14.601 0.009 14.923 0.011 14.186 0.010 13.441 0.007 132.885714 11.844603 2005.9314 3 1.8916241033416936e-05 1 132.802904 11.878436 08511269+1152423 8.650 0.018 8.122 0.018 7.976 0.018 AAA 111 111 000 0 0 1 2 6 J085112.68+115242.3 -- 15.155 0.008 11.021 0.001 10.173 0.001 12.134 0.003 10.411 0.001 132.802859 11.878427 2006.0161 3 4.49466797568373e-05 1 132.833993 11.879987 08512015+1152479 12.772 0.022 12.545 0.020 12.457 0.021 AAA 222 111 000 0 0 1 2 6 J085120.15+115247.9 * -- 15.296 0.004 14.071 0.003 13.699 0.003 14.988 0.012 13.600 0.003 132.833976 11.879981 2006.0161 3 1.7684806515114607e-05 1 132.856334 11.877449 08512552+1152388 12.150 0.022 11.897 0.020 11.810 0.021 AAA 222 111 000 0 0 1 2 6 J085125.51+115238.7 -- 14.970 0.004 14.346 0.011 13.956 0.010 15.090 0.011 13.223 0.004 132.856327 11.877434 2006.0161 3 1.649012706041507e-05 1 132.840646 11.877172 08512175+1152378 10.105 0.023 9.816 0.022 9.760 0.018 AAA 222 111 000 0 0 1 1 6 J085121.74+115237.8 -- 14.588 0.007 14.476 0.008 10.468 0.001 10.394 0.001 11.480 0.002 132.840623 11.877171 2005.9314 3 2.252979828595256e-05
In [ ]:
Content source: hombit/scientific_python
Similar notebooks: