Bin_info_ = plate_[5].data
Flux_ = plate_[0].data
primhdu_ = plate_[0]
PLATEIDs.append(primhdu_.header['PLATEID'])
ORMASK.append(plate_[3].data)
ANDMASK.append(plate_[2].data)
INVAR.append(plate_[1].data)
log_wavst.append(primhdu_.header['COEFF0'])
MJDs.append(primhdu_.header['MJD'])
BinInfos.append(Bin_info_)
Flux.append(Flux_)
list = fits.open('Superset_DR12Q.fits', memmap=True) #opening file
supers = list[1].data # storing BINTABLE extension data
print("Restoring Data")
Full_Data = storing(PLATEIDs, supers)
X, Y, Train_z, Train_mag, And, In, wavst, ID = MLAData(Full_Data, BinInfos,
Flux, log_wavst,
ANDMASK, INVAR)
wav_ratio = 10**0.0001
plate_no = 0
all_wav = []
while plate_no < len(wavst):
append_count = 0
cent_wav = 10**wavst[plate_no]
wavelengths = []
wavelengths.append(cent_wav)
current_wav = cent_wav
while append_count < (4600 - 1):
current_wav = current_wav * wav_ratio
Flux_ = plate_[0].data
primhdu_ = plate_[0]
PLATEIDs.append(primhdu_.header['PLATEID'])
ORMASK.append(plate_[3].data)
ANDMASK.append(plate_[2].data)
INVAR.append(plate_[1].data)
log_wavst.append(primhdu_.header['COEFF0'])
MJDs.append(primhdu_.header['MJD'])
BinInfos.append(Bin_info_)
Flux.append(Flux_)
list = fits.open(TrainingDir + slash + 'Superset_DR12Q.fits',
memmap=True) #opening file
print(len(log_wavst))
supers = list[1].data # storing BINTABLE extension data
Full_Data = storing(PLATEIDs, supers)
X, Y, Train_z, Train_mag, wavst = MLAData(Full_Data, BinInfos, Flux, log_wavst,
ANDMASK, INVAR)
i = 0
d = open(TrainingDir + slash + "Full.txt", 'w')
while i < len(PLATEIDs):
C_Plate = PLATEIDs[i]
a1 = X[i]
print(len(a1))
a2 = np.array(Y[i])
print(len(a2))
if len(a1) == 0 & len(a2) == 0:
i = i + 1
else:
Bin_info_ = plate_[5].data
Flux_ = plate_[0].data
primhdu_ = plate_[0]
PLATEIDs.append(primhdu_.header['PLATEID'])
ORMASK.append(plate_[3].data)
ANDMASK.append(plate_[2].data)
INVAR.append(plate_[1].data)
log_wavst.append(primhdu_.header['COEFF0'])
MJDs.append(primhdu_.header['MJD'])
BinInfos.append(Bin_info_)
Flux.append(Flux_)
list = fits.open('Superset_DR12Q.fits', memmap=True) #opening file
supers = list[1].data # storing BINTABLE extension data
print("Restoring Data")
Full_Data = storing(PLATEIDs, supers)
X, Y, Train_z, Train_mag, And, In, wavst, ID, MJ, MatchedPlates = MLAData(
Full_Data, BinInfos, Flux, log_wavst, ANDMASK, INVAR, MJDs, PLATEIDs)
plate_no = 0
all_Testwav = []
wav_ratio = 10**0.0001
while plate_no < len(X):
append_count = 0
cent_wav = 10**wavst[plate_no]
wavelengths = []
wavelengths.append(cent_wav)
current_wav = cent_wav
while append_count < (4600 - 1):
current_wav = current_wav * wav_ratio
wavelengths.append(current_wav)