Changeset 1090


Ignore:
Timestamp:
09.01.2016 16:34:14 (4 years ago)
Author:
klaus
Message:

bugfixes; added file output on fk figures

Location:
SHX/trunk/SeismicHandler/commands
Files:
4 edited

Legend:

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

    r1077 r1090  
    1313from obspy.core.util.geodetics import gps2DistAzimuth 
    1414from obspy.taup import TauPyModel 
     15from obspy.sh.core import fromUTCDateTime 
    1516 
    1617DEG_TO_KM = 111.19 
     
    113114            fdsnadr = self.qualifiers["ADDR"] 
    114115         
     116        # assume that then the first (and only) origin is the preferred one 
     117        evpar['includeallorigins'] = False 
     118        evpar['includeallmagnitudes'] = False 
     119         
    115120        client = Client( fdsnadr ) 
    116121        evlist = client.get_events( **evpar ) 
     
    118123            if len(ev.origins) < 0: 
    119124                continue 
    120             # take first origin, dont' know how to get preferred 
     125            # take first origin, shouldn't be any others 
    121126            orig = ev.origins[0] 
    122127            magn = ev.magnitudes[0] 
     
    132137                    continue 
    133138                offtime += arrivals[0].time 
    134             timestr = offtime.strftime( "%d-%h-%Y_%T" ) 
     139            timestr = fromUTCDateTime( offtime ) 
    135140            print >> Settings.echo_channel, "%s %7.3f %8.3f %5.1f %3.1f" % ( 
    136141                timestr,orig.latitude,orig.longitude,kmdepth,magn.mag) 
  • SHX/trunk/SeismicHandler/commands/fetch.py

    r1063 r1090  
    111111            log_message("warning.fetch", msg % ", ".join(stations)) 
    112112        add_traces(final, gain=True) 
     113 
     114        # reread DB due to bug in Db interface 
     115        Stations().read( clear=True ) 
    113116 
    114117    def from_cache(self, station, channel, component, start, duration): 
  • SHX/trunk/SeismicHandler/commands/fetch_area.py

    r1063 r1090  
    161161            log_message("warning.fetch", msg % ", ".join(chanlist)) 
    162162        add_traces(final, gain=True) 
     163 
     164        # reread DB due to bug in Db interface 
     165        Stations().read( clear=True ) 
     166 
    163167 
    164168    def _getHostPortAndUser( self ): 
  • SHX/trunk/SeismicHandler/commands/fk.py

    r1081 r1090  
    133133        plottype = None 
    134134        timeplot_possible = False 
     135        slownessmap_plot = None 
     136        polar_plot = None 
     137        scatter_plot = None 
    135138        if len(self.parameters) > 3 and self.parameters[3]: 
    136139            frqmin = float( self.parameters[3] ) 
     
    154157        if self.qualifiers["SCATTER"]: 
    155158            plottype = "TimeScatter" 
     159            if self.qualifiers["SCATTER"] is not True: 
     160                scatter_plot = self.qualifiers["SCATTER"] 
    156161        if self.qualifiers["POLAR"]: 
    157162            plottype = "TimePolar" 
     163            if self.qualifiers["POLAR"] is not True: 
     164                polar_plot = self.qualifiers["POLAR"] 
    158165        if self.qualifiers["SLOWNESSMAP"]: 
    159166            plottype = "SlownessMap" 
     167            if self.qualifiers["SLOWNESSMAP"] is not True: 
     168                slownessmap_plot = self.qualifiers["SLOWNESSMAP"] 
    160169        if self.qualifiers["KMUNITS"]: 
    161170            cnvfactor = 1. 
     
    230239                title = "slo:%5.2f  az:%5.1f" % (result_slo,result_az) 
    231240            self.plotSlowmap( np.fliplr(np.rot90(global_relpowmap*cnvfactor)), 
    232                 slomax*cnvfactor, title ) 
     241                slomax*cnvfactor, title, outfile=slownessmap_plot ) 
    233242        elif plottype == 'TimeScatter': 
    234             self.plotScatter( fkout ) 
     243            self.plotScatter( fkout, outfile=scatter_plot ) 
    235244        elif plottype == 'TimePolar': 
    236             self.plotPolar( fkout, slomax*cnvfactor ) 
     245            self.plotPolar( fkout, slomax*cnvfactor, outfile=polar_plot ) 
    237246        elif plottype == None: 
    238247            pass 
     
    241250         
    242251        
    243     def plotScatter( self, fkout ): 
     252    def plotScatter( self, fkout, outfile=None ): 
    244253        if len(fkout.T[0]) < 10: 
    245254            print "trace too short to plot" 
     
    262271        fig.autofmt_xdate() 
    263272        fig.subplots_adjust(left=0.15, top=0.95, right=0.95, bottom=0.2, hspace=0) 
    264         plt.show() 
     273        if outfile == None: 
     274            plt.show() 
     275        else: 
     276            plt.savefig( outfile.lower() ) 
    265277     
    266     def plotPolar( self, fkout, slomax ): 
     278    def plotPolar( self, fkout, slomax, outfile=None ): 
    267279        if len(fkout.T[0]) < 10: 
    268280            print "trace too short to plot" 
     
    313325        ColorbarBase(cax, cmap=cmap, 
    314326            norm=Normalize(vmin=hist.min(), vmax=hist.max())) 
    315  
    316         plt.show() 
    317  
    318     def plotSlowmap( self, powmap, slomax, title="" ): 
     327        if outfile == None: 
     328            plt.show() 
     329        else: 
     330            plt.savefig( outfile.lower() ) 
     331 
     332    def plotSlowmap( self, powmap, slomax, title="", outfile=None ): 
    319333        fig = plt.figure( facecolor='white' ) 
    320334        ax = fig.add_subplot(1, 1, 1) 
     
    330344        ax.set_title( title ) 
    331345        plt.draw() 
    332         plt.show() 
     346        if outfile == None: 
     347            plt.show() 
     348        else: 
     349            plt.savefig( outfile.lower() ) 
    333350         
    334351 
Note: See TracChangeset for help on using the changeset viewer.