m42 fits文件处理

数据来自：https://dr12.sdss.org/fields/name?name=m42

import numpy as np
import matplotlib.pyplot as plt
from astropy.io import fits
m42=fits.open('fits_file/frame-g-006073-4-0063.fits')
m42.info()
m42[0].header
data=m42[0].data
plt.imshow(data)
plt.show()
np.max(data)
histogram=plt.hist(data.flat,bins=210)
#bins 条形数，color='r,g,y,c',density是否以密度的形式显示，
#histtype线条的类型"bar""barstacked",align对齐方式"left""mid""right"
#orientation "horizontal","vertical"
#log 是否以科学计数法
from matplotlib.colors import LogNorm
plt.imshow(data,norm=LogNorm(),cmap='gray')
plt.show()

m42u=fits.open('fits_file/frame-u-006073-4-0063.fits')
m42i=fits.open('fits_file/frame-i-006073-4-0063.fits')
m42g=fits.open('fits_file/frame-g-006073-4-0063.fits')

from astropy.visualization import make_lupton_rgb
g=m42g[0].data
i=m42i[0].data
u=m42u[0].data
rgb_default=make_lupton_rgb(g,i,u,stretch=0.5,Q=5)
plt.imshow(rgb_default,origin='lower')
plt.savefig("rgbm422.jpg")