In [1]:
%matplotlib inline
import numpy as np
from astropy.io import fits as pf
import matplotlib.pyplot as plt
In [3]:
#open the candidate file
dfile = 'QSO_Candidates_Final_Clustering_Set.fits'
data = pf.open(dfile)[1].data
In [4]:
#Separate the 7 colors into arrays
u = data['ug'] + data ['gr'] + data['ri'] + data['imag']
g = data ['gr'] + data['ri'] + data['imag']
r = data['ri'] + data['imag']
i = data['imag']
z = data['imag'] - data['iz']
s1= data['imag'] - (data['zs1'] + data['iz'])
s2= data['imag'] - (data['s1s2'] + data['zs1'] + data['iz'])
In [25]:
u = u[:, np.newaxis]
g = g[:, np.newaxis]
r = r[:, np.newaxis]
i = i[:, np.newaxis]
z = z[:, np.newaxis]
s1= s1[:, np.newaxis]
s2= s2[:, np.newaxis]
In [34]:
colarray = np.hstack([u,g,r,i,z,s1,s2])
In [39]:
bandpass = [3543, 4770, 6231, 7625, 9134, 36000, 45000]
In [ ]: