Changeset 618

Show
Ignore:
Timestamp:
09/16/08 09:28:24 (2 months ago)
Author:
laidler
Message:

#149: Removed tapering from spectrum.py:.resample() methods

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved
  • trunk/pysynphot/lib/spectrum.py

    r617 r618  
    312312        hdu.writeto(filename) 
    313313 
    314             
    315  
    316314                                                  
    317315    def integrate(self,fluxunits='photlam'): 
     316        #Extract the flux in the desired units 
    318317        u=self.fluxunits 
    319318        self.convert(fluxunits) 
    320  
    321319        wave,flux=self.getArrays() 
    322320        self.convert(u) 
    323              
     321        #then do the integration     
    324322        return self.trapezoidIntegration(wave,flux) 
    325323     
     
    547545            newasc = False 
    548546 
    549         ## First need to pad the ends of the spectrum with zeros 
    550         tapered = self.taper() 
    551  
    552547        ## Use numpy interpolation function 
    553         if tapered._wavetable[0]<tapered._wavetable[-1]: 
     548        if self._wavetable[0]<self._wavetable[-1]: 
    554549            oldasc = True 
    555             ans = N.interp(newwave,tapered._wavetable, 
    556                            tapered._fluxtable) 
     550            ans = N.interp(newwave,self._wavetable, 
     551                           self._fluxtable) 
    557552        else: 
    558553            oldasc = False 
    559             rev = N.interp(newwave,tapered._wavetable[::-1], 
    560                            tapered._fluxtable[::-1]) 
     554            rev = N.interp(newwave,self._wavetable[::-1], 
     555                           self._fluxtable[::-1]) 
    561556            ans = rev[::-1] 
    562557 
     
    10081003            newasc = False 
    10091004 
    1010         #First need to pad the ends with zeros(?) 
    1011         tapered=self.taper() 
    10121005        ## Use numpy interpolation function 
    1013         if tapered._wavetable[0]<tapered._wavetable[-1]: 
     1006        if self._wavetable[0]<self._wavetable[-1]: 
    10141007            oldasc = True 
    1015             ans = N.interp(newwave,tapered._wavetable, 
    1016                            tapered._throughputtable) 
     1008            ans = N.interp(newwave,self._wavetable, 
     1009                           self._throughputtable) 
    10171010        else: 
    10181011            oldasc = False 
    1019             rev = N.interp(newwave,tapered._wavetable[::-1], 
    1020                            tapered._throughputtable[::-1]) 
     1012            rev = N.interp(newwave,self._wavetable[::-1], 
     1013                           self._throughputtable[::-1]) 
    10211014            ans = rev[::-1] 
    10221015