In [1]:
%autosave 60
%matplotlib inline
In [2]:
import matplotlib.pyplot as plt
import numpy as np
In [3]:
import astropy
In [4]:
import astropy.units as u
import astropy.constants as c
In [5]:
c.M_sun
Out[5]:
In [11]:
SOLLAR_MASS = c.M_sun.to_value('g')
type(SOLLAR_MASS), SOLLAR_MASS
Out[11]:
In [10]:
type(c.M_sun.cgs)
Out[10]:
In [14]:
c.M_sun.cgs.value, c.M_sun.cgs.unitt
Out[14]:
In [17]:
c.G.cgs.value
Out[17]:
In [20]:
x = 1 * u.km
x.to('au')
Out[20]:
In [21]:
x_val = x.to_value('pc')
type(x_val), x_val
Out[21]:
In [27]:
x = 10 * u.km
t = 1 * u.hour
v = x / t
v = v.to(u.cm / u.s)
v + 1 * u.mau / (10 * u.year) # не используйте mau
Out[27]:
In [30]:
(1 * u.kau).to(u.km)
Out[30]:
In [31]:
(10 * u.mpc).to(u.pc)
Out[31]:
In [32]:
10 * u.km + 7 * u.kg
In [33]:
v = 2 * np.pi * (1.0 * u.au) / (1.0 * u.year)
v.to(u.km / u.s)
Out[33]:
In [35]:
v = np.sqrt(c.G * c.M_sun / c.au)
v.to(u.km / u.s)
Out[35]:
In [39]:
x = (10 * u.kg) / (0.1 * u.M_jup)
print(x.to(u.dimensionless_unscaled))
np.exp(x)
Out[39]:
In [41]:
a = (10 * u.g) ** (1 / 3) # не делайте так
a
Out[41]:
In [45]:
(10 * u.g) ** (2 / 7 + 1 / 15) # так тем более
Out[45]:
In [47]:
from fractions import Fraction as Fr
In [48]:
Fr(2, 7) + Fr(1, 15)
Out[48]:
In [49]:
(10 * u.g) ** (Fr(2, 7) + Fr(1, 15))
Out[49]:
In [52]:
np.cos(100 * u.arcmin)
Out[52]:
In [57]:
delta_m = 1 * u.mag
mAB = 10 * u.ABmag
display(mAB.physical.to(u.erg / u.s / u.cm**2 / u.Hz))
m = mAB + delta_m
display(m.physical.to(u.erg / u.s / u.cm**2 / u.Hz))
In [61]:
m1 = u.Magnitude(10 * u.count / u.s)
m2 = u.Magnitude(1e4 * u.count / u.hour)
dm = (m1 - m2).decompose()
dm, dm.unit, dm.value
Out[61]:
In [63]:
a = np.arange(10) * u.cm
a
Out[63]:
In [62]:
from astropy import coordinates as coord
In [64]:
coord.SkyCoord('01h23m45.6s -05d23m17s')
Out[64]:
In [67]:
coord.SkyCoord(ra=10 * u.deg + 30 * u.arcmin, dec=-10 * u.deg - 30 * u.arcmin)
Out[67]:
In [69]:
coord.SkyCoord(ra=[1, 2, 3], dec=[5, 4, 3], unit=u.deg)
Out[69]:
In [75]:
coord.SkyCoord(ra=[1, 2, 3], dec=[5, 4, 3], unit=['h', 'deg'])
Out[75]:
In [77]:
coord.SkyCoord(l=5*u.deg, b=10*u.deg, frame='galactic')
Out[77]:
In [80]:
c = coord.SkyCoord('01h23m45.6s -05d23m17s')
c.galactic, c.ra, c.dec
Out[80]:
In [81]:
from astropy.time import Time
In [82]:
t = Time.now() - np.arange(5) * u.hour
t
Out[82]:
In [83]:
loc = coord.EarthLocation.of_site('subaru')
In [85]:
loc.height
Out[85]:
In [86]:
loc.geodetic
Out[86]:
In [88]:
alt_az_frame = coord.AltAz(obstime=t, location=loc)
In [89]:
alt_az_frame
Out[89]:
In [100]:
from astropy import visualization
display(c)
c.transform_to(alt_az_frame)
dt = np.arange(24) * u.hour
t = Time.now() + dt
moon = coord.get_moon(t)
display(moon)
alt_az_frame = coord.AltAz(obstime=t, location=loc)
moon_subaru = moon.transform_to(alt_az_frame)
with visualization.quantity_support():
plt.plot(dt, moon_subaru.alt.to(u.deg), 'x')
In [101]:
coord.SkyCoord.from_name('M31')
Out[101]:
In [103]:
from astroquery.simbad import Simbad
from astroquery.vizier import Vizier
In [104]:
m67 = Simbad.query_object('M67')
m67
Out[104]:
In [108]:
table = Simbad.query_region(coord.SkyCoord('08h51m18s +11d48m0s'),
radius=1 * u.arcmin)
table
Out[108]:
In [111]:
m67 = Simbad.query_object('M67')
# coord.SkyCoord.guess_from_table(m67)
center = coord.SkyCoord(ra=m67['RA'], dec=m67['DEC'],
unit=(u.hourangle, u.deg))
center
Out[111]:
In [116]:
vizier = Vizier(
row_limit=10000,
)
tmass_table = vizier.query_region(
center,
radius=5 * u.arcmin,
catalog='2MASS',
)[0]
display(tmass_table)
tmass = coord.SkyCoord(ra=tmass_table['RAJ2000'],
dec=tmass_table['DEJ2000'],
unit='deg')
sdss_table = vizier.query_region(
center,
radius=5 * u.arcmin,
catalog='SDSS7',
)[0]
display(sdss_table)
sdss = coord.SkyCoord(ra=sdss_table['RA_ICRS'],
dec=sdss_table['DE_ICRS'],
unit='deg')
In [137]:
tmass_close = tmass[center.separation(tmass) < 1 * u.arcmin]
plt.figure(figsize=(5,5))
with visualization.quantity_support():
plt.scatter(center.ra, center.dec, color='k', marker='*', s=96)
plt.scatter(tmass_close.ra, tmass_close.dec, marker='x',
color='red')
plt.scatter(sdss.ra, sdss.dec, marker='.', color='blue', s=4)
In [148]:
idx_tmass, idx_sdss, sep, dist = sdss.search_around_sky(
tmass_close,
1*u.arcsec,
)
with visualization.quantity_support():
for a, b in zip(tmass[idx_tmass], sdss[idx_sdss]):
plt.plot(u.Quantity([a.ra, b.ra]), u.Quantity([a.dec, b.dec]),
'-', color='blue')
In [151]:
idx_tmass, sep, dist = sdss.match_to_catalog_sky(tmass)
with visualization.quantity_support():
for a, b in zip(tmass[idx_tmass], sdss):
plt.plot(u.Quantity([a.ra, b.ra]), u.Quantity([a.dec, b.dec]),
'-', color='blue')
In [ ]: