Name = dataheader.get('OBJECT') headerKWs['OBJECT'] = Name headerKWs["DATE-OBS"] = dataheader.get('DATE-OBS') headerKWs["INSTRUMENT"] = dataheader.get('INSTRUME') headerKWs["OBSERVER"] = dataheader.get('OBSERVER') headerKWs["EXPTIME"] = dataheader.get('EXPTIME') Score = Moog960.Score() Melody = Moog960.ObservedMelody(Score=Score, name=Name) for wl, I, sn, v in zip(wlsol, data, snr, var): header = pyfits.Header() header.set('NAME', Name) header.set("WLSTART", numpy.min(wl)) header.set("WLSTOP", numpy.max(wl)) Melody.addPhrase(Moog960.ObservedPhrase(observedData=[SpectralTools.Spectrum(wl=wl, I=I, dI=I/sn, header=header, spectrum_type="OBSERVED")], Melody=Melody)) ID = Melody.phrases[-1].ID parameters = {} parameters["PHRASE"] = ID parameters["MELODY"] = Melody.ID parameters["SCORE"] = Melody.Score.ID parameters["WLSTART"] = header.get('WLSTART') parameters["WLSTOP"] = header.get('WLSTOP') parameters["SELECTED"] = False if not(Melody.ID in Melody.Score.observed_labels.keys()): Melody.Score.observed_labels[Melody.ID] = {} Melody.Score.observed_labels[Melody.ID][ID] = [] Melody.Score.observed_labels[Melody.ID][ID].append(Moog960.Label(parameters, reference=Melody.phrases[-1].observedData[-1], Phrase=Melody.phrases[-1],
headerKWs["OBSERVER"] = dataheader.get('OBSERVER') headerKWs["EXPTIME"] = dataheader.get('EXPTIME') Score = Moog960.Score() Melody = Moog960.ObservedMelody(Score=Score, name=Name) for wl, I, sn, v in zip(wlsol, data, snr, var): header = pyfits.Header() header.set('NAME', Name) header.set("WLSTART", numpy.min(wl)) header.set("WLSTOP", numpy.max(wl)) Melody.addPhrase( Moog960.ObservedPhrase(observedData=[ SpectralTools.Spectrum(wl=wl, I=I, dI=I / sn, header=header, spectrum_type="OBSERVED") ], Melody=Melody)) ID = Melody.phrases[-1].ID parameters = {} parameters["PHRASE"] = ID parameters["MELODY"] = Melody.ID parameters["SCORE"] = Melody.Score.ID parameters["WLSTART"] = header.get('WLSTART') parameters["WLSTOP"] = header.get('WLSTOP') parameters["SELECTED"] = False if not (Melody.ID in Melody.Score.observed_labels.keys()): Melody.Score.observed_labels[Melody.ID] = {}