Changeset 1081


Ignore:
Timestamp:
06.01.2016 18:15:48 (4 years ago)
Author:
klaus
Message:

better plot interface in fk; mplplot bugfix and other label type

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

Legend:

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

    r1079 r1081  
    1010""" 
    1111 
     12import sys 
    1213from SeismicHandler.basics.command import BaseCommand 
    1314from SeismicHandler.modules.traces import traces_from_list 
     
    9798        "SCATTER", 
    9899        "POLAR", 
     100        "SLOWNESSMAP", 
    99101        "KMUNITS", 
    100102    ] 
     
    125127        slomax = 15. 
    126128        winlen = timeb - timea 
    127         plot_slowmap = True 
    128129        winfrac = 0.05 
    129130        prewhiten=0 
    130131        slosymb = None 
    131132        azsymb = None 
    132         plotscatter = False 
     133        plottype = None 
     134        timeplot_possible = False 
    133135        if len(self.parameters) > 3 and self.parameters[3]: 
    134136            frqmin = float( self.parameters[3] ) 
     
    147149        if self.qualifiers["WINLEN"]: 
    148150            winlen = float( self.qualifiers["WINLEN"] ) 
    149             plot_slowmap = False 
     151            timeplot_possible = True 
    150152        if self.qualifiers["WINFRAC"]: 
    151153            winfrac = float( self.qualifiers["WINFRAC"] ) 
    152154        if self.qualifiers["SCATTER"]: 
    153             plotscatter = True 
     155            plottype = "TimeScatter" 
    154156        if self.qualifiers["POLAR"]: 
    155             plotscatter = False 
     157            plottype = "TimePolar" 
     158        if self.qualifiers["SLOWNESSMAP"]: 
     159            plottype = "SlownessMap" 
    156160        if self.qualifiers["KMUNITS"]: 
    157161            cnvfactor = 1. 
     
    196200                'elevation' : elev, 
    197201            }) 
    198             trc.stats 
    199          
    200         store = None 
    201         if plot_slowmap: 
     202         
     203        # switch to slowness map if no time plot possible 
     204        if plottype != None: 
     205            if not timeplot_possible: 
     206                plottype = 'SlownessMap' 
     207         
     208        if plottype == 'SlownessMap': 
    202209            store = dumpSlow 
     210        else: 
     211            store = dumpIndex 
    203212 
    204213        fkout = array_processing( traces, store=store, **fkargs ) 
     
    215224            self.symbols.set( azsymb, result_az ) 
    216225         
    217          
    218         if plot_slowmap: 
     226        if plottype == 'SlownessMap': 
    219227            if cnvfactor == 1.: 
    220228                title = "slo:%5.3f  az:%5.1f" % (result_slo,result_az) 
     
    223231            self.plotSlowmap( np.fliplr(np.rot90(global_relpowmap*cnvfactor)), 
    224232                slomax*cnvfactor, title ) 
    225         elif plotscatter: 
     233        elif plottype == 'TimeScatter': 
    226234            self.plotScatter( fkout ) 
     235        elif plottype == 'TimePolar': 
     236            self.plotPolar( fkout, slomax*cnvfactor ) 
     237        elif plottype == None: 
     238            pass 
    227239        else: 
    228             self.plotPolar( fkout, slomax*cnvfactor ) 
     240            raise ShxError( "Program Bug (illegal plot type)", status=1111 ) 
    229241         
    230242        
     
    325337    global_relpowmap = relpow 
    326338 
     339def dumpIndex( relpow, abspow, i ): 
     340    print "slowstep %3d\r" % i, 
     341    sys.stdout.flush() 
     342 
    327343# Standard main part (same for all commands). 
    328344if __name__ == "__main__": 
  • SHX/trunk/SeismicHandler/commands/mplplot.py

    r1080 r1081  
    186186        TIME_HH:MM              :  hour and minute of trace start 
    187187        TIME_HH:MM:SS           :  hour, minute, sec of trace start 
     188        DAY                     :  date of trace start 
    188189        """ 
    189190        if labelspec == None or labelspec == 'STATION-COMP': 
     
    200201                lab = "00:00" 
    201202            return lab 
     203        elif labelspec == 'DAY': 
     204            return trc.stats.starttime.strftime( "%d-%h-%Y" ) 
    202205        else: 
    203206            return "**ILLSPEC**" 
Note: See TracChangeset for help on using the changeset viewer.