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)