python
from astropy import coordinates as coord
from astropy.time import Time
from astropy.units import Unit
ra = coord.Angle('10h42m30s')
dec = coord.Angle('-30d00m00s')
c = coord.SkyCoord(ra, dec, frame='icrs')
galactic = c.transform_to(coord.Galactic)
print(galactic.l, galactic.b)
python
from astropy.io import fits
import matplotlib.pyplot as plt
hdulist = fits.open('image.fits')
image_data = hdulist[0].data
plt.imshow(image_data, cmap='gray')
plt.show()
python
from astropy.modeling.blackbody import blackbody_lambda
import numpy as np
flux = blackbody_lambda(wavelength, temperature)
plt.plot(wavelength, flux)
plt.show()
