|
21 | 21 | import numpy as np
|
22 | 22 |
|
23 | 23 |
|
24 |
| -from doppler.lsf import GaussianLsf, GaussHermiteLsf |
25 |
| -# LSF class dictionary |
26 |
| -lsfclass = {'gaussian': GaussianLsf, 'gauss-hermite': GaussHermiteLsf} |
27 | 24 | from apogee.apred import wave
|
28 | 25 | from apogee.apred import sincint
|
29 | 26 |
|
@@ -66,26 +63,6 @@ def mask(self,bdval) :
|
66 | 63 | """
|
67 | 64 | self.mask=(np.bitwise_and(self.bitmask,bdval)!=0) | (np.isfinite(self.flux)==False)
|
68 | 65 |
|
69 |
| - def Spec1D(self,res=22500) : |
70 |
| - """ Convert to a Nidever spec1D object, with axes flipped |
71 |
| - """ |
72 |
| - self.flux = self.flux.T |
73 |
| - if self.err is not None : self.err = self.err.T |
74 |
| - if self.bitmask is not None : self.bitmask = self.bitmask.T |
75 |
| - if self.mask is not None : self.mask = self.mask.T |
76 |
| - if self.sky is not None : self.sky = self.sky.T |
77 |
| - if self.skyerr is not None : self.skyerr = self.skyerr.T |
78 |
| - if self.telluric is not None : self.telluric = self.telluric.T |
79 |
| - if self.telerr is not None : self.telerr = self.telerr.T |
80 |
| - bd = np.where( (np.isfinite(self.flux)==False) | (self.err <= 0.0) )[0] |
81 |
| - if len(bd)>0: |
82 |
| - self.flux[bd] = 0.0 |
83 |
| - self.err[bd] = 1e30 |
84 |
| - self.mask[bd] = True |
85 |
| - self.lsftype = 'Gaussian' |
86 |
| - self.lsfsig = self.wave/res/2.354 |
87 |
| - self.lsf = lsfclass[self.lsftype.lower()](wave=self.wave,xtype='Wave',lsftype=self.lsftype,sigma=self.lsfsig) |
88 |
| - |
89 | 66 | def interp(self,new,nres) :
|
90 | 67 | """ Interpolate to new wavelengths
|
91 | 68 | """
|
|
0 commit comments