Changeset 1192 for SHX


Ignore:
Timestamp:
23.04.2016 12:00:53 (3 years ago)
Author:
klaus
Message:

command for plotting spectra; campat change

Location:
SHX/trunk/SeismicHandler/commands
Files:
1 added
1 edited

Legend:

Unmodified
Added
Removed
  • SHX/trunk/SeismicHandler/commands/psd.py

    r1191 r1192  
    1010from SeismicHandler.modules.screen import Screen 
    1111from obspy.signal import spectral_estimation 
     12from matplotlib import mlab 
    1213import pylab 
    1314 
     
    9394            else: 
    9495                loverlap = 0 
    95             pxx = spectral_estimation.psd(data, detrend=pylab.detrend_mean, 
    96                 NFFT=lnfft, noverlap=loverlap )[0] 
     96            #pxx = spectral_estimation.psd(data, detrend=pylab.detrend_mean, 
     97            #    NFFT=lnfft, noverlap=loverlap )[0] 
     98            pxx = mlab.psd( data, detrend=pylab.detrend_mean, 
     99                NFFT=lnfft, noverlap=loverlap, sides="onesided", 
     100                scale_by_freq=True )[0] 
    97101            pt = t.copy() 
    98102            pt.data = pxx[1:] 
Note: See TracChangeset for help on using the changeset viewer.